The Mean UV Spectrum of z\sim 4 LBGs

# [

Tucker Jones11affiliation: Astronomy Department, California Institute of Technology, MC249-17, Pasadena, CA 91125, USA , Daniel P. Stark22affiliation: Institute of Astronomy, Cambridge CB3 0HA, UK 33affiliation: Steward Observatory, University of Arizona, Tucson, AZ 85721 44affiliation: Hubble Fellow , Richard S. Ellis11affiliation: Astronomy Department, California Institute of Technology, MC249-17, Pasadena, CA 91125, USA
###### Abstract

We present and discuss the mean rest-frame ultraviolet spectrum for a sample of 81 Lyman Break Galaxies (LBGs) selected to be B-band dropouts with a mean redshift of and apparent magnitudes . Most of the individual spectra are drawn from our ongoing survey in the GOODS fields with the Keck DEIMOS spectrograph described in earlier papers in the series, and we have augmented our sample with published data taken with FORS2 on the VLT. In general we find similar trends in the spectral diagnostics to those found in the earlier, more extensive survey of LBGs at undertaken by Shapley et al. (2003). Specifically, we find low-ionization absorption lines which trace the presence of neutral outflowing gas are weaker in galaxies with stronger Ly emission, bluer UV spectral slopes, lower stellar masses, lower UV luminosities, and smaller half-light radii. This is consistent with a physical picture whereby star formation drives outflows of neutral gas which scatters Ly and gives rise to strong low-ionization absorption lines, while increasing the stellar mass, size, metallicity, and dust content of galaxies. Typical galaxies are thus expected to have stronger Ly emission and weaker low-ionization absorption at earlier times (higher redshifts). Indeed, our mean spectrum at shows somewhat weaker low-ionization absorption lines than at and available data at high redshift demonstrate that this evolutionary trend continues. Although the total absorption by low-ionization transitions weakens at high redshift, the fine structure emission lines are stronger suggesting a greater concentration of neutral gas at small galactocentric radius ( kpc). In conjunction with earlier results from our spectroscopic survey which demonstrated an increased fraction of LBGs with Ly emission at higher redshift, we argue that the reduced low-ionization absorption is likely caused by a decrease in the covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. At present we cannot distinguish between differences in the covering fraction, outflow kinematics and geometry as the underlying cause of this interesting evolution. However, our continuing survey will enable us to extend these diagnostics more reliably to higher redshift and determine the implications for the escape fraction of ionizing photons which governs the role of early galaxies in cosmic reionization.

galaxies: high redshift – galaxies: evolution

The Mean UV Spectrum of LBGs]Keck Spectroscopy of Faint Lyman Break Galaxies: III. The Mean Ultraviolet Spectrum at

## 1. Introduction

Considerable progress has been made over the past decade in charting the demographics of high redshift galaxies. Multi-wavelength surveys have defined the luminosity functions of UV and sub-mm selected star-forming sources (Reddy & Steidel, 2009; Wardlow et al., 2011) as well as the coeval population of quiescent massive red galaxies (Brammer et al., 2011). Spitzer data has revealed the time-dependent stellar mass density - a complementary quantity which represents the integral of the past star formation activity (e.g. Stark et al. 2009). Through these surveys a well-defined picture of the history of star formation and mass assembly over has been empirically determined (Hopkins & Beacom, 2006; Ellis, 2008; Robertson et al., 2010). The redshift range corresponds to the peak of star formation activity where the Hubble sequence starts to emerge, and the earlier era corresponding to is an even more formative one where mass assembly was particularly rapid.

Intermediate dispersion spectroscopy of carefully-selected Lyman break galaxies (LBGs) has been particularly important in defining population trends that cannot be identified from photometric data alone. A very influential study at was undertaken by Shapley et al. (2003) who used composite Keck LRIS spectra of various subsets of nearly 1000 LBGs to examine the role of hot stars, Hii regions and dust obscuration, as well as to measure the outflow kinematics and absorption line properties of neutral and ionized gas. Composite spectra are particularly useful for measuring weak lines which cannot be studied in detail for individual objects. Through these careful studies, a detailed picture of the mass-dependent evolution of LBGs has emerged (see Shapley 2011 for a recent review).

In earlier papers in this series (Stark et al. 2010, hereafter Paper I; Stark et al. 2011, hereafter Paper II), we introduced an equivalent spectroscopic survey of LBGs selected from a photometric catalog of more distant LBGs with in the Great Observatories Origins Deep Survey (GOODS) fields (Giavalisco et al., 2004; Stark et al., 2009). Whereas the Shapley et al. (2003) study targetted the study of LBGs close to the peak of activity in the overall cosmic star formation history, this earlier period corresponds to a less well-studied era when the rate of mass assembly is particularly rapid. From photometric data alone, Stark et al. (2009) deduced some significant changes in the characteristics of star formation at compared to later times, for example a shorter timescale of activity ( Myr). We considered it crucial to understand these changes in LBG properties if these galaxies are to be used as probes of cosmic reionization at higher redshifts.

At the time of writing, our Keck survey is continuing with increasing emphasis at high redshift. Paper I presented the first substantial results from a survey of LBGs at observed with the Keck/DEIMOS spectrograph. Paper II augmented this data with a further sample following more ambitious exposures focusing primarily on LBGs. Incorporating a sample of ESO VLT spectra, retrospectively selected using similar photometric criteria as those for the Keck sample from the FORS2 study of Vanzella et al. (2005, 2006, 2008, 2009), the current dataset amounts to a sample of 546 galaxies over the redshift range .

Our earlier papers in this series concentrated primarily on the rate of occurrence of Lyman (Ly ) emission in our spectra (the “Ly fraction”). The overall goal was to understand the significantly different evolutionary trends in the luminosity functions of LBGs and narrow-band selected Ly emitters (LAEs, Ouchi et al. 2008) prior to the use of the Ly fraction as a test of when reionization ended (Schenker et al., 2011). Paper I confirmed a result found by Shapley et al. (2003) at lower redshift, namely that Ly emission is more frequent in lower luminosity LBGs rising to a high proportion % at . More importantly, the Ly fraction was found to rise modestly with redshift over . By correlating the visibility of Ly emission with UV continuum slopes derived from the HST photometry, it was argued that these trends in the visibility of line emission most probably arise from different amounts of dust obscuration. Reduced dust extinction in lower luminosity LBGs and those at higher redshift has also been deduced from studies of larger photometric samples (Bouwens et al. 2009, see also Reddy & Steidel 2009).

Paper I also discussed the possibility that the covering fraction of hydrogen may be lower in low-luminosity LBGs. Strong Ly emission in luminous LBGs is often associated with low equivalent width interstellar absorption lines arising from a non-uniform covering fraction of neutral hydrogen (Quider et al., 2009; Shapley et al., 2003). This trend suggests that the high Ly fraction in faint LBGs is partially due to a lower covering fraction, which would imply that Lyman continuum photons may more easily escape from intrinsically faint galaxies. Such a result would have great importance in understanding the role of galaxies in maintaining cosmic reionization (Robertson et al., 2010).

The present paper represents our first analysis of the spectral properties of LBGs derived from composite spectra in the manner pioneered at by Shapley et al. (2003). The large database now amassed following the campaigns at Keck and the VLT makes a similar study now practical in the redshift range where there is evidence of increased short-term star formation and the mass assembly rate is particularly rapid. Via detailed studies of low-ionization absorption line and emission line profiles, we aim to examine possible changes in the kinematics and covering fraction of neutral gas, which affects the strength of Ly and the escape fraction of ionizing photons, as we approach the reionization era. As our redshift survey continues, in this paper we focus on a sample of galaxies with , selected as B-dropout LBGs. Combining our Keck sample with data from the VLT (see Paper I for details), herein we examine the spectral features and trends in composite spectra drawn from a sample of 131 galaxies.

