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.
- Award ID(s):
- 2045600
- 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 Reionization leads to large spatial fluctuations in the intergalactic temperature that can persist well after its completion. We study the imprints of such fluctuations on the $z$ ∼ 5 Ly α forest flux power spectrum using a set of radiation-hydrodynamic simulations that model different reionization scenarios. We find that large-scale coherent temperature fluctuations bring ${\sim}20\text{--}60{{\ \rm per\ cent}}$ extra power at k ∼ 0.002 s km−1, with the largest enhancements in the models where reionization is extended or ends the latest. On smaller scales (k ≳ 0.1 s km−1), we find that temperature fluctuations suppress power by ${\lesssim}10{{\ \rm per\ cent}}$. We find that the shape of the power spectrum is mostly sensitive to the reionization mid-point rather than temperature fluctuations from reionization’s patchiness. However, for all of our models with reionization mid-points of $z$ ≤ 8 ($z$ ≤ 12), the shape differences are ${\lesssim}20{{\ \rm per\ cent}}$ (${\lesssim}40{{\ \rm per\ cent}}$) because of a surprisingly well-matched cancellation between thermal broadening and pressure smoothing that occurs for realistic thermal histories. We also consider fluctuations in the ultraviolet background, finding their impact on the power spectrum to be much smaller than temperature fluctuations at k ≳ 0.01 s km−1. Furthermore, we compare our models to power spectrum measurements, finding that none of our models with reionization mid-points of $z$ < 8 is strongly preferred over another and that all of our models with mid-points of $z$ ≥ 8 are excluded at 2.5σ. Future measurements may be able to distinguish between viable reionization models if they can be performed at lower k or, alternatively, if the error bars on the high-k power can be reduced by a factor of 1.5.
-
more » « less
ABSTRACT We use deep imaging from the JWST Advanced Deep Extragalactic Survey (JADES) to study the evolution of the ionizing photon production efficiency, ξion. We estimate ξion for a sample of 677 galaxies at z ∼ 4–9 using NIRCam (Near-Infrared Camera) photometry. Specifically, combinations of the medium and wide bands F335M–F356W and F410M–F444W to constrain emission lines that trace ξion: Hα and [O iii]. Additionally, we use the spectral energy distribution fitting code prospector to fit all available photometry and infer galaxy properties. The flux measurements obtained via photometry are consistent with FRESCO (First Reionisation Epoch Spectroscopic Complete Survey) and NIRSpec-derived fluxes. Moreover, the emission-line-inferred measurements are consistent with the prospector estimates. We also confirm the observed ξion trend with redshift and MUV, and find: log ξion(z, MUV) = (0.05 ± 0.02)z + (0.11 ± 0.02)MUV + (27.33 ± 0.37). We use prospector to investigate correlations of ξion with other galaxy properties. We see a clear correlation between ξion and burstiness in the star formation history of galaxies, given by the ratio of recent to older star formation, where burstiness is more prevalent at lower stellar masses. We also convolve our ξion relations with luminosity functions from the literature, and constant escape fractions of 10 per cent and 20 per cent, to place constraints on the cosmic ionizing photon budget. By combining our results, we find that if our sample is representative of the faint low-mass galaxy population, galaxies with bursty star formation are efficient enough in producing ionizing photons and could be responsible for the reionization of the Universe.
-
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 Spitzer/Infrared Array Camera (IRAC) imaging has revealed that the brightest z ∼ 7−8 galaxies often exhibit young ages and strong nebular line emission, hinting at high ionizing efficiency among early galaxies. However, IRAC’s limited sensitivity has long hindered efforts to study the fainter, more numerous population often thought largely responsible for reionization. Here, we use Cosmic Evolution Early Release Science (CEERS) JWST/NIRCam data to characterize 116 ultraviolet (UV)-faint (median MUV = −19.5) z ∼ 6.5−8 galaxies. The spectral energy distributions are typically dominated by young (∼10–50 Myr), low-mass (M* ∼ 108 M⊙) stellar populations, and we find no need for extremely high stellar masses (∼1011 M⊙). Considering previous studies of UV-bright (MUV ∼ −22) z ∼ 7−8 galaxies, we find evidence for a strong (5–10 times) increase in specific star formation rate (sSFR) toward lower luminosities (median sSFR = 103 Gyr−1 in CEERS). The larger sSFRs imply a more dominant contribution from OB stars in the relatively numerous UV-faint population, perhaps suggesting that these galaxies are very efficient ionizing agents (median ξion = 1025.7 erg−1 Hz). In spite of the much larger sSFRs, we find little increase in [O iii] + H β equivalent widths towards fainter MUV (median ≈780 $\mathrm{\mathring{A}}$). If confirmed, this may indicate that a substantial fraction of our CEERS galaxies possess extremely low metallicities (≲3 per cent Z⊙) where [O iii] emission is suppressed. Alternatively, high ionizing photon escape fractions or bursty star formation histories can also weaken the nebular lines in a subset of our sample. While the majority of galaxies in our sample are very blue (median β = −2.0), we identify a significant tail of very dusty galaxies (β ∼ −1) at ≈0.5$L_\mathrm{UV}^\ast$ which may contribute significantly to the z ∼ 7−8 star formation rate density.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stad1057
