HI Reionization Research Articles

The ALPINE-ALMA [CII] survey

The ALMA-ALPINE [CII] survey is aimed at characterizing the properties of a sample of normal star-forming galaxies (SFGs). The ALMA Large Program to INvestigate (ALPINE) features 118 galaxies observed in the [CII]-158 μm line and far infrared (FIR) continuum emission during the period of rapid mass assembly, right after the end of the HI reionization, at redshifts of 4 < z < 6. We present the survey science goals, the observational strategy, and the sample selection of the 118 galaxies observed with ALMA, with an average beam minor axis of about 0.85″, or ∼5 kpc at the median redshift of the survey. The properties of the sample are described, including spectroscopic redshifts derived from the UV-rest frame, stellar masses, and star-formation rates obtained from a spectral energy distribution (SED) fitting. The observed properties derived from the ALMA data are presented and discussed in terms of the overall detection rate in [CII] and FIR continuum, with the observed signal-to-noise distribution. The sample is representative of the SFG population in the main sequence at these redshifts. The overall detection rate in [CII] is 64% for a signal-to-noise ratio (S/N) threshold larger than 3.5 corresponding to a 95% purity (40% detection rate for S/N > 5). Based on a visual inspection of the [CII] data cubes together with the large wealth of ancillary data, we find a surprisingly wide range of galaxy types, including 40% that are mergers, 20% extended and dispersion-dominated, 13% compact, and 11% rotating discs, with the remaining 16% too faint to be classified. This diversity indicates that a wide array of physical processes must be at work at this epoch, first and foremost, those of galaxy mergers. This paper sets a reference sample for the gas distribution in normal SFGs at 4 < z < 6, a key epoch in galaxy assembly, which is ideally suited for studies with future facilities, such as the James Webb Space Telescope (JWST) and the Extremely Large Telescopes (ELTs).

UV and Lyα luminosity functions of galaxies and star formation rate density at the end of HI reionization from the VIMOS UltraDeep Survey (VUDS)

Context. The star formation rate density (SFRD) evolution presents an area of great interest in the studies of galaxy evolution and reionization. The current constraints of SFRD at z > 5 are based on the rest-frame UV luminosity functions with the data from photometric surveys. The VIMOS UltraDeep Survey (VUDS) was designed to observe galaxies at redshifts up to ∼6 and opened a window for measuring SFRD at z > 5 from a spectroscopic sample with a well-controlled selection function. Aims. We establish a robust statistical description of the star-forming galaxy population at the end of cosmic HI reionization (5.0 ≤ z ≤ 6.6) from a large sample of 49 galaxies with spectroscopically confirmed redshifts. We determine the rest-frame UV and Lyα luminosity functions and use them to calculate SFRD at the median redshift of our sample z = 5.6. Methods. We selected a sample of galaxies at 5.0 ≤ zspec ≤ 6.6 from the VUDS. We cleaned our sample from low redshift interlopers using ancillary photometric data. We identified galaxies with Lyα either in absorption or in emission, at variance with most spectroscopic samples in the literature where Lyα emitters (LAE) dominate. We determined luminosity functions using the 1/Vmax method. Results. The galaxies in this redshift range exhibit a large range in their properties. A fraction of our sample shows strong Lyα emission, while another fraction shows Lyα in absorption. UV-continuum slopes vary with luminosity, with a large dispersion. We find that star-forming galaxies at these redshifts are distributed along the main sequence in the stellar mass vs. SFR plane, described with a slope α = 0.85 ± 0.05. We report a flat evolution of the specific SFR compared to lower redshift measurements. We find that the UV luminosity function is best reproduced by a double power law, while a fit with a Schechter function is only marginally inferior. The Lyα luminosity function is best fitted with a Schechter function. We derive a logSFRDUV(M⊙ yr−1 Mpc−3) = −1.45+0.06−0.08 and logSFRDLyα(M⊙ yr−1 Mpc−3) = −1.40+0.07−0.08. The SFRD derived from the Lyα luminosity function is in excellent agreement with the UV-derived SFRD after correcting for IGM absorption. Conclusions. Our new SFRD measurements at a mean redshift of z = 5.6 are ∼0.2 dex above the mean SFRD reported in Madau & Dickinson (2014, ARA&A, 52, 415), but in excellent agreement with results from Bouwens et al. (2015a, ApJ, 803, 34). These measurements confirm the steep decline of the SFRD at z > 2. The bright end of the Lyα luminosity function has a high number density, indicating a significant star formation activity concentrated in the brightest LAE at these redshifts. LAE with equivalent width EW > 25 Å contribute to about 75% of the total UV-derived SFRD. While our analysis favors low dust content in 5.0 < z < 6.6, uncertainties on the dust extinction correction and associated degeneracy in spectral fitting will remain an issue, when estimating the total SFRD until future surveys extending spectroscopy to the NIR rest-frame spectral domain, such as with JWST.

