For studies of galaxy formation and evolution, one of the major benefits of the James Webb Space Telescope is that space-based IFUs like those on its NIRSpec and MIRI instruments will enable spatially resolved spectroscopy of distant galaxies, including spectroscopy at the scale of individual star-forming regions in galaxies that have been gravitationally lensed. In the meantime, there is only a very small subset of lensed sources where work like this is possible even with the Hubble Space Telescope's Wide Field Camera 3 infrared channel grisms. We examine two of these sources, SDSS J1723+3411 and SDSS J2340+2947, using HST WFC3/IR grism data and supporting spatially-unresolved spectroscopy from several ground-based instruments to explore the size of spatial variations in observed strong emission line ratios like O32, R23, which are sensitive to ionization parameter and metallicity, and the Balmer decrement as an indicator of reddening. We find significant spatial variation in the reddening and the reddening-corrected O32 and R23 values which correspond to spreads of a few tenths of a dex in ionization parameter and metallicity. We also find clear evidence of a negative radial gradient in star formation in SDSS J2340+2947 and tentative evidence of one in SDSS J1723+3411, though its star formation is quite asymmetric. Finally, we find that reddening can vary enough spatially to make spatially-resolved reddening corrections necessary in order to characterize gradients in line ratios and the physical conditions inferred from them, necessitating the use of space-based IFUs for future work on larger, more statistically robust samples.

The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

We investigate a large sample of massive galaxies at $z\sim1$ with combined $HST$ broad-band and grism observations to constrain the star-formation histories of these systems as they transition from a star-forming state to quiescence. Among our sample of massive $(M_*>10^{10}~{\rm M_\odot})$ galaxies at $0.7<z<1.2$, dust-corrected H$\alpha$ and UV star-formation indicators agree with a small dispersion ($\sim0.2$~dex) for galaxies on the main sequence, but diverge and exhibit substantial scatter ($\sim0.7$~dex) once they drop significantly below the star-forming main sequence. Significant H$\alpha$ emission is present in galaxies with low dust-corrected UV SFR values as well as galaxies classified as quiescent using the $UVJ$ diagram. We compare the observed H$\alpha$ flux distribution to the expected distribution assuming bursty or smooth star-formation histories, and find that massive galaxies at $z\sim1$ are most consistent with a quick, bursty quenching process. This suggests that mechanisms such as feedback, stochastic gas flows, and minor mergers continue to induce low-level bursty star formation in massive galaxies at moderate redshift, even as they quench.

In the past decade, our understanding of how stars and galaxies formed during the first 5 billion years after the Big Bang has been revolutionized by observations that leverage gravitational lensing by intervening masses, which act as natural cosmic telescopes to magnify background sources. Previous studies have harnessed this effect to probe the distant universe at ultraviolet, optical, infrared and millimeter wavelengths. However, strong lensing studies of young, star-forming galaxies have never extended into X-ray wavelengths, which uniquely trace high-energy phenomena. Here we report an X-ray detection of star formation in a highly magnified, strongly lensed galaxy. This lensed galaxy, seen during the first third of the history of the Universe, is a low--mass, low--metallicity starburst with elevated X-ray emission, and is a likely analog to the first generation of galaxies. Our measurements yield insight into the role that X-ray emission from stellar populations in the first generation of galaxies may play in re-ionizing the Universe. This observation paves the way for future strong lensing-assisted X-ray studies of distant galaxies reaching orders of magnitude below the detection limits of current deep fields, and previews the depths that will be attainable with future X-ray observatories.

We report the serendipitous detection of two 3 mm continuum sources found in deep ALMA Band 3 observations to study intermediate redshift galaxies in the COSMOS field. One is near a foreground galaxy at 1.3", but is a previously unknown dust-obscured star-forming galaxy (DSFG) at probable $z_{CO}=3.329$, illustrating the risk of misidentifying shorter wavelength counterparts. The optical-to-mm spectral energy distribution (SED) favors a grey $\lambda^{-0.4}$ attenuation curve and results in significantly larger stellar mass and SFR compared to a Calzetti starburst law, suggesting caution when relating progenitors and descendants based on these quantities. The other source is missing from all previous optical/near-infrared/sub-mm/radio catalogs ("ALMA-only"), and remains undetected even in stacked ultradeep optical ($>29.6$ AB) and near-infrared ($>27.9$ AB) images. Using the ALMA position as a prior reveals faint $SNR\sim3$ measurements in stacked IRAC 3.6+4.5, ultradeep SCUBA2 850$\mu$m, and VLA 3GHz, indicating the source is real. The SED is robustly reproduced by a massive $M^*=10^{10.8}$M$_\odot$ and $M_{gas}=10^{11}$M$_\odot$, highly obscured $A_V\sim4$, star forming $SFR\sim300$ M$_{\odot}$yr$^{-1}$ galaxy at redshift $z=5.5\pm$1.1. The ultrasmall 8 arcmin$^{2}$ survey area implies a large yet uncertain contribution to the cosmic star formation rate density CSFRD(z=5) $\sim0.9\times10^{-2}$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$, comparable to all ultraviolet-selected galaxies combined. These results indicate the existence of a prominent population of DSFGs at $z>4$, below the typical detection limit of bright galaxies found in single-dish sub-mm surveys, but with larger space densities $\sim3 \times 10^{-5}$ Mpc$^{-3}$, higher duty cycles $50-100\%$, contributing more to the CSFRD, and potentially dominating the high-mass galaxy stellar mass function.

