Galaxy metallicity scaling relations provide a powerful tool for understanding galaxy evolution, but obtaining unbiased global galaxy gas-phase oxygen abundances requires proper treatment of the ...various line-emitting sources within spectroscopic apertures. We present a model framework that treats galaxies as ensembles of H ii and diffuse ionized gas (DIG) regions of varying metallicities. These models are based upon empirical relations between line ratios and electron temperature for H ii regions, and DIG strong-line ratio relations from SDSS-IV MaNGA IFU data. Flux-weighting effects and DIG contamination can significantly affect properties inferred from global galaxy spectra, biasing metallicity estimates by more than 0.3 dex in some cases. We use observationally motivated inputs to construct a model matched to typical local star-forming galaxies, and quantify the biases in strong-line ratios, electron temperatures, and direct-method metallicities as inferred from global galaxy spectra relative to the median values of the H ii region distributions in each galaxy. We also provide a generalized set of models that can be applied to individual galaxies or galaxy samples in atypical regions of parameter space. We use these models to correct for the effects of flux-weighting and DIG contamination in the local direct-method mass-metallicity and fundamental metallicity relations, and in the mass-metallicity relation based on strong-line metallicities. Future photoionization models of galaxy line emission need to include DIG emission and represent galaxies as ensembles of emitting regions with varying metallicity, instead of as single H ii regions with effective properties, in order to obtain unbiased estimates of key underlying physical properties.
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 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.
ABSTRACT
We present constraints on the massive star and ionized gas properties for a sample of 62 star-forming galaxies at z ∼ 2.3. Using BPASS stellar population models, we fit the rest-UV spectra ...of galaxies in our sample to estimate age and stellar metallicity which, in turn, determine the ionizing spectrum. In addition to the median properties of well-defined subsets of our sample, we derive the ages and stellar metallicities for 30 high-SNR individual galaxies – the largest sample of individual galaxies at high redshift with such measurements. Most galaxies in this high-SNR subsample have stellar metallicities of 0.001 < Z* < 0.004. We then use Cloudy + BPASS photoionization models to match observed rest-optical line ratios and infer nebular properties. Our high-SNR subsample is characterized by a median ionization parameter and oxygen abundance, respectively, of log (U)med = −2.98 ± 0.25 and 12 + log (O/H)med = 8.48 ± 0.11. Accordingly, we find that all galaxies in our sample show evidence for α-enhancement. In addition, based on inferred log (U) and 12 + log (O/H) values, we find that the local relationship between ionization parameter and metallicity applies at z ∼ 2. Finally, we find that the high-redshift galaxies most offset from the local excitation sequence in the BPT diagram are the most α-enhanced. This trend suggests that α-enhancement resulting in a harder ionizing spectrum at fixed oxygen abundance is a significant driver of the high-redshift galaxy offset on the BPT diagram relative to local systems. The ubiquity of α-enhancement among z ∼ 2.3 star-forming galaxies indicates important differences between high-redshift and local galaxies that must be accounted for in order to derive physical properties at high redshift.
Abstract
We investigate the nature of the relation among stellar mass, star formation rate, and gas-phase metallicity (the
–SFR–
Z
relation) at high redshifts using a sample of 260 star-forming ...galaxies at
z
∼ 2.3 from the MOSDEF survey. We present an analysis of the high-redshift
–SFR–
Z
relation based on several emission-line ratios for the first time. We show that a
–SFR–
Z
relation clearly exists at
z
∼ 2.3. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of
z
∼ 0 and
z
∼ 2.3 galaxies, we find that
z
∼ 2.3 galaxies have ∼0.1 dex lower metallicity at fixed
and SFR. In the context of chemical evolution models, this evolution of the
–SFR–
Z
relation suggests an increase with redshift of the mass-loading factor at fixed
, as well as a decrease in the metallicity of infalling gas that is likely due to a lower importance of gas recycling relative to accretion from the intergalactic medium at high redshifts. Performing this analysis simultaneously with multiple metallicity-sensitive line ratios allows us to rule out the evolution in physical conditions (e.g., N/O ratio, ionization parameter, and hardness of the ionizing spectrum) at fixed metallicity as the source of the observed trends with redshift and with SFR at fixed
at
z
∼ 2.3. While this study highlights the promise of performing high-order tests of chemical evolution models at high redshifts, detailed quantitative comparisons ultimately await a full understanding of the evolution of metallicity calibrations with redshift.
We use extensive spectroscopy from the MOSFIRE Deep Evolution Field survey to investigate the relationships between rest-frame optical emission line equivalent widths (W) and a number of galaxy and ...interstellar medium (ISM) characteristics for a sample of 1134 star-forming galaxies at redshifts 1.4 z 3.8. We examine how the equivalent widths of , , λλ4960, 5008, + Hβ, , and , depend on stellar mass, UV slope, age, star formation rate (SFR) and specific SFR (sSFR), ionization parameter and excitation conditions (O32 and /Hβ), gas-phase metallicity, and ionizing photon production efficiency ( ion). The trend of increasing W with decreasing stellar mass is strongest for (and +Hβ). More generally, the equivalent widths of all the lines increase with redshift at a fixed stellar mass or fixed gas-phase metallicity, suggesting that high equivalent width galaxies are common at high redshift. This redshift evolution in equivalent widths can be explained by the increase in SFR and decrease in metallicity with redshift at a fixed stellar mass. Consequently, the dependence of W on sSFR is largely invariant with redshift, particularly when examined for galaxies of a given metallicity. Our results show that high equivalent width galaxies, specifically those with high , have low stellar masses, blue UV slopes, young ages, high sSFRs, ISM line ratios indicative of high ionization parameters, high ion, and low metallicities. As these characteristics are often attributed to galaxies with high ionizing escape fractions, galaxies with high W are likely candidates for the population that dominates cosmic reionization.
Abstract
We examine the factors responsible for the variation in the ionization parameter (
U
) of high-redshift star-forming galaxies based on medium-resolution JWST/NIRSpec observations obtained by ...the Cosmic Evolution Early Release Science survey. The sample consists of 48 galaxies with redshifts
z
spec
= 2.7−6.3, which are largely representative of typical galaxies at these redshifts. The S
ii
λ
λ
6718, 6733 doublet is used to estimate electron densities (
n
e
), and dust-corrected H
α
luminosities are used to compute ionizing photon rates (
Q
). Using composite spectra of galaxies in bins of O
iii
λ
λ
4960, 5008/O
ii
λ
λ
3727, 3730 (O32) as a proxy for
U
, we determine that galaxies with higher O32 have 〈
n
e
〉 ≃ 500 cm
−3
that are ≳5 × larger than that of lower-O32 galaxies. We do not find a significant difference in 〈
Q
〉 between low- and high-O32 galaxies. Photoionization modeling indicates a large spread in
log
U
of ≈1.5 dex at a fixed
Z
neb
. On the other hand, the data indicate a highly significant correlation between
U
and star-formation-rate surface density (Σ
SFR
), which appears to be redshift invariant at
z
∼ 1.6−6.3, and possibly up to
z
∼ 9.5. We consider several avenues through which metallicity and Σ
SFR
(or gas density) may influence
U
, including variations in
n
e
and
Q
, internal dust extinction of ionizing photons, and the effects of gas density on the volume filling fraction. Based on these considerations, we conclude that gas density may play a more central role than metallicity in modulating
U
at these redshifts.
Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of active galactic nucleus (AGN)-driven ionized outflows in a sample of 159 AGNs at 1.4 ≤ z ≤ 3.8. The sample ...spans AGN bolometric luminosities of 1044-47 erg s−1 and includes both quiescent and star-forming galaxies extending across 3 orders of magnitude in stellar mass. We identify and characterize outflows from the Hβ, O iii, H , and N ii emission line spectra. We detect outflows in 17% of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of ∼400-3500 km s and spatial extents of 0.3-11.0 kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically 0.1-1 and increase with AGN luminosity, capable of exceeding unity at . In these more luminous sources, the ionized outflow alone is likely sufficient to regulate star formation and, when combined with outflowing neutral and molecular gas, may be able to quench star formation in their host galaxies.
ABSTRACT
We present detections of O iii λ4363 and direct-method metallicities for star-forming galaxies at z = 1.7–3.6. We combine new measurements from the MOSFIRE Deep Evolution Field (MOSDEF) ...survey with literature sources to construct a sample of 18 galaxies with direct-method metallicities at z > 1, spanning 7.5 < 12+log(O/H) < 8.2 and log(M*/M⊙) = 7–10. We find that strong-line calibrations based on local analogues of high-redshift galaxies reliably reproduce the metallicity of the z > 1 sample on average. We construct the first mass–metallicity relation at z > 1 based purely on direct-method O/H, finding a slope that is consistent with strong-line results. Direct-method O/H evolves by ≲0.1 dex at fixed M* and star formation rate from z ∼ 0 to 2.2. We employ photoionization models to constrain the ionization parameter and ionizing spectrum in the high-redshift sample. Stellar models with supersolar O/Fe and binary evolution of massive stars are required to reproduce the observed strong-line ratios. We find that the z > 1 sample falls on the z ∼ 0 relation between ionization parameter and O/H, suggesting no evolution of this relation from z ∼ 0 to z ∼ 2. These results suggest that the offset of the strong-line ratios of this sample from local excitation sequences is driven primarily by a harder ionizing spectrum at fixed nebular metallicity compared to what is typical at z ∼ 0, naturally explained by supersolar O/Fe at high redshift caused by rapid formation time-scales. Given the extreme nature of our z > 1 sample, the implications for representative z ∼ 2 galaxy samples at ∼1010 M⊙ are unclear, but similarities to z > 6 galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization.
Abstract
We utilize medium-resolution JWST/NIRSpec observations of 164 galaxies at
z
= 2.0–9.3 from the Cosmic Evolution Early Release Science (CEERS) survey to investigate the evolution of the ...excitation and ionization properties of galaxies at high redshifts. Our results represent the first statistical constraints on the evolution of the O
III
/H
β
versus N
II
/H
α
, S
II
/H
α
, and O
I
/H
α
“BPT” diagrams at
z
> 2.7, and the first analysis of the O
32
versus R
23
diagram at
z
> 4 with a large sample. We divide the sample into five redshift bins containing 30–40 galaxies each. The subsamples at
z
∼ 2.3,
z
∼ 3.3, and
z
∼ 4.5 are representative of the main-sequence star-forming galaxy population at these redshifts, while the
z
∼ 5.6 and
z
∼ 7.5 samples are likely biased toward high specific star formation rate, due to selection effects. Using composite spectra, we find that each subsample at
z
= 2.0–6.5 falls on the same excitation sequence in the N
II
and S
II
BPT diagrams and the O
32
–R
23
diagram on average, and is offset from the sequences followed by
z
= 0 H
II
regions in the same diagrams. The direction of these offsets are consistent with high-redshift star-forming galaxies uniformly having harder ionizing spectra than typical local galaxies at fixed nebular metallicity. The similarity of the average line ratios suggests that the ionization conditions of the interstellar medium do not strongly evolve between
z
∼ 2 and
z
∼ 6. Overall, the rest-optical line ratios suggest the
z
= 2.7–9.3 CEERS/NIRSpec galaxies at log(
M
*
/
M
⊙
) ∼ 7.5–10 have high degrees of ionization and moderately low oxygen abundances (∼0.1–0.3
Z
⊙
), but are not extremely metal-poor (<0.1
Z
⊙
) even at
z
> 6.5.