We study the incidence of nuclear activity in a large sample of massive post-starburst galaxies at z~0.7 selected from the Sloan Digital Sky Survey, and identify active galactic nuclei based on radio continuum and optical emission lines. Over our mass range of 10^10.6-10^11.5 Msun, the incidence of radio activity is weakly dependent on stellar mass and independent of stellar age, while radio luminosity depends strongly on stellar mass. Optical nuclear activity incidence depends most strongly on the Dn4000 line index, a proxy for stellar age, with an active fraction that is ~ten times higher in the youngest versus oldest post-starburst galaxies. Since a similar trend is seen between age and molecular gas fractions, we argue that, like in local galaxies, the age trend reflects a peak in available fueling rather than feedback from the central black hole on the surrounding galaxy.

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 a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.

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 use Hubble Space Telescope (HST) imaging and near-infrared spectroscopy from Keck/MOSFIRE to study the sub-structure around the progenitor of a Milky Way-mass galaxy in the Hubble Frontier Fields (HFF). Specifically, we study an $r_e = 40^{+70}_{-30}$pc, $M_{\star} \sim 10^{8.2} M_{\odot}$ rest-frame ultra-violet luminous "clump" at a projected distance of $\sim$100~pc from a $M_{\star} \sim 10^{9.8}$M$_{\odot}$ galaxy at $z = 2.36$ with a magnification $\mu = 5.21$. We measure the star formation history of the clump and galaxy by jointly modeling the broadband spectral energy distribution from HST photometry and H$\alpha$ from MOSFIRE spectroscopy. Given our inferred properties (e.g., mass, metallicity, dust) of the clump and galaxy, we explore scenarios in which the clump formed \emph{in-situ} (e.g., a star forming complex) or \emph{ex-situ} (e.g., a dwarf galaxy being accreted). If it formed \emph{in-situ}, we conclude that the clump is likely a single entity as opposed to a aggregation of smaller star clusters, making it one of the most dense star clusters cataloged. If it formed \emph{ex-situ}, then we are witnessing an accretion event with a 1:40 stellar mass ratio. However, our data alone are not informative enough to distinguish between \emph{in-situ} and \emph{ex-situ} scenarios to a high level of significance. We posit that the addition of high-fidelity metallicity information, such as [OIII]4363\AA, which can be detected at modest S/N with only a few hours of JWST/NIRSpec time, may be a powerful discriminant. We suggest that studying larger samples of moderately lensed sub-structures across cosmic time can provide unique insight into the hierarchical formation of galaxies like the Milky Way.

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.

We use high-resolution, multi-band imaging of ~16,500 galaxies in the CANDELS fields at 0 < z < 2.5 to study the evolution of color gradients and half-mass radii over cosmic time. We find that galaxy color gradients at fixed mass evolve rapidly between z~2.5 and z~1, but remain roughly constant below z~1. This result implies that the sizes of both star-forming and quiescent galaxies increase much more slowly than previous studies found using half-light radii. The half-mass radius evolution of quiescent galaxies is fully consistent with a model which uses observed minor merger rates to predict the increase in sizes due to the accretion of small galaxies. Progenitor bias may still contribute to the growth of quiescent galaxies, particularly if we assume a slower timescale for the minor merger growth model. The slower half-mass radius evolution of star-forming galaxies is in tension with cosmological simulations and semi-analytic galaxy models. Further detailed, consistent comparisons with simulations are required to place these results in context.

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.

The chemical composition of galaxies has been measured out to z~4. However, nearly all studies beyond z~0.7 are based on strong-line emission from HII regions within star-forming galaxies. Measuring the chemical composition of distant quiescent galaxies is extremely challenging, as the required stellar absorption features are faint and shifted to near-infrared wavelengths. Here, we present ultra-deep rest-frame optical spectra of five massive quiescent galaxies at z~1.4, all of which show numerous stellar absorption lines. We derive the abundance ratios [Mg/Fe] and [Fe/H] for three out of five galaxies; the remaining two galaxies have too young luminosity-weighted ages to yield robust measurements. Similar to lower-redshift findings, [Mg/Fe] appears positively correlated with stellar mass, while [Fe/H] is approximately constant with mass. These results may imply that the stellar mass-metallicity relation was already in place at z~1.4. While the [Mg/Fe]-mass relation at z~1.4 is consistent with the z<0.7 relation, [Fe/H] at z~1.4 is ~0.2 dex lower than at z<0.7. With a [Mg/Fe] of 0.44(+0.08,-0.07) the most massive galaxy may be more alpha-enhanced than similar-mass galaxies at lower redshift, but the offset is less significant than the [Mg/Fe] of 0.6 previously found for a massive galaxy at z=2.1. Nonetheless, these results combined may suggest that [Mg/Fe] in the most massive galaxies decreases over time, possibly by accreting low-mass, less alpha-enhanced galaxies. A larger galaxy sample is needed to confirm this scenario. Finally, the abundance ratios indicate short star-formation timescales of 0.2-1.0 Gyr.