This manuscript describes the public release of the Hubble Legacy Fields (HLF) project photometric catalog for the extended GOODS-South region from the Hubble Space Telescope (HST) archival program AR-13252. The analysis is based on the version 2.0 HLF data release that now includes all ultraviolet (UV) imaging, combining three major UV surveys. The HLF data combines over a decade worth of 7475 exposures taken in 2635 orbits totaling 6.3 Msec with the HST Advanced Camera for Surveys Wide Field Channel (ACS/WFC) and the Wide Field Camera 3 (WFC3) UVIS/IR Channels in the greater GOODS-S extragalactic field, covering all major observational efforts (e.g., GOODS, GEMS, CANDELS, ERS, UVUDF and many other programs; see Illingworth et al 2019, in prep). The HLF GOODS-S catalogs include photometry in 13 bandpasses from the UV (WFC3/UVIS F225W, F275W and F336W filters), optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters), to near-infrared (WFC3/IR F098M, F105W, F125W, F140W and F160W filters). Such a data set makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range from high resolution mosaics that are largely contiguous. Here, we describe a photometric analysis of 186,474 objects in the HST imaging at wavelengths 0.2--1.6$\mu$m. We detect objects from an ultra-deep image combining the PSF-homogenized and noise-equalized F850LP, F125W, F140W and F160W images, including Gaia astrometric corrections. SEDs were determined by carefully taking the effects of the point-spread function in each observation into account. All of the data presented herein are available through the HLF website (https://archive.stsci.edu/prepds/hlf/).

We present high spatial resolution imaging of the CO(1-0) line from the Karl G. Jansky Very Large Array (VLA) of COSMOS27289, a massive, compact star forming galaxy at z=2.234. This galaxy was selected to be structurally similar to z~2 passive galaxies. Our previous observations showed that it is very gas-poor with respect to typical star-forming galaxies at these redshifts, consistent with a rapid transition to quiescence as the molecular gas is depleted. The new data show that both the molecular gas fraction f_H2 = M_H2/M_star and the molecular gas depletion time t_dep = M_H2/SFR are lower in the central 1-2kpc of the galaxy and rise at larger radii ~2-4kpc. These observations are consistent with a scenario in which COSMOS27289 will imminently cease star formation in the inner regions before the outskirts, i.e. inside-out quenching, the first time this phenomenon has been seen via observations of molecular gas in the high-redshift universe. We find good qualitative and quantitative agreement with a hydrodynamical simulation of galaxy quenching, in which the central suppression of molecular gas arises due to rapid gas consumption and outflows that evacuate the central regions of gas. Our results provide independent evidence for inside-out quenching of star formation as a plausible formation mechanism for z~2 quiescent galaxies.

Galaxy observations are influenced by many physical parameters: stellar masses, star formation rates (SFRs), star formation histories (SFHs), metallicities, dust, black hole activity, and more. As a result, inferring accurate physical parameters requires high-dimensional models which capture or marginalize over this complexity. Here we re-assess inferences of galaxy stellar masses and SFRs using the 14-parameter physical model Prospector-$\alpha$ built in the Prospector Bayesian inference framework. We fit the photometry of 58,461 galaxies from the 3D-HST catalogs at $0.5 < z < 2.5$. The resulting stellar masses are $\sim0.1-0.3$ dex larger than the fiducial masses while remaining consistent with dynamical constraints. This change is primarily due to the systematically older SFHs inferred with Prospector. The SFRs are $\sim0.1-1+$ dex lower than UV+IR SFRs, with the largest offsets caused by emission from "old" ($t>100$ Myr) stars. These new inferences lower the observed cosmic star formation rate density by $\sim0.2$ dex and increase the observed stellar mass growth by $\sim 0.1$ dex, finally bringing these two quantities into agreement and implying an older, more quiescent Universe than found by previous studies at these redshifts. We corroborate these results by showing that the Prospector-$\alpha$ SFHs are both more physically realistic and are much better predictors of the evolution of the stellar mass function. Finally, we highlight examples of observational data which can break degeneracies in the current model; these observations can be incorporated into priors in future models to produce new & more accurate physical parameters.

We present strong gravitational lensing models for 37 galaxy clusters from the SDSS Giant Arcs Survey. We combine data from multi-band Hubble Space Telescope WFC3imaging, with ground-based imaging and spectroscopy from Magellan, Gemini, APO, and MMT, in order to detect and spectroscopically confirm new multiply-lensed background sources behind the clusters. We report spectroscopic or photometric redshifts of sources in these fields, including cluster galaxies and background sources. Based on all available lensing evidence, we construct and present strong lensing mass models for these galaxy clusters.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

We present $Hubble\ Space\ Telescope$ ($HST$) imaging and grism spectroscopy of a strongly lensed LIRG at $z=0.816$, SGAS 143845.1$+$145407, and use the magnification boost of gravitational lensing to study the distribution of star formation throughout this galaxy. Based on the $HST$ imaging data, we create a lens model for this system; we compute the mass distribution and magnification map of the $z=0.237$ foreground lens. We find that the magnification of the lensed galaxy ranges between $2$ and $10$, with a total magnification (measured over all the images of the source) of $\mu=11.8^{+4.6}_{-2.4}$. We find that the total projected mass density within $\sim34$ kpc of the brightest cluster galaxy is $6.0^{+0.3}_{-0.7}\times10^{12}\,M_{\odot}$. Using the lens model we create a source reconstruction for SGAS 143845.1$+$145407, which paired with a faint detection of H$\alpha$ in the grism spectroscopy, allows us to finally comment directly on the distribution of star formation in a $z\sim1$ LIRG. We find widespread star formation across this galaxy, in agreement with the current understanding of these objects. However, we note a deficit of H$\alpha$ emission in the nucleus of SGAS 143845.1$+$145407, likely due to dust extinction.

One of the greatest challenges to theoretical models of massive galaxy formation is the regulation of star formation at early times. The relative roles of molecular gas expulsion, depletion, and stabilization are uncertain as direct observational constraints of the gas reservoirs in quenched or quenching galaxies at high redshift are scant. We present ALMA observations of CO(2-1) in a massive ($\log M_{\star}/M_{\odot}=11.2$), recently quenched galaxy at $z=1.522$. The optical spectrum of this object shows strong Balmer absorption lines, which implies that star formation ceased $\sim$0.8 Gyr ago. We do not detect CO(2-1) line emission, placing an upper limit on the molecular $\mathrm{H_2}$ gas mass of 1.1$\times10^{10}\,M_{\odot}$. The implied gas fraction is $f_{\rm{H_2}}{\equiv M_{H_2}/M_{\star}}<7\%$, $\sim10\times$ lower than typical star forming galaxies at similar stellar masses at this redshift, among the lowest gas fractions at this specific star formation rate at any epoch, and the most stringent constraint on the gas contents of a $z>1$ passive galaxy to date. Our observations show that the depletion of $\mathrm{H_2}$ from the interstellar medium of quenched objects can be both efficient and fairly complete, in contrast to recent claims of significant cold gas in recently quenched galaxies. We explore the variation in observed gas fractions in high-$z$ galaxies and show that galaxies with high stellar surface density have low $f_{\rm{H_2}}$, similar to recent correlations between specific star formation rate and stellar surface density.

We investigate the median flattening of galaxies at $0.2<z<4.0$ in all five CANDELS/3D-HST fields via the apparent axis ratio $q$. We separate the sample into bins of redshift, stellar-mass, s\'ersic index, size, and UVJ determined star-forming state to discover the most important drivers of the median $q$ ($q_{med}$). Quiescent galaxies at $z<1$ and $M_{*}>10^{11}M_{\odot}$ are rounder than those at lower masses, consistent with the hypothesis that they have grown significantly through dry merging. The massive quiescent galaxies at higher redshift become flatter, and are as flat as star forming massive galaxies at $2.5<z<3.5$, consistent with formation through direct transformations or wet mergers. We find that in quiescent galaxies, correlations with $q_{med}$ and $M_{*}$, $z$ and $r_{e}$ are driven by the evolution in the s\'ersic index ($n$), consistent with the growing accumulation of minor mergers at lower redshift. Interestingly, $n$ does not drive these trends fully in star-forming galaxies. Instead, the strongest predictor of $q$ in star-forming galaxies is the effective radius, where larger galaxies are flatter. Our findings suggest that $q_{med}$ is tracing bulge-to-total ($B/T$) galaxy ratio which would explain why smaller/more massive star-forming galaxies are rounder than their extended/less massive analogues, although it is unclear why s\'ersic index correlates more weakly with flattening for star forming galaxies than for quiescent galaxies.

The IRAC mapping of the NMBS-II fields program is an imaging survey at 3.6 and 4.5$\mu$m with the Spitzer Infrared Array Camera (IRAC). The observations cover three Canada-France-Hawaii Telescope Legacy Survey Deep (CFHTLS-D) fields, including one also imaged by AEGIS, and two MUSYC fields. These are then combined with archival data from all previous programs into deep mosaics. The resulting imaging covers a combined area of about 3 $deg^2$, with at least $\sim$2 hr integration time for each field. In this work, we present our data reduction techniques and document the resulting coverage maps at 3.6 and 4.5$\mu$m. All of the images are W-registered to the reference image, which is either the z-band stack image of the 25\% best seeing images from the CFHTLS-D for CFHTLS-D1, CFHTLS-D3, and CFHTLS-D4, or the K-band images obtained at the Blanco 4-m telescope at CTIO for MUSYC1030 and MUSYC1255. We make all images and coverage maps described herein publicly available via the Spitzer Science Center.

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 present a targeted follow-up Hubble Space Telescope WFC3 F160W imaging study of very massive galaxies $(\log(M_{\rm{star}}/M_{\odot})> 11.2)$ selected from a combination of ground-based near-infrared galaxy surveys (UltraVISTA, NMBS-II, UKIDSS UDS) at $1.5<z<3$. We find that these galaxies are diverse in their structures, with $\sim1/3$ of the targets being composed of close pairs, and span a wide range in sizes. At $1.5<z<2.5$, the sizes of both star-forming and quiescent galaxies are consistent with the extrapolation of the stellar mass-size relations determined at lower stellar masses. At $2.5<z<3.0$, however, we find evidence that quiescent galaxies are systematically larger than expected based on the extrapolation of the relation derived using lower stellar mass galaxies. We used the observed light profiles of the blended systems to decompose their stellar masses and investigate the effect of the close pairs on the measured number densities of very massive galaxies in the early universe. We estimate correction factors to account for close-pair blends and apply them to the observed stellar mass functions measured using ground-based surveys. Given the large uncertainties associated with this extreme population of galaxies, there is currently little tension between the (blending-corrected) number density estimates and predictions from theoretical models. Although we currently lack the statistics to robustly correct for close-pair blends, we show that this is a systematic effect which can reduce the observed number density of very massive galaxies by up to a factor of $\sim1.5$, and should be accounted for in future studies of stellar mass functions.

A decade of study has established that the molecular gas properties of star-forming galaxies follow coherent scaling relations out to z~3, suggesting remarkable regularity of the interplay between molecular gas, star formation, and stellar growth. Passive galaxies, however, are expected to be gas-poor and therefore faint, and thus little is known about molecular gas in passive galaxies beyond the local universe. Here we present deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2-1) emission in 8 massive (Mstar ~ 10^11 Msol) galaxies at z~0.7 selected to lie a factor of 3-10 below the star-forming sequence at this redshift, drawn from the Large Early Galaxy Astrophysics Census (LEGA-C) survey. We significantly detect half the sample, finding molecular gas fractions <~0.1. We show that the molecular and stellar rotational axes are broadly consistent, arguing that the molecular gas was not accreted after the galaxies became quiescent. We find that scaling relations extrapolated from the star-forming population over-predict both the gas fraction and gas depletion time for passive objects, suggesting the existence of either a break or large increase in scatter in these relations at low specific star formation rate. Finally, we show that the gas fractions of the passive galaxies we have observed at intermediate redshifts are naturally consistent with evolution into local massive early-type galaxies by continued low-level star formation, with no need for further gas accretion or dynamical stabilization of the gas reservoirs in the intervening 6 billion years.

We present Hubble multi-wavelength photometric catalogs, including (up to) 17 filters with the Advanced Camera for Surveys and Wide Field Camera 3 from the ultra-violet to near-infrared for the Hubble Frontier Fields and associated parallels. We have constructed homogeneous photometric catalogs for all six clusters and their parallels. To further expand these data catalogs, we have added ultra-deep $K_{S}$-band imaging at 2.2~\micron\ from the Very Large Telescope HAWK-I and Keck-I MOSFIRE instruments. We also add post-cryogenic \spitzer\ imaging at 3.6~\micron\ and 4.5~\micron\ with the Infrared Array Camera (IRAC), as well as archival IRAC 5.8~\micron\ and 8.0~\micron\ imaging when available. We introduce the public release of the multi-wavelength (0.2--8~\micron) photometric catalogs, and we describe the unique steps applied for the construction of these catalogs. Particular emphasis is given to the source detection band, the contamination of light from the bright cluster galaxies and intra-cluster light. In addition to the photometric catalogs, we provide catalogs of photometric redshifts and stellar population properties. Furthermore, this includes all the images used in the construction of the catalogs, including the combined models of bright cluster galaxies and intra-cluster light, the residual images, segmentation maps and more. These catalogs are a robust data set of the Hubble Frontier Fields and will be an important aide in designing future surveys, as well as planning follow-up programs with current and future observatories to answer key questions remaining about first light, reionization, the assembly of galaxies and many more topics, most notably, by identifying high-redshift sources to target.

The multiplexing capability of slitless spectroscopy is a powerful asset in creating large spectroscopic datasets, but issues such as spectral confusion make the interpretation of the data challenging. Here we present a new method to search for emission lines in the slitless spectroscopic data from the 3D-HST survey utilizing the Wide-Field Camera 3 on board the Hubble Space Telescope. Using a novel statistical technique, we can detect compact (extended) emission lines at 90% completeness down to fluxes of 1.5 (3.0) times 10^{-17} erg/s/cm^2, close to the noise level of the grism exposures, for objects detected in the deep ancillary photometric data. Unlike previous methods, the Bayesian nature allows for probabilistic line identifications, namely redshift estimates, based on secondary emission line detections and/or photometric redshift priors. As a first application, we measure the comoving number density of Extreme Emission Line Galaxies (restframe [O III] 5007 equivalent widths in excess of 500 Angstroms). We find that these galaxies are nearly 10 times more common above z~1.5 than at z<0.5. With upcoming large grism surveys such as Euclid and WFIRST as well as grisms featuring prominently on the NIRISS and NIRCam instruments on James Webb Space Telescope, methods like the one presented here will be crucial for constructing emission line redshift catalogs in an automated and well-understood manner.

