Intergalactic Medium Neutral Fraction Research Articles

IGM damping wing constraints on the tail end of reionization from the enlarged XQR-30 sample

ABSTRACT The attenuation of Lyα photons by neutral hydrogen in the intergalactic medium (IGM) at z ≳ 5 continues to be a powerful probe for studying the epoch of reionization. Given a framework to estimate the intrinsic (true) Lyα emission of high-z sources, one can infer the ionization state of the IGM during reionization. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high signal-to-noise quasar spectra between $5.8\lesssim \, z\lesssim \, 6.6$ obtained with VLT/X-shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Lyα emission profiles (obtained via covariance matrix reconstruction of the Lyα and N v emission lines from unattenuated high-ionization emission line profiles; C iv, Si iv + O iv], and C iii]) and (ii) distribution of ionized regions within the IGM using synthetic damping wing profiles drawn from a 1.63 Gpc3 reionization simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionization. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14 -0.12}$ and $0.29\substack{+0.14 -0.13}$ at z = 6.15 and 6.35. Further, we also report 68th percentile upper limits of $\bar{x}_{\mathrm{H\, {\small I}}{}} \lt 0.21$, 0.20, 0.21, and 0.18 at z = 5.8, 5.95, 6.05, and 6.55. These results imply reionization is still ongoing at $5.8\lesssim \, z\lesssim \, 6.55$, consistent with previous results from XQR-30 (dark fraction and Lyα forest) along with other observational probes considered in the literature.

JADES: The emergence and evolution of Lyα emission and constraints on the intergalactic medium neutral fraction

The rest-frame UV recombination emission line Lyα can be powered by ionising photons from young massive stars in star-forming galaxies, but the fact that it can be resonantly scattered by neutral gas complicates its interpretation. For reionisation-era galaxies, a neutral intergalactic medium will scatter Lyα from the line of sight, making Lyα a useful probe of the neutral fraction evolution. Here, we explore Lyα in JWST/NIRSpec spectra from the ongoing JADES programme, which targets hundreds of galaxies in the well-studied GOODS-S and GOODS-N fields. These sources are UV-faint (−20.4 < MUV < −16.4) and thus represent a poorly explored class of galaxy. We fitted the low spectral resolution spectra (R ∼ 100) of a subset of 84 galaxies in GOODS-S with zspec > 5.6 (as derived with optical lines) with line and continuum models to search for significant line emission. Through exploration of the R100 data, we find evidence for Lyα in 17 sources. This sample allowed us to place observational constraints on the fraction of galaxies with Lyα emission in the redshift range 5.6 < z < 7.5, with a decrease from z = 6 to z = 7. We also find a positive correlation between the Lyα equivalent width and MUV, as seen in other samples. We used these results to estimate the neutral gas fraction at z ∼ 7, and our estimates are in agreement with previous results (XHI ∼ 0.5 − 0.9).

The Impact of Cosmic Variance on Inferences of Global Neutral Fraction Derived from Lyα Luminosity Functions during Reionization

We investigate the impact of field-to-field variation, deriving from cosmic variance, in measured Lyα emitter (LAE) luminosity functions (LFs) and this variation’s impact on inferences of the neutral fraction of the intergalactic medium (IGM) during reionization. We post-process a z = 7 IGM simulation to populate the dark matter halos with LAEs. These LAEs have realistic UV magnitudes, Lyα fluxes, and Lyα line profiles. We calculate the attenuation of Lyα emission in universes with varying IGM neutral fraction, . In a = 0.3 simulation, we perform 100 realizations of a mock 2 deg2 survey with a redshift window Δz = 0.5 and flux limit f Lyα > 1 × 10−17 erg s−1 cm−2; such a survey is typical in depth and volume of the largest LAE surveys conducted today. For each realization, we compute the LAE LF and use it to recover the input . Comparing the inferred values of across the ensemble of the surveys, we find that cosmic variance, deriving from large-scale structure and variation in the neutral gas along the sightline, imposes a floor in the uncertainty of when = 0.3. We explore mitigation strategies to decrease this uncertainty, such as increasing the volume, decreasing the flux limit, or probing the volume with many independent fields. Increasing the area and/or depth of the survey does not mitigate the uncertainty, but composing a survey with many independent fields is effective. This finding highlights the best strategy for LAE surveys aiming at constraining the of the universe during reionization.

