skip to main content


Title: Constraining reionization in progress at z  = 5.7 with Lyman-α emitters: voids, peaks, and cosmic variance
ABSTRACT A number of independent observations suggest that the intergalactic medium was significantly neutral at z = 7 and that reionization was, perhaps, still in progress at z = 5.7. The narrowband survey, SILVERRUSH, has mapped over 2000 Lyman-α emitters (LAEs) at these redshifts ( G58). Previous analyses have assumed that reionization was over by z = 5.7, but this data may actually sample the final stages of reionization when the last neutral islands were relegated to the cosmic voids. Motivated by these developments, we re-examine LAE void and peak statistics and their ability to constrain reionization. We construct models of the LAE distribution in (1 Gpc h−1)3 volumes, spanning a range of neutral fractions at z = 5.7 and 6.6. Models with a higher neutral fraction show an enhanced probability of finding holes in the LAE distribution. When comparing models at fixed mean surface density, however, LAEs obscured by neutral gas in the voids must be compensated by visible LAEs elsewhere. Hence, in these models, the likelihood of finding an overdense peak is also enhanced in the latter half of reionization. Compared to the widely used angular two-point correlation function (2PCF), we find that the void probability function (VPF) provides a more sensitive test of models during the latter half of reionization. By comparison, at neutral fractions $\sim 50{{\ \rm per\ cent}}$, the VPF and a simple peak thresholding statistic are both similar to the 2PCF in constraining power. Lastly, we find that the cosmic variance and large-scale asymmetries observed in the SILVERRUSH fields are consistent with large-scale structure in a ΛCDM universe.  more » « less
Award ID(s):
2007499
NSF-PAR ID:
10281495
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
501
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5294 to 5308
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. ABSTRACT

    The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here, we introduce an emissivity framework based on Lyα emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between production and escape of Lyα and LyC photons. Using constraints on the escape fractions of bright LAEs (LLyα > 0.2L*) at z ≈ 2 obtained from resolved Lyα profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at z ≈ 2–6, which is non-trivial given the strong evolution in both the star formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between z ≈ 6−9 under plausible assumptions. Within this framework, the >10 × rise in the UV population-averaged fesc between z ≈ 3–7 naturally arises due to the same phenomena that drive the growing LAE fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm {esc}}=50{{\ \rm per\ cent}}$, ξion = 1025.9 Hz erg−1 for half the bright LAEs), a highly ionizing minority of galaxies with MUV < −17 accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a ‘disco’ ionizing background. We conclude proposing tests to further develop our suggested Lyα-anchored formalism.

     
    more » « less
  3. Abstract

    We probe what constraints for the global ionized hydrogen fraction the void probability function (VPF) clustering can give for the Lyman Alpha Galaxies in the Epoch of Reionization (LAGER) narrowband survey as a function of area. Neutral hydrogen acts like a fog for Lyαemission, and measuring the drop in the luminosity function of Lyαemitters (LAEs) has been used to constrain the ionization fraction in narrowband surveys. However, the clustering of LAEs is independent of the luminosity function’s inherent evolution, and can offer additional constraints for reionization under different models. The VPF measures how likely a given circle is to be empty. It is a volume-averaged clustering statistic that traces the behavior of higher-order correlations, and its simplicity offers helpful frameworks for planning surveys. Using the Jensen et al. simulations of LAEs within various amounts of ionized intergalactic medium, we predict the behavior of the VPF in one (301 × 150.5 × 30 Mpc3), four (5.44 × 106Mpc3), or eight (1.1 × 107Mpc3) fields of LAGER imaging. We examine the VPF at 5′ and 13′, corresponding to the minimum scale implied by the LAE density and the separation of the 2D VPF from random, and the maximum scale from the eight-field 15.5 deg2LAGER area. We find that even a single DECam field of LAGER (2–3 deg2) could discriminate between mostly neutral versus ionized. Additionally, we find four fields allow for the distinction between 30%, 50%, and 95% ionized, and eight fields could even distinguish between 30%, 50%, 73%, and 95% ionized.

     
    more » « less
  4. 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
  5. Abstract The observed large-scale scatter in Ly α opacity of the intergalactic medium at z < 6 implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly α opacity and density. We present selections of z = 5.7 Ly α -emitting galaxies (LAEs) in the fields surrounding two highly opaque quasar sightlines with long Ly α troughs. The fields lie toward the z = 6.0 quasar ULAS J0148+0600, for which we reanalyze previously published results using improved photometric selection, and toward the z = 6.15 quasar SDSS J1250+3130, for which results are presented here for the first time. In both fields, we report a deficit of LAEs within 20 h −1 Mpc of the quasar. The association of highly opaque sightlines with galaxy underdensities in these two fields is consistent with models in which the scatter in Ly α opacity is driven by large-scale fluctuations in the ionizing UV background or by an ultra-late reionization that has not yet concluded at z = 5.7. 
    more » « less