We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various star formation history for each galaxy. Like others, we find that the slope of the Main Sequence (MS) of star formation in the $(M_*;SFR)$ plane bends at high mass. We observe clear morphological differences among galaxies across the MS, which also correlate with stellar age. At all redshifts, galaxies that are quenching or quenched, and thus old, have high $\Sigma_1$ (the projected density within the central 1 kpc), while younger, star-forming galaxies span a much broader range of $\Sigma_1$, which includes the high values observed for quenched galaxies, but also extends to much lower values. As galaxies age and quench, the stellar age and the dispersion of $\Sigma_1$ for fixed values of $M_{*}$ shows two different regimes, one, at the low--mass end, where quenching might be driven by causes external to the galaxies; the other, at the high--mass end, where quenching is driven by internal causes, very likely the mass given the low scatter of $\Sigma_1$ (mass quenching). We suggest that the monotonic increase of central density as galaxies grow is one manifestation of a more general phenomenon of structural transformation that galaxies undergo as they evolve.

The total star formation budget of galaxies consists of the sum of the unobscured star formation, as observed in the rest-frame ultraviolet (UV), together with the obscured component that is absorbed and re-radiated by dust grains in the infrared. We explore how the fraction of obscured star formation depends on stellar mass for mass-complete samples of galaxies at $0<z<2.5$. We combine GALEX and WISE photometry for SDSS-selected galaxies with the 3D-HST treasury program and Spitzer/MIPS 24$\mu$m photometry in the well-studied 5 extragalactic CANDELS fields. We find a strong dependence of the fraction of obscured star formation ($f_{\mathrm{obscured}}$=SFR$_{\mathrm{IR}}$/SFR$_{\mathrm{UV+IR}}$) on stellar mass, with remarkably little evolution in this fraction with redshift out to $z$=2.5. 50\% of star formation is obscured for galaxies with log(M/M$_{\odot}$)=9.4; although unobscured star formation dominates the budget at lower masses, there exists a tail of low mass extremely obscured star-forming galaxies at $z>1$. For log(M/M$_{\odot}$)$>$10.5, $>$90\% of star formation is obscured at all redshifts. We also show that at fixed total SFR, $f_{\mathrm{obscured}}$ is lower at higher redshift. At fixed mass, high-redshift galaxies are observed to have more compact sizes and much higher star formation rates, gas fractions and hence surface densities (implying higher dust obscuration), yet we observe no redshift evolution in $f_{\mathrm{obscured}}$ with stellar mass. This poses a challenge to theoretical models to reproduce, where the observed compact sizes at high redshift seem in tension with lower dust obscuration.

Using the combined resolving power of the Hubble Space Telescope and gravitational lensing, we resolve star-forming structures in a z~2.5 galaxy on scales much smaller than the usual kiloparsec diffraction limit of HST. SGAS J111020.0+645950.8 is a clumpy, star forming galaxy lensed by the galaxy cluster SDSS J1110+6459 at z = 0.659, with a total magnification ~30x across the entire arc. We use a hybrid parametric/non-parametric strong lensing mass model to compute the deflection and magnification of this giant arc, reconstruct the light distribution of the lensed galaxy in the source plane, and resolve the star formation into two dozen clumps. We develop a forward-modeling technique to model each clump in the source plane. We ray trace the model to the image plane, convolve with the instrumental point spread function (PSF), and compare with the GALFIT model of the clumps in the image plane, which decomposes clump structure from more extended emission. This technique has the advantage, over ray tracing, by accounting for the asymmetric lensing shear of the galaxy in the image plane and the instrument PSF. At this resolution, we can begin to study star formation on a clump-by-clump basis, toward the goal of understanding feedback mechanisms and the buildup of exponential disks at high redshift.

We present measurements of the surface density of star formation, the star-forming clump luminosity function, and the clump size distribution function, for the lensed galaxy SGAS J111020.0+645950.8 at a redshift of z=2.481. The physical size scales that we probe, radii r=30--50 pc, are considerably smaller scales than have yet been studied at these redshifts. The star formation surface density we find within these small clumps is consistent with surface densities measured previously for other lensed galaxies at similar redshift. Twenty-two percent of the rest-frame ultraviolet light in this lensed galaxy arises from small clumps, with r<100 pc. Within the range of overlap, the clump luminosity function measured for this lensed galaxy is remarkably similar to those of z~0 galaxies. In this galaxy, star-forming regions smaller than 100 pc---physical scales not usually resolved at these redshifts by current telescopes---are important locations of star formation in the distant universe. If this galaxy is representative, this may contradict the theoretical picture in which the critical size scale for star formation in the distant universe is of order 1 kiloparsec. Instead, our results suggest that current telescopes have not yet resolved the critical size scales of star-forming activity in galaxies over most of cosmic time.