A plan of the paper follows. We briefly review the spectroscopic observations and their data reduction in §2; much of the relevant discussion is contained in Paper I. In §3 we describe the selection of individual spectra that we consider appropriate for forming the composite mean spectrum at and the associated sources of uncertainty. §4 examines the mean spectrum in detail and introduces the various diagnostic features in the context of a physical model for LBGs of different masses and star-formation rates. In §5 we compare spectroscopic trends grouped by observable properties such as mass and luminosity with those found at by Shapley et al. (2003). In §6 we discuss those trends which appear to be redshift-dependent, discussing implications for the role of early star-forming galaxies in cosmic reionization. Finally we summarize our results in §7.

Throughout this paper, we adopt a flat CDM cosmology with , , and km sMpc. All magnitudes in this paper are quoted in the AB system (Oke, 1974).

## 2. Observations and Data Reduction

The rationale and procedures used to undertake our spectroscopic survey of LBGs over were introduced in detail in Paper I and the interested reader is referred to that paper for further detail. Here we recount only the basic details. Our target LBGs were selected as , or –band “dropouts” based on deep photometry in the two GOODS fields. The photometric catalog used in our analysis will be described in detail in Stark et al. (2011, in preparation). The selection and photometric approach is largely similar to that described in Stark et al. (2009) and Paper I, but we highlight two key updates. First the selection is performed on the v2 GOODS ACS catalogs. Second, we utilize deep ground-based near-IR imaging in GOODS-N obtained from WIRCAM on CFHT (Wang et al., 2010) and deep HST Wide Field Camera 3 / IR imaging of GOODS-S from CANDELS (Grogin et al., 2011; Koekemoer et al., 2011).

### 2.1. Keck/DEIMOS

The majority of spectra we present are taken from an ongoing survey with the DEep Imaging Multi-Object Spectrograph (DEIMOS; Faber et al. 2003) on the Keck II telescope. In this paper we have used observations taken in 2008 April and 2009 March (masks GN081, GN082, GN083, GN094, and GN095 in Paper I), which targeted a total of 261 -drop and 88 -drop galaxies. These data were taken with the 600 lines mm grating, covering the wavelength range 4850–10150 Å  with a resolution of Å. We do not use data at wavelengths Å  which are affected by strong and variable absorption by atmospheric water vapor. The seeing was typically 08 and ranged between 05 and 10 during the observations.

All data were reduced and calibrated using a modified version of the IDL pipeline spec2d, developed specifically for DEIMOS by the DEEP2 survey team (Davis et al., 2003). The data were reduced as described in Paper I, with the addition of two important modifications. First, the continuum traces of all target galaxies and other objects occupying the same slit were carefully masked to exclude object flux from the sky background model. Second, the b-spline fit to the sky background was modified to include a 2nd-order polynomial fit to the spatial dimension. These modifications significantly improved the sky subtraction, particularly at long wavengths Å where bright sky lines can be problematic.

The reduced one- and two-dimensional spectra were visually inspected using the IDL program SpecPro (Masters & Capak, 2011). Spectra which suffered from poor data quality were excluded from further analysis. The exclusions included unacceptable amounts of scattered light within the detector mosaic, defective CCD columns, contamination from bright nearby sources, or poor sky subtraction, for example arising from the location of the slit on the detector. A few low-redshift interlopers were also identified and excluded from further analysis. For the remaining spectra, redshifts were measured from either Ly emission (where present) or interstellar absorption lines. Galaxies identified as hosting strong AGN based on the presence of Civ or other strong emission lines were excluded. We identified 5 stars, 2 dusty low-redshift galaxies (), 4 AGN, and 134 star-forming galaxies with secure redshifts . After rejecting poor data, the final sample consists of 94 high-quality spectra with accurate redshifts. Examples of high redshift dropout spectra are shown in Figure 1.

### 2.2. Archival VLT/FORS2 Spectroscopy in GOODS-S

To augment our sample of high redshift spectra, we have also made use of data from the FORS2 program of Vanzella et al. (2005, 2006, 2008, 2009) which targeted dropouts in the GOODS-S field. The characteristics of that survey in terms of resolution and spectral coverage are very similar to that undertaken at Keck with DEIMOS and details can be found in Paper I. Using the coordinates provided in the published FORS2 database we queried the version 2.0 ACS catalogs for GOODS-S and undertook our own photometric measures and dropout selection criteria in an identical fashion to that used for our Keck survey. The magnitude distribution of the FORS2 sample is generally weighted towards sources brighter than those in the overall Keck survey, but for the purposes of constructing the mean spectra discussed in this paper, the bulk of the individual spectra are of comparable brightness.

### 2.3. Redshift measurements

Care is needed in deriving accurate systemic redshifts from rest-frame UV spectra since the strongest features trace the kinematics of outflowing gas rather than that of the stars. Stellar absorption lines are usually too faint to be measured precisely given the signal to noise of the spectra. Typically the only features detected in individual spectra are Ly and strong interstellar absorption lines such as Siii1260, Oi1302+Siii1304, Cii1334, Siiv1393,1402, Siii1526, and Civ 1548,1550 (Figure 1). Absorption by outflowing gas results in a blueshift of interstellar absorption lines. Outflowing neutral hydrogen along the line of sight leads to a Ly profile which displays very broad ( km s) blueshifted absorption and net redshifted emission. The high-ionization Siiv and Civ lines also arise in P-Cygni stellar winds with broad blueshifted absorption. The magnitude of these offsets has been well-quantified for star-forming galaxies at , and is typically km sfor interstellar absorption lines and km sfor Ly emission (Steidel et al., 2010).

To determine accurate redshifts for making composite spectra, we restrict our sample to those with redshifts measured from either Ly emission () or low-ionization interstellar absorption lines (). Although, as discussed above, these are not at the systemic redshift, we follow well-established techniques to correct for the typical offsets (Steidel et al., 2010). We do not use redshifts based on Ly absorption or high-ionization lines (Siiv and Civ) because of the complex and variable blueshifts of these features with respect to the systemic redshift. We define as the centroid of the emission line and consider only spectra in which the line is detected at significance. Interstellar absorption line redshifts require careful treatment to avoid spurious identification of sky line residuals. We consider only the low-ionization Siii1260, Oi1302+Siii1304, and Cii1334 features which are typically found in the highest signal-to-noise regions of our spectra. Absorption features at longer wavelengths are less reliable due to the higher density of strong night sky lines, while shorter wavelength transitions are lost in the Ly forest. To measure absorption line redshifts, we first estimate the redshift from Ly (either in emission or absorption). We then fit Gaussian profiles to the spectrum near the expected position of the three features above. We require the best fit of all three lines to be consistent to within km swith a combined significance of . If these conditions are met, we define as the weighted mean redshift of the three interstellar features.

With the criteria above, suitable redshifts are available for a total of 131 high-quality DEIMOS and FORS2 spectra. 91 redshifts are based on measures of emission only, 31 have only, and 9 have both and . The distribution of redshifts and absolute UV magnitudes for this sample is shown in Figure 2. Unless stated otherwise, further analysis in this paper is restricted to the 81 sources with redshift and apparent magnitude . We applied this additional magnitude criterion in order to ensure a well-defined continuum signal/noise in each individual spectrum. 53 of these 81 sources are drawn from the Keck survey and 28 from FORS2.

### 2.4. Sample bias

