We present an absorption-line survey of optically thick gas clouds -- Lyman Limit Systems (LLSs) -- observed at high dispersion with spectrometers on the Keck and Magellan telescopes. We measure column densities of neutral hydrogen NHI and associated
metal-line transitions for 157 LLSs at z=1.76-4.39 restricted to 10^17.3 < NHI < 10^20.3. An empirical analysis of ionic ratios indicates an increasing ionization state of the gas with decreasing NHI and that the majority of LLSs are highly ionized, confirming previous expectations. The Si^+/H^0 ratio spans nearly four orders-of-magnitude, implying a large dispersion in the gas metallicity. Fewer than 5% of these LLSs have no positive detection of a metal transition; by z~3, nearly all gas that is dense enough to exhibit a very high Lyman limit opacity has previously been polluted by heavy elements. We add new measurements to the small subset of LLS (~5-10) that may have super-solar abundances. High Si^+/Fe^+ ratios suggest an alpha-enhanced medium whereas the Si^+/C^+ ratios do not exhibit the super-solar enhancement inferred previously for the Lya forest.
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 $zsim 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.2pm 2.3, 15.1pm 1.8, 10.3pm 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=(37pm 2)h_{70}^{-1}$ Mpc and $eta = -5.4pm 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.
We present the first science results from our Hubble Space Telescope Survey for Lyman limit absorption systems (LLS) using the low dispersion spectroscopic modes of the Advanced Camera for Surveys and the Wide Field Camera 3. Through an analysis of 7
1 quasars, we determine the incidence frequency of LLS per unit redshift and per unit path length, l(z) and l(x) respectively, over the redshift range 1 < z< 2.6, and find a weighted mean of l(x)=0.29 +/-0.05 for 2.0 < z < 2.5 through a joint analysis of our sample and that of Ribaudo et al. (2011). Through stacked spectrum analysis, we determine a median (mean) value of the mean free path to ionizing radiation at z=2.4 of lambda_mfp = 243(252)h^(-1) Mpc, with an error on the mean value of +/- 43h^(-1) Mpc. We also re-evaluate the estimates of lambda_mfp from Prochaska et al. (2009) and place constraints on the evolution of lambda_mfp with redshift, including an estimate of the breakthrough redshift of z = 1.6. Consistent with results at higher z, we find that a significant fraction of the opacity for absorption of ionizing photons comes from systems with N_HI <= 10^{17.5} cm^(-2) with a value for the total Lyman opacity of tau_lyman = 0.40 +/- 0.15. Finally, we determine that at minimum, a 5-parameter (4 power-law) model is needed to describe the column density distribution function f(N_HI, X) at z sim 2.4, find that f(N_HI,X) undergoes no significant change in shape between z sim 2.4 and z sim 3.7, and provide our best fit model for f(N_HI,X).