For lensed galaxy SGAS J111020.0+645950.8 at redshift z=2.481, which is magnified by a factor of 28 +- 8, we analyze the morphology of star formation as traced by rest-frame ultraviolet emission, in both the highly-magnified source plane, and in simulations of how this galaxy would appear without lensing magnification. Were this galaxy not lensed but drawn from an HST deep field, we would conclude that almost all its star formation arises from an exponential disk (S\'ersic index of 1.0 +- 0.4) with an effective radius of r_e = 2.7 +- 0.3 kpc measured from two-dimensional fitting to F606W using Galfit, and r_e=1.9 +- 0.1 kpc measured by fitting a radial profile to F606W elliptical isophotes. At the normal spatial resolution of the deep fields, there is no sign of clumpy star formation within SGAS J111020.0+645950.8 . However, the enhanced spatial resolution enabled by gravitational lensing tells a very different story: much of the star formation arises in two dozen clumps with sizes of r=30--50 pc spread across the 7 kpc length of the galaxy. The color and spatial distribution of the diffuse component suggests that still smaller clumps are unresolved. Despite this clumpy, messy morphology, the radial profile is still well-characterized by an exponential profile. In this lensed galaxy, stars are forming in complexes with sizes well below 100 pc; such sizes are wholly unexplored by surveys of galaxy evolution at 1<z<3.

We have submitted to MAST the 1.5 version data release of the Hubble Legacy Fields (HLF) project covering a 25 x 25 arcmin area over the GOODS-S (ECDF-S) region from the HST archival program AR-13252. The release combines exposures from Hubble's two main cameras, the Advanced Camera for Surveys (ACS/WFC) and the Wide Field Camera 3 (WFC3/IR), taken over more than a decade between mid-2002 to the end of 2016. The HLF includes essentially all optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters) and infrared (WFC3/ IR F098M, F105W, F125W, F140W and F160W filters) data taken by Hubble over the original CDF-S region including the GOODS-S, ERS, CANDELS and many other programs (31 in total). The data has been released at https://archive.stsci.edu/prepds/hlf/ as images with a common astrometric reference frame, with corresponding inverse variance weight maps. We provide one image per filter of WFC3/IR images at 60 mas per pixel resolution and two ACS/WFC images per filter, at both 30 and 60 mas per pixel. Since this comprehensive dataset combines data from 31 programs on the GOODS-S/CDF-S, the AR proposal identified the MAST products by the global name "Hubble Legacy Field", with this region being identified by "HLF-GOODS-S". This dataset complements that of the Frontier Fields program. The total incorporated in the HLF-GOODS-S is 5.8 Msec in 7211 exposures from 2442 orbits. This is ~70% of a HST full cycle!

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

In this paper, we investigate the relationship between star formation and structure, using a mass-complete sample of 27,893 galaxies at $0.5<z<2.5$ selected from 3D-HST. We confirm that star-forming galaxies are larger than quiescent galaxies at fixed stellar mass (M$_{\star}$). However, in contrast with some simulations, there is only a weak relation between star formation rate (SFR) and size within the star-forming population: when dividing into quartiles based on residual offsets in SFR, we find that the sizes of star-forming galaxies in the lowest quartile are 0.27$\pm$0.06 dex smaller than the highest quartile. We show that 50% of star formation in galaxies at fixed M$_{\star}$ takes place within a narrow range of sizes (0.26 dex). Taken together, these results suggest that there is an abrupt cessation of star formation after galaxies attain particular structural properties. Confirming earlier results, we find that central stellar density within a 1 kpc fixed physical radius is the key parameter connecting galaxy morphology and star formation histories: galaxies with high central densities are red and have increasingly lower SFR/M$_{\star}$, whereas galaxies with low central densities are blue and have a roughly constant (higher) SFR/M$_{\star}$ at a given redshift. We find remarkably little scatter in the average trends and a strong evolution of $>$0.5 dex in the central density threshold correlated with quiescence from $z\sim0.7-2.0$. Neither a compact size nor high-$n$ are sufficient to assess the likelihood of quiescence for the average galaxy; rather, the combination of these two parameters together with M$_{\star}$ results in a unique quenching threshold in central density/velocity.

In this paper, we use stacking analysis to trace the mass-growth, colour evolution, and structural evolution of present-day massive galaxies ($\log(M_{*}/M_{\odot})=11.5$) out to $z=5$. We utilize the exceptional depth and area of the latest UltraVISTA data release, combined with the depth and unparalleled seeing of CANDELS to gather a large, mass-selected sample of galaxies in the NIR (rest-frame optical to UV). Progenitors of present-day massive galaxies are identified via an evolving cumulative number density selection, which accounts for the effects of merging to correct for the systematic biases introduced using a fixed cumulative number density selection, and find progenitors grow in stellar mass by $\approx1.5~\mathrm{dex}$ since $z=5$. Using stacking, we analyze the structural parameters of the progenitors and find that most of the stellar mass content in the central regions was in place by $z\sim2$, and while galaxies continue to assemble mass at all radii, the outskirts experience the largest fractional increase in stellar mass. However, we find evidence of significant stellar mass build up at $r<3~\mathrm{kpc}$ beyond $z>4$ probing an era of significant mass assembly in the interiors of present day massive galaxies. We also compare mass assembly from progenitors in this study to the EAGLE simulation and find qualitatively similar assembly with $z$ at $r<3~\mathrm{kpc}$. We identify $z\sim1.5$ as a distinct epoch in the evolution of massive galaxies where progenitors transitioned from growing in mass and size primarily through in-situ star formation in disks to a period of efficient growth in $r_{e}$ consistent with the minor merger scenario.

We present the results of a pilot near-infrared (NIR) spectroscopic campaign of five very massive galaxies ($\log(\text{M}_\star/\text{M}_\odot)>11.45$) in the range of $1.7<z<2.7$. We measure an absorption feature redshift for one galaxy at $z_\text{spec}=2.000\pm0.006$. For the remaining galaxies, we combine the photometry with the continuum from the spectra to estimate continuum redshifts and stellar population properties. We define a continuum redshift ($z_{\rm cont}$ ) as one in which the redshift is estimated probabilistically using EAZY from the combination of catalog photometry and the observed spectrum. We derive the uncertainties on the stellar population synthesis properties using a Monte Carlo simulation and examine the correlations between the parameters with and without the use of the spectrum in the modeling of the spectral energy distributions (SEDs). The spectroscopic constraints confirm the extreme stellar masses of the galaxies in our sample. We find that three out of five galaxies are quiescent (star formation rate of $\lesssim 1 M_\odot~yr^{-1}$) with low levels of dust obscuration ($A_{\rm V} < 1$) , that one galaxy displays both high levels of star formation and dust obscuration (${\rm SFR} \approx 300 M_\odot~{\rm yr}^{-1}$, $A_{\rm V} \approx 1.7$~mag), and that the remaining galaxy has properties that are intermediate between the quiescent and star-forming populations.

We make publicly available a catalog of calibrated environmental measures for galaxies in the five 3D-HST/CANDELS deep fields. Leveraging the spectroscopic and grism redshifts from the 3D-HST survey, multi wavelength photometry from CANDELS, and wider field public data for edge corrections, we derive densities in fixed apertures to characterize the environment of galaxies brighter than $JH_{140} < 24$ mag in the redshift range $0.5<z<3.0$. By linking observed galaxies to a mock sample, selected to reproduce the 3D-HST sample selection and redshift accuracy, each 3D-HST galaxy is assigned a probability density function of the host halo mass, and a probability that is a central or a satellite galaxy. The same procedure is applied to a $z=0$ sample selected from SDSS. We compute the fraction of passive central and satellite galaxies as a function of stellar and halo mass, and redshift, and then derive the fraction of galaxies that were quenched by environment specific processes. Using the mock sample, we estimate that the timescale for satellite quenching is $t_{\rm quench} \sim 2-5$ Gyr; longer at lower stellar mass or lower redshift, but remarkably independent of halo mass. This indicates that, in the range of environments commonly found within the 3D-HST sample, satellites are quenched by exhaustion of their gas reservoir in absence of cosmological accretion. We find that the quenching times can be separated into a delay phase during which satellite galaxies behave similarly to centrals at fixed stellar mass, and a phase where the star formation rate drops rapidly ($\sim 0.4-0.6$ Gyr), as shown previously at $z=0$. We conclude that this scenario requires satellite galaxies to retain a large reservoir of multi-phase gas upon accretion, even at high redshift, and that this gas sustains star formation for the long quenching times observed.

SDSS J2222+2745 is a galaxy cluster at z=0.49, strongly lensing a quasar at z=2.805 into six widely separated images. In recent HST imaging of the field, we identify additional multiply lensed galaxies, and confirm the sixth quasar image that was identified by Dahle et al. (2013). We used the Gemini North telescope to measure a spectroscopic redshift of z=4.56 of one of the secondary lensed galaxies. These data are used to refine the lens model of SDSS J2222+2745, compute the time delay and magnifications of the lensed quasar images, and reconstruct the source image of the quasar host and a second lensed galaxy at z=2.3. This second galaxy also appears in absorption in our Gemini spectra of the lensed quasar, at a projected distance of 34 kpc. Our model is in agreement with the recent time delay measurements of Dahle et al. (2015), who found tAB=47.7+/-6.0 days and tAC=-722+/-24 days. We use the observed time delays to further constrain the model, and find that the model-predicted time delays of the three faint images of the quasar are tAD=502+/-68 days, tAE=611+/-75 days, and tAF=415+/-72 days. We have initiated a follow-up campaign to measure these time delays with Gemini North. Finally, we present initial results from an X-ray monitoring program with Swift, indicating the presence of hard X-ray emission from the lensed quasar, as well as extended X-ray emission from the cluster itself, which is consistent with the lensing mass measurement and the cluster velocity dispersion.

Early quiescent galaxies at z~2 are known to be remarkably compact compared to their nearby counterparts. Possible progenitors of these systems include galaxies that are structurally similar, but are still rapidly forming stars. Here, we present Karl G. Jansky Very Large Array (VLA) observations of the CO(1-0) line towards three such compact, star-forming galaxies at z~2.3, significantly detecting one. The VLA observations indicate baryonic gas fractions >~5 times lower and gas depletion times >~10 times shorter than normal, extended massive star-forming galaxies at these redshifts. At their current star formation rates, all three objects will deplete their gas reservoirs within 100Myr. These objects are among the most gas-poor objects observed at z>2, and are outliers from standard gas scaling relations, a result which remains true regardless of assumptions about the CO-H2 conversion factor. Our observations are consistent with the idea that compact, star-forming galaxies are in a rapid state of transition to quiescence in tandem with the build-up of the z~2 quenched population. In the detected compact galaxy, we see no evidence of rotation or that the CO-emitting gas is spatially extended relative to the stellar light. This casts doubt on recent suggestions that the gas in these compact galaxies is rotating and significantly extended compared to the stars. Instead, we suggest that, at least for this object, the gas is centrally concentrated, and only traces a small fraction of the total galaxy dynamical mass.

The FourStar galaxy evolution survey (ZFOURGE) is a 45 night legacy program with the FourStar near-infrared camera on Magellan and one of the most sensitive surveys to date. ZFOURGE covers a total of $400\ \mathrm{arcmin}^2$ in cosmic fields CDFS, COSMOS and UDS, overlapping CANDELS. We present photometric catalogs comprising $>70,000$ galaxies, selected from ultradeep $K_s$-band detection images ($25.5-26.5$ AB mag, $5\sigma$, total), and $>80\%$ complete to $K_s<25.3-25.9$ AB. We use 5 near-IR medium-bandwidth filters ($J_1,J_2,J_3,H_s,H_l$) as well as broad-band $K_s$ at $1.05\ - 2.16\ \mu m$ to $25-26$ AB at a seeing of $\sim0.5$". Each field has ancillary imaging in $26-40$ filters at $0.3-8\ \mu m$. We derive photometric redshifts and stellar population properties. Comparing with spectroscopic redshifts indicates a photometric redshift uncertainty $\sigma_z={0.010,0.009}$, and 0.011 in CDFS, COSMOS, and UDS. As spectroscopic samples are often biased towards bright and blue sources, we also inspect the photometric redshift differences between close pairs of galaxies, finding $\sigma_{z,pairs}= 0.01-0.02$ at $1<z<2.5$. We quantify how $\sigma_{z,pairs}$ depends on redshift, magnitude, SED type, and the inclusion of FourStar medium bands. $\sigma_{z,pairs}$ is smallest for bright, blue star-forming samples, while red star-forming galaxies have the worst $\sigma_{z,pairs}$. Including FourStar medium bands reduces $\sigma_{z,pairs}$ by 50\% at $1.5<z<2.5$. We calculate SFRs based on ultraviolet and ultradeep far-IR $Spitzer$/MIPS and Herschel/PACS data. We derive rest-frame $U-V$ and $V-J$ colors, and illustrate how these correlate with specific SFR and dust emission to $z=3.5$. We confirm the existence of quiescent galaxies at $z\sim3$, demonstrating their SFRs are suppressed by $>\times15$.

Using the UltraVISTA DR1 and 3D-HST catalogs, we construct a stellar-mass-complete sample, unique for its combination of surveyed volume and depth, to study the evolution of the fractions of quiescent galaxies, moderately unobscured star-forming galaxies, and dusty star-forming galaxies as a function of stellar mass over the redshift interval $0.2 \le z \le 3.0$. We show that the role of dusty star-forming galaxies within the overall galaxy population becomes more important with increasing stellar mass, and grows rapidly with increasing redshift. Specifically, dusty star-forming galaxies dominate the galaxy population with $\log{(M_{\rm star}/M_{\odot})} \gtrsim 10.3$ at $z\gtrsim2$. The ratio of dusty and non-dusty star-forming galaxies as a function of stellar mass changes little with redshift. Dusty star-forming galaxies dominate the star-forming population at $\log{(M_{\rm star}/M_{\odot})} \gtrsim 10.0-10.5$, being a factor of $\sim$3-5 more common, while unobscured star-forming galaxies dominate at $\log{(M_{\rm star}/M_{\odot})} \lesssim 10$. At $\log{(M_{\rm star}/M_{\odot})} > 10.5$, red galaxies dominate the galaxy population at all redshift $z<3$, either because they are quiescent (at late times) or dusty star-forming (in the early universe).

We present an overview of the "KIFF" project, which provides ultra-deep Ks-band imaging of all six of the Hubble Frontier Fields clusters Abell 2744, MACS-0416, Abell S1063, Abell 370, MACS-0717 and MACS-1149. All of these fields have recently been observed with large allocations of Directors' Discretionary Time with the HST and Spitzer telescopes covering 0.4 < lambda < 1.6 microns and 3.6--4.5 microns, respectively. VLT/HAWK-I integrations of the first four fields reach 5-sigma limiting depths of Ks~26.0 (AB, point sources) and have excellent image quality (FWHM ~ 0."4). Shorter Keck/MOSFIRE integrations of the MACS-0717 (MACS-1149) field better observable in the north reach limiting depths Ks=25.5 (25.1) with seeing FWHM ~0."4 (0."5). In all cases the Ks-band mosaics cover the primary cluster and parallel HST/ACS+WFC3 fields. The total area of the Ks-band coverage is 490 arcmin^2. The Ks-band at 2.2 microns crucially fills the gap between the reddest HST filter (1.6 micron ~ H-band) and the IRAC 3.6 micron passband. While reaching the full depths of the space-based imaging is not currently feasible from the ground, the deep Ks-band images provide important constraints on both the redshifts and the stellar population properties of galaxies extending well below the characteristic stellar mass across most of the age of the universe, down to, and including, the redshifts of the targeted galaxy clusters (z < 0.5).

Recent surveys have identified a seemingly ubiquitous population of galaxies with elevated [OIII]/H$\beta$ emission line ratios at $z > 1$, though the nature of this phenomenon continues to be debated. The [OIII]/H$\beta$ line ratio is of interest because it is a main component of the standard diagnostic tools used to differentiate between active galactic nuclei (AGN) and star-forming galaxies, as well as the gas-phase metallicity indicators $O_{23}$ and $R_{23}$. Here, we investigate the primary driver of increased [OIII]/H$\beta$ ratios by median-stacking rest-frame optical spectra for a sample of star-forming galaxies in the 3D-HST survey in the redshift range $z\sim1.4-2.2$. Using $N = 4220$ star-forming galaxies, we stack the data in bins of mass and specific star formation rates (sSFR) respectively. After accounting for stellar Balmer absorption, we measure [OIII]$\lambda5007$\AA/H$\beta$ down to $\mathrm{M} \sim 10^{9.2} \ \mathrm{M_\odot}$ and sSFR $\sim 10^{-9.6} \ \mathrm{yr}^{-1}$, more than an order of magnitude lower than previous work at similar redshifts. We find an offset of $0.59\pm0.05$ dex between the median ratios at $z\sim2$ and $z\sim0$ at fixed stellar mass, in agreement with existing studies. However, with respect to sSFR, the $z \sim 2$ stacks all lie within 1$\sigma$ of the median SDSS ratios, with an average offset of only $-0.06\pm 0.05$. We find that the excitation properties of galaxies are tightly correlated with their sSFR at both $z\sim2$ and $z\sim0$, with a relation that appears to be roughly constant over the last 10 Gyr of cosmic time.

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 explore the presence of non-stellar rest-frame near-IR ($2-5 \ \mu \mathrm{m}$) emission in galaxies at $z \sim 1$. Previous studies identified this excess in relatively small samples and suggested that such non-stellar emission, which could be linked to the $3.3 \ \mu \mathrm{m}$ polycyclic aromatic hydrocarbons feature or hot dust emission, is associated with an increased star formation rate (SFR). In this Letter, we confirm and quantify the presence of an IR excess in a significant fraction of galaxies in the 3D-HST GOODS catalogs. By constructing a matched sample of galaxies with and without strong non-stellar near-IR emission, we find that galaxies with such emission are predominantly star-forming galaxies. Moreover, star-forming galaxies with an excess show increased mid- and far-IR and H$\alpha$ emission compared to other star-forming galaxies without. While galaxies with a near-IR excess show a larger fraction of individually detected X-ray active galactic nuclei (AGNs), an X-ray stacking analysis, together with the IR-colors and H$\alpha$ profiles, shows that AGNs are unlikely to be the dominant source of the excess in the majority of galaxies. Our results suggest that non-stellar near-IR emission is linked to increased SFRs and is ubiquitous among star-forming galaxies. As such, the near-IR emission might be a powerful tool to measure SFRs in the era of the James Webb Space Telescope.

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 derive average radial gradients in the dust attenuation towards HII regions in 609 galaxies at z~1.4, using measurements of the Balmer decrement out to r~3kpc. The Balmer decrements are derived from spatially resolved maps of Halpha and Hbeta emission from the 3D-HST survey. We find that with increasing stellar mass (M) both the normalization and strength of the gradient in dust attenuation increases. Galaxies with a mean mass of <log(M)> = 9.2Msun have little dust attenuation at all radii, whereas galaxies with <log(M)>= 10.2Msun have dust attenuation toward Halpha A(Halpha)~2mag in their central regions. We parameterize this as A(Halpha) = b + c log(r), with b = 0.9 + 1.0 log(M10), c = -1.9 - 2.2 log(M10), r in kpc, and M10 the stellar mass in units of 10^10Msun. This expression can be used to correct spatially resolved measurements of Halpha to radial distributions of star formation. When applied to our data, we find that the star formation rates in the central r<1kpc of galaxies in the highest mass bin are ~ 6 Msun/yr, six times higher than before correction and approximately half of the total star formation rate of these galaxies. If this high central star formation rate is maintained for several Gyr, a large fraction of the stars in present-day bulges likely formed in-situ.

We explore star-formation histories (SFHs) of galaxies based on the evolution of the star-formation rate stellar mass relation (SFR-M*). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M* relation at 0.5 < z < 4. Similar to recent works we find that the average infrared SEDs of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that the SFR-M* relation is not consistent with a single power-law of the form SFR ~ M*^a at any redshift; it has a power-law slope of a~1 at low masses, and becomes shallower above a turnover mass (M_0) that ranges from 10^9.5 - 10^10.8 Msol, with evidence that M_0 increases with redshift. We compare our measurements to results from state-of-the-art cosmological simulations, and find general agreement in the slope of the SFR-M* relation albeit with systematic offsets. We use the evolving SFR-M* sequence to generate SFHs, finding that typical SFRs of individual galaxies rise at early times and decline after reaching a peak. This peak occurs earlier for more massive galaxies. We integrate these SFHs to generate mass-growth histories and compare to the implied mass-growth from the evolution of the stellar mass function. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass-growth rates that are as much as 10x higher than inferred from the stellar mass function. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers.

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.

We present a study of photometric redshift accuracy in the 3D-HST photometric catalogs, using 3D-HST grism redshifts to quantify and dissect trends in redshift accuracy for galaxies brighter than $H_{F140W}<24$ with an unprecedented and representative high-redshift galaxy sample. We find an average scatter of $0.0197\pm0.0003(1+z)$ in the Skelton et al. (2014) photometric redshifts. Photometric redshift accuracy decreases with magnitude and redshift, but does not vary monotonically with color or stellar mass. The 1-$\sigma$ scatter lies between $0.01-0.03$(1+z) for galaxies of all masses and colors below $z<2.5$ (for $H_{F140W}{<}24$), with the exception of a population of very red ($U-V > 2$), dusty star-forming galaxies for which the scatter increases to $\sim0.1(1+z)$. Although the overall photometric redshift accuracy for quiescent galaxies is better than for star-forming galaxies, scatter depends more strongly on magnitude and redshift than on galaxy type. We verify these trends using the redshift distributions of close pairs and extend the analysis to fainter objects, where photometric redshift errors further increase to $\sim0.046(1+z)$ at $H_{F160W}=26$. We demonstrate that photometric redshift accuracy is strongly filter-dependent and quantify the contribution of multiple filter combinations. We evaluate the widths of redshift probability distribution functions and find that error estimates are underestimated by a factor of $\sim1.1-1.6$, but that uniformly broadening the distribution does not adequately account for fitting outliers. Finally, we suggest possible applications of these data in planning for current and future surveys and simulate photometric redshift performance in the LSST, DES, and combined DES and VHS surveys.

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.

We study the rest-frame ultra-violet sizes of massive (~0.8 x 10^11 M_Sun) galaxies at 3.4<z<4.2, selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles to HST/WFC3/F160W images from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). Massive quiescent galaxies are very compact, with a median circularized half-light radius r_e = 0.63 +/- 0.18 kpc. Removing 5/16 (31%) sources with signs of AGN activity does not change the result. Star-forming galaxies have r_e = 2.0 +/- 0.60 kpc, 3.2 +/- 1.3 x larger than quiescent galaxies. Quiescent galaxies at z~4 are on average 6.0 +\- 0.17 x smaller than at z~0 and 1.9 +/- 0.7 x smaller than at z~2. Star-forming galaxies of the same stellar mass are 2.4 +/- 0.7 x smaller than at z~0. Overall, the size evolution at 0<z<4 is well described by a powerlaw, with r_e = 5.08 +/- 0.28 (1+z)^(-1.44+/-0.08) kpc for quiescent and r_e = 6.02 +/- 0.28 (1+z)^(-0.72+/-0.05) kpc for star-forming galaxies. Compact star-forming galaxies are rare in our sample: we find only 1/14 (7%) with r_e / (M / 10^11 M_Sun)^0.75 < 1.5, whereas 13/16 (81%) of the quiescent galaxies is compact. The number density of compact quiescent galaxies at z~4 is 1.8 +/- 0.8 x 10^-5 Mpc^-3 and increases rapidly, by >5 x, between 2<z<4. The paucity of compact star-forming galaxies at z~4 and their large rest-frame ultra-violet median sizes suggest that the formation phase of compact cores is very short and/or highly dust obscured.

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 investigate the stellar populations of 25 massive, galaxies ($\log[M_\ast/M_\odot] \geq 10.9$) at $1.5 < z < 2$ using data obtained with the K-band Multi-Object Spectrograph (KMOS) on the ESO VLT. Targets were selected to be quiescent based on their broadband colors and redshifts using data from the 3D-HST grism survey. The mean redshift of our sample is $\bar{z} = 1.75$, where KMOS YJ-band data probe age- and metallicity-sensitive absorption features in the rest-frame optical, including the $G$ band, Fe I, and high-order Balmer lines. Fitting simple stellar population models to a stack of our KMOS spectra, we derive a mean age of $1.03^{+0.13}_{-0.08}$ Gyr. We confirm previous results suggesting a correlation between color and age for quiescent galaxies, finding mean ages of $1.22^{+0.56}_{-0.19}$ Gyr and $0.85^{+0.08}_{-0.05}$ Gyr for the reddest and bluest galaxies in our sample. Combining our KMOS measurements with those obtained from previous studies at $0.2 < z < 2$ we find evidence for a $2-3$ Gyr spread in the formation epoch of massive galaxies. At $z < 1$ the measured stellar ages are consistent with passive evolution, while at $1 < z \lesssim2$ they appear to saturate at $\sim$1 Gyr, which likely reflects changing demographics of the (mean) progenitor population. By comparing to star-formation histories inferred for "normal" star-forming galaxies, we show that the timescales required to form massive galaxies at $z \gtrsim 1.5$ are consistent with the enhanced $\alpha$-element abundances found in massive local early-type galaxies.

