We use photoionization models that are designed to reconcile the joint rest-UV-optical spectra of high-z star-forming galaxies to self-consistently infer the gas chemistry and nebular ionization and ...excitation conditions for ∼150 galaxies from the Keck Baryonic Structure Survey (KBSS), using only observations of their rest-optical nebular spectra. We find that the majority of z ∼ 2-3 KBSS galaxies are moderately O-rich, with an interquartile range in 12 + log(O/H) = 8.29-8.56, and have significantly sub-solar Fe enrichment, with an interquartile range of Fe/H = −0.79, −0.53, which contributes additional evidence in favor of super-solar O/Fe in high-z galaxies. The model-inferred ionization parameters and N/O are strongly correlated with common strong-line indices (such as O32 and N2O2), with the latter exhibiting similar behavior to local extragalactic H ii regions. In contrast, diagnostics commonly used for measuring gas-phase O/H (such as N2 and O3N2) show relatively large scatter with the overall amount of oxygen present in the gas and behave differently than observed at z ∼ 0. We provide a new calibration for using R23 to measure O/H in typical high-z galaxies, although it is most useful for relatively O-rich galaxies; combining O32 and R23 does not yield a more effective calibration. Finally, we consider the implications for the intrinsic correlations between physical conditions across the galaxy sample and find that N/O varies with O/H in high-z galaxies in a manner that is almost identical to local H ii regions. However, we do not find a strong anti-correlation between ionization parameter and metallicity (O/H or Fe/H) in high-z galaxies, which is one of the principal bases for using strong-line ratios to infer oxygen abundance.
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.
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 a detailed analysis of 317 2.0 ≤ z ≤ 2.7 star-forming galaxies from the Keck Baryonic Structure Survey. Using complementary spectroscopic observations with Keck/LRIS and Keck/MOSFIRE, as ...well as spectral energy distribution (SED) fits to broadband photometry, we examine the joint rest-UV and rest-optical properties of the same galaxies, including stellar and nebular dust attenuation, metallicity, and star formation rate (SFR). The inferred parameters of the stellar population (reddening, age, SFR, and stellar mass) are strongly dependent on the details of the assumed stellar population model and the shape of the attenuation curve. Nebular reddening is generally larger than continuum reddening, but with large scatter. Compared to local galaxies, high-redshift galaxies have lower gas-phase metallicities (and/or higher nebular excitation) at fixed nebular reddening, and higher nebular reddening at fixed stellar mass, consistent with gas fractions that increase with redshift. We find that continuum reddening is correlated with 12 + log(O/H)O3N2 at 3.0 significance, whereas nebular reddening is correlated with only 1.1 significance. This may reflect the dependence of both continuum reddening and O3N2 on the shape of the ionizing radiation field produced by the massive stars. Finally, we show that H -based and SED-based estimates of SFR exhibit significant scatter relative to one another, and on average agree only for particular combinations of spectral synthesis models and attenuation curves. We find that the SMC extinction curve predicts consistent SFRs if we assume the subsolar (0.14 Z ) binary star models that are favored for high-redshift galaxies.
Abstract
We present a detailed study of the rest-optical (3600–7000 Å) nebular spectra of ∼380 star-forming galaxies at
, obtained with Keck/Multi-object Spectrometer for Infrared Exploration ...(MOSFIRE) as part of the Keck Baryonic Structure Survey (KBSS). The KBSS-MOSFIRE sample is representative of star-forming galaxies at these redshifts, with stellar masses
and star formation rates SFR = 3–1000
yr
−1
. We focus on robust measurements of many strong diagnostic emission lines for individual galaxies: O
ii
λλ
3727, 3729, Ne
iii
λ
3869, H
β
, O
iii
4960, 5008, N
ii
λλ
6549, 6585, H
α
, and S
ii
λλ
6718, 6732. Comparisons with observations of typical local galaxies from the Sloan Digital Sky Survey and between subsamples of KBSS-MOSFIRE show that high-redshift galaxies exhibit a number of significant differences in addition to the well-known offset in log(O
iii
/H
β
) and log(N
ii
/H
α
). We argue that the primary difference between H
ii
regions in
galaxies and those at
is an enhancement in the degree of nebular excitation, as measured by O
iii
/H
β
and
(O
iii
+O
ii
)/H
β
. At the same time, KBSS-MOSFIRE galaxies are ∼10 times more massive than
galaxies with similar ionizing spectra and have higher N/O (likely accompanied by higher O/H) at fixed excitation. These results indicate the presence of harder ionizing radiation fields at fixed N/O and O/H relative to typical
galaxies, consistent with Fe-poor stellar population models that include massive binaries, and highlight a population of massive, high-specific star formation rate galaxies at high redshift with systematically different star formation histories than galaxies of similar stellar mass today.
We present a spectroscopic survey of 318 faint (R ~ 27, L ~ 0.1Llow *), Ly alpha -emission-selected galaxies (LAEs) in regions centered on the positions of hyperluminous QSOs (HLQSOs) at 2.5 < z < 3. ...A sample of 32 LAEs with rest-frame optical emission line spectra from Keck/Multi-Object Spectrometer For InfraRed Exploration (MOSFIRE) are used to interpret the LAE spectra in the context of their systemic redshifts. The fields are part of the Keck Baryonic Structure Survey, which includes substantial ancillary multi-wavelength imaging from both the ground and space. From a quantitative analysis of the diverse Ly alpha spectral morphologies, including line widths, asymmetries, and multi-peaked profiles, we find that peak widths and separations are typically smaller than among samples of more luminous continuum-selected galaxies (Lyman-break galaxies and their analogs; LBGs) at similar redshifts. We find tentative evidence for an association between Ly alpha spectral morphology and external illumination by the nearby HLQSO. Using the MOSFIRE subsample, we find that the peak of the resolved (R approx = 1300) Ly alpha line is shifted by +200 km s super(-1) with respect to systemic across a diverse set of galaxies including both LAEs and LBGs. We also find a small number of objects with significantly blueshifted Ly alpha emission, a potential indicator of accreting gas. The Ly alpha -to-H alpha line ratios measured for the MOSFIRE subset suggest that the LAEs in this sample have Ly alpha escape fractions functionof sub(esc,Ly alpha ) approx = 30%, significantly higher than typical LBG samples. Using redshifts calibrated by our MOSFIRE sample, we construct composite LAE spectra, finding the first evidence for metal-enriched outflows in such intrinsically faint high-redshift galaxies. These outflows have smaller continuum covering fractions (functionof sub(c) approx = 0.3) and velocities (v sub(ave) approx = 100-200 km s super(-1), v sub(max) approx = 500 km s super(-1)) than those associated with typical LBGs, suggesting that the gas covering fraction is a likely driver of the high Ly alpha and Ly-continuum escape fractions of LAEs with respect to LBGs. Our results suggest a similar scaling of outflow velocity with star formation rate (SFR) as is observed at lower redshifts (v sub(outflow) ~ SFR super(0.25)) and indicate that a substantial fraction of gas is ejected with v > v sub(esc). Further observations, including deep spectroscopy in the observed near-IR, will further probe the evolution and enrichment of these galaxies in the context of their gaseous environments.
ABSTRACT
We present the rest-frame optical spectroscopic properties of 60 faint (
R
AB
∼ 27;
L
∼ 0.1
L
*
) Ly
α
-selected galaxies (LAEs) at
z
≈ 2.56. These LAEs also have rest-UV spectra of their ...Ly
α
emission line morphologies, which trace the effects of interstellar and circumgalactic gas on the escape of Ly
α
photons. We find that the LAEs have diverse rest-optical spectra, but their average spectroscopic properties are broadly consistent with the extreme low-metallicity end of the populations of continuum-selected galaxies selected at
z
≈ 2–3. In particular, the LAEs have extremely high O
iii
λ
5008/H
β
ratios (log(O
iii
/H
β
) ∼ 0.8) and low N
ii
λ
6585/H
α
ratios (log(N
ii
/H
α
) < 1.15). Coupled with a detection of the O
iii
λ
4364 auroral line, these measurements indicate that the star-forming regions in faint LAEs are characterized by high electron temperatures (
T
e
≈ 1.8 × 10
4
K), low oxygen abundances (12 + log(O/H) ≈ 8.04,
Z
neb
≈ 0.22
Z
⊙
), and high excitations with respect to their more luminous continuum-selected analogs. Several of our faintest LAEs have line ratios consistent with even lower metallicities, including six with 12 + log(O/H) ≈ 6.9–7.4 (
Z
neb
≈ 0.02–0.05
Z
⊙
). We interpret these observations in light of new models of stellar evolution (including binary interactions) that have been shown to produce long-lived populations of hot, massive stars at low metallicities. We find that strong, hard ionizing continua are required to reproduce our observed line ratios, suggesting that faint galaxies are efficient producers of ionizing photons and important analogs of reionization-era galaxies. Furthermore, we investigate the physical trends accompanying Ly
α
emission across the largest current sample of combined Ly
α
and rest-optical galaxy spectroscopy, including both the 60 KBSS-Ly
α
LAEs and 368 more luminous galaxies at similar redshifts. We find that the net Ly
α
emissivity (parameterized by the Ly
α
equivalent width) is strongly correlated with nebular excitation and ionization properties and weakly correlated with dust attenuation, suggesting that metallicity plays a strong role in determining the observed properties of these galaxies by modulating their stellar spectra, nebular excitation, and dust content.
Abstract
The chemical abundance patterns of gas and stars in galaxies are powerful probes of galaxies’ star formation histories and the astrophysics of galaxy assembly but are challenging to measure ...with confidence in distant galaxies. In this paper, we report the first measurements of the correlation between stellar mass (
M
*
) and multiple tracers of chemical enrichment (including O, N, and Fe) in individual
z
∼ 2–3 galaxies, using a sample of 195 star-forming galaxies from the Keck Baryonic Structure Survey. The galaxies’ chemical abundances are inferred using photoionization models capable of reconciling high-redshift galaxies’ observed extreme rest-UV and rest-optical spectroscopic properties. We find that the
M
*
–O/H relation for our sample is relatively shallow, with moderately large scatter, and is offset ∼0.35 dex higher than the corresponding
M
*
–Fe/H relation. The two relations have very similar slopes, indicating a high level of
α
-enhancement—O/Fe ≈ 2.2 × (O/Fe)
⊙
—across two decades in
M
*
. The
M
*
–N/H relation has the steepest slope and largest intrinsic scatter, which likely results from the fact that many
z
∼ 2 galaxies are observed near or past the transition from “primary” to “secondary” N production, and may reflect uncertainties in the astrophysical origin of N. Together, these results suggest that
z
∼ 2 galaxies are old enough to have seen substantial enrichment from intermediate-mass stars, but are still young enough that Type Ia supernovae have not had time to contribute significantly to their enrichment.
Abstract
The resonantly scattered Ly
α
line illuminates the extended halos of neutral hydrogen in the circumgalactic medium of galaxies. We present integral field Keck Cosmic Web Imager observations ...of double-peaked, spatially extended Ly
α
emission in 12 relatively low-mass (
M
⋆
∼ 10
9
M
⊙
)
z
∼ 2 galaxies characterized by extreme nebular emission lines. Using individual spaxels and small bins as well as radially binned profiles of larger regions, we find that for most objects in the sample the Ly
α
blue-to-red peak ratio increases, the peak separation decreases, and the fraction of flux emerging at line center increases with radius. We use new radiative transfer simulations to model each galaxy with a clumpy, multiphase outflow with radially varying outflow velocity, and self-consistently apply the same velocity model to the low-ionization interstellar absorption lines. These models reproduce the trends of peak ratio, peak separation, and trough depth with radius, and broadly reconcile outflow velocities inferred from Ly
α
and absorption lines. The galaxies in our sample are well-described by a model in which neutral, outflowing clumps are embedded in a hotter, more highly ionized inter-clump medium (ICM), whose residual neutral content produces absorption at the systemic redshift. The peak ratio, peak separation, and trough flux fraction are primarily governed by the line-of-sight component of the outflow velocity, the H
i
column density, and the residual neutral density in the ICM respectively. The azimuthal asymmetries in the line profile further suggest nonradial gas motions at large radii and variations in the H
i
column density in the outer halos.
We present an analysis of the galaxy distribution surrounding 15 of the most luminous (> ~ 10 super(14) L sub(middot in circle); M sub(1450) Asymptotically = to -30) QSOs in the sky with z ...Asymptotically = to 2.7. Our data are drawn from the Keck Baryonic Structure Survey, which has been optimized to examine the small-scale interplay between galaxies and the intergalactic medium during the peak of the galaxy formation era at z ~ 2-3. In this work, we use the positions and spectroscopic redshifts of 1558 galaxies that lie within ~3' (4.2 h super(-1) comoving Mpc; cMpc) of the hyperluminous QSO (HLQSO) sight line in 1 of 15 independent survey fields, together with new measurements of the HLQSO systemic redshifts. By combining the spatial and redshift distributions, we measure the galaxy-HLQSO cross-correlation function, the galaxy-galaxy autocorrelation function, and the characteristic scale of galaxy overdensities surrounding the sites of exceedingly rare, extremely rapid, black hole accretion. On average, the HLQSOs lie within significant galaxy overdensities, characterized by a velocity dispersion sigmav Asymptotically = to 200 km s super(-1) and a transverse angular scale of ~25" (~200 physical kpc). We argue that such scales are expected for small groups with log (M sub(h)/M sub(middot in circle)) Asymptotically = to 13. The galaxy-HLQSO cross-correlation function has a best-fit correlation length r super(GQ) sub(0) = (7.3 + or - 1.3) h super(-1) cMpc, while the galaxy autocorrelation measured from the spectroscopic galaxy sample in the same fields has r super(GG) sub(0) = (6.0 + or - 0.5) h super(-1) cMpc. Based on a comparison with simulations evaluated at z ~ 2.6, these values imply that a typical galaxy lives in a host halo with log (M sub(h)/M sub(middot in circle)) = 11.9 + or - 0.1, while HLQSOs inhabit host halos of log (M sub(h)/M sub(middot in circle)) = 12.3 + or - 0.5. In spite of the extremely large black hole masses implied by their observed luminosities log (M sub(BH)/M sub(middot in circle)) > ~9.7, it appears that HLQSOs do not require environments very different from their much less luminous QSO counterparts. Evidently, the exceedingly low space density of HLQSOs (<, ~ 10 super(-9) cMpc super(-3)) results from a one-in-a-million event on scales < <1 Mpc, and not from being hosted by rare dark matter halos.
PDF
