We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mass. Massive galaxies at z > 3 have either a flat or declining star formation history (SFH), whereas similar stellar mass galaxies at z~2.0 exhibit a slight peak in the past 500 Myrs. Comparing with the IllustrisTNG cosmological simulation, we find that (i) the dynamical mass of massive galaxies in simulations (log(M_star) > 10.0) increases by ~0.1 dex at a fixed stellar mass between z=2.0-3.0, and (ii) dynamical mass growth is coupled with a rapid rise in the ex-situ stellar mass fraction (stars accreted from other galaxies) for massive galaxies at z < 3.5. We speculate that the rising contribution of ex-situ stellar mass to the total stellar mass growth of massive galaxies is driving the higher integrated velocity dispersion and rising SFHs of massive galaxies at z~2.0 compared to galaxies of similar stellar masses at z > 3.

We present spectra of the most massive quiescent galaxy yet discovered at $z>3$, spectroscopically confirmed via the detection of Balmer absorption features in the $H-$ and $K-$bands of Keck/MOSFIRE. The spectra confirm a galaxy with no significant ongoing star formation, consistent with the lack of rest-frame UV flux and overall photometric spectral energy distribution. With a stellar mass of $3.1^{+0.1}_{-0.2} \times 10^{11} ~\rm{M}_\odot$ at $z = 3.493$, this galaxy is nearly three times more massive than the highest redshift spectroscopically confirmed absorption-line identified galaxy known. The star-formation history of this quiescent galaxy implies that it formed $>1000 ~\rm{M}_\odot$/yr for almost 0.5 Gyr beginning at $z\sim7.2$, strongly suggestive that it is the descendant of massive dusty star-forming galaxies at $5<z<7$ recently observed with ALMA. While galaxies with similarly extreme stellar masses are reproduced in some simulations at early times, such a lack of ongoing star formation is not seen there. This suggests the need for a more rapid quenching process than is currently prescribed, challenging our current understanding of how ultra-massive galaxies form and evolve in the early Universe.

We present 16 new ultrabright $H_{AB}\lesssim25$ galaxy candidates at z~8 identified over the COSMOS/UltraVISTA field. The new search takes advantage of the deepest-available ground-based optical and near-infrared observations, including the DR3 release of UltraVISTA and full-depth Spitzer/IRAC observations from the SMUVS and SPLASH programs. Candidates are selected using Lyman-break criteria, combined with strict optical non-detection and SED-fitting criteria, minimizing contamination by low-redshift galaxies and low-mass stars. HST/WFC3 coverage from the DASH program reveals that one source evident in our ground-based near-IR data has significant substructure and may actually correspond to 3 separate z~8 objects, resulting in a sample of 18 galaxies, 10 of which seem to be fairly robust (with a >97% probability of being at z>7). The UV-continuum slope $\beta$ for the bright z~8 sample is $\beta=-2.2\pm0.6$, bluer but still consistent with that of similarly bright galaxies at z~6 ($\beta=-1.55\pm0.17$) and z~7 ($\beta=-1.75\pm0.18$). Their typical stellar masses are 10$^{9.1^{+0.5}_{-0.4}}M_{\odot}$, with the SFRs of $32^{+44}_{-32}M_{\odot}$/year, specific SFR of $4^{+8}_{-4}$ Gyr$^{-1}$, stellar ages of $\sim22^{+69}_{-22}$\,Myr, and low dust content A$_V=0.15^{+0.30}_{-0.15}$ mag. Using this sample we constrain the bright end of the z~8 UV luminosity function (LF). When combined with recent empty field LF estimates at z~8-9, the resulting z~8 LF can be equally well represented by either a Schechter or a double power-law (DPL) form. Assuming a Schechter parameterization, the best-fit characteristic magnitude is $M^*= -20.95^{+0.30}_{-0.35}$ mag with a very steep faint end slope $\alpha=-2.15^{+0.20}_{-0.19}$. These new candidates include amongst the brightest yet found at these redshifts, 0.5-1.0 mag brighter than found over CANDELS, providing excellent targets for follow-up studies.

