skip to main content

Explore Research Products in the NSF-PAR

Title: The thesan project: ionizing escape fractions of reionization-era galaxies
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
Award ID(s):
2007355 1814259
NSF-PAR ID:
10433779
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
520
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2757 to 2780
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT We present the escape fraction of hydrogen ionizing photons (fesc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z ≥ 5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 haloes in Mvir ∼ 108–$10^{12}\, M_{\odot }$ (M* ∼ 104–$10^{10}\, M_{\odot }$) at z = 5–12. We find the sample average 〈fesc〉increases with halo mass for Mvir ∼ 108–$10^{9.5}\, M_{\odot }$, becomes nearly constant for 109.5–$10^{11}\, M_{\odot }$, and decreases at ${\gtrsim}10^{11}\, M_{\odot }$. Equivalently, 〈fesc〉 increases with stellar mass up to $M_{\ast }\sim 10^8\, M_{\odot }$ and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate 〈fesc〉 ∼ 0.2 for galaxies at $M_{\ast }\sim 10^8\, M_{\odot }$. Nearly half of the escaped ionizing photons come from stars 1–3 Myr old and the rest from stars 3–10 Myr old. Binaries only have a modest effect, boosting 〈fesc〉 by ∼25–35 per cent and the number of escaped photons by 60–80 per cent. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of 〈fesc〉 at the high-mass end is due to dust attenuation, while at the low-mass end, 〈fesc〉 decreases owing to inefficient star formation and hence feedback. At fixed mass, 〈fesc〉 tends to increase with redshift. Although the absolute 〈fesc〉does not fully converge with resolution in our simulations, the mass- and redshift-dependence of 〈fesc〉 is likely robust. Our simulations produce sufficient ionizing photons for cosmic reionization. 
    more » « less
  2. ; ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    The observability of Lyα emitting galaxies (LAEs) during the Epoch of Reionization can provide a sensitive probe of the evolving neutral hydrogen gas distribution, thus setting valuable constraints to distinguish different reionization models. In this study, we utilize the new thesan suite of large-volume ($L_\text{box} = 95.5\, \text{cMpc}$) cosmological radiation-hydrodynamic simulations to directly model the Lyα emission from individual galaxies and the subsequent transmission through the intergalactic medium. thesan combines the arepo-rt radiation-hydrodynamic solver with the IllustrisTNG galaxy formation model and includes high- and medium-resolution simulations designed to investigate the impacts of halo-mass-dependent escape fractions, alternative dark matter models, and numerical convergence. We find important differences in the Lyα transmission based on reionization history, bubble morphology, frequency offset from line centre, and galaxy brightness. For a given global neutral fraction, Lyα transmission reduces when low-mass haloes dominate reionization over high-mass haloes. Furthermore, the variation across sightlines for a single galaxy is greater than the variation across all galaxies. This collectively affects the visibility of LAEs, directly impacting observed Lyα luminosity functions (LFs). We employ Gaussian Process Regression using SWIFTEmulator to rapidly constrain an empirical model for dust escape fractions and emergent spectral-line profiles to match observed LFs. We find that dust strongly impacts the Lyα transmission and covering fractions of MUV ≲ −19 galaxies in $M_\text{vir} \gtrsim 10^{11}\, \text{M}_{\bigodot }$ haloes, such that the dominant mode of removing Lyα photons in non-LAEs changes from low-IGM transmission to high dust absorption around z ∼ 7.

     
    more » « less
  3. ; ; ; ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    ABSTRACT

    We measure escape fractions, fesc, of ionizing radiation from galaxies in the sphinx suite of cosmological radiation-hydrodynamical simulations of reionization, resolving haloes with $M_{\rm vir}\gtrsim 7.5 \times 10^7 \ {\rm {M}_{\odot }}$ with a minimum cell width of ≈10 pc. Our new and largest 20 co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times 1011 M⊙. The simulated galaxies agree well with observational constraints of the ultraviolet (UV) luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over time-scales of a few Myr, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the intergalactic medium. Statistically, fesc peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies (M* ≈ 107 M⊙, M1500 ≈ −17, Z ≲ 5 × 10−3 Z⊙), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to sphinx, galaxies too dim to be yet observed, with ${M_{1500}}\gtrsim -17$, provide about 55 per cent of the photons contributing to reionization. The global averaged fesc naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.

     
    more » « less
  4. ; ; ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    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
  5. ; ; ; ; ( , Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT We study the escape fraction of ionizing photons (fesc) in two cosmological zoom-in simulations of galaxies in the reionization era with halo mass Mhalo ∼ 1010 and $10^{11}\, \mathrm{ M}_{\odot }$ (stellar mass M* ∼ 107 and $10^9\, \mathrm{ M}_{\odot }$) at z = 5 from the Feedback in Realistic Environments project. These simulations explicitly resolve the formation of proto-globular clusters (GCs) self-consistently, where 17–39 per cent of stars form in bound clusters during starbursts. Using post-processing Monte Carlo radiative transfer calculations of ionizing radiation, we compute fesc from cluster stars and non-cluster stars formed during a starburst over ∼100 Myr in each galaxy. We find that the averaged fesc over the lifetime of a star particle follows a similar distribution for cluster stars and non-cluster stars. Clusters tend to have low fesc in the first few Myr, presumably because they form preferentially in more extreme environments with high optical depths; the fesc increases later as feedback starts to destroy the natal cloud. On the other hand, some non-cluster stars formed between cluster complexes or in the compressed shells at the front of a superbubble can also have high fesc. We find that cluster stars on average have comparable fesc to non-cluster stars. This result is robust across several star formation models in our simulations. Our results suggest that the fraction of ionizing photons from proto-GCs to cosmic reionization is comparable to the cluster formation efficiencies in high-redshift galaxies and thus proto-GCs likely contribute an appreciable fraction of photons but are not the dominant sources for reionization. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email