Radial mass-to-light ratio gradients cause the half-mass and half-light radii of galaxies to differ, potentially biasing studies that use half-light radii. Here we present the largest catalog to date of galaxy half-mass radii at z > 1: 7,006 galaxies in the CANDELS fields at 1.0 < z < 2.5. The sample includes both star-forming and quiescent galaxies with stellar masses 9.0 < log(M_* / M_\odot) < 11.5. We test three methods for calculating half-mass radii from multi-band PSF-matched HST imaging: two based on spatially-resolved SED modeling, and one that uses a rest-frame color profile. All three methods agree, with scatter <~0.3 dex. In agreement with previous studies, most galaxies in our sample have negative color gradients (the centers are redder than the outskirts, and r_e,mass < r_e,light). We find that color gradient strength has significant trends with increasing stellar mass, half-light radius, U-V color, and stellar mass surface density. These trends have not been seen before at z>1. Furthermore, color gradients of star-forming and quiescent galaxies show a similar redshift evolution: they are flat at z>~2, then steeply decrease as redshift decreases. This affects the galaxy mass-size relation. The normalizations of the star-forming and quiescent r_mass-M_* relations are 10-40% smaller than the corresponding r_light-M_* relations; the slopes are ~0.1-0.3 dex shallower. Finally, the half-mass radii of star-forming and quiescent galaxies at M_* = 10^{10.5}M_\odot only grow by ~1%$ and ~8% between z~2.25 and z~1.25. This is significantly less than the ~37% and ~47% size increases found when using the half-light radius.

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 present COSMOS-Drift And SHift (DASH), a Hubble Space Telescope WFC3 imaging survey of the COSMOS field in the H_160 filter. The survey comprises 456 individual WFC3 pointings corresponding to an area of 0.49 deg^2 (0.66 deg^2 when including archival data) and reaches a 5 point-source limit of H_160 =25.1 (0".3 aperture). COSMOS-DASH is the widest HST/WFC3 imaging survey in H_160 filter, tripling the extragalactic survey area in the near-infrared at HST resolution. We make the reduced H_160 mosaic available to the community. We use this dataset to measure the sizes of 162 galaxies with log(M_star/M_sun) > 11.3 at 1.5 < z < 3.0, and augment this sample with 748 galaxies at 0.1 < z < 1.5 using archival ACS imaging. We find that the median size of galaxies in this mass range changes with redshift as r_eff = (10.4+/-0.4)(1 +z)^(0.65+/-0.05) kpc. Separating the galaxies into star forming and quiescent galaxies using their restframe U-V and V-J colors, we find no statistical difference between the median sizes of the most massive star-forming and quiescent galaxies at z = 2.5: they are 4.9+/-0.9 kpc and 4.3 +/-0.3 kpc respectively. However, we do find a significant difference in the S`ersic index between the two samples, such that massive quiescent galaxies have higher central densities than star forming galaxies. We extend the size-mass analysis to lower masses by combining it with the 3D-HST/CANDELS sample of van der Wel et al. (2014), and derive empirical relations between size, mass, and redshift. Fitting a relation of the form r_eff = A m_star^a, m_star = M_star/5x10^10 M_sun and r_eff in kpc, we find log A = -0.25 log (1 + z) + 0.79 and a = -0.13 log(1 + z) + 0.27. We also provide relations for the subsamples of star forming and quiescent galaxies. Our results confirm previous studies that were based on smaller samples or ground-based imaging.

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.

The process by which massive galaxies transition from blue, star-forming disks into red, quiescent galaxies remains one of the most poorly-understood aspects of galaxy evolution. In this investigation, we attempt to gain a better understanding of how star formation is quenched by focusing on a massive post-starburst galaxy at z = 0.747. The target has a high stellar mass and a molecular gas fraction of ~30% -- unusually high for its low star formation rate. We look for indicators of star formation suppression mechanisms in the stellar kinematics and age distribution of the galaxy obtained from spatially resolved Gemini Integral-Field spectra and in the gas kinematics obtained from ALMA. We find evidence of significant rotation in the stars, but we do not detect a stellar age gradient within 5 kpc. The molecular gas is aligned with the stellar component, and we see no evidence of strong gas outflows. Our target may represent the product of a merger-induced starburst or of morphological quenching; however, our results are not completely consistent with any of the prominent quenching models.

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 high spatial resolution (FWHM$\sim$0.14'') observations of the CO($8-7$) line in GDS-14876, a compact star-forming galaxy at $z=2.3$ with total stellar mass of $\log(M_{\star}/M_{\odot})=10.9$. The spatially resolved velocity map of the inner $r\lesssim1$~kpc reveals a continous velocity gradient consistent with the kinematics of a rotating disk with $v_{\rm rot}(r=1\rm kpc)=163\pm5$ km s$^{-1}$ and $v_{\rm rot}/\sigma\sim2.5$. The gas-to-stellar ratios estimated from CO($8-7$) and the dust continuum emission span a broad range, $f^{\rm CO}_{\rm gas}=M_{\rm gas}/M_{\star}=13-45\%$ and $f^{\rm cont}_{\rm gas}=50-67\%$, but are nonetheless consistent given the uncertainties in the conversion factors. The dynamical modeling yields a dynamical mass of$\log(M_{\rm dyn}/M_{\odot})=10.58^{+0.5}_{-0.2}$ which is lower, but still consistent with the baryonic mass, $\log$(M$_{\rm bar}$= M$_{\star}$ + M$^{\rm CO}_{\rm gas}$/M$_{\odot}$)$=11.0$, if the smallest CO-based gas fraction is assumed. Despite a low, overall gas fraction, the small physical extent of the dense, star-forming gas probed by CO($8-7$), $\sim3\times$ smaller than the stellar size, implies a strong concentration that increases the gas fraction up to $f^{\rm CO, 1\rm kpc}_{\rm gas}\sim 85\%$ in the central 1 kpc. Such a gas-rich center, coupled with a high star-formation rate, SFR$\sim$ 500 M$_{\odot}$ yr$^{-1}$, suggests that GDS-14876 is quickly assembling a dense stellar component (bulge) in a strong nuclear starburst. Assuming its gas reservoir is depleted without replenishment, GDS-14876 will quickly ($t_{\rm depl}\sim27$ Myr) become a compact quiescent galaxy that could retain some fraction of the observed rotational support.

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.

The physical mechanisms that quench star formation, turning blue star-forming galaxies into red quiescent galaxies, remain unclear. In this Letter, we investigate the role of gas supply in suppressing star formation by studying the molecular gas content of post-starburst galaxies. Leveraging the wide area of the SDSS, we identify a sample of massive intermediate-redshift galaxies that have just ended their primary epoch of star formation. We present ALMA CO(2-1) observations of two of these post-starburst galaxies at z~0.7 with M* ~ 2x10^11 Msun. Their molecular gas reservoirs of (6.4 +/- 0.8) x 10^9 Msun and (34.0 +/- 1.6) x 10^9 Msun are an order of magnitude larger than comparable-mass galaxies in the local universe. Our observations suggest that quenching does not require the total removal or depletion of molecular gas, as many quenching models suggest. However, further observations are required both to determine if these apparently quiescent objects host highly obscured star formation and to investigate the intrinsic variation in the molecular gas properties of post-starburst galaxies.

Accurate measurements of galaxy masses and sizes are key to tracing galaxy evolution over time. Cosmological zoom-in simulations provide an ideal test bed for assessing the recovery of galaxy properties from observations. Here, we utilize galaxies with $M_*\sim10^{10}-10^{11.5}M_{\odot}$ at z~1.7-2 from the MassiveFIRE cosmological simulation suite, part of the Feedback in Realistic Environments (FIRE) project. Using mock multi-band images, we compare intrinsic galaxy masses and sizes to observational estimates. We find that observations accurately recover stellar masses, with a slight average underestimate of ~0.06 dex and a ~0.15 dex scatter. Recovered half-light radii agree well with intrinsic half-mass radii when averaged over all viewing angles, with a systematic offset of ~0.1 dex (with the half-light radii being larger) and a scatter of ~0.2 dex. When using color gradients to account for mass-to-light variations, recovered half-mass radii also exceed the intrinsic half-mass radii by ~0.1 dex. However, if not properly accounted for, aperture effects can bias size estimates by ~0.1 dex. No differences are found between the mass and size offsets for star-forming and quiescent galaxies. Variations in viewing angle are responsible for ~25% of the scatter in the recovered masses and sizes. Our results thus suggest that the intrinsic scatter in the mass-size relation may have previously been overestimated by ~25%. Moreover, orientation-driven scatter causes the number density of very massive galaxies to be overestimated by ~0.5 dex at $M_*\sim10^{11.5}M_{\odot}$.

We investigate the optical and Wide-field Survey Explorer (WISE) colors of "E+A" identified post-starburst galaxies, including a deep analysis on 190 post-starbursts detected in the 2{\mu}m All Sky Survey Extended Source Catalog. The post-starburst galaxies appear in both the optical green valley and the WISE Infrared Transition Zone (IRTZ). Furthermore, we find that post-starbursts occupy a distinct region [3.4]-[4.6] vs. [4.6]-[12] WISE colors, enabling the identification of this class of transitioning galaxies through the use of broad-band photometric criteria alone. We have investigated possible causes for the WISE colors of post-starbursts by constructing a composite spectral energy distribution (SED), finding that mid-infrared (4-12{\mu}m) properties of post-starbursts are consistent with either 11.3{\mu}m polycyclic aromatic hydrocarbon emission, or Thermally Pulsating Asymptotic Giant Branch (TP-AGB) and post-AGB stars. The composite SED of extended post- starburst galaxies with 22{\mu}m emission detected with signal to noise >3 requires a hot dust component to produce their observed rising mid-infrared SED between 12 and 22{\mu}m. The composite SED of WISE 22{\mu}m non-detections (S/N<3), created by stacking 22{\mu}m images, is also flat, requiring a hot dust component. The most likely source of this mid-infrared emission of these E+A galaxies is a buried active galactic nucleus. The inferred upper limit to the Eddington ratios of post-starbursts are 1e-2 to 1e-4, with an average of 1e-3. This suggests that AGNs are not radiatively dominant in these systems. This could mean that including selections able to identify active galactic nuclei as part of a search for transitioning and post-starburst galaxies would create a more complete census of the transition pathways taken as a galaxy quenches its star formation.

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.

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 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 spatially-resolved Atacama Large Millimeter/sub-millimeter Array (ALMA) 870 $\mu$m dust continuum maps of six massive, compact, dusty star-forming galaxies (SFGs) at $z\sim2.5$. These galaxies are selected for their small rest-frame optical sizes ($r_{\rm e, F160W}\sim1.6$ kpc) and high stellar-mass densities that suggest that they are direct progenitors of compact quiescent galaxies at $z\sim2$. The deep observations yield high far-infrared (FIR) luminosities of L$_{\rm IR}=10^{12.3-12.8}$ L$_{\odot}$ and star formation rates (SFRs) of SFR$=200-700$ M$_{\odot}$yr$^{-1}$, consistent with those of typical star-forming "main sequence" galaxies. The high-spatial resolution (FWHM$\sim$0.12"-0.18") ALMA and HST photometry are combined to construct deconvolved, mean radial profiles of their stellar mass and (UV+IR) SFR. We find that the dusty, nuclear IR-SFR overwhelmingly dominates the bolometric SFR up to $r\sim5$ kpc, by a factor of over 100$\times$ from the unobscured UV-SFR. Furthermore, the effective radius of the mean SFR profile ($r_{\rm e, SFR}\sim1$ kpc) is $\sim$30% smaller than that of the stellar mass profile. The implied structural evolution, if such nuclear starburst last for the estimated gas depletion time of $\Delta t=\pm100$ Myr, is a 4$\times$ increase of the stellar mass density within the central 1 kpc and a 1.6$\times$ decrease of the half-mass radius. This structural evolution fully supports dissipation-driven, formation scenarios in which strong nuclear starbursts transform larger, star-forming progenitors into compact quiescent galaxies.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

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.

We present low-resolution near-infrared stacked spectra from the 3D-HST survey up to $z=2.0$ and fit them with commonly used stellar population synthesis models: BC03 (Bruzual & Charlot, 2003), FSPS10 (Flexible Stellar Population Synthesis, Conroy & Gunn 2010), and FSPS-C3K (Conroy, Kurucz, Cargile, Castelli, in prep). The accuracy of the grism redshifts allows the unambiguous detection of many emission and absorption features, and thus a first systematic exploration of the rest-frame optical spectra of galaxies up to $z=2$. We select massive galaxies ($\rm log(M_{*} / M_{\odot}) > 10.8$), we divide them into quiescent and star-forming via a rest-frame color-color technique, and we median-stack the samples in 3 redshift bins between $z=0.5$ and $z=2.0$. We find that stellar population models fit the observations well at wavelengths below $\rm 6500 \AA$ rest-frame, but show systematic residuals at redder wavelengths. The FSPS-C3K model generally provides the best fits (evaluated with a $\chi^2_{red}$ statistics) for quiescent galaxies, while BC03 performs the best for star-forming galaxies. The stellar ages of quiescent galaxies implied by the models, assuming solar metallicity, vary from 4 Gyr at $z \sim 0.75$ to 1.5 Gyr at $z \sim 1.75$, with an uncertainty of a factor of 2 caused by the unknown metallicity. On average the stellar ages are half the age of the Universe at these redshifts. We show that the inferred evolution of ages of quiescent galaxies is in agreement with fundamental plane measurements, assuming an 8 Gyr age for local galaxies. For star-forming galaxies the inferred ages depend strongly on the stellar population model and the shape of the assumed star-formation history.

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 build a set of composite galaxy SEDs by de-redshifting and scaling multi-wavelength photometry from galaxies in the ZFOURGE survey, covering the CDFS, COSMOS, and UDS fields. From a sample of ~4000 K_s-band selected galaxies, we define 38 composite galaxy SEDs that yield continuous low-resolution spectra (R~45) over the rest-frame range 0.1-4 um. Additionally, we include far infrared photometry from the Spitzer Space Telescope and the Herschel Space Observatory to characterize the infrared properties of our diverse set of composite SEDs. From these composite SEDs we analyze the rest-frame UVJ colors, as well as the ratio of IR to UV light (IRX) and the UV slope ($\beta$) in the IRX$-\beta$ dust relation at 1<z<3. Blue star-forming composite SEDs show IRX and $\beta$ values consistent with local relations; dusty star-forming galaxies have considerable scatter, as found for local IR bright sources, but on average appear bluer than expected for their IR fluxes. We measure a tight linear relation between rest-frame UVJ colors and dust attenuation for star-forming composites, providing a direct method for estimating dust content from either (U-V) or (V-J) rest-frame colors for star-forming galaxies at intermediate redshifts.

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 the relation between galaxy structure and spectral type, using a K-selected galaxy sample at 0.5<z<2.0. Based on similarities between the UV-to-NIR spectral energy distributions, we classify galaxies into 32 spectral types. The different types span a wide range in evolutionary phases, and thus -- in combination with available CANDELS/F160W imaging -- are ideal to study the structural evolution of galaxies. Effective radii (R_e) and Sersic parameters (n) have been measured for 572 individual galaxies, and for each type, we determine R_e at fixed stellar mass by correcting for the mass-size relation. We use the rest-frame U-V vs. V-J diagram to investigate evolutionary trends. When moving into the direction perpendicular to the star-forming sequence, in which we see the Halpha equivalent width and the specific star formation rate (sSFR) decrease, we find a decrease in R_e and an increase in n. On the quiescent sequence we find an opposite trend, with older redder galaxies being larger. When splitting the sample into redshift bins, we find that young post-starburst galaxies are most prevalent at z>1.5 and significantly smaller than all other galaxy types at the same redshift. This result suggests that the suppression of star formation may be associated with significant structural evolution at z>1.5. At z<1, galaxy types with intermediate sSFRs (10^{-11.5}-10^{-10.5} yr^-1) do not have post-starburst SED shapes. These galaxies have similar sizes as older quiescent galaxies, implying that they can passively evolve onto the quiescent sequence, without increasing the average size of the quiescent galaxy population.

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 present reduced data and data products from the 3D-HST survey, a 248-orbit HST Treasury program. The survey obtained WFC3 G141 grism spectroscopy in four of the five CANDELS fields: AEGIS, COSMOS, GOODS-S, and UDS, along with WFC3 $H_{140}$ imaging, parallel ACS G800L spectroscopy, and parallel $I_{814}$ imaging. In a previous paper (Skelton et al. 2014) we presented photometric catalogs in these four fields and in GOODS-N, the fifth CANDELS field. Here we describe and present the WFC3 G141 spectroscopic data, again augmented with data from GO-1600 in GOODS-N. The data analysis is complicated by the fact that no slits are used: all objects in the WFC3 field are dispersed, and many spectra overlap. We developed software to automatically and optimally extract interlaced 2D and 1D spectra for all objects in the Skelton et al. (2014) photometric catalogs. The 2D spectra and the multi-band photometry were fit simultaneously to determine redshifts and emission line strengths, taking the morphology of the galaxies explicitly into account. The resulting catalog has 98,663 measured redshifts and line strengths down to $JH_{IR}\leq 26$ and 22,548 with $JH_{IR}\leq 24$, where we comfortably detect continuum emission. Of this sample 5,459 galaxies are at $z>1.5$ and 9,621 are at $0.7<z<1.5$, where H$\alpha$ falls in the G141 wavelength coverage. Based on comparisons with ground-based spectroscopic redshifts, and on analyses of paired galaxies and repeat observations, the typical redshift error for $JH_{IR}\leq 24$ galaxies in our catalog is $\sigma_z \approx 0.003 \times (1+z)$, i.e., one native WFC3 pixel. The $3\sigma$ limit for emission line fluxes of point sources is $1.5\times10^{-17}$ ergs s$^{-1}$ cm$^{-2}$. We show various representations of the full dataset, as well as individual examples that highlight the range of spectra that we find in the survey.

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

In this paper we study a key phase in the formation of massive galaxies: the transition of star forming galaxies into massive (M_stars~10^11 Msun), compact (r_e~1 kpc) quiescent galaxies, which takes place from z~3 to z~1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star forming galaxies exists at z~2, using K-band spectroscopy of 25 of these objects at 2.0<z<2.5. They have a median NII/Halpha ratio of 0.6, are highly obscured with SFR(tot)/SFR(Halpha)~10, and have a large range of observed line widths. We infer from the kinematics and spatial distribution of Halpha that the galaxies have rotating disks of ionized gas that are a factor of ~2 more extended than the stellar distribution. By combining measurements of individual galaxies, we find that the kinematics are consistent with a nearly Keplerian fall-off from V_rot~500 km/s at 1 kpc to V_rot~250 km/s at 7 kpc, and that the total mass out to this radius is dominated by the dense stellar component. Next, we study the size and mass evolution of the progenitors of compact massive galaxies. Even though individual galaxies may have had complex histories with periods of compaction and mergers, we show that the population of progenitors likely followed a simple inside-out growth track in the size-mass plane of d(log r_e) ~ 0.3 d(log M_stars). This mode of growth gradually increases the stellar mass within a fixed physical radius, and galaxies quench when they reach a stellar density or velocity dispersion threshold. As shown in other studies, the mode of growth changes after quenching, as dry mergers take the galaxies on a relatively steep track in the size-mass plane.

It is well established that (1) star-forming galaxies follow a relation between their star formation rate (SFR) and stellar mass (M$_{\star}$), the "star-formation sequence", and (2) the SFRs of galaxies correlate with their structure, where star-forming galaxies are less concentrated than quiescent galaxies at fixed mass. Here, we consider whether the scatter and slope of the star-formation sequence is correlated with systematic variations in the Sersic indices, $n$, of galaxies across the SFR-M$_{\star}$ plane. We use a mass-complete sample of 23,848 galaxies at $0.5<z<2.5$ selected from the 3D-HST photometric catalogs. Galaxy light profiles parameterized by $n$ are based on Hubble Space Telescope CANDELS near-infrared imaging. We use a single SFR indicator empirically-calibrated from stacks of Spitzer/MIPS 24$\mu$m imaging, adding the unobscured and obscured star formation. We find that the scatter of the star-formation sequence is related in part to galaxy structure; the scatter due to variations in $n$ at fixed mass for star-forming galaxies ranges from 0.14$\pm$0.02 dex at $z\sim2$ to 0.30$\pm$0.04 dex at $z<1$. While the slope of the log(SFR)-log(M$_{\star}$) relation is of order unity for disk-like galaxies, galaxies with $n>2$ (implying more dominant bulges) have significantly lower SFR/M$_{\star}$ than the main ridgeline of the star-formation sequence. These results suggest that bulges in massive $z\sim2$ galaxies are actively building up, where the stars in the central concentration are relatively young. At $z<1$, the presence of older bulges within star-forming galaxies lowers global SFR/M$_{\star}$, decreasing the slope and contributing significantly to the scatter of the star-formation sequence.

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

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]

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

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.

We present results on the z~2.3 mass-metallicity relation (MZR) using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We use an initial sample of 87 star-forming galaxies with spectroscopic coverage of H\beta, [OIII]\lambda 5007, H\alpha, and [NII]\lambda 6584 rest-frame optical emission lines, and estimate the gas-phase oxygen abundance based on the N2 and O3N2 strong-line indicators. We find a positive correlation between stellar mass and metallicity among individual z~2.3 galaxies using both the N2 and O3N2 indicators. We also measure the emission-line ratios and corresponding oxygen abundances for composite spectra in bins of stellar mass. Among composite spectra, we find a monotonic increase in metallicity with increasing stellar mass, offset ~0.15-0.3 dex below the local MZR. When the sample is divided at the median star-formation rate (SFR), we do not observe significant SFR dependence of the z~2.3 MZR among either individual galaxies or composite spectra. We furthermore find that z~2.3 galaxies have metallicities ~0.1 dex lower at a given stellar mass and SFR than is observed locally. This offset suggests that high-redshift galaxies do not fall on the local "fundamental metallicity relation" among stellar mass, metallicity, and SFR, and may provide evidence of a phase of galaxy growth in which the gas reservoir is built up due to inflow rates that are higher than star-formation and outflow rates. However, robust conclusions regarding the gas-phase oxygen abundances of high-redshift galaxies await a systematic reappraisal of the application of locally calibrated metallicity indicators at high redshift.

We explore the relation between the dynamical mass-to-light ratio ($M/L$) and rest-frame color of massive quiescent galaxies out to z~2. We use a galaxy sample with measured stellar velocity dispersions in combination with Hubble Space Telescope and ground-based multi-band photometry. Our sample spans a large range in $\log M_{dyn}/L_{g}$ (of 1.6~dex) and $\log M_{dyn}/L_{K}$ (of 1.3~dex). There is a strong, approximately linear correlation between the $M/L$ for different wavebands and rest-frame color. The root-mean-scatter scatter in $\log~M_{dyn}/L$ residuals implies that it is possible to estimate the $M/L$ with an accuracy of ~0.25 dex from a single rest-frame optical color. Stellar population synthesis (SPS) models with a Salpeter stellar initial mass function (IMF) can not simultaneously match $M_{dyn}/L_{g}$ vs. $(g-z)_{rest-frame}$ and $M_{dyn}/L_{K}$ vs. $(g-K)_{rest-frame}$. By changing the slope of the IMF we are still unable to explain the M/L of the bluest and reddest galaxies. We find that an IMF with a slope between $\alpha=2.35$ and $\alpha=1.35$ provides the best match. We also explore a broken IMF with a Salpeter slope at $M<1M_{\odot}$ and $M>4M_{\odot}$ and a slope $\alpha$ in the intermediate region. The data favor a slope of $\alpha=1.35$ over $\alpha=2.35$. Nonetheless, our results show that variations between different SPS models are comparable to the IMF variations. In our analysis we assume that the variation in $M/L$ and color is driven by differences in age, and that other contributions (e.g., metallicity evolution, dark matter) are small. These assumptions may be an important source of uncertainty as galaxies evolve in more complex ways.

We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift z = 0 galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.

We investigate the star formation rates of quiescent galaxies at high redshift (0.3 < z < 2.5) using 3D-HST WFC3 grism spectroscopy and Spitzer mid-infrared data. We select quiescent galaxies on the basis of the widely used UVJ color-color criteria. Spectral energy distribution fitting (rest frame optical and near-IR) indicates very low star formation rates for quiescent galaxies (sSFR ~ 10^-12 yr^-1). However, SED fitting can miss star formation if it is hidden behind high dust obscuration and ionizing radiation is re-emitted in the mid-infrared. It is therefore fundamental to measure the dust-obscured SFRs with a mid-IR indicator. We stack the MIPS-24um images of quiescent objects in five redshift bins centered on z = 0.5, 0.9, 1.2, 1.7, 2.2 and perform aperture photometry. Including direct 24um detections, we find sSFR ~ 10^-11.9 * (1+z)^4 yr^-1. These values are higher than those indicated by SED fitting, but at each redshift they are 20-40 times lower than those of typical star forming galaxies. The true SFRs of quiescent galaxies might be even lower, as we show that the mid-IR fluxes can be due to processes unrelated to ongoing star formation, such as cirrus dust heated by old stellar populations and circumstellar dust. Our measurements show that star formation quenching is very efficient at every redshift. The measured SFR values are at z > 1.5 marginally consistent with the ones expected from gas recycling (assuming that mass loss from evolved stars refuels star formation) and well above that at lower redshifts.

The Fundamental Plane (FP) of early-type galaxies, relating the effective radius, velocity dispersion, and surface brightness, has long been recognized as a unique tool for analyzing galaxy structure and evolution. With the discovery of distant quiescent galaxies and the introduction of high sensitivity near-infrared spectrographs, it is now possible to explore the FP out to z~2. In this Letter we study the evolution of the FP out to z~2 using kinematic measurements of massive quiescent galaxies ($M_{*}>10^{11} M_{\odot}$). We find preliminary evidence for the existence of an FP out to z~2. The scatter of the FP, however, increases from z~0 to z~2, even when taking into account the larger measurement uncertainties at higher redshifts. We find a strong evolution of the zero point from z~2 to z~0: $\Delta\log_{10}M/L_g\propto(-0.49\pm0.03)~z$. In order to assess whether our spectroscopic sample is representative of the early-type galaxy population at all redshifts, we compare their rest-frame g-z colors with those from a larger mass complete sample of quiescent galaxies. At z>1 we find that the spectroscopic sample is bluer. We use the color offsets to estimate a mass-to-light ratio (M/L) correction. The implied FP zero point evolution after correction is significantly smaller: $\Delta\log_{10}M/L_g\propto(-0.39\pm 0.02)~z$. This is consistent with an apparent formation redshift of $z_{\rm{form}}=6.62^{+3.19}_{-1.44}$ for the underlying population, ignoring the effects of progenitor bias. A more complete spectroscopic sample is required at z~2 to properly measure the M/L evolution from the FP evolution.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

The dense interiors of massive galaxies are among the most intriguing environments in the Universe. In this paper we ask when these dense cores were formed and determine how galaxies gradually assembled around them. We select galaxies that have a stellar mass >3x10^10 Msun inside r=1 kpc out to z=2.5, using the 3D-HST survey and data at low redshift. Remarkably, the number density of galaxies with dense cores appears to have decreased from z=2.5 to the present. This decrease is probably mostly due to stellar mass loss and the resulting adiabatic expansion, with some contribution from merging. We infer that dense cores were mostly formed at z>2.5, consistent with their largely quiescent stellar populations. While the cores appear to form early, the galaxies in which they reside show strong evolution: their total masses increase by a factor of 2-3 from z=2.5 to z=0 and their effective radii increase by a factor of 5-6. As a result, the contribution of dense cores to the total mass of the galaxies in which they reside decreases from ~50% at z=2.5 to ~15% at z=0. Because of their early formation, the contribution of dense cores to the total stellar mass budget of the Universe is a strong function of redshift. The stars in cores with M_1kpc>3x10^10 Msun make up ~0.1% of the stellar mass density of the Universe today but 10%-20% at z~2, depending on their IMF. The formation of these cores required the conversion of ~10^11 Msun of gas into stars within ~1 kpc, while preventing significant star formation at larger radii.

The 3D-HST and CANDELS programs have provided WFC3 and ACS spectroscopy and photometry over ~900 square arcminutes in five fields: AEGIS, COSMOS, GOODS-North, GOODS-South, and the UKIDSS UDS field. All these fields have a wealth of publicly available imaging datasets in addition to the HST data, which makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range. In this paper we describe a photometric analysis of the CANDELS and 3D-HST HST imaging and the ancillary imaging data at wavelengths 0.3um to 8um. Objects were selected in the WFC3 near-IR bands, and their SEDs were determined by carefully taking the effects of the point spread function in each observation into account. A total of 147 distinct imaging datasets were used in the analysis. The photometry is made available in the form of six catalogs: one for each field, as well as a master catalog containing all objects in the entire survey. We also provide derived data products: photometric redshifts, determined with the EAZY code, and stellar population parameters determined with the FAST code. We make all the imaging data that were used in the analysis available, including our reductions of the WFC3 imaging in all five fields. 3D-HST is a spectroscopic survey with the WFC3 and ACS grisms, and the photometric catalogs presented here constitute a necessary first step in the analysis of these grism data. All the data presented in this paper are available through the 3D-HST website.

The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct comparison between dust around star-forming regions ($A_{V,\mathrm{HII}}$) and the integrated dust content ($A_{V,\mathrm{star}}$). We select a sample of 163 galaxies between $1.36\le{}z\le1.5$ with H$\alpha$ signal-to-noise ratio $\ge5$ and measure Balmer decrements from stacked spectra to calculate $A_{V,\mathrm{HII}}$. First, we stack spectra in bins of $A_{V,\mathrm{star}}$, and find that $A_{V,\mathrm{HII}}=1.86\,A_{V,\mathrm{star}}$, with a significance of $\sigma=1.7$. Our result is consistent with the two-component dust model, in which galaxies contain both diffuse and stellar birth cloud dust. Next, we stack spectra in bins of specific star formation rate ($\log\,\mathrm{SSFR}$), star formation rate ($\log\,\mathrm{SFR}$), and stellar mass ($\log{}M_*$). We find that on average $A_{V,\mathrm{HII}}$ increases with SFR and mass, but decreases with increasing SSFR. Interestingly, the data hint that the amount of extra attenuation decreases with increasing SSFR. This trend is expected from the two-component model, as the extra attenuation will increase once older stars outside the star-forming regions become more dominant in the galaxy spectrum. Finally, using Balmer decrements we derive dust-corrected H$\alpha$ SFRs, and find that stellar population modeling produces incorrect SFRs if rapidly declining star formation histories are included in the explored parameter space.

We have used near-ultraviolet (NUV) to mid-infrared (MIR) composite spectral energy distributions (SEDs) to simultaneously model the attenuated stellar and dust emission of 0.5 < z < 2.0 galaxies. These composite SEDs were previously constructed from the photometric catalogs of the NEWFIRM Medium-Band Survey, by stacking the observed photometry of galaxies that have similar rest-frame NUV-to-NIR SEDs. In this work, we include a stacked MIPS 24 micron measurement for each SED type to extend the SEDs to rest-frame MIR wavelengths. Consistent with previous studies, the observed MIR emission for most SED types is higher than expected from only the attenuated stellar emission. We fit the NUV-to-MIR composite SEDs by the Flexible Stellar Population Synthesis (SPS) models, which include both stellar and dust emission. We compare the best-fit star formation rates (SFRs) to the SFRs based on simple UV+IR estimators. Interestingly, the UV and IR luminosities overestimate SFRs - compared to the model SFRs - by more than ~ 1 dex for quiescent galaxies, while for the highest star-forming galaxies in our sample the two SFRs are broadly consistent. The difference in specific SFRs also shows a gradually increasing trend with declining specific SFR, implying that quiescent galaxies have even lower specific SFRs than previously found. Contributions from evolved stellar populations to both the UV and the MIR SEDs most likely explain the discrepancy. Based on this work, we conclude that SFRs should be determined from modeling the attenuated stellar and dust emission simultaneously, instead of employing simple UV+IR-based SFR estimators.

We examine how the total X-ray luminosity correlates with stellar mass, stellar population, and redshift for a K-band limited sample of ~3500 galaxies at 0.5<z<2.0 from the NEWFIRM Medium Band Survey in the COSMOS field. The galaxy sample is divided into 32 different galaxy types, based on similarities between the spectral energy distributions. For each galaxy type, we further divide the sample into bins of redshift and stellar mass, and perform an X-ray stacking analysis using the Chandra COSMOS (C-COSMOS) data. We find that full band X-ray luminosity is primarily increasing with stellar mass, and at similar mass and spectral type is higher at larger redshifts. When comparing at the same stellar mass, we find that the X-ray luminosity is slightly higher for younger galaxies (i.e., weaker 4000\AA breaks), but the scatter in this relation is large. We compare the observed X-ray luminosities to those expected from low and high mass X-ray binaries (XRBs). For blue galaxies, XRBs can almost fully account for the observed emission, while for older galaxies with larger 4000\AA breaks, active galactic nuclei (AGN) or hot gas dominate the measured X-ray flux. After correcting for XRBs, the X-ray luminosity is still slightly higher in younger galaxies, although this correlation is not significant. AGN appear to be a larger component of galaxy X-ray luminosity at earlier times, as the hardness ratio increases with redshift. Together with the slight increase in X-ray luminosity this may indicate more obscured AGNs or higher accretion rates at earlier times.

We present Keck/MOSFIRE K-band spectroscopy of the first mass-selected sample of galaxies at $z\sim2.3$. Targets are selected from the 3D-HST Treasury survey. The six detected galaxies have a mean [NII]$\lambda$6584/H$\alpha$ ratio of $0.27\pm0.01$, with a small standard deviation of 0.05. This mean value is similar to that of UV-selected galaxies of the same mass. The mean gas-phase oxygen abundance inferred from the [NII]/H$\alpha$ ratios depends on the calibration method, and ranges from 12+log(O/H)$_{gas}=8.57$ for the {Pettini} & {Pagel} (2004) calibration to 12+log(O/H)$_{gas}= 8.87$ for the {Maiolino} {et~al.} (2008) calibration. Measurements of the stellar oxygen abundance in nearby quiescent galaxies with the same number density indicate 12+log(O/H)$_{stars}= 8.95$, similar to the gas-phase abundances of the $z\sim2.3$ galaxies if the {Maiolino} {et~al.} (2008) calibration is used. This suggests that these high-redshift star forming galaxies may be progenitors of today's massive early-type galaxies. The main uncertainties are the absolute calibration of the gas-phase oxygen abundance and the incompleteness of the $z\sim2.3$ sample: the galaxies with detected H$\alpha$ tend to be larger and have higher star formation rates than the galaxies without detected H$\alpha$, and we may still be missing the most dust-obscured progenitors.

The Fundamental Plane (FP) is an empirical relation between the size, surface brightness, and velocity dispersion of early-type galaxies. This relation has been studied extensively for early-type galaxies in the local universe to constrain galaxy formation mechanisms. The evolution of the zeropoint of this plane has been extended to high redshifts to study the luminosity evolution of massive galaxies, under the assumption of structural homology. In this work, we assess this assumption by replacing surface brightness with stellar mass density and present the evolution of the "mass FP" for massive, quiescent galaxies since z~2. By accounting for stellar populations, we thereby isolate and trace structural and dynamical evolution. Despite the observed dramatic evolution in the sizes and morphologies of massive galaxies since z~3, we find that quiescent galaxies lie on the mass FP out to z~2. In contrast with ~1.4 dex evolution in the luminosity FP, average residuals from the z~0 mass FP are less than ~0.15 dex since z~2. Assuming the Hyde & Bernardi (2009) mass FP slope, we find that this minimal offset scales as (1+z)^{-0.095+/-0.043}. This result lends credence to previous studies that derived luminosity evolution from the FP. Therefore, despite their compact sizes and suggestions that massive galaxies are more disk-like at z~2, the relationship between their dynamics and structural properties are consistent with local early-type galaxies. Finally, we find no strong evidence for a tilt of the mass FP relative to the Virial plane, but emphasize the need for full models including selection biases to fully investigate this issue.

This letter utilizes composite spectral energy distributions (SEDs) constructed from NEWFIRM Medium-Band Survey photometry to constrain the dust attenuation curve in 0.5<z<2.0 galaxies. Based on similarities between the full SED shapes (0.3-8 micron), we have divided galaxies in 32 different spectral classes and stacked their photometry. As each class contains galaxies over a range in redshift, the resulting rest-frame SEDs are well-sampled in wavelength and show various spectral features including Halpha and the UV dust bump at 2175 Angstrom. We fit all composite SEDs with flexible stellar population synthesis models, while exploring attenuation curves with varying slopes and UV bump strengths. The Milky Way and Calzetti law provide poor fits at UV wavelengths for nearly all SEDs. Consistent with previous studies, we find that the best-fit attenuation law varies with spectral type. There is a strong correlation between the best-fit dust slope and UV bump strength, with steeper laws having stronger bumps. Moreover, the attenuation curve correlates with specific star formation rate (SFR), with more active galaxies having shallower dust curves and weaker bumps. There is also a weak correlation with inclination. The observed trends can be explained by differences in the dust-to-star geometry, a varying grain size distribution, or a combination of both. Our results have several implications for galaxy evolution studies. First, the assumption of a universal dust model leads to biases in derived galaxy properties. Second, the presence of a dust bump may result in underestimated values for the UV slope, used to correct SFRs of distant galaxies.

Galaxies with the mass of the Milky Way dominate the stellar mass density of the Universe but it is uncertain how and when they were assembled. Here we study progenitors of these galaxies out to z=2.5, using data from the 3D-HST and CANDELS Treasury surveys. We find that galaxies with present-day stellar masses of log(M)~10.7 built ~90% of their stellar mass since z=2.5, with most of the star formation occurring before z=1. In marked contrast to the assembly history of massive elliptical galaxies, mass growth is not limited to large radii: the mass in the central 2 kpc of the galaxies increased by a factor of 3.2+-0.8 between z=2.5 and z=1. We therefore rule out simple models in which bulges were fully assembled at high redshift and disks gradually formed around them. Instead, bulges (and black holes) likely formed in lockstep with disks, through bar instabilities, migration, or other processes. We find that after z=1 the growth in the central regions gradually stopped and the disk continued to build, consistent with recent studies of the gas distributions in z~1 galaxies and the properties of many spiral galaxies today.

Using stellar kinematics measurements, we investigate the growth of massive, quiescent galaxies from z~2 to today. We present X-Shooter spectra from the UV to NIR and dynamical mass measurements of 5 quiescent massive (>10^11 Msun) galaxies at z~2. This triples the sample of z>1.5 galaxies with well constrained (dsigma <100 km/s) velocity dispersion measurements. From spectral population synthesis modeling we find that these galaxies have stellar ages that range from 0.5-2 Gyr, with no signs of ongoing star formation. We measure velocity dispersions (290-450 km/s) from absorption lines and find that they are 1.6-2.1 times higher than those of galaxies in the SDSS at the same mass. Sizes are measured using GALFIT from HST-WFC3 H160 and UDS K-band images. The dynamical masses correspond well to the SED based stellar masses, with dynamical masses that are ~15% higher. We find that M_*/M_dyn may decrease slightly with time, which could reflect the increase of the dark matter fraction within an increasing effective radius. We combine different stellar kinematic studies from the literature, and examine the structural evolution from z~2 to z~0: we confirm that at fixed dynamical mass, the effective radius increases by a factor of ~2.8, and the velocity dispersion decreases by a factor of ~1.7. The mass density within one effective radius decreases by a factor of ~20, while within a fixed physical radius (1 kpc) it decreases only mildly (~2). When we allow for an evolving mass limit by selecting a population of galaxies at fixed number density, a stronger size growth with time is found (~4), velocity dispersion decreases by a factor of ~1.4, and interestingly, the mass density within 1 kpc is consistent with no evolution. This finding suggests that massive quiescent galaxies at ~2 grow inside-out, consistent with the expectations from minor mergers.