The ZFIRE survey has spectroscopically confirmed two proto-clusters using the MOSFIRE instrument on Keck 1: one at z=2.095 in COSMOS and another at z=1.62 in UDS. Here we use an updated ZFIRE dataset to derive the properties of ionized gas regions of proto-cluster galaxies by extracting fluxes from emission lines H$\beta$ 4861\AA, [OIII] 5007\AA, H$\alpha$ 6563\AA, [NII] 6585\AA, and [SII] 6716,6731\AA. We measure gas-phase metallicity of members in both proto-clusters using two indicators, including a strong-line indicator relatively independent of ionization parameter and electron density. Proto-cluster and field galaxies in both UDS and COSMOS lie on the same Mass-Metallicity Relation with both metallicity indicators. We compare our results to recent IllustrisTNG results, which reports no significant gas-phase metallicity offset between proto-cluster and field galaxies until z=1.5. This is in agreement with our observed metallicities, where no offset is measured between proto-cluster and field populations. We measure tentative evidence from stacked spectra that UDS high mass proto-cluster and field galaxies have differing [OIII]/H$\beta$ ratios, however these results are dependent on the sample size of the high mass stacks.

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 investigate the stellar and dust properties of massive (log$(M_*/M_\odot) \ge 10.5$) and dusty ($A_V \ge 1$) galaxies at $1 \le z \le 4$ by modeling their spectral energy distributions (SEDs) obtained from the combination of UltraVISTA DR3 photometry and \textit{Herschel} PACS-SPIRE data using MAGPHYS. Although the rest-frame U-V vs V-J (UVJ) diagram traces well the star-formation rates (SFR) and dust obscuration (A$_V$) out to $z \sim 3$, $\sim$15-20\% of the sample surprisingly resides in the quiescent region of the UVJ diagram, while $\sim50$\% at $3<z<4$ fall in the unobscured star-forming region. The median SED of massive dusty galaxies exhibits weaker MIR and UV emission, and redder UV slopes with increasing cosmic time. The IR emission for our sample has a significant contribution ($>20\%$) from dust heated by evolved stellar populations rather than star formation, demonstrating the need for panchromatic SED modeling. The local relation between dust mass and SFR is followed only by a sub-sample with cooler dust temperatures, while warmer objects have reduced dust masses at a given SFR. Most star-forming galaxies in our sample do not follow local IRX-$\beta$ relations, though IRX does strongly correlate with A$_V$. Our sample follows local relations, albeit with large scatter, between ISM diagnostics and sSFR. We show that FIR-detected sources represent the extreme of a continuous population of dusty galaxies rather than a fundamentally different population. Finally, using commonly adopted relations to derive SFRs from the combination of the rest-frame UV and the observed 24$\mu$m is found to overestimate the SFR by a factor of 3-5 for the galaxies in our sample.

The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$\beta$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H$\beta$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim8$ galaxies, such as the ionizing photon production efficiency with $\log(\xi_\mathrm{ion}/\mathrm{erg}^{-1}\hspace{1mm}\mathrm{Hz})\geq25.77$. We find a relatively tight correlation between the [OIII]+H$\beta$ and UV luminosity, which we use to derive for the first time the [OIII]+H$\beta$ luminosity function (LF) at $z\sim8$. The $z\sim8$ [OIII]+H$\beta$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H$\beta$ luminosity at a given UV luminosity from $z\sim3$ to $z\sim8$. Finally, using the [OIII]+H$\beta$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim30$ mag will be required.

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering.

