The Lyman-$\alpha$ (hereafter Ly-$\alpha$) forest is a probe of large-scale matter density fluctuations at high redshift, $z > 2.1$. It consists of HI absorption spectra along individual lines-of-sight. If the line-of-sight density is large enough, 3D maps of HI absorption can be inferred by tomographic reconstruction. In this article, we investigate the Ly-$\alpha$ forest available in the Stripe 82 field ($220\,\mathrm{deg^{2}}$), based on the quasar spectra from SDSS Data Release DR16. The density of observed quasar spectra is $37\,\mathrm{deg^{-2}}$ with a mean pixel signal-to-noise ratio of two per angstrom. This study provides an intermediate case between the average SDSS density and that of the much denser but smaller CLAMATO survey. We derive a 3D map of large-scale matter fluctuations from these data, using a Wiener filter technique. The total volume of the map is $0.94\,\mathrm{h^{-3} Gpc^{3}}$. Its resolution is $13\,\mathrm{h^{-1} Mpc}$, which is related to the mean transverse distance between nearest lines-of-sight. From this map, we provide a catalog of voids and protocluster candidates in the cosmic web. The map-making and void catalog are compared to simulated eBOSS Stripe 82 observations. A stack over quasar positions provides a visualization of the Ly-$\alpha$ quasar cross-correlation. This tomographic reconstruction constitutes the largest-volume high-redshift 3D map of matter fluctuations.

The Ly-$\alpha$ forest 1D flux power spectrum is a powerful probe of several cosmological parameters. Assuming a $\Lambda$CDM cosmology including massive neutrinos, we find that the latest SDSS DR14 BOSS and eBOSS Ly-$\alpha$ forest data is in very good agreement with current weak lensing constraints on $(\Omega_m, \sigma_8)$ and has the same small level of tension with Planck. We did not identify a systematic effect in the data analysis that could explain this small tension, but we show that it can be reduced in extended cosmological models where the spectral index is not the same on the very different times and scales probed by CMB and Ly-$\alpha$ data. A particular case is that of a $\Lambda$CDM model including a running of the spectral index on top of massive neutrinos. With combined Ly-$\alpha$ and Planck data, we find a slight (3$\sigma$) preference for negative running, $\alpha_s= -0.010 \pm 0.004$ (68% CL). Neutrino mass bounds are found to be robust against different assumptions. In the $\Lambda$CDM model with running, we find $\sum m_\nu <0.11$ eV at the 95% confidence level for combined Ly-$\alpha$ and Planck (temperature and polarisation) data, or $\sum m_\nu < 0.09$ eV when adding CMB lensing and BAO data. We further provide strong and nearly model-independent bounds on the mass of thermal warm dark matter. For a conservative configuration consisting of SDSS data restricted to $z<4.5$ combined with XQ-100 \lya data, we find $m_X > 5.3\;\mathrm{keV}$ (95\%CL).

The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparable to the Hubble time, memory of reionization heating is retained. Therefore, a late ending reionization has significant implications for the structure of the $z\sim5-6$ Lyman-$\alpha$ (ly$\alpha$) forest. Using state-of-the-art hydrodynamical simulations that allow us to vary the timing of reionization and its associated heat injection, we argue that extant thermal signatures from reionization can be detected via the Ly$\alpha$ forest power spectrum at $5< z<6$. This arises because the small-scale cutoff in the power depends not only the the IGM temperature at these epochs, but is also particularly sensitive to the pressure smoothing scale set by the IGM full thermal history. Comparing our different reionization models with existing measurements of the Ly$\alpha$ forest flux power spectrum at $z=5.0-5.4$, we find that models satisfying Planck's $\tau_e$, constraint favor a moderate amount of heat injection consistent with galaxies driving reionization, but disfavoring quasar-driven scenarios. We explore the impact of different reionization histories and heating models on the shape of the power spectrum, and find that they can produce similar effects, but argue that this degeneracy can be broken with high enough quality data. We study the feasibility of measuring the flux power spectrum at $z\simeq 6$ using mock quasar spectra and conclude that a sample of $\sim10$ high-resolution spectra with an attainable signal-to-noise ratio will allow distinguishing between different reionization scenarios.

We present new measurements of the flux power-spectrum P(k) of the $z<0.5$ HI Lyman-$\alpha$ forest spanning scales k ~ 0.001-0.1 s/km. These results were derived from 65 far ultraviolet quasar spectra (resolution R~18000) observed with the Cosmic Origin Spectrograph (COS) on board the Hubble Space Telescope. The analysis required careful masking of all contaminating, coincident absorption from HI and metal-line transitions of the Galactic interstellar medium and intervening absorbers as well as proper treatment of the complex COS line-spread function. From the P(k) measurements, we estimate the HI photoionization rate ($\Gamma_{\rm HI}$) in the z<0.5 intergalactic medium. Our results confirm most of the previous $\Gamma_{\rm HI}$ estimates. We conclude that previous concerns of a photon underproduction crisis are now resolved by demonstrating that the measured $\Gamma_{\rm HI}$ can be accounted for by ultraviolet emission from quasars alone. In a companion paper, we will present constraints on the thermal state of the $z<0.5$ intergalactic medium from the P(k) measurements presented here.

