publications

2024

1. Semi-analytic modelling of Pop. III star formation and metallicity evolution – I. Impact on the UV luminosity functions at z = 9–16

   Emanuele M Ventura, Yuxiang Qin, Sreedhar Balu, and 1 more author

   Monthly Notices of the Royal Astronomical Society, Feb 2024

   Abs DOI

   ABSTRACT We implemented Population III (Pop. III) star formation in mini-haloes within the Meraxes semi-analytic galaxy formation and reionization model, run on top of a N-body simulation with L = 10 h−1 cMpc with 20483 particles resolving all dark matter haloes down to the mini-haloes (∼105 M⊙). Our modelling includes the chemical evolution of the IGM, with metals released through supernova-driven bubbles that expand according to the Sedov–Taylor model. We found that SN-driven metal bubbles are generally small, with radii typically of 150 ckpc at z = 6. Hence, the majority of the first galaxies are likely enriched by their own star formation. However, as reionization progresses, the feedback effects from the UV background become more pronounced, leading to a halt in star formation in low-mass galaxies, after which external chemical enrichment becomes more relevant. We explore the sensitivity of the star formation rate density and stellar mass functions to the unknown values of free parameters. We also discuss the observability of Pop. III dominated systems with JWST, finding that the inclusion of Pop. III galaxies can have a significant effect on the total UV luminosity function at z = 12–16. Our results support the idea that the excess of bright galaxies detected with JWST might be explained by the presence of bright top-heavy Pop. III dominated galaxies without requiring an increased star formation efficiency.

2. The impact of ionization morphology and X-ray heating on the cosmological 21-cm skew spectrum

   Jadien H Cook, Sreedhar Balu, Bradley Greig, and 4 more authors

   Monthly Notices of the Royal Astronomical Society, Mar 2024

   Abs DOI

   ABSTRACT The cosmological 21-cm signal offers a potential probe of the early Universe and the first ionizing sources. Current experiments probe the spatially dependent variance (Gaussianity) of the signal through the power spectrum (PS). The signal, however, is expected to be highly non-Gaussian due to the complex topology of reionization and X-ray heating. We investigate the non-Gaussianities of X-ray heating and reionization, by calculating the skew spectrum (SS) of the 21-cm signal using meraxes, which couples a semi-analytic galaxy population with seminumerical reionization simulations. The SS is the cross-spectrum of the quadratic temperature brightness field with itself. We generate a set of seven simulations from z = 30 to 5, varying the halo mass threshold for hosting star formation, the X-ray luminosity per star formation rate, and the minimum X-ray energy escaping host galaxies. We find the SS is predominantly negative as a function of redshift, transitioning to positive towards the start of reionization, and peaking during the mid-point of reionization. We do not see a negative dip in the SS during reionization, likely due to the specifics of modelling ionization sources. We normalize the SS by the PS during reionization isolating the non-Gaussianities. We find a trough (\k}sim }, 0.1}, }textrm {Mpc}^{-1}\) and peak (\k}sim }, 0.4{-}1}, }textrm {Mpc}^{-1}\) in the normalized SS during the mid-to-late periods of reionization. These correlate to the ionization topology, and neutral islands in the IGM. We calculate the cosmic variance of the normalized SS, and find these features are detectable in the absence of foregrounds with the SKA_LOW.

2023

1. Thermal and reionization history within a large-volume semi-analytic galaxy formation simulation

   Sreedhar Balu, Bradley Greig, Yisheng Qiu, and 4 more authors

   Monthly Notices of the Royal Astronomical Society, Apr 2023

   Abs DOI

   We predict the 21-cm global signal and power spectra during the Epoch of Reionization using the meraxes semi-analytic galaxy formation and reionization model, updated to include X-ray heating and thermal evolution of the intergalactic medium. Studying the formation and evolution of galaxies together with the reionization of cosmic hydrogen using semi-analytic models (such as M eraxes) requires N-body simulations within large volumes and high-mass resolutions. For this, we use a simulation of side-length 210 h−1 Mpc with 43203 particles resolving dark matter haloes to masses of \5}times 10^8 }rm{ }h^{-1}}, }mathrm{M_}odot }\. To reach the mass resolution of atomically cooled galaxies, thought to be the dominant population contributing to reionization, at z = 20 of {}sim 2}times 10^7 }text{ }h^{-1}}, }mathrm{M_}odot } we augment this simulation using the darkforest Monte Carlo merger tree algorithm (achieving an effective particle count of ∼1012). Using this augmented simulation, we explore the impact of mass resolution on the predicted reionization history as well as the impact of X-ray heating on the 21-cm global signal and the 21-cm power spectra. We also explore the cosmic variance of 21-cm statistics within 703 h−3 Mpc3 sub-volumes. We find that the midpoint of reionization varies by Δz ∼ 0.8 and that the cosmic variance on the power spectrum is underestimated by a factor of 2–4 at k ∼ 0.1–0.4 Mpc−1 due to the non-Gaussian nature of the 21-cm signal. To our knowledge, this work represents the first model of both reionization and galaxy formation which resolves low-mass atomically cooled galaxies while simultaneously sampling sufficiently large scales necessary for exploring the effects of X-rays in the early Universe.