Using high resolution spectra from the VLT LEGA-C program, we reconstruct the star formation histories (SFHs) of 607 galaxies at redshifts $z = 0.6-1.0$ and stellar masses $\gtrsim10^{10}$M$_{\odot}$ using a custom full spectrum fitting algorithm that incorporates the emcee and FSPS packages. We show that the mass-weighted age of a galaxy correlates strongly with stellar velocity dispersion ($\sigma_*$) and ongoing star-formation (SF) activity, with the stellar content in higher-$\sigma_*$ galaxies having formed earlier and faster. The SFHs of quiescent galaxies are generally consistent with passive evolution since their main SF epoch, but a minority show clear evidence of a rejuvenation event in their recent past. The mean age of stars in galaxies that are star-forming is generally significantly younger, with SF peaking after $z<1.5$ for almost all star-forming galaxies in the sample: many of these still have either constant or rising SFRs on timescales $>100$Myrs. This indicates that $z>2$ progenitors of $z\sim1$ star-forming galaxies are generally far less massive. Finally, despite considerable variance in the individual SFHs, we show that the current SF activity of massive galaxies ($>$L$_*$) at $z\sim1$ correlates with SF levels at least $3$Gyrs prior: SFHs retain `memory' on a large fraction of the Hubble time. Our results illustrate a novel approach to resolve the formation phase of galaxies, and in identifying their individual evolutionary paths, connects progenitors and descendants across cosmic time. This is uniquely enabled by the high-quality continuum spectroscopy provided by the LEGA-C survey.

We perform a kinematic and morphological analysis of 44 star-forming galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared spectroscopy from Keck/MOSFIRE and F160W imaging from CANDELS/3D-HST as part of the ZFIRE survey. Our sample consists of cluster and field galaxies from $2.0 < z < 2.5$ with K band multi-object slit spectroscopic measurements of their H$\alpha$ emission lines. H$\alpha$ rotational velocities and gas velocity dispersions are measured using the Heidelberg Emission Line Algorithm (HELA), which compares directly to simulated 3D data-cubes. Using a suite of simulated emission lines, we determine that HELA reliably recovers input S$_{0.5}$ and angular momentum at small offsets, but $V_{2.2}/\sigma_g$ values are offset and highly scattered. We examine the role of regular and irregular morphology in the stellar mass kinematic scaling relations, deriving the kinematic measurement S$_{0.5}$, and finding $\log(S_{0.5}) = (0.38\pm0.07)\log(M/M_{\odot}-10) + (2.04\pm0.03)$ with no significant offset between morphological populations and similar levels of scatter ($\sim0.16$ dex). Additionally, we identify a correlation between M$_{\star}$ and $V_{2.2}/\sigma_g$ for the total sample, showing an increasing level of rotation dominance with increasing M$_{\star}$, and a high level of scatter for both regular and irregular galaxies. We estimate the specific angular momenta ($j_{disk}$) of these galaxies and find a slope of $0.36\pm0.12$, shallower than predicted without mass-dependent disk growth, but this result is possibly due to measurement uncertainty at M$_{\star}$ $<$ 9.5. However, through a K-S test we find irregular galaxies to have marginally higher $j_{disk}$ values than regular galaxies, and high scatter at low masses in both populations.

We present stellar rotation curves and velocity dispersion profiles for 104 quiescent galaxies at $z=0.6-1$ from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey. Rotation is typically probed across 10-20kpc, or to an average of 2.7${\rm R_e}$. Combined with central stellar velocity dispersions ($\sigma_0$) this provides the first determination of the dynamical state of a sample selected by a lack of star formation activity at large lookback time. The most massive galaxies ($M_{\star}>2\times10^{11}\,M_{\odot}$) generally show no or little rotation measured at 5kpc ($|V_5|/\sigma_0<0.2$ in 8 of 10 cases), while ${\sim}64\%$ of less massive galaxies show significant rotation. This is reminiscent of local fast- and slow-rotating ellipticals and implies that low- and high-redshift quiescent galaxies have qualitatively similar dynamical structures. We compare $|V_5|/\sigma_0$ distributions at $z\sim0.8$ and the present day by re-binning and smoothing the kinematic maps of 91 low-redshift quiescent galaxies from the CALIFA survey and find evidence for a decrease in rotational support since $z\sim1$. This result is especially strong when galaxies are compared at fixed velocity dispersion; if velocity dispersion does not evolve for individual galaxies then the rotational velocity at 5kpc was an average of ${94\pm22\%}$ higher in $z\sim0.8$ quiescent galaxies than today. Considering that the number of quiescent galaxies grows with time and that new additions to the population descend from rotationally-supported star-forming galaxies, our results imply that quiescent galaxies must lose angular momentum between $z\sim1$ and the present, presumably through dissipationless merging, and/or that the mechanism that transforms star-forming galaxies also reduces their rotational support.

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang.

We obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and Hyde, respectively. No dust emission is detected at the location of Jekyll, implying SFR < 13 Msun/yr, which is the most stringent upper limit ever obtained for a quiescent galaxy at these redshifts. The two sources are confirmed to lie at the same redshift thanks to the detection of [CII]158 in Hyde, which provides one the few robust redshifts for an "H-dropout" galaxy. The line has a rotating-disk velocity profile blueshifted from Jekyll by 549+/-60 km/s, demonstrating that it is produced by another galaxy. Careful de-blending of the Spitzer imaging confirms the existence of Hyde, and its non-detection with Hubble requires extreme attenuation by dust. Modeling the photometry of both galaxies shows that Jekyll has fully quenched >200 Myr prior to observation and still presents a challenge for models, while Hyde only harbors moderate star-formation (SFR<120 Msun/yr) and is located at least a factor 1.4 below the z~4 main sequence. Hyde could also have stopped forming stars <200 Myr before being observed, which would be consistent with its hight compactness similar to z~4 quiescent galaxies and its low [CII]/FIR ratio, but significant SF cannot be ruled out. Finally, we show that Hyde hosts a dense reservoir of gas comparable to that of extreme starbursts, suggesting that its SFR was reduced without expelling the gas outside of the galaxy. We argue that Jekyll and Hyde can be seen as two stages of quenching, and provide a unique laboratory to study this poorly understood phenomenon. [abridged]

In the early Universe finding massive galaxies that have stopped forming stars present an observational challenge as their rest-frame ultraviolet emission is negligible and they can only be reliably identified by extremely deep near-infrared surveys. These have revealed the presence of massive, quiescent early-type galaxies appearing in the universe as early as z$\sim$2, an epoch 3 Gyr after the Big Bang. Their age and formation processes have now been explained by an improved generation of galaxy formation models where they form rapidly at z$\sim$3-4, consistent with the typical masses and ages derived from their observations. Deeper surveys have now reported evidence for populations of massive, quiescent galaxies at even higher redshifts and earlier times, however the evidence for their existence, and redshift, has relied entirely on coarsely sampled photometry. These early massive, quiescent galaxies are not predicted by the latest generation of theoretical models. Here, we report the spectroscopic confirmation of one of these galaxies at redshift z=3.717 with a stellar mass of 1.7$\times$10$^{11}$ M$_\odot$ whose absorption line spectrum shows no current star-formation and which has a derived age of nearly half the age of the Universe at this redshift. The observations demonstrates that the galaxy must have quickly formed the majority of its stars within the first billion years of cosmic history in an extreme and short starburst. This ancestral event is similar to those starting to be found by sub-mm wavelength surveys pointing to a possible connection between these two populations. Early formation of such massive systems is likely to require significant revisions to our picture of early galaxy assembly.

We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5\sigma$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $\beta=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form.

We study galactic star-formation activity as a function of environment and stellar mass over 0.5<z<2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey. We estimate the galaxy environment using a Bayesian-motivated measure of the distance to the third nearest neighbor for galaxies to the stellar mass completeness of our survey, $\log(M/M_\odot)>9 (9.5)$ at z=1.3 (2.0). This method, when applied to a mock catalog with the photometric-redshift precision ($\sigma_z / (1+z) \lesssim 0.02$), recovers galaxies in low- and high-density environments accurately. We quantify the environmental quenching efficiency, and show that at z> 0.5 it depends on galaxy stellar mass, demonstrating that the effects of quenching related to (stellar) mass and environment are not separable. In high-density environments, the mass and environmental quenching efficiencies are comparable for massive galaxies ($\log (M/M_\odot)\gtrsim$ 10.5) at all redshifts. For lower mass galaxies ($\log (M/M)_\odot) \lesssim$ 10), the environmental quenching efficiency is very low at $z\gtrsim$ 1.5, but increases rapidly with decreasing redshift. Environmental quenching can account for nearly all quiescent lower mass galaxies ($\log(M/M_\odot) \sim$ 9-10), which appear primarily at $z\lesssim$ 1.0. The morphologies of lower mass quiescent galaxies are inconsistent with those expected of recently quenched star-forming galaxies. Some environmental process must transform the morphologies on similar timescales as the environmental quenching itself. The evolution of the environmental quenching favors models that combine gas starvation (as galaxies become satellites) with gas exhaustion through star-formation and outflows ("overconsumption"), and additional processes such as galaxy interactions, tidal stripping and disk fading to account for the morphological differences between the quiescent and star-forming galaxy populations.

We investigate the stellar kinematics and stellar populations of 58 radio-loud galaxies of intermediate luminosities (L$_{3 GHz}$ $>$ 10$^{23}$ W Hz$^{-1}$ ) at 0.6 < z < 1. This sample is constructed by cross-matching galaxies from the deep VLT/VIMOS LEGA-C spectroscopic survey with the VLA 3 GHz dataset. The LEGA-C continuum spectra reveal for the first time stellar velocity dispersions and age indicators of z $\sim$ 1 radio galaxies. We find that $z\sim 1$ radio-loud AGN occur exclusively in predominantly old galaxies with high velocity dispersions: $\sigma_*>$ 175 km s$^{-1}$, corresponding to black hole masses in excess of $10^8$ M$_{\odot}$. Furthermore, we confirm that at a fixed stellar mass the fraction of radio-loud AGN at z $\sim$ 1 is 5 - 10 times higher than in the local universe, suggesting that quiescent, massive galaxies at z $\sim$ 1 switch on as radio AGN on average once every Gyr. Our results strengthen the existing evidence for a link between high black-hole masses, radio loudness and quiescence at z $\sim$ 1.

We present the first determination of the galaxy luminosity function (LF) at z~4, 5, 6 and 7 in the rest-frame optical at lambda_rest~900 nm (z' band). The rest-frame optical light traces the content in low-mass evolved stars (~stellar mass - Mstar), minimizing potential measurement biases for Mstar. Moreover it is less affected by nebular line emission contamination and dust attenuation, is independent of stellar population models, and can be probed up to z~8 through Spitzer/IRAC. Our analysis leverages the unique full depth Spitzer/IRAC 3.6um-to-8.0um data over the CANDELS/GOODS-N, CANDELS/GOODS-S and COSMOS/UltraVISTA fields. We find that at absolute magnitudes M_z' fainter than >-23 mag, M_z' linearly correlates with M_UV,1600. At brighter M_z', M_UV,1600 presents a turnover, suggesting that the stellar mass-to-light ratio Mstar/L_UV,1600 could be characterised by a very broad range of values at high stellar masses. Median-stacking analysis recovers a Mstar/L_z' roughly independent on M_z' for M_z'>-23 mag, but exponentially increasing at brighter magnitudes. We find that the evolution of the LF marginally prefers a pure evolution in luminosity over a pure evolution in density, with the characteristic luminosity decreasing by a factor ~5x between z~4 and z~7. Direct application of the recovered Mstar/L_z' generates stellar mass functions consistent with average measurements from the literature. Measurements of the stellar-to-halo mass ratio at fixed cumulative number density show that it is roughly constant with redshift for Mh>10^12Msun. This is also supported by the fact that the evolution of the LF at 4<z<7 can be accounted for by a rigid displacement in luminosity corresponding to the evolution of the halo mass from abundance matching.

Using deep multi-wavelength photometry of galaxies from ZFOURGE, we group galaxies at $2.5<z<4.0$ by the shape of their spectral energy distributions (SEDs). We identify a population of galaxies with excess emission in the $K_s$-band, which corresponds to [OIII]+H$\beta$ emission at $2.95<z<3.65$. This population includes 78% of the bluest galaxies with UV slopes steeper than $\beta = -2$. We de-redshift and scale this photometry to build two composite SEDs, enabling us to measure equivalent widths of these Extreme [OIII]+H$\beta$ Emission Line Galaxies (EELGs) at $z\sim3.5$. We identify 60 galaxies that comprise a composite SED with [OIII]+H$\beta$ rest-frame equivalent width of $803\pm228$\AA\ and another 218 galaxies in a composite SED with equivalent width of $230\pm90$\AA. These EELGs are analogous to the `green peas' found in the SDSS, and are thought to be undergoing their first burst of star formation due to their blue colors ($\beta < -1.6$), young ages ($\log(\rm{age}/yr)\sim7.2$), and low dust attenuation values. Their strong nebular emission lines and compact sizes (typically $\sim1.4$ kpc) are consistent with the properties of the star-forming galaxies possibly responsible for reionizing the universe at $z>6$. Many of the EELGs also exhibit Lyman-$\alpha$ emission. Additionally, we find that many of these sources are clustered in an overdensity in the Chandra Deep Field South, with five spectroscopically confirmed members at $z=3.474 \pm 0.004$. The spatial distribution and photometric redshifts of the ZFOURGE population further confirm the overdensity highlighted by the EELGs.