Quiescent galaxies at z~2 have been identified in large numbers based on rest-frame colors, but only a small number of these galaxies have been spectroscopically confirmed to show that their rest-frame optical spectra show either strong Balmer or metal absorption lines. Here, we median stack the rest-frame optical spectra for 171 photometrically-quiescent galaxies at 1.4 < z < 2.2 from the 3D-HST grism survey. In addition to Hbeta (4861A), we unambiguously identify metal absorption lines in the stacked spectrum, including the G-band (4304A), Mg I (5175A), and Na I (5894A). This finding demonstrates that galaxies with relatively old stellar populations already existed when the universe was ~3 Gyr old, and that rest-frame color selection techniques can efficiently select them. We find an average age of 1.3^0.1_0.3 Gyr when fitting a simple stellar population to the entire stack. We confirm our previous result from medium-band photometry that the stellar age varies with the colors of quiescent galaxies: the reddest 80% of galaxies are dominated by metal lines and have a relatively old mean age of 1.6^0.5_0.4 Gyr, whereas the bluest (and brightest) galaxies have strong Balmer lines and a spectroscopic age of 0.9^0.2_0.1 Gyr. Although the spectrum is dominated by an evolved stellar population, we also find [OIII] and Hbeta emission. Interestingly, this emission is more centrally concentrated than the continuum with L_[OIII] = 1.7 +/- 0.3 x 10^40 erg s^-1, indicating residual central star formation or nuclear activity.

The assembly of galaxies can be described by the distribution of their star formation as a function of cosmic time. Thanks to the WFC3 grism on HST it is now possible to measure this beyond the local Universe. Here we present the spatial distribution of Halpha emission for a sample of 54 strongly star-forming galaxies at z~1 in the 3D-HST Treasury survey. By stacking the Halpha emission we find that star formation occurred in approximately exponential distributions at z~1, with median Sersic index of n=1.0+-0.2. The stacks are elongated with median axis ratios of b/a=0.58+-0.09 in Halpha, consistent with (possibly thick) disks at random orientation angles. Keck spectra obtained for a subset of eight of the galaxies show clear evidence for rotation, with inclination-corrected velocities of 90 to 330 km/s. The most straightforward interpretation of our results is that star formation in strongly star-forming galaxies at z~1 generally occurred in disks. The disks appear to be "scaled-up" versions of nearby spiral galaxies: they have EW(Halpha)~100 Angstroms out to the solar orbit and they have star formation surface densities above the threshold for driving galactic scale winds.

We present deep Keck/LRIS spectroscopy and HST/WFC3 imaging in the rest-frame optical for a sample of eight galaxies at z~1.5 with high photometrically-determined stellar masses. The data are combined with VLT/XShooter spectra of five galaxies from van de Sande et al. (2011, 2012 to be submitted). We find that these thirteen galaxies have high velocity dispersions, with a median of sigma=301 km s^{-1}. This high value is consistent with their relatively high stellar masses and compact sizes. We study their stellar populations using the strength of Balmer absorption lines, which are not sensitive to dust absorption. We find a large range in Balmer absorption strength, with many galaxies showing very strong lines indicating young ages. The median Hdelta_A equivalent width, determined directly or inferred from the H10 line, is 5.4 Angstroms, indicating a luminosity-weighted age of ~1 Gyr. Although this value may be biased towards higher values because of selection effects,high-dispersion galaxies with such young ages are extremely rare in the local Universe. Interestingly we do not find a simple correlation with rest-frame U-V color: some of the reddest galaxies have very strong Balmer absorption lines. These results demonstrate that many high-dispersion galaxies at z~1.5 were quenched recently. This implies that there must be a population of star-forming progenitors at z~2 with high velocity dispersions or linewidths, which are notoriously absent from CO/Halpha selected surveys.

We investigate the evolution of the Halpha equivalent width, EW(Halpha), with redshift and its dependence on stellar mass, taking advantage of the first data from the 3D-HST survey, a large spectroscopic Treasury program with the Hubble Space Telescope WFC3. Combining our Halpha measurements of 854 galaxies at 0.8<z<1.5 with those of ground based surveys at lower and higher redshift, we can consistently determine the evolution of the EW(Halpha) distribution from z=0 to z=2.2. We find that at all masses the characteristic EW(Halpha) is decreasing towards the present epoch, and that at each redshift the EW(Halpha) is lower for high-mass galaxies. We measure a slope of EW(Halpha) ~ (1+z)^(1.8) with little mass dependence. Qualitatively, this measurement is a model-independent confirmation of the evolution of star forming galaxies with redshift. A quantitative conversion of EW(Halpha) to sSFR is very model dependent, because of differential reddening corrections between the continuum SED and the Balmer lines. The observed EW(Halpha) can be reproduced with a simple model in which the SFR for galaxies rises to the epoch of z~2.5 and then decreases with time to z = 0. The model implies that the EW(Halpha) rises to 400 A at z=8. The sSFR evolves faster than EW(Halpha), as the mass-to-light ratio also evolves with redshift. In this context, we find that the sSFR evolves as (1+z)^(3.2), nearly independent of mass, consistent with previous reddening insensitive estimates. We confirm previous results that the observed slope of the sSFR-z relation is steeper than the one predicted by models, but models and observations agree in finding little mass dependence.

We present evidence of large-scale outflows from three low-mass (log(M/M_sun)~9.75) star-forming (SFR >4 M_sun/yr) galaxies observed at z=1.24, z=1.35 and z=1.75 in the 3D-HST Survey. Each of these galaxies is located within a projected physical distance of 60 kpc around the sight line to the quasar SDSS J123622.93+621526.6, which exhibits well-separated strong (W_r>0.8A) Mg II absorption systems matching precisely to the redshifts of the three galaxies. We derive the star formation surface densities from the H-alpha emission in the WFC3 G141 grism observations for the galaxies and find that in each case the star formation surface density well-exceeds 0.1 M_sun/yr/kpc^2, the typical threshold for starburst galaxies in the local Universe. From a small but complete parallel census of the 0.65<z<2.6 galaxies with H_140<24 proximate to the quasar sight line, we detect Mg II absorption associated with galaxies extending to physical distances of 130 kpc. We determine that the W_r>0.8A Mg II covering fraction of star-forming galaxies at 1<z<2 may be as large as unity on scales extending to at least 60 kpc, providing early constraints on the typical extent of starburst-driven winds around galaxies at this redshift. Our observations additionally suggest that the azimuthal distribution of W_r>0.4A Mg II absorbing gas around star-forming galaxies may evolve from z~2 to the present, consistent with recent observations of an increasing collimation of star formation-driven outflows with time from z~3.

We present 3D-HST, a near-infrared spectroscopic Treasury program with the Hubble Space Telescope for studying the processes that shape galaxies in the distant Universe. 3D-HST provides rest-frame optical spectra for a sample of ~7000 galaxies at 1<z<3.5, the epoch when 60% of all star formation took place, the number density of quasars peaked, the first galaxies stopped forming stars, and the structural regularity that we see in galaxies today must have emerged. 3D-HST will cover 3/4 (625 sq.arcmin) of the CANDELS survey area with two orbits of primary WFC3/G141 grism coverage and two to four parallel orbits with the ACS/G800L grism. In the IR these exposure times yield a continuum signal-to-noise of ~5 per resolution element at H~23.1 and a 5sigma emission line sensitivity of 5x10-17 erg/s/cm2 for typical objects, improving by a factor of ~2 for compact sources in images with low sky background levels. The WFC3/G141 spectra provide continuous wavelength coverage from 1.1-1.6 um at a spatial resolution of ~0."13, which, combined with their depth, makes them a unique resource for studying galaxy evolution. We present the preliminary reduction and analysis of the grism observations, including emission line and redshift measurements from combined fits to the extracted grism spectra and photometry from ancillary multi-wavelength catalogs. The present analysis yields redshift estimates with a precision of sigma(z)=0.0034(1+z), or sigma(v)~1000 km/s. We illustrate how the generalized nature of the survey yields near-infrared spectra of remarkable quality for many different types of objects, including a quasar at z=4.7, quiescent galaxies at z~2, and the most distant T-type brown dwarf star known. The CANDELS and 3D-HST surveys combined will provide the definitive imaging and spectroscopic dataset for studies of the 1<z<3.5 Universe until the launch of the James Webb Space Telescope.

We investigate the build-up of galaxies at z~1 using maps of Halpha and stellar continuum emission for a sample of 57 galaxies with rest-frame Halpha equivalent widths >100 Angstroms in the 3D-HST grism survey. We find that the Halpha emission broadly follows the rest-frame R-band light but that it is typically somewhat more extended and clumpy. We quantify the spatial distribution with the half-light radius. The median Halpha effective radius r_e(Halpha) is 4.2+-0.1 kpc but the sizes span a large range, from compact objects with r_e(Halpha) ~ 1.0 kpc to extended disks with r_e(Halpha) ~ 15 kpc. Comparing Halpha sizes to continuum sizes, we find <r_e(Halpha)/r_e(R)>=1.3+-0.1 for the full sample. That is, star formation, as traced by Halpha, typically occurs out to larger radii than the rest-frame R-band stellar continuum; galaxies are growing their radii and building up from the inside out. This effect appears to be somewhat more pronounced for the largest galaxies. Using the measured Halpha sizes, we derive star formation rate surface densities. We find that they range from ~0.05 Msun yr^{-1} kpc^{-2} for the largest galaxies to ~5 Msun yr^{-1} kpc^{-2} for the smallest galaxies, implying a large range in physical conditions in rapidly star-forming z~1 galaxies. Finally, we infer that all galaxies in the sample have very high gas mass fractions and stellar mass doubling times < 500 Myr. Although other explanations are also possible, a straightforward interpretation is that we are simultaneously witnessing the rapid formation of compact bulges and large disks at z~1.

We study the growth of the red sequence through the number density and structural evolution of a sample of young and old quiescent galaxies at 0<z<2. The galaxies are selected from the NEWFIRM Medium-Band Survey (NMBS) in the Cosmic Evolution Survey (COSMOS) field. We find a large population of massive young recently quenched ("post-starburst") galaxies at z>1 that are almost non-existent at z<1; their number density is 5 x 10^{-5} Mpc^{-3} at z=2, whereas it is a factor of 10 less at z=0.5. The observed number densities of young and old quiescent galaxies at z>1 are consistent with a simple model in which all old quiescent galaxies were once identified as post-starburst galaxies. We find that the overall population of quiescent galaxies have smaller sizes and slightly more elongated shapes at higher redshift, in agreement with other recent studies. Interestingly, the most recently quenched galaxies at 1<z<2 are not larger, and possibly even smaller, than older galaxies at those redshifts. This result is inconsistent with the idea that the evolution of the average size of quiescent galaxies is largely driven by continuous transformations of larger, star-forming galaxies: in that case, the youngest quiescent galaxies would also be the largest. Instead, mergers or other mechanisms appear to be required to explain the size growth of quiescent galaxies from z=2 to the present.

We present first results from the 3D-HST program, a near-IR spectroscopic survey performed with the Wide Field Camera 3 on the Hubble Space Telescope. We have used 3D-HST spectra to measure redshifts and Halpha equivalent widths for a stellar mass-limited sample of 34 galaxies at 1<z<1.5 with M(stellar)>10^11 M(sun) in the COSMOS, GOODS, and AEGIS fields. We find that a substantial fraction of massive galaxies at this epoch are forming stars at a high rate: the fraction of galaxies with Halpha equivalent widths >10 A is 59%, compared to 10% among SDSS galaxies of similar masses at z=0.1. Galaxies with weak Halpha emission show absorption lines typical of 2-4 Gyr old stellar populations. The structural parameters of the galaxies, derived from the associated WFC3 F140W imaging data, correlate with the presence of Halpha: quiescent galaxies are compact with high Sersic index and high inferred velocity dispersion, whereas star-forming galaxies are typically large two-armed spiral galaxies, with low Sersic index. Some of these star forming galaxies might be progenitors of the most massive S0 and Sa galaxies. Our results challenge the idea that galaxies at fixed mass form a homogeneous population with small scatter in their properties. Instead we find that massive galaxies form a highly diverse population at z>1, in marked contrast to the local Universe.

We measure spectral features of ~3500 K-selected galaxies at 0.5<z<2.0 from high quality medium-band photometry using a new technique. First, we divide the galaxy sample in 32 subsamples based on the similarities between the full spectral energy distributions (SEDs) of the galaxies. For each of these 32 galaxy types we construct a composite SED by de-redshifting and scaling the observed photometry. This approach increases the signal-to-noise ratio and sampling of galaxy SEDs and allows for model-independent stellar population studies. The composite SEDs are of spectroscopic quality, and facilitate -- for the first time -- Halpha measurement for a large magnitude-limited sample of distant galaxies. The linewidths indicate a photometric redshift uncertainty of dz<0.02x(1+z). The composite SEDs also show the Balmer and 4000 Angstrom breaks, MgII absorption at ~2800 Angstrom, the dust absorption feature at 2175 Angstrom, and blended [OIII]+Hbeta emission. We compare the total equivalent width of Halpha, [NII], and [SII] (W_Halpha+) with the strength of the 4000 Angstrom break (D(4000)) and the best-fit specific star formation rate, and find that all these properties are strongly correlated. This is a reassuring result, as currently most distant stellar population studies are based on just continuum emission. Furthermore, the relation between W_Halpha+ and Dn(4000) provides interesting clues to the SFHs of galaxies, as these features are sensitive to different stellar ages. We find that the correlation between W_Halpha+ and D(4000) at 0.5<z<2.0 is similar to z~0, and that the suppression of star formation in galaxies at z<2 is generally not abrupt, but a gradual process.

We present a study of the evolution of the galaxy Velocity Dispersion Function (VDF) from z=0 to z=1.5 using photometric data from the UKIDSS Ultra Deep Survey (UDS) and Newfirm Medium Band Survey (NMBS) COSMOS surveys. The VDF has been measured locally using direct kinematic measurements from the Sloan Digital Sky Survey, but direct studies of the VDF at high redshift are difficult as they require velocity dispersion measurements of many thousands of galaxies. Taylor et al. (2010) demonstrated that dynamical and stellar mass are linearly related when the structure of the galaxy is accounted for. We show that the stellar mass, size and S\'ersic index can reliably predict the velocity dispersions of SDSS galaxies. We apply this relation to galaxies at high redshift and determine the evolution of the inferred VDF. We find that the VDF at z~0.5 is very similar to the VDF at z=0. At higher redshifts, we find that the number density of galaxies with dispersions <~200 km/s is lower, but the number of high dispersion galaxies is constant or even higher. At fixed cumulative number density, the velocity dispersions of galaxies with logN[Mpc^{-3}]<-3.5 increase with time by a factor of ~1.4 from z~1.5-0, whereas the dispersions of galaxies with lower number density are approximately constant or decrease with time. The VDF appears to show less evolution than the stellar mass function, particularly at the lowest number densities. We note that these results are still somewhat uncertain and we suggest several avenues for further calibrating the inferred velocity dispersions.

