skip to main content


Title: Inferring the intergalactic medium neutral fraction at z ∼ 6–8 with low-luminosity Lyman break galaxies
ABSTRACT

We present a Bayesian inference on the neutral hydrogen fraction of the intergalactic medium (IGM), $\overline{x}_{\small HI}$, at z ∼ 6–8 using the properties of Lyman break galaxies (LBGs) during the epoch of reionization. We use large samples of LBG candidates at 5.5 ≤ z ≤ 8.2 with spectroscopy from Keck/DEIMOS and Keck/MOSFIRE. For each galaxy, we incorporate either the Lyman-α (Lyα) equivalent width (EW) for detections or the EW limit spectrum for non-detections to parametrize the EW distribution at various ultraviolet brightnesses for a given redshift. Using our reference sample of galaxy candidates from the ionized universe at z ∼ 6.0, we are able to infer $\overline{x}_{\small HI}$ at two redshifts: z ∼ 6.7 and z ∼ 7.6. This work includes intrinsically faint, gravitationally lensed galaxies at z ∼ 6.0 in order to constrain the intrinsic faint-end Lyα EW distribution and provide a comparable population of galaxies to counterparts in our sample that are at higher redshift. The inclusion of faint galaxy candidates, in addition to a more sophisticated modelling framework, allows us to better isolate effects of the interstellar medium and circumgalactic medium on the observed Lyα distribution from those of the IGM. We infer an upper limit of $\overline{x}_{\small HI}$ ≤ 0.25 (0.44) at z = 6.7 ± 0.2 and a neutral fraction of $\overline{x}_{\small HI}$ = $0.83^{+0.08}_{-0.11}$ (0.83$^{+0.11}_{-0.21}$) at z = 7.6 ± 0.6, both within 68 per cent (95 per cent) uncertainty, results that favour a moderately late and fairly rapid reionization.

 
more » « less
NSF-PAR ID:
10377235
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
517
Issue:
3
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 3263-3274
Size(s):
["p. 3263-3274"]
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The reionization of hydrogen is closely linked to the first structures in the Universe, so understanding the timeline of reionization promises to shed light on the nature of these early objects. In particular, transmission of Lyman alpha (Ly α) from galaxies through the intergalactic medium (IGM) is sensitive to neutral hydrogen in the IGM, so can be used to probe the reionization timeline. In this work, we implement an improved model of the galaxy UV luminosity to dark matter halo mass relation to infer the volume-averaged fraction of neutral hydrogen in the IGM from Ly α observations. Many models assume that UV-bright galaxies are hosted by massive dark matter haloes in overdense regions of the IGM, so reside in relatively large ionized regions. However, observations and N-body simulations indicate that scatter in the UV luminosity–halo mass relation is expected. Here, we model the scatter (though we assume the IGM topology is unaffected) and assess the impact on Ly α visibility during reionization. We show that UV luminosity–halo mass scatter reduces Ly α visibility compared to models without scatter, and that this is most significant for UV-bright galaxies. We then use our model with scatter to infer the neutral fraction, $\overline{x}_{\mathrm{ H}\,{\small I}}$, at z ∼ 7 using a sample of Lyman-break galaxies in legacy fields. We infer $\overline{x}_{\mathrm{ H}\,{\small I}} = 0.55_{-0.13}^{+0.11}$ with scatter, compared to $\overline{x}_{\mathrm{ H}\,{\small I}} = 0.59_{-0.14}^{+0.12}$ without scatter, a very slight decrease and consistent within the uncertainties. Finally, we place our results in the context of other constraints on the reionization timeline and discuss implications for future high-redshift galaxy studies. 
    more » « less
  2. ABSTRACT

    We describe JWST/NIRSpec prism measurements of Ly α emission in z ≳ 5 galaxies. We identify Ly α detections in 10 out of 69 galaxies with robust rest-optical emission-line redshift measurements at 5 ≤ z < 7 in the Cosmic Evolution Early Release Science (CEERS) and DDT-2750 observations of the Extended Groth Strip field. Galaxies at z ≃ 6 with faint continuum (F150W=27–29 mag) are found with extremely large rest-frame Ly α equivalent widths (EWs; ranging up to 286 Å). Likely Ly α detections are also seen in two new z > 7 galaxies (z = 7.49 and 7.17) from the second epoch of CEERS observations, both showing large Ly α EWs that likely indicate significant transmission through the intergalactic medium (IGM). We measure high Ly α escape fractions in the 12 Ly α emitters in our sample (median 0.28), two of which show $f_{\rm esc}^{ {\rm Ly}\alpha }$ near unity (>0.80). We find that $50_{-11}^{+11}$ per cent of z ≃ 6 galaxies with [O iii] + H β EW>1000 Å have $f_{\rm esc}^{ {\rm Ly}\alpha }$ >0.2, consistent with the fractions found in lower redshift samples with matched [O iii] + H β EWs. While uncertainties are still significant, we find that only $10_{-5}^{+9}$ per cent of z > 7 galaxies with similarly strong rest optical emission lines show such large $f_{\rm esc}^{ {\rm Ly}\alpha }$, as may be expected if IGM attenuation of Ly α increases towards higher redshifts. We identify photometric galaxy overdensities near the z ≳ 7 Ly α emitters, potentially providing the ionizing flux necessary to create large ionized sightlines that facilitate Ly α transmission. Finally, we investigate the absence of Ly α emission in a comparable (and spectroscopically confirmed) galaxy overdensity at z = 7.88 in the Abell 2744 field, discussing new prism spectra of the field obtained with the UNCOVER program.

     
    more » « less
  3. ABSTRACT

    The ionizing photon escape fraction [Lyman continuum (LyC) fesc] of star-forming galaxies is the single greatest unknown in the reionization budget. Stochastic sightline effects prohibit the direct separation of LyC leakers from non-leakers at significant redshifts. Here we circumvent this uncertainty by inferring fesc using resolved (R > 4000) Lyman α (Lyα) profiles from the X-SHOOTER Lyα survey at z = 2 (XLS-z2). With empirically motivated criteria, we use Lyα profiles to select leakers ($f_{\mathrm{ esc}} > 20{{\ \rm per\ cent}}$) and non-leakers ($f_{\mathrm{ esc}} < 5{{\ \rm per\ cent}}$) from a representative sample of >0.2L* Lyman α emitters (LAEs). We use median stacked spectra of these subsets over λrest ≈ 1000–8000 Å to investigate the conditions for LyC fesc. Our stacks show similar mass, metallicity, MUV, and βUV. We find the following differences between leakers versus non-leakers: (i) strong nebular C iv and He ii emission versus non-detections; (ii) [O iii]/[O ii] ≈ 8.5 versus ≈3; (iii) Hα/Hβ indicating no dust versus E(B − V) ≈ 0.3; (iv) Mg ii emission close to the systemic velocity versus redshifted, optically thick Mg ii; and (v) Lyα fesc of ${\approx} 50{{\ \rm per\ cent}}$ versus ${\approx} 10{{\ \rm per\ cent}}$. The extreme equivalent widths (EWs) in leakers ([O iii]+$\mathrm{ H}\beta \approx 1100$ Å rest frame) constrain the characteristic time-scale of LyC escape to ≈3–10 Myr bursts when short-lived stars with the hardest ionizing spectra shine. The defining traits of leakers – extremely ionizing stellar populations, low column densities, a dust-free, high-ionization state interstellar medium (ISM) – occur simultaneously in the $f_{\rm esc} > 20{{\ \rm per\ cent}}$ stack, suggesting they are causally connected, and motivating why indicators like [O iii]/[O ii] may suffice to constrain fesc at z > 6 with the James Webb Space Telescope (JWST). The leakers comprise half of our sample, have a median LyC$f_{\rm esc} \approx 50{{\ \rm per\ cent}}$ (conservative range: $20\!-\!55{{\ \rm per\ cent}}$), and an ionizing production efficiency $\log ({\xi _{\rm {ion}}/\rm {Hz\ erg^{-1}}})\approx 25.9$ (conservative range: 25.7–25.9). These results show LAEs – the type of galaxies rare at z ≈ 2, but that become the norm at higher redshift – are highly efficient ionizers, with extreme ξion and prolific fesc occurring in sync.

     
    more » « less
  4. 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
  5. ABSTRACT The mean free path of ionizing photons, λmfp, is a key factor in the photoionization of the intergalactic medium (IGM). At z ≳ 5, however, λmfp may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of λmfp that address this bias and extend up to z ∼ 6 for the first time. Our measurements at z ∼ 5 are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At z ∼ 6 we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $\lambda _{\rm mfp} = 9.09^{+1.62}_{-1.28}$ proper Mpc and $0.75^{+0.65}_{-0.45}$ proper Mpc (68 per cent confidence) at z = 5.1 and 6.0, respectively. The results at z = 5.1 are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At z = 6.0, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at z = 6.0 falls well below extrapolations from lower redshifts, indicating rapid evolution in λmfp over 5 < z < 6. This evolution disfavours models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by z = 6, but is qualitatively consistent with models wherein reionization completed at z = 6 or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20 per cent neutral at z = 6, although our measurement at z = 6.0 is even lower than these models prefer. 
    more » « less