MOSFIRE Deep Evolution Field Survey Research Articles

Ly α Halo Properties and Dust in the Circumgalactic Medium of z ∼ 2 Star-forming Galaxies

We present Keck Cosmic Web Imager integral-field unit observations around extended Lyα halos of 27 typical star-forming galaxies with redshifts 2.0 < z < 3.2 drawn from the MOSFIRE Deep Evolution Field survey. We examine the average Lyα surface brightness profiles in bins of star formation rate (SFR), stellar mass (M *), age, stellar continuum reddening, SFR surface density (ΣSFR), and ΣSFR normalized by stellar mass (ΣsSFR). The scale lengths of the halos correlate with stellar mass, age, and stellar continuum reddening and anticorrelate with SFR, ΣSFR, and ΣsSFR. These results are consistent with a scenario in which the down-the-barrel fraction of Lyα emission is modulated by the low-column-density channels in the interstellar medium, and in which the neutral gas covering fraction is related to the physical properties of the galaxies. Specifically, we find that this covering fraction increases with stellar mass, age, and E(B − V) and decreases with SFR, ΣSFR, and ΣsSFR. We also find that the resonantly scattered Lyα emission suffers greater attenuation than the (nonresonant) stellar continuum emission, and that the difference in attenuation increases with stellar mass, age, and stellar continuum reddening, and decreases with ΣsSFR. These results imply that more reddened galaxies have more dust in their circumgalactic medium.

The MOSDEF survey: properties of warm ionized outflows at z = 1.4–3.8

ABSTRACT We use the large spectroscopic data set of the MOSFIRE Deep Evolution Field survey to investigate the kinematics and energetics of ionized gas outflows. Using a sample of 598 star-forming galaxies at redshift 1.4 &amp;lt; z &amp;lt; 3.8, we decompose [O iii] and $\rm {H}\,\alpha$ emission lines into narrow and broad components, finding significant detections of broad components in 10 per cent of the sample. The ionized outflow velocity from individual galaxies appears independent of galaxy properties, such as stellar mass, star formation rate (SFR), and SFR surface density (ΣSFR). Adopting a simple outflow model, we estimate the mass-, energy-, and momentum-loading factors of the ionized outflows, finding modest values with averages of 0.33, 0.04, and 0.22, respectively. The larger momentum- than energy-loading factors, for the adopted physical parameters, imply that these ionized outflows are primarily momentum driven. We further find a marginal correlation (2.5σ) between the mass-loading factor and stellar mass in agreement with predictions by simulations, scaling as ηm$\propto M_{\star }^{-0.45}$. This shallow scaling relation is consistent with these ionized outflows being driven by a combination of mechanical energy generated by supernovae explosions and radiation pressure acting on dusty material. In a majority of galaxies, the outflowing material does not appear to have sufficient velocity to escape the gravitational potential of their host, likely recycling back at later times. Together, these results suggest that the ionized outflows traced by nebular emission lines are negligible, with the bulk of mass and energy carried out in other gaseous phases.

The Impact of Star-formation-rate Surface Density on the Electron Density and Ionization Parameter of High-redshift Galaxies* *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 use the large spectroscopic data set of the MOSFIRE Deep Evolution Field survey to investigate some of the key factors responsible for the elevated ionization parameters (U) inferred for high-redshift galaxies, focusing in particular on the role of star-formation-rate surface density (ΣSFR). Using a sample of 317 galaxies with spectroscopic redshifts z spec ≃ 1.9–3.7, we construct composite rest-frame optical spectra in bins of ΣSFR and infer electron densities, n e, using the ratio of the [O ii] λ λ3727, 3730 doublet. Our analysis suggests a significant (≃3σ) correlation between n e and ΣSFR. We further find significant correlations between U and ΣSFR for composite spectra of a subsample of 113 galaxies, and for a smaller sample of 25 individual galaxies with inferences of U. The increase in n e—and possibly also the volume filling factor of dense clumps in H ii regions—with ΣSFR appear to be important factors in explaining the relationship between U and ΣSFR. Further, the increase in n e and SFR with redshift at a fixed stellar mass can account for most of the redshift evolution of U. These results suggest that the gas density, which sets n e and the overall level of star formation activity, may play a more important role than metallicity evolution in explaining the elevated ionization parameters of high-redshift galaxies.

