No Arabic abstract
We report new measurements of the intergalactic medium (IGM) Ly$alpha$ and Ly$beta$ effective optical depth at $5.3<z<6.5$, using a new sample of quasar sightlines including 32 quasars at $6.308le zle7.00$. These quasars provide a large statistical sample to measure the IGM evolution during the transition phase of the reionization epoch. We construct a data set of deep optical spectra of these quasars using VLT, Keck, Gemini, LBT, and MMT. We measure the Ly$alpha$ effective optical depth at $5.36<z<6.57$ using the Ly$alpha$ forests of both individual spectra and the stacked spectrum. The large scatter of individual measurements is consistent with previous work, suggesting an inhomogeneous reionization process. Combining our new measurements and previous results, we obtain a best-fit for the Ly$alpha$ effective optical depth evolution at $z>5.3$, $taupropto(1+z)^{8.6pm1.0}$. We then estimate the observed Ly$beta$ effective optical depth using Ly$beta$ forests and convert them to Ly$alpha$ optical depth for comparison, which provides additional constraints on the evolution of the IGM optical depth. The Ly$beta$-based measurements are generally in agreement with the best-fit evolution obtained from Ly$alpha$ forests. Using this new sample, we identify 389 Ly$alpha$ and 50 Ly$beta$ transmission spikes at $5.5<z<6.3$. The upper limits of Ly$alpha$ optical depth estimated using transmission spikes are well consistent with our best-fit evolution. The evolution in number density of these high-redshift transmission spikes suggests a rapid transition phase at the end of the reionization. Comparison of our optical depth measurements with hydrodynamical simulations indicates a IGM neutral hydrogen fraction $langle f_{rm HI}ranglegtrsim10^{-4}$ at $z=6$.
Luminous $zge7$ 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$alpha$ 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 towards that line of sight. However, damping wing features have only been detected in two $z>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/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $rm (1.39pm0.16) times10^{9} ~M_odot$ SMBH accreting at an Eddington ratio of $lambda_{rm Edd}=0.7pm0.1$, consistent with the values seen in other luminous $zsim 7$ quasars. Furthermore, the Ly$alpha$ 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 $langle x_{rm HI} rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3le t_{rm Q}le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.
We test the galactic outflow model by probing associated galaxies of four strong intergalactic CIV absorbers at $z=5$--6 using the Hubble Space Telescope (HST) ACS ramp narrowband filters. The four strong CIV absorbers reside at $z=5.74$, $5.52$, $4.95$, and $4.87$, with column densities ranging from $N_{rm{CIV}}=10^{13.8}$ cm$^{-2}$ to $10^{14.8}$ cm$^{-2}$. At $z=5.74$, we detect an i-dropout Ly$alpha$ emitter (LAE) candidate with a projected impact parameter of 42 physical kpc from the CIV absorber. This LAE candidate has a Ly$alpha$-based star formation rate (SFR$_{rm{Lyalpha}}$) of 2 $M_odot$ yr$^{-1}$ and a UV-based SFR of 4 $M_odot$ yr$^{-1}$. Although we cannot completely rule out that this $i$-dropout emitter may be an [OII] interloper, its measured properties are consistent with the CIV powering galaxy at $z=5.74$. For CIV absorbers at $z=4.95$ and $z=4.87$, although we detect two LAE candidates with impact parameters of 160 kpc and 200 kpc, such distances are larger than that predicted from the simulations. Therefore we treat them as non-detections. For the system at $z=5.52$, we do not detect LAE candidates, placing a 3-$sigma$ upper limit of SFR$_{rm{Lyalpha}}approx 1.5 M_odot$ yr$^{-1}$. In summary, in these four cases, we only detect one plausible CIV source at $z=5.74$. Combining the modest SFR of the one detection and the three non-detections, our HST observations strongly support that smaller galaxies (SFR$_{rm{Lyalpha}} lesssim 2 M_odot$ yr$^{-1}$) are main sources of intergalactic CIV absorbers, and such small galaxies play a major role in the metal enrichment of the intergalactic medium at $zgtrsim5$.
We compare a sample of five high-resolution, high S/N Ly$alpha$ forest spectra of bright $6<z lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 pm 1600,10500pm 2100,12000 pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $zsim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.
We explore the kinematics of 27 z~6 quasar host galaxies observed in [CII]-158 micron ([CII]) emission with the Atacama Large Millimeter/sub-millimeter Array at a resolution of ~0.25. We find that nine of the galaxies show disturbed [CII] emission, either due to a close companion galaxy or recent merger. Ten galaxies have smooth velocity gradients consistent with the emission arising from a gaseous disk. The remaining eight quasar host galaxies show no velocity gradient, suggesting that the gas in these systems is dispersion-dominated. All galaxies show high velocity dispersions with a mean of 129+-10 km/s. To provide an estimate of the dynamical mass within twice the half-light radius of the quasar host galaxy, we model the kinematics of the [CII] emission line using our publicly available kinematic fitting code, qubefit. This results in a mean dynamical mass of 5.0+-0.8(+-3.5) x 10^10 Msun. Comparison between the dynamical mass and the mass of the supermassive black hole reveals that the sample falls above the locally derived bulge mass--black hole mass relation at 2.4sigma significance. This result is robust even if we account for the large systematic uncertainties. Using several different estimators for the molecular mass, we estimate a gas mass fraction of >10%, indicating gas makes up a large fraction of the baryonic mass of z~6 quasar host galaxies. Finally, we speculate that the large variety in [CII] kinematics is an indication that gas accretion onto z~6 super massive black holes is not caused by a single precipitating factor.
We present the discovery of one or two extremely faint z~6 quasars in 6.5 deg^2 utilizing a unique capability of the wide-field imaging of the Subaru/Suprime-Cam. The quasar selection was made in (i-z_B) and (z_B-z_R) colors, where z_B and z_R are bandpasses with central wavelengths of 8842A and 9841A, respectively. The color selection can effectively isolate quasars at z~6 from M/L/T dwarfs without the J-band photometry down to z_R<24.0, which is 3.5 mag. deeper than SDSS. We have selected 17 promising quasar candidates. The follow-up spectroscopy for seven targets identified one apparent quasar at z=6.156 with M_1450=-23.10. We also identified one possible quasar at z=6.041 with a faint continuum of M_1450=-22.58 and a narrow Lyman-alpha emission with HWHM=427 km/s, which cannot be distinguished from Lyman-alpha emitters. We derive the quasar luminosity function at z~6 by combining our faint quasar sample with the bright quasar samples by SDSS and CFHQS. Including our data points invokes a higher number density in the faintest bin of the quasar luminosity function than the previous estimate employed. This suggests a steeper faint-end slope than lower-z, though it is yet uncertain based on a small number of spectroscopically identified faint quasars and several quasar candidates are still remain to be diagnosed. The steepening of the quasar luminosity function at the faint-end does increase the expected emission rate of the ionizing photon, however, it only changes by a factor of ~2-6. This was found to be still insufficient for the required photon budget of reionization at z~6.