• We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$\alpha$ equivalent width (EW), and various galaxy properties, we study the evolutionary trends among Ly$\alpha$, low- and high-ionization interstellar (LIS and HIS) absorption features, and integrated galaxy properties. We also examine the redshift evolution of Ly$\alpha$ and LIS absorption kinematics, and fine-structure emission EWs. The connections among the strengths of Ly$\alpha$, LIS lines, and dust extinction are redshift-independent, as is the decoupling of Ly$\alpha$ and HIS line strengths, and the bulk outflow kinematics as traced by LIS lines. Stronger Ly$\alpha$ emission is observed at higher redshift at fixed UV luminosity, stellar mass, SFR, and age. Much of this variation in average Ly$\alpha$ strength with redshift, and the variation in Ly$\alpha$ strength at fixed redshift, can be explained in terms of variations in neutral gas covering fraction and/or dust content in the ISM and CGM. However, based on the connection between Ly$\alpha$ and CIII] emission strengths, we additionally find evidence for variations in the intrinsic production rate of Ly$\alpha$ photons at the highest Ly$\alpha$ EWs. The challenge now is to understand the observed evolution in neutral gas covering fraction and dust extinction within a coherent model for galaxy formation, and make robust predictions for the escape of ionizing radiation at z > 6.
  • 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 investigate the nature and evolution of large-scale structure within the SSA22 protocluster region at $z=3.09$ using cosmological simulations. A redshift histogram constructed from current spectroscopic observations of the SSA22 protocluster reveals two separate peaks at $z = 3.065$ (blue) and $z = 3.095$ (red). Based on these data, we report updated overdensity and mass calculations for the SSA22 protocluster. We find $\delta_{b,gal}=4.8 \pm 1.8$, $\delta_{r,gal}=9.5 \pm 2.0$ for the blue and red peaks, respectively, and $\delta_{t,gal}=7.6\pm 1.4$ for the entire region. These overdensities correspond to masses of $M_b = (0.76 \pm 0.17) \times 10^{15} h^{-1} M_{\odot}$, $M_r = (2.15 \pm 0.32) \times 10^{15} h^{-1} M_{\odot}$, and $M_t=(3.19 \pm 0.40) \times 10^{15} h^{-1} M_{\odot}$ for the red, blue, and total peaks, respectively. We use the Small MultiDark Planck (SMDPL) simulation to identify comparably massive $z\sim 3$ protoclusters, and uncover the underlying structure and ultimate fate of the SSA22 protocluster. For this analysis, we construct mock redshift histograms for each simulated $z\sim 3$ protocluster, quantitatively comparing them with the observed SSA22 data. We find that the observed double-peaked structure in the SSA22 redshift histogram corresponds not to a single coalescing cluster, but rather the proximity of a $\sim 10^{15}h^{-1} M_{\odot}$ protocluster and at least one $>10^{14} h^{-1} M_{\odot}$ cluster progenitor. Such associations in the SMDPL simulation are easily understood within the framework of hierarchical clustering of dark matter halos. We finally find that the opportunity to observe such a phenomenon is incredibly rare, with an occurrence rate of $7.4h^3 \mbox{ Gpc}^{-3}$.
  • Galaxy metallicity scaling relations provide a powerful tool for understanding galaxy evolution, but obtaining unbiased global galaxy gas-phase oxygen abundances requires proper treatment of the various line-emitting sources within spectroscopic apertures. We present a model framework that treats galaxies as ensembles of HII and diffuse ionized gas (DIG) regions of varying metallicities. These models are based upon empirical relations between line ratios and electron temperature for HII regions, and DIG strong-line ratio relations from SDSS-IV MaNGA IFU data. Flux-weighting effects and DIG contamination can significantly affect properties inferred from global galaxy spectra, biasing metallicity estimates by more than 0.3 dex in some cases. We use observationally-motivated inputs to construct a model matched to typical local star-forming galaxies, and quantify the biases in strong-line ratios, electron temperatures, and direct-method metallicities as inferred from global galaxy spectra relative to the median values of the HII region distributions in each galaxy. We also provide a generalized set of models that can be applied to individual galaxies or galaxy samples in atypical regions of parameter space. We use these models to correct for the effects of flux-weighting and DIG contamination in the local direct-method mass-metallicity and fundamental metallicity relations, and in the mass-metallicity relation based on strong-line metallicities. Future photoionization models of galaxy line emission need to include DIG emission and represent galaxies as ensembles of emitting regions with varying metallicity, instead of as single HII regions with effective properties, in order to obtain unbiased estimates of key underlying physical properties.
  • 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.
  • 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.
  • The CIII]{\lambda}{\lambda}1907, 1909 rest-frame UV emission doublet has recently been detected in galaxies during the epoch of reionization (z > 6), with high equivalent width (>10 {\AA}, rest frame). Currently, it is possible to obtain much more detailed information for star-forming galaxies at significantly lower redshift. Accordingly, studies of their far-UV spectra are useful for understanding the factors modulating the strength of CIII] emission. We present the first statistical sample of CIII] emission measurements in star-forming galaxies at z ~ 1. Our sample is drawn from the DEEP2 survey and spans the redshift 0.64 < z < 1.35 (<z> = 1.08). We find that the median equivalent width (EW) of individual CIII] detections in our sample (1.30 {\AA}) is much smaller than the typical value observed thus far at z > 6. Furthermore, out of 184 galaxies with coverage of CIII], only 40 have significant detections. Galaxies with individual CIII] detections have bluer colors and lower luminosities on average than those without, implying that strong CIII] emitters are in general young and low-mass galaxies without significant dust extinction. Using stacked spectra, we further investigate how CIII] strength correlates with multiple galaxy properties (M_B, U-B, M*, star-formation rate, specific star-formation rate) and rest-frame near-UV (FeII* and MgII) and optical [OIII] and H{\beta}) emission line strengths. These results provide a detailed picture of the physical environment in star-forming galaxies at z ~ 1, and motivate future observations of strong CIII] emitters at similar redshifts.
  • Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.
  • We present observations of Q1549-C25, an ~L* star-forming galaxy at z=3.15 for which Lyman-continuum (LyC) radiation is significantly detected in deep Keck/LRIS spectroscopy. We find no evidence for contamination from a lower-redshift interloper close to the line of sight in the high signal-to-noise spectrum of Q1549-C25. Furthermore, the morphology of Q1549-C25 in V_606, J_125, and H_160 Hubble Space Telescope (HST) imaging reveals that the object consists of a single, isolated component within 1". In combination, these data indicate Q1549-C25 as a clean spectroscopic detection of LyC radiation, only the second such object discovered to date at z~3. We model the spectral energy distribution (SED) of Q1549-C25, finding evidence for negligible dust extinction, an age (assuming continuous star formation) of ~1 Gyr, and a stellar mass of M_*=7.9x10^9 M_sun. Although it is not possible to derive strong constraints on the absolute escape fraction of LyC emission, f_esc(LyC), from a single object, we use simulations of intergalactic and circumgalactic absorption to infer f_esc(LyC)>=0.51 at 95% confidence. The combination of deep Keck/LRIS spectroscopy and HST imaging is required to assemble a larger sample of objects like Q1549-C25, and obtain robust constraints on the average f_esc(LyC) at z~3 and beyond.
  • We present the first statistical sample of rest-frame far-UV spectra of star-forming galaxies at z~1. These spectra are unique in that they cover the high-ionization CIV{\lambda}{\lambda}1548, 1550 doublet. We also detect low-ionization features such as SiII{\lambda}1527, FeII{\lambda}1608, AlII{\lambda}1670, NiII{\lambda}{\lambda}1741, 1751 and SiII{\lambda}1808, and intermediate-ionization features from AlIII{\lambda}{\lambda}1854, 1862. Comparing the properties of absorption lines of lower- and higher- ionization states provides a window into the multi-phase nature of circumgalactic gas. Our sample is drawn from the DEEP2 survey and spans the redshift range 1.01 < z < 1.35 (<z> = 1.25). By isolating the interstellar CIV absorption from the stellar P-Cygni wind profile we find that 69% of the CIV profiles are blueshifted with respect to the systemic velocity. Furthermore, CIV shows a small but significant blueshift relative to FeII (offset of the best-fit linear regression -76 $\pm$ 26 km/s). At the same time, the CIV blueshift is on average comparable to that of MgII{\lambda}{\lambda}2796, 2803. At this point, in explaining the larger blueshift of CIV absorption at the ~ 3-sigma level, we cannot distinguish between the faster motion of highly-ionized gas relative to gas traced by FeII, and filling in on the red side from resonant CIV emission. We investigate how far-UV interstellar absorption kinematics correlate with other galaxy properties using stacked spectra. These stacking results show a direct link between CIV absorption and the current SFR, though we only observe small velocity differences among different ionization states tracing the outflowing ISM.
  • We present measurements of the electron-temperature based oxygen abundance for a highly star-forming galaxy at z=3.08, COSMOS-1908. This is the highest redshift at which [OIII]$\lambda$4363 has been detected, and the first time that this line has been measured at z>2. We estimate an oxygen abundance of 12+log(O/H)$=8.00^{+0.13}_{-0.14}$. This galaxy is a low-mass ($10^{9.3}$ M$_{\odot}$), highly star-forming ($\sim50$ M$_{\odot}$ yr$^{-1}$) system that hosts a young stellar population ($\sim160$ Myr). We investigate the physical conditions of the ionized gas in COSMOS-1908 and find that this galaxy has a high ionization parameter, little nebular reddening ($E(B-V)_{\rm gas}<0.14$), and a high electron density ($n_e\sim500$ cm$^{-3}$). We compare the ratios of strong oxygen, neon, and hydrogen lines to the direct-method oxygen abundance for COSMOS-1908 and additional star-forming galaxies at z=0-1.8 with [OIII]$\lambda$4363 measurements, and show that galaxies at z$\sim$1-3 follow the same strong-line correlations as galaxies in the local universe. This agreement suggests that the relationship between ionization parameter and O/H is similar for z$\sim$0 and high-redshift galaxies. These results imply that metallicity calibrations based on lines of oxygen, neon, and hydrogen do not strongly evolve with redshift and can reliably estimate abundances out to z$\sim$3, paving the way for robust measurements of the evolution of the mass-metallicity relation to high redshift.
  • Star-forming galaxies form a sequence in the [OIII]/H-beta vs. [NII]/H-alpha diagnostic diagram, with low metallicity, highly ionized galaxies falling in the upper left corner. Drawing from a large sample of UV-selected star-forming galaxies at z~2 with rest-frame optical nebular emission line measurements from Keck-MOSFIRE, we select the extreme ~5% of the galaxies lying in this upper left corner, requiring log([NII]/H-alpha) <= 1.1 and log([OIII]/H-beta) >= 0.75. These cuts identify galaxies with 12 + log(O/H) <~ 8.0, when oxygen abundances are measured via the O3N2 diagnostic. We study the Ly-alpha properties of the resulting sample of 14 galaxies. The mean (median) rest-frame Ly-alpha equivalent width is 39 (36) A, and 11 of the 14 objects (79%) are Ly-alpha-emitters (LAEs) with W_Lya > 20 A. We compare the equivalent width distribution of a sample of 522 UV-selected galaxies at 2.0<z<2.6 identified without regard to their optical line ratios; this sample has mean (median) Ly-alpha equivalent width -1 (-4) A, and only 9% of these galaxies qualify as LAEs. The extreme galaxies typically have lower attenuation at Ly-alpha than those in the comparison sample, and have ~50% lower median oxygen abundances. Both factors are likely to facilitate the escape of Ly-alpha: in less dusty galaxies Ly-alpha photons are less likely to be absorbed during multiple scatterings, while the harder ionizing spectrum and higher ionization parameter associated with strong, low metallicity star formation may reduce the covering fraction or column density of neutral hydrogen, further easing Ly-alpha escape. The use of nebular emission line ratios may prove useful in the identification of galaxies with low opacity to Ly-alpha photons across a range of redshifts.
  • We present the results of a densely sampled spectroscopic survey of the SSA22 protocluster at $z\approx 3.09$. Our sample with Keck/LRIS spectroscopy includes 106 Ly$\alpha$ Emitters (LAEs) and 40 Lyman Break Galaxies (LBGs) at $z=3.05-3.12$. These galaxies are contained within the $9'\times9'$ region in which the protocluster was discovered, which also hosts the maximum galaxy overdensity in the SSA22 region. The redshift histogram of our spectroscopic sample reveals two distinct peaks, at $z=3.069$ (blue, 43 galaxies) and $z=3.095$ (red, 103 galaxies). Furthermore, objects in the blue and red peaks are segregated on the sky, with galaxies in the blue peak concentrating towards the western half of the field. These results suggest that the blue and red redshift peaks represent two distinct structures in physical space. Although the double-peaked redshift histogram is traced in the same manner by LBGs and LAEs, and brighter and fainter galaxies, we find that nine out of 10 X-ray AGNs in SSA22, and all seven spectroscopically-confirmed giant Ly$\alpha$ "blobs," reside in the red peak. We combine our dataset with sparsely sampled spectroscopy from the literature over a significantly wider area, finding preliminary evidence that the double-peaked structure in redshift space extends beyond the region of our dense spectroscopic sampling. In order to fully characterize the three-dimensional structure, dynamics, and evolution of large-scale structure in the SSA22 overdensity, we require the measurement of large samples of LAE and LBG redshifts over a significantly wider area, as well as detailed comparisons with cosmological simulations of massive cluster formation.
  • 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.
  • High redshift star-forming galaxies are likely responsible for the reionization of the Universe, yet direct detection of their escaping ionizing (Lyman continuum) photons has proven to be extremely challenging. In this study, we search for escaping Lyman continuum of the Cosmic Horseshoe, a gravitationally lensed, star-forming galaxy at z=2.38 with a large magnification of $\sim24$. Transmission at wavelengths of low ionization interstellar absorption lines in the rest-frame ultraviolet suggest a patchy, partially transparent interstellar medium. This makes it an ideal candidate for direct detection of the Lyman continuum. We obtained a 10-orbit Hubble near-UV image using the WFC3/UVIS F275W filter that probes wavelengths just below the Lyman limit at the redshift of the Horseshoe in an attempt to detect escaping Lyman continuum radiation. After fully accounting for the uncertainties in the opacity of the intergalactic medium as well as accounting for the charge transfer inefficiency in the WFC3 CCDs, we find a $3 \sigma$ upper-limit for the relative escape fraction of $f_{esc,rel}<0.08$. This value is a factor of five lower than the value (0.4) predicted by the 40\% transmission in the low-ion absorption lines. We discuss the possible causes for this discrepancy and consider the implications for future attempts at both direct Lyman continuum detection as well as indirect estimates of the escape fraction.
  • 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.
  • Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the physical conditions of star-forming regions in $z\sim2.3$ galaxies, specifically the electron density and ionization state. From measurements of the [O II]$\lambda\lambda$3726,3729 and [S II]$\lambda\lambda$6716,6731 doublets, we find a median electron density of $\sim250$ cm$^{-3}$ at $z\sim2.3$, an increase of an order of magnitude compared to measurements of galaxies at $z\sim0$. While $z\sim2.3$ galaxies are offset towards significantly higher O$_{32}$ values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O$_{32}$ vs. R$_{23}$ and O3N2 diagrams. Based on these results, we propose that $z\sim2.3$ star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and $z\sim2.3$ galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from $z\sim0$ to $z\sim2.3$. We find that $z\sim2.3$ galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving [N II]$\lambda$6584. We conclude that the offset of $z\sim2.3$ galaxies from the local star-forming sequence in the [N II] BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at $z\sim0$ and $z\sim2$.
  • We present $U_{336}V_{606}J_{125}H_{160}$ follow-up $HST$ observations of 16 $z\sim3$ candidate LyC emitters in the HS1549+1919 field. With these data, we obtain high spatial-resolution photometric redshifts of all sub-arcsecond components of the LyC candidates in order to eliminate foreground contamination and identify robust candidates for leaking LyC emission. Of the 16 candidates, we find one object with a robust LyC detection that is not due to foreground contamination. This object (MD5) resolves into two components; we refer to the LyC-emitting component as MD5b. MD5b has an observed 1500\AA\ to 900\AA\ flux-density ratio of $(F_{UV}/F_{LyC})_{obs}=4.0\pm2.0$, compatible with predictions from stellar population synthesis models. Assuming minimal IGM absorption, this ratio corresponds to a relative (absolute) escape fraction of $f_{esc,rel}^{MD5b}=75-100$% ($f_{esc,abs}^{MD5b}=14-19$%). The stellar population fit to MD5b indicates an age of $\lesssim50$Myr, which is in the youngest 10% of the $HST$ sample and the youngest third of typical $z\sim3$ Lyman break galaxies, and may be a contributing factor to its LyC detection. We obtain a revised, contamination-free estimate for the comoving specific ionizing emissivity at $z=2.85$, indicating (with large uncertainties) that star-forming galaxies provide roughly the same contribution as QSOs to the ionizing background at this redshift. Our results show that foreground contamination prevents ground-based LyC studies from obtaining a full understanding of LyC emission from $z\sim3$ star-forming galaxies. Future progress in direct LyC searches is contingent upon the elimination of foreground contaminants through high spatial-resolution observations, and upon acquisition of sufficiently deep LyC imaging to probe ionizing radiation in high-redshift galaxies.
  • We use Hubble Space Telescope Wide-Field Camera 3 (HST/WFC3) rest-frame optical imaging to select a pilot sample of star-forming galaxies in the redshift range z = 2.00-2.65 whose multi-component morphologies are consistent with expectations for major mergers. We follow up this sample of major merger candidates with Keck/NIRSPEC longslit spectroscopy obtained in excellent seeing conditions (FWHM ~ 0.5 arcsec) to obtain Halpha-based redshifts of each of the morphological components in order to distinguish spectroscopic pairs from false pairs created by projection along the line of sight. Of six pair candidates observed, companions (estimated mass ratios 5:1 and 7:1) are detected for two galaxies down to a 3sigma limiting emission-line flux of ~ 10^{-17} erg/s/cm2. This detection rate is consistent with a ~ 50% false pair fraction at such angular separations (1-2 arcsec), and with recent claims that the star-formation rate (SFR) can differ by an order of magnitude between the components in such mergers. The two spectroscopic pairs identified have total SFR, SFR surface densities, and stellar masses consistent on average with the overall z ~ 2 star forming galaxy population.
  • 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.
  • We present results on the dust attenuation curve of z~2 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies with nebular spectroscopic redshifts in the range z= 1.36-2.59 and high S/N measurements of, or upper limits on, the H-alpha and H-beta emission lines obtained with Keck/MOSFIRE. We construct composite SEDs of galaxies in bins of specific SFR and Balmer optical depth in order to directly constrain the dust attenuation curve from the UV through near-IR for typical star-forming galaxies at high redshift. Our results imply an attenuation curve that is very similar to the SMC extinction curve at wavelengths redward of 2500 Angstroms. At shorter wavelengths, the shape of the curve is identical to that of the Calzetti relation, but with a lower normalization (R_V). Hence, the new attenuation curve results in SFRs that are ~20% lower, and log stellar masses that are 0.16 dex lower, than those obtained with the Calzetti attenuation curve. Moreover, we find that the difference in the reddening---and the total attenuation---of the ionized gas and stellar continuum correlates strongly with SFR, such that for dust-corrected SFRs larger than 20 Msun/yr assuming a Chabrier IMF, the nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple model that can account for these trends is one in which the UV through optical stellar continuum is dominated by a population of less reddened stars, while the nebular line and bolometric luminosities become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based SFRs may underestimate the total SFR at even modest levels of ~20 Msun/yr. [Abridged]
  • We use a sample of 262 spectroscopically confirmed star-forming galaxies at redshifts $2.08\leq z\leq 2.51$ to compare H$\alpha$, UV, and IR star-formation-rate diagnostics and to investigate the dust properties of the galaxies. At these redshifts, the H$\alpha$ line shifts to the $K_{s}$-band. By comparing $K_{s}$-band photometry to underlying stellar population model fits to other UV, optical, and near-infrared data, we infer the H$\alpha$ flux for each galaxy. We obtain the best agreement between H$\alpha$- and UV-based SFRs if we assume that the ionized gas and stellar continuum are reddened by the same value and that the Calzetti attenuation curve is applied to both. Aided with MIPS 24$\mu$m data, we find that an attenuation curve steeper than the Calzetti curve is needed to reproduce the observed IR/UV ratios of galaxies younger than 100 Myr. Furthermore, using the bolometric star-formation rate inferred from the UV and mid-IR data (SFR$_{IR}$+SFR$_{UV}$), we calculated the conversion between the H$\alpha$ luminosity and SFR to be $(7.5\pm1.3) \times 10^{-42}$ for a Salpeter IMF, which is consistent with the Kennicutt (1998) conversion. The derived conversion factor is independent of any assumption of the dust correction and is robust to stellar population model uncertainties.
  • We present results on the excitation properties of z~2.3 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) Survey. With its coverage of the full suite of strong rest-frame optical emission lines, MOSDEF provides an unprecedented view of the rest-frame optical spectra of a representative sample of distant star-forming galaxies. We investigate the locations of z~2.3 MOSDEF galaxies in multiple emission-line diagnostic diagrams. These include the [OIII]/Hb vs. [NII]/Ha and [OIII]/Hb vs. [SII]/Ha "BPT" diagrams, as well as the O_32 vs. R_23 excitation diagram. We recover the well-known offset in the star-forming sequence of high-redshift galaxies in the [OIII]/Hb vs. [NII]/Ha BPT diagram relative to SDSS star-forming galaxies. However, the shift for our rest-frame optically selected sample is less significant than for rest-frame-UV selected and emission-line selected galaxies at z~2. Furthermore, we find that the offset is mass-dependent, only appearing within the low-mass half of the z~2.3 MOSDEF sample, where galaxies are shifted towards higher [NII]/Ha at fixed [OIII]/Hb. Within the [OIII]/Hb vs. [SII]/Ha and O_32 vs. R_23 diagrams, we find that z~2.3 galaxies are distributed like local ones, and therefore attribute the shift in the [OIII]/Hb vs. [NII]/Ha BPT diagram to elevated N/O abundance ratios among lower-mass (M_*<10^10 M_sun) high-redshift galaxies. The variation in N/O ratios calls into question the use at high redshift of oxygen abundance indicators based on nitrogen lines, but the apparent invariance with redshift of the excitation sequence in the O_32 vs. R_23 diagram paves the way for using the combination of O_32 and R_23 as an unbiased metallicity indicator over a wide range in redshift. This indicator will allow for an accurate characterization of the shape and normalization of the mass-metallicity relationship over more than 10 Gyr.
  • We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on rest-frame optical AGN identification and completeness at z~2.3. With our sample of 50 galaxies and 10 X-ray and IR-selected AGN with measured H-beta, [OIII], H-alpha, and [NII] emission lines, we investigate the location of AGN in the BPT, MEx (mass-excitation), and CEx (color-excitation) diagrams. We find that the BPT diagram works well to identify AGN at z~2.3 and that the z~0 AGN/star-forming galaxy classifications do not need to shift substantially at z~2.3 to robustly separate these populations. However, the MEx diagram fails to identify all of the AGN identified in the BPT diagram, and the CEx diagram is substantially contaminated at high redshift. We further show that AGN samples selected using the BPT diagram have selection biases in terms of both host stellar mass and stellar population, in that AGN in low mass and/or high specific star formation rate galaxies are difficult to identify using the BPT diagram. These selection biases become increasingly severe at high redshift, such that optically-selected AGN samples at high redshift will necessarily be incomplete. We also find that the gas in the narrow-line region appears to be more enriched than gas in the host galaxy for at least some MOSDEF AGN. However, AGN at z~2 are generally less enriched than local AGN with the same host stellar mass.