LIMFAST. I. A Seminumerical Tool for Line Intensity Mapping

We present LIMFAST, a seminumerical code for simulating high-redshift galaxy formation and cosmic reionization as revealed by multitracer line intensity mapping (LIM) signals. LIMFAST builds upon and extends the 21cmFAST code widely used for 21 cm cosmology by implementing state-of-the-art models of galaxy formation and evolution. The metagalactic radiation background, including the production of various star formation lines, together with the 21 cm line signal tracing the neutral intergalactic medium (IGM), is self-consistently described by photoionization modeling and stellar population synthesis coupled to the galaxy formation model. We introduce basic structure and functionalities of the code, and demonstrate its validity and capabilities by showing broad agreements between the predicted and observed evolution of cosmic star formation, IGM neutral fraction, and metal enrichment. We also present the LIM signals of 21 cm, Lyα, Hα, Hβ, [O ii], and [O iii] lines simulated by LIMFAST, and compare them with results from the literature. We elaborate on how several major aspects of our modeling framework, including models of star formation, chemical enrichment, and photoionization, may impact different LIM observables and thus become testable once applied to observational data. LIMFAST aims at being an efficient and resourceful tool for intensity mapping studies in general, exploring a wide range of scenarios of galaxy evolution and reionization and frequencies over which useful cosmological signals can be measured.

(Nearly) Model-independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at z ∼ 5–7 Using Dark Pixel Fractions in Lyα and Lyβ Forests

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at z > 7 and largely ionized by z ∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low () or close to unity (). In particular, the neutral fraction evolution of the IGM at the critical redshift range of z = 6–7 is poorly constrained. We present new constraints on at z ∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 < z < 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyα and Lyβ forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at z ∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of (1σ), (1σ), and (1σ). The dark pixel fractions at z > 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

Crucial Factors for Lyα Transmission in the Reionizing Intergalactic Medium: Infall Motion, H ii Bubble Size, and Self-shielded Systems

Abstract Using the CoDa II simulation, we study the Lyα transmissivity of the intergalactic medium (IGM) during reionization. At z &gt; 6, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Lyα, suppressing the transmission up to roughly the circular velocity of the galaxy. In some sight lines, an optically thin blob generated by a supernova in a neighboring galaxy results in a peak feature, which can be mistaken for a blue peak. Redward of the resonance absorption, the damping-wing opacity correlates with the global IGM neutral fraction and the UV magnitude of the source galaxy. Brighter galaxies tend to suffer lower opacity because they tend to reside in larger H ii regions, and the surrounding IGM transmits redder photons, which are less susceptible to attenuation, owing to stronger infall velocity. The H ii regions are highly nonspherical, causing both sight-line-to-sight-line and galaxy-to-galaxy variation in opacity. Also, self-shielded systems within H ii regions strongly attenuate the emission for certain sight lines. All these factors add to the transmissivity variation, requiring a large sample size to constrain the average transmission. The variation is largest for fainter galaxies at higher redshift. The 68% range of the transmissivity is similar to or greater than the median for galaxies with M UV ≥ −21 at z ≥ 7, implying that more than a hundred galaxies would be needed to measure the transmission to 10% accuracy.

The Evolution of the Lyman-alpha Luminosity Function during Reionization