Consistent modelling of the meta-galactic UV background and the thermal/ionization history of the intergalactic medium

Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$\alpha$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium.

Cosmic reionization after Planck II: contribution from quasars

In the light of the recent Planck downward revision of the electron scattering optical depth, and of the discovery of a faint AGN population at $z > 4$, we reassess the actual contribution of quasars to cosmic reionization. To this aim, we extend our previous MCMC-based data-constrained semi-analytic reionization model and study the role of quasars on global reionization history. We find that, the quasars can alone reionize the Universe only for models with very high AGN emissivities at high redshift. These models are still allowed by the recent CMB data and most of the observations related to HI reionization. However, they predict an extended and early HeII reionization ending at $z\gtrsim4$ and a much slower evolution in the mean HeII Ly-$\alpha$ forest opacity than what the actual observation suggests. Thus when we further constrain our model against the HeII Ly-$\alpha$ forest data, this AGN-dominated scenario is found to be clearly ruled out at 2-$\sigma$ limits. The data seems to favour a standard two-component picture where quasar contributions become negligible at $z\gtrsim6$ and a non-zero escape fraction of $\sim10\%$ is needed from early-epoch galaxies. For such models, mean neutral hydrogen fraction decreases to $\sim10^{-4}$ at $z=6.2$ from $\sim0.8$ at $z=10.0$ and helium becomes doubly ionized at much later time, $z\sim3$. We find that, these models are as well in good agreement with the observed thermal evolution of IGM as opposed to models with very high AGN emissivities.

Constraining the UV emissivity of AGN throughout cosmic time via X-ray surveys

The cosmological process of hydrogen (HI) reionization in the intergalactic medium is thought to be driven by UV photons emitted by star-forming galaxies and ionizing active galactic nuclei (AGN). The contribution of QSOs to HI reionization at $z>4$ has been traditionally believed to be quite modest. However, this view has been recently challenged by new estimates of a higher faint-end UV luminosity function (LF). To set firmer constraints on the emissivity of AGN at $z<6$, we here make use of complete X-ray selected samples including deep Chandra and new COSMOS data, capable to efficiently measure the 1 ryd comoving AGN emissivity up to $z\sim5-6$ and down to five magnitudes fainter than probed by current optical surveys, without any luminosity extrapolation. We find good agreement between the logN$\rm{_H}\lesssim21-22$ cm$^{-2}$ X-ray LF and the optically-selected QSO LF at all redshifts for $M_{1450}\leq -23$. The full range of the logN$\rm{_H}\lesssim21-22$ cm$^{-2}$ LF ($M_{1450} \leq -17$) was then used to quantify the contribution of AGN to the photon budget critical value needed to keep the Universe ionized. We find that the contribution of ionizing AGN at $z = 6$ is as small as 1\% - 7\%, and very unlikely to be greater than 30\%, thus excluding an AGN-dominated reionization scenario.

The Giant Gemini GMOS survey of zem &gt; 4.4 quasars – I. Measuring the mean free path across cosmic time

We have obtained spectra of 163 quasars at $z_\mathrm{em}>4.4$ with the Gemini Multi Object Spectrometers on the Gemini North and South telescopes, the largest publicly available sample of high-quality, low-resolution spectra at these redshifts. From this homogeneous data set, we generated stacked quasar spectra in three redshift intervals at $z\sim 5$. We have modelled the flux below the rest-frame Lyman limit ($\lambda_\mathrm{r}<912$\AA) to assess the mean free path $\lambda_\mathrm{mfp}^{912}$ of the intergalactic medium to HI-ionizing radiation. At mean redshifts $z_\mathrm{q}=4.56$, 4.86 and 5.16, we measure $\lambda_\mathrm{mfp}^{912}=(22.2\pm 2.3, 15.1\pm 1.8, 10.3\pm 1.6)h_{70}^{-1}$ proper Mpc with uncertainties dominated by sample variance. Combining our results with $\lambda_\mathrm{mfp}^{912}$ measurements from lower redshifts, the data are well modelled by a simple power-law $\lambda_\mathrm{mfp}^{912}=A[(1+z)/5]^\eta$ with $A=(37\pm 2)h_{70}^{-1}$ Mpc and $\eta = -5.4\pm 0.4$ between $z=2.3$ and $z=5.5$. This rapid evolution requires a physical mechanism -- beyond cosmological expansion -- which reduces the cosmic effective Lyman limit opacity. We speculate that the majority of HI Lyman limit opacity manifests in gas outside galactic dark matter haloes, tracing large-scale structures (e.g. filaments) whose average density (and consequently neutral fraction) decreases with cosmic time. Our measurements of the strongly redshift-dependent mean free path shortly after the completion of HI reionization serve as a valuable boundary condition for numerical models thereof. Having measured $\lambda_\mathrm{mfp}^{912}\approx 10$ Mpc at $z=5.2$, we confirm that the intergalactic medium is highly ionized by that epoch and that the redshift evolution of the mean free path does not show a break that would indicate a recent end to HI reionization.