skip to main content


Title: Effects of Photoionization and Photoheating on Lyα Forest Properties from Cholla Cosmological Simulations
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
Award ID(s):
1828315
NSF-PAR ID:
10313043
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
912
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The presence of excess scatter in the Ly-α forest at z ∼ 5.5, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionization. However, low data quality and systematic uncertainties complicate the use of Ly-α transmission as a precision probe of reionization’s end stages. In this paper, we assemble a sample of 67 quasar sightlines at z > 5.5 with high signal-to-noise ratios of >10 per ≤15 km s−1 spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR > 20 per pixel) spectra of 30 quasars at z > 5.7. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-α systems. We present improved measurements of the mean Ly-α transmission over 4.9 < z < 6.1. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-α transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to z ≤ 5.2 (<1σ), and mild tension (∼2.5σ) at z = 5.3. Homogeneous UVB models are ruled out by excess Ly-α transmission scatter at z ≥ 5.4 with high confidence (>3.5σ). Our results indicate that reionization-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least z = 5.3 (t = 1.1 Gyr after the big bang). This is further evidence for a late end to reionization. 
    more » « less
  2. ABSTRACT

    Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly α transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly α profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at z ∼ 7 remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [C ii]-based redshifts from ALMA. Here, we report MMT/Binospec Ly α spectroscopy of eight UV-bright (MUV ∼ −22) galaxies at z ≃ 7 selected from the ALMA REBELS survey. We detect Ly α in four of eight galaxies and use the [C ii] systemic redshifts to investigate the Ly α velocity profiles. The Ly α lines are significantly redshifted from systemic (average velocity offset = 223 km s–1) and broad (FWHM ≈ 300–650 km s−1), with two sources showing emission extending to ≈750 km s−1. We find that the broadest Ly α profiles are associated with the largest [C ii] line widths, suggesting a potential link between the Ly α FWHM and the dynamical mass. Since Ly α photons at high velocities transmit efficiently through the z = 7 IGM, our data suggest that velocity profiles play a significant role in boosting the Ly α visibility of the most UV-luminous reionization-era galaxies.

     
    more » « less
  3. Abstract We present a new investigation of the intergalactic medium (IGM) near the end of reionization using “dark gaps” in the Ly α forest. Using spectra of 55 QSOs at z em > 5.5, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly α forest where the transmission averaged over 1 comoving h −1 Mpc bins falls below 5%. Nine ultralong ( L > 80 h −1 Mpc) dark gaps are identified at z < 6. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than 30 h −1 Mpc, F 30 , as a function of redshift. We measure F 30 ≃ 0.9, 0.6, and 0.15 at z = 6.0, 5.8, and 5.6, respectively, with the last of these long dark gaps persisting down to z ≃5.3. Comparing our results with predictions from hydrodynamical simulations, we find that the data are consistent with models wherein reionization extends significantly below redshift six. Models wherein the IGM is essentially fully reionized that retain large-scale fluctuations in the ionizing UV background at z ≲6 are also potentially consistent with the data. Overall, our results suggest that signatures of reionization in the form of islands of neutral hydrogen and/or large-scale fluctuations in the ionizing background remain present in the IGM until at least z ≃ 5.3. 
    more » « less
  4. ABSTRACT

    Revealing the cosmic hydrogen reionization history is one of the main goals of the modern cosmology. z > 5 quasars (QSOs) have been used as back-lights to investigate the evolution of the intervening intergalactic medium (IGM) during the cosmic reionization since their first discovery. However, due to the small population of luminous QSOs (∼130 QSOs known to date), a tight constraint on the reionization history has not yet been placed. In this work, we aim to tighten the constraint using the 93 QSOs (5.5 < z < 7.1) recently discovered in the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQS) project. This is the largest QSO sample used to constrain the epoch of reionization. We measure the mean IGM Ly α transmission and the QSO near-zone size using the UV spectra of these QSOs. The mean IGM Ly α transmission rises above zero at z ≲ 6, indicating the end of the reionization. The near-zone sizes of the SHELLQs QSOs are consistent with sizes spanned by QSOs of lifetime tq ∼ 1–100 Myr in simulations. Due to the scatter created by the low signal-to-noise spectra and large Ly α redshift uncertainty, we cannot conclude whether the redshift evolution of the near-zone size is affected by the reionization effect.

     
    more » « less
  5. 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