Abstract The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Lyα) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Lyα photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Lyα visibility. The galaxy Lyα luminosity function (LF) is thus a useful tool to constrain the timeline of reionization. In this work, we model the Lyα LF as a function of redshift, z = 5–10, and average IGM neutral hydrogen fraction, x ¯ Hɪ . We combine the Lyα luminosity probability distribution obtained from inhomogeneous reionization simulations with a model for the UV LF to model the Lyα LF. As the neutral fraction increases, the average number density of Lyα emitting galaxies decreases, and are less luminous, though for x ¯ Hɪ ≲ 0.4 there is only a small decrease in the Lyα LF. We use our model to infer the IGM neutral fraction at z = 6.6, 7.0, and 7.3 from observed Lyα LFs. We conclude that there is a significant increase in the neutral fraction with increasing redshift: x ¯ Hɪ ( z = 6.6 ) = 0.08 − 0.05 + 0.08 , x ¯ Hɪ ( z = 7.0 ) = 0.28 ± 0.05 and x ¯ Hɪ ( z = 7.3 ) = 0.83 − 0.07 + 0.06 . We predict trends in the Lyα luminosity density and Schechter parameters as a function of redshift and the neutral fraction. We find that the Lyα luminosity density decreases as the universe becomes more neutral. Furthermore, as the neutral fraction increases, the faint-end slope of the Lyα LF steepens, and the characteristic Lyα luminosity shifts to lower values; hence, we conclude that the evolving shape of the Lyα LF—not just its integral—is an important tool to study reionization.

Searching for the shadows of giants – II. The effect of local ionization on the Ly α absorption signatures of protoclusters at redshift z ∼ 2.4

ABSTRACT Local variations in the intergalactic medium (IGM) neutral hydrogen fraction will affect the Ly α absorption signature of protoclusters identified in tomographic surveys. Using the IllustrisTNG simulations, we investigate how the AGN proximity effect and hot, collisionally ionized gas arising from gravitational infall and black hole feedback changes the Ly α absorption associated with $M_{z=0}\simeq 10^{14}\, {\rm M}_\odot$ protoclusters at z ≃ 2.4. We find that protocluster galaxy overdensities exhibit a weak anticorrelation with Ly α transmission in IGM transmission maps, but local H$\, \rm \scriptstyle I$ ionization enhancements due to hot $T\gt 10^{6}\rm \, K$ gas or nearby AGN can disrupt this relationship within individual protoclusters. On average, however, we find that strong reductions in the IGM neutral fraction are limited to within $\lesssim 5h^{-1}\, \textrm {cMpc}$ of the dark matter haloes. Local ionization enhancements will therefore have a minimal impact on the completeness of protocluster identification in tomographic surveys if smoothing Ly α transmission maps over scales of $\sim 4 h^{-1}\, \textrm {cMpc}$, as is typically done in observations. However, if calibrating the relationship between the matter density and Ly α transmission in tomographic maps using simple analytical models for the Ly α forest opacity, the presence of hot gas around haloes can still result in systematically lower estimates of Mz = 0 for the most massive protoclusters.

Ly α as a tracer of cosmic reionization in the SPHINX radiation-hydrodynamics cosmological simulation

ABSTRACT The Ly α emission line is one of the most promising probes of cosmic reionization but isolating the signature of a change in the ionization state of the intergalactic medium (IGM) is challenging because of intrinsic evolution and internal radiation transfer effects. We present the first study of the evolution of Ly α emitters (LAE) during the epoch of reionization based on a full radiation-hydrodynamics cosmological simulation that is able to capture both the large-scale process of reionization and the small-scale properties of galaxies. We predict the Ly α emission of galaxies in the 103 cMpc3sphinx simulation at 6 ≤ z ≤ 9 by computing the full Ly α radiation transfer from interstellar medium (ISM) to IGM scales. sphinx is able to reproduce many observational constraints such as the UV/Ly α luminosity functions and stellar mass functions at z ≳ 6 for the dynamical range probed by our simulation (M1500 ≳ −18, LLy α ≲ 1042 erg s−1, M⋆ ≲ 109 M⊙). As intrinsic Ly α emission and internal Ly α escape fractions barely evolve from z = 6–9, the observed suppression of Ly α luminosities with increasing redshift is fully attributed to IGM absorption. For most observable galaxies (M1500 ≲ −16), the Ly α line profiles are slightly shifted to the red due to internal radiative transfer effects that mitigates the effect of IGM absorption. Overall, the enhanced Ly α suppression during reionization traces the IGM neutral fraction $x_{{\rm H\, {i}}}$ well, but the predicted amplitude of this reduction is a strong function of the Ly α peak shift, which is set at ISM/circumgalactic medium scales. We find that a large number of LAEs could be detectable in very deep surveys during reionization when $x_{{\rm H\, {i}}}$ is still $\approx 50{{\ \rm per\ cent}}$.