While the sample is 90% complete for dropout selected galaxies to an apparent magnitude , at fainter magnitudes Ly emission or strong interstellar absorption features are required for reliable redshifts. This results in a bias towards stronger low-ionization absorption lines for the fraction (%) without Ly . We can quantify this bias through the detectability of the average low-ionization absorption line strength shown in the composite spectrum in Figure 3 (see §3). We find that we can measure redshifts for this average absorption line strength with 90% completeness at , declining through 50% at to zero at . Of the galaxies with detected absorption line redshifts, those with have line strengths only 10% stronger than for those with . Clearly this is a small effect.

The bias toward stronger Ly emission for fainter galaxies is manifest in Figure 2 where a paucity of objects with Å can be seen at faint magnitudes. This bias was fully quantified in Paper I using Monte Carlo simulations. For galaxies at , a sample completeness of 95% is reached at Å  for , and Å  for . For the sample presented in this paper, the least biased subset is that with strong Ly emission, followed by that with bright apparent magnitudes .

A final issue in considering composite spectra is that these are comprised of individual spectra across with different rest-frame wavelength ranges. As our spectra generally cover the wavelength range Å, galaxies at contribute to the rest frame Å, while those at contribute to Å. At longer rest-frame wavelengths the composite spectrum is therefore largely contributed by galaxies at lower redshift. In addition there are wavelength-dependent sources of noise discussed in §3.2. In summary, the data used to construct the composite spectrum in Figure 3 correspond to a mean redshift at Å, at 1200 Å, at 1500 Å, and at 1650 Å. For the wavelength range of interest in this work, the redshift bias is not particularly troublesome.

## 3. Composite Spectra

We now turn to the presentation of the composite spectra. We need to account for the difference between the redshifts determined using Ly on the one hand and the low-ionization interstellar lines on the other hand and the systemic redshift prior to co-addition. We also seek to understand the signal/noise of the composite in terms of the statistical uncertainties and the variance among the individual spectra used to construct the composite.

Composite spectra are constructed by shifting the individual spectra into the rest frame according to a deduced systemic redshift and then averaging the set. In general terms we will first identify a subsample based on their observable properties. Each spectrum in the sample is shifted to the adopted rest frame and interpolated to a common wavelength scale with a dispersion of 0.12 Å. All spectra are normalized to have a median in the range 1250–1500 Å. Spectra taken with DEIMOS are smoothed to a resolution of 1.9 Åto match the lower resolution of FORS2 data. The spectra are then averaged at each wavelength using a -clipped mean to reject outliers arising from sky subtraction residuals and cosmetic defects. An equal number of positive and negative outliers are rejected at each wavelength, totaling at most 30% of the data. The remaining data are averaged with an arithmetic mean.

Uncertainty in a composite spectrum will arise from both the finite signal to noise and the variance of the individual galaxies. For example, the variance in Ly equivalent widths in our sample is much greater than the uncertainty measured in the individual spectra. It is especially important to quantify the sample variance for weak features that are generally not detected in individual spectra. We account for sample variance with a bootstrap technique. For each composite spectrum we create 100 alternate composites using the same number of spectra but drawn at random from the parent sample. Each alternate has an average 63% of the sample represented with 37% duplicates. Every measurement made on the composite spectrum is repeated for each of the 100 alternates. We then take the measurement error to be the standard deviation of the 100 alternate measurements, which reflects both the sample variance and finite signal to noise.

As discussed in §2.3 the most challenging issue is to determine the systemic redshift prior to shifting to the rest-frame. Here we follow the approach used by Shapley et al. (2003). As a first approximation we use the value of (where available) to construct a composite spectrum. This enables us to locate the stellar photospheric line Ciii1176 in the composite where we detect a velocity difference of km s  with respect to Ly . We can thus infer that Ly emission in our sample is redshifted on average by km s. In a similar fashion, stacking spectra using the redshift based on low-ionization interstellar absorption results in a detection of Ciii1176 with a velocity offset of km s. For comparison, at 3 Steidel et al (2010) find km s  and km s. To construct composite spectra, we use either shifted by km sor shifted by km sto approximate the systemic redshift of each galaxy. We use the Ly -based redshift when available since it is typically determined with greater precision than . Figure 3 shows the composite spectrum of 81 galaxies in our sample with and apparent magnitude using this method.

### 3.1. Uncertainties in the Systemic Redshift

A natural concern is the extent to which these applied shifts might vary within the sample used to make the composite. This can be estimated from observations of higher signal to noise from spectra taken at lower redshift. Steidel et al. (2010) quantify the offset between , , and the systemic redshift in a sample of 89 galaxies at . They find velocity offsets km s  and km srelative to H. Assuming this is representative of our data, the uncertainty in our systemic redshift is therefore likely to be km s. An upper limit on can be estimated from the width of spectral lines in the composite spectrum. In particular the rest-wavelength of the stellar line Ciii1176 in the composite provides a valuable measurement of the average offset from the systemic velocity and its width provides an upper limit on the effective spectral resolution. We measure a systemic velocity of km s  in the composite spectrum (Figure 3) and a FWHM km s  (deconvolved from the instrumental resolution km s). The uncertainty in the adopted redshifts about the true systemic stellar value is therefore km s  FWHM or equivalently km s, comparable in fact to what was achieved for individual spectra at lower redshift by Steidel et al. (2010).

### 3.2. Error spectrum

Figure 3 shows the composite spectrum of our sample as well as the error spectrum derived using the bootstrap technique discussed above. The error at each pixel is calculated as the standard deviation of all averaged data points (excluding outliers), divided by the square root of the number of data points. There are several wavelength-dependent factors contributing to the error spectrum in addition to the finite signal-to-noise of individual spectra. One factor is the intrinsic sample variance, seen clearly as a noise spike at the position of Ly , and also evident for Cii1334 and other absorption features. Another is the decreased instrument throughput at Å, leading to higher noise at shorter wavelengths. Similarly, stronger sky line emission causes increasing noise at longer wavelengths. Finally, the number of contributing spectra peaks at rest-frame Å , with increased noise at higher and lower wavelengths where fewer spectra are included. Our observed wavelength range Å  corresponds to a rest-frame Å  at the mean redshift of our magnitude-limited sample. Ultimately we acheive a signal to noise ratio in the continuum of between the Ly line and rest-frame 1800 Å  and in the Ly forest, with a peak S/N at 1350 Å.

## 4. Features in the Composite Spectrum

We now discuss the composite spectrum (Figure 3) in more detail, focusing on the strong spectral features at 1215–1550 Å  where we have the best signal to noise. In this wavelength range we detect Ly , Nv1239,1243, Siii1260, Siii*1265, Oi1302+Siii1304 (blended), Siii*1309, Cii1334, Siiv1394,1403, Siii 1527, Siii*1533, and Civ1548,1550 at high significance. Our composite spectrum is very similar to the composite of LBGs presented in Shapley et al. (2003), which we show in Figure 4 for comparison. The absolute magnitude distribution of our sample is broadly similar to that at ; both cover the range . Our discussion below follows closely that presented originally by Shapley et al. (2003) but we are also interested in whether there are differences seen over the redshift range . We will discuss these possible evolutionary trends in §5.

### 4.1. Lyman Break Galaxies: a physical picture