Interpreting observations of distant galaxies in terms of constraints on physical parameters - such as stellar mass, star-formation rate (SFR) and dust optical depth - requires spectral synthesis modelling. We analyse the reliability of these physical parameters as determined under commonly adopted `classical' assumptions: star-formation histories assumed to be exponentially declining functions of time, a simple dust law and no emission-line contribution. Improved modelling techniques and data quality now allow us to use a more sophisticated approach, including realistic star-formation histories, combined with modern prescriptions for dust attenuation and nebular emission (Pacifici et al. 2012). We present a Bayesian analysis of the spectra and multi-wavelength photometry of 1048 galaxies from the 3D-HST survey in the redshift range 0.7<z<2.8 and in the stellar mass range 9<log(M/Mo)<12. We find that, using the classical spectral library, stellar masses are systematically overestimated (~0.1 dex) and SFRs are systematically underestimated (~0.6 dex) relative to our more sophisticated approach. We also find that the simultaneous fit of photometric fluxes and emission-line equivalent widths helps break a degeneracy between SFR and optical depth of the dust, reducing the uncertainties on these parameters. Finally, we show how the biases of classical approaches can affect the correlation between stellar mass and SFR for star-forming galaxies (the `Star-Formation Main Sequence'). We conclude that the normalization, slope and scatter of this relation strongly depend on the adopted approach and demonstrate that the classical, oversimplified approach cannot recover the true distribution of stellar mass and SFR.