Can Galaxy Evolution Mimic Cosmic Reionization?

Abstract Lyα emitting galaxies are powerful tools to probe the late stages of cosmic reionization. The observed sudden drop in Lyα fraction at z &gt; 6 is often interpreted as a sign of reionization, since the intergalactic medium (IGM) is more neutral and opaque to Lyα photons. Crucially, this interpretation of the observations is only valid under the assumption that galaxies themselves experience a minimal evolution at these epochs. By modeling Lyα radiative transfer effects in and around galaxies, we examine whether a change in the galactic properties can reproduce the observed drop in the Lyα fraction. We find that an increase in the galactic neutral hydrogen content or a reduction in the outflow velocity toward higher redshift both lead to a lower Lyα escape fraction, and can thus mimic an increasing neutral fraction of the IGM. We furthermore find that this change in galactic properties leads to systematically different Lyα spectra which can be used to differentiate the two competing effects. Using the CANDELSz7 survey measurements which indicate slightly broader lines at z ∼ 6, we find that the scenario of a mere increase in the galactic column density toward higher z is highly unlikely. We also show that a decrease in outflow velocity is not ruled out by existing data but leads to more prominent blue peaks at z &gt; 6. Our results caution using Lyα observations to estimate the IGM neutral fraction without accounting for the potential change in the galactic properties, e.g., by mapping out the evolution of Lyα spectral characteristics.

Interpreting LOFAR 21-cm signal upper limits at z ≈ 9.1 in the context of high-z galaxy and reionization observations

ABSTRACT Using the latest upper limits on the 21-cm power spectrum at z ≈ 9.1 from the Low Frequency Array (LOFAR), we explore the regions of parameter space which are inconsistent with the data. We use 21cmmc, a Monte Carlo Markov chain sampler of 21cmfast which directly forward models the three dimensional (3D) cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parametrization from 21cmfast, which includes mass-dependent star formation rates and ionizing escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-z galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionization such as ultra-violet (UV) luminosity functions, background UV photoionization rate, intergalactic medium (IGM) neutral fraction, and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at z ≈ 9.1 are excluded at ≳2σ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below ∼2.6 K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{\rm H\, \rm {\small I}} \lesssim 0.6$. Note, these limits are only obtained from 141 h of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21cmmc in combining 21-cm observations with other high-z probes to constrain the astrophysics of galaxies.

Exploring reionization and high-z galaxy observables with recent multiredshift MWA upper limits on the 21-cm signal

ABSTRACT We use the latest multiredshift (z = 6.5−8.7) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. We explore these limits in the context of reionization astrophysics by using 21CMMC to connect the disfavoured regions of parameter space to existing observational constraints on reionization such as high-z galaxy ultraviolet (UV) luminosity functions, the background UV photoionization rate, the intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of disfavoured models to already be inconsistent with existing observational constraints. These can be broadly classified into two types of models: (i) ‘cold’ reionization and (ii) pure matter density fluctuations in a cold, neutral IGM (i.e. no reionization). Interestingly, a small subsample of models inconsistent with the MWA is consistent with the aforementioned constraints (excluding the X-ray emissivity). This implies that the current MWA limits are already providing unique information to disfavour models of reionization, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding 1σ lower limits of ϵX, 0.5−2 keV ≳ 1034.5 erg s−1 Mpc−3. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, and 2.4 K at z = 6.5, 6.8, 7.1, 7.8, 8.2, and 8.7. With this, we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on ϵX, 0.5–2 keV and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction.

A Significantly Neutral Intergalactic Medium Around the Luminous z = 7 Quasar J0252–0503

Abstract Luminous z ≥ 7 quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Lyα damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction toward that line of sight. However, damping wing features have only been detected in two z &gt; 7 quasars in previous studies. In this paper, we present new high-quality optical and near-infrared spectroscopy of the z = 7.00 quasar DES J025216.64–050331.8 obtained with Keck/Near-Infrared Echellette Spectrometer and Gemini/GMOS. By using the Mg ii single-epoch virial method, we find that it hosts a SMBH accreting at an Eddington ratio of λ Edd = 0.7 ± 0.1, consistent with the values seen in other luminous z ∼ 7 quasars. Furthermore, the Lyα region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at z = 7 of within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of . This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

The spin–temperature dependence of the 21-cm–LAE cross-correlation

ABSTRACT Cross-correlating 21 cm with known cosmic signals will be invaluable proof of the cosmic origin of the first 21-cm detections. As some of the widest fields available, comprising thousands of sources with reasonably known redshifts, narrow-band Lyman-α emitter (LAE) surveys are an obvious choice for such cross-correlation. Here, we revisit the 21-cm–LAE cross-correlation, relaxing the common assumption of reionization occurring in a pre-heated intergalactic medium (IGM). Using specifications from the Square Kilometre Array and the Subaru Hyper Supreme-Cam, we present new forecasts of the 21-cm–LAE cross-correlation function at z ∼ 7. We sample a broad parameter space of the mean IGM neutral fraction and spin temperature, ($\bar{x}_{\rm H\,{\small I}}$, $\bar{T}_{\rm S}$). The sign of the cross-correlation roughly follows the sign of the 21-cm signal: Ionized regions that surround LAEs correspond to relative hot spots in the 21-cm signal when the neutral IGM is colder than the CMB, and relative cold spots when the neutral IGM is hotter than the CMB. The amplitude of the cross-correlation function generally increases with increasing $\bar{x}_{\rm H\,{\small I}}$, following the increasing bias of the cosmic H ii regions. As is the case for 21 cm, the strongest cross signal occurs when the IGM is colder than the CMB, providing a large contrast between the neutral regions and the ionized regions, which host LAEs. We also vary the topology of reionization and the epoch of X-ray heating. The cross-correlation during the first half of reionization is sensitive to these topologies, and could thus be used to constrain them.

Inferences on the timeline of reionization at z ∼ 8 from the KMOS Lens-Amplified Spectroscopic Survey

Abstract Detections and non-detections of Lyman alpha (Lyα) emission from z &amp;gt; 6 galaxies (&amp;lt;1 Gyr after the big bang) can be used to measure the timeline of cosmic reionization. Of key interest to measuring reionization’s mid-stages, but also increasing observational challenge, are observations at z &amp;gt; 7, where Lyα redshifts to near infra-red wavelengths. Here we present a search for z &amp;gt; 7.2 Lyα emission in 53 intrinsically faint Lyman Break Galaxy candidates, gravitationally lensed by massive galaxy clusters, in the KMOS Lens-Amplified Spectroscopic Survey (KLASS). With integration times of ∼7–10 h, we detect no Lyα emission with signal-to-noise ratio (S/N) &amp;gt; 5 in our sample. We determine our observations to be 80 per cent complete for 5σ spatially and spectrally unresolved emission lines with integrated line flux &amp;gt;5.7 × 10−18 erg s−1 cm−2. We define a photometrically selected sub-sample of 29 targets at z = 7.9 ± 0.6, with a median 5σ Lyα EW limit of 58 Å. We perform a Bayesian inference of the average intergalactic medium (IGM) neutral hydrogen fraction using their spectra. Our inference accounts for the wavelength sensitivity and incomplete redshift coverage of our observations, and the photometric redshift probability distribution of each target. These observations, combined with samples from the literature, enable us to place a lower limit on the average IGM neutral hydrogen fraction of $\gt 0.76 \,\, (68{{\ \rm per\ cent}}), \,\, \gt 0.46 \,\, (95{{\ \rm per\ cent}})$ at z ∼ 8, providing further evidence of rapid reionization at z ∼ 6–8. We show that this is consistent with reionization history models extending the galaxy luminosity function to $M_ \rm {\small UV}\lesssim -12$, with low ionizing photon escape fractions, $f_\textrm{esc} \lesssim 15{{\ \rm per\ cent}}$.

The Universe Is Reionizing at z ∼ 7: Bayesian Inference of the IGM Neutral Fraction Using Lyα Emission from Galaxies

Abstract We present a new flexible Bayesian framework for directly inferring the fraction of neutral hydrogen in the intergalactic medium (IGM) during the Epoch of Reionization (EoR, z ∼ 6–10) from detections and non-detections of Lyman Alpha (Lyα) emission from Lyman Break galaxies (LBGs). Our framework combines sophisticated reionization simulations with empirical models of the interstellar medium (ISM) radiative transfer effects on Lyα. We assert that the Lyα line profile emerging from the ISM has an important impact on the resulting transmission of photons through the IGM, and that these line profiles depend on galaxy properties. We model this effect by considering the peak velocity offset of Lyα lines from host galaxies’ systemic redshifts, which are empirically correlated with UV luminosity and redshift (or halo mass at fixed redshift). We use our framework on the sample of LBGs presented in Pentericci et al. and infer a global neutral fraction at z ∼ 7 of , consistent with other robust probes of the EoR and confirming that reionization is ongoing ∼700 Myr after the Big Bang. We show that using the full distribution of Lyα equivalent width detections and upper limits from LBGs places tighter constraints on the evolving IGM than the standard Lyα emitter fraction, and that larger samples are within reach of deep spectroscopic surveys of gravitationally lensed fields and James Webb Space Telescope NIRSpec.

VLT/FORS2 view at z ~ 6: Lyman-α emitter fraction and galaxy physical properties at the edge of the epoch of cosmic reionization

The fraction of Lyman-α emitters (LAEs) among the galaxy population has been found to increase from z ~ 0 to z ~ 6 and drop dramatically at z> 6. This drop has been interpreted as an effect of an increasingly neutral intergalactic medium (IGM) with increasing redshift, while a Lyman continuum escape fraction evolving with redshift and/or a sudden change of galaxy physical properties can also contribute to the decreasing LAE fraction. We report the result of a large VLT/FORS2 program aiming to confirm spectroscopically a large galaxy sample at z ≥ 6 that has been selected in several independent fields through the Lyman break technique. Combining those data with archival data, we create a large and homogeneous sample of z ~ 6 galaxies (N = 127), complete in terms of Lyα detection at > 95% for Lyα equivalent width EW(Lyα) ≥ 25 Å. We use this sample to derive a new measurement of the LAE fraction at z ~ 6 and derive the physical properties of these galaxies through spectral energy distribution (SED) fitting. We find a median LAE fraction at z ~ 6 lower than in previous studies, while our sample exhibits typical properties for z ~ 6 galaxies in terms of UV luminosity and UVβ slope. The comparison of galaxy physical properties between LAEs and non-LAEs is comparable to results at lower redshift: LAEs with the largest EW(Lyα) exhibit bluer UV slopes, are slightly less massive and less star-forming. The main difference between LAEs and non-LAEs is that the latter are significantly dustier. Using predictions of our SED fitting code accounting for nebular emission, we find an effective Lyα escape fraction fesceff(Lyα) = 0.23-0.17+0.36 remarkably consistent with the value derived by comparing UV luminosity function with Lyα luminosity function. We conclude that the drop in the LAE fraction from z ~ 6 to z> 6 is less dramatic than previously found and the effect of an increasing IGM neutral fraction is possibly observed at 5 <z< 6. The processes driving the escape of Lyα photons at z ~ 6 are similar to those at lower redshifts and based on our derived fesceff(Lyα), we find that the IGM has a relatively small impact on Lyα photon visibility at z ~ 6, with a lower limit for the IGM transmission to Lyα photons, TIGM ≳ 0.20, likely due to the presence of outflows.

First Results from the Lyman Alpha Galaxies in the Epoch of Reionization (LAGER) Survey: Cosmological Reionization at z ∼ 7