It is helpful to begin by describing a possible physical picture of LBGs based on many analyses of the extensive observations at (see Shapley 2011 for a recent review). Typical * LBGs at have ultraviolet half-light radii kpc, stellar masses M, and star formation rates Myr(Bouwens et al., 2004; Ferguson et al., 2004; Shapley et al., 2001). The star formation surface density is sufficient to drive “superwinds” of outflowing gas, similar to those seen in local galaxies where Myrkpc (Heckman, 2002). Indeed, blueshifted interstellar absorption lines confirm there are outflows with typical velocities of km s, and in some cases as high as 800 km s(Shapley et al., 2003; Pettini et al., 2002; Quider et al., 2009, 2010). However, the physical origin of these outflows remains unclear (e.g. Murray et al. 2010). It is thought these outflows produce an extended circumgalactic medium (CGM) of ejected material. Outflowing gas is found in both low and high ionization states (e.g. Siii and Siiv). Low-ionization transitions such as Siii and Cii are mostly associated with neutral hydrogen whereas high ionization lines occur in the fully ionized component. Based on trends in the strength of the various species, Shapley et al. (2003) suggest a geometry in which discrete clouds of neutral gas are embedded in a halo of ionized gas. Steidel et al. (2010) show that the neutral and ionized CGM components both extend to radii of at least 125 kpc.

### 4.2. Lyα

Ly is the most prominent and diverse feature in our individual spectra. The line originates from hydrogen recombination in H II regions photoionized by massive stars and, as these stars dominate the adjacent continuum, its intrinsic equivalent width should be within the range Å  for nearly all stellar populations (e.g. Forero-Romero et al. 2011). Our composite spectrum reveals a complex line profile with strong absorption extending blueward to km s, and redshifted emission to km s(Figure 5) with a peak offset of km srelative to the systemic velocity. The net equivalent width Å  is only % of the expected intrinsic value.

The form of this Ly profile has been readily understood in terms of the physical picture discussed in §4.1. Ly emission produced at the systemic velocity can escape only along a line of sight free of neutral hydrogen or if it is shifted in velocity far from resonance. As photons escape they encounter blueshifted clouds of partially neutral gas, which absorb and re-emit isotropically. Photons backscattered from neutral clouds at small radii will appear redshifted and have a higher chance of escaping. Additionally, photons scattered from neutral gas at the edge of the CGM can escape into the ionized intergalactic medium and will be observed as blueshifted emission (Steidel et al., 2010). Indeed, Figure 5 shows such a weak blueshifted emission peak in the composite spectrum. On the other hand, the reduced equivalent width compared to that expected intrinsically could be due to many effects including dust extinction, scattering at large radii where emission falls outside the spectroscopic slits and a non-zero escape fraction of ionizing radiation. Since these effects cannot easily be disentangled, the absolute strength of the Ly line must be interpreted with caution (Paper II, Schenker et al. 2011).

We can estimate the spatial extent of the neutral CGM (i.e. the radius at which escaping Ly is last scattered) based on the strength of the diffuse blueshifted Ly emission. The spatial profile of extended Ly emission has been well-quantified at by Steidel et al. (2011) who find that diffuse haloes are a generic property of star-forming galaxies leading to a total Ly flux greater than that measured within the spectroscopic slits.

Average Ly surface brightness profiles are well-fit at large radii by an exponential form, viz.

 ΣLyα(b)=Σ0exp(b/bl) (1)

where kpc. The blueshifted emission shown in Figure 5 has an equivalent width Å  measured within a slit aperture of , corresponding to kpc at . Assuming this approximates the peak value, we estimate per 50 kpc and integrating Equation 1 yields a total equivalent width of . Using the measured value of , the scale length is given by

 bl=(2.3±0.2)√Whalo/\AAkpc. (2)

Although the total Ly flux is not measured, we can estimate its value from theoretical expectations as well as observations at lower redshift. The equivalent width in Figure 5 is Å  (excluding the blueshifted emission component), thus we can write Å. As discussed above, the observed is diminished by dust and the escape of ionizing radiation. Indeed, Steidel et al. (2011) measure Å  in various subsamples of their data, all lower than the expected intrinsic value, =100-200 Å . Considering the composite DEIMOS spectra shown in Figure 5 and assuming an intrinsic Å (Forero-Romero et al., 2011) we derive upper limits of kpc for galaxies with , and kpc for those with Å. This latter constraint is the most stringent, and also likely closest to the true value of . Steidel et al. (2011) find that their sample of Ly emitters (LAEs, defined as Å) has the highest ( Å) and largest scale length ( kpc) of any sub-sample that they analyze. Our constraint of kpc for the LAEs therefore suggests that the characteristic size of Ly haloes may be smaller at than at . Due to the large inherent uncertainties, direct measurements of low surface brightness Ly emission are needed to confirm this possibility.

### 4.3. Low-ionization metal transitions

According to the physical picture in §4.1, absorption in low-ionization transitions occurs in both the interstellar medium and outflowing clouds of cool gas. The former is at the systemic velocity while the latter is blueshifted. Such transitions are generally saturated (Shapley et al., 2003; Pettini et al., 2002) so the line depth at a given velocity provides a measure of the areal covering fraction of O and B stars by neutral gas along the line of sight. Specifically, the line profile is given by where is the continuum flux. While typical LBGs are too faint for detailed line profiles, high-resolution spectra of a few bright lensed sources at (Pettini et al., 2002; Quider et al., 2009, 2010) have revealed absorption velocities ranging from to km s  with the highest covering fraction at km s. The mean low-ionization absorption line velocity offset in our composite is km sin good agreement with that at (Shapley et al., 2003; Steidel et al., 2010).

Of particular diagnostic value are the two relatively unblended transitions of Siii at 1260 and 1527 Å whose equivalent width ratio is a valuable tracer of the optical depth (Shapley et al., 2003). The profiles of both lines and their ratio are shown in Figure 6. Gas near the systemic velocity ( km s) is clearly optically thick, while for km s  it is intermediate in optical depth. This difference is also seen in high-resolution spectra of lensed galaxies (Pettini et al., 2002; Quider et al., 2010), and optically thin gas is seen at large galactocentric radii (Steidel et al., 2010) suggesting that smaller, optically thin clouds are more easily accelerated to high velocity and large distances. We note that the column density at which both Siii transitions become optically thick () is cm. Assuming as measured for the lensed galaxy cB58 (Pettini et al., 2002), we estimate that optically thin absorption occurs in clouds with hydrogen column densities cm.

#### 4.3.1 Fine structure transitions

A satisfying aspect of our composite spectrum is the successful identification of fine structure emission lines of Siii (see marked features in Figure 3). Siii ions in the CGM absorb photons in the resonance transitions and immediately re-emit a photon at approximately the same velocity. A photon absorbed at 1260 Å  will be re-emitted at either the same wavelength or to the fine structure transition Siii*1265. Likewise a photon absorbed at 1527 Å  will be re-emitted as either Siii1527 or Siii*1533, and an absorbed Siii1304 photon can be re-emitted as Siii*1309. In all three cases the probability of emission in the resonant and fine structure transitions is approximately equal.

Since no absorption is seen in the fine structure transitions, we infer that atoms in the excited ground state will typically decay to the ground state before absorbing another photon. Since every absorbed photon is re-emitted, the net equivalent width of the resonant and fine structure transitions is . The precise equivalent widths depend on the initial absorption and the optical depth. If the gas is optically thin, re-emitted photons will immediately escape giving . In the limit of optically thick gas, resonant photons will be continuously scattered until they emerge as Siii* (after 2 scatterings on average). In this case, and . Since the majority of absorption occurs in optically thick gas (§4.3), we expect the equivalent width of Siii absorption lines to reflect the kinematics and covering fraction of neutral gas with minimal contamination from Siii re-emission.

#### 4.3.2 The spatial extent of low-ionization absorption