Using observations made with MOSFIRE on Keck I as part of the ZFIRE survey, we present the stellar mass Tully-Fisher relation at 2.0 < z < 2.5. The sample was drawn from a stellar mass limited, Ks-band selected catalog from ZFOURGE over the CANDELS area in the COSMOS field. We model the shear of the Halpha emission line to derive rotational velocities at 2.2X the scale radius of an exponential disk (V2.2). We correct for the blurring effect of a two-dimensional PSF and the fact that the MOSFIRE PSF is better approximated by a Moffat than a Gaussian, which is more typically assumed for natural seeing. We find for the Tully-Fisher relation at 2.0 < z < 2.5 that logV2.2 =(2.18 +/- 0.051)+(0.193 +/- 0.108)(logM/Msun - 10) and infer an evolution of the zeropoint of Delta M/Msun = -0.25 +/- 0.16 dex or Delta M/Msun = -0.39 +/- 0.21 dex compared to z = 0 when adopting a fixed slope of 0.29 or 1/4.5, respectively. We also derive the alternative kinematic estimator S0.5, with a best-fit relation logS0.5 =(2.06 +/- 0.032)+(0.211 +/- 0.086)(logM/Msun - 10), and infer an evolution of Delta M/Msun= -0.45 +/- 0.13 dex compared to z < 1.2 if we adopt a fixed slope. We investigate and review various systematics, ranging from PSF effects, projection effects, systematics related to stellar mass derivation, selection biases and slope. We find that discrepancies between the various literature values are reduced when taking these into account. Our observations correspond well with the gradual evolution predicted by semi-analytic models.

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 introduce the ALMA Redshift 4 Survey (AR4S), a systematic ALMA survey of all the known galaxies with stellar mass (M*) larger than 5e10 Msun at 3.5<z<5 in the GOODS--south, UDS and COSMOS CANDELS fields. The sample we have analyzed in this paper is composed of 96 galaxies observed with ALMA at 890um (180um rest-frame) with an on-source integration time of 1.3 min per galaxy. We detected 32% of the sample at more than 3 sigma significance. Using the stacked ALMA and Herschel photometry, we derived an average dust temperature of 40+/-2 K for the whole sample, and extrapolate the Lir and SFR for all our galaxies based on their ALMA flux. We then used a forward modeling approach to estimate their intrinsic sSFR distribution, deconvolved of measurement errors and selection effects: we find a linear relation between SFR and M*, with a median sSFR=2.8+/-0.8 Gyr and a dispersion around that relation of 0.28+/-0.13 dex. This latter value is consistent with that measured at lower redshifts, which is proof that the main sequence of star-forming galaxies was already in place at z=4, at least among massive galaxies. These new constraints on the properties of the main sequence are in good agreement with the latest predictions from numerical simulations, and suggest that the bulk of star formation in galaxies is driven by the same mechanism from z=4 to the present day, that is, over at least 90% of the cosmic history. We also discuss the consequences of our results on the population of early quiescent galaxies. This paper is part of a series that will employ these new ALMA observations to explore the star formation and dust properties of the massive end of the z=4 galaxy population.