Recent photometric studies have shown that early-type galaxies at fixed stellar mass were smaller and denser at earlier times. In this paper we assess that finding by deriving the dynamical mass of such a compact quiescent galaxy at z=1.8. We have obtained a high-quality spectrum with full UV-NIR wavelength coverage of galaxy NMBS-C7447 using X-Shooter on the VLT. We determined a velocity dispersion of 294 +- 51 km/s. Given this velocity dispersion and the effective radius of 1.64 +- 0.15 kpc (as determined from HST-WFC3 F160W observations) we derive a dynamical mass of 1.7 +- 0.5 x 10^11 Msun. Comparison of the full spectrum with stellar population synthesis models indicates that NMBS-C774 has a relatively young stellar population (0.40 Gyr) with little or no star formation and a stellar mass of ~1.5 x 10^11 Msun. The dynamical and photometric stellar mass are in good agreement. Thus, our study supports the conclusion that the mass densities of quiescent galaxies were indeed higher at earlier times, and this earlier result is not caused by systematic measurement errors. By combining available spectroscopic measurements at different redshifts, we find that the velocity dispersion at fixed dynamical mass was a factor of ~1.8 higher at z=1.8 compared to z=0. Finally, we show that the apparent discrepancies between the few available velocity dispersion measurements at z>1.5 are consistent with the intrinsic scatter of the mass-size relation.

We present deep near-infrared (NIR) medium-bandwidth photometry over the wavelength range 1-1.8 microns in the All-wavelength Extended Groth strip International Survey (AEGIS) and Cosmic Evolution Survey (COSMOS) fields. The observations were carried out as part of the NEWFIRM Medium-Band Survey (NMBS), an NOAO survey program on the Mayall 4m telescope on Kitt Peak using the NOAO Extremely Wide-Field Infrared Imager (NEWFIRM). In this paper, we describe the full details of the observations, data reduction and photometry for the survey. We also present a public K-selected photometric catalog, along with accurate photometric redshifts. The redshifts are computed with 37 (20) filters in the COSMOS (AEGIS) fields, combining the NIR medium-bandwidth data with existing ultraviolet (UV; Galaxy Evolution Explorer), visible and NIR (Canada-France-Hawaii Telescope and Subaru) and mid-IR (Spitzer/IRAC) imaging. We find excellent agreement with publicly available spectroscopic redshifts, with sigma_z/(1+z)~1-2% for ~4000 galaxies at z=0-3. The NMBS catalogs contain ~13,000 galaxies at z>1.5 with accurate photometric redshifts and rest-frame colors. Due to the increased spectral resolution obtained with the five NIR medium-band filters, the median 68% confidence intervals of the photometric redshifts of both quiescent and star-forming galaxies are a factor of ~2 times smaller when comparing catalogs with medium-band NIR photometry to NIR broadband photometry. We show evidence for a clear bimodal color distribution between quiescent and star-forming galaxies that persists to z~3, a higher redshift than has been probed so far.

We study the build-up of the bimodal galaxy population using the NEWFIRM Medium-Band Survey, which provides excellent redshifts and well-sampled spectral energy distributions of ~27,000 galaxies with K<22.8 at 0.4 < z < 2.2. We first show that star-forming galaxies and quiescent galaxies can be robustly separated with a two-color criterion over this entire redshift range. We then study the evolution of the number density and mass density of quiescent and star-forming galaxies, extending the results of the COMBO-17, DEEP2, and other surveys to z=2.2. The mass density of quiescent galaxies with M > 3 10^11 solar masses increases by a factor of ~10 from z=2 to the present day, whereas the mass density in star-forming galaxies is flat or decreases over the same time period. Modest mass growth by a factor of 2 of individual quiescent galaxies can explain roughly half of the strong density evolution at M>10^11 solar masses, due to the steepness of the exponential tail of the mass function. The rest of the density evolution of massive, quiescent galaxies is likely due to transformation (e.g. quenching) of the massive star-forming population, a conclusion which is consistent with the density evolution we observe for the star-forming galaxies themselves, which is flat or decreasing with cosmic time. Modest mass growth does not explain the evolution of less massive quiescent galaxies (~10^10.5 solar masses), which show a similarly steep increase in their number densities. The less massive quiescent galaxies are therefore continuously formed by transforming galaxies from the star-forming population

We derive stellar masses from SED fitting to rest-frame optical and UV fluxes for 401 star-forming galaxies at z 4, 5, and 6 from Hubble-WFC3/IR observations of the ERS combined with the deep GOODS-S Spitzer/IRAC data (and include a previously-published z 7 sample). A mass-luminosity relation with strongly luminosity-dependent M/Luv ratios is found for the largest sample (299 galaxies) at z 4. The relation M \propto L_{UV,1500}^(1.7+/-0.2) has a well-determined intrinsic sample variance of 0.5 dex. This relation is also consistent with the more limited samples at z 5-7. This z 4 mass-luminosity relation, and the well-established faint UV luminosity functions at z 4-7, are used to derive galaxy mass functions (MF) to masses M~10^8 at z 4-7. A bootstap approach is used to derive the MFs to account for the large scatter in the M--Luv relation and the luminosity function uncertainties, along with an analytical crosscheck. The MFs are also corrected for the effects of incompleteness. The incompleteness-corrected MFs are steeper than previously found, with slopes \alpha_M-1.4 to -1.6 at low masses. These slopes are, however, still substantially flatter than the MFs obtained from recent hydrodynamical simulations. We use these MFs to estimate the stellar mass density (SMD) of the universe to a fixed M_{UV,AB}<-18 as a function of redshift and find a SMD growth \propto(1+z)^{-3.4 +/-0.8} from z 7 to z 4. We also derive the SMD from the completeness-corrected MFs to a mass limit M~10^{8} Msun. Such completeness-corrected MFs and the derived SMDs will be particularly important for model comparisons as future MFs reach to lower masses.

We present Hubble Space Telescope NIC2 morphologies of a spectroscopic sample of massive galaxies at z~2.3, by extending our sample of 9 compact quiescent galaxies (r_e~0.9 kpc) with 10 massive emission-line galaxies. The emission-line galaxies are classified by the nature of their ionized emission; there are six star-forming galaxies and four galaxies hosting an active galactic nucleus (AGN). The star-forming galaxies are the largest among the emission-line galaxies, with a median size of re = 2.8 kpc. The three galaxies with the highest star formation rates (> 100M \odot/yr) have irregular and clumpy morphologies. The AGN host galaxies are more similar to the compact quiescent galaxies in terms of their structures (re~1.1 kpc) and spectral energy distributions. The total sample clearly separates into two classes in a color-mass diagram: the large star-forming galaxies that form the blue cloud, and the compact quiescent galaxies on the red sequence. However, it is unclear how or even if the two classes are evolutionary related. Three out of six massive star-forming galaxies have dense cores and thus may passively evolve into compact galaxies due to fading of outer star-forming regions. For these galaxies a reverse scenario, in which compact galaxies grow inside-out by star formation is also plausible. We do caution though that the sample is small. Nonetheless, it is evident that a Hubble sequence of massive galaxies with strongly correlated galaxy properties is already in place at z > 2.

Stellar population synthesis (SPS) models are a key ingredient of many galaxy evolution studies. Unfortunately, the models are still poorly calibrated for certain stellar evolution stages. Of particular concern is the treatment of the thermally-pulsing asymptotic giant branch (TP-AGB) phase, as different implementations lead to systematic differences in derived galaxy properties. Post-starburst galaxies are a promising calibration sample, as TP-AGB stars are thought to be most prominently visible during this phase. Here, we use post-starburst galaxies in the NEWFIRM medium-band survey to assess different SPS models. The available photometry allows the selection of a homogeneous and well-defined sample of 62 post-starburst galaxies at 0.7<z<2.0, from which we construct a well-sampled composite spectral energy distribution (SED) over the range 1200-40 000 Angstrom. The SED is well-fit by the Bruzual & Charlot SPS models, while the Maraston models do not reproduce the rest-frame optical and near-infrared parts of the SED simultaneously. When the fitting is restricted to lambda < 6000 Angstrom, the Maraston models overpredict the near-infrared luminosity, implying that these models give too much weight to TP-AGB stars. Using the flexible SPS models by Conroy et al, and assuming solar metallicity, we find that the contribution of TP-AGB stars to the integrated SED is a factor of ~3 lower than predicted by the latest Padova TP-AGB models. Whether this is due to lower bolometric luminosities, shorter lifetimes, and/or heavy dust obscuration of TP-AGB stars remains to be addressed. Altogether, our data demand a low contribution from TP-AGB stars to the SED of post-starburst galaxies.