From the picture above (§4.3.1), we expect that the equivalent width of the fine structure emission lines should be equal and opposite to the resonant line equivalent width. This is not the case: and . ( is more difficult to quantify since Siii1304 is blended with Oi1302). There are two possible explanations. One is that scattered photons have larger path lengths and are subject to greater dust attenuation. Since the UV continuum slopes imply little differential extinction (mean E(B-V) ) and we expect only 2 scatterings for the average Siii* photon, this seems unlikely. More reasonably, the emitting region could be larger than that sampled by our slits. Although the slit aperture samples most of the continuum light and its line of sight absorption, the CGM which absorbs and isotropically re-emits scattered photons extends to much larger radii (Steidel et al., 2010, 2011). The strength of fine structure emission provides a direct measurement of the amount of absorption at small radii contained within the slit. The fine structure to resonant absorption line ratio suggests that a fraction (combining both 1265 and 1527 measures) of the total fine structure emission is contained within the extraction aperture of 1 arcsec. The half-light radius of Siii* emission, and hence Siii absorption, thus corresponds to kpc at . This scale is only slightly larger than the median half-light radius of galaxies in our sample, and would be reached in only 20 Myr at the typical outflow velocity (190 km s).

It is instructive to reconsider the composite spectrum from Shapley et al. (2003) where the fine structure lines are clearly seen (but were not interpreted fully along the discussion above at the time). The Siii* lines are noticeably stronger in our composite, with an average at c.f. Å  at . Furthermore, the Siii absorption lines are weaker in the composite, suggesting the absorption takes place at larger radii at . This difference cannot be explained through instrumental differences between the two surveys. Smoothing our composite to match the 3.25 Å  resolution of Shapley et al. (2003) reduces the equivalent widths by %. Likewise the wider slits used for the data (14, 11 kpc) c.f. the data (10; 7 kpc) is not the cause. If the CGM properties are similar, we expect the fine structure strength to constitute a larger fraction of the resonant absorption line strength in the data, contrary to observations. Defining , we find for the composite spectrum and for , indicating a smaller characteristic radius of fine structure emission at .

Inescapably, therefore, we conclude the circumgalactic gas around LBGs at differs in two important ways from that at : there is greater low-ionization absorption at small radii at and less total low-ionization absorption. In addressing the origin of this effect, Figure 10 shows that the difference in equivalent width is at least partially linked to the kinematic offset from the Ly emission. Thus it seems a higher covering fraction at small radii is required to produce the stronger fine structure emission. Higher resolution observations of individual bright, or perhaps gravitationally-lensed, galaxies at 4 will ultimately be required to separate the relative contributions of covering fraction and kinematics in explaining this result.

### 4.4. High-ionization lines

The high-ionization lines Siiv and Civ arise both in interstellar gas and in stellar P-Cygni winds. The velocity centroid of Siiv, km s, is consistent with the low-ionization interstellar absorption lines suggesting that most absorption is interstellar in origin. This is supported by the absence of a significant redshifted emission component expected for a P-Cygni profile. In contrast, Civ is broader with a larger absorption velocity offset km sand redshifted P-Cygni emission indicating a large contribution from stellar winds. We also detect the Nv1239,1243 P-Cygni feature, although the proximity to Ly makes this feature difficult to study in detail.

We can determine the optical depth of highly ionized outflowing gas from the ratio of Siiv absorption lines. The ratio for optically thin absorption, and 1.0 in the optically thick case. Shapley et al. (2003) find optically thin absorption in composite spectra of LBGs at , whereas we measure a ratio in the composite spectrum shown in Figure 3 indicating a significant contribution of optically thick absorption at . The total equivalent width of the Siiv doublet is weaker by a factor in Figure 3 compared to the composite of Shapley et al. (2003). The combination of higher optical depth and lower equivalent width suggests that the velocity range and/or covering fraction of the ionized gas traced by Siiv is lower at higher redshift.

#### 4.4.1 Metallicity

The P-Cygni profile of Civ is sensitive to metallicity, and the combination of Civ and Heii equivalent widths constrains both the age and metallicity. Here we compare the equivalent widths measured in the composite spectrum (Figure 3) with theoretical models in order to estimate the typical metallicity of galaxies in our sample. We note that the equivalent width of Civ contains significant interstellar absorption, so the value reported in Table 1 should be treated as an upper limit on the P-Cygni component. Assuming the interstellar absorption component of Civ is similar to that of Siiv ( Å), we take the P-Cygni absorption component to have equivalent width Å  with a conservative uncertainty. We compare this estimate and the measured equivalent width of Heii (reported in Table 1) with predictions from the stellar population synthesis code BPASS presented in Eldridge & Stanway (2009). We consider BPASS models which include binary evolution with continuous star formation rate, and determine the difference between observed and predicted equivalent width as a function of metallicity and relative carbon abundance. We restrict the stellar population age to be within the 1 scatter of the median value for galaxies in our sample, determined to be years from spectral energy density fits assuming constant star formation with a Kroupa (2002) initial mass function. Figure 7 shows the resulting relative error between BPASS models and measured equivalent width, which we define as

 Relativeerror=0.5⎡⎣(ΔWC{\sc iv}σC{\sc iv})2+(ΔWHe{\sc ii}σHe{\sc ii})2⎤⎦. (3)

Relative error values are thus consistent with the data.

Models with and somewhat depleted carbon abundance are in good agreement with the data. We note that observations of both local and high-redshift galaxies indicate typical relative carbon abundances (C/O) lower than the solar value (equivalent to in Figure 7) for metallicities (Kobulnicky & Skillman, 1998; Shapley et al., 2003; Erb et al., 2010; Eldridge & Stanway, 2011). Solar metallicity models () do not fit the data. We conclude that the typical metallicity of galaxies in our sample is or solar metallicity. Bright galaxies at in the FORS2 sample have gas-phase metallicity or about the solar value, in reasonable agreement (Maiolino et al., 2008). These values are consistent with the allowed metallicity range inferred for the Shapley et al. (2003) composite using the same BPASS models (Eldridge & Stanway, 2011). Given that the equivalent widths of Civ and Heii in Shapley et al. (2003) agree well with those in Table 1, we expect the typical metallicity of galaxies to be similar in both samples.

## 5. Spectroscopic trends

We now turn to an analysis of how the various spectroscopic features discussed in §4 are related to observable properties of LBGs as a prelude to considering how they might evolve with redshift. We will begin with the dependence of low-ionization absorption line strength with Ly equivalent width. This is motivated by the common physical dependence of these features – both are governed by the kinematics and covering fraction of neutral circumgalactic gas – and also because this has been examined in detail for LBGs at (Shapley et al., 2003).

We construct composite spectra of three sub-samples of galaxies divided according to their Ly equivalent width. Defining as the average equivalent width of 1260, 1303, 1334, and 1527 Å features, Figure 8 shows weaker for galaxies with stronger Ly . This was also noted for LBGs by Vanzella et al. (2009). Both the =3 and =4 samples show this trend, although is weaker at for galaxies with strong Ly emission. This difference can be attributed to the luminosity-dependent trend of stronger and weaker in fainter galaxies (Paper I; Shapley et al. 2003; Vanzella et al. 2009). If we consider a subset of galaxies in our sample with similar absolute magnitudes to those observed at (), we recover the same normalization (Figure 8). Our sample is 90% complete at the corresponding apparent magnitudes, so we expect a negligible bias.

Having established the correlation of with , we now examine the dependence with other demographic properties. We divide the full sample into two bins of equal size according to each property of interest. The results are shown in Figure 9. We briefly review the trend of each property with and and discuss the physical origin. Many of these trends were seen in Shapley et al. (2003) and Vanzella et al. (2009) and are clearly inter-related due to correlations between the demographic properties.

