An official website of the United States government
Abstract A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work, we calculate the ionizing escape fractions (fesc) of reionization-era galaxies from the state-of-the-art thesan simulations, which combine an accurate radiation-hydrodynamic solver (arepo-rt) with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($$L_\text{box} = 95.5\, \text{cMpc}$$). This allows the formation of numerous massive haloes ($$M_\text{halo} \gtrsim 10^{10}\, {\text{M}_{\odot }}$$), which are often statistically underrepresented in previous studies but are believed to be important to achieving rapid reionization. We find that low-mass galaxies ($$M_\text{stars} \lesssim 10^7\, {\text{M}_{\odot }}$$) are the main drivers of reionization above z ≳ 7, while high-mass galaxies ($$M_\text{stars} \gtrsim 10^8\, {\text{M}_{\odot }}$$) dominate the escaped ionizing photon budget at lower redshifts. We find a strong dependence of fesc on the effective star formation rate (SFR) surface density defined as the SFR per gas mass per escape area, i.e. $$\bar{\Sigma }_\text{SFR} = \text{SFR}/M_\text{gas}/R_{200}^2$$. The variation in halo escape fractions decreases for higher mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. Dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global fesc depends on the dust model. Finally, active galactic nuclei are unimportant for reionization in thesan and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells.
more »
« less
Introducing the thesan project: radiation-magnetohydrodynamic simulations of the epoch of reionization
ABSTRACT We introduce the thesan project, a suite of large volume ($$L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $$3.1 \times 10^6\, {\rm M_\odot }$$ and $$5.8 \times 10^5\, {\rm M_\odot }$$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.
more »
« less
- Award ID(s):
- 2307699
- PAR ID:
- 10363140
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 511
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 4005-4030
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($$M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$$), run to z = 0 with $$30\, \mathrm{M}_{\odot }$$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $$M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $$M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($$M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core $${\gt}200\rm \, pc$$, while lower mass UFDs exhibit cusps down to $${\lesssim}100\rm \, pc$$, as expected from energetic arguments. Our dwarfs with $$M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).more » « less
-
ABSTRACT We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (∼1000 galaxies with $$M_{\rm star}\gt 10^8\, M_\odot$$ at z = 0) at a resolution ($$\Delta {}x\sim {}20\, {\rm pc}$$ , $$m_{\rm b}\sim {}6\times {}10^4\, M_\odot$$ ) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L = 22.1 Mpc) thanks to its exceptionally high dynamic range (≳106) and the inclusion of multichannel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that star formation rates, gas masses, and metallicities of simulated galaxies with $$M_{\rm star}\lt 10^{10.5-11}\, M_\odot$$ broadly agree with observations. These galaxy scaling relations extend to low masses ($$M_{\rm star}\sim {}10^7\, M_\odot$$ ) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z ∼ 0–5. At low z , FIREbox predicts a peak in the stellar-mass–halo-mass relation but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a ΛCDM Universe while also highlighting modelling challenges to be addressed in next-generation galaxy simulations.more » « less
-
ABSTRACT The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at $$z \gtrsim 6$$ is an open question for understanding which galaxies primarily drove cosmic reionization. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful ‘near-far’ technique for studying the reionization-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass ($$M_{\rm star}\lesssim 10^{9}\, \rm {M_\odot }$$) galaxies across 11 Milky Way (MW)- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterize their progenitors at $$z=6\!-\!9$$, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with $$M_{\rm star}\sim 10^5\, \rm {M_\odot }$$ ($$\sim 10^9\, \rm {M_\odot }$$) had $$\approx 1$$ ($$\approx 30$$) progenitors at $$z\approx 7$$, and its main progenitor comprised $$\approx 100~{{\ \rm per\ cent}}$$ ($$\approx 10~{{\ \rm per\ cent}}$$) of the total stellar mass of all its progenitors at $$z\approx 7$$. We show that although only $$\sim 15~{{\ \rm per\ cent}}$$ of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at $$z \sim 6 \!-\! 9$$. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just low-mass galaxies at $z = 0$ and not the MW/M31 hosts correctly recovers the slope of the SMF down to $$M_{\rm star} \sim 10^{4.5}\, \rm {{\rm M}_{\odot }}$$ at $$z \gtrsim 6$$. Thus, we validate the ‘near-far’ approach as an unbiased tool for probing low-mass reionization-era galaxies.more » « less
-
ABSTRACT The connection between the escape fraction of ionizing radiation (fesc) and the properties of galaxies, such as stellar mass ($$\rm M_{\rm *}$$), age, star-formation rate (SFR), and dust content, are key inputs for reionization models, but many of these relationships remain untested at high redshift. We present an analysis of a sample of 96 $$z$$ ∼ 3 galaxies from the Keck Lyman Continuum Spectroscopic Survey (KLCS). These galaxies have both sensitive Keck/LRIS spectroscopic measurements of the Lyman continuum (LyC) region, and multiband photometry that places constraints on stellar population parameters. We construct composite spectra from subsamples binned as a function of galaxy property and quantify the ionizing-photon escape for each composite. We find a significant anti-correlation between fesc and $$\rm M_{\rm *}$$, consistent with predictions from cosmological zoom-in simulations. We also find significant anti-correlation between fesc and E(B−V), encoding the underlying physics of LyC escape in our sample. We also find no significant correlation between fesc and either stellar age or specific SFR (= SFR/$$\rm M_{\rm *}$$), challenging interpretations that synchronize recent star formation and favorable conditions for ionizing escape. The galaxy properties now shown to correlate with fesc in the KLCS are Lyα equivalent width, UV Luminosity, $$\rm M_{\rm *}$$, SFR, and E(B−V), but not age or sSFR. This comprehensive analysis of galaxy properties and LyC escape at high redshift will be used to guide future models and observations of the reionization epoch.more » « less
