No Arabic abstract
We present a new investigation of the intergalactic medium (IGM) near the end of reionization using dark gaps in the Lyman-alpha (Ly$alpha$) forest. Using spectra of 55 QSOs at $z_{rm em}>5.5$, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly$alpha$ forest where the transmission averaged over 1 comoving $h^{-1},{rm Mpc}$ bins falls below 5%. Nine ultra-long ($L > 80~h^{-1},{rm Mpc}$) dark gaps are identified at $z<6$. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than $30~h^{-1},{rm Mpc}$, $F_{30}$, as a function of redshift. We measure $F_{30} simeq 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 simeq 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 lesssim 6$ are also potentially consistent with the data. Overall, our results suggest that signature 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 simeq 5.3$.
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 $zsim5-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 Plancks $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 $zsimeq 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 compare a sample of five high-resolution, high S/N Ly$alpha$ forest spectra of bright $6<z lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics 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 characterise 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 prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 pm 1600,10500pm 2100,12000 pm 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 $zsim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.
The transmission of Lyman-{alpha} (Ly{alpha}) in the spectra of distant quasars depends on the density, temperature, and ionization state of the intergalactic medium (IGM). Therefore, high-redshift (z > 5) Ly{alpha} forests could be invaluable in studying the late stages of the epoch of reionization (EoR), as well as properties of the sources that drive it. Indeed, high-quality quasar spectra have now firmly established the existence of large-scale opacity fluctuations at z > 5, whose physical origins are still debated. Here we introduce a Bayesian framework capable of constraining the EoR and galaxy properties by forward-modelling the high-z Ly{alpha} forest. Using priors from galaxy and CMB observations, we demonstrate that the final overlap stages of the EoR (when >95% of the volume was ionized) should occur at z < 5.6, in order to reproduce the large-scale opacity fluctuations seen in forest spectra. However, it is the combination of patchy reionization and the inhomogeneous UV background that produces the longest Gunn-Peterson troughs. Ly{alpha} forest observations tighten existing constraints on the characteristic ionizing escape fraction of galaxies, with the combined observations suggesting f_{rm esc} approx 7^4_3%, and disfavoring a strong evolution with the galaxys halo (or stellar) mass.
Conventional wisdom was that thermal relics from the epoch of reionization (EOR) would vanish swiftly. Recently, however, it was shown that these relics can survive to lower redshifts ($z sim 2$) than previously thought, due to gas at mean density being heated to $T sim 3 times 10^4$ K by reionization, which is inhomogeneous, and shocks. Given the high sensitivities of upcoming Ly$alpha$ forest surveys, this effect will be a novel broadband systematic for cosmological application. From the astrophysical point of view, however, the imprint of inhomogeneous reionization can shed light on the EOR and cosmic dawn. We utilize a hybrid method -- which includes two different simulation codes capable of handling the huge dynamical range -- to show the impact of patchy reionization on the Ly$alpha$ forest and its dependence on different astrophysical scenarios. We found statistically significant deviations in the 1D Ly$alpha$ power spectrum at $k = 0.14$ cMpc$^{-1}$ that range from $sim 1%$ at $z = 2$ up to almost $sim 20%$ at $z = 4$. The deviations in the 3D Ly$alpha$ power spectrum, at the same wavenumber, are large and range from a few per cent at $z = 2$ up to $sim 50%$ at $z = 4$, although these deviations ignore the effect of He II reionization and AGN feedback at $z<4$. By exploiting different $k$-dependence of power spectrum among various astrophysical scenarios, the effect of patchy reionization on the Ly$alpha$ forest power spectrum can open a new window into cosmic reionization and possibly cosmic dawn.
We use a set of AMR hydrodynamic simulations post-processed with the radiative-transfer code RADAMESH to study how inhomogeneous HeII reionization affects the intergalactic medium (IGM). We propagate radiation from active galactic nuclei (AGNs) considering two scenarios for the time evolution of the ionizing sources. We find that HeII reionization takes place in a very inhomogeneous fashion, through the production of well separated bubbles of the ionized phase that subsequently percolate. Overall, the reionization process is extended in time and lasts for a redshift interval Delta z>1. At fixed gas density, the temperature distribution is bimodal during the early phases of HeII reionization and cannot be described by a simple effective equation of state. When HeII reionization is complete, the IGM is characterized by a polytropic equation of state with index gamma~1.20. This relation is appreciably flatter than at the onset of the reionization process (gamma=1.56) and also presents a much wider dispersion around the mean. We extract HI and HeII Ly-alpha absorption spectra from the simulations and fit Voigt profiles to them. We find that the regions where helium is doubly ionized are characterized by different probability density functions of the curvature and of the Doppler b parameters of the HI Ly-alpha forest as a consequence of the bimodal temperature distribution during the early phases of HeII reionization. The column-density ratio in HeII and HI shows a strong spatial variability. Its probability density function rapidly evolves with time reflecting the increasing volume fraction in which ionizing radiation is harder due to the AGN contribution. Finally we show that the number density of the flux-transmission windows per unit redshift and the mean size of the dark gaps in the HeII spectra have the potential to distinguish between different reionization scenarios. (abridged)