Less luminous galaxies have stronger and weaker illustrating that higher star formation rates drive larger amounts of neutral gas into the CGM with higher velocity and/or covering fraction. Defining the ultraviolet spectral slope as for B-drops (Bouwens et al., 2009), we also find that bluer galaxies with lower have stronger and weaker . This trend was also noted in Paper I, which showed that LBGs with strong Ly emission have systematically bluer . Since neutral gas also presumably contains dust, the same gas which gives rise to also reddens the continuum. Stellar masses are measured for 60% of our sample for which there is unconfused Spitzer/IRAC photometry. As expected from the trends with luminosity, lower mass galaxies have stronger and weaker . Finally, measuring half-light radii from the GOODS ACS data with Sextractor, we find smaller galaxies have stronger and weaker .

Could these trends could be due to selection effects or sample bias? The trend of stronger Ly for less luminous LBGs is of particular concern, since fainter galaxies will require a larger for detection. To address this, we consider only the 32 galaxies with apparent magnitudes for which the spectroscopic sample is 90% complete. This sub-sample is divided into two equal bins and the results confirm that trends seen in the larger sample also hold in brighter galaxies unaffected by sample bias.

The trends shown in Figure 9 are generally consistent with the overall trend of with (Figure 8). The variations mostly arise from the distribution of neutral gas (traced by ), with relatively little effect from other demographic properties examined (, , , ). However, demographic properties do have some effect. The strongest deviation seen in Figure 9 is that with , which shows stronger than would be expected at a given for galaxies with red UV slopes. Noted also by Shapley et al. (2003), this suggests that outflowing neutral gas contains dust which reddens the continuum. Also, more luminous galaxies (i.e. those with higher star formation rates) have stronger at a given (Figure 8). This is likely due to increased absorption at large velocities in galaxies with higher SFR, as observed at (Weiner et al., 2009). For all other demographics, the composite spectra have values of within of that expected purely based on .

In summary, the trends seen in Figure 9 arise almost entirely because of variations in the neutral gas covering fraction and/or kinematics, which are themselves correlated with the demographic properties.

### 5.1. Kinematics

We have now established that outflowing neutral gas is the dominant factor in determining both and . However, variations from the trend of with are apparent, particularly with and as discussed above. The most obvious mechanism for this behaviour is a systematic difference in the covering fraction and kinematics of neutral gas. Assuming (where is the outflow velocity), a higher and lower velocity range can conspire to give a constant . However, changing and will also affect the transmission of Ly photons resulting in a different . It is therefore of interest to consider how to distinguish between the covering fraction and kinematics of neutral gas.

The kinematics of neutral gas can be roughly parameterized by the velocity dispersion and centroid of low-ionization interstellar absorption lines. We have measure the velocity of low-ionization absorption lines with respect to Ly as a proxy for outflow velocity. is strongly correlated with and in LBGs at (Shapley et al., 2003) in the sense that larger velocities are associated with stronger interstellar absorption and weaker Ly emission. Figure 10 shows measured from the same composite spectra used to determine demographic trends (e.g. Figure 9), as well as the results from Shapley et al. (2003). All composites are consistent (within ) with the relation measured at as well as the mean measured for B-dropout galaxies by Vanzella et al. (2009). We measure a velocity dispersion for each composite, and find that each is within of the effective spectral resolution km s(measured for the stellar [Ciii] feature in Figure 3). We are therefore unable to detect trends in outflow kinematics with demographic properties in composite spectra. Higher signal-to-noise data, higher spectral resolution, or detailed studies of individual galaxies are required to address trends in kinematics and covering fraction of neutral gas at .

## 6. The evolving CGM

We now turn to a discussion of the redshift evolution of the neutral CGM surrounding typical LBGs. The most useful probes are Ly and low-ionization absorption lines which trace the kinematics and covering fraction of neutral gas, and fine structure emission lines which provide a constraint on the spatial extent of the absorbing gas (§4.3.2).

The common dependence of Ly and low-ionization absorption lines on the neutral CGM results in a strong correlation between and . Various physical properties of LBGs are correlated with both and , but in such a way that the relation between and remains nearly constant (Figure 9). Furthermore, the relation at fixed does not change significantly with redshift between and (Figure 8).

To further examine evolutionary trends with redshift, we now divide our spectroscopic sample into two bins of redshift at fixed , now including galaxies at all redshifts (no longer restricted to as in previous sections). We consider galaxies with absolute magnitude , chosen to be representative of the sample in Shapley et al. (2003). There are 64 galaxies in our sample within this range (see Figure 2). We construct composite spectra of galaxies with redshift above and below the median and measure the equivalent width of and low-ionization lines (both resonant absorption and fine structure emission) in each composite. The results are given in Table 2 along with the demographic properties of galaxies in each sub-sample. The quantities , , and are measured from photometry while is determined from a direct fit to the ultraviolet continuum in the composite spectrum. We define and fit the rest frame Å to determine , with uncertainty quantified using the bootstrap method described in §3. Aside from redshift, the demographics of each sub-sample are quite similar. The higher redshift galaxies have slightly higher average , smaller , and smaller (bluer) . is consistent for both to within the sample variance, and is also consistent with the value Å  measured for the composite spectrum of LBGs in Shapley et al (2003). The most striking difference is in the strength of the low-ionization absorption lines, which are significantly weaker at higher redshifts (Figure 11). The variation in is not explained by systematic differences in or demographic properties, hence we seek an alternate explanation.

We first examine whether the evolution in could arise as a result of different equivalent width distributions for Ly . Although the mean across our two redshift subsamples is similar, the lower redshift subsample has a broader distribution and contains more galaxies with Ly in absorption (). This is reflected in the larger sample variance in at lower redshift (Table 2). We evaluate the effect of this potential bias on by constructing a composite spectrum from a subset of the galaxies with intermediate Å, resulting in a consistent mean with sample variance reduced by a factor of 2.5. This composite has Å, 0.2 Å higher than when the full range of is used, but still considerably lower than the value Å measured for the higher redshift galaxies. We therefore conclude that differences in the distribution are insufficient to explain the observed variation in absorption line strength with redshift.

There are several possible physical explanations for the evolution of with redshift shown in Figure 11. One possibility is that the kinematics and/or covering fraction of neutral gas are systematically different. For example, an outflowing wind with fixed input energy and momentum will reach higher velocity at lower redshifts due to lower density of the IGM. However, this effect should be stronger between than from , whereas the decrement in is much stronger from (Figure 11). Furthermore, we measure a higher offset between the velocity centroid of Ly emission and low-ionization absorption for galaxies () than for (). Both measurements are consistent with the trend shown in Figure 10. Based on the trend of with seen at (Shapley et al., 2003), we would then expect stronger at higher redshift. Kinematics are thus unable to explain the difference in , at least with the information currently available. Data with higher spectral resolution and signal-to-noise are needed to fully address differences in the covering fraction and kinematics of neutral gas.

Another obvious potential cause of weaker , or equivalently decreased at a fixed , is the presence of neutral hydrogen with very low column density of heavy elements. In this case we would expect optically thin absorption line profiles. This could be directly tested as shown in §4.3.1, but the signal-to-noise of Siii in the high redshift composite is too low to constrain the optical depth. We do see some evidence that the column density of neutral gas is lower at higher redshifts based on Ly line profiles in the composite spectra, shown in Figure 12. The Ly absorption trough at Å  is significantly weaker in the higher redshift composite. This absorption arises at least in part from damping wings of high column density gas with cm associated with the low-ionization metal absorption lines (e.g. Pettini et al. 2000, 2002). The higher redshift galaxies are therefore characterized by lower typical or/and a lower covering factor of high column density gas. Higher signal-to-noise data at is required to determine whether this affects the optical depth of low-ionization metal transitions.