2. Fisher matrix forecasts on the astrophysics of galaxies during the epoch of reionization from the 21-cm power spectra

   Sreedhar Balu, Bradley Greig, and J Stuart B Wyithe

   Monthly Notices of the Royal Astronomical Society, Aug 2023

   Abs DOI

   ABSTRACT The hyperfine 21-cm transition of neutral hydrogen from the early Universe (z > 5) is a sensitive probe of the formation and evolution of the first luminous sources. Using the Fisher matrix formalism we explore the complex and degenerate high-dimensional parameter space associated with the high-z sources of this era and forecast quantitative constraints from a future 21-cm power spectrum (21-cm PS) detection. This is achieved using {}rm {}small ERAXES} a coupled semi-analytic galaxy formation model and reionization simulation, applied to an N-body halo merger tree with a statistically complete population of all atomically cooled galaxies out to z ∼ 20. Our mock observation assumes a 21-cm detection spanning z ∈ [5, 24] from a 1000 h mock observation with the forthcoming Square Kilometre Array, and is calibrated with respect to ultraviolet luminosity functions (UV LFs) at z ∈ [5, 10], the optical depth of CMB photons to Thompson scattering from Planck, and various constraints on the IGM neutral fraction at z > 5. In this work, we focus on the X-ray luminosity, ionizing UV photon escape fraction, star formation, and supernova feedback of the first galaxies. We demonstrate that it is possible to recover five of the eight parameters describing these properties with better than 50 per cent precision using just the 21-cm PS. By combining with UV LFs, we are able to improve our forecast, with five of the eight parameters constrained to better than 10 per cent (and all below 50 per cent).

3. Implications of z ≳ 12 JWST galaxies for galaxy formation at high redshift

   Yuxiang Qin, Sreedhar Balu, and J Stuart B Wyithe

   Monthly Notices of the Royal Astronomical Society, Sep 2023

   Abs DOI

   ABSTRACT Using a semi-analytic galaxy formation model, we study analogues of eight z ≳ 12 galaxies recently discovered by James Webb Space Telescope (JWST). We select analogues from a cosmological simulation with a (311 cMpc)3 volume and an effective particle number of 1012 enabling the resolution of every atomic-cooling galaxy at z ≤ 20. We vary model parameters to reproduce the observed ultraviolet (UV) luminosity function at 5 < z < 13, aiming for a statistically representative high-redshift galaxy mock catalogue. Using the forward-modelled JWST photometry, we identify analogues from this catalogue and study their properties as well as possible evolutionary paths and local environment. We find faint JWST galaxies (MUV ≳ − 19.5) to remain consistent with the standard galaxy formation model and that our fiducial catalogue includes large samples of their analogues. The properties of these analogues broadly agree with conventional spectral energy distribution-fitting results, except for having systematically lower redshifts due to the evolving ultraviolet luminosity function, and for having higher specific star formation rates as a result of burstier histories in our model. On the other hand, only a handful of bright galaxy analogues can be identified for the observed z ∼ 12 galaxies. Moreover, in order to reproduce the z ≳ 16 JWST galaxy candidates, boosting star-forming efficiencies through reduced feedback regulation and increased gas depletion rate is necessary relative to models of lower redshift populations. This suggests star formation in the first galaxies could differ significantly from their lower redshift counterparts. We also find that these candidates are subject to low-redshift contamination, which is present in our fiducial results as both the dusty or quiescent galaxies at z ∼ 5.