An Updated Dust-to-Star Geometry: Dust Attenuation Does Not Depend on Inclination in 1.3 ≤z ≤2.6 Star-forming Galaxies from MOSDEF

We investigate dust attenuation and its dependence on viewing angle for 308 star-forming galaxies at 1.3 ≤ z ≤ 2.6 from the MOSFIRE Deep Evolution Field survey. We divide galaxies with a detected Hα emission line and coverage of Hβ into eight groups by stellar mass, star formation rate (SFR), and inclination (i.e., axis ratio), and we then stack their spectra. From each stack, we measure the Balmer decrement and gas-phase metallicity, and then we compute the median A V and UV continuum spectral slope (β). First, we find that none of the dust properties (Balmer decrement, A V, or β) varies with the axis ratio. Second, both stellar and nebular attenuation increase with increasing galaxy mass, showing little residual dependence on SFR or metallicity. Third, nebular emission is more attenuated than stellar emission, and this difference grows even larger at higher galaxy masses and SFRs. Based on these results, we propose a three-component dust model in which attenuation predominantly occurs in star-forming regions and large, dusty star-forming clumps, with minimal attenuation in the diffuse ISM. In this model, nebular attenuation primarily originates in clumps, while stellar attenuation is dominated by star-forming regions. Clumps become larger and more common with increasing galaxy mass, creating the above mass trends. Finally, we argue that a fixed metal yield naturally leads to mass regulating dust attenuation. Infall of low-metallicity gas increases the SFR and lowers the metallicity, but leaves the dust column density mostly unchanged. We quantify this idea using the Kennicutt–Schmidt and fundamental metallicity relations, showing that galaxy mass is indeed the primary driver of dust attenuation.

The MOSDEF survey: probing resolved stellar populations at z ∼ 2 Using a new bayesian-defined morphology metric called patchiness

ABSTRACT We define a new morphology metric called ‘patchiness’ (P) that is sensitive to deviations from the average of a resolved distribution, does not require the galaxy centre to be defined, and can be used on the spatially resolved distribution of any galaxy property. While the patchiness metric has a broad range of applications, we demonstrate its utility by investigating the distribution of dust in the interstellar medium (ISM) of 310 star-forming galaxies at spectroscopic redshifts 1.36 &amp;lt; z &amp;lt; 2.66 observed by the MOSFIRE Deep Evolution Field survey. The stellar continuum reddening distribution, derived from high-resolution multiwaveband CANDELS/3D-HST imaging, is quantified using the patchiness, Gini, and M20 coefficients. We find that the reddening maps of high-mass galaxies, which are dustier and more metal-rich on average, tend to exhibit patchier distributions (high P) with the reddest components concentrated within a single region (low M20). Our results support a picture where dust is uniformly distributed in low-mass galaxies (≲1010 M⊙), implying efficient mixing of dust throughout the ISM. On the other hand, the dust distribution is patchier in high-mass galaxies (≳1010 M⊙). Dust is concentrated near regions of active star formation and dust mixing time-scales are expected to be longer in high-mass galaxies, such that the outskirt regions of these physically larger galaxies remain relatively unenriched. This study presents direct evidence for patchy dust distributions on scales of a few kpc in high-redshift galaxies, which previously has only been suggested as a possible explanation for the observed differences between nebular and stellar continuum reddening, star formation rate indicators, and dust attenuation curves.

Infrared Spectral Energy Distributions and Dust Masses of Sub-solar Metallicity Galaxies at z ∼ 2.3