We compare galaxy scaling relations as a function of environment at $z\sim2$ with our ZFIRE survey where we have measured H$\alpha$ fluxes for 90 star-forming galaxies selected from a mass-limited [$\log(M_{\star}/M_{\odot})>9$] sample based on ZFOURGE. The cluster galaxies (37) are part of a confirmed system at z=2.095 and the field galaxies (53) are at $1.9<z<2.4$; all are in the COSMOS legacy field. There is no statistical difference between H$\alpha$-emitting cluster and field populations when comparing their star formation rate (SFR), stellar mass ($M_{\star}$), galaxy size ($r_{eff}$), SFR surface density [$\Sigma$(H$\alpha_{star}$)], and stellar age distributions. The only difference is that at fixed stellar mass, the H$\alpha$-emitting cluster galaxies are $\log(r_{eff})\sim0.1$ larger than in the field. Approximately 19% of the H$\alpha$-emitters in the cluster and 26% in the field are IR-luminous ($L_{IR}>2\times10^{11} L_{\odot}$). Because the LIRGs in our combined sample are $\sim5$ times more massive than the low-IR galaxies, their radii are $\sim70$% larger. To track stellar growth, we separate galaxies into those that lie above, on, and below the H$\alpha$ star-forming main sequence (SFMS) using $\Delta$SFR$(M_{\star})=\pm0.2$ dex. Galaxies above the SFMS (starbursts) tend to have higher H$\alpha$ SFR surface densities and younger light-weighted stellar ages compared to galaxies below the SFMS. Our results indicate that starbursts (+SFMS) in the cluster and field at $z\sim2$ are growing their stellar cores. Lastly, we compare to the (SFR-$M_{\star}$) relation from RHAPSODY cluster simulations and find the predicted slope is nominally consistent with the observations. However, the predicted cluster SFRs tend to be too low by a factor of $\sim2$ which seems to be a common problem for simulations across environment.

