skip to main content


Title: Constraining the Temperature-density Relation of the Inter-galactic Medium from Analytically Modeling Lyα Forest Absorbers
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
Award ID(s):
2007499
NSF-PAR ID:
10420272
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Research in Astronomy and Astrophysics
Volume:
23
Issue:
4
ISSN:
1674-4527
Page Range / eLocation ID:
045007
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. ABSTRACT We compare a sample of five high-resolution, high S/N  Ly α forest spectra of bright 6 < z < ∼6.5 QSOs aimed at spectrally resolving the last remaining transmission spikes at z > 5 with those obtained from mock absorption spectra from the Sherwoodand Sherwood–Relics simulation suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile-fitting procedure for the inverted transmitted flux, 1 − F, similar to the widely used Voigt profile fitting of the transmitted flux F at lower redshifts, to characterize the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominent in low temperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution (≤ 8  km s−1) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of T0 = (11 000 ± 1600, 10 500 ± 2100, 12 000 ± 2200) K at z = (5.4, 5.6, 5.8). This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to z ∼ 5.3, the width distribution of the observed transmission spikes is consistent with the range of T0 caused by spatial fluctuations in the temperature–density relation. 
    more » « less
  3. Abstract The Ly-α forest is the large-scale structure probe for which we appear to have modeling control to the highest wavenumbers. This makes the Ly-α forest of great interest for constraining the warmness/fuzziness of dark matter and the timing of reionization processes. However, the standard statistic, the Ly-α forest power spectrum, is unable to strongly constrain the IGM temperature-density relation, and this inability further limits how well other high wavenumber-sensitive parameters can be constrained. With the aim of breaking these degeneracies, we measure the power spectrum of the Ly-β forest and its cross correlation with the coeval Ly-α forest using the one hundred spectra of z = 3.5 − 4.5 quasars in the VLT/X-Shooter XQ-100 Legacy Survey, motivated by the Ly-β transition’s smaller absorption cross section that makes it sensitive to somewhat higher densities relative to the Ly-α transition. Our inferences from this measurement for the IGM temperature-density relation appear to latch consistently onto the recent tight lower-redshift Ly-α forest constraints. The z = 3.4 − 4.7 trends we find using the Ly-α–Ly-β cross correlation show a flattening of the slope of the temperature-density relation with decreasing redshift. This is the trend anticipated from ongoing He ii reionization and there being sufficient time to reach the asymptotic temperature-density slope after hydrogen reionization completes. Furthermore, our measurements provide a consistency check on IGM models that explain the Ly-α forest, with the cross correlation being immune to systematics that are uncorrelated between the two forests, such as metal line contamination. 
    more » « less
  4. 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
  5. Abstract

    Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low (x¯HI103) or close to unity (x¯HI1). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints onx¯HIatz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits ofx¯HI(z=6.3)<0.79±0.04(1σ),x¯HI(z=6.5)<0.87±0.03(1σ), andx¯HI(z=6.7)<0.940.09+0.06(1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

     
    more » « less