We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$\alpha$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0 (\rho/\rho_0)^{\gamma-1}$, and previous work has inferred $T_0$ and $\gamma$ by fitting its low-$b$ cutoff. This approach discards the majority of available data, and is susceptible to systematics related to cutoff determination. We present a method that exploits all information encoded in the $b$-$N_{\text{HI}}$ distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, then use principal component decomposition to create an emulator which can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D $(T_0, \gamma)$ measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at $z=2$. Using 200 absorbers, equivalent in pathlength to a single Ly$\alpha$ forest spectrum, we measure $\log T_0 =4.092^{+0.050}_{-0.055}$ and $\gamma=1.49^{+0.073}_{-0.074}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of $\log T_0$ and $\gamma$ with precision on $\log T_0$ ($\gamma$) nearly two (three) times higher for current dataset sizes.

We determine the thermal evolution of the intergalactic medium (IGM) over $3\, \mathrm{Gyr}$ of cosmic time $1.8<z<5.4$ by comparing measurements of the Ly{\alpha} forest power spectrum to a suite of $\sim70$ hydrodynamical simulations. We conduct Bayesian inference of IGM thermal parameters using an end-to-end forward modeling framework whereby mock spectra generated from our simulation grid are used to build a custom emulator which interpolates the power spectrum between thermal grid points. The temperature at mean density $T_0$ rises steadily from $T_0\sim 6000\, \mathrm{K}$ at $z=5.4$, peaks at $14000\, \mathrm{K}$ for $z\sim 3.4$, and decreases at lower redshift reaching $T_0\sim 7000\, \mathrm{K}$ by $z\sim1.8$. This evolution provides conclusive evidence for photoionization heating resulting from the reionization of He II, as well as the subsequent cooling of the IGM due to the expansion of the Universe after all reionization events are complete. Our results are broadly consistent with previous measurements of thermal evolution based on a variety of approaches, but the sensitivity of the power spectrum, the combination of high precision BOSS measurements of large-scale modes ($k\lesssim 0.02\, \mathrm{s/km}$) with our recent determination of the small-scale power, our large grid of models, and our careful statistical analysis allow us to break the well known degeneracy between the temperature at mean density $T_0$ and the slope of the temperature density relation $\gamma$ that has plagued previous analyses. At the highest redshifts $z\geq5$ we infer lower temperatures than expected from the standard picture of IGM thermal evolution leaving little room for additional smoothing of the Ly{\alpha} forest by free streaming of warm dark matter.

We decompose the Lyman-{\alpha} (Ly{\alpha}) forest of an extensive sample of 74 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$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range $2.0\leq z \leq 3.4$, which allows us to measure the evolution of the IGM temperature-density ($T= T_0 (\rho/ \rho_0)^{\gamma-1}$) relation parameters $T_0$ and $\gamma$. We calibrate our measurements against Ly$\alpha$ forest simulations, using 21 different thermal models of the IGM at each redshift, also allowing for 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 $T_0$ shows a suggestive peak at $z=2.8$, while our evolution of $\gamma$ is consistent with $\gamma\simeq 1.4$ and disfavors inverted temperature-density relations. Our measured evolution of $T_0$ and $\gamma$ are generally in good agreement with previous determinations in the literature. Both the peak in the evolution of $T_0$ at $z = 2.8$, as well as the high temperatures $T_0\simeq 15000-20000\,$K that we observe at $2.4 < z < 3.4$, strongly suggest that a significant episode of heating occurred after the end of HI reionization, which was most likely the cosmic reionization of HeII.

We present a new measurement of the Ly${\alpha}$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N 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, DLAs, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes $k > 0.02\,\mathrm{s\,km^{-1}}$, unaccessible 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, however 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 < 0.1\,\mathrm{s\,km^{-1}} $. As a result, by masking metals and restricting to $k < 0.1\,\mathrm{s\,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${\alpha}$ forest simulations. These mocks are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform MCMC fits. Our results agree well with BOSS on scales $k < 0.02\,\mathrm{s\,km^{-1}}$ where the measurements overlap. The combination of BOSS' percent level low-$k$ precision with our $5-15\%$ high-$k$ measurements, results in a powerful new dataset 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-$\alpha$ 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.