We examine the connection between the observed star-forming sequence (SFR $\propto$ $M^{\alpha}$) and the observed evolution of the stellar mass function between $0.2 < z < 2.5$. We find the star-forming sequence cannot have a slope $\alpha$ $\lesssim$ 0.9 at all masses and redshifts, as this would result in a much higher number density at $10 < \log(\mathrm{M/M_{\odot}}) < 11$ by $z=1$ than is observed. We show that a transition in the slope of the star-forming sequence, such that $\alpha=1$ at $\log(\mathrm{M/M_{\odot}})<10.5$ and $\alpha=0.7-0.13z$ ({Whitaker} {et~al.} 2012) at $\log(\mathrm{M/M_{\odot}})>10.5$, greatly improves agreement with the evolution of the stellar mass function. We then derive a star-forming sequence which reproduces the evolution of the mass function by design. This star-forming sequence is also well-described by a broken-power law, with a shallow slope at high masses and a steep slope at low masses. At $z=2$, it is offset by $\sim$0.3 dex from the observed star-forming sequence, consistent with the mild disagreement between the cosmic SFR and recent observations of the growth of the stellar mass density. It is unclear whether this problem stems from errors in stellar mass estimates, errors in SFRs, or other effects. We show that a mass-dependent slope is also seen in other self-consistent models of galaxy evolution, including semi-analytical, hydrodynamical, and abundance-matching models. As part of the analysis, we demonstrate that neither mergers nor hidden low-mass quiescent galaxies are likely to reconcile the evolution of the mass function and the star-forming sequence. These results are supported by observations from {Whitaker} {et~al.} (2014).

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.

We constrain the slope of the star formation rate ($\log\Psi$) to stellar mass ($\log\mathrm{M_{\star}}$) relation down to $\log(\mathrm{M_{\star}/M_{\odot}})=8.4$ ($\log(\mathrm{M_{\star}/M_{\odot}})=9.2$) at $z=0.5$ ($z=2.5$) with a mass-complete sample of 39,106 star-forming galaxies selected from the 3D-HST photometric catalogs, using deep photometry in the CANDELS fields. For the first time, we find that the slope is dependent on stellar mass, such that it is steeper at low masses ($\log\mathrm{\Psi}\propto\log\mathrm{M_{\star}}$) than at high masses ($\log\mathrm{\Psi}\propto(0.3-0.6)\log\mathrm{M_{\star}}$). These steeper low mass slopes are found for three different star formation indicators: the combination of the ultraviolet (UV) and infrared (IR), calibrated from a stacking analysis of Spitzer/MIPS 24$\mu$m imaging; $\beta$-corrected UV SFRs; and H$\alpha$ SFRs. The normalization of the sequence evolves differently in distinct mass regimes as well: for galaxies less massive than $\log(\mathrm{M_{\star}/M_{\odot}})<10$ the specific SFR ($\Psi/\mathrm{M_{\star}}$) is observed to be roughly self-similar with $\Psi/\mathrm{M_{\star}}\propto(1+z)^{1.9}$, whereas more massive galaxies show a stronger evolution with $\Psi/\mathrm{M_{\star}}\propto(1+z)^{2.2-3.5}$ for $\log(\mathrm{M_{\star}/M_{\odot}})=10.2-11.2$. The fact that we find a steep slope of the star formation sequence for the lower mass galaxies will help reconcile theoretical galaxy formation models with the observations. The results of this study support the analytical conclusions of Leja et al. (2014).

We determine the intrinsic, 3-dimensional shape distribution of star-forming galaxies at 0<z<2.5, as inferred from their observed projected axis ratios. In the present-day universe star-forming galaxies of all masses 1e9 - 1e11 Msol are predominantly thin, nearly oblate disks, in line with previous studies. We now extend this to higher redshifts, and find that among massive galaxies (M* > 1e10 Msol) disks are the most common geometric shape at all z < 2. Lower-mass galaxies at z>1 possess a broad range of geometric shapes: the fraction of elongated (prolate) galaxies increases toward higher redshifts and lower masses. Galaxies with stellar mass 1e9 Msol (1e10 Msol) are a mix of roughly equal numbers of elongated and disk galaxies at z~1 (z~2). This suggests that galaxies in this mass range do not yet have disks that are sustained over many orbital periods, implying that galaxies with present-day stellar mass comparable to that of the Milky Way typically first formed such sustained stellar disks at redshift z~1.5-2. Combined with constraints on the evolution of the star formation rate density and the distribution of star formation over galaxies with different masses, our findings imply that, averaged over cosmic time, the majority of stars formed in disks.

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.

We present a detailed analysis of Hubble Space Telescope (HST), Wide Field Camera 3 (WFC3) G141 grism spectroscopy for seven star-forming regions of the highly magnified lensed starburst galaxy RCSGA 032727-132609 at z=1.704. We measure the spatial variations of the extinction in RCS0327 through the observed H$\gamma$/H$\beta$ emission line ratios, finding a constant average extinction of $\mathrm{E(B-V)_{gas}}=0.40\pm0.07$. We infer that the star formation is enhanced as a result of an ongoing interaction, with measured star formation rates derived from demagnified, extinction-corrected H$\beta$ line fluxes for the individual star-forming clumps falling >1-2 dex above the star formation sequence. When combining the HST/WFC3 [OIII]$\lambda$5007/H$\beta$ emission line ratio measurements with [NII]/H$\alpha$ line ratios from Wuyts et a. (2014), we find that the majority of the individual star-forming regions fall along the local "normal" abundance sequence. With the first detections of the He I $\lambda$5876 and He II $\lambda$4686 recombination lines in a distant galaxy, we probe the massive-star content of the star-forming regions in RCS0327. The majority of the star-forming regions have a He I $\lambda$5876 to H$\beta$ ratio consistent with the saturated maximum value, which is only possible if they still contain hot O-stars. Two regions have lower ratios, implying that their last burst of new star formation ended ~5 Myr ago. Together, the He I $\lambda$5876 and He II $\lambda$4686 to H$\beta$ line ratios provide indirect evidence for the order in which star formation is stopping in individual star-forming knots of this high redshift merger. We place the spatial variations of the extinction, star formation rate and ionization conditions in the context of the star formation history of RCS0327.

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 0<z<3. Separating early- and late-type galaxies on the basis of star-formation activity, we confirm that early-type galaxies are on average smaller than late-type galaxies at all redshifts, and find a significantly different rate of average size evolution at fixed galaxy mass, with fast evolution for the early-type population, R_eff ~ (1+z)^-1.48, and moderate evolution for the late-type population, R_eff ~ (1+z)^-0.75. The large sample size and dynamic range in both galaxy mass and redshift, in combination with the high fidelity of our measurements due to the extensive use of spectroscopic data, not only fortify previous results, but also enable us to probe beyond simple average galaxy size measurements. At all redshifts the slope of the size-mass relation is shallow, R_eff ~ M_star^0.22, for late-type galaxies with stellar mass >3x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times.