We present an analysis of the clustering of galaxies as a function of their stellar mass at 1 < z < 2 using data from the NEWFIRM Medium Band Survey (NMBS). The precise photometric redshifts and stellar masses that the NMBS produces allows us to define a series of mass limited samples of galaxies more massive than 0.7, 1 and 3x10^10 Msun in redshift intervals centered on z = 1.1, 1.5 and 1.9 respectively. In each redshift interval we show that there exists a strong dependence of clustering strength on the stellar mass limit of the sample, with more massive galaxies showing a higher clustering amplitude on all scales. We further interpret our clustering measurements in the LCDM cosmological context using the halo model of galaxy clustering. We show that the typical halo mass of central and satellite galaxies increases with stellar mass, whereas the satellite fraction decreases with stellar mass, qualitatively the same as is seen at z < 1. We see little evidence of any redshift dependence in the stellar mass-to-halo mass relationship over our narrow redshift range. However, when we compare with similar measurements at z~0, we see clear evidence for a change in this relation. If we assume a universal baryon fraction, the ratio of stellar mass to halo mass reveals the fraction of baryons that have been converted to stars. We see that the peak in this star formation efficiency for central galaxies shifts to higher halo masses at higher redshift, moving from ~7x10^11 Msun at z~0 to ~3x10^12 Msun at z~1.5, revealing evidence of `halo downsizing'. Finally we show that for highly biased galaxy populations at z > 1 there may be a discrepancy between the measured space density and clustering and that predicted by the halo model. This could imply that there is a problem with one or more ingredients of the halo model at these redshifts, for instance the halo bias relation or the halo profile.

We study the evolution of the size - stellar mass relation for a large spectroscopic sample of galaxies in the GOODs North field up to $z \sim 3.5$. The sizes of the galaxies are measured from $\textit{K}_{s}$-band images (corresponding to rest-frame optical/NIR) from the Subaru 8m telescope. We reproduce earlier results based on photometric redshifts that the sizes of galaxies at a given mass evolve with redshift. Specifically, we compare sizes of UV-bright galaxies at a range of redshifts: Lyman break galaxies (LBGs) selected through the U-drop technique ($z \sim 2.5-3.5$), BM/BX galaxies at $z \sim 1.5-2.5$, and GALEX LBGs at low redshift ($z \sim 0.6-1.5$). The median sizes of these UV-bright galaxies evolve as $(1+z)^{-1.11\pm0.13}$ between $z \sim 0.5-3.5$. The UV-bright galaxies are significantly larger than quiescent galaxies at the same mass and redshift by $0.45\pm0.09$ dex. We also verify the correlation between color and stellar mass density of galaxies to high redshifts. The sizes of sub-mm galaxies in the same field are measured and compared with BM/BX galaxies. We find that median half-light radii of SMGs is $2.90 \pm 0.45$ kpc and there is little difference in their size distribution to the UV-bright star forming galaxies.

We present the Spitzer IRAC/MUSYC Public Legacy Survey in the Extended CDF-South (SIMPLE), which consists of deep IRAC observations covering the ~1,600 arcmin^2 area surrounding GOODS-S. The limiting magnitudes of the SIMPLE IRAC mosaics typically are 23.8, 23.6, 21.9, and 21.7, at 3.6 um, 4.5 um, 5.8 um, and 8.0 um, respectively (5-sigma total point source magnitudes in AB). The SIMPLE IRAC images are combined with the 10'x15' GOODS IRAC mosaics in the center. We give detailed descriptions of the observations, data reduction, and properties of the final images, as well as the detection and photometry methods used to build a catalog. Using published optical and near-infrared data from the Multiwavelength Survey by Yale-Chile (MUSYC), we construct an IRAC-selected catalog, containing photometry in UBVRIz'JHK, [3.6 um], [4.5 um], [5.8 um], and [8.0 um]. The catalog contains 43,782 sources with S/N > 5 at 3.6 um, 19,993 of which have 13-band photometry. We compare this catalog to the publicly available MUSYC and FIREWORKS catalogs and discuss the differences. Using a high signal-to-noise sub-sample of 3,391 sources with ([3.6] + [4.5])/2 < 21.2, we investigate the star formation rate history of massive galaxies out to z ~ 1.8. We find that at z ~ 1.8 at least 30% +/-7% of the most massive galaxies (Mstar > 10^11 Msol) are passively evolving, in agreement with earlier results from surveys covering less area.

We present an analysis of the far-infrared (FIR) spectral energy distributions (SEDs) of two massive K-selected galaxies at z = 2.122 and z = 2.024 detected at 24um, 70um, 160um by Spitzer, 250um, 350um, 500um by BLAST, and 870um by APEX. The large wavelength range of these observations and the availability of spectroscopic redshifts allow us to unambiguously identify the peak of the redshifted thermal emission from dust at ~ 300um. The SEDs of both galaxies are reasonably well fit by synthetic templates of local galaxies with L_IR ~ 10^{11} L_{sun} -- 10^{12} L_{sun} yet both galaxies have L_{IR} ~ 10^{13} L_{sun}. This suggests that these galaxies are not high redshift analogues of the Hyper-LIRGs/ULIRGs used in local templates, but are instead "scaled up" versions of local ULIRGs/LIRGs. For both galaxies the star formation rate determined from the best-fit FIR SEDs (SFR(L_{IR})) agrees with the SFR determined from the dust corrected H-alpha luminosity (SFR(H-alpha)) to within a factor of ~ 2; however, when the SFR of these galaxies is estimated using only the observed 24um flux and the standard luminosity-dependent template method (SFR(24um)), it systematically overestimates the SFR by as much as a factor of 6. A larger sample of 24 K-selected galaxies at z ~ 2.3 drawn from the Kriek et al. (2008) GNIRS sample shows the same trend between SFR(24um) and SFR(H-alpha). Using that sample we show that SFR(24um) and SFR(H-alpha) are in better agreement when SFR(24um) is estimated using the log average of local templates rather than selecting a single luminosity-dependent template, because this incorporates lower luminosity templates. The better agreement between SFRs from lower luminosity templates suggests that the FIR SEDs of the BLAST-detected galaxies may be typical for massive galaxies at z ~ 2, and that the majority are scaled up versions of lower luminosity local galaxies.

[Abridged] We use the NEWFIRM Medium-Band Survey (NMBS) to characterize the properties of a mass-complete sample of 14 galaxies at 3.0<z<4.0 with M_star>2.5x10^11 Msun, and to derive more accurate measurements of the high-mass end of the stellar mass function (SMF) of galaxies at z=3.5, with significantly reduced contributions from photometric redshift errors and cosmic variance to the total error budget of the SMF. The typical very massive galaxy at z=3.5 is red and faint in the observer's optical, with median r=26.1, and rest-frame U-V=1.6. About 60% of the sample have optical colors satisfying either the U- or the B-dropout color criteria, although ~50% of these galaxies have r>25.5. About 30% of the sample has SFRs from SED modeling consistent with zero. However, >80% of the sample is detected at 24 micron, with total infrared luminosities in the range (0.5-4.0)x10^13 Lsun. This implies the presence of either dust-enshrouded starburst activity (with SFRs of 600-4300 Msun/yr) and/or highly-obscured active galactic nuclei (AGN). The contribution of galaxies with M_star>2.5x10^11 Msun to the total stellar mass budget at z=3.5 is ~8%. We find an evolution by a factor of 2-7 and 3-22 from z~5 and z~6, respectively, to z=3.5. The previously found disagreement at the high-mass end between observed and model-predicted SMFs is now significant at the 3sigma level. However, systematic uncertainties dominate the total error budget, with errors up to a factor of ~8 in the densities, bringing the observed SMF in marginal agreement with the predicted SMF. Additional systematic uncertainties on the high-mass end could be introduced by either 1) the intense star-formation and/or the very common AGN activities as inferred from the MIPS 24 micron detections, and/or 2) contamination by a significant population of massive, old, and dusty galaxies at z~2.6.

With a complete, mass-selected sample of quiescent galaxies from the NEWFIRM Medium-Band Survey (NMBS), we study the stellar populations of the oldest and most massive galaxies (>10^11 Msun) to high redshift. The sample includes 570 quiescent galaxies selected based on their extinction-corrected U-V colors out to z=2.2, with accurate photometric redshifts, sigma_z/(1+z)~2%, and rest-frame colors, sigma_U-V~0.06 mag. We measure an increase in the intrinsic scatter of the rest-frame U-V colors of quiescent galaxies with redshift. This scatter in color arises from the spread in ages of the quiescent galaxies, where we see both relatively quiescent red, old galaxies and quiescent blue, younger galaxies towards higher redshift. The trends between color and age are consistent with the observed composite rest-frame spectral energy distributions (SEDs) of these galaxies. The composite SEDs of the reddest and bluest quiescent galaxies are fundamentally different, with remarkably well-defined 4000A- and Balmer-breaks, respectively. Some of the quiescent galaxies may be up to 4 times older than the average age- and up to the age of the universe, if the assumption of solar metallicity is correct. By matching the scatter predicted by models that include growth of the red sequence by the transformation of blue galaxies to the observed intrinsic scatter, the data indicate that most early-type galaxies formed their stars at high redshift with a burst of star formation prior to migrating to the red sequence. The observed U-V color evolution with redshift is weaker than passive evolution predicts; possible mechanisms to slow the color evolution include increasing amounts of dust in quiescent galaxies towards higher redshift, red mergers at z<1, and a frosting of relatively young stars from star formation at later times.

We investigate the Spitzer/IRAC properties of 36 z~7 z-dropout galaxies and 3 z~8 Y-dropout galaxies derived from deep/wide-area WFC3/IR data of the Early Release Science, the ultradeep HUDF09, and wide-area NICMOS data. We fit stellar population synthesis models to the SEDs to derive mean redshifts, stellar masses, and ages. The z~7 galaxies are best characterized by substantial ages (>100 Myr) and M/L_V~0.2. The main trend with decreasing luminosity is that of bluing of the far-UV slope from beta~-2.0 to beta~-3.0. This can be explained by decreasing metallicity, except for the lowest luminosity galaxies (0.1 L*{z=3}), where low metallicity and smooth SFHs fail to match the blue far-UV and moderately red H-[3.6] color. Such colors may require episodic SFHs with short periods of activity and quiescence ("on-off" cycles) and/or a contribution from emission lines. The stellar mass of our sample of z~7 star forming galaxies correlates with SFR according to log M* = 8.70 (+-0.09) + 1.06(+-0.10) log SFR, implying star formation may have commenced at z>10. No galaxies are found with SFRs much higher or lower than the past averaged SFR suggesting that the typical star formation timescales are probably a substantial fraction of the Hubble time. We report the first IRAC detection of Y-dropout galaxies at z~8. The average rest-frame U-V=0.3 (AB) of the 3 galaxies are similar to faint z~7 galaxies, implying similar M/L. The stellar mass density to M_{UV,AB}<-18 is rho*(z=8) = 1.8^{+0.7}_{-1.0} x 10^6 M_sun Mpc^-3, following log rho*(z)= 10.6(+-0.6) - 4.4(+-0.7) log(1+z) [M_sun Mpc^-3] over 3<z< 8.

The presence of massive, compact, quiescent galaxies at z>2 presents a major challenge for theoretical models of galaxy formation and evolution. Using one of the deepest large public near-IR surveys to date, we investigate in detail the correlations between star formation and galaxy structural parameters (size, stellar mass, and surface density) from z=2 to the present. At all redshifts, massive quiescent galaxies (i.e. those with little or no star formation) occupy the extreme high end of the surface density distribution and follow a tight mass-size correlation, while star-forming galaxies show a broad range of both densities and sizes. Conversely, galaxies with the highest surface densities comprise a nearly-homogeneous population with little or no ongoing star formation, while less dense galaxies exhibit high star-formation rates and varying levels of dust obscuration. Both the sizes and surface densities of quiescent galaxies evolve strongly from z=2-0; we parameterize this evolution for both populations with simple power law functions and present best-fit parameters for comparison to future theoretical models. Higher-mass quiescent galaxies undergo faster structural evolution, consistent with previous results. Interestingly, star-forming galaxies' sizes and densities evolve at rates similar to those of quiescent galaxies. It is therefore possible that the same physical processes drive the structural evolution of both populations, suggesting that "dry mergers" may not be the sole culprit in this size evolution.

We use a robust sample of 11 z~7 galaxies (z-dropouts) to estimate the stellar mass density of the universe when it was only ~750 Myr old. We combine the very deep optical to near-Infrared photometry from the HST ACS and NICMOS cameras with mid-Infrared Spitzer IRAC imaging available through the GOODS program. After carefully removing the flux from contaminating foreground sources we have obtained reliable photometry in the 3.6 and 4.5 micron IRAC channels. The spectral shapes of these sources, including their rest frame optical colors, strongly support their being at z~7 with a mean photometric redshift of <z>=7.2+/-0.5. We use Bruzual & Charlot (2003) synthetic stellar population models to constrain their stellar masses and star formation histories. We find stellar masses that range over 0.1 -12x10^9 M_sol and average ages from 20 Myr to up to 425 Myr with a mean of ~300 Myr, suggesting that in some of these galaxies most of the stars were formed at z>8 (and probably at z>~10). The best fits to the observed SEDs are consistent with little or no dust extinction, in agreement with recent results at z~4-8. The star formation rates (SFR) are in the range from 5-20 M_sol/yr. From this sample we measure a stellar mass density of 6.6_{-3.3}^{+5.4}x10^5 M_sol/Mpc^3 to a limit of M_{UV,AB}<-20 (or 0.4 L*(z=3)). Combined with a fiducial lower limit for their ages (80 Myr) this implies a maximum SFR density of 0.008 M_sol/yr/Mpc^3. This is well below the critical level needed to reionize the universe at z~8 using standard assumptions. However, this result is based on luminous sources (>L*) and does not include the dominant contribution of the fainter galaxies. Strikingly, we find that the specific SFR is constant from z~7 to z~2 but drops substantially at more recent times.

We study the Spitzer Infrared Array Camera (IRAC) mid-infrared (rest-frame optical) fluxes of 14 newly WFC3/IR-detected z=7 z_{850}-dropout galaxies and 5 z=8 Y_{105}-dropout galaxies. The WFC3/IR depth and spatial resolution allow accurate removal of contaminating foreground light, enabling reliable flux measurements at 3.6 micron and 4.5 micron. None of the galaxies are detected to [3.6]=26.9 (AB, 2 sigma), but a stacking analysis reveals a robust detection for the z_{850}-dropouts and an upper limit for the Y_{105}-dropouts. We construct average broadband SEDs using the stacked ACS, WFC3, and IRAC fluxes and fit stellar population synthesis models to derive mean redshifts, stellar masses, and ages. For the z_{850}-dropouts, we find z=6.9^{+0.1}_{-0.1}, (U-V)_{rest}=0.4, reddening A_V=0, stellar mass M*=1.2^{+0.3}_{-0.6} x 10^9 M_sun (Salpeter IMF). The best-fit ages ~300Myr, M/L_V=0.2, and SSFR=1.7Gyr^{-1} are similar to values reported for luminous z=7 galaxies, indicating the galaxies are smaller but not younger. The sub-L* galaxies observed here contribute significantly to the stellar mass density and under favorable conditions may have provided enough photons for sustained reionization at 7<z<11. In contrast, the z=8.3^{+0.1}_{-0.2} Y_{105}-dropouts have stellar masses that are uncertain by 1.5 dex due to the near-complete reliance on far-UV data. Adopting the 2 sigma upper limit on the M/L(z=8), the stellar mass density to M_{UV,AB} < -18 declines from rho*(z=7)=3.7^{+1.0}_{-1.8} x 10^6 M_sun Mpc^{-3} to rho*(z=8) < 8 x 10^5 M_sun Mpc^{-3}, following (1+z)^{-6} over 3<z<8. Lower masses at z=8 would signify more dramatic evolution, which can be established with deeper IRAC observations, long before the arrival of the James Webb Space Telescope.

We study the growth of massive galaxies from z=2 to the present using data from the NEWFIRM Medium Band Survey. The sample is selected at a constant number density of n=2x10^-4 Mpc^-3, so that galaxies at different epochs can be compared in a meaningful way. We show that the stellar mass of galaxies at this number density has increased by a factor of ~2 since z=2, following the relation log(M)=11.45-0.15z. In order to determine at what physical radii this mass growth occurred we construct very deep stacked rest-frame R-band images at redshifts z=0.6, 1.1, 1.6, and 2.0. These image stacks of typically 70-80 galaxies enable us to characterize the stellar distribution to surface brightness limits of ~28.5 mag/arcsec^2. We find that massive galaxies gradually built up their outer regions over the past 10 Gyr. The mass within a radius of r=5 kpc is nearly constant with redshift whereas the mass at 5-75 kpc has increased by a factor of ~4 since z=2. Parameterizing the surface brightness profiles we find that the effective radius and Sersic n parameter evolve as r_e~(1+z)^-1.3 and n~(1+z)^-1.0 respectively. The data demonstrate that massive galaxies have grown mostly inside-out, assembling their extended stellar halos around compact, dense cores with possibly exponential radial density distributions. Comparing the observed mass evolution to the average star formation rates of the galaxies we find that the growth is likely dominated by mergers, as in-situ star formation can only account for ~20% of the mass build-up from z=2 to z=0. The main uncertainties in this study are possible redshift-dependent systematic errors in the total stellar masses and the conversion from light-weighted to mass-weighted radial profiles.

We select 25,000 galaxies from the NEWFIRM Medium Band Survey (NMBS) to study the rest-frame U-V color distribution of galaxies at 0 < z < 2.5. The five unique NIR filters of the NMBS enable the precise measurement of photometric redshifts and rest-frame colors for 9,900 galaxies at 1 < z < 2.5. The rest-frame U-V color distribution at all z<~2.5 is bimodal, with a red peak, a blue peak, and a population of galaxies in between (the green valley). Model fits to the optical-NIR SEDs and the distribution of MIPS-detected galaxies indicate that the colors of galaxies in the green valley are determined largely by the amount of reddening by dust. This result does not support the simplest interpretation of green valley objects as a transition from blue star-forming to red quiescent galaxies. We show that correcting the rest-frame colors for dust reddening allows a remarkably clean separation between the red and blue sequences up to z~2.5. Our study confirms that dusty starburst galaxies can contribute a significant fraction to red sequence samples selected on the basis of a single rest-frame color (i.e. U-V), so extra care must be taken if samples of truly "red and dead" galaxies are desired. Interestingly, of galaxies detected at 24 microns, 14% remain on the red sequence after applying the reddening correction.

Using a sample of nine massive compact galaxies at z ~ 2.3 with rest-frame optical spectroscopy and comprehensive U through 8um photometry we investigate how assumptions in SED modeling change the stellar mass estimates of these galaxies, and how this affects our interpretation of their size evolution. The SEDs are fit to Tau-models with a range of metallicities, dust laws, as well as different stellar population synthesis codes. These models indicate masses equal to, or slightly smaller than our default masses. The maximum difference is 0.16 dex for each parameter considered, and only 0.18 dex for the most extreme combination of parameters. Two-component populations with a maximally old stellar population superposed with a young component provide reasonable fits to these SEDs using the models of Bruzual & Charlot (2003); however, using models with updated treatment of TP-AGB stars the fits are poorer. The two-component models predict masses that are 0.08 to 0.22 dex larger than the Tau-models. We also test the effect of a bottom-light IMF and find that it would reduce the masses of these galaxies by 0.3 dex. Considering the range of allowable masses from the Tau-models, two-component fits, and IMF, we conclude that on average these galaxies lie below the mass-size relation of galaxies in the local universe by a factor of 3-9, depending on the SED models used.

We present a detailed study of the spatially resolved kinematics, star-formation and stellar mass in a highly amplified galaxy at z=4.92 behind the lensing cluster MS1358+62. We use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the velocity structure of the galaxy from the nebular [OII] emission. Using a detailed gravitational lens model, we account for lensing amplification factor 12.+/-2.0 and find that this intrinsically L* galaxy has a stellar mass of M*=7+/-2x10^8Mo, a dynamical mass of Mdyn=3+/-1x10^9csc^2(i)Mo (within of 2kpc) and a star-formation rate of 42+/-8Mo/yr. The source-plane UV/optical morphology of this galaxy is dominated by five discrete star-forming regions. Exploiting the dynamical information we derive masses for individual star-forming regions of Mcl~10^(8-9)Mo with sizes of ~200pc. We find that, at a fixed size, the star-formation rate density within these HII regions is approximately two orders of magnitude greater than those observed in local spiral/starburst galaxies, but consistent with the most massive HII regions in the local Universe such as 30Dor. Finally, we compare the spatially resolved nebular emission line velocity with the Ly-alpha and UV-ISM lines and find that this galaxy is surrounded by a galactic scale outflow in which the Ly-alpha appears redshifted by ~150km/s and the UV-ISM lines blue-shifted by -200km/s from the (systemic) nebular emission. The velocity structure of the outflow mirrors that of the nebular emission suggesting the outflow is young (~15Myr), and has yet to burst out of the system. Taken together, these results suggest that this young galaxy is undergoing its first major epoch of mass assembly.

Several recent studies have shown that about half of the massive galaxies at z~2 are in a quiescent phase. Moreover, these galaxies are commonly found to be ultra-compact with half-light radii of ~1 kpc. We have obtained a ~29 hr spectrum of a typical quiescent, ultra-dense galaxy at z=2.1865 with the Gemini Near-Infrared Spectrograph. The spectrum exhibits a strong optical break and several absorption features, which have not previously been detected in z>2 quiescent galaxies. Comparison of the spectral energy distribution with stellar population synthesis models implies a low star formation rate (SFR) of 1-3 Msol/yr, an age of 1.3-2.2 Gyr, and a stellar mass of ~2x10^11 Msol. We detect several faint emission lines, with emission-line ratios of [NII]/Halpha, [SII]/Halpha and [OII]/[OIII] typical of low-ionization nuclear emission-line regions. Thus, neither the stellar continuum nor the nebular emission implies active star formation. The current SFR is <1% of the past average SFR. If this galaxy is representative of compact quiescent galaxies beyond z=2, it implies that quenching of star formation is extremely efficient and also indicates that low luminosity active galactic nuclei (AGNs) could be common in these objects. Nuclear emission is a potential concern for the size measurement. However, we show that the AGN contributes <8% to the rest-frame optical emission. A possible post-starburst population may affect size measurements more strongly; although a 0.5 Gyr old stellar population can make up <10% of the total stellar mass, it could account for up to ~40% of the optical light. Nevertheless, this spectrum shows that this compact galaxy is dominated by an evolved stellar population.

Recent studies have found that the oldest and most luminous galaxies in the early Universe are surprisingly compact, having stellar masses similar to present-day elliptical galaxies but much smaller sizes. This finding has attracted considerable attention as it suggests that massive galaxies have grown by a factor of ~five in size over the past ten billion years. A key test of these results is a determination of the stellar kinematics of one of the compact galaxies: if the sizes of these objects are as extreme as has been claimed, their stars are expected to have much higher velocities than those in present-day galaxies of the same mass. Here we report a measurement of the stellar velocity dispersion of a massive compact galaxy at redshift z=2.186, corresponding to a look-back time of 10.7 billion years. The velocity dispersion is very high at 510 (+165, -95) km/s, consistent with the mass and compactness of the galaxy inferred from photometric data and indicating significant recent structural and dynamical evolution of massive galaxies. The uncertainty in the dispersion was determined from simulations which include the effects of noise and template mismatch. However, we caution that some subtle systematic effect may influence the analysis given the low signal-to-noise ratio of our spectrum.

We present SED modeling of a sample of 34 K-selected galaxies at z ~ 2.3. These galaxies have NIR spectroscopy that samples the rest-frame Balmer/4000A break as well as deep photometry in thirteen broadband filters. New to our analysis is IRAC data that extend the SEDs into the rest-frame NIR. Comparing parameters determined from SED fits with and without the IRAC data we find that the IRAC photometry significantly improves the confidence intervals of Tau, A_v, stellar mass, and SFR for individual galaxies, but does not systematically alter the mean parameters of the sample. We use the IRAC data to assess how well current stellar population synthesis codes describe the rest-frame NIR SEDs of young galaxies where discrepancies between treatments of the TP-AGB phase of stellar evolution are most pronounced. The models of Bruzual & Charlot (2003), Maraston (2005), and Charlot & Bruzual (2008) all successfully reproduce the SEDs of our galaxies with < 5% differences in the quality of fit; however, the best-fit masses from each code differ systematically by as much as a factor of 1.5, and other parameters vary more, up to factors of 2-3. A comparison of best-fit stellar population parameters from different SPS codes, dust laws, and metallicities shows that the choice of SPS code is the largest systematic uncertainty in most parameters, and that systematic uncertainties are typically larger than the formal random uncertainties. The SED fitting confirms our previous result that galaxies with strongly suppressed star formation account for ~ 50% of the K-bright population at z ~ 2.3; however, the uncertainty in this fraction is large due to systematic differences in the SSFRs derived from the three SPS models.

Over the last decade optical spectroscopic surveys have characterized the low redshift galaxy population and uncovered populations of star-forming galaxies back to z ~ 7. This work has shown that the primary epoch of galaxy building and black hole growth occurs at redshifts of 2 to 3. The establishment of the concordance LCDM cosmology shifted the focus of galaxy population studies from constraining cosmological parameters to characterizing the processes which regulate the formation and evolution of galaxies.In the next decade, high redshift observers will attempt to formulate a coherent evolutionary picture connecting galaxies in the high redshift Universe to galaxies today. In order to link galaxy populations at different redshifts, we must not only characterize their evolution in a systematic way, we must establish which physical processes are responsible for it. Considerable progress has already been made in understanding how galaxies evolved from z ~ 1 to the present day. Large spectroscopic surveys in the near infrared are required to push these studies back towards the main epoch of galaxy building. Only then will we understand the full story of the formation of L* galaxies like our own Milky Way. A large near-IR spectroscopic survey will also provide the calibration needed to avoid systematics in the large photometric programs proposed to study the nature of dark matter and dark energy. We provide an outline design for a multi-object 0.4 to 1.8 micron spectrograph, which could be placed on an existing telescope, and which would allow a full characterization of the galaxy population out to z ~ 2. We strongly recommend a serious further study to design a real instrument, which will be required for galaxy formation studies to advance to the next frontier.

Recent studies have shown that massive quiescent galaxies at high redshift are much more compact than present-day galaxies of the same mass. Here we compare the radial stellar density profiles and the number density of a sample of massive galaxies at z ~ 2.3 to nearby massive elliptical galaxies. We confirm that the average stellar densities of the z ~ 2.3 galaxies within the effective radius, rho(<r_e), are two orders of magnitude higher than those of local elliptical galaxies of the same stellar mass. However, we also find that the densities measured within a constant physical radius of 1 kpc, rho(<1 kpc), are higher by a factor of 2-3 only. This suggests that inside-out growth scenarios are plausible, in which the compact high redshift galaxies make up the centers of normal nearby ellipticals. The compact galaxies are common at high redshift, which enables us to further constrain their evolution by requiring that the number density of their descendants does not exceed constraints imposed by the z=0 galaxy mass function. We infer that size growth must be efficient, with (r_{1+2}/r_1) ~ (M_{1+2}/M_1)^2. A simple model where compact galaxies with masses ~ 10^11 Msun primarily grow through minor mergers produces descendants with the approximate sizes, stellar densities, and number density of elliptical galaxies with masses 2-3 x10^11 Msun in the local Universe. We note that this model also predicts evolution in the M_BH - sigma relation, such that the progenitors of elliptical galaxies have lower black hole masses at fixed velocity dispersion. The main observational uncertainty is the conversion from light to mass; measurements of kinematics are needed to calibrate the masses and stellar densities of the high redshift galaxies.

We present a new K-selected, optical-to-near-infrared photometric catalog of the Extended Chandra Deep Field South (ECDFS), making it publicly available to the astronomical community. The dataset is founded on publicly available imaging, supplemented by original zJK imaging data obtained as part of the MUltiwavelength Survey by Yale-Chile (MUSYC). The final photometric catalog consists of photometry derived from nine band U-K imaging covering the full 0.5x0.5 sq. deg. of the ECDFS, plus H band data for approximately 80% of the field. The 5sigma flux limit for point-sources is K = 22.0 (AB). This is also the nominal completeness and reliability limit of the catalog: the empirical completeness for 21.75 < K < 22.00 is 85+%. We have verified the quality of the catalog through both internal consistency checks, and comparisons to other existing and publicly available catalogs. As well as the photometric catalog, we also present catalogs of photometric redshifts and restframe photometry derived from the ten band photometry. We have collected robust spectroscopic redshift determinations from published sources for 1966 galaxies in the catalog. Based on these sources, we have achieved a (1sigma) photometric redshift accuracy of Dz/(1+z) = 0.036, with an outlier fraction of 7.8%. Most of these outliers are X-ray sources. Finally, we describe and release a utility for interpolating restframe photometry from observed SEDs, dubbed InterRest. Particularly in concert with the wealth of already publicly available data in the ECDFS, this new MUSYC catalog provides an excellent resource for studying the changing properties of the massive galaxy population at z < 2. (Abridged)

Deep near-infrared imaging surveys allow us to select and study distant galaxies in the rest-frame optical, and have transformed our understanding of the early Universe. As the vast majority of K- or IRAC-selected galaxies is too faint for spectroscopy, the interpretation of these surveys relies almost exclusively on photometric redshifts determined from fitting templates to the broad-band photometry. The best-achieved accuracy of these redshifts Delta(z)/(1+z) ~ 0.06 at z>1.5, which is sufficient for determining the broad characteristics of the galaxy population but not for measuring accurate rest-frame colors, stellar population parameters, or the local galaxy density. We have started a near-infrared imaging survey with the NEWFIRM camera on the Kitt Peak 4m telescope to greatly improve the accuracy of photometric redshifts in the range 1.5<z<3.5. The survey uses five medium-bandwidth filters, which provide crude "spectra" over the wavelength range 1-1.8 micron for all objects in the 27.6 x 27.6 arcmin NEWFIRM field. In this first paper, we illustrate the technique by showing medium band NEWFIRM photometry of several galaxies at 1.7<z<2.7 from the near-infrared spectroscopic sample of Kriek et al. (2008). The filters unambiguously pinpoint the location of the redshifted Balmer break in these galaxies, enabling very accurate redshift measurements. The full survey will provide similar data for ~8000 faint K-selected galaxies at z>1.5 in the COSMOS and AEGIS fields. The filter set also enables efficient selection of exotic objects such as high redshift quasars, galaxies dominated by emission lines, and very cool brown dwarfs; we show that late T and candidate "Y" dwarfs could be identified using only two of the filters.

The existence of massive galaxies with strongly suppressed star formation at z~2.3, identified in a previous paper, suggests that a red sequence may already be in place beyond z=2. In order to test this hypothesis, we study the rest-frame U-B color distribution of massive galaxies at 2<z<3. The sample is drawn from our near-infrared spectroscopic survey for massive galaxies. The color distribution shows a statistically significant (>3 sigma) red sequence, which hosts ~60% of the stellar mass at the high-mass end. The red-sequence galaxies have little or no ongoing star formation, as inferred from both emission-line diagnostics and stellar continuum shapes. Their strong Balmer breaks and their location in the rest-frame U-B, B-V plane indicate that they are in a post-starburst phase, with typical ages of ~0.5-1.0 Gyr. In order to study the evolution of the red sequence, we compare our sample with spectroscopic massive galaxy samples at 0.02<z<0.045 and 0.6<z<1.0. The rest-frame U-B color reddens by ~0.25 mag from z~2.3 to the present at a given mass. Over the same redshift interval, the number and stellar mass density on the high-mass end (>10^11 Msol) of the red sequence grow by factors of ~8 and ~6, respectively. We explore simple models to explain the observed evolution. Passive evolution models predict too strong d(U-B), and produce z~0 galaxies that are too red. More complicated models that include aging, galaxy transformations, and red mergers can explain both the number density and color evolution of the massive end of the red sequence between z~2.3 and the present.

Using deep near-infrared spectroscopy Kriek et al. (2006) found that ~45% of massive galaxies at z~2.3 have evolved stellar populations and little or no ongoing star formation. Here we determine the sizes of these quiescent galaxies using deep, high-resolution images obtained with HST/NIC2 and laser guide star-assisted Keck/AO. Considering that their median stellar mass is 1.7x10^11 Solar masses the galaxies are remarkably small, with a median effective radius of 0.9 kpc. Galaxies of similar mass in the nearby Universe have sizes of ~5 kpc and average stellar densities which are two orders of magnitude lower than the z~2.3 galaxies. These results extend earlier work at z~1.5 and confirm previous studies at z>2 which lacked spectroscopic redshifts and imaging of sufficient resolution to resolve the galaxies. Our findings demonstrate that fully assembled early-type galaxies make up at most ~10% of the population of K-selected quiescent galaxies at z~2.3, effectively ruling out simple monolithic models for their formation. The galaxies must evolve significantly after z~2.3, through dry mergers or other processes, consistent with predictions from hierarchical models.

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have completed a near-infrared spectroscopic survey for K-bright galaxies at z~2.3, selected from the MUSYC survey. We derived spectroscopic redshifts from emission lines or from continuum features and shapes for all 36 observed galaxies. The continuum redshifts are driven by the Balmer/4000 Angstrom break, and have an uncertainty in dz/(1+z) of <0.019. We use this unique sample to determine, for the first time, how accurately redshifts and other properties of massive high-redshift galaxies can be determined from broadband photometric data alone. We find that the photometric redshifts of the galaxies in our sample have a systematic error of 0.08 and a random error of 0.13 in dz/(1+z). The systematic error can be reduced by using optimal templates and deep photometry; the random error, however, will be hard to reduce below 5%. The spectra lead to significantly improved constraints for stellar population parameters. For most quantities this improvement is about equally driven by the higher spectral resolution and by the much reduced redshift uncertainty. Properties such as the age, A_V, current star formation rate, and the star formation history are generally very poorly constrained with broadband data alone. Interestingly stellar masses and mass-to-light ratios are among the most stable parameters from broadband data. Nevertheless, photometric studies may overestimate the number of massive galaxies at 2<z<3, and thus underestimate the evolution of the stellar mass density. Finally,the spectroscopy supports our previous finding that red galaxies dominate the high-mass end of the galaxy population at z=2-3.

We present HST NICMOS+ACS and Spitzer IRAC+MIPS observations of 41 galaxies at 2<z<3.5 in the FIRES MS1054 field with red and blue rest-frame optical colors. About half of the galaxies are very compact (effective radii r_e < 1 kpc) at rest-frame optical wavelengths, the others are extended (1< r_e < 10 kpc). For reference, 1 kpc corresponds to 0.12 arcsec at z=2.5 in the adopted cosmology. We separate actively star forming galaxies from quiescent galaxies by modeling their rest-frame UV-NIR SEDs. The star forming galaxies span the full range of sizes, while the quiescent galaxies all have r_e<2kpc. In the redshift range where MIPS 24 micron imaging is a sensitive probe of re-radiated dust emission (z<2.5), the 24 micron fluxes confirm that the light of the small quiescent galaxies is dominated by old stars, rather than dust-enshrouded star formation or AGN activity. The inferred surface mass densities and velocity dispersions for the quiescent galaxies are very high compared to those in local galaxies. The galaxies follow a Kormendy relation (between surface brightness and size) with approximately the same slope as locally, but shifted to brighter surface brightnesses, consistent with a mean stellar formation redshift of z_f~5. This paper demonstrates a direct relation between star formation activity and size at z~2.5, and the existence of a significant population of massive, extremely dense, old stellar systems without readily identifiable counterparts in the local universe.

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have assembled a complete sample of 20 K-selected galaxies at 2.0<z<2.7 with high quality near-infrared spectra. As described in a previous paper, 9 of these 20 galaxies have strongly suppressed star formation and no detected emission lines. The present paper concerns the 11 galaxies with detected Halpha emission, and studies the origin of the line emission using the GNIRS spectra and follow-up observations with SINFONI on the VLT. Based on their [NII]/Halpha ratios, the spatial extent of the line emission and several other diagnostics, we infer that four of the eleven emission-line galaxies host narrow line active galactic nuclei (AGNs). The AGN host galaxies have stellar populations ranging from evolved to star-forming. Combining our sample with a UV-selected galaxy sample at the same redshift that spans a broader range in stellar mass, we find that black-hole accretion is more effective at the high-mass end of the galaxy distribution (~2.9x10^11 Msun) at z~2.3. Furthermore, by comparing our results with SDSS data, we show that the AGN activity in massive galaxies has decreased significantly between z~2.3 and z~0. AGNs with similar normalized accretion rates as those detected in our K-selected galaxies reside in less massive galaxies (~4.0x10^10 Msun) at low redshift. This is direct evidence for downsizing of AGN host galaxies. Finally, we speculate that the typical stellar mass-scale of the actively accreting AGN host galaxies, both at low and at high redshift, might be similar to the mass-scale at which star-forming galaxies seem to transform into red, passive systems.

We study the clustering properties of K-selected galaxies at 2<z<3.5 using deep multiwavelength imaging in three fields from the MUSYC survey. These are the first measurements to probe the spatial correlation function of K-selected galaxies in this redshift range on large scales, allowing for robust conclusions about the dark matter halos that host these galaxies. K-selected galaxies with K<21 have a correlation length r_0 ~ 6 h^-1 Mpc, larger than typical values found for optically-selected galaxies. The correlation length does not depend on K-band magnitude, but it does increase strongly with color; the J-K>2.3 distant red galaxies (DRGs) have r_0 ~ 11 h^-1 Mpc. Contrary to findings for optically-selected galaxies, K-selected galaxies that are faint in the R-band cluster more strongly than brighter galaxies. These results suggest that a color-density relation was in place at z>2. Our results indicate that K-bright blue galaxies and K-bright red galaxies are fundamentally different, as they have different clustering properties. Using a simple model of one galaxy per halo, we infer halo masses ~ 5 x 10^12 M_sun for K<21 galaxies and ~ 2 x 10^13 M_sun for DRGs. A comparison of the observed space density of DRGs to the density of their host halos suggests large halo occupation numbers; however, this result is at odds with the lack of a strong small-scale excess in the angular correlation function. We find that the z=0 descendants of the galaxies considered here reside primarily in groups and clusters. [abridged]

We present deep near-infrared JHK imaging of four 10'x10' fields. The observations were carried out as part of the Multiwavelength Survey by Yale-Chile (MUSYC) with ISPI on the CTIO 4m telescope. The typical point source limiting depths are J~22.5, H~21.5, and K~21 (5sigma; Vega). The effective seeing in the final images is ~1.0". We combine these data with MUSYC UBVRIz imaging to create K-selected catalogs that are unique for their uniform size, depth, filter coverage, and image quality. We investigate the rest-frame optical colors and photometric redshifts of galaxies that are selected using common color selection techniques, including distant red galaxies (DRGs), star-forming and passive BzKs, and the rest-frame UV-selected BM, BX, and Lyman break galaxies (LBGs). These techniques are effective at isolating large samples of high redshift galaxies, but none provide complete or uniform samples across the targeted redshift ranges. The DRG and BM/BX/LBG criteria identify populations of red and blue galaxies, respectively, as they were designed to do. The star-forming BzKs have a very wide redshift distribution, a wide range of colors, and may include galaxies with very low specific star formation rates. In comparison, the passive BzKs are fewer in number, have a different distribution of K magnitudes, and have a somewhat different redshift distribution. By combining these color selection criteria, it appears possible to define a reasonably complete sample of galaxies to our flux limit over specific redshift ranges. However, the redshift dependence of both the completeness and sampled range of rest-frame colors poses an ultimate limit to the usefulness of these techniques.

We present first results of a spectroscopic survey targeting K-selected galaxies at z=2.0-2.7 using the GNIRS instrument on Gemini-South. We obtained near-infrared spectra with a wavelength coverage of 1.0-2.5 micron for 26 K-bright galaxies (K<19.7) selected from the MUSYC survey using photometric redshifts. We successfully derived spectroscopic redshifts for all 26 galaxies using rest-frame optical emission lines or the redshifted Balmer/4000 Angstrom break. Twenty galaxies have spectroscopic redshifts in the range 2.0<z<2.7, for which bright emission lines like Halpha and [OIII] fall in atmospheric windows. Surprisingly, we detected no emission lines for nine of these 20 galaxies. The median 2 sigma upper limit on the rest-frame equivalent width of Halpha for these nine galaxies is ~10 Angstrom. The stellar continuum emission of these same nine galaxies is best fitted by evolved stellar population models. The best-fit star formation rate (SFR) is zero for five out of nine galaxies, and consistent with zero within 1 sigma for the remaining four. Thus, both the Halpha measurements and the independent stellar continuum modeling imply that 45% of our K-selected galaxies are not forming stars intensely. This high fraction of galaxies without detected line emission and low SFRs may imply that the suppression of star formation in massive galaxies occurs at higher redshift than is predicted by current CDM galaxy formation models. However, obscured star formation may have been missed, and deep mid-infrared imaging is needed to clarify this situation.

(Abridged) We present the evolution of the volume averaged properties of the rest-frame optically luminous galaxy population to z~3, determined from four disjoint deep fields with optical to near-infrared wavelength coverage. We select galaxies above a rest-frame V-band luminosity of 3x10^10 Lsol and characterize their rest-frame UV through optical properties via the mean spectral energy distribution (SED). To measure evolution we apply the same selection criteria to a sample of galaxies from the Sloan Digital Sky Survey and COMBO-17. The mean rest-frame 2200Ang through V-band SED becomes steadily bluer with increasing redshift but at z<3 the mean SED falls within the range defined by ``normal'' galaxies in the nearby Universe. We measure stellar mass-to-light ratios (Mstar/L) by fitting models to the rest-frame UV-optical SEDs and derive the stellar mass density. The stellar mass density in luminous galaxies has increased by a factor of 3.5-7.9 from z=3 to z=0.1, including field-to-field variance uncertainties. After correcting to total, the measured mass densities at z<2 lie below the integral of the star formation rate (SFR) density as a function of redshift as derived from UV selected samples. This may indicate a systematic error in the mass densities or SFR(z) estimates. We find large discrepancies between recent model predictions for the evolution of the mass density and our results, even when our observational selection is applied to the models. Finally we determine that Distant Red Galaxies (selected to have J_s - K_s>2.3) in our LV selected samples contribute 30% and 64% of the stellar mass budget at z~2 and z~ 2.8 respectively. These galaxies are largely absent from UV surveys and this result highlights the need for mass selection of high redshift galaxies.

We use near-infrared (NIR) spectroscopy obtained with GNIRS on Gemini, NIRSPEC on KECK, and ISAAC on the VLT to study the rest-frame optical continua of three `Distant Red Galaxies' (having Js - Ks > 2.3) at z>2. All three galaxy spectra show the Balmer/4000 Angstrom break in the rest-frame optical. The spectra allow us to determine spectroscopic redshifts from the continuum with an estimated accuracy dz/(1+z) ~ 0.001-0.04. These redshifts agree well with the emission line redshifts for the 2 galaxies with Halpha emission. This technique is particularly important for galaxies that are faint in the rest-frame UV, as they are underrepresented in high redshift samples selected in optical surveys and are too faint for optical spectroscopy. Furthermore, we use the break, continuum shape, and equivalent width of Halpha together with evolutionary synthesis models to constrain the age, star formation timescale, dust content, stellar mass and star formation rate of the galaxies. Inclusion of the NIR spectra in the stellar population fits greatly reduces the range of possible solutions for stellar population properties. We find that the stellar populations differ greatly among the three galaxies, ranging from a young dusty starburst with a small break and strong emission lines to an evolved galaxy with a strong break and no detected line emission. The dusty starburst galaxy has an age of 0.3 Gyr and a stellar mass of 1*10^11 Msun. The spectra of the two most evolved galaxies imply ages of 1.3-1.4 Gyr and stellar masses of 4*10^11 Msun. The large range of properties seen in these galaxies strengthens our previous much more uncertain results from broadband photometry. Larger samples are required to determine the relative frequency of dusty starbursts and (nearly) passively evolving galaxies at z~2.5.