Abstract We present the first results from the ongoing Lyman Alpha Galaxies in the Epoch of Reionization (LAGER) project, which is the largest narrowband survey for z ∼ 7 galaxies to date. Using a specially built narrowband filter NB964 for the superb large-area Dark Energy Camera (DECam) on the NOAO/CTIO 4 m Blanco telescope, LAGER has collected 34 hr NB964 narrowband imaging data in the 3 deg2 COSMOS field. We have identified 23 Lyα Emitter candidates at z = 6.9 in the central 2-deg2 region, where DECam and public COSMOS multi-band images exist. The resulting luminosity function (LF) can be described as a Schechter function modified by a significant excess at the bright end (four galaxies with L Lyα ∼ 1043.4±0.2 erg s−1). The number density at L Lyα ∼ 1043.4±0.2 erg s−1 is little changed from z = 6.6, while at fainter L Lyα it is substantially reduced. Overall, we see a fourfold reduction in Lyα luminosity density from z = 5.7 to z = 6.9. Combined with a more modest evolution of the continuum UV luminosity density, this suggests a factor of ∼3 suppression of Lyα by radiative transfer through the z ∼ 7 intergalactic medium (IGM). It indicates an IGM neutral fraction of x H i ∼ 0.4–0.6 (assuming Lyα velocity offsets of 100–200 km s−1). The changing shape of the Lyα LF between z ≲ 6.6 and z = 6.9 supports the hypothesis of ionized bubbles in a patchy reionization at z ∼ 7.

Are we witnessing the epoch of reionisation at z = 7.1 from the spectrum of J1120+0641?

We quantify the presence of Ly\alpha\ damping wing absorption from a partially-neutral intergalactic medium (IGM) in the spectrum of the $z=7.08$ QSO, ULASJ1120+0641. Using a Bayesian framework, we simultaneously account for uncertainties in: (i) the intrinsic QSO emission spectrum; and (ii) the distribution of cosmic HI patches during the epoch of reionisation (EoR). For (i) we use a new intrinsic Ly\alpha\ emission line reconstruction method (Greig et al.), sampling a covariance matrix of emission line properties built from a large database of moderate-$z$ QSOs. For (ii), we use the Evolution of 21-cm Structure (EOS; Mesinger et al.) simulations, which span a range of physically-motivated EoR models. We find strong evidence for the presence of damping wing absorption redward of Ly\alpha\ (where there is no contamination from the Ly\alpha\ forest). Our analysis implies that the EoR is not yet complete by $z=7.1$, with the volume-weighted IGM neutral fraction constrained to $\bar{x}_{\rm H\,{\scriptsize I}} = 0.40\substack{+0.21 -0.19}$ at $1\sigma$ ($\bar{x}_{\rm H\,{\scriptsize I}} = 0.40\substack{+0.41 -0.32}$ at $2\sigma$). This result is insensitive to the EoR morphology. Our detection of significant neutral HI in the IGM at $z=7.1$ is consistent with the latest Planck 2016 measurements of the CMB Thompson scattering optical depth (Planck Collaboration XLVII).

Constraints on the temperature of the intergalactic medium atz= 8.4 with 21-cm observations

We compute robust lower limits on the spin temperature, $T_{\rm S}$, of the $z=8.4$ intergalactic medium (IGM), implied by the upper limits on the 21-cm power spectrum recently measured by PAPER-64. Unlike previous studies which used a single epoch of reionization (EoR) model, our approach samples a large parameter space of EoR models: the dominant uncertainty when estimating constraints on $T_{\rm S}$. Allowing $T_{\rm S}$ to be a free parameter and marginalizing over EoR parameters in our Markov Chain Monte Carlo code 21CMMC, we infer $T_{\rm S}\ge3 {\rm K}$ (corresponding approximately to $1\sigma$) for a mean IGM neutral fraction of $\bar{x}_{\rm H{\scriptsize I}}\gtrsim0.1$. We further improve on these limits by folding-in additional EoR constraints based on: (i) the dark fraction in QSO spectra, which implies a strict upper limit of $\bar{x}_{\rm H{\scriptsize I}}[z=5.9]\leq 0.06+0.05 \,(1\sigma)$; and (ii) the electron scattering optical depth, $\tau_{\rm e}=0.066\pm0.016\,(1\sigma)$ measured by the Planck satellite. By restricting the allowed EoR models, these additional observations tighten the approximate $1\sigma$ lower limits on the spin temperature to $T_{\rm S} \ge 6$ K. Thus, even such preliminary 21-cm observations begin to rule out extreme scenarios such as `cold reionization', implying at least some prior heating of the IGM. The analysis framework developed here can be applied to upcoming 21-cm observations, thereby providing unique insights into the sources which heated and subsequently reionized the very early Universe.