We derive H{\alpha} fluxes for a large spectroscopic and photometric-redshift-selected sample of sources over GOODS-North and South in the redshift range z = 3.8-5.0 with deep HST, Spitzer/IRAC, and ground-based observations. The H{\alpha} flux is inferred based on the offset between the IRAC 3.6 {\mu}m flux and that predicted from the best-fit SED. We demonstrate that the H{\alpha} flux correlates well with dust- corrected UV star-formation rate (SFR) and therefore can serve as an independent SFR indicator. However, we also find a systematic offset in the SFR_H{\alpha}/SFR_UV ratios for z ~ 4-5 galaxies relative to local relations (assuming the same dust corrections for nebular regions and stellar light). We show that we can resolve the modest tension in the inferred SFRs by assuming bluer intrinsic UV slopes (increasing the dust correction), a rising star-formation history or assuming a low metallicity stellar population with a hard ionizing spectrum (increasing the L_H{\alpha}/SFR ratio). Using H{\alpha} as a SFR indicator, we find a normalization of the star formation main sequence in good agreement with recent SED-based determinations and also derive the SFR functions at z ~ 4-8. In addition, we assess for the first time the burstiness of star formation in z ~ 4 galaxies on <100 Myr time scales by comparing UV and H{\alpha}-based sSFRs; their one-to-one relationship argues against significantly bursty star-formation histories. Further progress will be made on these results, by incorporating new results from ALMA to constrain the dust-obscured star formation in high-redshift UV-selected samples.

The gas accretion and star-formation histories of galaxies like the Milky Way remain an outstanding problem in astrophysics. Observations show that 8 billion years ago, the progenitors to Milky Way-mass galaxies were forming stars 30 times faster than today and predicted to be rich in molecular gas, in contrast with low present-day gas fractions ($<$10%). Here we show detections of molecular gas from the CO(J=3-2) emission (rest-frame 345.8 GHz) in galaxies at redshifts z=1.2-1.3, selected to have the stellar mass and star-formation rate of the progenitors of today's Milky Way-mass galaxies. The CO emission reveals large molecular gas masses, comparable to or exceeding the galaxy stellar masses, and implying most of the baryons are in cold gas, not stars. The galaxies' total luminosities from star formation and CO luminosities yield long gas-consumption timescales. Compared to local spiral galaxies, the star-formation efficiency, estimated from the ratio of total IR luminosity to CO emission,} has remained nearly constant since redshift z=1.2, despite the order of magnitude decrease in gas fraction, consistent with results for other galaxies at this epoch. Therefore the physical processes that determine the rate at which gas cools to form stars in distant galaxies appear to be similar to that in local galaxies.

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

We perform a kinematic analysis of galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared (NIR) spectroscopy from Keck/MOSFIRE as part of the ZFIRE survey. Our sample consists of 75 Ks-band selected star-forming galaxies from the ZFOURGE survey with stellar masses ranging from log(M$_{\star}$/M$_{\odot}$)$=9.0-11.0$, 28 of which are members of a known overdensity at $z=2.095$. We measure H$\alpha$ emission-line integrated velocity dispersions ($\sigma_{\rm int}$) from 50$-$230 km s$^{-1}$, consistent with other emission-line studies of $z\sim2$ field galaxies. From these data we estimate virial, stellar, and gas masses and derive correlations between these properties for cluster and field galaxies at $z\sim2$. We find evidence that baryons dominate within the central effective radius. However, we find no statistically significant differences between the cluster and the field, and conclude that the kinematics of star-forming galaxies at $z\sim2$ are not significantly different between the cluster and field environments.

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

We investigate active galactic nuclei (AGN) candidates within the FourStar Galaxy Evolution Survey (ZFOURGE) to determine the impact they have on star-formation in their host galaxies. We first identify a population of radio, X-ray, and infrared-selected AGN by cross-matching the deep $K_{s}$-band imaging of ZFOURGE with overlapping multi-wavelength data. From this, we construct a mass-complete (log(M$_{*}$/M$_{\odot}$) $\ge$ 9.75), AGN luminosity limited sample of 235 AGN hosts over z = 0.2 - 3.2. We compare the rest-frame U - V versus V - J (UVJ) colours and specific star-formation rates (sSFRs) of the AGN hosts to a mass-matched control sample of inactive (non-AGN) galaxies. UVJ diagnostics reveal AGN tend to be hosted in a lower fraction of quiescent galaxies and a higher fraction of dusty galaxies than the control sample. Using 160{\mu}m Herschel PACS data, we find the mean specific star-formation rate of AGN hosts to be elevated by 0.34$\pm$0.07 dex with respect to the control sample across all redshifts. This offset is primarily driven by infrared-selected AGN, where the mean sSFR is found to be elevated by as much as a factor of ~5. The remaining population, comprised predominantly of X-ray AGN hosts, is found mostly consistent with inactive galaxies, exhibiting only a marginal elevation. We discuss scenarios that may explain these findings and postulate that AGN are less likely to be a dominant mechanism for moderating galaxy growth via quenching than has previously been suggested.

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3\sigma$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

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.

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.

