Lyα Forest Research Articles

Ionising properties of galaxies in JADES for a stellar mass complete sample: resolving the cosmic ionising photon budget crisis at the Epoch of Reionisation

Abstract We use NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) to study the ionising properties of a sample of 14652 galaxies at 3 ≤ zphot ≤ 9, 90% complete in stellar mass down to log(M⋆/[M⊙]) ≈ 7.5. Out of the full sample, 1620 of the galaxies have spectroscopic redshift measurements from the literature. We use the spectral energy distribution fitting code Prospector to fit all available photometry and infer galaxy properties. We find a significantly milder evolution of the ionising photon production efficiency (ξion) with redshift and UV magnitude than previously reported. Interestingly, we observe two distinct populations in ξion, distinguished by their burstiness (given by SFR10/SFR100). Both populations show the same evolution with z and MUV, but have a different ξion normalisation. We convolve the more representative log (ξion(z, MUV)) relations (accounting for ∼97% of the sample), with luminosity functions from literature, to place constraints on the cosmic ionising photon budget. By combining our results, we find that one of our models can match the observational constraints from the Lyα forest at z ≲ 6. We conclude that galaxies with MUV between −16 and −20, adopting a reasonable escape fraction, can produce enough ionising photons to ionise the Universe, without exceeding the required ionising photon budget.

Measurements of the thermal and ionization state of the intergalactic medium during the cosmic afternoon

Abstract We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 Lyα absorption lines fitted from 12 archival HST STIS quasar spectra. We employ the machine-learning-based inference method that uses joint Doppler parameter - column density ($b-N_{{\rm {H\,{\small I}}}{}}$) distributions obtained from Lyα forest decomposition. Our results show that the H i photoionization rates, $\Gamma _{{\rm{H\,{\small I}}}{}}$, agree with recent UV background synthesis models, with $\log (\Gamma _{{\rm{H\,{\small I}}}}/\text{s}^{-1})={-11.79}^{+0.18}_{-0.15}$, ${-11.98}^{+0.09}_{-0.09}$, and ${-12.32}^{+0.10}_{-0.12}$ at z = 1.4, 1.2, and 1 respectively. We obtain the IGM temperature at the mean density, T0, and the adiabatic index, γ, as [log (T0/K), γ] = [${4.13}^{+0.12}_{-0.10}$, ${1.34}^{+0.10}_{-0.15}$], $[{3.79}^{+0.11}_{-0.11}$, ${1.70}^{+0.09}_{-0.09}]$ and $[{4.12}^{+0.15}_{-0.25}$, ${1.34}^{+0.21}_{-0.26}]$ at z = 1.4, 1.2 and 1. Our measurements of T0 at z = 1.4 and 1.2 are consistent with the trend predicted from previous z < 3 temperature measurements and theoretical expectations, where the IGM cools down after He ii reionization in the absence of any non-standard heating. However, our T0 measurement at z = 1 shows unexpectedly high IGM temperature. Given the relatively large uncertainty in these measurements, where $\sigma _{T_0} \sim 5000$ K, mostly emanating from the limited size of our dataset, we can not conclude whether the IGM cools down as expected. Lastly, we generate mock datasets to test the constraining power of future measurement with larger datasets. The results demonstrate that, with redshift pathlength Δz ∼ 2 for each redshift bin, three times the current dataset, we can constrain the T0 of IGM within 1500K, which would be sufficient to constrain the IGM thermal history at z < 1.5 conclusively.

High-redshift, Small-scale Tests of Ultralight Axion Dark Matter Using Hubble and Webb Galaxy UV Luminosities

Abstract We calculate the abundance of UV-bright galaxies in the presence of ultralight axion (ULA) dark matter (DM), finding that axions suppress their formation with a non-trivial dependence on redshift and luminosity. We set limits on axion DM using UV luminosity function (UVLF) data, excluding a single axion as all the DM for m ax < 10−21.6 eV and limiting axions with − 26 ≤ log ( m ax / eV ) ≤ − 23 to be less than 22% of the DM (both at 95% credibility). These limits use UVLF measurements from 24,000 sources from the Hubble Space Telescope (HST) at redshifts 4 ≤ z ≤ 10. We marginalize over a parametric model connecting halo mass and UV luminosity. Our results bridge a window in axion mass and fraction previously unconstrained by cosmological data, between large-scale cosmic microwave background and galaxy clustering and the small-scale Lyα forest. These high-z measurements provide a powerful consistency check of low-z tests of axion DM, including the recent hint for a sub-dominant ULA DM fraction in Lyα forest data. We also consider a sample of 25 spectroscopically confirmed high-z galaxies from the James Webb Space Telescope (JWST), finding these data to be consistent with HST. Combining HST and JWST UVLF data does not improve our constraints beyond HST alone, but future JWST measurements have the potential to improve these results. We also find an excess of low-mass halos (<109 M ⊙) at z < 3, which could be probed by subgalactic structure probes (e.g., stellar streams, satellite galaxies, and strong lensing).

