We decompose the Lyman- (Ly ) forest of an extensive sample of 75 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in ...the peculiar velocity field have the smallest Doppler parameters, resulting in a low-b cutoff in the b- distribution set primarily by the thermal state of the intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range , which allows us to measure the evolution of the IGM temperature-density ( ) relation parameters T0 and γ. We calibrate our measurements against mock Ly forest data generated using 26 hydrodynamic simulations with different thermal histories from the THERMAL suite, also encompassing different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of T0 (γ) shows a suggestive peak (dip) at z = 2.9 (z = 3). Our measured evolution of T0 and γ is generally in good agreement with previous determinations in the literature. Both the peak in the evolution of T0 at z = 2.8, as well as the high temperatures that we observe at 2.4 < z < 3.4, strongly suggest that a significant episode of heating occurred after the end of H i reionization, which was most likely the cosmic reionization of He ii.
Abstract
We present two measurements of the temperature–density relationship (TDR) of the intergalactic medium (IGM) in the redshift range 2.55 < z < 2.95 using a sample of 13 high-quality quasar ...spectra and high resolution numerical simulations of the IGM. Our approach is based on fitting the neutral hydrogen column density $N_{\rm H\, \small {I}}$ and the Doppler parameter b of the absorption lines in the Lyα forest. The first measurement is obtained using a novel Bayesian scheme that takes into account the statistical correlations between the parameters characterizing the lower cut-off of the $b\hbox{--}N_{\rm H\, \small {I}}$ distribution and the power-law parameters T0 and γ describing the TDR. This approach yields T0/103 K = 15.6 ± 4.4 and γ = 1.45 ± 0.17 independent of the assumed pressure smoothing of the small-scale density field. In order to explore the information contained in the overall $b\hbox{--}N_{\rm H\, \small {I}}$ distribution rather than only the lower cut-off, we obtain a second measurement based on a similar Bayesian analysis of the median Doppler parameter for separate column-density ranges of the absorbers. In this case, we obtain T0/103 K = 14.6 ± 3.7 and γ = 1.37 ± 0.17 in good agreement with the first measurement. Our Bayesian analysis reveals strong anticorrelations between the inferred T0 and γ for both methods as well as an anticorrelation of the inferred T0 and the pressure smoothing length for the second method, suggesting that the measurement accuracy can in the latter case be substantially increased if independent constraints on the smoothing are obtained. Our results are in good agreement with other recent measurements of the thermal state of the IGM probing similar (over-)density ranges.
ABSTRACT The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Ly forest is degenerate with the 3D ...structure of the IGM on ∼100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T 104 K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected filtering scale cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Ly flux, Freal, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The Freal power spectrum is accurately described by a simple fitting function with cutoff at λF, allowing us to rigorously quantify the pressure smoothing scale for the first time: we find λF = 79 kpc (comoving) at z = 3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Ly forest absorption in close quasar pairs, recently proposed as a method to measure the pressure smoothing scale. Our results enable one to quantify the pressure smoothing scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude T0 and slope γ of the temperature-density relation should be augmented with a third pressure smoothing scale parameter λF.
Although the baryons in the intergalactic medium (IGM) trace dark matter fluctuations on megaparsec scales, on smaller scales ~100 kpc, fluctuations are suppressed because the finite temperature gas ...is pressure supported against gravity, analogous to the classical Jeans argument. This Jeans filtering scale, which quantifies the small-scale structure of the IGM, has fundamental cosmological implications. First, it provides a thermal record of heat injected by ultraviolet photons during cosmic reionization events, and thus constrains the thermal and reionization history of the universe. Second, the Jeans scale determines the clumpiness of the IGM, a critical ingredient in models of cosmic reionization. Third, it sets the minimum mass scale for gravitational collapse from the IGM, and hence plays a pivotal role in galaxy formation. Unfortunately, it is extremely challenging to measure the Jeans scale via the standard technique of analyzing purely longitudinal Lyalpha forest spectra, because the thermal Doppler broadening of absorption lines along the line-of-sight, is highly degenerate with Jeans smoothing. In this work, we show that the Jeans filtering scale can be directly measured by characterizing the coherence of correlated Lyalpha forest absorption in close quasar pairs, with separations small enough ~100 kpc to resolve it. We present a novel technique for this purpose, based on the probability density function (PDF) of phase angle differences of homologous longitudinal Fourier modes in close quasar pair spectra. A Bayesian formalism is introduced based on the phase angle PDF, and Markov Chain Monte Carlo techniques are used to characterize the precision of a hypothetical Jeans scale measurement, and explore degeneracies with other thermal parameters governing the IGM. A semi-analytical model of the Lyalpha forest is used to generate a large grid (500) of thermal models from a dark matter only simulation. Our full parameter study indicates that a realistic sample of only 20 close quasar pair spectra can pinpoint the Jeans scale to Asymptotically = to5% precision, independent of the amplitude T sub(0) and slope gamma of the temperature-density relation of the IGM T = T sub(0)(rho/rho) super(gamma-1). This exquisite sensitivity arises because even long-wavelength one-dimensional Fourier modes ~10 Mpc, i.e., two orders of magnitude larger than the Jeans scale, are nevertheless dominated by projected small-scale three-dimensional (3D) power. Hence phase angle differences between all modes of quasar pair spectra actually probe the shape of the 3D power spectrum on scales comparable to the pair separation. We show that this new method for measuring the Jeans scale is unbiased and is insensitive to a battery of systematics that typically plague Lyalpha forest measurements, such as continuum fitting errors, imprecise knowledge of the noise level and/or spectral resolution, and metal-line absorption.
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.
Abstract
We present a new measurement of the Ly
α
forest power spectrum at 1.8 <
z
< 3.4 using 74 Keck/HIRES and VLT/UVES high-resolution, high-signal-to-noise-ratio quasar spectra. We developed a ...custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise and corrects for the bias induced by masking missing data, damped Ly
α
absorption systems, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes
, inaccessible to previous SDSS/BOSS analyses. It is well known that these high-
k
modes are highly sensitive to the thermal state of the intergalactic medium, but contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power and find a modest effect on modes with
. As a result, by masking metals and restricting to
, their impact is completely mitigated. We present an end-to-end Bayesian forward-modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Ly
α
forest simulations. These mock spectra are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform Markov chain Monte Carlo fits. Our results agree well with BOSS on scales
, where the measurements overlap. The combination of the percent-level low-
k
precision of BOSS with our 5%–15% high-
k
measurements results in a powerful new data set for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables.
The distribution of diffuse gas in the intergalactic medium (IGM) imprints a series of hydrogen absorption lines on the spectra of distant background quasars known as the Lyman-α forest. Cosmological ...hydrodynamical simulations predict that IGM density fluctuations are suppressed below a characteristic scale where thermal pressure balances gravity. We measured this pressure-smoothing scale by quantifying absorption correlations in a sample of close quasar pairs. We compared our measurements to hydrodynamical simulations, where pressure smoothing is determined by the integrated thermal history of the IGM. Our findings are consistent with standard models for photoionization heating by the ultraviolet radiation backgrounds that reionized the universe.
We present a new measurement of the Ly forest power spectrum at 1.8 < z < 3.4 using 74 Keck/HIRES and VLT/UVES high-resolution, high-signal-to-noise-ratio quasar spectra. We developed a custom ...pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise and corrects for the bias induced by masking missing data, damped Ly absorption systems, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes , inaccessible to previous SDSS/BOSS analyses. It is well known that these high-k modes are highly sensitive to the thermal state of the intergalactic medium, but contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power and find a modest effect on modes with . As a result, by masking metals and restricting to , their impact is completely mitigated. We present an end-to-end Bayesian forward-modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Ly forest simulations. These mock spectra are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform Markov chain Monte Carlo fits. Our results agree well with BOSS on scales , where the measurements overlap. The combination of the percent-level low-k precision of BOSS with our 5%-15% high-k measurements results in a powerful new data set for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables.
We present a new measurement of the Lyα forest power spectrum at 1.8 < z < 3.4 using 74 Keck/HIRES and VLT/UVES high-resolution, high-signal-to-noise-ratio quasar spectra. We developed a custom ...pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise and corrects for the bias induced by masking missing data, damped Lyα absorption systems, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes $k\gt 0.02\,{\rm{s}}\,{\mathrm{km}}^{-1}$, inaccessible to previous SDSS/BOSS analyses. It is well known that these high-k modes are highly sensitive to the thermal state of the intergalactic medium, but contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power and find a modest effect on modes with $k\lt 0.1\,{\rm{s}}\,{\mathrm{km}}^{-1}$. As a result, by masking metals and restricting to $k\lt 0.1\,{\rm{s}}\,{\mathrm{km}}^{-1}$, their impact is completely mitigated. We present an end-to-end Bayesian forward-modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Lyα forest simulations. These mock spectra are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform Markov chain Monte Carlo fits. Our results agree well with BOSS on scales $k\lt 0.02\,{\rm{s}}\,{\mathrm{km}}^{-1}$, where the measurements overlap. The combination of the percent-level low-k precision of BOSS with our 5%–15% high-k measurements results in a powerful new data set for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. Lastly, the power spectra and their covariance matrices are provided as electronic tables.