No Arabic abstract
The lya forest at high redshifts is a powerful probe of reionization. Modeling and observing this imprint comes with significant technical challenges: inhomogeneous reionization must be taken into account while simultaneously being able to resolve the web-like small-scale structure prior to reionization. In this work we quantify the impact of inhomogeneous reionization on the lya forest at lower redshifts ($2 < z < 4$), where upcoming surveys such as DESI will enable precision measurements of the flux power spectrum. We use both small box simulations capable of handling the small-scale structure of the lya forest and semi-numerical large box simulations capable of representing the effects of inhomogeneous reionization. We find that inhomogeneous reionization could produce a measurable effect on the lya forest power spectrum. The deviation in the 3D power spectrum at $z_{rm obs} = 4$ and $k = 0.14 rm{Mpc}^{-1}$ ranges from $19 - 36%$, with a larger effect for later reionization. The corrections decrease to $2.0 - 4.1%$ by $z_{rm obs} = 2$. The impact on the 1D power spectrum is smaller, and ranges from $3.3 - 6.5%$ at $z_{rm obs}=4$ to $0.35 - 0.75%$ at $z_{rm obs}=2$, values which are comparable to the statistical uncertainties in current and upcoming surveys. Furthermore, we study how can this systematic be constrained with the help of the quadrupole of the 21 cm power spectrum.
In La Plante et al. (2017), we presented a new suite of hydrodynamic simulations with the aim of accurately capturing the process of helium II reionization. In this paper, we discuss the observational signatures present in the He II Ly$alpha$ forest. We show that the effective optical depth of the volume $tau_mathrm{eff}$ is not sufficient for capturing the ionization state of helium II, due to the large variance inherent in sightlines. However, the He II flux PDF can be used to determine the timing of helium II reionization. The amplitude of the one-dimensional flux power spectrum can also determine the ionization state of helium II. We show that even given the currently limited number of observations ($sim$50 sightlines), measurements of the flux PDF can yield information about helium II reionization. Further, measurements using the one-dimensional power spectrum can provide clear indications of the timing of reionization, as well as the relative bias of sources of ionizing radiation.
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)
One well-known way to constrain the hydrogen neutral fraction, x_H, of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lya absorption line. We examine this methods effectiveness in light of recent models showing that the IGM neutral fraction is highly inhomogeneous on large scales during reionization. Using both analytic models and semi-numeric simulations, we show that the picket-fence absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of x_H. In particular, we show that simple fits to the damping wing tend to overestimate the true neutral fraction in a partially ionized universe, with a fractional error of ~ 30% near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate x_H if the observational sample only probes a subset of the entire halo population, such as quasars with large HII regions. We also find that the damping wing absorption profile is generally steeper than one would naively expect in a homogeneously ionized universe. The profile steepens and the sightline-to-sightline scatter increases as reionization progresses. Of course, the bias and scatter also depend on x_H and so can, at least in principle, be used to constrain it. Damping wing constraints must therefore be interpreted by comparison to theoretical models of inhomogeneous reionization.
Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${left<n^2right>/left<nright>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.
The angular positions of quasars are deflected by the gravitational lensing effect of foreground matter. The Lyman-alpha forest seen in the spectra of these quasars is therefore also lensed. We propose that the signature of weak gravitational lensing of the forest could be measured using similar techniques that have been applied to the lensed Cosmic Microwave Background, and which have also been proposed for application to spectral data from 21cm radio telescopes. As with 21cm data, the forest has the advantage of spectral information, potentially yielding many lensed slices at different redshifts. We perform an illustrative idealized test, generating a high resolution angular grid of quasars (of order arcminute separation), and lensing the Lyman-alphaforest spectra at redshifts z=2-3 using a foreground density field. We find that standard quadratic estimators can be used to reconstruct images of the foreground mass distribution at z~1. There currently exists a wealth of Lya forest data from quasar and galaxy spectral surveys, with smaller sightline separations expected in the future. Lyman-alpha forest lensing is sensitive to the foreground mass distribution at redshifts intermediate between CMB lensing and galaxy shear, and avoids the difficulties of shape measurement associated with the latter. With further refinement and application of mass reconstruction techniques, weak gravitational lensing of the high redshift Lya forest may become a useful new cosmological probe.