The COS CGM Compendium. V. The Dichotomy of O vi Associated with Low- and High-metallicity Cool Gas at z < 1

We analyze the O vi content and kinematics for 126 H i-selected absorbers at 0.14 ≲ z ≲ 0.73 for which the metallicities of their cool photoionized phase have been determined. We separate the absorbers into 100 strong Lyα forest systems (SLFSs with 15 ≲ log N(H i) < 16.2) and 26 partial Lyman Limit systems (pLLSs with 16.2 ≤ log N(H i) ≤ 17.2). The sample is drawn from the COS CGM Compendium (CCC) and has O vi coverage in signal-to-noise ratio ≥ 8 Hubble Space Telescope/COS G130M/G160M QSO spectra, yielding a 2σ completeness level of logN (O vi) ≥ 13.6. The O vi detection rates differ substantially between low-metallicity (LM; [X/H] ≤ −1.4) and high-metallicity (HM; [X/H] > −1.4) SLFSs, with 20% and 60% detection rates, respectively. The O vi detection frequency for the HM and LM pLLSs is, however, similar at ∼60%. The SLFSs and pLLSs without detected O vi are consistent with the absorbing gas being in a single phase, while those with O vi trace multiphase gas. We show that the O vi velocity widths and column densities have different distributions in LM and HM gas. We find a strong correlation between O vi column density and metallicity. The strongest ( logN (O vi) ≳ 14) and broadest O vi absorbers are nearly always associated with HM absorbers, while weaker O vi absorbers are found in both LM and HM absorbers. From comparisons with galaxy-selected and blind O vi surveys, we conclude absorbers with logN (O vi) ≳ 14 most likely arise in the CGM of star-forming galaxies. Absorbers with weak O vi likely trace the extended CGM or intergalactic medium, while those without O vi likely originate in the intergalactic medium.

Deep learning the intergalactic medium using lyman-alpha forest at 4 ≤ z ≤ 5

Abstract Unveiling the thermal history of the intergalactic medium (IGM) at 4 ≤ z ≤ 5 holds the potential to reveal early onset He ii reionization or lingering thermal fluctuations from H i reionization. We set out to reconstruct the IGM gas properties along simulated Lyman-alpha forest data on pixel-by-pixel basis, employing deep neural networks. Our approach leverages the Sherwood-Relics simulation suite, consisting of diverse thermal histories, to generate mock spectra. Our convolutional and residual networks with likelihood metric predicts the Lyα optical depth-weighted density or temperature for each pixel in the Lyα forest skewer. We find that our network can successfully reproduce IGM conditions with high fidelity across range of instrumental signal-to-noise. These predictions are subsequently translated into the temperature-density plane, facilitating the derivation of reliable constraints on thermal parameters. This allows us to estimate temperature at mean cosmic density, T0 , with one sigma confidence, $\delta {\rm T_{\rm 0}}\,$≲ 1000K, using only one 20h−1cMpc sightline (Δz ≃ 0.04) with a typical reionization history. Existing studies utilize redshift pathlength comparable to Δz ≃ 4 for similar constraints. We can also provide more stringent constraints on the slope (1σ confidence interval, δγ ≲ 0.1) of the IGM temperature-density relation as compared to other traditional approaches. We test the reconstruction on a single high signal-to-noise observed spectrum (20h−1cMpc segment), and recover thermal parameters consistent with current measurements. This machine learning approach has the potential to provide accurate yet robust measurements of IGM thermal history at the redshifts in question.

A measurement of the escaping ionising efficiency of galaxies at redshift 5

The escaping ionising efficiency from galaxies, $f_ esc ion $, is a crucial ingredient for understanding their contribution to hydrogen reionisation, but both of its components, $f_ esc $ and $ ion $, are extremely difficult to measure. We measured the average escaping ionising efficiency $ esc ion of galaxies at $z=5$ implied by the mean level of ionisation in the intergalactic medium via the Lyman-alpha forest. We used the fact that $ N ion UV f_ esc ion $, the product of the ionising output and the UV density $ UV $, can be calculated from the known average strength of the UV background and the mean free path of ionising photons. These quantities, as well as $ UV $, are robustly measured at $z We calculated the missing factor of $ esc ion at $z=5$ during a convenient epoch after hydrogen reionisation had been completed and the intergalactic medium had reached ionisation equilibrium but before bright quasars began to dominate the ionising photon production. Intuitively, our constraint corresponds to the required escaping ionising production from galaxies in order to avoid over- or under-ionising the Lyman-alpha forest. We obtained a measurement of $ esc ion /$erg $ at $z=5$ when integrating the $ UV $ down to a limiting magnitude $M_ lim =-11$. Our measurement of the escaping ionising efficiency of galaxies is in rough agreement with both observations of early galaxies and with most models.

The PAU Survey: The quasar and UV luminosity functions at 2.7&lt;z&lt;5.3

We present the Lyman-alpha ( and ultraviolet (UV) luminosity function (LF) in bins of redshift for quasars selected in the Physics of the Accelerating Universe Survey (PAUS). A sample of 915 objects was selected at $2.7&lt;z&lt;5.3$ within an effective area of $ 36$ deg$^2$ observed in 40 narrow-band (NB) filters (FWHM $ We cover the intermediate--bright luminosity regime of the LF erg\,s $; $-29&lt;M_ UV &lt;-24$). The continuous wavelength coverage of the PAUS NB set allows very efficient target identification and precise redshift measurements. We show that our method is able to retrieve a relatively complete ($C 85&lt;!PCT!&gt;$) and pure ($P 90&lt;!PCT!&gt;$) sample of quasars for $ $. In order to obtain corrections for the LF estimation, and assess the accuracy of our selection method, we produced mock catalogs of $0&lt;z&lt;4.3$ quasars and galaxies that mimic our target population and their main contaminants. Our results show a clear evolution of the and UV LFs, with a declining tendency in the number density of quasars toward increasing redshifts. In addition, the faint-end power-law slope of the LF becomes steeper with redshift, suggesting that the number density of quasars declines faster than that of fainter emitters. By integrating the LF, we find that the total emitted by bright quasars per unit volume rapidly declines with increasing redshift, being subdominant to that of star-forming galaxies by several orders of magnitude by $z 4$. Finally, we stack the NB pseudo-spectra of a visually selected ``golden sample'' of 591 quasars to obtain photometric composite SEDs in bins of redshift, enabling us to measure the mean intergalactic medium absorption using the Lyman-alpha forest as a function of redshift, yielding results consistent with previous spectroscopic determinations.

LYαNNA: A deep learning field-level inference machine for the Lyman-α forest

The inference of astrophysical and cosmological properties from the Lyman-α forest conventionally relies on summary statistics of the transmission field that carry useful but limited information. We present a deep learning framework for inference from the Lyman-α forest at the field level. This framework consists of a 1D residual convolutional neural network (ResNet) that extracts spectral features and performs regression on thermal parameters of the intergalactic medium that characterize the power-law temperature-density relation. We trained this supervised machinery using a large set of mock absorption spectra from NYX hydrodynamic simulations at z = 2.2 with a range of thermal parameter combinations (labels). We employed Bayesian optimization to find an optimal set of hyperparameters for our network, and then employed a committee of 20 neural networks for increased statistical robustness of the network inference. In addition to the parameter point predictions, our machine also provides a self-consistent estimate of their covariance matrix with which we constructed a pipeline for inferring the posterior distribution of the parameters. We compared the results of our framework with the traditional summary based approach, namely the power spectrum and the probability density function (PDF) of transmission, in terms of the area of the 68% credibility regions as our figure of merit (FoM). In our study of the information content of perfect (noise- and systematics-free) Lyα forest spectral datasets, we find a significant tightening of the posterior constraints – factors of 10.92 and 3.30 in FoM over the power spectrum only and jointly with PDF, respectively – which is the consequence of recovering the relevant parts of information that are not carried by the classical summary statistics.

The 3D Lyman-α forest power spectrum from eBOSS DR16

Abstract We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly-α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ∼ 205,000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp (ik · r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artifacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line-of-sight of Ly-α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly-α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range (0.02 ≤ k ≤ 0.35) h Mpc−1 at the 1 − 2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly-α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO and 4MOST.

The Negative Baryon Acoustic Oscillation Shift in the Lyα Forest from Cosmological Simulations

We present the first measurement of the Lyα forest baryon acoustic oscillations (BAO) shift parameter from cosmological simulations. In particular, we generate a suite of 1000 accurate effective field-level bias-based Lyα forest simulations of volume V=(1h−1Gpc)3 at z = 2, both in real and redshift space, calibrated upon two fixed-and-paired cosmological hydrodynamic simulations. To measure the BAO, we stack the 3D power spectra of the 1000 different realizations, compute the average, and use a model accounting for a proper smooth-peak component decomposition of the power spectrum, to fit it via an efficient Markov Chain Monte Carlo scheme estimating the covariance matrices directly from the simulations. We report the BAO shift parameters to be α=0.9969−0.0014+0.0014 and α=0.9905−0.0027+0.0027 in real and redshift space, respectively. We also measure the bias b lya and the BAO broadening parameter Σnl, finding blya=−0.1786−0.0001+0.0001 and Σnl=3.87−0.20+0.20 in real space, and blya=−0.073−0.004+0.005 and Σnl=6.55−0.22+0.23 in redshift space. Moreover, we measure the linear Kaiser factor βlya=1.39−0.18+0.24 from the isotropic redshift space fit. Overall, we find evidence for a negative shift of the BAO peak at the ∼2.2σ and ∼3.5σ levels in real and redshift space, respectively. This work sets new important theoretical constraints on the Lyα forest BAO scale and offers a potential solution to the tension emerging from previous observational analysis, in light of ongoing and upcoming Lyα forest spectroscopic surveys, such as DESI, the Prime Focus Spectrograph Survey, and WEAVE-QSO.

Cosmic Reionization in the JWST Era: Back to AGNs?

Deep surveys with the James Webb Space Telescope (JWST) have revealed an emergent population of moderate-luminosity, broad-line active galactic nuclei (AGNs) at 4 ≲ z ≲ 13 powered by accretion onto early massive black holes. The high number densities reported, together with the large Lyman-continuum (LyC) production efficiency and leakiness into the intergalactic medium that are typical of UV-selected AGNs, lead us to reassess a scenario where AGNs are the sole drivers of the cosmic hydrogen/helium reionization process. Our approach is based on the assumptions, grounded in recent observations, that (a) the fraction of broad-line AGNs among galaxies is around 10%–15%; (b) the mean escape fraction of hydrogen LyC radiation is high, at ≳80%, in AGN hosts and is negligible otherwise; and (c) internal absorption at 4 ryd or a steep ionizing EUV spectrum delay full reionization of He ii until z ≃ 2.8–3.0, in agreement with observations of the He ii Lyα forest. In our fiducial models, (1) hydrogen reionization is 99% completed by redshift z ≃ 5.3–5.5 and reaches its midpoint at z ≃ 6.5–6.7; (2) the integrated Thomson scattering optical depth to reionization is ≃0.05, consistent with constraints from cosmic microwave background anisotropy data; and (3) the abundant AGN population detected by JWST does not violate constraints on the unresolved X-ray background.

TensorFlow Hydrodynamics Analysis for Ly-[formula omitted] Simulations

We introduce the Python program THALAS (TensorFlow Hydrodynamics Analysis for Lyman-Alpha Simulations), which maps baryon fields (baryon density, temperature, and velocity) to Lyα optical depth fields in both real space and redshift space. Unlike previous Lyα codes, THALAS is fully differentiable, enabling a wide variety of potential applications for general analysis of hydrodynamical simulations and cosmological inference. To demonstrate THALAS’s capabilities, we apply it to the Lyα forest inversion problem: given a Lyα optical depth field, we reconstruct the corresponding real-space dark matter density field. Such applications are relevant to both cosmological and three-dimensional tomographic analyses of Lyman Alpha forest data.

Impact of Self-shielding Minihalos on the Lyα Forest at High Redshift

Dense gas in minihalos with masses of 106−108 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Lyα forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Lyα opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10−12 s−1, as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Lyα absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Lyα flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc−1 (or 10−3 km−1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10−12 s−1) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Lyα forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.

Measurement of the small-scale 3D Lyman-α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D, of the Lyα forest at wavenumbers k corresponding to small, ∼ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D. This methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D.

Simulating the epoch of helium reionization in photon-conserving semi-numerical code SCRIPT

The reionization of the second electron of helium (HeII) leaves important imprints on the thermal and ionization state of the intergalactic medium (IGM). Observational evidence suggests that HeII reionization ended at z ≃ 3 due to ionizing photons emitted predominantly by quasars. We present efficient semi-numerical simulations of helium reionization in a 230 h-1 Mpc box, that takes into account the spatial patchiness of reionization coupled with photoheating of the IGM. Dark matter haloes are assigned quasars using empirical measurements of the quasar luminosity function, assuming a universal quasar lifetime consistent with duty cycle values inferred from measurements of the quasar clustering. The ionizing photon field from quasars is then included in the semi-numerical Code for ReionIzation with PhoTon conservation (SCRIPT), which was originally developed for modeling hydrogen reionization. In this work, we make appropriate modifications to SCRIPT for modeling inhomogenous HeII reionization and the corresponding thermal history of the IGM is modelled via a subgrid prescription. Our model has three main free parameters i.e. the global clumping factor 𝒞HeIII, the temperature increase due to photoheating T re He and the quasar spectral energy distribution (SED) index, α UV. Our fiducial model with 𝒞HeIII = 15.6 and T re He ∼ 6000 K gives reasonable values for the empirical measurements of the temperature density equation of state at these redshifts, assuming that quasars brighter than M1450 < -21 and having α UV = 1.7 contribute to HeII reionization. The efficiency of our code shows promising prospects for performing parameter estimation in future, for models of HeII reionization using observations of the Lyα forest.

IGM damping wing constraints on the tail end of reionization from the enlarged XQR-30 sample

ABSTRACT The attenuation of Lyα photons by neutral hydrogen in the intergalactic medium (IGM) at z ≳ 5 continues to be a powerful probe for studying the epoch of reionization. Given a framework to estimate the intrinsic (true) Lyα emission of high-z sources, one can infer the ionization state of the IGM during reionization. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high signal-to-noise quasar spectra between $5.8\lesssim \, z\lesssim \, 6.6$ obtained with VLT/X-shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Lyα emission profiles (obtained via covariance matrix reconstruction of the Lyα and N v emission lines from unattenuated high-ionization emission line profiles; C iv, Si iv + O iv], and C iii]) and (ii) distribution of ionized regions within the IGM using synthetic damping wing profiles drawn from a 1.63 Gpc3 reionization simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionization. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14 -0.12}$ and $0.29\substack{+0.14 -0.13}$ at z = 6.15 and 6.35. Further, we also report 68th percentile upper limits of $\bar{x}_{\mathrm{H\, {\small I}}{}} \lt 0.21$, 0.20, 0.21, and 0.18 at z = 5.8, 5.95, 6.05, and 6.55. These results imply reionization is still ongoing at $5.8\lesssim \, z\lesssim \, 6.55$, consistent with previous results from XQR-30 (dark fraction and Lyα forest) along with other observational probes considered in the literature.

