We present an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFRs) inferred from mid-infrared (MIR) data that are substantially elevated relative to those estimated from dust-corrected UV data. We use a sample of ∼200 galaxies and active galactic nuclei (AGNs) at 1.40 < z < 2.61 with 24 μm detections (rest-frame 8 μm) from MIPS/Spitzer. We find that the identification of MIR-excess galaxies strongly depends on the methodologies used to estimate IR luminosity (L1R) and to correct the UV light for dust attenuation. We find that extrapolations of the SFR from the observed 24 μm flux, using luminosity-dependent templates based on local galaxies, substantially overestimate in z ∼ 2 galaxies. By including Herschel observations and using a stellar-mass-dependent, luminosity-independent L1R, we obtain more reliable estimates of the SFR and a lower fraction of MIR-excess galaxies. Once stellar-mass selection biases are taken into account, we identify ∼24% of our galaxies as MIR excess. However, is not elevated in MIR-excess galaxies compared to MIR-normal galaxies, indicating that the intrinsic fraction of MIR excess may be lower. Using X-ray-, IR-, and optically selected AGNs in MOSDEF, we do not find a higher prevalence for AGNs in MIR-excess galaxies relative to MIR-normal galaxies. A stacking analysis of X-ray-undetected galaxies does not reveal a harder spectrum in MIR-excess galaxies relative to MIR-normal galaxies. Our analysis indicates that AGN activity does not contribute substantially to the MIR excess and instead implies that it is likely due to the enhanced polycyclic aromatic hydrocarbon emission.

Using observations from the first 2 yr of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we study 13 active galactic nucleus (AGN) driven outflows detected from a sample of 67 X-ray, IR, and/or optically selected AGNs at. The AGNs have bolometric luminosities of, including both quasars and moderate-luminosity AGNs. We detect blueshifted, ionized gas outflows in the Hβ, [O iii], Hα, and/or [N ii] emission lines of 19% of the AGNs, while only 1.8% of the MOSDEF galaxies have similarly detected outflows. The outflow velocities span ∼300 to 1000 km s-1. Eight of the 13 outflows are spatially extended on similar scales to the host galaxies, with spatial extents of 2.5-11.0 kpc. Outflows are detected uniformly across the star-forming main sequence, showing little trend with the host galaxy star formation rate. Line ratio diagnostics indicate that the outflowing gas is photoionized by the AGNs. We do not find evidence for positive AGN feedback, in either our small MOSDEF sample or a much larger Sloan Digital Sky Survey sample, using the BPT diagram. Given that a galaxy with an AGN is 10 times more likely to have a detected outflow, the outflowing gas is photoionized by the AGNs, and estimates of the mass and energy outflow rates indicate that stellar feedback is insufficient to drive at least some of these outflows; they are very likely to be AGN driven. The outflows have mass-loading factors of the order of unity, suggesting that they help regulate star formation in their host galaxies, though they may be insufficient to fully quench it.

We present the first spectroscopic measurement of multiple rest-frame optical emission lines at z > 4. During the MOSFIRE Deep Evolution Field survey, we observed the galaxy GOODSN-17940 with the Keck I/MOSFIRE spectrograph. The K-band spectrum of GOODSN-17940 includes significant detections of the [O ii]λλ3726,3729, [Ne iii]λ3869, and Hγ emission lines and a tentative detection of Hδ, indicating z spec = 4.4121. GOODSN-17940 is an actively star-forming z > 4 galaxy based on its K-band spectrum and broadband spectral energy distribution. A significant excess relative to the surrounding continuum is present in the Spitzer/IRAC channel 1 photometry of GOODSN-17940, due primarily to strong Hα emission with a rest-frame equivalent width of EW(Hα) = 1200. Based on the assumption of 0.5 Z o models and the Calzetti attenuation curve, GOODSN-17940 is characterized by . The Balmer decrement inferred from Hα/Hγ is used to dust correct the Hα emission, yielding . These M ∗ and star formation rate (SFR) values place GOODSN-17940 an order of magnitude in SFR above the z ∼ 4 star-forming "main sequence." Finally, we use the observed ratio of [Ne iii]/[O ii] to estimate the nebular oxygen abundance in GOODSN-17940, finding O/H ∼ 0.2 (O/H)o. Combining our new [Ne iii]/[O ii] measurement with those from stacked spectra at z ∼ 0, 2, and 3, we show that GOODSN-17940 represents an extension to z > 4 of the evolution toward higher [Ne iii]/[O ii] (i.e., lower O/H) at fixed stellar mass. It will be possible to perform the measurements presented here out to z ∼ 10 using the James Webb Space Telescope.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on the identification, selection biases, and host galaxy properties of 55 X-ray, IR, and optically selected active galactic nuclei (AGNs) at 1.4 < z < 3.8. We obtain rest-frame optical spectra of galaxies and AGNs and use the BPT diagram to identify optical AGNs. We examine the uniqueness and overlap of the AGNs identified at different wavelengths. There is a strong bias against identifying AGNs at any wavelength in low-mass galaxies, and an additional bias against identifying IR AGNs in the most massive galaxies. AGN hosts span a wide range of star formation rates (SFRs), similar to inactive galaxies once stellar mass selection effects are accounted for. However, we find (at ∼ 2-3σ significance) that IR AGNs are in less dusty galaxies with relatively higher SFR and optical AGNs in dusty galaxies with relatively lower SFR. X-ray AGN selection does not display a bias with host galaxy SFR. These results are consistent with those from larger studies at lower redshifts. Within star-forming galaxies, once selection biases are accounted for, we find AGNs in galaxies with similar physical properties as inactive galaxies, with no evidence for AGN activity in particular types of galaxies. This is consistent with AGNs being fueled stochastically in any star-forming host galaxy. We do not detect a significant correlation between SFR and AGN luminosity for individual AGN hosts, which may indicate the timescale difference between the growth of galaxies and their supermassive black holes.

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-1 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 at LAGN≳ 1046 erg s-1. 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.

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.

