Abstract
We present results of a deep spectroscopic survey quantifying the statistics of the escape of ionizing radiation from star-forming galaxies at
z
∼ 3. We measure the ratio of ionizing to ...non-ionizing UV flux density
, where
f
900
is the mean flux density evaluated over the range 880, 910 Å. We quantify the emergent ratio of ionizing to non-ionizing UV flux density by analyzing high signal-to-noise ratio composite spectra formed from subsamples with common observed properties and numbers sufficient to reduce the statistical uncertainty in the modeled IGM+CGM correction to obtain precise values of
, including a full-sample average
= 0.057 ± 0.006. We show that
increases monotonically with
, inducing an inverse correlation with UV luminosity as a by-product. We fit the composite spectra using stellar spectral synthesis together with models of the ISM in which a fraction
f
c
of the stellar continuum is covered by gas with column density
. We show that the composite spectra simultaneously constrain the intrinsic properties of the stars (
L
900
/
L
1500
)
int
along with
f
c
,
,
, and
f
esc,abs
, the absolute escape fraction of ionizing photons. We find a sample-averaged
f
esc,abs
= 0.09 ± 0.01, with subsamples falling along a linear relation
. Using the far-UV luminosity function, the distribution function
n
(
W
(Ly
α
)), and the relationship between
and
, we estimate the total ionizing emissivity of
z
∼ 3 star-forming galaxies with
M
uv
≤ −19.5, which exceeds the contribution of quasi-stellar objects by a factor of ∼3, and accounts for ∼50% of the total
ϵ
LyC
at
z
∼ 3 estimated using indirect methods.
ABSTRACT
Observations of reionization-era analogues at z ∼ 3 are a powerful tool for constraining reionization. Rest-ultraviolet observations are particularly useful, in which both direct and ...indirect tracers of ionizing-photon production and escape can be observed. We analyse a sample of 124 z ∼ 3 galaxies from the Keck Lyman Continuum Spectroscopic Survey, with sensitive spectroscopic measurements of the Lyman continuum region. We present a method of removing foreground contamination from our sample using high-resolution, multiband Hubble Space Telescope imaging. We re-measure the global properties of the cleaned sample of 13 individually detected Lyman continuum sources and 107 individually undetected sources, including a sample-averaged absolute escape fraction of fesc, abs = 0.06 ± 0.01 and a sample-averaged ratio of ionizing to non-ionizing ultraviolet flux density of 〈f900/f1500〉out = 0.040 ± 0.006, corrected for attenuation from the intergalactic and circumgalactic media. Based on composite spectra, we also recover a strong positive correlation between 〈f900/f1500〉out and Lyα equivalent width (Wλ(Ly$\rm \alpha$)) and a negative correlation between 〈f900/f1500〉out and UV luminosity. As in previous work, we interpret the relationship between 〈f900/f1500〉out and Wλ(Ly$\rm \alpha$) in terms of the modulation of the escape of ionizing radiation from star-forming galaxies based on the covering fraction of neutral gas. We also use a Wλ(Ly$\rm \alpha$)-weighted 〈f900/f1500〉out to estimate an ionizing emissivity from star-forming galaxies at z ∼ 3 as ϵLyC ≃ 5.5 × 1024 erg s−1 Hz−1 Mpc−3. This estimate, evaluated using the uncontaminated sample of this work, affirms that the contribution of galaxies to the ionizing background at z ∼ 3 is comparable to that of active galactic nuclei.
We present a combined analysis of rest-frame far-UV (FUV; 1000-2000 A) and rest-frame optical (3600-7000 A) composite spectra formed from very deep Keck/LRIS and Keck/MOSFIRE observations of a sample ...of 30 star-forming galaxies with z= 2.40 + or - 0.11, selected to be broadly representative of the full KBSS-MOSFIRE spectroscopic survey. Since the same massive stars are responsible for the observed FUV continuum and for the excitation of the observed nebular emission, a self-consistent stellar population synthesis model should simultaneously match the details of the FUV stellar+nebular continuum and-when inserted as the excitation source in photoionization models-predict all observed nebular emission line ratios. We find that only models including massive star binaries, having low stellar metallicity ( Zlow */Z sub(middot in circle)Asymptotically = to 0.1) but relatively high nebular (ionized gas-phase) abundances ( Z sub(nch)/Zmid dot in circle sub(middot in circle)Asymptotically = to 0.5), can successfully match all of the observational constraints. We show that this apparent discrepancy is naturally explained by highly super-solar O/Fe (Asymptotically = to4-5 (O/Fe)Asymptotically = to), expected for a gas whose enrichment is dominated by the products of core-collapse supernovae. While O dominates the physics of the ionized gas (and thus the nebular emission lines), Fe dominates the extreme-UV (EUV) and FUV opacity and controls the mass-loss rate from massive stars, resulting in particularly dramatic effects for massive stars in binary systems. This high nebular excitation-caused by the hard EUV spectra of Fe-poor massive stars-is much more common at high redshift (z> ~ 2) than low redshift due to systematic differences in the star formation history of typical galaxies.
We present the first spectroscopic measurements of the shape of the far-ultraviolet (far-UV; lambda = 950-1500 A) dust attenuation curve at high redshift (z~ 3). Our analysis employs rest-frame UV ...spectra of 933 galaxies at z~ 3, 121 of which have very deep spectroscopic observations (> ~7 hr) at lambda = 850-1300 A, with the Low Resolution Imaging Spectrograph on the Keck Telescope. By using an iterative approach in which we calculate the ratios of composite spectra in different bins of continuum color excess, E(B- V), we derive a dust curve that implies a lower attenuation in the far-UV for a given E(B- V) than those obtained with standard attenuation curves. We demonstrate that the UV composite spectra of z~ 3 galaxies can be modeled well by assuming our new attenuation curve, a high covering fraction of H I, and absorption from the Lyman-Werner bands of H sub(2) with a small (< ~20%) covering fraction. The low covering fraction of H sub(2) relative to that of the H sub(1) and dust suggests that most of the dust in the ISM of typical galaxies at z~ 3 is unrelated to the catalysis of H sub(2), and is associated with other phases of the ISM (i.e., the ionized and neutral gas). The far-UV dust curve implies a factor of thickapproximate2 lower dust attenuation of Lyman continuum (ionizing) photons relative to those inferred from the most commonly assumed attenuation curves for L* galaxies at z~ 3. Our results may be utilized to assess the degree to which ionizing photons are attenuated in H II regions or, more generally, in the ionized or low column density (N(H I) < ~ 10 super(17,2) cm super(-2)) neutral ISM of high-redshift galaxies.
Abstract
We investigate the evolution of galaxy gas-phase metallicity (O/H) over the range
z
= 0–3.3 using samples of ∼300 galaxies at
z
∼ 2.3 and ∼150 galaxies at
z
∼ 3.3 from the MOSDEF survey. ...This analysis crucially utilizes different metallicity calibrations at
z
∼ 0 and
z
> 1 to account for evolving interstellar medium (ISM) conditions. We find significant correlations between O/H and stellar mass (
M
*
) at
z
∼ 2.3 and
z
∼ 3.3. The low-mass power-law slope of the mass–metallicity relation (MZR) is remarkably invariant over
z
= 0–3.3, such that O/H ∝
at all redshifts in this range. At fixed
M
*
, O/H decreases with increasing redshift as
d
log(O/H)/
dz
= −0.11 ± 0.02. We find no evidence that the fundamental metallicity relation between
M
*
, O/H, and star formation rate evolves out to
z
∼ 3.3. We employ analytic chemical evolution models to place constraints on the mass and metal loading factors of galactic outflows. The efficiency of metal removal increases toward lower
M
*
at fixed redshift and toward higher redshift at fixed
M
*
. These models suggest that the slope of the MZR is primarily set by the scaling of the outflow metal loading factor with
M
*
, not by the change in gas fraction as a function of
M
*
. The evolution toward lower O/H at fixed
M
*
with increasing redshift is driven by both higher gas fraction (leading to stronger dilution of ISM metals) and higher metal removal efficiency. These results suggest that the processes governing the smooth baryonic growth of galaxies via gas flows and star formation hold in the same form over at least the past 12 Gyr.
Using a large sample of spectroscopically confirmed z~ 3 galaxies, we establish an empirical relationship between reddening (E(B- V)), neutral gas covering fraction (f sub(cov)(H I)), and the escape ...of ionizing (Lyman continuum, LyC) photons. Our sample includes 933 galaxies at z~ 3, 121 of which have deep spectroscopic observations (> ~7 hr) at 850 < ~ lambda sub(rest)< ~ 1300 A with the Low Resolution Imaging Spectrograph on Keck. The high covering fraction of outflowing optically thick H i indicated by the composite spectra of these galaxies implies that photoelectric absorption, rather than dust attenuation, dominates the depletion of LyC photons. By modeling the composite spectra as the combination of an unattenuated stellar spectrum including nebular continuum emission with one that is absorbed by H 1 and reddened by a line-of-sight extinction, we derive an empirical relationship between E(B- V) and f sub(cov)(H 1). Galaxies with redder UV continua have larger covering fractions of H 1 characterized by higher line-of-sight extinctions. We develop a model which connects the ionizing escape fraction with E(B- V), and which may be used to estimate the ionizing escape fraction for an ensemble of galaxies. Alternatively, direct measurements of the escape fraction for our sample allow us to constrain the intrinsic LyC-to-UV flux density ratio to be left angle bracket S(900 A)/S(1500 A)right angle bracket sub(int)> ~ 0.20, a value that favors stellar population models that include weaker stellar winds, a flatter initial mass function, and/or binary evolution. Last, we demonstrate how the framework discussed here may be used to assess the pathways by which ionizing radiation escapes from high-redshift galaxies.
Using observations from the MOSFIRE Deep Evolution Field survey, we investigate the physical conditions of star-forming regions in z ~ 2.3 galaxies, specifically the electron density and ionization ...state. From measurements of the OIIlambdalambda3726,3729 and SIIlambdalambda6716,6731 doublets, we find a median electron density of ~250 cm super(-3) at z ~ 2.3, an increase of an order of magnitude compared to measurements of galaxies at z ~ 0. While z ~ 2.3 galaxies are offset toward significantly higher O sub(32) values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O sub(32) versus R sub(23) and O3N2 diagrams. Based on these results, we propose that z ~ 2.3 star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and z ~ 2.3 galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from z ~ 0 to z ~ 2.3. We find that z ~ 2.3 galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving NIIlambda6584. We conclude that the offset of z ~ 2.3 galaxies from the local star-forming sequence in the NII BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at z ~ 0 and z ~ 2.
We use a newly assembled sample of 3545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z = 1.5-2.5 to constrain the relationship between UV slope (β) and dust ...attenuation (LIR/LUV IRX). Our sample significantly extends the range of LUV and β probed in previous UV-selected samples, including those as faint as M1600 = −17.4 ( ) and −2.6 β 0.0. IRX is measured using stacks of deep Herschel data, and the results are compared with predictions of the IRX−β relation for different assumptions of the stellar population model and obscuration curve. We find that z = 1.5-2.5 galaxies have an IRX−β relation that is consistent with the predictions for an SMC curve if we invoke subsolar-metallicity models currently favored for high-redshift galaxies, while the commonly assumed starburst curve overpredicts the IRX at a given β by a factor of 3. IRX is roughly constant with LUV for LUV 3 × 109 L . Thus, the commonly observed trend of fainter galaxies having bluer β may simply reflect bluer intrinsic slopes for such galaxies, rather than lower obscurations. The IRX−β relation for young/low-mass galaxies at z 2 implies a dust curve that is steeper than the SMC. The lower attenuations and higher ionizing photon output for low-metallicity stellar populations point to Lyman continuum production efficiencies, ion, that may be elevated by a factor of 2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z ∼ 2.
We present observations of Q1549-C25, an ~L* star-forming galaxy at z= 3.15 for which Lyman-continuum (LyC) radiation is significantly detected in deep Keck/LRIS spectroscopy. We find no evidence of ...contamination from a lower-redshift interloper close to the line of sight in the high signal-to-noise spectrum of Q1549-C25. Furthermore, the morphology of Q1549-C25 in V sub(606), J sub(125), and H sub(160)Hubble Space Telescope(HST) imaging reveals that the object consists of a single, isolated component within 1''. In combination, these data indicate Q1549-C25 as a clean spectroscopic detection of LyC radiation, only the second such object discovered to date at z~ 3. We model the spectral energy distribution of Q1549-C25, finding evidence of negligible dust extinction, an age (assuming continuous star formation) of ~1 Gyr, and a stellar mass of Mlow * = 7.9 x 10 super(9)Mmiddot in circle. Although it is not possible to derive strong constraints on the absolute escape fraction of LyC emission, f sub(esc)(LyC), from a single object, we use simulations of intergalactic and circumgalactic absorption to infer f sub(esc)(LyC) > or =, slanted 0.51 at 95% confidence. The combination of deep Keck/LRIS spectroscopy and Hubble Space Telescope imaging is required to assemble a larger sample of objects like Q1549-C25, and obtain robust constraints on the average f sub(esc)(LyC) at z~ 3 and beyond.
We use the deep ground-based optical photometry of the Lyman Break Galaxy (LBG) Survey to derive robust measurements of the faint-end slope ( alpha ) of the UV luminosity function (LF) at redshifts ...1.9 <= z <= 3.4. Our sample includes >2000 spectroscopic redshifts and 31000 LBGs in 31 spatially independent fields over a total area of 3261 arcmin2. These data allow us to select galaxies to 0.07L* and 0.10L* at z ~ 2 and z ~ 3, respectively. A maximum-likelihood analysis indicates steep values of alpha (z = 2) = -1.73 ± 0.07 and alpha (z = 3) = -1.73 ± 0.13. This result is robust to luminosity-dependent systematics in the Ly alpha equivalent width and reddening distributions, and is similar to the steep values advocated at z 4, and implies that 93% of the unobscured UV luminosity density at z ~ 2-3 arises from sub-L* galaxies. With a realistic luminosity-dependent reddening distribution, faint to moderately luminous galaxies account for 70% and 25% of the bolometric luminosity density and present-day stellar mass density, respectively, when integrated over 1.9 <= z < 3.4. We find a factor of 8-9 increase in the star-formation rate density between z ~ 6 and z ~ 2, due to both a brightening of L* and an increasing dust correction proceeding to lower redshifts. Combining the UV LF with stellar mass estimates suggests a relatively steep low-mass slope of the stellar mass function at high redshift. The previously observed discrepancy between the integral of the star-formation history and stellar mass density measurements at z ~ 2 may be reconciled by invoking a luminosity-dependent reddening correction to the star-formation history combined with an accounting for the stellar mass contributed by UV-faint galaxies. The steep and relatively constant faint-end slope of the UV LF at z 2 contrasts with the shallower slope inferred locally, suggesting that the evolution in the faint-end slope may be dictated simply by the availability of low-mass halos capable of supporting star formation at z 2.