Weaker low-ionization lines could also be caused by a systematically higher ionization state at higher redshifts resulting in lower column densities of low-ionization gas. This would also explain the weaker damped Ly absorption trough seen at higher redshifts (Figure 12). This scenario can be tested by measuring the equivalent widths of higher-ionization silicon transitions, in particular Siiii1206 and Siiv1393,1402. We find that the Siiv lines are weaker at higher redshift by a factor of in , consistent with the decrement in . We do not detect the Siiii1206 line in the high redshift composite, with a 1 upper limit of times that of the lower redshift composite. Siiii1206 is marked in Figure 12 and is clearly stronger at lower redshift. We therefore find no evidence of a significant change in the ionization state of outflowing gas at higher redshifts.

A final possibility for lower column density of heavy elements is that Ly is scattered by “cold streams” of nearly metal-free neutral gas, predicted by simulations to accrete onto galaxies at small radii and with increasing rates at higher redshift (e.g. Dekel et al. 2009; Faucher-Giguère et al. 2011). If they exist, such streams could be identified as Hi absorption systems with no corresponding absorption from heavy elements. Steidel et al. (2010) find little evidence for the presence of such streams around LBGs at . A similar study at higher redshift would be very interesting but also extremely challenging with current observational facilities.

Could the difference in with redshift be caused in part by a different physical extent of circumgalactic gas? In §4.3.2 we showed that fine structure emission lines can be used to measure the amount of absorption taking place within the size of the spectroscopic aperture. This quantity does appear to change with redshift, in the sense that absorption takes place at smaller characteristic radius in LBGs at higher redshift. The composite spectra of LBGs at has weaker absorption lines and stronger fine structure emission than the composite of Shapley et al. (2003) (Figure 4), in qualitative agreement with trends seen at (Shapley et al. 2003; see their Figure 9). We have only a weak constraint on the redshift evolution at due to limited signal-to-noise in the high redshift composite. The ratio of total Siii fine structure emission in the high- and low-redshift composites is (taken as

 W1265+W1309+W1533W1260+W1303+W1527

with values given in Table 2), while the corresponding ratio of absorption line strength is . The inferred fraction of low-ionization absorption within the slit aperture is thus a factor larger in the higher redshift sample, consistent with no evolution. Higher signal-to-noise is required to accurately constrain the spatial extent of low-ionization absorption at .

To summarize, we find weaker low-ionization absorption lines in LBGs at higher redshift with fixed . The difference in is not consistent with the trends observed for Ly and various demographic properties examined in §5. We have discussed several possible causes for this discrepancy including variation in the kinematics, covering fraction, optical depth, ionization state, and spatial extent of absorbing gas. The data show that the ionization state has no significant effect, but we are unable to conclusively address other possible causes. Spectra of galaxies at with higher signal-to-noise are required to more accurately constrain the optical depth and spatial extent. Additionally, spectra of individual galaxies taken with higher spectral resolution will be required to independently determine the covering fraction and kinematic structure of absorbing gas.

### 6.1. Galaxy evolution

We argue in §5 that Ly equivalent width is determined primarily by the neutral CGM, which is correlated with various demographic galaxy properties (Figure 9). It is well established that these demographic properties vary with redshift and we expect to vary accordingly. Large photometric surveys have shown that Lyman break galaxies at increasingly higher redshifts have lower luminosities (Bouwens et al., 2007, 2011a), bluer UV spectral slopes (Bouwens et al., 2009, 2011b), smaller stellar masses (Stark et al., 2009; González et al., 2011), and smaller sizes (Bouwens et al., 2004; Ferguson et al., 2004). This is in accordance with inside-out galaxy growth: galaxies increase in size and stellar mass as they evolve with time, while increasing metallicity and dust content reddens the ultraviolet continuum. Simultaneously, star formation drives large-scale outflows of gas which reach larger distances and are accelerated to larger velocities at later times (e.g. Murray et al. 2010). Galaxies which are more evolved (i.e. larger, more massive, redder) should therefore have a CGM characterized by larger spatial extent, larger velocity range, and higher covering fraction of neutral gas. Observationally this results in weaker , stronger , and relatively weaker fine structure emission. These are precisely the trends observed at both (Shapley et al., 2003) and (this work).

### 6.2. Lyα in the epoch of reionization

We reiterate that Lyman break galaxies at increasingly higher redshifts have lower luminosities, bluer UV spectral slopes, smaller stellar masses, and smaller sizes. Notably, all trends in the demographics of galaxies at higher redshift are correlated with stronger and weaker (Figure 9). We therefore expect typical galaxies at higher redshifts to have, on average, stronger Ly emission. Earlier results from this survey confirm that strong Ly emission is more frequent in galaxies at higher redshift (Paper I; Paper II). We find no evidence of this trend reversing. In fact, galaxies with extremely small size, low mass, and blue tend to be the strongest Ly emitters (e.g. Erb et al. 2010). We do, however, expect the average Ly emission strength to decrease significantly at increasing redshifts in the epoch of reionization due to neutral hydrogen in the IGM (e.g. Haiman & Spaans 1999).

We have devoted considerable discussion to the properties of Ly in part because Ly is of great interest as a tracer of cosmic reionization. Several authors have now presented evidence that reionization was incomplete at based on a rapidly decreasing fraction of galaxies with strong Ly emission at (Schenker et al., 2011; Ono et al., 2011; Pentericci et al., 2011). Although trends at suggest that the galaxies observed by these authors should have a higher fraction of strong Ly emission, we have presented evidence that the relation between and is systematically different at (Figure 11). This implies a systematic difference in the spatial, kinematic, or optical depth structure of neutral circumgalactic gas compared to galaxies at lower redshift. The physical origin of this evolution and its effect on at will need to be understood in order to fully interpret the results of Ly surveys at higher redshifts in the context of reionization.

## 7. Summary

The rest-frame ultraviolet spectra of star-forming galaxies contains a wealth of information about the properties of the circumgalactic medium. In this paper we have presented an analysis of several features which trace the CGM with a focus on the properties of neutral gas. We find that the trends observed at lower redshift (; Shapley et al. 2003) also hold at with approximately the same normalization (Figure 8). However, we find evidence for rapid evolution at with lower at fixed and , suggesting a systematic difference in the spatial distribution, kinematics, ionization state, or optical depth of circumgalactic gas at higher redshifts. We determine that the ionization state is not responsible for the observed evolution but are unable to distinguish between kinematics, covering fraction, optical depth, or the spatial extent of neutral gas as the likely cause. We are collecting additional spectra of LBGs at with our ongoing survey, including high spectral resolution observations of bright lensed galaxies from which we can disentangle the kinematic profile and covering fraction of neutral gas. These data will allow us to address the precise magnitude and physical origin of evolution in circumgalactic gas properties.

As a final note, we emphasize that the neutral CGM is of great interest in the context of reionization of the universe. Neutral gas in the circumgalactic medium absorbs ionizing radiation, thereby inhibiting the ability of galaxies to reionize the universe. The escape fraction of ionizing photons is one of the most important and uncertain factors in determining the contribution of star-forming galaxies to reionization (Robertson et al., 2010). We have shown that typical galaxies at higher redshift have weaker low-ionization absorption lines based on their demographic trends, and presented new evidence that absorption lines are systematically weaker at even for fixed demographic properties. This is likely caused by a lower covering fraction and/or velocity range of neutral gas, and we will address the physical origin of this evolution with future data from our ongoing survey. Determining the redshift evolution of neutral gas covering fraction in LBGs will be of great interest for interpreting surveys of Ly emission in the context of reionization and addressing the role of star-forming galaxies in reionizing the universe.