We present measurements of [SIII]$\lambda\lambda$9069,9531 for a sample of $z\sim1.5$ star-forming galaxies, the first sample with measurements of these lines at z>0.1. We employ the line ratio S$_{32}$$\equiv$[SIII]$\lambda\lambda$9069,9531/[SII]$\lambda\lambda$6716,6731 as a novel probe of evolving ISM conditions. Since this ratio includes the low-ionization line [SII], it is crucial that the effects of diffuse ionized gas (DIG) on emission-line ratios be accounted for in $z\sim0$ integrated galaxy spectra, or else that comparisons be made to samples of local HII regions in which DIG emission is not present. We find that S$_{32}$ decreases with increasing stellar mass at both $z\sim1.5$ and $z\sim0$, but that the dependence is weak suggesting S$_{32}$ has a very shallow anticorrelation with metallicity, in contrast with O$_{32}$ that displays a strong metallicity dependence. As a result, S$_{32}$ only mildly evolves with redshift at fixed stellar mass. The $z\sim1.5$ sample is systematicallty offset towards lower S$_{32}$ and higher [SII]/H$\alpha$ at fixed [OIII]/H$\beta$ relative to $z=0$ HII regions. By comparing to photoionization model grids, we find that such trends can be explained by a scenario in which the ionizing spectrum is harder at fixed O/H with increasing redshift, but are inconsistent with an increase in ionization parameter at fixed O/H. This analysis demonstrates the advantages of expanding beyond the strongest rest-optical lines for evolutionary studies, and the particular utility of [SIII] for characterizing evolving ISM conditions and stellar compositions. These measurements provide a basis for estimating [SIII] line strengths for high-redshift galaxies, a line that the James Webb Space Telescope will measure out to z~5.5.

We study the properties of 30 spectroscopically identified pairs of galaxies observed during the peak epoch of star formation in the universe. These systems are drawn from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at 1.4 ≤ z ≤ 3.8, and are interpreted as early-stage galaxy mergers. Galaxy pairs in our sample are identified as two objects whose spectra were collected on the same Keck/MOSFIRE spectroscopic slit. Accordingly, all pairs in the sample have projected separations R proj ≤ 60 kpc. The velocity separation for pairs was required to be Δv ≤ 500 km s -1 , which is a standard threshold for defining interacting galaxy pairs at low redshift. Stellar mass ratios in our sample range from 1.1 to 550, with 12 ratios closer than or equal to 3:1, the common definition of a "major merger." Studies of merging pairs in the local universe indicate an enhancement in star formation activity and deficit in gas-phase oxygen abundance relative to isolated galaxies of the same mass. We compare the MOSDEF pairs sample to a control sample of isolated galaxies at the same redshift, finding no measurable SFR enhancement or metallicity deficit at fixed stellar mass for the pairs sample. The lack of significant difference between the average properties of pairs and control samples appears in contrast to results from low-redshift studies, although the small sample size and lower signal-to-noise of the high-redshift data limit definitive conclusions on redshift evolution. These results are consistent with some theoretical works, suggesting a reduced differential effect of precoalescence mergers on galaxy properties at high redshift - specifically that precoalescence mergers do not drive strong starbursts.

We present ionized gas kinematics for 681 galaxies at z &tild; 1.4-3.8 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 (Σbar,e). Intrinsic dispersion correlates most with Σbar,e, and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at z &tild; 1.4-3.8 are baryon dominated within their effective radii (RE), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within RE, fbar, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, fbar does not appear to vary with stellar mass but increases with decreasing RE and increasing Σbar,e. For galaxies at z ≥ 2, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low-metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity (L X) per star formation rate (SFR) with redshift due to the decrease of metallicity (Z) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 z ∼ 2 star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for ∼1500 galaxies in the CANDELS fields at 1.37 < z < 3.80. Using Chandra data from the Chandra AEGIS-X Deep, Chandra Deep Field North, and Chandra Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and Z because the galaxies are too faint to be individually detected. In agreement with previous studies, the average L X/SFR of our z ∼ 2 galaxy sample is enhanced by ≈0.4-0.8 dex relative to local HMXB L X-SFR scaling relations. Splitting our sample by Z, we find that L X/SFR and Z are anticorrelated with 97% confidence. This observed Z dependence for HMXB-dominated galaxies is consistent with both the local L X-SFR-Z relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak owing to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of L X/SFR and the Z dependence of HMXBs.