Using the deep multi-wavelength MUSYC, GOODS, and FIRES surveys we construct a stellar mass-limited sample of galaxies at 2<z<3. The sample comprises 294 galaxies with M>10^11 Solar masses distributed over four independent fields with a total area of almost 400 sq arcmin. The mean number density of massive galaxies in this redshift range is (2.2+-0.6) x 10^-4 Mpc^-3. We present median values and 25th and 75th percentiles for the distributions of observed R mags, observed J-K colors, and rest-frame UV continuum slopes, M/L(V) ratios, and U-V colors. The galaxies show a large range in all these properties. The ``median galaxy'' is faint in the observer's optical (R=25.9), red in the observed near-IR (J-K=2.48), has a rest-frame UV spectrum which is relatively flat (beta=-0.4), and rest-frame optical colors resembling those of nearby spiral galaxies (U-V=0.62). We determine which galaxies would be selected as Lyman break galaxies (LBGs) or Distant Red Galaxies (DRGs, having J-K>2.3) in this mass-limited sample. By number DRGs make up 69% of the sample and LBGs 20%, with a small amount of overlap. By mass DRGs make up 77% and LBGs 17%. Neither technique provides a representative sample of massive galaxies at 2<z<3 as they only sample the extremes of the population. As we show here, multi-wavelength surveys with high quality photometry are essential for an unbiased census of massive galaxies in the early Universe. The main uncertainty in this analysis is our reliance on photometric redshifts; confirmation of the results presented here requires extensive near-infrared spectroscopy of optically-faint samples.