## Acknowledgments

D.P.S. acknowledges support from NASA through Hubble Fellowship grant #HST-HF-51299.01 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc, for NASA under contract NAS5-265555. R.S.E. acknowledges the hospitality of Piero Madau and colleagues at the University of California, Santa Cruz where this work was completed. We thank Masami Ouchi, Chuck Steidel, Max Pettini, Anna Quider, Crystal Martin, Alice Shapley, and Gwen Rudie for helpful discussions. The analysis pipeline used to reduce the DEIMOS data was developed at UC Berkeley with support from NSF grant AST-0071048. Most of the data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain.

## References

• Balestra et al. (2010) Balestra, I., Mainieri, V., Popesso, P., et al. 2010, A&A, 512, A12
• Bouwens et al. (2004) Bouwens, R. J., Illingworth, G. D., Blakeslee, J. P., Broadhurst, T. J., & Franx, M. 2004, ApJ, 611, L1
• Bouwens et al. (2007) Bouwens, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, ApJ, 670, 928
• Bouwens et al. (2009) Bouwens, R. J., Illingworth, G. D., Franx, M., et al. 2009, ApJ, 705, 936
• Bouwens et al. (2011a) Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2011, ApJ, 737, 90
• Bouwens et al. (2011b) Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2011, arXiv:1109.0994
• Brammer et al. (2011) Brammer, G. B., Whitaker, K. E., van Dokkum, P. G., et al. 2011, ApJ, 739, 24
• Davis et al. (2003) Davis, M., Faber, S. M., Newman, J., et al. 2003, Proc. SPIE, 4834, 161
• Dekel et al. (2009) Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 451
• Eldridge & Stanway (2009) Eldridge, J. J., & Stanway, E. R. 2009, MNRAS, 400, 1019
• Eldridge & Stanway (2011) Eldridge, J. J., & Stanway, E. R. 2011, MNRAS, 1595
• Ellis (2008) Ellis, R. S. 2008, Saas-Fee Advanced Course 36: First Light in the Universe, 259
• Erb et al. (2010) Erb, D. K., Pettini, M., Shapley, A. E., et al. 2010, ApJ, 719, 1168
• Faber et al. (2003) Faber, S. M., Phillips, A. C., Kibrick, R. I., et al. 2003, Proc. SPIE, 4841, 1657
• Faucher-Giguère et al. (2011) Faucher-Giguère, C.-A., Kereš, D., & Ma, C.-P. 2011, MNRAS, 417, 2982
• Ferguson et al. (2004) Ferguson, H. C., Dickinson, M., Giavalisco, M., et al. 2004, ApJ, 600, L107
• Forero-Romero et al. (2011) Forero-Romero, J. E., Yepes, G., Gottloeber, S., & Prada, F. 2011, arXiv:1109.0228
• Giavalisco et al. (2004) Giavalisco, M., Ferguson, H. C., Koekemoer, A. M., et al. 2004, ApJ, 600, L93
• González et al. (2011) González, V., Labbé, I., Bouwens, R. J., et al. 2011, ApJ, 735, L34
• Grogin et al. (2011) Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, arXiv:1105.3753
• Haiman & Spaans (1999) Haiman, Z., & Spaans, M. 1999, ApJ, 518, 138
• Heckman (2002) Heckman, T. M. 2002, Extragalactic Gas at Low Redshift, 254, 292
• Hopkins & Beacom (2006) Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142
• Koekemoer et al. (2011) Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, arXiv:1105.3754
• Kobulnicky & Skillman (1998) Kobulnicky, H. A., & Skillman, E. D. 1998, ApJ, 497, 601
• Kroupa (2002) Kroupa, P. 2002, Science, 295, 82
• Maiolino et al. (2008) Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463
• Masters & Capak (2011) Masters, D., & Capak, P. 2011, PASP, 123, 638
• Murray et al. (2010) Murray, N., Quataert, E., & Thompson, T. A. 2010, ApJ, 709, 191
• Oke (1974) Oke, J. B. 1974, ApJS, 27, 21
• Ono et al. (2011) Ono, Y., Ouchi, M., Mobasher, B., et al. 2011, arXiv:1107.3159
• Ouchi et al. (2008) Ouchi, M., Shimasaku, K., Akiyama, M., et al. 2008, ApJS, 176, 301
• Pentericci et al. (2011) Pentericci, L., Fontana, A., Vanzella, E., et al. 2011, arXiv:1107.1376
• Pettini et al. (2000) Pettini, M., Steidel, C. C., Adelberger, K. L., Dickinson, M., & Giavalisco, M. 2000, ApJ, 528, 96
• Pettini et al. (2002) Pettini, M., Rix, S. A., Steidel, C. C., et al. 2002, ApJ, 569, 742
• Quider et al. (2009) Quider, A. M., Pettini, M., Shapley, A. E., & Steidel, C. C. 2009, MNRAS, 398, 1263
• Quider et al. (2010) Quider, A. M., Shapley, A. E., Pettini, M., Steidel, C. C., & Stark, D. P. 2010, MNRAS, 402, 1467
• Reddy & Steidel (2009) Reddy, N. A., & Steidel, C. C. 2009, ApJ, 692, 778
• Robertson et al. (2010) Robertson, B. E., Ellis, R. S., Dunlop, J. S., McLure, R. J., & Stark, D. P. 2010, Nature, 468, 49
• Schenker et al. (2011) Schenker, M. A., Stark, D. P, Ellis, R. S., et al. 2011, arXiv:1107.1261
• Shapley et al. (2001) Shapley, A. E., Steidel, C. C., Adelberger, K. L., et al. 2001, ApJ, 562, 95
• Shapley et al. (2003) Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 588, 65
• Shapley (2011) Shapley, A. E. 2011, ARA&A, 49, 525
• Stark et al. (2009) Stark, D. P., Ellis, R. S., Bunker, A., et al. 2009, ApJ, 697, 1493
• Stark et al. (2010) Stark, D. P., Ellis, R. S., Chiu, K., Ouchi, M., & Bunker, A. 2010, MNRAS, 408, 1628
• Stark et al. (2011) Stark, D. P., Ellis, R. S., & Ouchi, M. 2011, ApJ, 728, L2
• Steidel et al. (2003) Steidel, C. C., Adelberger, K. L., Shapley, A. E., et al. 2003, ApJ, 592, 728
• Steidel et al. (2010) Steidel, C. C., Erb, D. K., Shapley, A. E., et al. 2010, ApJ, 717, 289
• Steidel et al. (2011) Steidel, C. C., Bogosavljević, M., Shapley, A. E., et al. 2011, ApJ, 736, 160
• Vanzella et al. (2005) Vanzella, E., Cristiani, S., Dickinson, M., et al. 2005, A&A, 434, 53
• Vanzella et al. (2006) Vanzella, E., Cristiani, S., Dickinson, M., et al. 2006, A&A, 454, 423
• Vanzella et al. (2008) Vanzella, E., Cristiani, S., Dickinson, M., et al. 2008, A&A, 478, 83
• Vanzella et al. (2009) Vanzella, E., Giavalisco, M., Dickinson, M., et al. 2009, ApJ, 695, 1163
• Wang et al. (2010) Wang, W.-H., Cowie, L. L., Barger, A. J., Keenan, R. C., & Ting, H.-C. 2010, ApJS, 187, 251
• Wardlow et al. (2011) Wardlow, J. L., Smail, I., Coppin, K. E. K., et al. 2011, MNRAS, 415, 1479
• Weiner et al. (2009) Weiner, B. J., Coil, A. L., Prochaska, J. X., et al. 2009, ApJ, 692, 187

## References

You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters