Detection of hot, metal-enriched outflowing gas

Detection of hot, metal-enriched outflowing gas around z2.3 star-forming galaxies in the Keck Baryonic Structure Survey


We use quasar absorption lines to study the physical conditions in the circumgalactic medium of redshift star-forming galaxies taken from the Keck Baryonic Structure Survey (KBSS). In Turner et al. (2014) we used the pixel optical depth technique to show that absorption by \hone and the metal ions \osix, \nfive, \cfour, \cthree and \sifour is strongly enhanced within  \kmps and projected distances proper kpc from sightlines to the background quasars. Here we demonstrate that the \osix absorption is also strongly enhanced at fixed \hone, \cfour, and \sifour optical depths, and that this enhancement extends out to  \kmps. At fixed \hone the increase in the median \osix optical depth near galaxies is 0.3–0.7 dex and is detected at 2–3- confidence for all seven \hone bins that have . We use ionization models to show that the observed strength of \osix as a function of \hone is consistent with enriched, photoionized gas for pixels with . However, for pixels with this would lead to implausibly high metallicities at low densities if the gas were photoionized by the background radiation. This indicates that the galaxies are surrounded by gas that is sufficiently hot to be collisionally ionized (K) and that a substantial fraction of the hot gas has a metallicity of solar. Given the high metallicity and large velocity extent (out to ) of this gas, we conclude that we have detected hot, metal enriched outflows arising from star-forming galaxies.

galaxies: formation – intergalactic medium – quasars: absorption lines

1 Introduction

Galaxy formation models predict that massive galaxies are surrounded by haloes of hot, chemically enriched gas, which may be penetrated by accreting streams of cold and largely metal-poor gas (e.g. Kereš et al., 2005; Dekel et al., 2009; van de Voort et al., 2011; van de Voort & Schaye, 2012). The hot component is heated through shocks associated with galactic winds and gas accretion, while the cold component is photo-heated to temperatures of K. The circumgalactic medium thus results from the inflow of gas into the potential well set by the dark matter halo, and the outflows driven by feedback from star formation and/or active galactic nuclei (AGN). Hence, the gas around galaxies holds valuable clues to the fuelling and feedback processes that currently limit our understanding of galaxy evolution.

This circumgalactic gas is very diffuse, making it difficult to detect in emission. Instead, the gas can be studied in absorption using high-quality spectra of background quasars. However, the low number density of bright quasars and the difficulty of obtaining accurate galaxy redshifts make it challenging to apply this technique to systematic surveys of the circumgalactic medium.

The advent of the Cosmic Origin Spectrograph on the Hubble Space Telescope and the MOSFIRE spectrograph on the Keck telescope have recently improved the situation by significantly increasing the quality of low-redshift quasar spectra and the accuracy with which redshifts can be measured for high-redshift galaxies, respectively. In particular, Tumlinson et al. (2011) found \osix to be ubiquitous within 150 proper kiloparsec (pkpc) of star-forming galaxies with median halo masses   (Werk et al., 2014). At , Turner et al. (2014) found strong enhancements of \osix, \nfive, \cfour, \cthree, \sifour, and \hone within 180 pkpc and 240 kms of star-forming galaxies that are also thought to be hosted by haloes with masses . These studies confirmed and extended earlier work on the association between galaxies and absorbers.

The ion \osix is of particular interest, because its absorption is relatively strong and because it can trace gas at  K, similar to the temperatures to which the gas is expected to be heated in shocks associated with winds and accretion events. Furthermore, it is the temperature range for which gas cooling from much higher temperatures is most likely to be detectable. Indeed, a number of simulations predict the presence of large amounts of \osix around galaxies (e.g. Stinson et al., 2012; Ford et al., 2013; Shen et al., 2013). In general, these studies find that \osix resides in a collisionally ionized gas phase for impact parameters  pkpc, outside of which \osix tends to be found in a cooler, photoionized gas. Tepper-García et al. (2011) estimate that two-thirds of \osix absorbers at in their simulations trace gas with  K (see also Oppenheimer et al. 2012).

At low redshifts, many observational programmes have found that \osix absorbers tend to be located within 300 pkpc of galaxy positions (e.g. Stocke et al., 2006; Chen & Mulchaey, 2009; Prochaska et al., 2011; Tumlinson et al., 2011). Although the presence of \osix around galaxies is well established, the properties of the gas that it traces are still under debate, as \osix systems are observed to arise in both photoionized and collisionally ionized phases (e.g, Danforth & Shull, 2008; Savage et al., 2014). While some studies find that \osix in galaxy haloes is consistent with being photoionized (e.g. Prochaska et al., 2011), many circumgalactic \osix absorbers are thought to be in a phase distinct from that of both \hone and lower ions at  K, due to their often complex and differing kinematic structures (e.g. Tripp et al., 2008) and/or inferences from ionization modelling (e.g. Werk et al., 2013). Using \hone \lya absorber stacking2 Pieri et al. (2014) found that for high ionization lines including \osix, their observations were most consistent with  K, with a possible contribution from warmer  K gas (inferred from N/O measurements) .

Because of its near-coincidence with the \hone \lyb forest, studying \osix absorption using line-fitting techniques becomes progressively more difficult with increasing redshift, although individual line fitting is still possible and has been done extensively at (e.g. Bergeron et al., 2002; Carswell et al., 2002; Simcoe et al., 2002, 2004, 2006; Lopez et al., 2007; Schaye et al., 2007; Fox et al., 2008). Rather than looking by eye for individual \osix absorption lines, \osix has also been studied using the automatic pixel optical depth technique to correlate the absorption from \osix with that of other ions (Cowie & Songaila, 1998; Schaye et al., 2000a; Aguirre et al., 2008; Pieri et al., 2010), optionally correcting for much of the contamination by \hone Lyman series lines (Aguirre et al., 2002, 2008; Turner et al., 2014).

Aracil et al. (2004) divided a sample of pixel optical depths into those near and far from strong \hone \lya absorption. They found that for gas at the same \hone optical depth, the amount of \osix was enhanced near strong \hone \lya absorption for . A similar study was then undertaken by Pieri et al. (2006) for the quasar Q (). Motivated by the very strong correlation between galaxies and strong \cfour absorbers found by Adelberger et al. (2003, 2005a), they used \cfour optical depth as a galaxy proxy, and found enhancements in both \cfour and \osix at fixed \hone close to galaxy positions.

In this work we extend the technique used in Pieri et al. (2006) and apply it to a much larger sample of 15 quasars and 854 spectroscopically confirmed galaxies (with impact parameters as small as 40 pkpc) taken from the Keck Baryonic Structure Survey (KBSS, Steidel et al. 2014). Rakic et al. (2012) and Rudie et al. (2012b) have already used an earlier version of the survey to study the distribution of neutral hydrogen around the galaxies, while Turner et al. (2014) measured the distribution of metal ions using the same data as is analysed here. Using a galactocentric approach and the pixel optical depth technique, Turner et al. (2014) found metal-line absorption to be strongly enhanced with respect to random regions for impact parameters  pkpc and line of sight (LOS) distances within  \kmps of the galaxy positions (or 1 pMpc in the case of pure Hubble flow). Furthermore, thanks to observations using MOSFIRE, the elongation of enhancement along the LOS was determined to be largely caused by gas peculiar velocities (rather than redshift errors).

While Turner et al. (2014) studied optical depth as a function of galaxy distance, here we will measure the enhancement of metal-line absorption at fixed \hone optical depth. This enables us to tell whether the enhancement in metal absorption near galaxies found by Turner et al. (2014) merely reflects the higher gas densities implied by the observed increase in \hone absorption, or whether it indicates that the circumgalactic gas has a higher metallicity or a different temperature compared to random regions with the same \hone optical depth.

This paper is organized as follows. In Section 2 we briefly review the properties of the galaxy and quasar samples, the galaxy redshifts, and the pixel optical depth technique. In Section 3 we present the principal observational results of this paper, which is that we find a strong and significant enhancement of \osix at fixed \hone, \cfour and \sifour for impact parameters  pkpc and velocities  \kmps. In section 4 we consider whether the observational result can be explained (1) if the gas near galaxies is photoionized and metal rich, (2) if the gas is photoionized by radiation from stars in the nearby galaxies, or (3) if the enriched gas is collisionally ionized. We find that small galactocentric distance pixels with are in agreement with scenarios (1) and (3), while for those that have , only scenario (3), i.e. the presence of hot, collisionally ionized gas, provides a consistent explanation. Finally, in Section 5 we summarize and discuss our main results.

Throughout the paper, we use proper rather than comoving units (denoted as pkpc and pMpc), and employ cosmological parameters as measured from the Planck mission (Planck Collaboration et al., 2013), i.e.  \kmps Mpc, , , and .

2 Observations and method

This work makes use of procedures that were detailed in Turner et al. (2014), which we will briefly outline here.

Ion \taumedk S/N
Table 1: The log of the median optical depth, and the median continuum S/N of all pixels (with normalized flux ) in the redshift range considered for the particular ion and recovery method.

2.1 Galaxy & QSO samples

The KBSS is centred around the fields of 15 hyper-luminous QSOs, all of which have been observed extensively with Keck/HIRES and therefore have very high quality spectra. Details about the QSO reduction and analysis, including the fitting out of DLA wings, are described in Rudie et al. (2012b). The HIRES spectra typically have (which corresponds to a  \kmps), and S/N ranging from to  pixel (we give the median S/N for the spectral regions covered by each ion in Table 1).

The survey focuses on obtaining spectroscopic redshifts for the galaxies in the above QSO fields. The full sample currently consists of galaxies at , that were chosen to have redshifts in the range probed by the QSO spectra using UV-colour selection techniques from Steidel et al. (2003, 2004) and Adelberger et al. (2004). Spectroscopic follow-up using the instruments LRIS, NIRSPEC and/or MOSFIRE was then performed on galaxies with apparent magnitudes (see Rudie et al. 2012b for more information about the galaxy follow-up strategy). The above selection typically results in galaxies with halo masses  \msol(Adelberger et al., 2005b; Conroy et al., 2008; Trainor & Steidel, 2012; Rakic et al., 2013), which corresponds to virial radii and circular velocities of  pkpc and  \kmps, respectively. Furthermore, these galaxies tend to have dynamical masses  \msol(Erb et al., 2006b), median star formation rates  \msol yr (Erb et al., 2006a; Steidel et al., 2014), gas-phase metallicities (Steidel et al., 2014), and stellar ages  Gyr (Erb et al., 2006b).

In this work, we focus on the subset of 21 KBSS galaxies that satisfy the following two constraints. First, we limit the sample to galaxies with impact parameters  pkpc. Although we also considered galaxies with impact parameters up to  pMpc as was done in Turner et al. (2014), we did not find any differences in the results from using impact parameter bins  pkpc compared to those from random regions. The choice of 180 pkpc can be further motivated by the fact that Turner et al. (2014) found a strong metal optical depth enhancement above the median value for random regions up to the same impact parameter values.

The second constraint concerns the velocity direction, where we consider only galaxies that lie within the \lya forest of the background QSO, defined as:


The left-hand side corresponds to the beginning of the \lyb forest, while the right-hand side is set by the QSO emission redshift less an offset factor to avoid proximity effects. This second criterion is required so that the pixel optical depth analysis (as described in § 2.3) can be applied to the same redshift range for different ions.

2.2 Galaxy redshifts

NIRSPEC and MOSFIRE are both near-IR spectrographs, with galaxy redshift measurement uncertainties estimated to be and  \kmps, respectively. With these instruments, we are able to probe the rest-frame optical wavelengths for the galaxies in this sample, and to observe the nebular emission lines H, H, and [\othree] 4959,5007. Since these lines arise in \htwo regions of galaxies and, unlike Ly, are not subject to significant resonant scattering, they are thought to be robust tracers of the systemic redshift, and for our analysis we therefore take the galaxy redshift to be equal to that of the nebular emission lines, .

Although we are continuously working to increase the number of galaxies that have MOSFIRE observations, currently 4 out of the 21 galaxies in our sub-sample have been observed only with LRIS, which probes the rest-frame UV. From these data, we can measure the galaxy redshifts from either interstellar absorption lines () or \lya emission lines (). However, these lines tend to be systematically offset from the systemic galaxy redshifts (Shapley et al., 2003; Adelberger et al., 2003; Steidel et al., 2010; Rakic et al., 2011). To correct for this, we consider the galaxies that have both a nebular and rest-frame UV redshift measurement, and use the average difference between or and to determine a bulk correction value. The specifics of the correction implementation, as well as the latest offset values, can be found in § 2.2 of Turner et al. (2014). Our final sample of 21 galaxies contains 17 and 4 galaxies with redshifts measured from nebular and rest-frame UV features, respectively. Their minimum, median, and maximum impact parameters are 35, 118, and 177 pkpc, respectively.

2.3 Pixel optical depths

In this work, we use the pixel optical depth method (Cowie & Songaila, 1998; Ellison et al., 2000; Schaye et al., 2000a; Aguirre et al., 2002; Schaye et al., 2003; Turner et al., 2014) to study how absorption varies with galaxy proximity. Because of the statistical nature of our approach, we lose information about individual systems. However, by studying correlations between the pixel optical depths of different transitions, we are able to probe gas to lower densities, even in the presence of strong contamination, in a fast and fully automated manner.

A complete description of the method used is given in § 3 and Appendix A of Turner et al. (2014); here we give a short summary. We initially define the optical depth for ion and multiplet component as


where is the normalized flux at the pixel with a wavelength given by where is the rest-wavelength of the th transition of ion .

Next, we correct the absorption by each ion for saturation and for possible sources of contamination. Beginning with \hone, if \lya is saturated, we use the optical depths of higher-order transitions (\lyb, \lyg etc.) at the same redshift. Of all the higher-order pixels at a single redshift, we take the minimum of the optical depths that are not saturated (if there are any), scaled to the \lya transition, and use it to replace the saturated \lya value.

The recovered \hone \lya is then used to subtract 5 orders of the Lyman series of \hone (starting from \lyb) from the optical depths of ions that have rest-wavelengths in the \lyb forest: \osix  and \cthree . We note that for the recovery of ions that require \hone subtraction, we do not mask the DLAs or fit out the wings for the \hone recovery (see Turner et al. 2014).

For ions that have a closely-spaced doublet (\osix  and \sifour ) we use the relative oscillator strengths to scale the optical depth of the weaker component to that of the stronger one, and take the minimum of the two optical depths at every redshift in order to correct for contamination.

Although \cfour  also has a closely-spaced doublet, due to its strength and the fact that it is located redwards of the \lya forest, most contamination comes from other \cfour systems. We therefore use a procedure where we iteratively subtract the optical depth of the weaker component at the position of the stronger component.

Although we attempt to correct for contamination, these corrections will generally not be perfect. In particular, we cannot correct for contamination of \osix by \hone \lya, which will cause our measurements of the \osix optical depth to be overestimates (however, on average it will affect all pixels in the forest equally). This contaminating \lya absorption, as well as residual absorption from other contaminating lines, set the limit down to which we can detect enhancements in the \osix optical depth.

3 Results

The first step in our analysis is to compute the median as a function of , which we will often denote \osix(\hone) for brevity. To do this, we take pixel pairs of and at each redshift , divide the pixels into 0.5 dex sized bins of , and compute the median of in each bin. We do this for both the full pixel sample (i.e. every pixel pair available from all 15 QSOs in the redshift range given by eq. 1, irrespective of the locations of the galaxies), and for pixels known to be located at small galactocentric distances, defined as those within3  \kmpsof the redshifts of galaxies with impact parameters  pkpc.

To avoid effects due to small number statistics, we do the following. For the full pixel sample, we divide the spectra into chunks of 5Å. We then compute the number of chunks that are sampled by each bin, and discard any bins that draw from fewer than five different chunks. For the small galactocentric distance pixel sample, we remove bins that do not have pixel contributions from at least five different galaxies. Finally, for both pixel samples we discard any bins containing fewer than 25 pixels in total.

Figure 1: Left panel: The \osix(\hone) relation for the full sample of pixels (black circles) and for pixels located at small galactocentric distances (red stars). Error bars indicate the uncertainty determined by bootstrap resampling the galaxies. The median values for the full sample are denoted by the horizontal and vertical dotted black lines. The \osix flat level, defined as the median value for all \osix pixels that have associated , is represented by the horizontal dashed grey line, and the shaded region shows its 1- error. The dashed line shows the optical depth relation for the full pixel sample after masking out regions within  \kmps of known galaxy redshifts, where we have considered galaxies with impact parameters of up to 2 pMpc. For a given \hone optical depth, we find a strong enhancement for the median \osix optical depth for pixels at small galactocentric distances. Centre and right panels: The same as the left panel except for \osix(\cfour) and \osix(\sifour). The fact that we observe the same trend when binning by different ions along the x-axis provides additional evidence that the effect we are seeing is real.
All pixels 0.04–0.18 pMpc
Table 2: The median pixel optical depth values plotted in Figure 1.

The derived \osix(\hone) relations are shown in the left panel of Fig. 1, where the black points and red stars indicate the relations for the full pixel sample and the small galactocentric distance sample, respectively. We also provide the data values shown in this figure in Table 2. The behaviour of the \osix(\hone) relation for the full pixel sample is consistent with previous measurements (e.g. Schaye et al., 2000a; Aguirre et al., 2008) and we briefly explain the observed characteristics here.

Focusing on the black circles, there are two distinct regions in \hone optical depth, separated by4 . For , increases with . This relation arises because a large number of these pixel pairs are probing regions which have been enriched by oxygen, and the value at each pixel is set by the median number density ratio of \osix to \hone.

For lower values of (below ), stays approximately constant, which indicates that the measured value is determined by residual contamination or noise and that the true, median value of is below this detection limit. In general, this asymptotic value of is slightly less than the median \osix optical depth of the full sample of pixels (the value is given in Table 1, and is indicated by the horizontal dotted black line in Fig. 1). We attempt to measure the constant value to which the \osix pixel optical depths asymptote, which we call the \osix flat level. As in Aguirre et al. (2008), we take this flat level to be the median of all pixels associated with .

Furthermore, we estimate a 1- error on this quantity by dividing each spectrum into 5 Å chunks and creating 1000 bootstrap resampled spectra. The flat level and the associated 1 error are denoted by the black dashed line and grey region, respectively. It is difficult to probe \osix optical depths below the flat level, as one becomes limited by contamination. Specifically, absorption from ions with rest wavelength less than that of \hone \lya are found bluewards of the QSO \lya emission line, and their recovery is limited by contamination from \hone and other metal lines. These metal transitions are less affected by the quality (S/N) of the spectra and are more sensitive to the method of recovery. On the other hand, for ions redwards of the QSO’s \lya emission, the median pixel optical depths are set mainly by the S/N and/or shot-noise, since the majority of pixels do not have detectable metal absorption. The two right panels of Figs. 4 and 5 in Aguirre et al. (2002) show the changes on the median levels for different recovery methods and S/N ratios for \osix (bluewards of the QSO \lya emission) and \cfour (redwards of the QSO \lya emission) using simulated spectra. These figures clearly demonstrate how the median level for \osix is more sensitive to recovery method while for \cfour changing the S/N has a greater effect.

Turning next to the red stars in the left panel of Fig. 1, we see a significant enhancement of at fixed for the small galactocentric distance pixels compared to the full pixel sample. Such a difference is not present if we consider larger impact parameter bins (not shown). We emphasize that although the full pixel sample is representative of random regions, it is composed of pixels both near and far from galaxies (and many of these galaxies are likely not detected in our survey). To demonstrate this, we examine the effect of masking out the regions in the spectra that are known to be near galaxies. We consider all galaxies in our sample with impact parameters  pMpc and with redshifts satisfying eq. 1, and mask out regions of  \kmps around these redshifts in all 15 of our QSO spectra. The resulting optical depth relation is shown as the dashed line in the left panel of Fig. 1, and we conclude that the full pixel sample relation is largely independent of the presence of regions proximate to detected galaxies.

Next, we vary the ion plotted along the x-axis. In the centre and right panels of Fig. 1, we show the relations for \osix(\cfour) and \osix(\sifour), respectively (the values are also given in Table 2). Just as for \osix(\hone), we observe significantly enhanced \osix optical depths at both fixed and fixed . The persistent enhancement of OVI absorption at fixed optical depth of three distinct ions (\hone, \cfour and \sifour) makes the individual detections with respect to each ion still more significant. Here we note that although the \osix(\hone) and \osix(\cfour) optical depth relations for which the galaxy positions have been masked out of the spectra show significant signal, this is not the case for \osix(\sifour). This suggests that strong \sifour absorption arises primarily near galaxies.

It is instructive to assess how sensitive the enhancement is to the chosen velocity range. We have examined the \osix(\hone) relation for velocity bins starting with  \kmps and increasing both velocity limits by increments of 10 \kmps(such that each cut spanned the same total velocity range of 340 \kmps). We found that for optical depth bins with , the enhancement is present up to a velocity range of  \kmps, and we take the midpoint of this bin,  \kmps, as an upper limit to the extent of the \osix enhancement. In Fig. 2 we show the relation for the  \kmps velocity bin (central panel), as well as higher and lower velocity cuts (right and left panels, respectively). For most of the optical depth bins with , the enhancement in \osix is only significant for velocities within  \kmps of the galaxy positions.

Although for the remainder of the analysis we will continue using the smallest velocity cut of  \kmps  we emphasize that for the lowest \hone optical depth bins, \osix is enhanced out to velocities of  \kmps, corresponding to times the typical circular velocities of the galaxies in our sample. While simulations predict that absorbers from galaxies below the detection limit can be projected to such velocities around their more massive counterparts (e.g., Rahmati & Schaye, 2014), in § 4 we determine that only a hot, collisionally ionized gas phase is observed out to these large velocities, which is certainly suggestive of outflows. Further comparisons with simulations will be required to fully disentangle outflow and clustering effects.

Figure 2: The same as Fig. 1 but we have overplotted the result of taking different velocity cuts when considering the small galactocentric distance points (blue squares). We have chosen velocity bins such that each cut spans the same velocity range of 340 \kmps, and find that the enhancement of \osix at fixed is present for out to the  \kmps velocity cut. Using the central value of this velocity range as our upper limit, we conclude that the enhancement persists out to  \kmps, which is times the typical circular velocities of the galaxies in our sample.

3.1 Are the observed differences in optical depth ratios real?

Figure 3: The probability density function of \hone (left), \cfour (centre) and \sifour (right) pixel optical depths for both the small galactocentric distance sample (red) and the full sample (black). For ease of comparison, we have shaded regions along the x-axis that are not included in the previous figures.
Figure 4: The same as Fig. 1 but instead of taking pixels within  \kmpsof the galaxy redshifts, the red stars show the effect of choosing regions within  \kmps of the galaxy redshifts. By excluding the positions directly around the galaxies, but still using pixels within  \kmps of the galaxy redshifts, we remove physical effects caused by the presence of the galaxy while probing the same spectral properties such as S/N, resolution, and distance from the QSO. Contrary to Fig. 1, the \osix absorption is not enhanced for the small galactocentric distance sample, which implies that the enhancement visible in that figure is not due to systematic differences in the spectral properties between the two samples.
Figure 5: The same as Fig. 1, except that we have overplotted the results of randomising the redshifts of the small galactocentric distance galaxies 1000 times. For each of these randomizations, we recompute \osix(\hone), \osix(\cfour) and \osix(\sifour), and in this figure we show the median, and the 84.1, 98.7, and 99.8 percentiles (corresponding to 0, 1, 2, and 3, respectively) from the full distribution of realizations. We find that the enhancement of \osix for fixed , and is approximately a 2–3- effect per bin.

Fig. 1 suggests that the gas near the sample galaxies has properties different from that in random regions with the same strength of \hone, \cfour, or \sifour optical depth. However, it is important to verify that the difference between the two pixel samples is not driven primarily by chance or systematic errors. In particular, limiting ourselves to only a few regions of the spectra could skew the results if, for example, these regions have different S/N, or inconsistent contamination levels due to being located at different redshifts.

First, we would like to be sure that the enhancement in the median \osix optical depths is not due to a small number of pixels. In Fig. 3 we show the pixel optical depth probability density functions (PDFs) for \hone, \cfour and \sifour, for both the small galactocentric distance and the full pixel samples (note that the x-axis ranges are larger than those shown in Fig. 1). Focusing on the small galactocentric distance pixels (red lines), one sees that the bins with enhanced median \osix (; and ) comprise the majority of the pixels. Thus, the enhancement in the median \osix optical depth cannot be attributed to a small number of pixels.

As a further test, in Fig. 4 we plot versus for the same shown in Fig. 1, but instead of using pixels within  \kmps of the galaxy redshifts, we look at regions further away, i.e.  \kmps. This velocity cut is near enough to the galaxy redshifts that within each impact parameter bin we are still looking at areas of the spectra with the same S/N and contamination characteristics, but far enough to avoid regions that may be associated with the galaxies. If the enhancement of the \osix optical depth at fixed \hone, \cfour, or \sifour that we detect near galaxies were caused by systematic differences in the spectral properties of the two samples, then we would expect to see the enhancement to a similar significance in these figures. However, in every case where we previously saw an enhancement in the median \osix optical depth for the pixels at small galactocentric distance, the effect is completely removed in Fig. 4.

Next, we examine whether the optical depth differences are consistent with random fluctuations. To do this, we take the galaxies from the small galactocentric distance sample, randomize their redshifts 1000 times, and calculate \osix(\hone), \osix(\cfour) and \osix(\sifour) for every random realization. In Fig. 5 we show the one-, two-, and three- percentiles that result from this procedure. From this, we measure the enhancements seen in \osix(\hone), \osix(\cfour), and \osix(\sifour) to be approximately a 2–3- effect per bin.

Additionally, because we do not observe differences in the \osix(\hone), \osix(\cfour) and \osix(\sifour) relations for galaxies with impact parameters  pkpc and the full pixel sample, here we check whether the \osix enhancement for the galaxies with impact parameters  pkpc can be attributed to properties other than the galaxy distance to the QSO sightline. We directly test whether the following three characteristics in every impact parameter bins are consistent with the full sample: (1) the galaxy redshifts (lower galaxy redshifts mean that \osix is more contaminated by \hone), (2) the galaxy velocity distance from the QSO \lya emission (to rule out QSO proximity effects) and (3) the S/N of the spectral regions (smaller S/N will bias the optical depth estimation high).

We have measured the -values resulting from a 2-sample Kolmogorov-Smirnoff test, comparing the galaxy redshift, median S/N within  \kmps from the redshift of the galaxy, and redshift difference between the galaxy and the QSO, between the galaxies in the small galactocentric distance sample and the full galaxy sample. In every instance we find -values greater than 0.1, which is consistent with the null hypothesis that the two samples are drawn from the same distribution.

Finally, in Appendix A we have tested how the \osix contamination correction affects our results. We determined that although the inferred values of the \osix optical depths are sensitive to changes in the correction procedure, the enhancement of \osix optical depths at fixed \hone the \osix for pixels at small galactocentric distances compared to random regions is unchanged, and remains significant irrespective of the contamination correction details. Thus, we conclude that the observed optical depth differences are neither due to chance nor to systematic variations in the spectral characteristics in either the QSOs or the galaxies.

Figure 6: Theoretical optical depth ratios () as a function of temperature and hydrogen density from CLOUDY modelling. For this figure we have assumed , however we allow this value to vary according to eq. 3 when comparing with the observations. The vertical dashed line denotes  K, which is the temperature that we assume in the case of PIE. Furthermore, under the assumption of PIE we can convert hydrogen number densities to \hone optical depths with eq. 5, and we show corresponding values on the right-hand y-axis. For CIE, we set can use the ratios to set lower limits on the metallicity by using the maximum ratio from this temperature-density plane. This maximum occurs at  K and is marked by the small black circle at .

4 Ionization models

In this section, we will investigate the physical origin of the difference in \osix observed between regions near our galaxies and in random locations. We focus on \osix(\hone) because binning by enables a more straightforward physical interpretation than binning by metal ion optical depth, due to the tight correlation between \hone absorption and gas density in photoionized gas. We consider three possible scenarios and examine their plausibilities. The first hypothesis we consider is enriched photoionized ( K) gas, where the enhanced \osix(\hone) near galaxies might be explained by higher gas metallicities. Second, we test the idea that ionizing radiation from the nearby galaxy could be responsible for the increase in \osix(\hone). Finally, we consider whether our observations can be explained by a hot, collisionally ionized enriched gas phase near galaxies. We argue that of these three, the first scenario (enriched photoionized gas) can account for pixels with , while only the third explanation (the presence of hot, enriched gas) is plausible for pixels with . Of course, it is important to note that every \hone bin likely contains a mixture of pixels from different gas phases and ionization sources; the final behaviour is simply determined by the dominant phase.

4.1 Photoionization by the background radiation

The optical depth of \hone is believed to be a good tracer of the photoionized gas density, even on an individual pixel basis (Aguirre et al., 2002). Hence, if the gas probed were predominantly photoionized and if the abundance of oxygen depended only on gas density, we would not expect to see any difference between the \osix(\hone) relations of all pixels and those known to lie near galaxies in Fig. 1. Since a clear difference is observed, we postulate that this could be caused by an increase in the oxygen abundance near galaxies at fixed gas densities.

To test this idea, we turn to ionization modelling using CLOUDY (Ferland et al., 2013, version 13.03). Our setup involves a plane-parallel slab illuminated uniformly by an ionizing background, along a grid of varying temperatures and hydrogen densities, that covers both photoionization equilibrium (PIE) and collisional ionization equilibrium (CIE, discussed in § 4.3). For the fiducial case we use the ionizing background from Haardt & Madau (2001) including contributions from both quasars and galaxies, normalized to match the metagalactic \hone photoionization rate, , from Becker et al. (2007).5 However, the shape and normalization of the background is subject to large uncertainties, and we explore the effect of varying them in Appendix C.6

Element Ref.
C 1.
O 2.
Si 3.
Table 3: Solar abundances used in this work, taken from CLOUDY 13. References are 1. Allende Prieto et al. (2002); 2. Allende Prieto et al. (2001); 3. Holweger (2001).

After obtaining individual ionization fractions from the above models as a function of density, , and temperature, , we can relate metallicity to optical depths using the following equation:


where and are the oscillator strengths and rest wavelengths of the transitions, and we assume solar abundances from CLOUDY 13 (listed in Table 3). The resulting optical depth ratio contours, as a function of temperature and hydrogen number density and assuming solar metallicity, are shown in Fig. 6. At temperatures  K, photoionization dominates and the contours are only weakly dependent on temperature, while at higher temperatures and for sufficiently high densities collisional ionization dominates and the contours are independent of the density.

Figure 7: Left panel: The ratio of \osix to \hone optical depths as a function of . To compute these points, we used the data from Fig. 1. First, we subtracted the \osix flat level from the median \osix optical depths, in order to correct for contamination. After dividing the resulting values by their corresponding \hone optical depths, we obtain the points plotted here. The black vertical dotted line corresponds to the median value of all \hone pixel optical depths, while the diagonal line shows the median value of all \osix pixel optical depths divided by the given value of along the x-axis. Right panel: Theoretical ratios of \osix to \hone as a function of taken from CLOUDY modelling (see Fig 6). First, in the case of PIE, we assume that the gas has a temperature of  K, and convert \hone optical depths to a hydrogen number density using eq. 5 (the corresponding values are shown on the upper y-axis). The blue curves show profiles taken along the vertical black dashed line in Fig. 6, where the different line styles demonstrate the effect of varying the metallicity. We also consider CIE, in which case cannot estimate the density from . Instead, we can use the maximum theoretical ratio (indicated by the black circle at Figure 6) to set a lower limit on the metallicity, and these are shown by the horizontal orange lines.

We would like to compare the observed optical depth ratios with those predicted from CLOUDY in order to estimate the gas metallicity. To compute the observed ratios, we begin with the points in the left panel Fig. 1. In order to correct for the presence of residual contamination from absorption by species other than \osix, we subtract the flat level (i.e. the asymptotic value of reached for small and indicated by the horizontal, dashed line in Fig. 1) from all values. Following Schaye et al. (2003), to be conservative we added the error on the flat level linearly (rather than in quadrature) to the errors on , and finally we divide every point by its corresponding . The results of this calculation are plotted in the left panel Fig. 7.

Next, since we are assuming that the gas is photoionized, we can transform the \hone optical depths to gas densities, as was for example done in Rakic et al. (2012). The relation between column density and pixel optical depth at the line centre, , is:


Here, is the oscillator strength, is the rest-wavelength of the transition, and is the Doppler line width, which we set to  \kmps, the typical value measured by Rudie et al. (2012b). In order to apply this to our data, we must assume that the \hone pixel optical depths are close to the line centre . This is not a bad approximation, because the steepness of the \hone column density distribution function implies that it is more likely to be close to the maximum optical depth of a weaker absorber than it is to be in the wing of a stronger absorber. Indeed, Fig. 11 of Rakic et al. (2012) shows that the above methodology leads to very close agreement with the simulation result of Aguirre et al. (2002).

To convert to a hydrogen number density, we turn to Schaye (2001), who derive a relation between density and column density by assuming the absorbers are gravitationally confined gas clouds, which implies that they have sizes on the order of the local Jeans length. Combining their eq. 8 with eq. 4, we obtain:


where we have assumed the metagalactic photoionization rate to be , the value from Becker et al. (2007), and we have taken the temperature to be  K, typical for a moderately overdense IGM region (e.g. Schaye et al., 2000b; Lidz et al., 2010; Becker et al., 2011; Rudie et al., 2012a). The result is insensitive to the precise value of the temperature, as long as the gas is predominantly photoionized, as assumed by Schaye (2001). We have also assumed a gas fraction close to the universal value of .

Thus, we can use eq. 5 to convert between \hone optical depth and hydrogen number density. Next, we consider the contour along a constant temperature of  K in Fig. 6, from which we obtain theoretical optical depth ratios as a function of . With eq. 5, we can convert the dependence to corresponding values of . The result of doing so is shown by the blue lines in Fig. 6, where the different line styles demonstrate the effect of varying the metallicity using eq. 3.

Figure 8: Left panel: Metallicity as a function of \hone pixel optical depth. The metallicity was inferred from the ratio of to (Fig. 7) assuming photoionization equilibrium (PIE,  K) as described in the text. On the upper axis we show converted to an overdensity using the relation for photoionized gas that we assumed in the ionization models. Random regions (black circles) show a positive correlation between metallicity and density, while at small galactocentric distances (red stars) for the metallicity increases inversely with density. Such an inverted metallicity-density relation is not physically expected and indicates that the assumption of PIE is probably incorrect. Right panel: Lower limits on the metallicity as a function , inferred again from the measured optical depth ratios but now assuming collisional ionization equilibrium (CIE). The size of each arrow is proportional to the size of the error bar on the measured optical depth ratio. The inferred limits on the metallicity are lower than for PIE, but still imply that at least some of the gas near galaxies is substantially enriched.

Hence, by interpolating between the blue curves in the right panel of Fig. 7, we can estimate a metallicity (assuming PIE) for each optical depth ratio in the left panel. The outcome of this procedure is plotted in the left panel of Fig. 8. To aid the interpretation, we indicate the baryon overdensity inferred from along the upper x-axis (again, we can only convert \hone optical depth to overdensity under the assumption of PIE).

Focusing first on the black points that correspond to the full pixel sample, we find that metallicity increases with overdensity , and is in agreement with previous pixel optical depth studies (e.g. Schaye et al., 2003; Aguirre et al., 2008). Next, examining the values derived from pixels with small galactocentric distance (red stars), for pixels with we infer the same metallicity-density relation as for the full pixel sample, but with metallicities that are  dex higher. We conclude that for these \hone optical depths, photoionized gas that is enriched with respect to random regions of the same density is consistent with our observations of \osix(\hone) at small galactocentric distances.

However, turning to pixels with , it is clear that the relation between metallicity and density is inverted. In this regime, the metallicity appears to increase with decreasing overdensity, with the highest metallicities found at the lowest overdensities, even reaching supersolar values. We do not consider this enrichment pattern to be physically plausible. Although some regions of relatively low density can be highly enriched, they usually arise in hot superbubbles. Furthermore, an underdense gas phase close to galaxies with  K and supersolar metallicities is not consistent with predictions from cosmological hydrodynamical simulations (e.g. Ford et al., 2013; Shen et al., 2013). The above suggests that one or more of our assumptions must be incorrect for the small galactocentric distance pixel sample.

4.2 Photoionization by stars in nearby galaxies

Figure 9: The data points are identical to those in the left panel of Fig. 1. The curves show the result for the full pixel sample (black circles) after dividing the \hone optical depths by the boost factor in the \hone photoionization rate expected due to ionizing radiation from the nearby galaxy for the relative escape fractions, , indicated in the legend and for a distance to the source equal to the median impact parameter of the pixels in the small galactocentric distance sample. Even for %, which is unrealistically high, the curve falls short of the red stars, implying that photoionization by radiation from the observed galaxies cannot explain the observed enhancement in \osix/\hone near galaxies.

In this section, we test whether the enhancement of for a fixed near galaxies can be explained by increased photoionization by radiation from the stars in the nearby galaxies. At sufficiently small galactocentric distances, we would expect the mean hydrogen ionization rate to be dominated by galactic radiation rather than by the metagalactic background radiation (e.g. Schaye, 2006; Rahmati et al., 2013), which would reduce the \hone optical depth at fixed density. Furthermore, since stars emit very little radiation with photon energies  Ryd, we do not expect the \osix optical depth to be modified by local stellar sources, given that the ionization energy for \ofive is about 8.4 Ryd. Hence, if photoionization by nearby galaxies were important, we would expect the red stars in Fig. 1 to be shifted upwards relative to the black points and the red stars in Fig. 7 to be shifted towards the top-left relative to the black points. This is qualitatively consistent with the observations, which suggests that the difference between the two samples may be due to photoionization by local sources rather than due to the presence of hot, enriched gas. However, as we will show below, quantitatively this scenario does not work out.

We first need to evaluate the photoionization rate due to local sources, , at  Ryd,


where is the Planck constant, is the incident flux, and is the photoionization cross section, taken to be where  cm. To estimate the incident flux from a nearby galaxy, we assume that the galaxy is an isotropically emitting point source:


where is the galaxy luminosity, is the escape fraction of photons at frequency , and is the proper distance from the point source. For the luminosity, we use measurements of the KBSS galaxy luminosity function from Reddy & Steidel (2009), who obtained  mag. This corresponds to  erg s Hz at a rest wavelength of  Å, where is the escape fraction at this wavelength and is the Hubble constant in units of 70 kmsMpc.

We then use the fact that spectral synthesis models corresponding to typical KBSS galaxy metallicities of can be well-approximated by a blackbody curve with  K (Steidel et al., 2014). Integrating this function numerically, we obtain:


where we take  pkpc because it is the median impact parameter of the galaxies that comprise our small galactocentric distance sample. Note that this choice is conservative since it assumes that the absorption within of the galaxy arises in gas that is precisely at the distance of the galaxy. In reality, some of this gas will be in front or behind the galaxy and will be greater than the impact parameter.

The remaining unknown in our estimate of the photoionization rate from local galaxies comes from the ratio of the escape fractions at 1 Ryd and at the observed wavelength of 1700 Å, . The literature contains many values for the escape fraction of Lyman continuum photons relative to that of non-ionizing UV continuum photons, . Values range from % to % (e.g. Steidel et al., 2001; Shapley et al., 2006; Iwata et al., 2009; Nestor et al., 2013; Mostardi et al., 2013; Cooke et al., 2014), although typical values are closer to the lower limit of this range.

Since we already expect  s from the extragalactic background (Becker et al., 2007), adding a contribution from nearby galaxies would boost the \hone photoionization rate by a factor of for an escape fraction % at the median impact parameter and would hence suppress the \hone optical depths by the same amount. To test whether ionizing radiation from nearby galaxies can explain the difference between the \osix(\hone) relations of the small galactocentric radii and the random pixel samples, we have re-calculated the \osix(\hone) relation for random regions after dividing all the \hone pixel optical depths by the factor by which local sources boost the \hone ionization rate for various values of the relative escape fraction. The results are compared with the observed relation for the small galactocentric radii sample in Fig. 9. For % the effect of local sources is not strong enough to reproduce the observations. Even for the highly unrealistic case that %, the boost in the \hone ionization rate is insufficient to completely account for the observed enhancement of \osix small galactocentric distances for .

On the other hand, while for it appears that relative escape fractions of % may be able to explain our observations, we reiterate that average observed relative escape fraction values are usually much lower than this for galaxies such as ours (e.g., Shapley et al., 2006). Furthermore, we note that by considering only the transverse rather than 3-dimensional distance to the galaxy, our estimate of the \hone photoionization rate boost is strictly an upper limit, and the true strength of the effect is almost certainly smaller.

As a final test, we have taken the blackbody spectrum used to approximate the local galaxy radiation and used it as input, along with the extragalactic background, into CLOUDY. We generated temperature-density planes showing optical depth ratios (as in the right panel of Fig. 7), as well as and , and compared these with the output using only the extragalactic background. For the \osix(\hone) relation, we indeed found that the predicted ratios for photoionized gas at  K were higher. However, as expected due to the fact that our assumed spectrum does not contain many photons with energies above 1 Ryd, for the observed ratio values there was almost no discernible difference between the \osix(\cfour) and \osix(\sifour) relations derived with and without the contribution from nearby galaxies.7 Hence, if local sources were responsible for the enhancement in \osix(\hone), we would not expect to see any difference between the small galactocentric distance and random pixel samples for the \osix(\cfour) and \osix(\sifour) relations. This is not consistent with the fact that we find significantly enhanced \osix for both fixed \cfour and fixed \sifour much like we do for the \osix(\hone) relation (see the centre and right panels of Fig. 1). Thus, we conclude that enhanced photoionization due to the proximity of the observed galaxies cannot explain the enhancement of \osix at fixed \hone that we observe for the small galactocentric distance pixel sample.

In the above scenario, we only consider photoionization by stars in the nearby galaxies because there is no evidence that our galaxies contain AGN. However, it was recently pointed out by Oppenheimer & Schaye (2013a) that because of the long recombination times for metals at densities typical of the circumgalactic medium, ions like \osix may remain out of ionization equilibrium long after the AGN episode is over. In fact, for reasonable AGN duty cycles, the authors argue that much of the \osix detected in quasar spectra resides in such fossil AGN proximity zones. On the other hand, in regions where the equilibrium neutral hydrogen fraction is low, \hone will equilibrate nearly instantaneously after the AGN turns off. Hence, in this scenario \hone would be photoionized and in equilibrium, while \osix is out of equilibrium and enhanced because of the fossil AGN proximity effect. We would then expect the enhancement in \osix to be strongest in regions of low \hone absorption, since such regions correspond to low densities and thus long recombination times. This is in qualitative agreement with our observations, although Oppenheimer & Schaye (2013a) show that at very low density oxygen will become so highly ionized that \osix is suppressed, which would be inconsistent with our findings.

More generally, since (fossil) AGN proximity and higher gas temperatures both tend to increase the abundance of more highly ionized species relative to ions with lower ionization potentials, many of the predictions of these two scenarios will be qualitatively similar, although it remains to be seen whether the fossil proximity effect can work quantitatively. One important difference, however, is the expected widths of the absorption lines. If the oxygen is hot enough for \osix to be collisionally ionized, then the absorption lines will be broader than if the gas were photoionized. We intend to measure the widths of the metal absorption lines and to model the fossil AGN proximity effect in future work.

4.3 Collisionally ionized gas

It may be that at small galactocentric distances, much of the \osix absorption arises in collisionally ionized rather than photoionized gas. Indeed, the behaviour of the \osix(\hone) relation for the near-galaxy pixels with points to collisional ionization playing a role. In general, we would not expect to see very much \hone associated with collisionally ionized \osix, as the \hone fraction is very low at  K. The behaviour of the [O/H]() relation that we inferred under the assumption of photoionization (Fig. 8) is rather suggestive, as the relation displays larger departures from the random regions with decreasing .

If we are in fact probing collisionally ionized gas, then we cannot measure its metallicity as was done in § 4.1 – first because we are unable to estimate the gas density from the \hone optical depth, and second because there is no natural equilibrium temperature. However, it is possible to put a lower limit on the metallicity by selecting the temperature and density that would maximize . In Fig. 6, we have marked the temperature and density where this theoretical maximum is reached by a small black circle. Then, in the right-hand panel of Fig. 7, the horizontal orange lines show the ratio at this point for different metallicities.

Thus, irrespective of \hone optical depth (which has no correspondence with density under the assumption of CIE) we can convert the observed optical depth ratios in the left-hand panel of Fig. 7 to lower-limits on the metallicity, and we show the result in the right panel of Fig. 8. As expected, the limits obtained are lower than the metallicities inferred when assuming photoionization. Nevertheless, some of the lower limits are sufficiently high to be interesting. The lower limit approaches at the lowest bins in \hone optical depth, where .

Given the reasonable metallicity limits in Fig. 8 combined with the fact that the previous two scenarios (photoionization by the extragalactic background or by local stellar radiation) cannot account for the \osix enhancement for , we conclude that collisional ionization is the dominant ionization source for a significant fraction of the small galactocentric distance pixels. We emphasize that in collisionally ionized gas, is no longer a good estimator of density, and we could potentially be probing very high overdensities.

We also examined the \osix/\cfour and \osix/\sifour ratios in Appendix B. We conclude that \osix pixels near galaxies are more consistent with a hotter gas phase than the full pixel sample, and furthermore we find that the strength of this trend increases inversely with \cfour and \sifour optical depths, as might be expected for collisionally ionized gas at  K where these other ions are not abundant. However, since the results are not constraining, the details are left to the Appendix.

5 Discussion and conclusions

We have used absorption by \hone and metals to study the physical conditions near star-forming galaxies in the fields of 15 hyper-luminous background QSOs that have been observed with Keck/HIRES as part of the Keck Baryonic Structure Survey (KBSS). We focused on 21 galaxies with impact parameters  pkpc and isolated the pixels of the quasar spectra that are within  \kmps of the accurate galaxy redshifts provided by the KBSS. In Turner et al. (2014) we showed that the absorption by \hone and various metal ions is strongly enhanced in these circumgalactic regions. The fact that both \hone and the metals are enhanced raises the question whether the observed increase in the metal absorption merely reflects the presence of higher gas densities near galaxies or whether it implies that the gas near galaxies has a higher metallicity at fixed density or a different temperature from gas in random regions.

To address this question, we measured the pixel optical depths of \osix as a function of \hone, \cfour and \sifour, and compared the results for the pixels located at small galactocentric distances to the full pixel sample, which is representative of random regions. Our main result is the detection of a 0.3–0.7 dex enhancement (which reaches its maximum at ) in the median optical depth of \osix at fixed for the small galactocentric distance sample when compared with the full pixel sample (Fig. 1). We verified that this enhancement, which we detected at 2–3 confidence per logarithmic bin for , , and , is not due to differences in the redshift distribution or the quality of the quasar spectra between the small galactocentric distance and full pixel samples.

We proposed and tested three different hypotheses that may explain the observed enhancement of \osix(\hone) near galaxies: (1) the gas is photoionized by the extragalactic background but has a higher metallicity at fixed density; (2) the gas is more highly ionized at fixed density by radiation from stars in the nearby galaxy; and (3) the enriched gas is hot and collisionally ionized.

To test scenario (1), we used CLOUDY ionization models and the relation between and density for photoionized, self-gravitating clouds from Schaye (2001) and Rakic et al. (2012) to convert the observed optical depth ratios into a metallicity-density relation, assuming K, as expected for photoionized gas. We found that the full pixel sample gave a physically plausible metallicity-density relation that is consistent with previous studies which did not have information on the locations of galaxies. Furthermore, the same metallicity-density relation (but shifted up by  dex in metallicity) is also present for small galactocentric distance pixels that have . Thus, for high \hone optical depths, the enhancement of \osix at fixed is consistent with arising predominantly from enriched, photoionized gas. However, making this same assumption of PIE for resulted in an [O/H] versus overdensity relation that cannot be easily explained, as [O/H] was found to increase strongly with decreasing overdensity, implying supersolar metallicities for underdense gas (left panel of Fig. 8). We therefore concluded that while photoionization by the background radiation is a plausible scenario for the random regions, it cannot account for the observed enhancement of \osix(\hone) near galaxies.

In scenario (2) \osix is enhanced at fixed \hone because radiation from stars in the nearby galaxy suppresses \hone, while leaving \osix unchanged. To test this explanation, we estimated the \hone photoionization rate due to the galaxies. However, we found that only under the unrealistic assumptions that the relative escape fraction is 100% and that the 3-D distance between the gas and the galaxy is equal to the impact parameter, can the flux of ionizing radiation from the galaxies explain the observed increase in near galaxies. Reducing the relative escape fraction to a still conservative value of 10% rules out this scenario as a dominant ionization source (Fig. 9). Furthermore, we found that such a galaxy proximity effect is predicted to have a much smaller effect on \osix(\cfour) and \osix(\sifour) than on \osix(\hone), which is inconsistent with our observations.

Contrary to photoionization by either the extragalactic background or local stellar radiation, scenario (3) can explain the enhancement in \osix(\hone) near galaxies for pixels with . If a substantial fraction of the enriched gas near galaxies is sufficiently hot for \osix to be collisionally ionized, i.e. K, then we can account for the observations. By assuming the maximum \osix/\hone ratio reached in collisional ionization equilibrium, we converted the observed optical depth ratios into lower limits on the metallicity, finding for gas with weak \hone absorption (right panel of Fig. 8). Indeed, this measurement is supported by other characterizations of KBSS galaxy properties. Rudie et al. (2012b) found higher temperatures for fixed near the KBSS galaxies, while Steidel et al. (2014) measured \htwo region metallicities of  Z, which could serve as a possible upper limit to the value presented here (although if metal-enriched winds drive most of the metals out of the star-forming regions, it is possible that the circumgalactic gas may have a higher metallicity than galactic \htwo regions). Furthermore, the inferred metallicities and temperatures of  K are in agreement with the predictions of van de Voort & Schaye (2012) for galaxies with the masses and redshifts similar to ours.

In summary, we favour the conclusion that our galaxies are surrounded by hot (K) gas of which a substantial fraction must have metallicity of solar. Furthermore, we find that this metal-enriched phase extends out to  \kmps of the galaxy positions (Fig. 2), which corresponds to times the halo circular velocities. Because of the relatively high temperature that requires shock-heating, the large velocity range extending far outside the haloes, and high metallicity, we conclude that we have detected hot, metal-enriched outflowing gas. Future comparisons with hydrodynamical simulations, considering ion ratios as well as the kinematics and line widths, will provide strong constraints on models of galaxy formation and may provide further insight into the interpretation of our observations.


We are very grateful to Milan Bogosavljevic, Alice Shapley, Dawn Erb, Naveen Reddy, Max Pettini, Ryan Trainor, and David Law for their invaluable contributions to the Keck Baryonic Structure Survey, without which the results presented here would have not been possible. We also thank Ryan Cooke for his help with the continuum fitting of QSO spectra. We gratefully acknowledge support from Marie Curie Training Network CosmoComp (PITN-GA-2009- 238356) and the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC Grant agreement 278594-GasAroundGalaxies. CCS, GCR, ALS acknowledge support from grants AST-0908805 and AST-13131472 from the US National Science Foundation. This work is based on data 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 NASA, and was made possible by the generous financial support of the W.M. Keck Foundation. We thank the W. M. Keck Observatory staff for their assistance with the observations. We also thank the Hawaiian people, as without their hospitality the observations presented here would have not been possible.

Appendix A Correction of OVI contamination

Figure 10: The same as the left panel Fig. 1, but varying the method of \osix contamination correction. The modifications shown here are the result of not subtracting any \hone (left panel), not taking the doublet minimum (centre panel), and not performing any correction at all (right panel). This figure demonstrates that the qualitative behaviour of the enhancement of \osix optical depth at fixed \hone for pixels at small galactocentric distance is entirely independent of the \osix contamination correction procedure.