Abstract We present results from Atacama Large Millimeter/submillimeter Array (ALMA) 1.2 mm continuum observations of a sample of 27 star-forming galaxies at z = 2.1–2.5 from the MOSFIRE Deep Evolution Field survey with metallicity and star formation rate measurements from optical emission lines. Using stacks of Spitzer, Herschel, and ALMA photometry (rest frame ∼8–400 μm), we examine the infrared (IR) spectral energy distributions (SED) of z ∼ 2.3 subsolar-metallicity (∼0.5 Z ⊙) luminous infrared galaxies (LIRGs). We find that the data agree well with an average template of higher-luminosity local low-metallicity dwarf galaxies (reduced χ 2 = 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ultraluminous IR galaxies, even in the most favorable case (with reduced χ 2 = 2.8), the templates are rejected at &gt;98% confidence. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar-metallicity galaxies may result from different grain properties or a harder/more intense ionizing radiation field that increases the dust temperature. The obscured star formation rate (SFR) indicated by the far-IR emission of the subsolar-metallicity galaxies is only ∼60% of the total SFR, considerably lower than that of the local LIRGs with ∼96%–97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data overestimate the Rayleigh–Jeans tail measurements by a factor of 2–20. These templates underestimate IR luminosities if fit to the observed ALMA fluxes by &gt;0.4 dex. At a given stellar mass or metallicity, dust masses at z ∼ 2.3 are an order of magnitude higher than z ∼ 0. Given the predicted molecular gas fractions, the observed z ∼ 2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies with similar metallicities.

The MOSFIRE Deep Evolution Field Survey: Implications of the Lack of Evolution in the Dust Attenuation–Mass Relation to z ∼ 2*

Abstract We investigate the relationship between dust attenuation and stellar mass (M *) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field survey at z ∼ 2.3 and the Sloan Digital Sky Survey (SDSS) at z ∼ 0, augmenting the latter optical data set with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of Hα and Hβ emission lines, and photometric measurements of the rest-UV stellar continuum. The Hα/Hβ ratio is used to determine the magnitude of attenuation at the wavelength of Hα, A Hα . Rest-UV colors and spectral energy distribution fitting are used to estimate A 1600, the magnitude of attenuation at a rest wavelength of 1600 Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and M * over the redshift range z ∼ 0 to z ∼ 2.3. Folding in the latest estimates of the evolution of M dust, (M dust/M gas), and gas surface density at fixed M *, we find that the expected M dust and dust mass surface density are both significantly higher at z ∼ 2.3 than at z ∼ 0. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation versus M *, it is essential to determine the relationship between metallicity and (M dust/M gas), the dust mass absorption coefficient and dust geometry, and the evolution of these relations and quantities from z ∼ 0 to z ∼ 2.3.

The MOSDEF survey: the dependence of H α-to-UV SFR ratios on SFR and size at z ∼ 2

ABSTRACT We perform an aperture-matched analysis of dust-corrected H α and UV star formation rates (SFRs) using 303 star-forming galaxies with spectroscopic redshifts 1.36 &amp;lt; zspec &amp;lt; 2.66 from the MOSFIRE Deep Evolution Field survey. By combining H α and H β emission line measurements with multiwaveband resolved Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey/3D-HST imaging, we directly compare dust-corrected H α and UV SFRs, inferred assuming a fixed attenuation curve shape and constant SFHs, within the spectroscopic aperture. Previous studies have found that H α and UV SFRs inferred with these assumptions generally agree for typical star-forming galaxies, but become increasingly discrepant for galaxies with higher SFRs (≳100 M⊙ yr−1), with H α-to-UV SFR ratios being larger for these galaxies. Our analysis shows that this trend persists even after carefully accounting for the apertures over which H α and UV-based SFRs (and the nebular and stellar continuum reddening) are derived. Furthermore, our results imply that H α SFRs may be higher in the centres of large galaxies (i.e. where there is coverage by the spectroscopic aperture) compared to their outskirts, which could be indicative of inside-out galaxy growth. Overall, we suggest that the persistent difference between nebular and stellar continuum reddening and high H α-to-UV SFR ratios at the centres of large galaxies may be indicative of a patchier distribution of dust in galaxies with high SFRs.

The MOSDEF Survey: Environmental Dependence of the Gas-phase Metallicity of Galaxies at 1.4 ≤ z ≤ 2.6*

Abstract Using the near-IR spectroscopy of the MOSFIRE Deep Evolution Field survey, we investigate the role of the local environment in the gas-phase metallicity of galaxies. The local environment measurements are derived from accurate and uniformly calculated photometric redshifts with well-calibrated probability distributions. Based on rest-frame optical emission lines, [N ii]λ6584 and Hα, we measure gas-phase oxygen abundances of 167 galaxies at 1.37 ≤ z ≤ 1.7 and 303 galaxies at 2.09 ≤ z ≤ 2.61, located in diverse environments. We find that at z ∼ 1.5, the average metallicity of galaxies in overdensities with M * ∼ 109.8 M ⊙, 1010.2 M ⊙, and 1010.8 M ⊙ is higher relative to their field counterparts by 0.094 ± 0.051, 0.068 ± 0.028, and 0.052 ± 0.043 dex, respectively. However, this metallicity enhancement does not exist at higher redshift, z ∼ 2.3, where, compared to the field galaxies, we find 0.056 ± 0.043, 0.056 ± 0.028, and 0.096 ± 0.034 dex lower metallicity for galaxies in overdense environments with M * ∼ 109.8 M ⊙, 1010.2 M ⊙ and 1010.7 M ⊙, respectively. Our results suggest that, at 1.37 ≤ z ≤ 2.61, the variation of mass–metallicity relation with local environment is small (&lt;0.1 dex), and reverses at z ∼ 2. Our results support the hypothesis that, at the early stages of cluster formation, owing to efficient gas cooling, galaxies residing in overdensities host a higher fraction of pristine gas with prominent primordial gas accretion, which lowers their gas-phase metallicity compared to their coeval field galaxies. However, as the universe evolves to lower redshifts (z ≲ 2), the shock-heated gas in overdensities cannot cool down efficiently, and galaxies become metal-rich rapidly due to the suppression of pristine gas inflow and re-accretion of metal-enriched outflows in overdensities.

The MOSDEF survey: differences in SFR and metallicity for morphologically selected mergers at z ∼ 2

ABSTRACT We study the properties of 55 morphologically-identified merging galaxy systems at z ∼ 2. These systems are flagged as mergers based on features such as tidal tails, double nuclei, and asymmetry. Our sample is drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey, along with a control sample of isolated galaxies at the same redshift. We consider the relationships between stellar mass, star formation rate (SFR), and gas-phase metallicity for both merging and non-merging systems. In the local universe, merging systems are characterized by an elevated SFR and depressed metallicity compared to isolated systems at a given mass. Our results indicate SFR enhancement and metallicity deficit for merging systems relative to non-merging systems for a fixed stellar mass at z ∼ 2, though larger samples are required to establish these preliminary results with higher statistical significance. In future work, it will be important to establish if the enhanced SFR and depressed metallicity in high-redshift mergers deviate from the ‘fundamental metallicity relation,’ as is observed in mergers in the local universe, and therefore shed light on gas flows during galaxy interactions.

The MOSDEF Survey: The First Direct Measurements of the Nebular Dust Attenuation Curve at High Redshift*

Abstract We use a sample of 532 star-forming galaxies at redshifts z ≃ 1.4–2.6 with deep rest-frame optical spectra from the MOSFIRE Deep Evolution Field (MOSDEF) survey to place the first constraints on the nebular attenuation curve at high redshift. Based on the first five low-order Balmer emission lines detected in the composite spectra of these galaxies (Hα through ), we derive a nebular attenuation curve that is similar in shape to that of the Galactic extinction curve, suggesting that the dust covering fraction and absorption/scattering properties along the lines of sight to massive stars at high redshift are similar to those of the average Milky Way sight line. The curve derived here implies nebular reddening values that are, on average, systematically larger than those derived for the stellar continuum. In the context of stellar population synthesis models that include the effects of stellar multiplicity, the difference in reddening of the nebular lines and stellar continuum may imply molecular cloud crossing timescales that are a factor of longer than those inferred for local molecular clouds, star formation rates that are constant or increasing with time such that newly formed and dustier OB associations always dominate the ionizing flux, and/or that the dust responsible for reddening the nebular emission may be associated with nonmolecular (i.e., ionized and neutral) phases of the interstellar medium. Our analysis points to a variety of investigations of the nebular attenuation curve that will be enabled with the next generation of ground- and space-based facilities.

The MOSDEF Survey: Neon as a Probe of ISM Physical Conditions at High Redshift*

Abstract We present results on the properties of neon emission in z ∼ 2 star-forming galaxies drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Doubly ionized neon ([Ne iii]λ3869) is detected at ≥3σ in 61 galaxies, representing ∼25% of the MOSDEF sample with Hα, Hβ, and [O iii]λ5007 detections at similar redshifts. We consider the neon emission-line properties of both individual galaxies with [Ne iii]λ3869 detections and composite z ∼ 2 spectra binned by stellar mass. With no requirement of [Ne iii]λ3869 detection, the latter provide a more representative picture of neon emission-line properties in the MOSDEF sample. The [Ne iii]λ3869/[O ii]λ3727 ratio (Ne3O2) is anticorrelated with stellar mass in z ∼ 2 galaxies, as expected based on the mass–metallicity relation. It is also positively correlated with the [O iii]λ5007/[O ii]λ3727 ratio (O32), but z ∼ 2 line ratios are offset toward higher Ne3O2 at fixed O32, compared with both local star-forming galaxies and individual H ii regions. Despite the offset toward higher Ne3O2 at fixed O32 at z ∼ 2, biases in inferred Ne3O2-based metallicity are small. Accordingly, Ne3O2 may serve as an important metallicity indicator deep into the reionization epoch. Analyzing additional rest-optical line ratios including [Ne iii]λ3869/[O iii]λ5007 (Ne3O3) and [O iii]λ5007/Hβ (O3Hβ), we conclude that the nebular emission-line ratios of z ∼ 2 star-forming galaxies suggest a harder ionizing spectrum (lower stellar metallicity, i.e., Fe/H) at fixed gas-phase oxygen abundance, compared to systems at z ∼ 0. These new results based on neon lend support to the physical picture painted by oxygen, nitrogen, hydrogen, and sulfur emission of an ionized interstellar medium in high-redshift star-forming galaxies irradiated by chemically young, α-enhanced massive stars.

The MOSDEF Survey: calibrating the relationship between H α star formation rate and radio continuum luminosity at 1.4 &amp;lt; z &amp;lt; 2.6

ABSTRACT The observed empirical relation between the star formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high-redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H α SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at 1.4 &amp;lt; z &amp;lt; 2.6 with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis, we compare the observed radio continuum luminosities with those predicted from the dust-corrected H α SFR assuming a range of z ∼ 0 relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H α SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to z ∼ 2.

The MOSDEF Survey: The Variation of the Dust Attenuation Curve with Metallicity

Abstract We derive a UV-optical stellar dust attenuation curve of galaxies at z = 1.4–2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field survey with Hα, Hβ, and [N ii]λ 6585 spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., Hα-to-Hβ luminosities. The average attenuation curve for the high-metallicity sample ( , corresponding to ) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175 Å extinction bump that is ∼0.5× the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample ( ) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the 3σ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B − V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star-forming regions but also reddening the continuum in the high-metallicity galaxies.

The MOSDEF Survey: Kinematic and Structural Evolution of Star-forming Galaxies at 1.4 ≤ z ≤ 3.8

Abstract We present ionized gas kinematics for 681 galaxies at from the MOSFIRE Deep Evolution Field survey, measured using models that account for random galaxy-slit misalignments together with structural parameters derived from CANDELS Hubble Space Telescope (HST) imaging. Kinematics and sizes are used to derive dynamical masses. Baryonic masses are estimated from stellar masses and inferred gas masses from dust-corrected star formation rates (SFRs) and the Kennicutt–Schmidt relation. We measure resolved rotation for 105 galaxies. For the remaining 576 galaxies we use models based on HST imaging structural parameters together with integrated velocity dispersions and baryonic masses to statistically constrain the median ratio of intrinsic ordered to disordered motion, . We find that increases with increasing stellar mass and decreasing specific SFR (sSFR). These trends may reflect marginal disk stability, where systems with higher gas fractions have thicker disks. For galaxies with detected rotation we assess trends between their kinematics and mass, sSFR, and baryon surface density ( ). Intrinsic dispersion correlates most with , and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at are baryon dominated within their effective radii ( ), with / increasing over time. The inferred baryon fractions within , , decrease over time, even at fixed mass, size, or surface density. At fixed redshift, does not appear to vary with stellar mass but increases with decreasing and increasing . For galaxies at , the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

The MOSDEF Survey: A Census of AGN-driven Ionized Outflows at z = 1.4–3.8

Abstract 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.

The MOSDEF survey: direct-method metallicities and ISM conditions at z ∼ 1.5–3.5

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 &amp;gt; 1, spanning 7.5 &amp;lt; 12+log(O/H) &amp;lt; 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 &amp;gt; 1 sample on average. We construct the first mass–metallicity relation at z &amp;gt; 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 &amp;gt; 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 &amp;gt; 1 sample, the implications for representative z ∼ 2 galaxy samples at ∼1010 M⊙ are unclear, but similarities to z &amp;gt; 6 galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization.

The MOSDEF Survey: Sulfur Emission-line Ratios Provide New Insights into Evolving Interstellar Medium Conditions at High Redshift∗

Abstract We present results on the emission-line properties of 1.3 ≤ z ≤ 2.7 galaxies drawn from the complete the MOSFIRE Deep Evolution Field (MOSDEF) survey. Specifically, we use observations of the emission-line diagnostic diagram of [O iii]λ 5007/Hβ versus [S ii]λλ6717,6731/Hα, i.e., the “[S ii] BPT diagram,” to gain insight into the physical properties of high-redshift star-forming regions. High-redshift MOSDEF galaxies are offset toward lower [S ii]λλ6717,6731/Hα at fixed [O iii]λ5007/Hβ, relative to local galaxies from the Sloan Digital Sky Survey (SDSS). Furthermore, at fixed [O iii]λ5007/Hβ, local SDSS galaxies follow a trend of decreasing [S ii]λλ6717,6731/Hα as the surface density of star formation (ΣSFR) increases. We explain this trend in terms of the decreasing fractional contribution from diffuse ionized gas (f DIG) as ΣSFR increases in galaxies, which causes galaxy-integrated line ratios to shift toward the locus of pure H ii-region emission. The z ∼ 0 relationship between f DIG and ΣSFR implies that high-redshift galaxies have lower f DIG values than typical local systems, given their significantly higher typical ΣSFR. When an appropriate low-redshift benchmark with zero or minimal f DIG is used, high-redshift MOSDEF galaxies appear offset toward higher [S ii]λλ6717,6731/Hα and/or [O iii]λ 5007/Hβ. The joint shifts of high-redshift galaxies in the [S ii] and [N ii] BPT diagrams are best explained in terms of the harder spectra ionizing their star-forming regions at fixed nebular oxygen abundance (expected for chemically young galaxies), as opposed to large variations in N/O ratios or higher ionization parameters. The evolving mixture of H ii regions and diffuse ionized gas is an essential ingredient of our description of the interstellar medium over cosmic time.

The MOSDEF Survey: Broad Emission Lines at z = 1.4–3.8*

Abstract We present results from the MOSFIRE Deep Evolution Field survey on broad flux from the nebular emission lines Hα, [N ii], [O iii], Hβ, and [S ii]. The sample consists of 127 star-forming galaxies at 1.37 &lt; z &lt; 2.61 and 84 galaxies at 2.95 &lt; z &lt; 3.80. We decompose the emission lines using narrow and broad Gaussian components that we define as having FWHM &lt; 275 km s−1 and FWHM &gt; 300 km s−1, respectively, for both individual galaxies and stacks. For individual galaxies, broad emission is detected at &gt;3σ in &lt;10% of galaxies and the broad flux accounts for 10%–70% of the total flux. In the stacks, we find a slight increase in broad to narrow flux ratio with mass but note that we cannot reliably detect broad emission with FWHM &lt; 275 km s−1, which may be significant at low masses. When placed on the N2-BPT diagram ([O iii]/Hβ versus [N ii]/Hα), the broad components of the stacks are shifted toward higher [O iii]/Hβ and [N ii]/Hα ratios compared to the narrow component. We compare the location of the broad components to shock models and find that the broad component could be explained as a shocked outflow, but we do not rule out other possibilities, such as the presence of an AGN. We discuss the possible consequences of shocked emission on the galaxy location in emission line diagnostic diagrams and calculation of SFR. We attempt to estimate the mass outflow rate/star formation rate, but caution that our results strongly rely on the assumptions regarding the physical properties of the outflow.

The MOSDEF Survey: Significant Evolution in the Rest-frame Optical Emission Line Equivalent Widths of Star-forming Galaxies at z = 1.4–3.8

Abstract 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.