skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The hydrodynamic response of small-scale structure to reionization drives large IGM temperature fluctuations that persist to z = 4
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
Award ID(s):
2045600
PAR ID:
10616018
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Monthly Notices of the Royal Astronomical Society
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society: Letters
Volume:
533
Issue:
1
ISSN:
1745-3925
Page Range / eLocation ID:
L100 to L106
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. null (Ed.)
    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
  3. ABSTRACT We explore the use of deep learning to infer the temperature of the intergalactic medium from the transmitted flux in the high-redshift Ly α forest. We train neural networks on sets of simulated spectra from redshift z = 2–3 outputs of cosmological hydrodynamic simulations, including high-temperature regions added in post-processing to approximate bubbles heated by He ii reionization. We evaluate how well the trained networks are able to reconstruct the temperature from the effect of Doppler broadening in the simulated input Ly α forest absorption spectra. We find that for spectra with high resolution (10 $$\, {\rm km}\, {\rm s}^{-1}$$ pixel) and moderate signal-to-noise ratio (20–50), the neural network is able to reconstruct the intergalactic medium temperature smoothed on scales of $$\sim 6 \, h^{-1}\, {\rm Mpc}$$ quite well. Concentrating on discontinuities, we find that high-temperature regions of width $$25 \, h^{-1}\, {\rm Mpc}$$ and temperature $$20\, 000$$ K can be fairly easily detected and characterized. We show an example where multiple sightlines are combined to yield tomographic images of hot bubbles. Deep learning techniques may be useful in this way to help us understand the complex temperature structure of the intergalactic medium around the time of helium reionization. 
    more » « less
  4. Abstract The absorption by neutral hydrogen in the intergalactic medium (IGM) produces the Ly α forest in the spectra of quasars. The Ly α forest absorbers have a broad distribution of neutral hydrogen column density N H I and Doppler b parameter. The narrowest Ly α absorption lines (of lowest b ) with neutral hydrogen column density above ∼10 13 cm −2 are dominated by thermal broadening, which can be used to constrain the thermal state of the IGM. Here we constrain the temperature-density relation T = T 0 ( ρ / ρ ¯ ) γ − 1 of the IGM at 1.6 < z < 3.6 by using N H I and b parameters measured from 24 high-resolution and high-signal-to-noise quasar spectra and by employing an analytic model to model the N H I -dependent low- b cutoff in the b distribution. In each N H I bin, the b cutoff is estimated using two methods, one non-parametric method from computing the cumulative b distribution and a parametric method from fitting the full b distribution. We find that the IGM temperature T 0 at the mean gas density ρ ¯ shows a peak of ∼1.5 × 10 4 K at z ∼ 2.7–2.9. At redshift higher than this, the index γ approximately remains constant, and it starts to increase toward lower redshifts. The evolution in both parameters is in good agreement with constraints from completely different approaches, which signals that He ii reionization completes around z ∼ 3. 
    more » « less
  5. Abstract We investigate the metal species associated with the Ly α forest in eBOSS quasar spectra. Metal absorption lines are revealed in stacked spectra from cross-correlating the selected Ly α absorbers in the forest and the flux fluctuation field. Up to 13 metal species are identified as being associated with relatively strong Ly α absorbers (those with flux fluctuations − 1.0 < δ Ly α < − 0.6 and with a neutral hydrogen column density of ∼ 10 15−16 cm −2 ) over the absorber redshift range of 2 < z abs < 4. The column densities of these species decrease toward higher redshift and for weaker Ly α absorbers. From modeling the column densities of various species, we find that the column density pattern suggests contributions from multiple gas components, both in the circumgalactic medium (CGM) and the intergalactic medium (IGM). While the low-ionization species (e.g., C ii , Si ii , and Mg ii ) can be explained by high-density, cool gas ( T ∼ 10 4 K) from the CGM, the high-ionization species may reside in low-density or high-temperature gas in the IGM. The measurements provide inputs for modeling the metal contamination in the Ly α forest baryon acoustic oscillation measurements. Comparisons with metal absorptions in high-resolution quasar spectra and hydrodynamic galaxy formation simulations can further elucidate the physical conditions of these Ly α absorbers. 
    more » « less