Here we examine how modifying the \osix contamination correction affects the resulting \osix(\hone) relation. It is possible that the small galactocentric distance pixel sample suffers from more \hone contamination than the full pixel sample, given the larger on average \hone optical depth values (see the left panel of Figure 3) and the proximity of \lyb to \osix ( \kmps). We postulate that if the enhancement of \osix for fixed \hone optical depths at small galactocentric distances were largely due to uncorrected \hone \lyb contamination, modification of the \osix contamination would have a notable effect on this enhancement.

To test this, in Fig. 10 we show the \osix(\hone) relation as in the left panel of Fig. 1, but with changes to the \osix contamination correction. In the fiducial case, we subtract 5 orders of the Lyman series of \hone beginning with \lyb, and then take the minimum of the \osix doublet optical depths. In figure 10, we first show \osix(\hone) after not subtracting any \hone but still taking the doublet minimum (left panel), not taking the doublet minimum but still subtracting \hone (centre panel), and not performing any correction at all (right panel). Although the absolute values of the apparent \osix optical depths increase when fewer contamination corrections are performed, the significance of \osix optical depth at fixed \hone for pixels at small galactocentric distance compared to the full pixel sample remains unchanged, indicating that our results are not sensitive to the \osix contamination correction procedure.

Finally, we perform a simple calculation to demonstrate that the \hone \lyb contamination is not predicted to have a large effect. To begin, we would like to estimate the typical \lyb strength at a distance  \kmps away from galaxies. For this, we use our measurement of the median optical depth as a function of distance along the galaxy LOS, as was done in the right panels of Fig. 6 from Turner et al. (2014), but extending out to  \kmps. We find that the \hone \lya optical depths asymptote to the median value of all pixels (0.051), and using the relative oscillator strengths of \lya and \lyb we can convert this to an expected median value of \hone \lyb (0.0082). Adding this to the median value of all \osix pixels (0.023), we obtain 0.031. This corresponds to in the left panel of Figure 1, and is significantly less than the observed \osix optical depth enhancement.

Appendix B OVI(CIV) and OVI(SiIV)

Figure 11: Left column: Median optical depth ratios derived from the \osix(\cfour) (top) and \osix(\sifour) (bottom) relations in Fig. 1. The diagonal lines show the median value of all \osix pixel optical depths divided by the given value of and along the x-axis, and demarcates a conservative detection level. Right column: Predicted optical depth ratios from CLOUDY.

In the centre and right panels of Fig. 1, we investigated whether the \osix(\cfour) and \osix(\sifour) relations depend on galactocentric distance. First, we note that we see similar behaviour for \osix(\cfour) and \osix(\sifour) as for \osix(\hone) – that is, we find a higher at fixed for pixels near galaxies compared to random locations, and the magnitude of this enhancement increases with decreasing .

As in § 3, we can use CLOUDY to model the optical depth ratios as a function of temperature and density. These temperature-density planes, along with the observed ratios, are shown in Fig. 11. We note that if we assume solar relative abundances, then the optical depth ratios determined from the CLOUDY models are fixed and do not depend on metallicity (unlike for optical depth ratios with in the denominator).

Looking first at \osix(\cfour), for pixels near galaxies and at random locations have similar optical depth ratios ranging from to , which correspond to a maximum temperature of  K. On the other hand, at lower \cfour optical depths, the ratios derived from pixels near galaxies are as high as , and are not inconsistent with the collisionally ionized region at  K.

Similarly, for \osix(\sifour) the optical depth ratios derived from pixels with small galactocentric distances have higher values, and therefore higher temperature upper limits, than in the random regions. They reach values as high as , which again corresponds to a maximum temperature of  K.

Taken together, the sample of pixels at small galactocentric distance show the same trend for \osix(\cfour) and \osix(\sifour), discrepant from that of random locations: increases inversely with and . If the conditions at near galaxies are most favourable to the collisional ionization of \osix (i.e. temperatures  K), then we would expect to find this oxygen in regions with less \cfour and \sifour, and such a scenario is certainly not inconsistent with the above relations.

Appendix C Ionization background

The spectral shape of the ionizing background radiation is a very large source of uncertainty for ionization modelling. In this section, we explore the impact that changing this background has on our results. We experiment with four different models, shown in Fig. 12, which have all been normalized to have the same \hone photoionization rate, which we take from (Becker et al., 2007) to be  s at . For our fiducial background, we use the Haardt & Madau (2001) quasar+galaxy model (HM01) and vary it in two ways. First, we use the quasar-only model, which increases the intensity above 1 Ryd (HM01 Q-only). We also invoke a model where the intensity above 4 Ryd is reduced by 1 dex, in order to simulate the absence of helium reionization (HM01 cut). Finally, we compare with the updated background from Haardt & Madau (2012, HM12).

Figure 12: Comparison of different extragalactic ionization backgrounds.
Figure 13: [O/H] derived from the \osix(\hone) relation, for four different extragalactic backgrounds (as indicated in each panel), under the assumption of PIE.

As a first test, we examine the values of [O/H] derived when we assume PIE with  K for the four different backgrounds in Fig. 13. We find that for the first three backgrounds (HM01, HM12, and Q-only) the results are in fairly good agreement. In general, the metallicity-density relation is steepest for HM01, and spans the largest range of [O/H] values.

Of the four backgrounds considered, the one salient outlier is the HM01 cut background (bottom-right panel). First, for the small galactocentric distance points, invoking this background produces reasonable (i.e. at least less than solar) metallicities for . However, the qualitative behaviour of inferred [O/H] as density decreases is still present. Furthermore, for , the inferred metallicities of the full pixel sample become unrealistically high (up to [O/H], while Steidel et al. 2014 find  Z for this galaxy sample), suggesting that this background is not a realistic choice.

Figure 14: The predicted values of for the HM01 (black solid contours) and HM01 cut (blue dashed contours) backgrounds.

Another significant change that would result from using the HM01 cut model is shown in Fig. 14, where we compare the the values of as a function of temperature and density for both the HM01 (solid black contours) and HM01 cut (dashed blue contours) backgrounds.

