skip to main content

Title: 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.
Authors:
; ;
Award ID(s):
1817256
Publication Date:
NSF-PAR ID:
10289365
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
503
Issue:
2
Page Range or eLocation-ID:
1668 to 1679
ISSN:
0035-8711
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Planck data provide precise constraints on cosmological parameters when assuming the base ΛCDM model, including a 0.17 per cent measurement of the age of the Universe, $t_0=13.797 \pm 0.023\, {\rm Gyr}$. However, the persistence of the ‘Hubble tension’ calls the base ΛCDM model’s completeness into question and has spurred interest in models such as early dark energy (EDE) that modify the assumed expansion history of the Universe. We investigate the effect of EDE on the redshift–time relation z↔t and find that it differs from the base ΛCDM model by at least ${\approx } 4{{\ \rm per\ cent}}$ at all t and z. Asmore »long as EDE remains observationally viable, any inferred t ← z or z ← t quoted to a higher level of precision do not reflect the current status of our understanding of cosmology. This uncertainty has important astrophysical implications: the reionization epoch – 10 > z > 6 – corresponds to disjoint lookback time periods in the base ΛCDM and EDE models, and the EDE value of t0 = 13.25 ± 0.17 Gyr is in tension with published ages of some stars, star clusters, and ultrafaint dwarf galaxies. However, most published stellar ages do not include an uncertainty in accuracy (due to, e.g. uncertain distances and stellar physics) that is estimated to be $\sim 7\!-\!10{{\ \rm per\ cent}}$, potentially reconciling stellar ages with $t_{0,\rm EDE}$. We discuss how the big data era for stars is providing extremely precise ages ($\lt 1{{\ \rm per\ cent}}$) and how improved distances and treatment of stellar physics such as convection could result in ages accurate to $4\!-\!5{{\ \rm per\ cent}}$, comparable to the current accuracy of t↔z. Such precise and accurate stellar ages can provide detailed insight into the high-redshift Universe independent of a cosmological model.« less
  2. 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 wheremore »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.« less
  3. 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 andmore »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.« less
  4. 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
  5. ABSTRACT We present a catalogue of 22 755 objects with slitless, optical, Hubble Space Telescope (HST) spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS). The data cover ∼220 sq. arcmin to 7-orbit (∼10 ks) depth in 20 parallel pointings of the Advanced Camera for Survey’s G800L grism. The fields are located 6 arcmin away from 10 massive galaxy clusters in the HFF and CLASH footprints. 13 of the fields have ancillary HST imaging from these or other programs to facilitate a large number of applications, from studying metal distributions at z ∼ 0.5, to quasars at z ∼ 4, to the star formationmore »histories of hundreds of galaxies in between. The spectroscopic catalogue has a median redshift of 〈z〉 = 0.60 with a median uncertainty of $\Delta z / (1+z)\lesssim 2{{\ \rm per\ cent}}$ at $F814\mathit{ W}\lesssim 23$ AB. Robust continuum detections reach a magnitude fainter. The 5 σ limiting line flux is $f_{\rm lim}\approx 5\times 10^{-17}\rm ~erg~s^{-1}~cm^{-2}$ and half of all sources have 50 per cent of pixels contaminated at ≲1 per cent. All sources have 1D and 2D spectra, line fluxes/uncertainties and identifications, redshift probability distributions, spectral models, and derived narrow-band emission-line maps from the Grism Redshift and Line Analysis tool (grizli). We provide other basic sample characterizations, show data examples, and describe sources and potential investigations of interest. All data and products will be available online along with software to facilitate their use.« less