We derive the 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 (MOSDEF) survey with H$\alpha$, H$\beta$, and [NII]$\lambda 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$\alpha$-to-H$\beta$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3\sigma$ 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.

We present ionized gas kinematics for 681 galaxies at $z\sim 1.4-3.8$ from the MOSFIRE Deep Evolution Field survey, measured using models which 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, $V/\sigma_{V,0}$. We find that $V/\sigma_{V,0}$ 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 ($\Sigma_{\mathrm{bar},e}$). Intrinsic dispersion correlates most with $\Sigma_{\mathrm{bar},e}$ and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at $z\sim 1.4-3.8$ are baryon dominated within their effective radii ($R_E$), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within $R_E$, $f_{\mathrm{bar}}$, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, $f_{\mathrm{bar}}$ does not appear to vary with stellar mass but increases with decreasing $R_E$ and increasing $\Sigma_{\mathrm{bar},e}$. For galaxies at $z\geq2$, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

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 present detections of [OIII]$\lambda$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<1$2+log(O/H$)<8.2$ and log(M$_*$/M$_{\odot})=7-10$. We find that strong-line calibrations based on local analogs 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 $\lesssim0.1$ dex at fixed M$_*$ and SFR from $z\sim0-2.2$. We employ photoionization models to constrain the ionization parameter and ionizing spectrum in the high-redshift sample. Stellar models with super-solar 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\sim0$ relation between ionization parameter and O/H, suggesting no evolution of this relation from $z\sim0$ to $z\sim2$. 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\sim0$, naturally explained by super-solar O/Fe at high redshift caused by rapid formation timescales. Given the extreme nature of our $z>1$ sample, the implications for representative $z\sim2$ galaxy samples at $\sim10^{10}$ M$_{\odot}$ are unclear, but similarities to $z>6$ galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization.

Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] 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 $\sim 400-3500 ~\mathrm{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 $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~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.

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_{\mathrm{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\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and 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_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.

We present results on the emission-line properties of 1.3<=z<=2.7 galaxies drawn from the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. Specifically, we use observations of the emission-line diagnostic diagram of [OIII]5007/Hb vs. [SII]6717,6731/Ha, i.e., the "[SII] BPT diagram," to gain insight into the physical properties of high-redshift star-forming regions. High-redshift MOSDEF galaxies are offset towards lower [SII]6717,6731/Ha at fixed [OIII]5007/Hb, relative to local galaxies from the Sloan Digital Sky Survey (SDSS). Furthermore, at fixed [OIII]5007/Hb, local SDSS galaxies follow a trend of decreasing [SII]6717,6731/Ha as the surface density of star formation (Sigma_SFR) increases. We explain this trend in terms of the decreasing fractional contribution from diffuse ionized gas (f_DIG) as Sigma_SFR increases in galaxies, which causes galaxy-integrated line ratios to shift towards the locus of pure HII-region emission. The z~0 relationship between f_DIG and Sigma_SFR implies that high-redshift galaxies have lower f_DIG values than typical local systems, given their significantly higher typical Sigma_SFR. When an appropriate low-redshift benchmark with zero or minimal f_DIG is used, high-redshift MOSDEF galaxies appear offset towards higher [SII]6717,6731/Ha and/or [OIII]5007/Hb. The joint shifts of high-redshift galaxies in the [SII] and [NII] 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 HII regions and DIG is an essential ingredient to our description of the ISM over cosmic time.

We investigate the infrared contribution from supermassive black hole activity versus host galaxy emission in the mid to far-infrared (IR) spectrum for a large sample of X-ray bright active galactic nuclei (AGN) residing in dusty, star-forming host galaxies. We select 703 AGN with L_X = 10^42-46 ergs/s at 0.1 < z < 5 from the Chandra XBootes X-ray Survey with rich multi-band observations in the optical to far-IR. This is the largest sample to date of X-ray AGN with mid and far-IR detections that uses spectral energy distribution (SED) decomposition to determine intrinsic AGN and host galaxy infrared luminosities. We determine weak or nonexistent relationships when averaging star-formation activity as a function of AGN activity, but see stronger positive trends when averaging L_X in bins of star-forming activity for AGN at low redshifts. We estimate an average dust covering factor of 33% based on infrared SEDs and bolometric AGN luminosity, corresponding to a Type 2 AGN population of roughly a third. We also see a population of AGN that challenge the inclination based unification model with individual dust covering factors that contradict the nuclear obscuration expected from observed X-ray hardness ratios. We see no strong connection between AGN fractions in the IR and corresponding total infrared, 24 um, or X-ray luminosities. The average rest-frame AGN contribution as a function of IR wavelength shows significant (~80%) contributions in the mid-IR that trail off at lambda > 30 um. Additionally, we provide a relation between observed L_X and pure AGN IR output for high-z AGN allowing future studies to estimate AGN infrared contribution using only observed X-ray flux density estimates.

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 \leq z \leq 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_{{\rm proj}}\leq 60$ kpc. The velocity separation for pairs was required to be $\Delta v \leq 500 \mbox{ 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 anenhancement 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 pre-coalescence mergers on galaxy properties at high redshift -- specifically that pre-coalescence mergers do not drive strong starbursts.

We use extensive spectroscopy from the MOSFIRE Deep Evolution Field (MOSDEF) survey to investigate the relationships between rest-frame optical emission line equivalent widths ($W$) and a number of galaxy and ISM characteristics for a sample of $1134$ star-forming galaxies at redshifts $1.4\lesssim z\lesssim 3.8$. We examine how the equivalent widths of [OII]$\lambda\lambda 3727, 3730$, H$\beta$, [OIII]$\lambda\lambda 4960, 5008$, [OIII]$+$H$\beta$, H$\alpha$, and H$\alpha$+[NII]$\lambda\lambda 6550, 6585$, depend on stellar mass, UV slope, age, star-formation rate (SFR) and specific SFR (sSFR), ionization parameter and excitation conditions (O32 and [OIII]/H$\beta$), gas-phase metallicity, and ionizing photon production efficiency ($\xi_{\rm ion}$). The trend of increasing $W$ with decreasing stellar mass is strongest for [OIII] (and [OIII]+H$\beta$). 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 $W({\rm [OIII]})$, have low stellar masses, blue UV slopes, young ages, high sSFRs, ISM line ratios indicative of high ionization parameters, high $\xi_{\rm 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.

Using the MOSDEF rest-frame optical spectroscopic survey, we investigate the star-formation histories (SFHs) of different galaxy types, ranging from actively star forming to quiescent at $1.4\leq~z\leq2.6$. SFHs are constrained utilizing stellar continuum spectroscopy, specifically through a combination of Balmer absorption lines, the 4000$\AA$ break, and the equivalent width of the H$\alpha$ emission line. To attain a sufficiently high signal-to-noise ratio (S/N) to conduct these measurements we stack spectra of galaxies with similar spectral types, as determined from their rest-frame $U-V$ and $V-J$ colors. We bin the MOSDEF sample into five spectral types, subdividing the quiescent and star-forming bins to better explore galaxies transitioning between the two. We constrain the average SFHs for each type, finding that quiescent and transitional galaxies in the MOSDEF sample are dominated by an SFH with an average star-formation timescale of $\tau\sim0.1-0.2$~Gyr. These findings contrast with measurements from the low-redshift Universe where, on average, galaxies form their stars over a more extended time period ($\tau>1$~Gyr). Furthermore, our spectral index measurements correlate with mass surface density for all spectral types. Finally, we compare the average properties of the galaxies in our transitional bins to investigate possible paths to quiescence, and speculate on the viability of a dusty post-starburst phase.

We investigate the nature of the relation among stellar mass, star-formation rate, and gas-phase metallicity (the M$_*$-SFR-Z relation) at high redshifts using a sample of 260 star-forming galaxies at $z\sim2.3$ from the MOSDEF survey. We present an analysis of the high-redshift M$_*$-SFR-Z relation based on several emission-line ratios for the first time. We show that a M$_*$-SFR-Z relation clearly exists at $z\sim2.3$. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of $z\sim0$ and $z\sim2.3$ galaxies, we find that $z\sim2.3$ galaxies have $\sim0.1$ dex lower metallicity at fixed M$_*$ and SFR. In the context of chemical evolution models, this evolution of the M$_*$-SFR-Z relation suggests an increase with redshift of the mass-loading factor at fixed M$_*$, 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 M$_*$ at $z\sim2.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 present an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFR) inferred from mid-infrared (MIR) data that is substantially elevated relative to that estimated from dust-corrected UV data. We use a sample of $\sim$200 galaxies and AGN at $1.40<z<2.61$ with 24 $\mu$m detections (rest-frame 8$\mu$m) from MIPS/\textit{Spitzer}. We find that the identification of MIR-excess galaxies strongly depends on the methodologies used to estimate IR luminosity ($\rm L_{IR}$) and to correct the UV light for dust attenuation. We find that extrapolations of the SFR from the observed 24 $\mu$m flux, using luminosity-dependent templates based on local galaxies, substantially overestimate $\rm L_{IR}$ in $z\sim2$ galaxies. By including \textit{Herschel} observations and using a stellar mass-dependent, luminosity-independent $\rm L_{IR}$, 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 $\sim24\%$ of our galaxies as MIR-excess. However, $\rm SFR_{H\alpha}$ 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 AGN in MOSDEF, we do not find a higher prevalence for AGN 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 PAH emission.

We combine spectroscopic measurements of H$\alpha$ and H$\beta$ and UV continuum photometry for a sample of 673 galaxies from the MOSFIRE Deep Evolution Field survey to constrain hydrogen ionizing photon production efficiencies ($\xi_{\rm ion}$, xi_ion) at z=1.4-2.6. We find average log(xi_ion/[Hz erg$^{-1}$])=25.06 (25.34), assuming the Calzetti (SMC) curve for the UV dust correction and a scatter of 0.28 dex in xi_ion distribution. After accounting for observational uncertainties and variations in dust attenuation, we conclude that the remaining scatter in xi_ion is likely dominated by galaxy-to-galaxy variations in stellar populations, including the slope and upper-mass cutoff of the initial mass function, stellar metallicity, star-formation burstiness, and stellar evolution (e.g., single/binary star evolution). Moreover, xi_ion is elevated in galaxies with high ionization states (high [OIII]/[OII]) and low oxygen abundances (low [NII]/H$\alpha$ and high [OIII]/H$\beta$) in the ionized ISM. However, xi_ion does not correlate with the offset from the z~0 star-forming locus in the BPT diagram, suggesting no change in the hardness of ionizing radiation accompanying the offset from the z~0 sequence. We also find that galaxies with blue UV spectral slopes ($\langle\beta\rangle$=-2.1) have elevated xi_ion by a factor of ~2 relative to the average xi_ion of the sample ($\langle\beta\rangle$=-1.4). If these blue galaxies are similar to those at z > 6, our results suggest that a lower Lyman continuum escape fraction is required for galaxies to maintain reionization, compared to the canonical xi_ion predictions from stellar population models. Furthermore, we demonstrate that even with robustly dust-corrected H$\alpha$, the UV dust attenuation can cause on average a ~0.3dex systematic uncertainty in xi_ion calculations.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on broad flux from the nebular emission lines H$\alpha$, [NII], [OIII], H$\beta$, and [SII]. The sample consists of 127 star-forming galaxies at $1.37 < z < 2.61$ and 84 galaxies at $2.95 < z < 3.80$. We decompose the emission lines using narrow ($\text{FWHM} < 275 \ \text{km s}^{-1}$) and broad ($\text{FWHM} > 300 \ \text{km s}^{-1}$) Gaussian components for individual galaxies and stacks. Broad emission is detected at $>3\sigma$ in $<10$% of galaxies and the broad flux accounts for 10-70% of the total flux. We find a slight increase in broad to narrow flux ratio with mass but note that we cannot reliably detect broad emission with $\text{FWHM} < 275 \ \text{km s}^{-1}$, which may be significant at low masses. Notably, there is a correlation between higher signal-to-noise (S/N) spectra and a broad component detection indicating a S/N dependence in our ability to detect broad flux. When placed on the N2-BPT diagram ([OIII]/H$\beta$ vs. [NII]/H$\alpha$) the broad components of the stacks are shifted towards higher [OIII]/H$\beta$ and [NII]/$\alpha$ 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 due to shocks, but we do not rule out other possibilities such as the presence of an AGN. We estimate the mass loading factor (mass outflow rate/star formation rate) assuming the broad component is a photoionized outflow and find that the mass loading factor increases as a function of mass which agrees with previous studies. We show that adding emission from shocked gas to $z\sim0$ SDSS spectra shifts galaxies towards the location of $z\sim2$ galaxies on several emission line diagnostic diagrams.

Using observations from the first two years of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we study 13 AGN-driven outflows detected from a sample of 67 X-ray, IR and/or optically-selected AGN at $z \sim 2$. The AGN have bolometric luminosities of $\sim10^{44}-10^{46} ~\mathrm{erg~s^{-1}}$, including both quasars and moderate-luminosity AGN. We detect blueshifted, ionized gas outflows in the H$\beta$ , [OIII], H$\alpha$ ~and/or [NII] emission lines of $19\%$ of the AGN, while only 1.8\% of the MOSDEF galaxies have similarly-detected outflows. The outflow velocities span $\sim$300 to 1000 km s$^{-1}$. Eight of the 13 outflows are spatially extended on similar scales as the host galaxies, with spatial extents of 2.5 to 11.0 kpc. Outflows are detected uniformly across the star-forming main sequence, showing little trend with the host galaxy SFR. Line ratio diagnostics indicate that the outflowing gas is photoionized by the AGN. We do not find evidence for positive AGN feedback, in either our small MOSDEF sample or a much larger SDSS sample, using the BPT diagram. Given that a galaxy with an AGN is ten times more likely to have a detected outflow, the outflowing gas is photoionzed by the AGN, 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 (MOSDEF) 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 [OII]$\lambda\lambda 3726,3729$, [NeIII]$\lambda3869$, and H$\gamma$ emission lines and a tentative detection of H$\delta$, indicating $z_{\rm{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$\alpha$ emission with a rest-frame equivalent width of $\mbox{EW(H}\alpha)=1200$ \AA. Based on the assumption of $0.5 Z_{\odot}$ models and the Calzetti attenuation curve, GOODSN-17940 is characterized by $M_*=5.0^{+4.3}_{-0.2}\times 10^9 M_{\odot}$. The Balmer decrement inferred from H$\alpha$/H$\gamma$ is used to dust correct the H$\alpha$ emission, yielding $\mbox{SFR(H}\alpha)=320^{+190}_{-140} M_{\odot}\mbox{ yr}^{-1}$. These $M_*$ and SFR values place GOODSN-17940 an order of magnitude in SFR above the $z\sim 4$ star-forming "main sequence." Finally, we use the observed ratio of [NeIII]/[OII] to estimate the nebular oxygen abundance in GOODSN-17940, finding $\mbox{O/H}\sim 0.2 \mbox{ (O/H)}_{\odot}$. Combining our new [NeIII]/[OII] measurement with those from stacked spectra at $z\sim 0, 2, \mbox{ and } 3$, we show that GOODSN-17940 represents an extension to $z>4$ of the evolution towards higher [NeIII]/[OII] (i.e., lower $\mbox{O/H}$) at fixed stellar mass. It will be possible to perform the measurements presented here out to $z\sim 10$ using the James Webb Space Telescope.

We present results on the variation of 7.7 micron Polycyclic Aromatic Hydrocarbon (PAH) emission in galaxies spanning a wide range in metallicity at z ~ 2. For this analysis, we use rest-frame optical spectra of 476 galaxies at 1.37 < z < 2.61 from the MOSFIRE Deep Evolution Field (MOSDEF) survey to infer metallicities and ionization states. Spitzer/MIPS 24 micron and Herschel/PACS 100 and 160 micron observations are used to derive rest-frame 7.7 micron luminosities (L(7.7)) and total IR luminosities (L(IR)), respectively. We find significant trends between the ratio of L(7.7) to L(IR) (and to dust-corrected SFR) and both metallicity and [OIII]/[OII] (O32) emission-line ratio. The latter is an empirical proxy for the ionization parameter. These trends indicate a paucity of PAH emission in low metallicity environments with harder and more intense radiation fields. Additionally, L(7.7)/L(IR) is significantly lower in the youngest quartile of our sample (ages of 500 Myr) compared to older galaxies, which may be a result of the delayed production of PAHs by AGB stars. The relative strength of L(7.7) to L(IR) is also lower by a factor of ~ 2 for galaxies with masses $M_* < 10^{10}M_{\odot}$, compared to the more massive ones. We demonstrate that commonly-used conversions of L(7.7) (or 24 micron flux density; f(24)) to L(IR) underestimate the IR luminosity by more than a factor of 2 at $M_*$ ~ $10^{9.6-10.0} M_{\odot}$. We adopt a mass-dependent conversion of L(7.7) to L(IR) with L(7.7)/L(IR)= 0.09 and 0.22 for $M_* < 10^{10}$ and $> 10^{10} M_{\odot}$, respectively. Based on the new scaling, the SFR-$M_*$ relation has a shallower slope than previously derived. Our results also suggest a higher IR luminosity density at z ~ 2 than previously measured, corresponding to a ~ 30% increase in the SFR density.

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 (AGN) at $1.4 < z < 3.8$. We obtain rest-frame optical spectra of galaxies and AGN and use the BPT diagram to identify optical AGN. We examine the uniqueness and overlap of the AGN identified at different wavelengths. There is a strong bias against identifying AGN at any wavelength in low mass galaxies, and an additional bias against identifying IR AGN in the most massive galaxies. AGN hosts span a wide range of star formation rate (SFR), similar to inactive galaxies once stellar mass selection effects are accounted for. However, we find (at $\sim 2-3\sigma$ significance) that IR AGN are in less dusty galaxies with relatively higher SFR and optical AGN 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 AGN in galaxies with similar physical properties as inactive galaxies, with no evidence for AGN activity in particular types of galaxies. This is consistent with AGN 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.

We present the first direct comparison between Balmer line and panchromatic SED-based SFRs for z~2 galaxies. For this comparison we used 17 star-forming galaxies selected from the MOSFIRE Deep Evolution Field (MOSDEF) survey, with $3\sigma$ detections for H$\alpha$ and at least two IR bands (Spitzer/MIPS 24$\mu$m and Herschel/PACS 100 and 160$\mu$m, and in some cases Herschel/SPIRE 250, 350, and 500$\mu$m). The galaxies have total IR (8-1000$\mu$m) luminosities of $\sim10^{11.4}-10^{12.4}\,\textrm{L}_\odot$ and star-formation rates (SFRs) of $\sim30-250\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We fit the UV-to-far-IR SEDs with flexible stellar population synthesis (FSPS) models - which include both stellar and dust emission - and compare the inferred SFRs with the SFR(H$\alpha$,H$\beta$) values corrected for dust attenuation using Balmer decrements. The two SFRs agree with a scatter of 0.17 dex. Our results imply that the Balmer decrement accurately predicts the obscuration of the nebular lines and can be used to robustly calculate SFRs for star-forming galaxies at z~2 with SFRs up to $\sim200\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We also use our data to assess SFR indicators based on modeling the UV-to-mid-IR SEDs or by adding SFR(UV) and SFR(IR), for which the latter is based on the mid-IR only or on the full IR SED. All these SFRs show a poorer agreement with SFR(H$\alpha$,H$\beta$) and in some cases large systematic biases are observed. Finally, we show that the SFR and dust attenuation derived from the UV-to-near-IR SED alone are unbiased when assuming a delayed exponentially declining star-formation history.

We present H$\alpha$ gas kinematics for 178 star-forming galaxies at z~2 from the MOSFIRE Deep Evolution Field survey. We have developed models to interpret the kinematic measurements from fixed-angle multi-object spectroscopy, using structural parameters derived from CANDELS HST/F160W imaging. For 35 galaxies we measure resolved rotation with a median $(V/\sigma_{V,0})_{R_E}=2.1$. We derive dynamical masses from the kinematics and sizes and compare them to baryonic masses, with gas masses estimated from dust-corrected H$\alpha$ star formation rates (SFRs) and the Kennicutt-Schmidt relation. When assuming that galaxies with and without observed rotation have the same median $(V/\sigma_{V,0})_{R_E}$, we find good agreement between the dynamical and baryonic masses, with a scatter of $\sigma_{RMS}=0.34$dex and a median offset of $\Delta\log_{10}M=0.04$dex. This comparison implies a low dark matter fraction (8% within an effective radius) for a Chabrier initial mass function (IMF), and disfavors a Salpeter IMF. Moreover, the requirement that $M_{dyn}/M_{baryon}$ should be independent of inclination yields a median value of $(V/\sigma_{V,0})_{R_E}=2.1$ for galaxies without observed rotation. If instead we treat the galaxies without detected rotation as early-type galaxies, the masses are also in reasonable agreement ($\Delta\log_{10}M=-0.07$dex, $\sigma_{RMS}=0.37$dex). The inclusion of gas masses is critical in this comparison; if gas masses are excluded there is an increasing trend of $M_{dyn}/M_{*}$ with higher specific SFR (SSFR). Furthermore, we find indications that $V/\sigma$ decreases with increasing H$\alpha$ SSFR for our full sample, which may reflect disk settling. We also study the Tully-Fisher relation and find that at fixed stellar mass $S_{0.5}=(0.5V_{2.2}^2+\sigma_{V,0}^2)^{1/2}$ was higher at earlier times. At fixed baryonic mass, we observe the opposite trend. [abridged]

We present results on the star-formation rate (SFR) versus stellar mass ($M_*$) relation (i.e., the "main sequence") among star-forming galaxies at 1.37$\leq$$z$$\leq$2.61 using the MOSFIRE Deep Evolution Field (MOSDEF) survey. Based on a sample of 261 galaxies with H$\alpha$ and H$\beta$ spectroscopy, we have estimated robust dust-corrected instantaneous SFRs over a large range in $M_*$ ($\sim10^{9.5}-10^{11.5}M_\odot$). We find a correlation between log(SFR(H$\alpha$)) and log($M_*$) with a slope of 0.65$\pm$0.08 (0.58$\pm$0.10) at 1.4<z<2.6 (2.1<z<2.6). We find that different assumptions for the dust correction, such as using the color-excess of the stellar continuum to correct the nebular lines, sample selection biases against red star-forming galaxies, and not accounting for Balmer absorption can yield steeper slopes of the log(SFR)-log($M_*$) relation. Our sample is immune from these biases as it is rest-frame optically selected, H$\alpha$ and H$\beta$ are corrected for Balmer absorption, and the H$\alpha$ luminosity is dust-corrected using the nebular color-excess computed from the Balmer decrement. The scatter of the log(SFR(H$\alpha$))-log($M_*$) relation, after accounting for the measurement uncertainties, is 0.31 dex at 2.1<z<2.6, which is 0.05 dex larger than the scatter in log(SFR(UV))-log($M_*$). Based on comparisons to a simulated SFR-$M_*$ relation with some intrinsic scatter, we argue that in the absence of direct measurements of galaxy-to-galaxy variations in the attenuation/extinction curves and the IMF, one cannot use the difference in the scatter of the SFR(H$\alpha$)- and SFR(UV)-$M_*$ relations to constrain the stochasticity of star formation in high-redshift galaxies.

We study the evidence for a connection between active galactic nuclei (AGN) fueling and star formation by investigating the relationship between the X-ray luminosities of AGN and the star formation rates (SFRs) of their host galaxies. We identify a sample of 309 AGN with $10^{41}<L_\mathrm{X}<10^{44} $ erg s$^{-1}$ at $0.2 < z < 1.2$ in the PRIMUS redshift survey. We find AGN in galaxies with a wide range of SFR at a given $L_X$. We do not find a significant correlation between SFR and the observed instantaneous $L_X$ for star forming AGN host galaxies. However, there is a weak but significant correlation between the mean $L_\mathrm{X}$ and SFR of detected AGN in star forming galaxies, which likely reflects that $L_\mathrm{X}$ varies on shorter timescales than SFR. We find no correlation between stellar mass and $L_\mathrm{X}$ within the AGN population. Within both populations of star forming and quiescent galaxies, we find a similar power-law distribution in the probability of hosting an AGN as a function of specific accretion rate. Furthermore, at a given stellar mass, we find a star forming galaxy $\sim2-3$ more likely than a quiescent galaxy to host an AGN of a given specific accretion rate. The probability of a galaxy hosting an AGN is constant across the main sequence of star formation. These results indicate that there is an underlying connection between star formation and the presence of AGN, but AGN are often hosted by quiescent galaxies.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.