No Arabic abstract
The unprecedentedly bright afterglow of Swift GRB 130606A at z = 5.91 gave us a unique opportunity to probe the reionization era by high precision analyses of the redward damping wing of Ly alpha absorption, but the reported constraints on the neutral hydrogen fraction (f_HI) in intergalactic medium (IGM) derived from spectra taken by different telescopes are in contradiction. Here we examine the origin of this discrepancy by analyzing the spectrum taken by VLT with our own analysis code previously used to fit the Subaru spectrum. Though the VLT team reported no evidence for IGM HI using the VLT spectrum, we confirmed our previous result of preferring non-zero IGM HI (the best-fit f_HI ~ 0.06, when IGM HI extends to the GRB redshift). The fit residuals of the VLT spectrum by the model without IGM HI show the same systematic trend as the Subaru spectrum. We consider that the likely origin of the discrepancy between the two teams is the difference of the wavelength ranges adopted in the fittings; our wavelength range is wider than that of the VLT team, and also we avoided the shortest wavelength range of deep Ly alpha absorption (lambda_obs < 8426 A), because this region is dominated by HI in the host galaxy and the systematic uncertainty about host HI velocity distribution is large. We also study the sensitivity of these results to the adopted Ly alpha cross section formulae, ranging from the classical Lorentzian function to the most recent one taking into account fully quantum mechanical scattering. It is found that the preference for non-zero IGM HI is robust against the choice of the cross section formulae, but it is quantitatively not negligible and hence one should be careful in future analyses.
The unprecedentedly bright optical afterglow of GRB 130606A located by Swift at a redshift close to the reionization era (z = 5.913) provides a new opportunity to probe the ionization status of intergalactic medium (IGM). Here we present an analysis of the red Ly alpha damping wing of the afterglow spectrum taken by Subaru/FOCAS during 10.4-13.2 hr after the burst. We find that the minimal model including only the baseline power-law and HI absorption in the host galaxy does not give a good fit, leaving residuals showing concave curvature in 8400-8900 A with an amplitude of about 0.6% of the flux. Such a curvature in the short wavelength range cannot be explained either by extinction at the host with standard extinction curves, intrinsic curvature of afterglow spectra, or by the known systematic uncertainties in the observed spectrum. The red damping wing by intervening HI gas outside the host can reduce the residual by about 3 sigma statistical significance. We find that a damped Ly alpha system is not favored as the origin of this intervening HI absorption, from the observed Ly beta and metal absorption features. Therefore absorption by diffuse IGM remains as a plausible explanation. A fit by a simple uniform IGM model requires HI neutral fraction of f_HI ~ 0.1-0.5 depending on the distance to the GRB host, implying high f_HI IGM associated with the observed dark Gunn-Peterson (GP) troughs. This gives a new evidence that the reionization is not yet complete at z ~ 6.
We present the discovery and subsequent spectroscopy with Gemini-North of the optical afterglow of the Swift gamma-ray burst (GRB) 140515A. The spectrum exhibits a well-detected continuum at wavelengths longer than 8915 Angs with a steep decrement to zero flux blueward of 8910 Angs due to Ly-alpha absorption at redshift z~6.33. Some transmission through the Lyman-alpha forest is present at 5.2<z<5.733, but none is detected at higher redshift, consistent with previous measurements from quasars and GRB 130606A. We model the red damping wing of Lyman-alpha in three ways that provide equally good fits to the data: (a) a single host galaxy absorber at z=6.327 with log(N_HI)=18.62+/-0.08; (b) pure intergalactic medium (IGM) absorption from z=6.0 to z=6.328 with a constant neutral hydrogen fraction of x_HI=0.056+0.011-0.027; and (c) a hybrid model with a host absorber located within an ionized bubble of radius 10 comoving Mpc in an IGM with x_HI=0.12+/-0.05 (x_HI<0.21 at the 2-sigma level). Regardless of the model, the sharpness of the dropoff in transmission is inconsistent with a substantial neutral fraction in the IGM at this redshift. No narrow absorption lines from the host galaxy are detected, indicating a host metallicity of [Z/H]<~ -0.8. Even if we assume that all of the hydrogen absorption is due to the host galaxy, the column is unusually low for a GRB sightline, similar to two out of the other three highest-redshift bursts with measured log(N_HI). This is possible evidence that the escape fraction of ionizing photons from normal star-forming galaxies increases at z>~6.
One well-known way to constrain the hydrogen neutral fraction, x_H, of the high-redshift intergalactic medium (IGM) is through the shape of the red damping wing of the Lya absorption line. We examine this methods effectiveness in light of recent models showing that the IGM neutral fraction is highly inhomogeneous on large scales during reionization. Using both analytic models and semi-numeric simulations, we show that the picket-fence absorption typical in reionization models introduces both scatter and a systematic bias to the measurement of x_H. In particular, we show that simple fits to the damping wing tend to overestimate the true neutral fraction in a partially ionized universe, with a fractional error of ~ 30% near the middle of reionization. This bias is generic to any inhomogeneous model. However, the bias is reduced and can even underestimate x_H if the observational sample only probes a subset of the entire halo population, such as quasars with large HII regions. We also find that the damping wing absorption profile is generally steeper than one would naively expect in a homogeneously ionized universe. The profile steepens and the sightline-to-sightline scatter increases as reionization progresses. Of course, the bias and scatter also depend on x_H and so can, at least in principle, be used to constrain it. Damping wing constraints must therefore be interpreted by comparison to theoretical models of inhomogeneous reionization.
The observed deficit of strongly Lyman-alpha emitting galaxies at z>6.5 is attributed to either increasing neutral hydrogen in the intergalactic medium (IGM) and/or to the evolving galaxy properties. To investigate this, we have performed very deep near-IR spectroscopy of z>7 galaxies using MOSFIRE on the Keck-I Telescope. We measure the Lyman-alpha fraction at z~8 (combined photometric redshift peak at z=7.7) using two methods. First, we derived NLy{alpha}/Ntot directly using extensive simulations to correct for incompleteness. Second, we used a Bayesian formalism (introduced by Treu et al. 2012) that compares the z>7 galaxy spectra to models of the Lyman-alpha equivalent width (WLy{alpha}) distribution at z~6. We explored two simple evolutionary scenarios: smooth evolution where Lyman-alpha is attenuated in all galaxies by a constant factor (perhaps owing to processes from galaxy evolution or a slowly increasing IGM opacity), and patchy evolution where Lyman-alpha is blocked in some fraction of galaxies (perhaps due to the IGM being opaque along only some fraction of sightlines). The Bayesian formalism places stronger constraints compared with the direct method. Combining our data with that in the literature we find that at z~8 the Lyman-alpha fraction has dropped by a factor >3(84% confidence interval) using both the smooth and patchy scenarios compared to the z~6 values. Furthermore, we find a tentative evidence that the data favor the patchy scenario over smooth (with positive Bayesian evidence), extending trends observed at z~7 to higher redshift. If this decrease is a result of reionization as predicted by theory, then our data imply the volume averaged neutral hydrogen fraction in the IGM to be >0.3 suggesting that the reionization of the universe is in progress at z~8.
Events such as GRB130606A at z=5.91, offer an exciting new window into pre-galactic metal enrichment in these very high redshift host galaxies. We study the environment and host galaxy of GRB 130606A, a high-z event, in the context of a high redshift population of GRBs. We have obtained multiwavelength observations from radio to gamma-ray, concentrating particularly on the X-ray evolution as well as the optical photometric and spectroscopic data analysis. With an initial Lorentz bulk factor in the range Gamma_0 ~ 65-220, the X-ray afterglow evolution can be explained by a time-dependent photoionization of the local circumburst medium, within a compact and dense environment. The host galaxy is a sub-DLA (log N (HI) = 19.85+/-0.15), with a metallicity content in the range from ~1/7 to ~1/60 of solar. Highly ionized species (N V and Si IV) are also detected. This is the second highest redshift burst with a measured GRB-DLA metallicity and only the third GRB absorber with sub-DLA HI column density. GRB lighthouses therefore offer enormous potential as backlighting sources to probe the ionization and metal enrichment state of the IGM at very high redshifts for the chemical signature of the first generation of stars.