No Arabic abstract
At redshifts beyond z>6, as the mean fraction of neutral hydrogen x_HI in the intergalactic medium (IGM) increases, the line flux of Lyman alpha (Lya) emitters can be significantly suppressed, which can result in a decrease in the observed number of emitters at a given Lya flux. However, cosmological HII regions surrounding the Lya emitting galaxies alleviate these effects. We use simple models of the Lya line suppression that incorporate the presence of HII regions to predict the overall effect of the Lya absorption on the Lya luminosity function. We find, in agreement with other recent studies, that unless ionizing sources are unusually strongly clustered, a fully neutral IGM may be inconsistent with the large abundance of confirmed z=6.5 Lya emitters. However, the presence of local HII regions prohibits placing a tight constraint on the mean neutral fraction. We find x_HI < 0.25; the presence of strong winds and/or the clustering of ionizing sources would further weaken this constraint. We conclude that the evolution of the Lya LF is consistent with reionization occurring near this redshift. Finally, we suggest that a measurement of observed Lya line width as a function of the Lya luminosity, in a future, larger sample of Lya emitters, may serve as a robust diagnostic of the neutral fraction in the IGM.
The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Ly$alpha$) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Ly$alpha$~photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Ly$alpha$~visibility. The galaxy Ly$alpha$~luminosity function (LF) is thus a useful tool to constrain the timeline of reionization. In this work, we model the Ly$alpha$~LF as a function of redshift, $z=5-10$, and average IGM neutral hydrogen fraction, $overline{x}_textsc{hi}$. We combine the Ly$alpha$~luminosity probability distribution obtained from inhomogeneous reionization simulations with a model for the UV LF to model the Ly$alpha$~LF. As the neutral fraction increases, the average number density of Ly$alpha$~emitting galaxies decreases, and are less luminous, though for $overline{x}_textsc{hi} lesssim 0.4$ there is only a small decrease of the Ly$alpha$~LF. We use our model to infer the IGM neutral fraction at $z=6.6, 7.0, 7.3$ from observed Ly$alpha$~LFs. We conclude that there is a significant increase in the neutral fraction with increasing redshift: $overline{x}_textsc{hi}(z=6.6)=0.08^{+ 0.08}_{- 0.05}, , overline{x}_textsc{hi}(z=7.0)=0.28 pm 0.05$ and $overline{x}_textsc{hi}(z=7.3)=0.83^{+ 0.06}_{- 0.07}$. We predict trends in the Ly$alpha$~luminosity density and Schechter parameters as a function of redshift and the neutral fraction. We find that the Ly$alpha$~luminosity density decreases as the universe becomes more neutral. Furthermore, as the neutral fraction increases, the faint-end slope of the Ly$alpha$~LF steepens, and the characteristic Ly$alpha$~luminosity shifts to lower values, concluding that the evolving shape of the Ly$alpha$~LF -- not just its integral -- is an important tool to study reionization.
We build a physical model for high-redshift Lyman Alpha emitters (LAEs) by coupling state of the art cosmological simulations (GADGET-2) with a dust model and a radiative transfer code (pCRASH). We post-process the cosmological simulation with pCRASH using five different values of the escape fraction of hydrogen ionizing photons (f_esc=0.05,0.25,0.5,0.75,0.95) until reionization is complete, i.e. the average neutral hydrogen fraction drops to <X_HI>~10^-4. Then, the only free-parameter left to match model results to the observed Lya and UV luminosity functions of LAEs at z~6.6 is the relative escape of Lyman Alpha (Lya) and continuum photons from the galactic environment (f_alpha/f_c). We find a three-dimensional degeneracy such that the theoretical model can be reconciled with observations for an IGM Lya transmission <T_alpha>_LAE~38-50% (which translates to <X_HI>~0.5-10^-4 for Gaussian emission lines), f_esc~0.05-0.50 and f_alpha/f_c~0.6-1.8.
We present the luminosity function (LF) for ultraluminous Ly$alpha$ emitting galaxies (LAEs) at z = 6.6. We define ultraluminous LAEs (ULLAEs) as galaxies with logL(Ly$alpha$) > 43.5 erg s$^{-1}$. We select our main sample using the g, r, i, z, and NB921 observations of a wide-area (30 deg$^2$) Hyper Suprime-Cam survey of the North Ecliptic Pole (NEP) field. We select candidates with g, r, i > 26, NB921 $leq$ 23.5, and NB921 - z $leq$ 1.3. Using the DEIMOS spectrograph on Keck II, we confirm 9 of our 14 candidates as ULLAEs at z = 6.6 and the remaining 5 as an AGN at z = 6.6, two [OIII]$lambda$5007 emitting galaxies at z = 0.84 and z = 0.85, and two non-detections. This emphasizes the need for full spectroscopic follow-up to determine accurate LFs. In constructing the ULLAE LF at z = 6.6, we combine our 9 NEP ULLAEs with two previously discovered and confirmed ULLAEs in the COSMOS field: CR7 and COLA1. We apply rigorous corrections for incompleteness based on simulations. We compare our ULLAE LF at z = 6.6 with LFs at z = 5.7 and z = 6.6 from the literature. Our data reject some previous LF normalizations and power law indices, but they are broadly consistent with others. Indeed, a comparative analysis of the different literature LFs suggests that none is fully consistent with any of the others, making it critical to determine the evolution from z = 5.7 to z = 6.6 using LFs constructed in exactly the same way at both redshifts.
We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$alpha$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$beta$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$beta$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$alpha$ and Lyman-$gamma$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively.
Lyman alpha (Lya) emission lines should be attenuated in a neutral intergalactic medium (IGM). Therefore the visibility of Lya emitters at high redshifts can serve as a valuable probe of reionization at about the 50% level. We present an imaging search for z=7.7 Lya emitting galaxies using an ultra-narrowband filter (filter width= 9A) on the NEWFIRM imager at the Kitt Peak National Observatory. We found four candidate Lya emitters in a survey volume of 1.4 x 10^4 Mpc^3, with a line flux brighter than 6x10^-18 erg/cm^2/s (5 sigma in 2 aperture). We also performed a detailed Monte-Carlo simulation incorporating the instrumental effects to estimate the expected number of Lya emitters in our survey, and found that we should expect to detect one Lya emitter, assuming a non-evolving Lya luminosity function (LF) between z=6.5 and z=7.7. Even if one of the present candidates is spectroscopically confirmed as a z~8 Lya emitter, it would indicate that there is no significant evolution of the Lya LF from z=3.1 to z~8. While firm conclusions would need both spectroscopic confirmations and larger surveys to boost the number counts of galaxies, we successfully demonstrate the feasibility of sensitive near-infrared (1.06 um) narrow-band searches using custom filters designed to avoid the OH emission lines that make up most of the sky background.