The observed ratios span values that range from to for the small galactocentric distance sample. For the HM01 background, this corresponds to gas that has either to  cm at  K, or  cm at  K. However, for HM01 cut, at the typical photoionized gas temperature of  K,  cm for the range of observed optical depths, which corresponds to underdense values (i.e. ). Although the small galactocentric distance pixel sample could be reconciled with temperatures of  K, it is unlikely that the full pixel sample is typically probing such hot gas, or densities below the cosmic mean. We conclude that the HM01 cut background is in tension with ratios of derived from the full pixel sample.


  1. pagerange: Detection of hot, metal-enriched outflowing gas around z2.3 star-forming galaxies in the Keck Baryonic Structure SurveyC
  2. Although this work does not use direct galaxy detections, the authors inferred that their strongest absorber sample probes regions defined as circumgalactic in Rudie et al. (2012b) in percent of cases. This was determined by using high-resolution QSO spectra from fields where the redshifts of Lyman break galaxies with impact parameters  kpc are already known. The QSO spectra were degraded to match the resolution of the spectra used in Pieri et al. (2014), and a correction for the expected volume density of Lyman break galaxies from Reddy et al. (2008) was applied.
  3. This velocity interval was chosen because (as explained in more detail in Rakic et al. 2012 and Turner et al. 2014) it is the scale over which the optical depths are smooth in the LOS direction.
  4. As noted in Aguirre et al. (2008), because the \osix pixel relations are not as strong as for \cfour, we fix the value of by hand rather than using functional fits. We use the same values as in Aguirre et al. (2008) of = 0.1 when pixel pairs are binned based on the optical depths of transitions that fall blueward of the QSO’s \lya emission (\osix, \hone) and 0.01 for those falling redward (\cfour, \sifour).
  5. Measurements of at vary between in Faucher-Giguère et al. (2008) up to in Becker & Bolton (2013). We choose the intermediate value of at taken from the fitting formulas of Becker et al. (2007).
  6. We note that the assumption of PIE may break down as non-equilibrium effects become important in cooling gas at temperatures  K. The presence of these effects causes collisional ionization to occur at lower temperatures compared to in equilibrium, although once an extragalactic background is included the impact of the effect of non-equilibrium cooling become less important (Oppenheimer & Schaye, 2013b).
  7. Invoking galaxy spectra with a soft X-ray component could potentially lead to more oxygen ionized to \osix (e.g., Cantalupo, 2010), however these spectral models remain uncertain and we leave testing of this scenario to a future work.


  1. Adelberger, K. L., Shapley, A. E., Steidel, C. C., et al. 2005a, \apj, 629, 636
  2. Adelberger, K. L., Steidel, C. C., Pettini, M., et al. 2005b, \apj, 619, 697
  3. Adelberger, K. L., Steidel, C. C., Shapley, A. E., et al. 2004, \apj, 607, 226
  4. Adelberger, K. L., Steidel, C. C., Shapley, A. E., & Pettini, M. 2003, \apj, 584, 45
  5. Aguirre, A., Dow-Hygelund, C., Schaye, J., & Theuns, T. 2008, \apj, 689, 851
  6. Aguirre, A., Schaye, J., & Theuns, T. 2002, \apj, 576, 1
  7. Allende Prieto, C., Lambert, D. L., & Asplund, M. 2001, \apjl, 556, L63
  8. —. 2002, \apjl, 573, L137
  9. Aracil, B., Petitjean, P., Pichon, C., & Bergeron, J. 2004, \aap, 419, 811
  10. Becker, G. D., & Bolton, J. S. 2013, \mnras, 436, 1023
  11. Becker, G. D., Bolton, J. S., Haehnelt, M. G., & Sargent, W. L. W. 2011, \mnras, 410, 1096
  12. Becker, G. D., Rauch, M., & Sargent, W. L. W. 2007, \apj, 662, 72
  13. Bergeron, J., Aracil, B., Petitjean, P., & Pichon, C. 2002, \aap, 396, L11
  14. Cantalupo, S. 2010, \mnras, 403, L16
  15. Carswell, B., Schaye, J., & Kim, T.-S. 2002, \apj, 578, 43
  16. Chen, H.-W., & Mulchaey, J. S. 2009, \apj, 701, 1219
  17. Conroy, C., Shapley, A. E., Tinker, J. L., Santos, M. R., & Lemson, G. 2008, \apj, 679, 1192
  18. Cooke, J., Ryan-Weber, E. V., Garel, T., & Díaz, C. G. 2014, \mnras, 441, 837
  19. Cowie, L. L., & Songaila, A. 1998, \nat, 394, 44
  20. Danforth, C. W., & Shull, J. M. 2008, \apj, 679, 194
  21. Dekel, A., Birnboim, Y., Engel, G., et al. 2009, \nat, 457, 451
  22. Ellison, S. L., Songaila, A., Schaye, J., & Pettini, M. 2000, \aj, 120, 1175
  23. Erb, D. K., Steidel, C. C., Shapley, A. E., et al. 2006a, \apj, 647, 128
  24. —. 2006b, \apj, 646, 107
  25. Faucher-Giguère, C.-A., Lidz, A., Hernquist, L., & Zaldarriaga, M. 2008, \apj, 688, 85
  26. Ferland, G. J., Porter, R. L., van Hoof, P. A. M., et al. 2013, \rmxaa, 49, 137
  27. Ford, A. B., Oppenheimer, B. D., Davé, R., et al. 2013, \mnras, 432, 89
  28. Fox, A. J., Bergeron, J., & Petitjean, P. 2008, \mnras, 388, 1557
  29. Haardt, F., & Madau, P. 2001, in Clusters of Galaxies and the High Redshift Universe Observed in X-rays, ed. D. M. Neumann & J. T. V. Tran
  30. Haardt, F., & Madau, P. 2012, \apj, 746, 125
  31. Holweger, H. 2001, in American Institute of Physics Conference Series, Vol. 598, Joint SOHO/ACE workshop ”Solar and Galactic Composition”, ed. R. F. Wimmer-Schweingruber, 23–30
  32. Iwata, I., Inoue, A. K., Matsuda, Y., et al. 2009, \apj, 692, 1287
  33. Kereš, D., Katz, N., Weinberg, D. H., & Davé, R. 2005, \mnras, 363, 2
  34. Lidz, A., Faucher-Giguère, C.-A., Dall’Aglio, A., et al. 2010, \apj, 718, 199
  35. Lopez, S., Ellison, S., D’Odorico, S., & Kim, T.-S. 2007, \aap, 469, 61
  36. Mostardi, R. E., Shapley, A. E., Nestor, D. B., et al. 2013, \apj, 779, 65
  37. Nestor, D. B., Shapley, A. E., Kornei, K. A., Steidel, C. C., & Siana, B. 2013, \apj, 765, 47
  38. Oppenheimer, B. D., Davé, R., Katz, N., Kollmeier, J. A., & Weinberg, D. H. 2012, \mnras, 420, 829
  39. Oppenheimer, B. D., & Schaye, J. 2013a, \mnras, 434, 1063
  40. —. 2013b, \mnras, 434, 1043
  41. Pieri, M. M., Frank, S., Mathur, S., et al. 2010, \apj, 716, 1084
  42. Pieri, M. M., Schaye, J., & Aguirre, A. 2006, \apj, 638, 45
  43. Pieri, M. M., Mortonson, M. J., Frank, S., et al. 2014, \mnras, 441, 1718
  44. Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2013, ArXiv e-prints, arXiv:1303.5076
  45. Prochaska, J. X., Weiner, B., Chen, H.-W., Mulchaey, J., & Cooksey, K. 2011, \apj, 740, 91
  46. Rahmati, A., & Schaye, J. 2014, \mnras, 438, 529
  47. Rahmati, A., Schaye, J., Pawlik, A. H., & Raicevic, M. 2013, \mnras, 431, 2261
  48. Rakic, O., Schaye, J., Steidel, C. C., et al. 2013, \mnras, 433, 3103
  49. Rakic, O., Schaye, J., Steidel, C. C., & Rudie, G. C. 2011, \mnras, 414, 3265
  50. —. 2012, \apj, 751, 94
  51. Reddy, N. A., & Steidel, C. C. 2009, \apj, 692, 778
  52. Reddy, N. A., Steidel, C. C., Pettini, M., et al. 2008, \apjs, 175, 48
  53. Rudie, G. C., Steidel, C. C., & Pettini, M. 2012a, \apjl, 757, L30
  54. Rudie, G. C., Steidel, C. C., Trainor, R. F., et al. 2012b, \apj, 750, 67
  55. Savage, B. D., Kim, T.-S., Wakker, B. P., et al. 2014, \apjs, 212, 8
  56. Schaye, J. 2001, \apj, 559, 507
  57. —. 2006, \apj, 643, 59
  58. Schaye, J., Aguirre, A., Kim, T.-S., et al. 2003, \apj, 596, 768
  59. Schaye, J., Carswell, R. F., & Kim, T.-S. 2007, \mnras, 379, 1169
  60. Schaye, J., Rauch, M., Sargent, W. L. W., & Kim, T.-S. 2000a, \apjl, 541, L1
  61. Schaye, J., Theuns, T., Rauch, M., Efstathiou, G., & Sargent, W. L. W. 2000b, \mnras, 318, 817
  62. Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, \apj, 588, 65
  63. Shapley, A. E., Steidel, C. C., Pettini, M., Adelberger, K. L., & Erb, D. K. 2006, \apj, 651, 688
  64. Shen, S., Madau, P., Guedes, J., et al. 2013, \apj, 765, 89
  65. Simcoe, R. A., Sargent, W. L. W., & Rauch, M. 2002, \apj, 578, 737
  66. —. 2004, \apj, 606, 92
  67. Simcoe, R. A., Sargent, W. L. W., Rauch, M., & Becker, G. 2006, \apj, 637, 648
  68. Steidel, C. C., Adelberger, K. L., Shapley, A. E., et al. 2003, \apj, 592, 728
  69. Steidel, C. C., Erb, D. K., Shapley, A. E., et al. 2010, \apj, 717, 289
  70. Steidel, C. C., Pettini, M., & Adelberger, K. L. 2001, \apj, 546, 665
  71. Steidel, C. C., Shapley, A. E., Pettini, M., et al. 2004, \apj, 604, 534
  72. Steidel, C. C., Rudie, G. C., Strom, A. L., et al. 2014, ArXiv e-prints, arXiv:1405.5473
  73. Stinson, G. S., Brook, C., Prochaska, J. X., et al. 2012, \mnras, 425, 1270
  74. Stocke, J. T., Penton, S. V., Danforth, C. W., et al. 2006, \apj, 641, 217
  75. Tepper-García, T., Richter, P., Schaye, J., et al. 2011, \mnras, 413, 190
  76. Trainor, R. F., & Steidel, C. C. 2012, \apj, 752, 39
  77. Tripp, T. M., Sembach, K. R., Bowen, D. V., et al. 2008, \apjs, 177, 39
  78. Tumlinson, J., Thom, C., Werk, J. K., et al. 2011, Science, 334, 948
  79. Turner, M. L., Schaye, J., Steidel, C. C., Rudie, G. C., & Strom, A. L. 2014, \mnras, 445, 794
  80. van de Voort, F., & Schaye, J. 2012, \mnras, 423, 2991
  81. van de Voort, F., Schaye, J., Booth, C. M., Haas, M. R., & Dalla Vecchia, C. 2011, \mnras, 414, 2458
  82. Werk, J. K., Prochaska, J. X., Thom, C., et al. 2013, \apjs, 204, 17
  83. Werk, J. K., Prochaska, J. X., Tumlinson, J., et al. 2014, \apj, 792, 8
Comments 0
Request Comment
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 minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description