Search

We aim to construct a machine-learning approach that allows for a pixel-by-pixel reconstruction of the intergalactic medium (IGM) density field for various warm dark matter (WDM) models using the Lyman-alpha forest. With this regression machinery, we constrain the mass of a potential WDM particle from observed Lyman-alpha sightlines directly from the density field. We design and train a Bayesian neural network on the supervised regression task of recovering the optical depth-weighted density field $\Delta_\tau$ as well as its reconstruction uncertainty from the Lyman-alpha forest flux field. We utilise the Sherwood-Relics simulation suite at $4.1\leq z \leq 5.0$ as the main training and validation dataset. Leveraging the density field recovered by our neural network, we construct an inference pipeline to constrain the WDM particle masses based on the probability distribution function of the density fields. We find that our trained Bayesian neural network can accurately recover within a $1\sigma$ error $\geq 85\%$ of the density field pixels from a validation simulated dataset that encompasses multiple WDM and thermal models of the IGM. When predicting on Lyman-alpha skewers generated using the alternative hydrodynamical code Nyx not included in the training data, we find a $1\sigma$ accuracy rate $\geq 75\%$. We consider 2 samples of observed Lyman-alpha spectra from the UVES and GHOST instruments, at $z=4.4$ and $z=4.9$ respectively and fit the density fields recovered by our Bayesian neural network to constrain WDM masses. We find lower bounds on the WDM particle mass of $m_{\mathrm{WDM}} \geq 3.8$ KeV and $m_{\mathrm{WDM}} \geq 2.2$ KeV at $2\sigma$ confidence, respectively. We are able to match current state-of-the-art WDM particle mass constraints using up to 40 times less observational data than Markov Chain Monte Carlo techniques based on the Lyman-alpha forest power spectrum., Comment: 12 pages, Submitted to A&A

Recent quasar absorption line observations suggest that reionization may end as late as $z \approx 5.3$. As a means to search for large neutral hydrogen islands at $z<6$, we revisit long dark gaps in the Ly$\beta$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$\alpha$ forest corresponding to both edges of these Ly$\beta$ dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is $z=5.8$. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be $\langle x_{\rm HI} \rangle \geq 6.1 \pm 3.9\%$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$., Comment: 8 pages, 5 figures, 1 table; accepted for publication in MNRAS Letters

Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$\alpha$ forest have not yet been able to discern whether reionisation occurred slowly and late, with substantial neutral hydrogen persisting at redshift $\sim 6$, or rapidly and earlier, with the apparent late end driven by the fluctuating UV background. Gunn-Peterson (GP) absorption troughs are solid indicators that reionisation is not complete until $z=5.3$, but whether they contain significantly neutral gas has not yet been proven. We aim to answer this question by directly measuring, for the first time, the neutral hydrogen fraction ($x_\mathrm{HI}$) at the end of the EoR ($5 \lesssim z \lesssim 6$) in high-redshift quasars spectra. For high neutral fractions $x_\mathrm{HI}\gtrsim0.1$, GP troughs exhibit damping wing (DW) absorption extending over $1000$ km s$^{-1}$ beyond the troughs. While conclusively detected in Lyman-$\alpha$ emission lines of quasars at $z\geq7$, DWs are challenging to observe in the general Lyman-$\alpha$ forest due to absorption complexities and small-scale stochastic transmission features. We report the first successful identification of the stochastic DW signal adjacent to GP troughs at redshifts $z=5.6$ through careful stacking of the dark gaps in Lyman-$\alpha$ forest. We use the signal to present a measurement of the corresponding global $x_\mathrm{HI}=0.19\pm0.07$ $(_{-0.16}^{+0.11})$ at $1\sigma$ $(2\sigma)$ at $z=5.6$ and a limit $x_\mathrm{HI}<0.44$ at $z=5.9$. The detection of this signal demonstrates the existence of substantially neutral islands near the conclusion of the EoR, unequivocally signaling a late-and-slow reionization scenario., Comment: 7 pages, 4 figures, published in A&A Letters, Appendix B added to address the comments of the reviewer

Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope, SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the THESAN simulation project at https://thesan-project.com. This project focuses on the $z \geq 5.5$ Universe, combining a radiation-hydrodynamics solver (AREPO-RT), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The THESAN suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope, affirming the accuracy and applicability of THESAN. The examples also serve as prototypes for how to utilise the released dataset to perform comparisons between predictions and observations., Comment: Data and documentation at https://www.thesan-project.com, comments and requests welcome, paper accepted for publication in MNRAS

The mean and the scatter of the HI content of a dark-matter halo as a function of the halo mass are useful statistics that can be used to test models of structure and galaxy formation. We investigate the possibility of constraining this HI-halo mass relation (HIHMR) from intensity maps of the redshifted 21 cm line. In particular, we use the geometry and topology of the brightness-temperature isocontours in a single frequency channel as quantified by the Minkowski functionals. First, we generate mock maps from a large N-body simulation considering the impact of thermal noise and foreground removal. We then use the Fisher information formalism to forecast constraints on a parametric model for the HIHMR. We consider a 20,000 deg$^2$ survey (originally proposed for dark-energy science) conducted with the Square Kilometre Array Phase 1 (SKA-1) MID observatory operating in single-dish mode. For a channel bandwidth of 2 MHz, we show that an integration time of a few$\,\times\,10^4$ s per pointing is sufficient to image the smoothed HI distribution at redshift $z \simeq 1$ and to measure the HIHMR in a nearly optimal way from the Minkowski functionals. Tighter constraints on some of the parameters can be obtained by using also an independent measurement of the mean HI density. Combining the results from different frequency channels provides exquisite constraints on the evolution of the HIHMR, especially in the central frequency range of the data cube., Comment: 13 pages, 9 figures. Section 6 added to address the comments of the reviewer, no changes in the results