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, more » 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. « less
Authors:
; ; ;
Award ID(s):
1828315
Publication Date:
NSF-PAR ID:
10313043
Journal Name:
The Astrophysical Journal
Volume:
912
Issue:
2
ISSN:
0004-637X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 ablemore »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.« less
  2. 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 deepmore »(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.« less
  3. Abstract We explore how the assumption of ionization equilibrium modulates the modeled intergalactic medium at the end of the hydrogen epoch of reionization using the cosmological radiation hydrodynamic Technicolor Dawn simulation. In neutral and partially ionized regions where the metagalactic ultraviolet background is weak, the ionization timescale t ion ≡ Γ −1 exceeds the Hubble time. Assuming photoionization equilibrium in such regions artificially boosts the ionization rate, accelerating reionization. By contrast, the recombination time t rec < t ion in photoionized regions, with the result that assuming photoionization equilibrium artificially increases the neutral hydrogen fraction. Using snapshots in the rangemore »8 ≥ z ≥ 5, we compare the predicted Ly α forest (LAF) flux power spectrum with and without the assumption of ionization equilibrium. Small scales ( k > 0.1 rad s km −1 ) exhibit reduced power from 7 ≤ z ≤ 5.5 in the ionization equilibrium case, while larger scales are unaffected. This occurs for the same reasons: ionization equilibrium artificially suppresses the neutral fraction in self-shielded gas and boosts ionizations in voids, suppressing small-scale fluctuations in the ionization field. When the volume-averaged neutral fraction drops below 10 −4 , the signature of nonequilibrium ionizations on the LAF disappears. Comparing with recent observations indicates that these nonequilibrium effects are not yet observable in the LAF flux power spectrum.« less
  4. 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, Fmore »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.« less
  5. 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 thatmore »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.

    « less