Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, 〈n2〉/〈n〉2, which corrects for unresolved fluctuations in gas density n on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution N-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale N-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally averaged, time-evolving clumping factor from the same high-resolution N-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20 per cent error, suggesting it may be a useful approximation.
This content will become publicly available on April 21, 2024
- Award ID(s):
- 2045600
- NSF-PAR ID:
- 10429294
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 522
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2047 to 2064
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT The ionizing photon escape fraction [Lyman continuum (LyC) fesc] of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshifts. Here we circumvent this uncertainty by inferring fesc using resolved (R > 4000) Lyman α (Lyα) profiles from the X-SHOOTER Lyα survey at z = 2 (XLS-z2). With empirically motivated criteria, we use Lyα profiles to select leakers ($f_{\mathrm{ esc}} > 20{{\ \rm per\ cent}}$) and non-leakers ($f_{\mathrm{ esc}} < 5{{\ \rm per\ cent}}$) from a representative sample of >0.2L* Lyman α emitters (LAEs). We use median stacked spectra of these subsets over λrest ≈ 1000–8000 Å to investigate the conditions for LyC fesc. Our stacks show similar mass, metallicity, MUV, and βUV. We find the following differences between leakers versus non-leakers: (i) strong nebular C iv and He ii emission versus non-detections; (ii) [O iii]/[O ii] ≈ 8.5 versus ≈3; (iii) Hα/Hβ indicating no dust versus E(B − V) ≈ 0.3; (iv) Mg ii emission close to the systemic velocity versus redshifted, optically thick Mg ii; and (v) Lyα fesc of ${\approx} 50{{\ \rm per\ cent}}$ versus ${\approx} 10{{\ \rm per\ cent}}$. The extreme equivalent widths (EWs) in leakers ([O iii]+$\mathrm{ H}\beta \approx 1100$ Å rest frame) constrain the characteristic time-scale of LyC escape to ≈3–10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers – extremely ionizing stellar populations, low column densities, a dust-free, high-ionization state interstellar medium (ISM) – occur simultaneously in the $f_{\rm esc} > 20{{\ \rm per\ cent}}$ stack, suggesting they are causally connected, and motivating why indicators like [O iii]/[O ii] may suffice to constrain fesc at z > 6 with the James Webb Space Telescope (JWST). The leakers comprise half of our sample, have a median LyC$f_{\rm esc} \approx 50{{\ \rm per\ cent}}$ (conservative range: $20\!-\!55{{\ \rm per\ cent}}$), and an ionizing production efficiency $\log ({\xi _{\rm {ion}}/\rm {Hz\ erg^{-1}}})\approx 25.9$ (conservative range: 25.7–25.9). These results show LAEs – the type of galaxies rare at z ≈ 2, but that become the norm at higher redshift – are highly efficient ionizers, with extreme ξion and prolific fesc occurring in sync.
-
more » « less
ABSTRACT The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here, we introduce an emissivity framework based on Lyα emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between production and escape of Lyα and LyC photons. Using constraints on the escape fractions of bright LAEs (LLyα > 0.2L*) at z ≈ 2 obtained from resolved Lyα profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at z ≈ 2–6, which is non-trivial given the strong evolution in both the star formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between z ≈ 6−9 under plausible assumptions. Within this framework, the >10 × rise in the UV population-averaged fesc between z ≈ 3–7 naturally arises due to the same phenomena that drive the growing LAE fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm {esc}}=50{{\ \rm per\ cent}}$, ξion = 1025.9 Hz erg−1 for half the bright LAEs), a highly ionizing minority of galaxies with MUV < −17 accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a ‘disco’ ionizing background. We conclude proposing tests to further develop our suggested Lyα-anchored formalism.
-
more » « less
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.
-
Abstract Becker et al. measured the mean free path of Lyman-limit photons in the intergalactic medium (IGM) at z = 6. The short value suggests that absorptions may have played a prominent role in reionization. Here we study physical properties of ionizing photon sinks in the wake of ionization fronts (I-fronts) using radiative hydrodynamic simulations. We quantify the contributions of gaseous structures to the Lyman-limit opacity by tracking the column-density distributions in our simulations. Within Δ t = 10 Myr of I-front passage, we find that self-shielding systems ( N H I > 10 17.2 cm −2 ) are comprised of two distinct populations: (1) overdensity Δ ∼ 50 structures in photoionization equilibrium with the ionizing background, and (2) Δ ≳ 100 density peaks with fully neutral cores. The self-shielding systems contribute more than half of the opacity at these times, but the IGM evolves considerably in Δ t ∼ 100 Myr as structures are flattened by pressure smoothing and photoevaporation. By Δ t = 300 Myr, they contribute ≲10% to the opacity in an average 1 Mpc 3 patch of the universe. The percentage can be a factor of a few larger in overdense patches, where more self-shielding systems survive. We quantify the characteristic masses and sizes of self-shielding structures. Shortly after I-front passage, we find M = 10 4 –10 8 M ⊙ and effective diameters d eff = 1–20 ckpc h −1 . These scales increase as the gas relaxes. The picture herein presented may be different in dark matter models with suppressed small-scale power.more » « less
