skip to main content


Title: The thesan project: Lyman-α emission and transmission during the Epoch of Reionization
ABSTRACT

The visibility of high-redshift Lyman-alpha emitting galaxies (LAEs) provides important constraints on galaxy formation processes and the Epoch of Reionization (EoR). However, predicting realistic and representative statistics for comparison with observations represents a significant challenge in the context of large-volume cosmological simulations. The thesan project offers a unique framework for addressing such limitations by combining state-of-the-art galaxy formation (IllustrisTNG) and dust models with the arepo-rt radiation-magnetohydrodynamics solver. In this initial study, we present Lyman-alpha centric analysis for the flagship simulation that resolves atomic cooling haloes throughout a $(95.5\, \text{cMpc})^3$ region of the Universe. To avoid numerical artefacts, we devise a novel method for accurate frequency-dependent line radiative transfer in the presence of continuous Hubble flow, transferable to broader astrophysical applications as well. Our scalable approach highlights the utility of LAEs and red damping-wing transmission as probes of reionization, which reveal nontrivial trends across different galaxies, sightlines, and frequency bands that can be modelled in the framework of covering fractions. In fact, after accounting for environmental factors influencing large-scale ionized bubble formation such as redshift and UV magnitude, the variation across galaxies and sightlines mainly depends on random processes including peculiar velocities and self-shielded systems that strongly impact unfortunate rays more than others. Throughout the EoR local and cosmological optical depths are often greater than or less than unity such that the exp (− τ) behaviour leads to anisotropic and bimodal transmissivity. Future surveys will benefit by targeting both rare bright objects and Goldilocks zone LAEs to infer the presence of these (un)predictable (dis)advantages.

 
more » « less
Award ID(s):
1814259
NSF-PAR ID:
10364560
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
512
Issue:
3
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 3243-3265
Size(s):
["p. 3243-3265"]
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    The observability of Lyα emitting galaxies (LAEs) during the Epoch of Reionization can provide a sensitive probe of the evolving neutral hydrogen gas distribution, thus setting valuable constraints to distinguish different reionization models. In this study, we utilize the new thesan suite of large-volume ($L_\text{box} = 95.5\, \text{cMpc}$) cosmological radiation-hydrodynamic simulations to directly model the Lyα emission from individual galaxies and the subsequent transmission through the intergalactic medium. thesan combines the arepo-rt radiation-hydrodynamic solver with the IllustrisTNG galaxy formation model and includes high- and medium-resolution simulations designed to investigate the impacts of halo-mass-dependent escape fractions, alternative dark matter models, and numerical convergence. We find important differences in the Lyα transmission based on reionization history, bubble morphology, frequency offset from line centre, and galaxy brightness. For a given global neutral fraction, Lyα transmission reduces when low-mass haloes dominate reionization over high-mass haloes. Furthermore, the variation across sightlines for a single galaxy is greater than the variation across all galaxies. This collectively affects the visibility of LAEs, directly impacting observed Lyα luminosity functions (LFs). We employ Gaussian Process Regression using SWIFTEmulator to rapidly constrain an empirical model for dust escape fractions and emergent spectral-line profiles to match observed LFs. We find that dust strongly impacts the Lyα transmission and covering fractions of MUV ≲ −19 galaxies in $M_\text{vir} \gtrsim 10^{11}\, \text{M}_{\bigodot }$ haloes, such that the dominant mode of removing Lyα photons in non-LAEs changes from low-IGM transmission to high dust absorption around z ∼ 7.

     
    more » « less
  2. ABSTRACT

    We present the spectroscopic confirmation of the brightest known gravitationally lensed Lyman-break galaxy in the Epoch of Reionization (EoR), A1703-zD1, through the detection of [C ii] 158 $\mu$m at a redshift of z = 6.8269 ± 0.0004. This source was selected behind the strong lensing cluster Abell 1703, with an intrinsic luminosity and a very blue Spitzer/Infrared Array Camera (IRAC) [3.6]–[4.5] colour, implying high equivalent width line emission of [O iii] + Hβ. [C ii] is reliably detected at 6.1σ cospatial with the rest-frame ultraviolet (UV) counterpart, showing similar spatial extent. Correcting for the lensing magnification, the [C ii] luminosity in A1703-zD1 is broadly consistent with the local $L_{\rm [C\, {\small II}]}$–star formation rate (SFR) relation. We find a clear velocity gradient of 103 ± 22 km $\rm s^{-1}$ across the source that possibly indicates rotation or an ongoing merger. We furthermore present spectral scans with no detected [C ii] above 4.6σ in two unlensed Lyman-break galaxies in the Extended Groth Strip (EGS)-Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field at z ∼ 6.6–6.9. This is the first time that the Northern Extended Millimeter Array (NOEMA) has been successfully used to observe [C ii] in a ‘normal’ star-forming galaxy at z > 6, and our results demonstrate its capability to complement the Atacama Large Millimeter/submillimeter Array (ALMA) in confirming galaxies in the EoR.

     
    more » « less
  3. Abstract A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work, we calculate the ionizing escape fractions (fesc) of reionization-era galaxies from the state-of-the-art thesan simulations, which combine an accurate radiation-hydrodynamic solver (arepo-rt) with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($L_\text{box} = 95.5\, \text{cMpc}$). This allows the formation of numerous massive haloes ($M_\text{halo} \gtrsim 10^{10}\, {\text{M}_{\odot }}$), which are often statistically underrepresented in previous studies but are believed to be important to achieving rapid reionization. We find that low-mass galaxies ($M_\text{stars} \lesssim 10^7\, {\text{M}_{\odot }}$) are the main drivers of reionization above z ≳ 7, while high-mass galaxies ($M_\text{stars} \gtrsim 10^8\, {\text{M}_{\odot }}$) dominate the escaped ionizing photon budget at lower redshifts. We find a strong dependence of fesc on the effective star formation rate (SFR) surface density defined as the SFR per gas mass per escape area, i.e. $\bar{\Sigma }_\text{SFR} = \text{SFR}/M_\text{gas}/R_{200}^2$. The variation in halo escape fractions decreases for higher mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. Dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global fesc depends on the dust model. Finally, active galactic nuclei are unimportant for reionization in thesan and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells. 
    more » « less
  4. ABSTRACT

    The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extragalactic astronomy. We investigate the IGM properties and their connection to galaxies at z ≥ 5.5 under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel thesan radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic ‘late’ reionization histories, match available constraints on global IGM properties, and reproduce the recently observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-α transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterizing the spectral shape and separation of Lyman-α transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, thesan reproduces the modulation of the Lyman-α flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, thesan produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.

     
    more » « less
  5. ABSTRACT

    Several studies have detected Lyman-alpha (Ly α) from bright ($M_{\small UV}\lesssim -21.5$) galaxies during the early stages of reionization despite the significantly neutral intergalactic medium. To explain these detections, it has been suggested that z > 7 Ly α emitters (LAEs) inhabit physical Mpc (pMpc)-scale ionized regions powered by overdensities of faint galaxies; however, systematic searches for these overdensities near LAEs have been challenging. Here, we use Cosmic Evolution Early Release Science JWST/Near Infrared Camera imaging to search for large-scale galaxy overdensities near two very ultraviolet (UV)-bright, z = 8.7 LAEs in the Extended Groth Strip (EGS) field. We colour select 27 z = 8.4–9.1 candidates, including the one LAE in the footprint (EGSY8p7). From spectral energy distribution models, we infer moderately faint UV luminosities ($-21.2\lesssim {M_{\small UV}}\lesssim -19.1$) and stellar masses of M* ≈ 107.5–8.8 M⊙. All are efficient ionizing agents ($\xi _{\text{ion}}^{*}\approx 10^{25.5-26.0}$ Hz erg−1) and are generally morphologically simple with only one compact (re ≲ 140 to ∼650 pc) star-forming component. 13 candidates lie within 5 arcmin of EGSY8p7, leading to a factor-of-four galaxy overdensity at ≲5 arcmin (∼1.4 projected pMpc at z ∼ 8.7) separations from EGSY8p7. Separations of 10–15 arcmin (∼2.7–4.1 projected pMpc) are consistent with an average field. The spatial distribution of our sample may qualitatively suggest an R ≥ 2 pMpc ionized bubble encompassing both LAEs in EGS, which is theoretically unexpected but may be possible for a galaxy population four times more numerous than the average to create with moderate escape fractions (fesc ≳ 0.15) over long times (≳ 200 Myr). Upcoming spectroscopic follow-up will characterize the size of any ionized bubble that may exist and the properties of the galaxies powering such a bubble.

     
    more » « less