-
Damping Wing-Like Features in the Stacked Ly$α$ Forest: Potential Neutral Hydrogen Islands at $z<6$
Authors:
Yongda Zhu,
George D. Becker,
Sarah E. I. Bosman,
Christopher Cain,
Laura C. Keating,
Fahad Nasir,
Valentina D'Odorico,
Eduardo Bañados,
Fuyan Bian,
Manuela Bischetti,
James S. Bolton,
Huanqing Chen,
Anson D'Aloisio,
Frederick B. Davies,
Rebecca L. Davies,
Anna-Christina Eilers,
Xiaohui Fan,
Prakash Gaikwad,
Bradley Greig,
Martin G. Haehnelt,
Girish Kulkarni,
Samuel Lai,
Ewald Puchwein,
Yuxiang Qin,
Emma V. Ryan-Weber
, et al. (6 additional authors not shown)
Abstract:
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$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ dark gaps and identify a damping wing-like extended absorption profi…
▽ More
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$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ 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$.
△ Less
Submitted 28 June, 2024; v1 submitted 20 May, 2024;
originally announced May 2024.
-
Damping wings in the Lyman-α forest: a model-independent measurement of the neutral fraction at 5.4<z<6.1
Authors:
Benedetta Spina,
Sarah E. I. Bosman,
Frederick B. Davies,
Prakash Gaikwad,
Yongda Zhu
Abstract:
Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$α$ 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…
▽ More
Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$α$ 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-$α$ emission lines of quasars at $z\geq7$, DWs are challenging to observe in the general Lyman-$α$ 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-$α$ 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σ$ $(2σ)$ 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.
△ Less
Submitted 21 August, 2024; v1 submitted 20 May, 2024;
originally announced May 2024.
-
The thesan project: public data release of radiation-hydrodynamic simulations matching reionization-era JWST observations
Authors:
Enrico Garaldi,
Rahul Kannan,
Aaron Smith,
Josh Borrow,
Mark Vogelsberger,
Rüdiger Pakmor,
Volker Springel,
Lars Hernquist,
Daniela Galárraga-Espinosa,
Jessica Y. -C. Yeh,
Xuejian Shen,
Clara Xu,
Meredith Neyer,
Benedetta Spina,
Mouza Almualla,
Yu Zhao
Abstract:
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…
▽ More
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.
△ Less
Submitted 21 March, 2024; v1 submitted 12 September, 2023;
originally announced September 2023.
-
The HI-halo mass relation at redshift $z \sim 1$ from the Minkowski functionals of 21 cm intensity maps
Authors:
Benedetta Spina,
Cristiano Porciani,
Carlo Schimd
Abstract:
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 isoco…
▽ More
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.
△ Less
Submitted 18 June, 2021; v1 submitted 22 January, 2021;
originally announced January 2021.