There is ongoing debate regarding the extent that environment affects galaxy size growth beyond z>1. To investigate the differences in star-forming and quiescent galaxy properties as a function of environment at z=2.1, we create a mass-complete sample of 59 cluster galaxies Spitler et al. (2012) and 478 field galaxies with log(M)>9 using photometric redshifts from the ZFOURGE survey. We compare the mass-size relation of field and cluster galaxies using measured galaxy semi-major axis half-light radii ($r_{1/2,maj}$) from CANDELS HST/F160W imaging. We find consistent mass normalized (log(M)=10.7) sizes for quiescent field galaxies ($r_{1/2,maj}=1.81\pm0.29$ kpc) and quiescent cluster galaxies ($r_{1/2,maj}=2.17\pm0.63$ kpc). The mass normalized size of star-forming cluster galaxies ($r_{1/2,maj}=4.00\pm0.26$ kpc ) is 12% larger (KS test $2.1\sigma$) than star-forming field galaxies ($r_{1/2,maj}=3.57\pm0.10$ kpc). From the mass-color relation we find that quiescent field galaxies with 9.7<log(M)<10.4 are slightly redder (KS test $3.6\sigma$) than quiescent cluster galaxies, while cluster and field quiescent galaxies with log(M)>10.4 have consistent colors. We find that star-forming cluster galaxies are on average 20% redder than star-forming field galaxies at all masses. Furthermore, we stack galaxy images to measure average radial color profiles as a function of mass. Negative color gradients are only present for massive star-forming field and cluster galaxies with log(M)>10.4, the remaining galaxy masses and types have flat profiles. Our results suggest given the observed differences in size and color of star-forming field and cluster galaxies, that the environment has begun to influence/accelerate their evolution. However, the lack of differences between field and cluster quiescent galaxies indicates that the environment has not begun to significantly influence their evolution at z~2.

We investigate the environmental dependence of the mass-metallicity relation at z=2 with MOSFIRE/Keck as part of the ZFIRE survey. Here, we present the chemical abundance of a Virgo-like progenitor at z=2.095 that has an established red sequence. We identified 43 cluster ($<z>=2.095\pm0.004$) and 74 field galaxies ($<z>=2.195\pm0.083$) for which we can measure metallicities. For the first time, we show that there is no discernible difference between the mass-metallicity relation of field and cluster galaxies to within 0.02dex. Both our field and cluster galaxy mass-metallicity relations are consistent with recent field galaxy studies at z~2. We present hydrodynamical simulations for which we derive mass-metallicity relations for field and cluster galaxies. We find at most a 0.1dex offset towards more metal-rich simulated cluster galaxies. Our results from both simulations and observations are suggestive that environmental effects, if present, are small and are secondary to the ongoing inflow and outflow processes that are governed by galaxy halo mass.

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

We present the spectroscopic confirmation of a galaxy cluster at $z=2.095$ in the COSMOS field. This galaxy cluster was first reported in the ZFOURGE survey as harboring evolved massive galaxies using photometric redshifts derived with deep near-infrared (NIR) medium-band filters. We obtain medium resolution ($R \sim$ 3600) NIR spectroscopy with MOSFIRE on the Keck 1 telescope and secure 180 redshifts in a $12'\times12'$ region. We find a prominent spike of 57 galaxies at $z=2.095$ corresponding to the galaxy cluster. The cluster velocity dispersion is measured to be $\sigma_{\rm v1D}$ = 552 $\pm$ 52 km/s. This is the first study of a galaxy cluster in this redshift range ($z \gt 2.0$) with the combination of spectral resolution ($\sim$26 km/s) and the number of confirmed members (${>}50$) needed to impose a meaningful constraint on the cluster velocity dispersion and map its members over a large field of view. Our $\Lambda$CDM cosmological simulation suggests that this cluster will most likely evolve into a Virgo-like cluster with ${\rm M_{vir}}{=}10^{14.4\pm0.3} {\rm M_\odot}$ ($68\%$ confidence) at $z\sim$ 0. The theoretical expectation of finding such a cluster is $\sim$ $4\%$. Our results demonstrate the feasibility of studying galaxy clusters at $z > 2$ in the same detailed manner using multi-object NIR spectrographs as has been done in the optical in lower redshift clusters.

