An official website of the United States government
ABSTRACT The thermal history and structure of the intergalactic medium (IGM) at $$z \ge 4$$ is an important boundary condition for reionization, and a key input for studies using the Ly $$\alpha$$ forest to constrain the masses of alternative dark matter candidates. Most such inferences rely on simulations that lack the spatial resolution to fully resolve the hydrodynamic response of IGM filaments and minihaloes to H i reionization heating. In this letter, we use high-resolution hydrodynamic + radiative transfer simulations to study how these affect the IGM thermal structure. We find that the adiabatic heating and cooling driven by the expansion of initially cold gas filaments and minihaloes sources significant small-scale temperature fluctuations. These likely persist in much of the IGM until $$z \le 4$$. Capturing this effect requires resolving the clumping scale of cold, pre-ionized gas, demanding spatial resolutions of $${\le} 2$$ $$h^{-1}$$kpc. Pre-heating of the IGM by X-rays can slightly reduce the effect. Our preliminary estimate of the effect on the Ly $$\alpha$$ forest finds that, at $$\log (k /[{\rm km^{-1} s}]) = -1.0$$, the Ly $$\alpha$$ forest flux power (at fixed mean flux) can increase $${\approx} 10~{{\ \rm per\ cent}}$$ going from 8 and 2 $$h^{-1}$$kpc resolution at $$z = 4{\!-\!}5$$ for gas ionized at $$z \ \lt\ 7$$. These findings motivate more careful analyses of how the effects studied here affect the Ly $$\alpha$$ forest.
more »
« less
Effect of separate initial conditions on the lyman-α forest in simulations
ABSTRACT Using a set of high resolution simulations, we quantify the effect of species-specific initial transfer functions on probes of the intergalactic medium (IGM) via the Lyman-α forest. We focus on redshifts 2–6, after H i reionization. We explore the effect of these initial conditions on measures of the thermal state of the low density IGM: the curvature, Doppler width cutoff, and Doppler width distribution. We also examine the matter and flux power spectrum, and potential consequences for constraints on warm dark matter models. We find that the curvature statistic is at most affected at the $$\approx 2{{\ \rm per\ cent}}$$ level at z = 6. The Doppler width cutoff parameters are affected by $$\approx 5{{\ \rm per\ cent}}$$ for the intercept, and $$\approx 8{{\ \rm per\ cent}}$$ for the fit slope, though this is subdominant to sample variation. The Doppler width distribution shows a $$\approx 30{{\ \rm per\ cent}}$$ effect at z = 3, however the distribution is not fully converged with simulation box size and resolution. The flux power spectrum is at most affected by $$\approx 5{{\ \rm per\ cent}}$$ at high redshift and small scales. We discuss numerical convergence with simulation parameters.
more »
« less
- Award ID(s):
- 1817256
- PAR ID:
- 10289365
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 503
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 1668 to 1679
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Recent observations from the EIGER JWST program have measured for the first time the quasar–galaxy cross-correlation function at $$z\approx 6$$. The autocorrelation function of faint $$z\approx 6$$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter haloes. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $$z\approx 2\!-\!4$$. We use a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $$z\approx 6$$ quasars and galaxies (i.e. their clustering properties and luminosity functions), and infer key quantities such as their luminosity–halo mass relation, the mass function of their host haloes, and their duty cycle/occupation fraction. Our key findings are (i) quasars reside on average in $$\approx 10^{12.5}\, {\rm M}_{\odot }$$ haloes (corresponding to $$\approx 5\sigma$$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($$\approx 1~{{\ \rm per\ cent}}$$); (iii) galaxies (that are bright in [O iii]) live in much smaller haloes ($$\approx 10^{10.9}\, {\rm M}_{\odot }$$) and have a larger duty cycle (occupation fraction) of $$\approx 13~{{\ \rm per\ cent}}$$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $$z\approx 2\!-\!6$$, and calls for new investigations of the role of quasar activity across cosmic time.more » « less
-
Abstract The density and temperature properties of the intergalactic medium (IGM) reflect the heating and ionization history during cosmological structure formation, and are primarily probed by the Ly α forest of neutral hydrogen absorption features in the observed spectra of background sources. We present the methodology and initial results from the Cholla IGM Photoheating Simulation (CHIPS) suite performed with the graphics process unit–accelerated Cholla code to study the IGM at high, uniform spatial resolution maintained over large volumes. In this first paper, we examine the IGM structure in CHIPS cosmological simulations that include IGM uniform photoheating and photoionization models where hydrogen reionization is completed early or by redshift z ∼ 6. Comparing with observations of the large- and small-scale Ly α transmitted flux power spectra P ( k ) at redshifts 2 ≲ z ≲ 5.5, the relative agreement of the models depends on scale, with the self-consistent Puchwein et al. IGM photoheating and photoionization model in good agreement with the flux P ( k ) at k ≳ 0.01 s km −1 at redshifts 2 ≲ z ≲ 3.5. On larger scales, the P ( k ) measurements increase in amplitude from z ∼ 4.6 to z ∼ 2.2, faster than the models, and lie in between the model predictions at 2.2 ≲ z ≲ 4.6 for k ≈ 0.002–0.01 s km −1 . We argue that the models could improve by changing the He ii photoheating rate associated with active galactic nuclei to reduce the IGM temperature at z ∼ 3. At higher redshifts, z ≳ 4.5, the observed flux P ( k ) amplitude increases at a rate intermediate between the models, and we argue that for models where hydrogen reionization is completed late ( z ∼ 5.5–6), resolving this disagreement will require inhomogeneous or “patchy” reionization. We then use an additional set of simulations to demonstrate that our results have numerically converged and are not strongly affected by varying cosmological parameters.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 mean free path of ionizing photons, λmfp, is a key factor in the photoionization of the intergalactic medium (IGM). At z ≳ 5, however, λmfp may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of λmfp that address this bias and extend up to z ∼ 6 for the first time. Our measurements at z ∼ 5 are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At z ∼ 6 we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $$\lambda _{\rm mfp} = 9.09^{+1.62}_{-1.28}$$ proper Mpc and $$0.75^{+0.65}_{-0.45}$$ proper Mpc (68 per cent confidence) at z = 5.1 and 6.0, respectively. The results at z = 5.1 are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At z = 6.0, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at z = 6.0 falls well below extrapolations from lower redshifts, indicating rapid evolution in λmfp over 5 < z < 6. This evolution disfavours models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by z = 6, but is qualitatively consistent with models wherein reionization completed at z = 6 or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20 per cent neutral at z = 6, although our measurement at z = 6.0 is even lower than these models prefer.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/mnras/stab555