Constraints on the Evolution of the Ionizing Background and Ionizing Photon Mean Free Path at the End of Reionization

The variations in Lyα forest opacity observed at z > 5.3 between lines of sight to different background quasars are too strong to be caused by fluctuations in the density field alone. The leading hypothesis for the cause of this excess variance is a late, ongoing reionization process at redshifts below six. Another model proposes strong ionizing background fluctuations coupled to a short, spatially varying mean free path of ionizing photons, without explicitly invoking incomplete reionization. With recent observations suggesting a short mean free path at z ∼ 6, and a dramatic improvement in z > 5 Lyα forest data quality, we revisit this latter possibility. Here, we apply the likelihood-free inference technique of approximate Bayesian computation (ABC) to jointly constrain the hydrogen photoionization rate ΓHI and the mean free path of ionizing photons λ mfp from the effective optical depth distributions at z = 5.0–6.1 from XQR-30. We find that the observations are well-described by fluctuating mean free path models with average mean free paths that are consistent with the steep trend implied by independent measurements at z ∼ 5–6, with a concomitant rapid evolution of the photoionization rate.

Covariance matrices for the Lyman-α forest using the lognormal approximation

We investigate the nature of correlations in the small-scale flux statistics of the Lyman-α (Lyα) forest across redshift bins. Understanding and characterising these correlations is important for unbiased cosmological and astrophysical parameter inference using the Lyα forest. We focus on the 1-dimensional flux power spectrum (FPS) and mean flux (F̅) simulated using the semi-numerical lognormal model we developed in earlier work. The lognormal model can capture the effects of long wavelength modes with relative ease as compared to full smoothed particle hydrodynamical (SPH) simulations that are limited by box volume. For a single redshift bin of size Δz ≃ 0.1, we show that the lognormal model predicts positive cross-correlations between k-bins in the FPS, and a negative correlation for F̅ × FPS, in qualitative agreement with SPH simulations and theoretical expectations. For measurements across two neighbouring redshift bins of width Δ z each (using long flux skewers of length 2Δ z that are 'split' in half), the lognormal model predicts an anti-correlation for FPS × FPS and a positive correlation for F̅ × FPS, caused by modes with the longest wavelengths. This is in contrast to SPH simulations which predict a negligible magnitude for cross-redshift correlations derived from such `split' skewers, and we discuss possible reasons for this difference. Finally, we perform a preliminary test of the impact of neglecting long wavelength modes on parameter inference, finding that whereas the correlation structure of neighbouring redshift bins has relatively little impact, the absence of long wavelength modes in the model can lead to ≳ 2 - σ biases in the inference of astrophysical parameters. Our results motivate a more careful treatment of long wavelength modes in analyses that rely on the small-scale Lyα forest for parameter inference.

The long-lasting effect of X-ray pre-heating in the post-reionization intergalactic medium

ABSTRACT X-ray photons can penetrate deep into the intergalactic medium (IGM), leading to pre-heating of the IGM prior to cosmic reionization. X-ray pre-heating wipes out some of the small-scale structures that would otherwise be present prior to the passage of an ionization front. Accurate modelling of the small-scale structure is vital to the post-reionization IGM, since the small-scale structure is ultimately the dominant source of long-lasting relics from hydrogen reionization. However, the precise impact of X-ray pre-heating in the fossils from hydrogen reionization is highly uncertain. In this work, we explore and establish for the first time, the long-lasting impact of X-ray pre-heating in the post-reionization IGM via hydrodynamic simulations with high-mass resolution. We find that the addition of X-ray pre-heating astrophysics leads to an overall lesser impact of the effect of inhomogeneous reionization in the Lyα forest – depending on specific X-ray prescription – at low redshifts (z ∼ 2) with respect to a model with no X-ray pre-heating. However, at high redshifts (z ∼ 4), our results indicate a strengthening of the relics of reionization in the Lyα forest because the IGM becomes more transparent compared to the scenario with no pre-heating. Thus, the absence of X-ray pre-heating in Lyα modelling can lead to a biased inference of cosmological parameters. Nevertheless, optimistically, the inclusion of X-ray pre-heating emerges as a promising novel avenue to probe the astrophysics of cosmic dawn.

Large-scale excess H i absorption around z ≈ 4 galaxies detected in a background galaxy spectrum in the MUSE eXtremely deep field

ABSTRACT Observationally mapping the relation between galaxies and the intergalactic medium (IGM) is of key interest for studies of cosmic reionization. Diffuse hydrogen gas has typically been observed in H i Lyman-α (Lyα) absorption in the spectra of bright background quasars. However, it is important to extend these measurements to background galaxies as quasars become increasingly rare at high redshift and rarely probe closely separated sight lines. Here, we use deep integral field spectroscopy in the MUSE eXtremely Deep Field to demonstrate the measurement of the Lyα transmission at z ≈ 4 in absorption to a background galaxy at z = 4.77. The H i transmission is consistent with independent quasar sight lines at similar redshifts. Exploiting the high number of spectroscopic redshifts of faint galaxies (500 between z = 4.0–4.7 within a radius of 8 arcmin) that are tracers of the density field, we show that Lyα transmission is inversely correlated with galaxy density, i.e. transparent regions in the Lyα forest mark underdense regions at z ≈ 4. Due to large-scale clustering, galaxies are surrounded by excess H i absorption over the cosmic mean out to 4 cMpc/h70. We also find that redshifts from the peak of the Lyα line are typically offset from the systemic redshift by +170 km s−1. This work extends results from z ≈ 2–3 to higher redshifts and demonstrates the power of deep integral field spectroscopy to simultaneously measure the ionization structure of the IGM and the large-scale density field in the early Universe.