We report the likely identification of a substantial population of massive M~10^11M_Sun galaxies at z~4 with suppressed star formation rates (SFRs), selected on rest-frame optical to near-IR colors from the FourStar Galaxy Evolution Survey. The observed spectral energy distributions show pronounced breaks, sampled by a set of near-IR medium-bandwidth filters, resulting in tightly constrained photometric redshifts. Fitting stellar population models suggests large Balmer/4000AA breaks, relatively old stellar populations, large stellar masses and low SFRs, with a median specific SFR of 2.9+/-1.8 x 10^-11/yr. Ultradeep Herschel/PACS 100micron, 160micron and Spitzer/MIPS 24micron data reveal no dust-obscured SFR activity for 15/19 (79%) galaxies. Two far-IR detected galaxies are obscured QSOs. Stacking the far-IR undetected galaxies yields no detection, consistent with the SED fit, indicating independently that the average specific SFR is at least 10x smaller than of typical star-forming galaxies at z~4. Assuming all far-IR undetected galaxies are indeed quiescent, the volume density is 1.8+/-0.7 x 10^-5Mpc^-3 to a limit of log10M/M_Sun>10.6, which is 10x and 80x lower than at z = 2 and z = 0.1. They comprise a remarkably high fraction (~35%) of z~4 massive galaxies, suggesting that suppression of star formation was efficient even at very high redshift. Given the average stellar age of 0.8Gyr and stellar mass of 0.8x10^11M_Sun, the galaxies likely started forming stars before z =5, with SFRs well in excess of 100M_Sun/yr, far exceeding that of similarly abundant UV-bright galaxies at z>4. This suggests that most of the star-formation in the progenitors of quiescent z~4 galaxies was obscured by dust.

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 combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-established z~4-7 UV LFs, UV-continuum slope trends with redshift and luminosity, and IRX-beta relations. They are well-described by Schechter relations. We extend the comparison baseline for SFR functions to z~2 by considering recent determinations of the H{\alpha} and mid-IR luminosity functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, \Delta\alpha\sim+0.13, and show no clear evolution with cosmic time (z~0-7). In addition, we find that the characteristic value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth as a function of redshift from ~5 M_sun/yr at z~8, 650 Myr after the Big Bang, to ~100 M_sun/yr at z~2, ~2.5 Gyr later. Recent results at z~10, close to the onset of galaxy formation, are consistent with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the redshift range z~2-8, providing validation for our dust corrections. These results provide strong evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time.

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 report the first results from the Z-FOURGE survey: the discovery of a candidate galaxy cluster at z = 2.2 consisting of two compact overdensities with red galaxies detected at >20{\sigma} above the mean surface density. The discovery was made possible by a new deep (Ks < 24.8 AB 5{\sigma}) Magellan/FourStar near-IR imaging survey with 5 custom medium-bandwidth filters. The filters pinpoint the location of the Balmer/4000\AA\ break in evolved stellar populations at 1.5 < z < 3.5, yielding significantly more accurate photometric redshifts than possible with broadband imaging alone. The overdensities are within 1' of each other in the COSMOS field and appear to be embedded in a larger structure that contains at least one additional overdensity (\sim10{\sigma}). Considering the global properties of the overdensities, the z = 2.2 system appears to be the most distant example of a galaxy cluster with a population of red galaxies. A comparison to a large {\Lambda}CDM simulation suggests that the system may consist of merging subclusters, with properties in between those of z > 2 protoclusters with more diffuse distributions of blue galaxies and the lower-redshift galaxy clusters with prominent red sequences. The structure is completely absent in public optical catalogs in COSMOS and only weakly visible in a shallower near-IR survey. The discovery showcases the potential of deep near-IR surveys with medium-band filters to advance the understanding of environment and galaxy evolution at z > 1.5.

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

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

We investigate the star formation history of the universe using FIREWORKS, a multiwavelength survey of the CDFS. We study the evolution of the specific star formation rate (sSFR) with redshift in different mass bins from z = 0 to z ~ 3. We find that the sSFR increases with redshift for all masses. The logarithmic increase of the sSFR with redshift is nearly independent of mass, but this cannot yet be verified at the lowest-mass bins at z > 0.8, due to incompleteness. We convert the sSFRs to a dimensionless growth rate to facilitate a comparison with a semi-analytic galaxy formation model that was implemented on the Millennium Simulation. The model predicts that the growth rates and sSFRs increase similarly with redshift for all masses, consistent with the observations. However, we find that for all masses, the inferred observed growth rates increase more rapidly with redshift than the model predictions. We discuss several possible causes for this discrepancy, ranging from field-to-field variance, conversions to SFR, and shape of the IMF. We find that none of these can solve the discrepancy completely. We conclude that the models need to be adapted to produce the steep increase in growth rate between redshift z=0 and z=1.

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

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

We present the mid-infrared colors of X-ray-detected AGN and explore mid-infrared selection criteria. Using a statistical matching technique, the likelihood ratio, over 900 IRAC counterparts were identified with a new MUSYC X-ray source catalog that includes ~1000 published X-ray sources in the Chandra Deep Field-South and Extended Chandra Deep Field-South. Most X-ray-selected AGN have IRAC spectral shapes consistent with power-law slopes, f_{nu} ~ nu^{alpha}, and display a wide range of colors, -2 < alpha < 2. Although X-ray sources typically fit to redder (more negative alpha) power-laws than non-X-ray detected galaxies, more than 50% do have flat or blue (galaxy-like) spectral shapes in the observed 3-8 micron band. Only a quarter of the X-ray selected AGN detected at 24 micron are well fit by featureless red power laws in the observed 3.6-24 micron, likely the subset of our sample whose infrared spectra are dominated by emission from the central AGN region. Most IRAC color-selection criteria fail to identify the majority of X-ray-selected AGN, finding only the more luminous AGN, the majority of which have broad emission lines. In deep surveys, these color-selection criteria select 10-20% of the entire galaxy population and miss many moderate luminosity AGN.