We take advantage of the Hubble Ultra Deep Field (UDF) data to study the restframe optical and ultra violet (UV) morphologies of the novel population of Distant Red Galaxies (DRGs). Six galaxies with J-Ks > 2.3 are found to Ks=21.5, five of which have photometric redshifts z_phot > 2, corresponding to a surface density of 0.9/arcmin^2. The surface brightness distributions of the z_phot > 2 galaxies are better represented by exponential disks than R^{1/4}-laws. Two of the z_phot > 2 galaxies are extended, while three have compact morphologies. The restframe optical morphology of the z_phot > 2 galaxies is quite different from the restframe UV morphology: all the galaxies have red central components which dominate in the NICMOS H_{160}-band images, and distinct off-center blue features which show up in (and often dominate) the ACS images. The mean measured effective radius of the z_phot > 2 galaxies is <r_e> =1.9+/-1.4 kpc, similar (within the errors) to the mean size of LBGs at similar redshifts. All the DRGs are resolved in the ACS images, while four are resolved in the NICMOS images. Two of the z_phot > 2 galaxies are bright X-ray sources and hence host AGN. The diverse restframe optical and UV morphological properties of DRGs derived here suggest that they have complex stellar populations, consisting of both evolved populations that dominate the mass and the restframe optical light, and younger populations, which show up as patches of star formation in the restframe UV light; in many ways resembling the properties of normal local galaxies. This interpretation is supported by fits to the broadband SEDs, which for all five z_phot > 2 are best represented by models with extended star formation histories and substantial amounts of dust.

Recent studies have shown that K-luminous galaxies at 2<z<2.5 have high Halpha line widths and high [NII]/Halpha ratios. If these lines arise from photo-ionization by massive stars in HII regions these results imply that massive, metal rich galaxies exist at high redshift. Here we investigate the ionization mechanism in a galaxy with Ks=19.1 and z=2.225 in the Chandra Deep Field South, using the new Gemini Near Infrared Spectrograph (GNIRS). GNIRS' cross-dispersed mode gives simultaneous access to the entire 1-2.5 micron wavelength range, allowing accurate measurements of line ratios of distant galaxies. From the ratio of Halpha/Hbeta we infer that the line emitting gas is heavily obscured, with E(B-V)=0.8. The reddening is higher than that inferred previously from the UV-optical continuum emission, consistent with findings for nearby star burst galaxies. We find that the galaxy has Seyfert-like line ratios, [NII]/Halpha ~ 0.6 and [OIII]/Hbeta ~ 6, which can be caused by photo-ionization by an active galactic nucleus (AGN) or shock ionization due to a strong galactic wind. Although we cannot exclude the presence of an AGN, the lack of AGN spectral features in the rest-frame ultraviolet, the consistency of radio, X-ray, and rest-frame UV star formation indicators, the fact that the [OIII]/Hbeta ratio remains high out to ~10 kpc from the nucleus, and the observed gas kinematics all argue for the wind hypothesis. Similar shock-induced ionization is seen in nearby star burst galaxies with strong winds. The evidence for shock ionization implies that measurements of metallicities and dynamical masses of star forming z>2 galaxies should be regarded with caution, especially since the existence of strong galactic winds in these objects is well established. [ABRIDGED]

We investigate the nature of the substantial population of high-z galaxies with Js-Ks>2.3 discovered as part of our FIRES survey. This colour cut efficiently isolates z>2 galaxies with red rest-frame optical colors ("Distant Red Galaxies" or DRGs). We select objects in the 2.5'x2.5' HDF-South (HDF-S) and 5'x5' field around the MS1054-03 cluster; the surface densities at Ks<21 are 1.6+-0.6 and 1.0+-0.2 arcmin^-2. We discuss the 34 DRGs at 2<z<3.5: 11 at Ks<22.5 in HDF-S and 23 at Ks<21.7 in the MS1054-03 field. We analyze the SEDs constructed from our deep near-infrared (NIR) and optical imaging from the ESO VLT and HST. We develop diagnostics involving I-Js, Js-H, and H-Ks to argue that the red NIR colors of DRGs cannot be attributed solely to extinction and require for many an evolved stellar population with prominent Balmer/4000A break. In the rest-frame, the optical colours of DRGs fall within the envelope of normal nearby galaxies and the UV colours suggest a wide range in star formation activity and/or extinction. This contrasts with the much bluer and more uniform SEDs of Lyman break galaxies (LBGs). From evolutionary synthesis models with constant star formation, solar metallicity, Salpeter IMF, and Calzetti et al. extinction law, we derive for the HDF-S (MS1054-03 field) DRGs median ages of 1.7(2.0) Gyr, A_V = 2.7(2.4) mag, stellar masses 0.8(1.6)x10^11 Msun, M/L_V = 1.2(2.3) Msun/LVsun, and SFR = 120(170) Msun/yr. Models assuming declining SFRs with e-folding timescales of 10Myr-1Gyr generally imply younger ages, lower A_V's and SFRs, but similar stellar masses within a factor of two. Compared to LBGs at similar redshifts and rest-frame L_V's, DRGs are older, more massive, and more obscured for any given star formation history. [ABRIDGED]