We present near-infrared spectroscopy of a sample of 22 Extreme Emission Line Galaxies at redshifts 1.3 < z < 2.3, confirming that these are low-mass (M* = 10^8 - 10^9 M_sun) galaxies undergoing intense starburst episodes (M*/SFR ~ 10-100 Myr). The sample is selected by [O III] or H{\alpha} emission line flux and equivalent width using near-infrared grism spectroscopy from the 3D-HST survey. High-resolution NIR spectroscopy is obtained with LBT/LUCI and VLT/X-SHOOTER. The [O III]/H{\beta} line ratio is high (> 5) and [N II]/H{\alpha} is always significantly below unity, which suggests a low gas-phase metallicity. We are able to determine gas-phase metallicities for 7 of our objects using various strong-line methods, with values in the range 0.05-0.30 Z_sun and with a median of 0.15 Z_sun; for 3 of these objects we detect [O III]{\lambda}4363 which allows for a direct constraint on the metallicity. The velocity dispersion, as measured from the nebular emission lines, is typically ~50 km/s. Combined with the observed star-forming activity, the Jeans and Toomre stability criteria imply that the gas fraction must be large (> 2/3), consistent with the difference between our dynamical and stellar mass estimates. The implied gas depletion time scale (several hundred Myr) is substantially longer than the inferred mass-weighted ages (~50 Myr), which further supports the emerging picture that most stars in low-mass galaxies form in short, intense bursts of star formation.

Most massive galaxies are thought to have formed their dense stellar cores at early cosmic epochs. However, cores in their formation phase have not yet been observed. Previous studies have found galaxies with high gas velocity dispersions or small apparent sizes but so far no objects have been identified with both the stellar structure and the gas dynamics of a forming core. Here we present a candidate core in formation 11 billion years ago, at z=2.3. GOODS-N-774 has a stellar mass of 1.0x10^11 Msun, a half-light radius of 1.0 kpc, and a star formation rate of 90[+45-20]Msun/yr. The star forming gas has a velocity dispersion 317+-30 km/s, amongst the highest ever measured. It is similar to the stellar velocity dispersions of the putative descendants of GOODS-N-774, compact quiescent galaxies at z~2 and giant elliptical galaxies in the nearby Universe. Galaxies such as GOODS-N-774 appear to be rare; however, from the star formation rate and size of the galaxy we infer that many star forming cores may be heavily obscured, and could be missed in optical and near-infrared surveys.

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.

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.

We present direct observational evidence for star formation quenching in galaxy groups in the redshift range 0<z<2.5. We utilize a large sample of nearly 6000 groups, selected by fixed cumulative number density from three photometric catalogs, to follow the evolving quiescent fractions of central and satellite galaxies over roughly 11 Gyr. At z~0, central galaxies in our sample range in stellar mass from Milky Way/M31 analogs (M=6.5x10^10 M\solar) to nearby massive ellipticals (M=1.5x10^11 M\solar). Satellite galaxies in the same groups reach masses as low as twice that of the Large Magellanic Cloud (M=6.5x10^9 M\solar). Using statistical background subtraction, we measure the average rest-frame colors of galaxies in our groups and calculate the evolving quiescent fractions of centrals and satellites over seven redshift bins. Our analysis shows clear evidence for star formation quenching in group halos, with a different quenching onset for centrals and their satellite galaxies. Using halo mass estimates for our central galaxies, we find that star formation shuts off in centrals when typical halo masses reach between 10^12 and 10^13 M\solar, consistent with predictions from the halo quenching model. In contrast, satellite galaxies in the same groups most likely undergo quenching by environmental processes, whose onset is delayed with respect to their central galaxy. Although star formation is suppressed in all galaxies over time, the processes that govern quenching are different for centrals and satellites. While mass plays an important role in determining the star formation activity of central galaxies, quenching in satellite galaxies is dominated by the environment in which they reside.

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.

Spectroscopic observations from the Large Binocular Telescope and the Very Large Telescope reveal kinematically narrow lines (~50 km/s) for a sample of 14 Extreme Emission Line Galaxies (EELGs) at redshifts 1.4 < z < 2.3. These measurements imply that the total dynamical masses of these systems are low (< 3x10^9 M_sun). Their large [O III] 5007 equivalent widths (500-1100 Angstroms) and faint blue continuum emission imply young ages of 10-100 Myr and stellar masses of 10^8-10^9 M_sun, confirming the presence of a violent starburst. The dynamical masses represent the first such determinations for low-mass galaxies at z > 1. The stellar mass formed in this vigorous starburst phase represents a large fraction of the total (dynamical) mass, without a significantly massive underlying population of older stars. The occurrence of such intense events in shallow potentials strongly suggests that supernova-driven winds must be of critical importance in the subsequent evolution of these systems.

We follow the structural evolution of star forming galaxies (SFGs) like the Milky Way by selecting progenitors to z~1.3 based on the stellar mass growth inferred from the evolution of the star forming sequence. We select our sample from the 3D-HST survey, which utilizes spectroscopy from the HST WFC3 G141 near-IR grism and enables precise redshift measurements for our sample of SFGs. Structural properties are obtained from Sersic profile fits to CANDELS WFC3 imaging. The progenitors of z=0 SFGs with stellar mass M=10^{10.5} Msun are typically half as massive at z~1. This late-time stellar mass assembly is consistent with recent studies that employ abundance matching techniques. The descendant SFGs at z~0 have grown in half-light radius by a factor of ~1.4 since z~1. The half-light radius grows with stellar mass as r_e M^{0.29}. While most of the stellar mass is clearly assembling at large radii, the mass surface density profiles reveal ongoing mass growth also in the central regions where bulges and pseudobulges are common features in present day late-type galaxies. Some portion of this growth in the central regions is due to star formation as recent observations of H-alpha maps for SFGs at z~1 are found to be extended but centrally peaked. Connecting our lookback study with galactic archeology, we find the stellar mass surface density at R=8 kpc to have increased by a factor of ~2 since z~1, in good agreement with measurements derived for the solar neighborhood of the Milky Way.

We analyze the resolved stellar populations of 473 massive star-forming galaxies at 0.7 < z < 1.5, with multi-wavelength broad-band imaging from CANDELS and Halpha surface brightness profiles at the same kiloparsec resolution from 3D-HST. Together, this unique data set sheds light on how the assembled stellar mass is distributed within galaxies, and where new stars are being formed. We find the Halpha morphologies to resemble more closely those observed in the ACS I band than in the WFC3 H band, especially for the larger systems. We next derive a novel prescription for Halpha dust corrections, which accounts for extra extinction towards HII regions. The prescription leads to consistent SFR estimates and reproduces the observed relation between the Halpha/UV luminosity ratio and visual extinction, both on a pixel-by-pixel and on a galaxy-integrated level. We find the surface density of star formation to correlate with the surface density of assembled stellar mass for spatially resolved regions within galaxies, akin to the so-called 'main sequence of star formation' established on a galaxy-integrated level. Deviations from this relation towards lower equivalent widths are found in the inner regions of galaxies. Clumps and spiral features, on the other hand, are associated with enhanced Halpha equivalent widths, bluer colors, and higher specific star formation rates compared to the underlying disk. Their Halpha/UV luminosity ratio is lower than that of the underlying disk, suggesting the ACS clump selection preferentially picks up those regions of elevated star formation activity that are the least obscured by dust. Our analysis emphasizes that monochromatic studies of galaxy structure can be severely limited by mass-to-light ratio variations due to dust and spatially inhomogeneous star formation histories.

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.

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.

3D-HST is a 248-orbit Treasury program to provide WFC3 and ACS grism spectroscopy over four extra-galactic fields (AEGIS, COSMOS, GOODS-South, and UDS), augmented with previously obtained data in GOODS-North. We present a new data release of the 3D-HST survey, version v3.0. This release follows the initial v0.5 release that accompanied the survey description paper (Brammer et al. 2012). The new v3.0 release includes the deepest near-IR HST grism spectra currently in existence, extracted from the 8-17 orbit depth observations in the Hubble Ultra Deep Field. Contamination-corrected 2D and 1D spectra, as well as derived redshifts, are made available for >250 objects in this 2'x2'field. The spectra are of extraordinary quality, and show emission features in many galaxies as faint as F140W=26-27, absorption features in quiescent galaxies at z~2, and several active galactic nuclei. In addition to these extremely deep grism data we provide reduced WFC3 F125W, F140W, and F160W image mosaics of all five 3D-HST/CANDELS fields.

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.

We present a statistical study of the environments of massive galaxies in four redshift bins between z=0.04 and z=1.6, using data from the Sloan Digital Sky Survey (SDSS) and the NEWFIRM Medium Band Survey (NMBS). We measure the projected radial distribution of galaxies in cylinders around a constant number density selected sample of massive galaxies and utilize a statistical subtraction of contaminating sources. Our analysis shows that massive primary galaxies typically live in group halos and are surrounded by 2 to 3 satellites with masses more than one-tenth of the primary galaxy mass. The cumulative stellar mass in these satellites roughly equals the mass of the primary galaxy itself. We further find that the radial number density profile of galaxies around massive primaries has not evolved significantly in either slope or overall normalization in the past 9.5 Gyr. A simplistic interpretation of this result can be taken as evidence for a lack of mergers in the studied groups and as support for a static evolution model of halos containing massive primaries. Alternatively, there exists a tight balance between mergers and accretion of new satellites such that the overall distribution of galaxies in and around the halo is preserved. The latter interpretation is supported by a comparison to a semi-analytic model, which shows a similar constant average satellite distribution over the same redshift range.

[Abridged] Using public data from the NMBS and CANDELS surveys, we study the population of massive galaxies at z>3 to identify the potential progenitors of z~2 compact, massive, quiescent (CMQ) galaxies, furthering our understanding of the evolution of massive galaxies. Our work is enabled by high-resolution CANDELS images and accurate photometric redshifts, stellar masses and star formation rates (SFRs) from 37-band NMBS photometry. The total number of z>3 massive galaxies is consistent with the number of massive quiescent (MQ) galaxies at z~2, implying that the SFRs for all of these galaxies must be much lower by z~2. We discover 4 CMQ galaxies at z>3, pushing back the time for which such galaxies have been observed. However, the volume density for these galaxies is significantly less than that of galaxies at z<2 with similar masses, SFRs, and sizes, implying that additional CMQ galaxies must be created in the ~1 Gyr between z=3 and z=2. We find 5 star-forming galaxies at z~3 that are compact (Re<1.4 kpc) and have stellar mass M*>10^(10.6)Msun, likely to become members of the CMQ galaxy population at z~2. We evolve the stellar masses and SFRs of each individual z>3 galaxy adopting 5 different star formation histories (SFHs) and studying the resulting population of massive galaxies at z=2.3. We find that declining or truncated SFHs are necessary to match the observed number density of MQ galaxies at z~2, whereas a constant SFH results in a number density significantly smaller than observed. All of our assumed SFHs imply number densities of CMQ galaxies at z~2 that are consistent with the observed number density. Better agreement with the observed number density of CMQ galaxies at z~2 is obtained if merging is included in the analysis and better still if star formation quenching is assumed to shortly follow the merging event, as implied by recent models of formation of MQ galaxies.

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 Hubble Space Telescope (HST) imaging and spectroscopy of the gravitational lens SL2SJ02176-0513, a cusp arc at z=1.847. The UV continuum of the lensed galaxy is very blue, which is seemingly at odds with its redder optical colors. The 3D-HST WFC3/G141 near-infrared spectrum of the lens reveals the source of this discrepancy to be extremely strong [OIII]5007 and H-beta emission lines with rest-frame equivalent widths of 2000 +/- 100 and 520 +/- 40 Angstroms, respectively. The source has a stellar mass ~10^8 Msun, sSFR \sim 100/Gyr, and detection of [OIII]4363 yields a metallicity of 12 + log(O/H) = 7.5 +/- 0.2. We identify local blue compact dwarf analogs to SL2SJ02176-0513, which are among the most metal-poor galaxies in the SDSS. The local analogs resemble the lensed galaxy in many ways, including UV/optical SED, spatial morphology and emission line equivalent widths and ratios. Common to SL2SJ02176-0513 and its local counterparts is an upturn at mid-IR wavelengths likely arising from hot dust heated by starbursts. The emission lines of SL2SJ02176-0513 are spatially resolved owing to the combination of the lens and the high spatial resolution of HST. The lensed galaxy is composed of two clumps with combined size r_e \sim 300 pc, and we resolve significant differences in UV color and emission line equivalent width between them. Though it has characteristics occasionally attributed to active galactic nuclei, we conclude that SL2SJ02176-0513 is a low-metallicity star-bursting dwarf galaxy. Such galaxies will be found in significant numbers in the full 3D-HST grism survey.

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 study the star formation rate (SFR) - stellar mass (M*) relation in a self-consistent manner from 0 < z < 2.5 with a sample of galaxies selected from the NEWFIRM Medium-Band Survey. We find a significant non-linear slope of the relation, SFR \propto M*^0.6, and a constant observed scatter of 0.34 dex, independent of redshift and M*. However, if we select only blue galaxies we find a linear relation SFR \propto M*, similar to previous results at z = 0 by Peng et al. (2010). This selection excludes red, dusty, star-forming galaxies with higher masses, which brings down the slope. By selecting on L_IR/L_UV (a proxy for dust obscuration) and the rest-frame U-V colors, we show that star-forming galaxies fall in three distinct regions of the log(SFR)-log(M*) plane: 1) actively star-forming galaxies with "normal" dust obscuration and associated colors (54% for log(M*) > 10 at 1 < z < 1.5), 2) red star-forming galaxies with low levels of dust obscuration and low specific SFRs (11%), and 3) dusty, blue star-forming galaxies with high specific SFRs (7%). The remaining 28% comprises quiescent galaxies. Galaxies on the "normal" star formation sequence show strong trends of increasing dust attenuation with stellar mass and a decreasing specific SFR, with an observed scatter of 0.25 dex (0.17 dex intrinsic scatter). The dusty, blue galaxies reside in the upper envelope of the star formation sequence with remarkably similar spectral shapes at all masses, suggesting that the same physical process is dominating the stellar light. The red, low-dust star-forming galaxies may be in the process of shutting off and migrating to the quiescent population.

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 present the rest-frame V-band luminosity function (LF) of galaxies at 0.4<z<4.0, measured from a near-infrared selected sample constructed from the NMBS, the FIRES, the FIREWORKS, and the ultra-deep NICMOS and WFC3 observations in the HDFN, HUDF, and GOODS-CDFS, all having high-quality optical to mid-infrared data. This unique sample combines data from surveys with a large range of depths and areas in a self-consistent way, allowing us to (1) minimize the uncertainties due to cosmic variance; and (2) simultaneously constrain the bright and faint ends with unprecedented accuracy over the targeted redshift range, probing the LF down to 0.1 L* at z~3.9. We find that (1) the faint end is fairly flat and with a constant slope from z=4, with alpha= -1.27 +/- 0.05; (2) the characteristic magnitude has dimmed by 1.3 mag from z~3.7 to z=0.1; (3) the characteristic density has increased by a factor of ~8 from z~3.7 to z=0.1, with 50% of this increase from z~4 to z~1.8; and (4) the luminosity density peaks at z~1-1.5, increasing by a factor of ~4 from z=4.0 to z~1-1.5, and subsequently decreasing by a factor of ~1.5 by z=0.1. We find no evidence for a steepening of the faint-end slope with redshift out to z=4, in contrast with previous observational claims and theoretical predictions. The constant faint-end slope suggests that the efficiency of stellar feedback may evolve with redshift. Alternative interpretations are discussed, such as different masses of the halos hosting faint galaxies at low and high redshifts and/or environmental effects.

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

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.

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

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.

Mergers between red galaxies are observed to be common in the nearby Universe, and are thought to be the dominant mechanism by which massive galaxies grow their mass at late times. These ``dry'' mergers can be readily identified in very deep ground based images, thanks to their extended low surface brightness tidal features. However, ground-based images lack the required resolution to determine the morphologies of the merging galaxies, and to measure the amount of dust and associated gas. We present HST/ACS and WFPC2 observations of a sample of 31 bulge-dominated red-sequence galaxies at z~0.1, comprised of ongoing mergers, merger remnants, and undisturbed galaxies. Nearly all galaxies have early-type morphologies and most are well-fit by r^1/4 law surface brightness profiles. We find that only 10% of the galaxies show evidence for the presence of dust. The amount of cold gas (or its upper limit) is calculated from the mean color-excess, assuming a simple relation between gas mass and dust mass. The gas mass is low for all galaxies, and we find that Mgas/Mstellar < 3x10^-4. We infer that red mergers in the nearby Universe mostly involve early-type galaxies containing little cold gas and dust. This may imply that the progenitors were mostly devoid of gas and/or that feedback mechanisms are very effective in preventing the gas to cool. The lack of gas in these objects may also imply a relatively large fraction of binary black holes in the centers of massive ellipticals.

We detect 353 X-ray point sources, mostly low-mass X-ray binaries (LMXBs), in four Chandra observations of Centaurus A (NGC 5128), the nearest giant early-type galaxy, and correlate this point source population with the largest available ensemble of confirmed and likely globular clusters associated with this galaxy. Of the X-ray sources, 31 are coincident with 30 globular clusters that are confirmed members of the galaxy by radial velocity measurement (2 X-ray sources match one globular cluster within our search radius), while 1 X-ray source coincides with a globular cluster resolved by HST images. Another 36 X-ray point sources match probable, but spectroscopically unconfirmed, globular cluster candidates. The color distribution of globular clusters and cluster candidates in Cen A is bimodal, and the probability that a red, metal rich GC candidate contains an LMXB is at least 1.7 times that of a blue, metal poor one. If we consider only spectroscopically confirmed GCs, this ratio increases to ~3. We find that LMXBs appear preferentially in more luminous (massive) GCs. These two effects are independent, and the latter is likely a consequence of enhanced dynamical encounter rates in more massive clusters which have on average denser cores. The X-ray luminosity functions of the LMXBs found in GCs and of those that are unmatched with GCs reveal similar underlying populations, though there is some indication that fewer X-ray faint LMXBs are found in globular clusters than X-ray bright ones. Our results agree with previous observations of the connection of GCs and LMXBs in early-type galaxies and extend previous work on Centaurus A.

We present an X-ray study of our nearest starburst galaxy IC 10, based on XMM-Newton and Chandra observations. A list of 73 XMM-Newton and 28 Chandra detections of point-like X-ray sources in the field is provided; a substantial fraction of them are likely stellar objects in the Milky Way due to the low Galactic latitude location of IC 10. The brightest source in the IC 10 field, X-1, shows a large variation by a factor of up to $\sim 6$ on time scales during the XMM-Newton observation. The X-ray spectra of the source indicate the presence of a multi-color blackbody accretion disk with an inner disk temperature T_{in} \approx 1.1 keV. These results are consistent with the interpretation of the source as a stellar mass black hole, probably accreting from a Wolf-Rayet star companion. We infer the mass of this black hole to be about 4 Msun if it is not spinning, or a factor of up to about 6 higher if there is significant spinning. We also detect an apparent diffuse soft X-ray emission component of IC 10. An effective method is devised to remove the X-ray CCD-readout streaks of X-1 that strongly affect the study of the diffuse component in the XMM-Newton and Chandra observations. We find that the diffuse X-ray morphology is oriented along the optical body of the galaxy and is chiefly associated with starburst regions. The diffuse component can be characterized by an optically thin thermal plasma with a mean temperature of $\sim 4 \times 10^6$ K and a 0.5-2 keV luminosity of $\sim 8 \times 10^{37} {\rm erg s^{-1}}$, representing only a small fraction of the expected mechanical energy inputs from massive stars in the galaxy. There is evidence that the hot gas is driving outflows from the starburst regions; therefore, the bulk of the energy inputs may be released in a galactic wind.