We use deep \textit{Hubble Space Telescope} spectroscopy to constrain the metallicities and (\editone{light-weighted}) ages of massive ($\log M_\ast/M_\odot\gtrsim10$) galaxies selected to have quiescent stellar populations at $1.0<z<1.8$. The data include 12--orbit depth coverage with the WFC3/G102 grism covering $\sim$ $8,000<\lambda<11,500$~\AA\, at a spectral resolution of $R\sim 210$ taken as part of the CANDELS Lyman-$\alpha$ Emission at Reionization (CLEAR) survey. At $1.0<z<1.8$, the spectra cover important stellar population features in the rest-frame optical. We simulate a suite of stellar population models at the grism resolution, fit these to the data for each galaxy, and derive posterior likelihood distributions for metallicity and age. We stack the posteriors for subgroups of galaxies in different redshift ranges that include different combinations of stellar absorption features. Our results give \editone{light-weighted ages of $t_{z \sim 1.1}= 3.2\pm 0.7$~Gyr, $t_{z \sim 1.2}= 2.2\pm 0.6$~Gyr, $t_{z\sim1.3}= 3.1\pm 0.6$~Gyr, and $t_{z\sim1.6}= 2.0 \pm 0.6$~Gyr, \editone{for galaxies at $z\sim 1.1$, 1.2, 1.3, and 1.6. This} implies that most of the massive quiescent galaxies at $1<z<1.8$ had formed $>68$\% of their stellar mass by a redshift of $z>2$}. The posteriors give metallicities of \editone{$Z_{z\sim1.1}=1.16 \pm 0.29$~$Z_\odot$, $Z_{z\sim1.2}=1.05 \pm 0.34$~$Z_\odot$, $Z_{z\sim1.3}=1.00 \pm 0.31$~$Z_\odot$, and $Z_{z\sim1.6}=0.95 \pm 0.39$~$Z_\odot$}. This is evidence that massive galaxies had enriched rapidly to approximately Solar metallicities as early as $z\sim3$.

We study the effects of galaxy environment on the evolution of the stellar-mass function (SMF) over 0.2 < z < 2.0 using the FourStar Galaxy Evolution (ZFOURGE) survey and NEWFIRM Medium-Band Survey (NMBS) down to the stellar-mass completeness limit, log M / Msun > 9.0 (9.5) at z = 1.0 (2.0). We compare the SMFs for quiescent and star-forming galaxies in the highest and lowest environments using a density estimator based on the distance to the galaxies' third-nearest neighbors. For star-forming galaxies, at all redshifts there are only minor differences with environment in the shape of the SMF. For quiescent galaxies, the SMF in the lowest densities shows no evolution with redshift, other than an overall increase in number density (phi*) with time. This suggests that the stellar-mass dependence of quenching in relatively isolated galaxies is both universal and does not evolve strongly. While at z >~ 1.5 the SMF of quiescent galaxies is indistinguishable in the highest and lowest densities, at lower redshifts it shows a rapidly increasing number density of lower-mass galaxies, log M / Msun ~= 9-10. We argue this evolution can account for all the redshift evolution in the shape of the total quiescent-galaxy SMF. This evolution in the quiescent-galaxy SMF at higher redshift (z > 1) requires an environmental-quenching efficiency that decreases with decreasing stellar mass at 0.5 < z < 1.5 or it would overproduce the number of lower-mass quiescent galaxies in denser environments. This requires a dominant environment process such as starvation combined with rapid gas depletion and ejection at z > 0.5 - 1.0 for galaxies in our mass range. The efficiency of this process decreases with redshift allowing other processes (such as galaxy interactions and ram-pressure stripping) to become more important at later times, z < 0.5.

We investigate the relationship between the black hole accretion rate (BHAR) and star-formation rate (SFR) for Milky Way (MW) and Andromeda (M31)-mass progenitors from z = 0.2 - 2.5. We source galaxies from the Ks-band selected ZFOURGE survey, which includes multi-wavelenth data spanning 0.3 - 160um. We use decomposition software to split the observed SEDs of our galaxies into their active galactic nuclei (AGN) and star-forming components, which allows us to estimate BHARs and SFRs from the infrared (IR). We perform tests to check the robustness of these estimates, including a comparison to BHARs and SFRs derived from X-ray stacking and far-IR analysis, respectively. We find as the progenit- ors evolve, their relative black hole-galaxy growth (i.e. their BHAR/SFR ratio) increases from low to high redshift. The MW-mass progenitors exhibit a log-log slope of 0.64 +/- 0.11, while the M31-mass progenitors are 0.39 +/- 0.08. This result contrasts with previous studies that find an almost flat slope when adopting X-ray/AGN-selected or mass-limited samples and is likely due to their use of a broad mixture of galaxies with different evolutionary histories. Our use of progenitor-matched samples highlights the potential importance of carefully selecting progenitors when searching for evolutionary relationships between BHAR/SFRs. Additionally, our finding that BHAR/SFR ratios do not track the rate at which progenitors quench casts doubts over the idea that the suppression of star-formation is predominantly driven by luminous AGN feedback (i.e. high BHARs).

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 use Hubble Space Telescope imaging to study the structural properties of ten of the most massive ($M \geq 10^{11.25}$ Msun) quiescent galaxies (QGs) in the UKIDSS UDS at $2.5<z<3.0$. The low spatial density of these galaxies required targeted WFC3 $H_{160}$ imaging, as such systems are rare in existing surveys like CANDELS. We fit Sersic models to the 2D light profiles and find that the median half-light radius is $R_e \sim 3$ kpc, a factor of $\sim 3$ smaller than QGs with similar masses at $z \sim 0$. Complementing our sample with similarly massive QGs at lower redshifts, we find that the median size evolves as $R_e \propto H(z)^{-0.85 \pm 0.12}$ (or alternatively, $R_e \propto (1+z)^{-0.90 \pm 0.12}$). This rate of evolution is slower than that for lower mass QGs. When compared to low redshift QGs, the axis ratio distribution for our high redshift massive QG sample is most consistent with those in which spheroids are dominant. These observations point to earlier size growth among massive QGs that also resulted in spheroidal systems. Finally, we measured residual-corrected surface brightness profiles for our sample. These show that the Sersic parameterization is generally representative out to several effective radii and does not miss excess low surface brightness light. The sizes inferred from the light profiles therefore confirm the compactness of these most massive high redshift QGs.

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.

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.

For the first time, we present the size evolution of a mass-complete (log(M*/Msol)>10) sample of star-forming galaxies over redshifts z=1-7, selected from the FourStar Galaxy Evolution Survey (ZFOURGE). Observed H-band sizes are measured from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) Hubble Space Telescope (HST)/F160W imaging. Distributions of individual galaxy masses and sizes illustrate that a clear mass-size relation exists up to z~7. At z~7, we find that the average galaxy size from the mass-size relation is more compact at a fixed mass of log(M*/Msol)=10.1, with r_1/2,maj=1.02+/-0.29 kpc, than at lower redshifts. This is consistent with our results from stacking the same CANDELS HST/F160W imaging, when we correct for galaxy position angle alignment. We find that the size evolution of star-forming galaxies is well fit by a power law of the form r_e = 7.07(1 + z)^-0.89 kpc, which is consistent with previous works for normal star-formers at 1<z<4. In order to compare our slope with those derived Lyman break galaxy studies, we correct for different IMFs and methodology and find a slope of -0.97+/-0.02, which is shallower than that reported for the evolution of Lyman break galaxies at z>4 (r_e\propto(1 +z)^-1.2+/-0.06). Therefore, we conclude the Lyman break galaxies likely represent a subset of highly star-forming galaxies that exhibit rapid size growth at z>4.

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.

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.

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

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

In order to reproduce the high-mass end of the galaxy mass-distribution, some process must be responsible for the suppression of star-formation in the most massive of galaxies. Commonly Active Galactic Nuclei (AGN) are invoked to fulfil this role, but the exact means by which they do so is still the topic of much debate, with studies finding evidence for both the suppression and enhancement of star-formation in AGN hosts. Using the ZFOURGE and NMBS galaxy surveys, we investigate the host galaxy properties of a mass-limited (M$_{\odot}$ $\ge$ 10$^{10.5}$ M$_{\odot}$), high-luminosity (L$_{1.4}$ $>$ 10$^{24}$ W Hz$^{-1}$) sample of radio-loud Active Galactic Nuclei to a redshift of z = 2.25. In contrast to low redshift studies, which associate radio-AGN activity with quiescent hosts, we find that the majority of z $>$ 1.5 radio-AGN are hosted by star-forming galaxies. Indeed, the stellar populations of radio-AGN are found to evolve with redshift in a manner that is consistent with the non-AGN mass-similar galaxy population. Interestingly, we find the radio-AGN fraction is constant across a redshift range of 0.25 $\le$ z $<$ 2.25, perhaps indicating that the radio-AGN duty cycle has little dependence on redshift or galaxy type. We do however see a strong relation between the radio-AGN fraction and stellar mass, with radio-AGN becoming rare below $\sim$ 10$^{10.5}$ M$_{\odot}$ or a halo-mass of 10$^{12}$ M$_{\odot}$. This halo-mass threshold is in good agreement with simulations that initiate radio-AGN feedback at this mass limit. Despite this we find that radio-AGN host star-formation rates are consistent with the non-AGN mass-similar galaxy sample, suggesting that while radio-AGN are in the right place to suppress star-formation in massive galaxies they are not necessarily responsible for doing so.

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 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 spectroscopically survey the galaxy cluster XMM-LSS J02182-05102 (hereafter IRC 0218) using LRIS (optical) and MOSFIRE (near-infrared) on Keck I as part of the ZFIRE survey. IRC 0218 has a narrow redshift range of $1.612<z_{\rm spec}<1.635$ defined by 33 members of which 20 are at R$_{\rm proj}<1$ Mpc. The cluster redshift and velocity dispersion are $z_{\rm cl}=1.6233\pm0.0003$ and $\sigma_{\rm cl}=254\pm50$ km s$^{-1}$. We reach NIR line sensitivities of $\sim0.3\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ that, combined with multi-wavelength photometry, provide extinction-corrected H$\alpha$ star formation rates (SFR), gas phase metallicities from [NII]/H$\alpha$, and stellar masses. We measure an integrated H$\alpha$ SFR of $\sim325{\rm M}_{\odot}$ yr$^{-1}$ (26 members; R$_{\rm proj}<2$ Mpc) and show that the elevated star formation in the cluster core (R$_{\rm proj}<0.25$ Mpc) is driven by the concentration of star-forming members, but the average SFR per H$\alpha$-detected galaxy is half that of members at R$_{\rm proj}\sim1$ Mpc. However, we do not detect any environmental imprint when comparing attenuation and gas phase metallicities: the cluster galaxies show similar trends with M$_{\star}$ as to the field, e.g. more massive galaxies have larger stellar attenuation. IRC 0218's gas phase metallicity-M$_{\star}$ relation (MZR) is offset to lower metallicities relative to $z\sim0$ and has a slope of $0.13\pm0.10$. Comparing the MZR in IRC 0218 to the COSMOS cluster at $z=2.1$ shows no evolution ($\Delta t\sim1$ Gyr): the MZR for both galaxy clusters are remarkably consistent with each other and virtually identical to several field surveys at $z\sim2$.

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.

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.

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 use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$).

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 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 study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$\sigma$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

We present an observational study of the stellar mass function of satellite galaxies around central galaxies at 0.2<z<1.2. Using statistical background subtraction of contaminating sources we derive satellite stellar mass distributions in four bins of central galaxy mass in three redshift ranges. Our results show that the stellar mass function of satellite galaxies increases with central galaxy mass, and that the distribution of satellite masses at fixed central mass is at most weakly dependent on redshift. We conclude that the average mass distribution of galaxies in groups is remarkably universal even out to z=1.2 and that it can be uniquely characterized by the group central galaxy mass. This further suggests that as central galaxies grow in stellar mass, they do so in tandem with the mass growth of their satellites. Finally, we classify all galaxies as either star forming or quiescent, and derive the mass functions of each subpopulation separately. We find that the mass distribution of both star forming and quiescent satellites show minimal redshift dependence at fixed central mass. However, while the fraction of quiescent satellite galaxies increases rapidly with increasing central galaxy mass, that of star forming satellites decreases.

Our understanding of the redshift $z>3$ galaxy population relies largely on samples selected using the popular "dropout" technique, typically consisting of UV-bright galaxies with blue colors and prominent Lyman breaks. As it is currently unknown if these galaxies are representative of the massive galaxy population, we here use the FourStar Galaxy Evolution (ZFOURGE) Survey to create a stellar mass-limited sample at $z=3-4$. Uniquely, ZFOURGE uses deep near-infrared medium-bandwidth filters to derive accurate photometric redshifts and stellar population properties. The mass-complete sample consists of 57 galaxies with log M $>10.6$, reaching below $M^{\star}$ at $z=3-4$. On average, the massive $z=3-4$ galaxies are extremely faint in the observed optical with median $R_{tot}^{AB}=27.48\pm0.41$ (restframe $M_{1700}=-18.05\pm0.37$). They lie far below the UV luminosity-stellar mass relation for Lyman break galaxies and are about $\sim100\times$ fainter at the same mass. The massive galaxies are red ($R-Ks_{AB}=3.9\pm0.2$; restframe UV-slope $\beta=-0.2\pm0.3$) likely from dust or old stellar ages. We classify the galaxy SEDs by their restframe $U-V$ and $V-J$ colors and find a diverse population: $46^{+6+10}_{-6-17}$% of the massive galaxies are quiescent, $54^{+8+17}_{-8-10}$% are dusty star-forming galaxies, and only $14^{+3+10}_{-3-4}$% resemble luminous blue star forming Lyman break galaxies. This study clearly demonstrates an inherent diversity among massive galaxies at higher redshift than previously known. Furthermore,we uncover a reservoir of dusty star-forming galaxies with $4\times$ lower specific star-formation rates compared to submillimeter-selected starbursts at $z>3$. With $5\times$ higher numbers, the dusty galaxies may represent a more typical mode of star formation compared to submillimeter-bright starbursts.

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 describe the Carnegie-Spitzer-IMACS (CSI) Survey, a wide-field, near-IR selected spectrophotometric redshift survey with IMACS on Magellan-Baade. CSI uses a flux-limited sample of galaxies in Spitzer IRAC 3.6micron imaging of SWIRE fields to efficiently trace the stellar mass of average galaxies to z~1.5. This paper provides an overview of the survey selection, observations, and processing of the photometry and spectrophotometry. We also describe the analysis of the data: new methods of fitting synthetic SEDs are used to derive redshifts, stellar masses, emission line luminosities, and coarse information on recent star-formation. Our unique methodology for analyzing low-dispersion spectra taken with multilayer prisms in IMACS, combined with panchromatic photometry from the ultraviolet to the IR, has yielded high quality redshifts for 43,347 galaxies in our first 5.3 sq. degs of the SWIRE XMM-LSS field. A new approach to assessing data quality is also described, and three different approaches are used to estimate our redshift errors, with robust agreement. Over the full range of 3.6micron fluxes of our selection, we find typical redshift uncertainties of sigma_z/(1+z) < 0.015. In comparisons with previously published spectroscopic redshifts we find scatters of sigma_z/(1+z) = 0.011 for galaxies at 0.7< z< 0.9, and sigma_z/(1+z) = 0.014 for galaxies at 0.9< z< 1.2. For galaxies brighter and fainter than i=23 mag, we find sigma_z/(1+z) = 0.008 and sigma_z/(1+z) = 0.022, respectively. Notably, our low-dispersion spectroscopy and analysis yields comparable redshift uncertainties and success rates for both red and blue galaxies, largely eliminating color-based systematics that can seriously bias observed dependencies of galaxy evolution on environment.

Using observations from the FourStar Galaxy Evolution Survey (ZFOURGE), we obtain the deepest measurements to date of the galaxy stellar mass function at 0.5 < z < 2.5. ZFOURGE provides well-constrained photometric redshifts made possible through deep medium-bandwidth imaging at 1-2um . We combine this with HST imaging from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS), allowing for the efficient selection of both blue and red galaxies down to stellar masses ~10^9.5 Msol at z ~ 2.5. The total surveyed area is 316 arcmin^2 distributed over three independent fields. We supplement these data with the wider and shallower NEWFIRM Medium-Band Survey (NMBS) to provide stronger constraints at high masses. Several studies at z<=1 have revealed a steepening of the slope at the low-mass end of the stellar mass function (SMF), leading to an upturn at masses <10^10 Msol that is not well-described by a standard single-Schechter function. We find evidence that this feature extends to at least z ~ 2, and that it can be found in both the star-forming and quiescent populations individually. The characteristic mass (M*) and slope at the lowest masses (alpha) of a double-Schechter function fit to the SMF stay roughly constant at Log(M/Msol) ~ 10.65 and ~-1.5 respectively. The SMF of star-forming galaxies has evolved primarily in normalization, while the change in shape is relatively minor. This is not the case for quiescent galaxies: the depth of our imaging allows us to show for the first time significantly more evolution at Log(M/Msol) < 10.5 than at higher masses. We find that the total mass density (down to 10^9 Msol) in star-forming galaxies has increased by a factor of ~2.2 since z ~ 2.5, whereas in quiescent galaxies it has increased by a factor of ~12 .

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 present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin^2) K-selected sources (K_AB < 24.0), split according to their rest-frame U - V vs. V - J colors into 72,216 star-forming and 9,034 quiescent galaxies, on maps from Spitzer/MIPS (24um), Herschel/SPIRE (250, 350, 500um), Herschel/PACS (100, 160um), and AzTEC (1100um). The fraction of the CIB resolved by our catalog is (69 $\pm$ 15)% at 24um, (78 $\pm$ 17)% at 70um, (58 $\pm$ 13)% at 100um, (78 $\pm$ 18)% at 160um, (80 $\pm$ 17)% at 250um, (69 $\pm$ 14)% at 350um, (65 $\pm$ 12)% at 500um, and (45 $\pm$ 8)% at 1100um. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at $\lambda$ < 200um is sourced predominantly from galaxies at z < 1, while at $\lambda$ > 200um the bulk originates from 1 < z < 2. Galaxies with stellar masses log(M/ M_sun)=9.5-11 are responsible for the majority of the CIB, with those in the log(M/ M_sun)=9.5-10 contributing mostly at $\lambda$ < 250um, and those in the log(M/ M_sun)=10.5-11 bin dominating at $\lambda$ > 350um. The contribution from galaxies in the log(M/ M_sun)=9.0-9.5 and log(M/ M_sun)=11.0-12.0 stellar mass bins contribute the least, both of order 5%, although the highest stellar-mass bin is a significant contributor to the luminosity density at z > 2. The luminosities of the galaxies responsible for the CIB shifts from a combination of "normal" and luminous infrared galaxies (LIRGs) at $\lambda$ < 160um, to LIRGs at 160um < $\lambda$ < 500um, to finally LIRGs and ultra-luminous infrared galaxies (ULIRGs) at $\lambda$ > 500um. Stacking analyses were performed with SIMSTACK (available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html) which accounts for possible biases due to clustering.

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.

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

Star-forming galaxies at redshifts z>6 are likely responsible for the reionization of the universe, and it is important to study the nature of these galaxies. We present three candidates for z~7 Lyman-break galaxies (LBGs) from a 155 arcmin^2 area in the CANDELS/COSMOS field imaged by the deep FourStar Galaxy Evolution (zFourGE) survey. The FourStar medium-band filters provide the equivalent of R~10 spectroscopy, which cleanly distinguishes between z~7 LBGs and brown dwarf stars. The distinction between stars and galaxies based on an object's angular size can become unreliable even when using HST imaging; there exists at least one very compact z~7 candidate (FWHM~0.5-1 kpc) that is indistinguishable from a point source. The medium-band filters provide narrower redshift distributions compared with broad-band-derived redshifts. The UV luminosity function derived using the three z~7 candidates is consistent with previous studies, suggesting an evolution at the bright end (MUV -21.6 mag) from z~7 to z~5. Fitting the galaxies' spectral energy distributions, we predict Lyman-alpha equivalent widths for the two brightest LBGs, and find that the presence of a Lyman-alpha line affects the medium-band flux thereby changing the constraints on stellar masses and UV spectral slopes. This illustrates the limitations of deriving LBG properties using only broad-band photometry. The derived specific star-formation rates for the bright LBGs are ~13 per Gyr, slightly higher than the lower-luminosity LBGs, implying that the star-formation rate increases with stellar mass for these galaxies.

We study the structural evolution of massive galaxies by linking progenitors and descendants at a constant cumulative number density of n_c=1.4x10^{-4} Mpc^{-3} to z~3. Structural parameters were measured by fitting Sersic profiles to high resolution CANDELS HST WFC3 J_{125} and H_{160} imaging in the UKIDSS-UDS at 1<z<3 and ACS I_{814} imaging in COSMOS at 0.25<z<1. At a given redshift, we selected the HST band that most closely samples a common rest-frame wavelength so as to minimize systematics from color gradients in galaxies. At fixed n_c, galaxies grow in stellar mass by a factor of ~3 from z~3 to z~0. The size evolution is complex: galaxies appear roughly constant in size from z~3 to z~2 and then grow rapidly to lower redshifts. The evolution in the surface mass density profiles indicates that most of the mass at r<2 kpc was in place by z~2, and that most of the new mass growth occurred at larger radii. This inside-out mass growth is therefore responsible for the larger sizes and higher Sersic indices of the descendants toward low redshift. At z<2, the effective radius evolves with the stellar mass as r_e M^{2.0}, consistent with scenarios that find dissipationless minor mergers to be a key driver of size evolution. The progenitors at z~3 were likely star-forming disks with r_e~2 kpc, based on their low Sersic index of n~1, low median axis ratio of b/a~0.52, and typical location in the star-forming region of the U-V versus V-J diagram. By z~1.5, many of these star-forming disks disappeared, giving rise to compact quiescent galaxies. Toward lower redshifts, these galaxies continued to assemble mass at larger radii and became the local ellipticals that dominate the high mass end of the mass function at the present epoch.

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

We describe the Carnegie-Spitzer-IMACS (CSI) Survey, a wide-field, near-IR selected spectrophotometric redshift survey with the Inamori Magellan Areal Camera and Spectrograph (IMACS) on Magellan-Baade. By defining a flux-limited sample of galaxies in Spitzer 3.6micron imaging of SWIRE fields, the CSI Survey efficiently traces the stellar mass of average galaxies to z~1.5. This first paper provides an overview of the survey selection, observations, processing of the photometry and spectrophotometry. We also describe the processing of the data: new methods of fitting synthetic templates of spectral energy distributions are used to derive redshifts, stellar masses, emission line luminosities, and coarse information on recent star-formation. Our unique methodology for analyzing low-dispersion spectra taken with multilayer prisms in IMACS, combined with panchromatic photometry from the ultraviolet to the IR, has yielded 37,000 high quality redshifts in our first 5.3 sq.degs of the SWIRE XMM-LSS field. We use three different approaches to estimate our redshift errors and find robust agreement. Over the full range of 3.6micron fluxes of our selection, we find typical uncertainties of sigma_z/(1+z) < 0.015. In comparisons with previously published VVDS redshifts, for example, we find a scatter of sigma_z/(1+z) = 0.012 for galaxies at 0.8< z< 1.2. For galaxies brighter and fainter than i=23 mag, we find sigma_z/(1+z) = 0.009 and sigma_z/(1+z) = 0.025, respectively. Notably, our low-dispersion spectroscopy and analysis yields comparable redshift uncertainties and success rates for both red and blue galaxies, largely eliminating color-based systematics that can seriously bias observed dependencies of galaxy evolution on environment.

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.

Using mass-selected galaxy samples from deep multiwavelength data we investigate the incidence of close galaxy pairs between z=0.4-2. Many such close pairs will eventually merge, and the pair fraction is therefore related to the merger rate. Over this redshift range the mean pair fraction is essentially constant (evolving as f_pair (1+z)^{-0.4 +/- 0.6}) with about 6+/-1% of massive galaxies having a 1:4 or greater companion within 30h^-1 kpc. Assuming the timescale over which pairs merge is not a strong function of redshift, this implies a similarly constant merger rate (per unit time) out to z=2. Since about three times as much cosmic time passes at z<1 as between z=1-2, this implies that correspondingly more mergers occur in the low-redshift universe. When minor companions (1:10 mass ratio or greater) are included, the pair fraction increases to ~20% and still does not evolve strongly with redshift. We also use a rest-frame color criterion to select pairs containing only quiescent galaxies (major "dry merger" progenitors), and find them to be similarly rare and constant with 4-7% of massive quiescent galaxies exhibiting a nearby companion. Thus, even though other studies find major mergers to be relatively uncommon since z=1, our results suggest that few additional mergers occur in the 1<z<2 range and other mechanisms may be required to explain the mass and size growth of galaxies over this epoch.

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.

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

Recent work has shown that the star formation-density relation -- in which galaxies with low star formation rates are preferentially found in dense environments -- is still in place at z~1, but the situation becomes less clear at higher redshifts. We use mass-selected samples drawn from the UKIDSS Ultra-Deep Survey to show that galaxies with quenched star formation tend to reside in dense environments out to at least z~1.8. Over most of this redshift range we are able to demonstrate that this star formation-density relation holds even at fixed stellar mass. The environmental quenching of star formation appears to operate with similar efficiency on all galaxies regardless of stellar mass. Nevertheless, the environment plays a greater role in the build-up of the red sequence at lower masses, whereas other quenching processes dominate at higher masses. In addition to a statistical analysis of environmental densities, we investigate a cluster at z=1.6, and show that the central region has an elevated fraction of quiescent objects relative to the field. Although the uncertainties are large, the environmental quenching efficiency in this cluster is consistent with that of galaxy groups and clusters at z~0. In this work we rely on photometric redshifts, and describe some of the pitfalls that large redshift errors can present.

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.

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

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

Much of the science that is made possible by multiwavelength redshift surveys requires the use of photometric redshifts. But as these surveys become more ambitious, and as we seek to perform increasingly accurate measurements, it becomes crucial to take proper account of the photometric redshift uncertainties. Ideally the uncertainties can be directly measured using a comparison to spectroscopic redshifts, but this may yield misleading results since spectroscopic samples are frequently small and not representative of the parent photometric samples. We present a simple and powerful empirical method to constrain photometric redshift uncertainties in the absence of spectroscopic redshifts. Close pairs of galaxies on the sky have a significant probability of being physically associated, and therefore of lying at nearly the same redshift. The difference in photometric redshifts in close pairs is therefore a measure of the redshift uncertainty. Some observed close pairs will arise from chance projections along the line of sight, but it is straightforward to perform a statistical correction for this effect. We demonstrate the technique using both simulated data and actual observations, and discuss how its usefulness can be limited by the presence of systematic photometric redshift errors. Finally, we use this technique to show how photometric redshift accuracy can depend on galaxy type.

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.

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

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

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

We investigate the properties of quiescent and star-forming galaxy populations to z~2 with purely photometric data, employing a novel rest-frame color selection technique. From the UKIDSS Ultra-Deep Survey Data Release 1, with matched optical and mid-IR photometry taken from the Subaru XMM Deep Survey and Spitzer Wide-Area Infrared Extragalactic Survey respectively, we construct a K-selected galaxy catalog and calculate photometric redshifts. Excluding stars, objects with uncertain z_phot solutions, those that fall in bad or incomplete survey regions, and those for which reliable rest-frame colors could not be derived, 30108 galaxies with K<22.4 (AB) and z<2.5 remain. The galaxies in this sample are found to occupy two distinct populations in the rest-frame U-V vs. V-J color space: a clump of red, quiescent galaxies (analogous to the red sequence) and a track of star-forming galaxies extending from blue to red U-V colors. This bimodal behavior is seen up to z~2. Due to a combination of measurement errors and passive evolution, the color-color diagram is not suitable to distinguish the galaxy bimodality at z>2 for this sample, but we show that MIPS 24um data suggest that a significant population of quiescent galaxies exists even at these higher redshifts. At z=1-2, the most luminous objects in the sample are divided roughly equally between star-forming and quiescent galaxies, while at lower redshifts most of the brightest galaxies are quiescent. Moreover, quiescent galaxies at these redshifts are clustered more strongly than those actively forming stars, indicating that galaxies with early-quenched star formation may occupy more massive host dark matter halos. This suggests that the end of star formation is associated with, and perhaps brought about by, a mechanism related to halo mass.

Recent studies have shown that distant red galaxies (DRGs), which dominate the high-mass end of the galaxy population at z~2.5, are more strongly clustered than the population of blue star-forming galaxies at similar redshifts. However these studies have been severely hampered by the small sizes of fields having deep near-infrared imaging. Here we use the large UKIDSS Ultra Deep Survey to study the clustering of DRGs. The size and depth of this survey allows for an unprecedented measurement of the angular clustering of DRGs at 2<z_phot<3 and K<21. The correlation function shows the expected power law behavior, but with an apparent upturn at theta<~10". We deproject the angular clustering to infer the spatial correlation length, finding 10.6+-1.6 h^-1 Mpc. We use the halo occupation distribution framework to demonstrate that the observed strong clustering of DRGs is not consistent with standard models of galaxy clustering, confirming previous suggestions that were based on smaller samples. Inaccurate photometric redshifts could artificially enhance the observed clustering, however significant systematic redshift errors would be required to bring the measurements into agreement with the models. Another possibility is that the underlying assumption that galaxies interact with their large-scale environment only through halo mass is not valid, and that other factors drive the evolution of the oldest, most massive galaxies at z~2.

We present K-band imaging of two ~30'x30' fields covered by the MUSYC Wide NIR Survey. The 1030 and 1255 fields were imaged with ISPI on the 4m Blanco telescope at CTIO to a 5 sigma point-source limiting depth of K~20 (Vega). Combining this data with the MUSYC Optical UBVRIz imaging, we created multi-band K-selected source catalogs for both fields. These catalogs, together with the MUSYC K-band catalog of the ECDF-S field, were used to select K<20 BzK galaxies over an area of 0.71 deg^2. This is the largest area ever surveyed for BzK galaxies. We present number counts, redshift distributions and stellar masses for our sample of 3261 BzK galaxies (2502 star-forming (sBzK) and 759 passively evolving (pBzK)), as well as reddening and star formation rate estimates for the star-forming BzK systems. We also present 2-point angular correlation functions and spatial correlation lengths for both sBzK and pBzK galaxies and show that previous estimates of the correlation function of these galaxies were affected by cosmic variance due to the small areas surveyed. We have measured correlation lengths r_0 of 8.89+/-2.03 Mpc and 10.82+/-1.72 Mpc for sBzK and pBzK galaxies respectively. This is the first reported measurement of the spatial correlation function of passive BzK galaxies. In the LambdaCDM scenario of galaxy formation, these correlation lengths at z~2 translate into minimum masses of ~4x10^{12} M_sun and ~9x10^{12} M_sun for the dark matter (DM) halos hosting sBzK and pBzK galaxies respectively. The clustering properties of the galaxies in our sample are consistent with them being the descendants of bright LBG at z~3, and the progenitors of present-day >1L* galaxies.

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

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

We present the first clustering results of X-ray selected AGN at z~3. Using Chandra X-ray imaging and UVR optical colors from MUSYC photometry in the ECDF-S field, we selected a sample of 58 z~3 AGN candidates. From the optical data we also selected 1385 LBG at 2.8<z< 3.8 with R<25.5. We performed auto-correlation and cross-correlation analyses, and here we present results for the clustering amplitudes and dark matter halo masses of each sample. For the LBG we find a correlation length of r_0,LBG = 6.7 +/- 0.5 Mpc, implying a bias value of 3.5 +/- 0.3 and dark matter (DM) halo masses of log(Mmin/Msun) = 11.8 +/- 0.1. The AGN-LBG cross-correlation yields r_0,AGN-LBG = 8.7 +/- 1.9 Mpc, implying for AGN at 2.8<z<3.8 a bias value of 5.5 +/- 2.0 and DM halo masses of log(Mmin/Msun) = 12.6 +0.5/-0.8. Evolution of dark matter halos in the Lambda CDM cosmology implies that today these z~3 AGN are found in high mass galaxies with a typical luminosity of 7+4/-2 L*.

We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host AGN. The LAEs exhibit a moderate spatial correlation length of r_0=3.6+0.8-1.0 Mpc, corresponding to a bias factor b=1.7+0.3-0.4, which implies median dark matter halo masses of log10(M_med) = 10.9+0.5-0.9 M_sun. Comparing the number density of LAEs, (1.5+-0.3)E-3/Mpc^3, with the number density of these halos finds a mean halo occupation ~1-10%. The evolution of galaxy bias with redshift implies that most z=3.1 LAEs evolve into present-day galaxies with L<2.5L*, whereas other z>3 galaxy populations typically evolve into more massive galaxies. Halo merger trees show that z=0 descendants occupy halos with a wide range of masses, with a median descendant mass close to that of L*. Only 30% of LAEs have sufficient stellar mass (>~3E9 M_sun) to yield detections in deep Spitzer-IRAC imaging. A two-population SED fit to the stacked UBVRIzJK+[3.6,4.5,5.6,8.0]micron fluxes of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0+0.6-0.4 E9 M_sun), moderate star formation rate (2+-1 M_sun/yr), a young component age of 20+30-10 Myr, and little dust (A_V<0.2). The best fit model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed.

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

The properties of galaxies are strongly correlated with their environment, with red galaxies dominating galaxy clusters and blue galaxies dominating the general field. However, not all field galaxies are young: studies of the colors, line strengths, and M/L ratios of massive early-type galaxies at 0<z<1.3 show that the most massive galaxies do not seem to care about their surroundings, and have very similar ages irrespective of their environment. There is good evidence that the growth of these galaxies does continue longer in the field than in clusters, via (nearly) dissipationless mergers of already old galaxies. These results are consistent with predictions of recent galaxy formation models, which incorporate AGN feedback to suppress star formation in the most massive halos. Systematic studies of the relation of galaxies with their environment beyond z=1 are difficult, and still somewhat contradictory. Intriguingly both the DEEP2 and VVDS surveys find that the color-density relation disappears at z~1.3, unfortunately just at the point where both surveys become highly incomplete. On the other hand, clustering studies at z~2.5 have shown that red galaxies cluster more strongly than blue galaxies, implying that the color-density relation was already in place at that redshift.

We present the rest-frame optical (B, V, and R band) luminosity functions (LFs) of galaxies at 2<z<3.5, measured from a K-selected sample constructed from the deep NIR MUSYC, the ultradeep FIRES, and the GOODS-CDFS. This sample is unique for its combination of area and range of luminosities. The faint-end slopes of the LFs at z>2 are consistent with those at z~0. The characteristic magnitudes are significantly brighter than the local values (e.g., ~1.2 mag in the R band), while the measured values for Phi_star are typically ~5 times smaller. The B-band luminosity density at z~2.3 is similar to the local value, and in the R band it is ~2 times smaller than the local value. We present the LF of Distant Red Galaxies (DRGs), which we compare to that of non-DRGs. While DRGs and non-DRGs are characterized by similar LFs at the bright end, the faint-end slope of the non-DRG LF is much steeper than that of DRGs. The contribution of DRGs to the global densities down to the faintest probed luminosities is 14%-25% in number and 22%-33% in luminosity. From the derived rest-frame U-V colors and stellar population synthesis models, we estimate the mass-to-light ratios (M/L) of the different subsamples. The M/L ratios of DRGs are ~5 times higher (in the R and V bands) than those of non-DRGs. The global stellar mass density at 2<z<3.5 appears to be dominated by DRGs, whose contribution is of order ~60%-80% of the global value. Qualitatively similar results are obtained when the population is split by rest-frame U-V color instead of observed J-K color.

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

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

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

We selected 40 candidate Lyman Alpha Emitting galaxies (LAEs) at z ~=3.1 with observed frame equivalent widths >150A and inferred emission line fluxes >2.5x10^-17 ergs/cm^2/s from deep narrow-band and broad-band MUSYC images of the Extended Chandra Deep Field South. Covering 992 sq. arcmin, this is the largest ``blank field'' surveyed for LAEs at z ~3, allowing an improved estimate of the space density of this population of 3+-1x10^-4 h_70^3/Mpc^3. Spectroscopic follow-up of 23 candidates yielded 18 redshifts, all at z ~=3.1. Over 80% of the LAEs are dimmer in continuum magnitude than the typical Lyman break galaxy spectroscopic limit of R= 25.5 (AB), with a median continuum magnitude R ~=27 and very blue continuum colors, (V-z) ~=0. Over 80% of the LAEs have the right UVR colors to be selected as Lyman break galaxies, but only 10% also have R<=25.5. Stacking the UBVRIzJK fluxes reveals that LAEs have stellar masses ~=5x10^8 h_70^-2 M_sun and minimal dust extinction, A_V < ~ 0.1. Inferred star formation rates are ~=6 h_70^-2 M_sun/yr, yielding a cosmic star formation rate density of 2x10^-3 h_70 M_sun/yr/Mpc^3. None of our LAE candidates show evidence for rest-frame emission line equivalent widths EW_rest>240A which might imply a non-standard IMF. One candidate is detected by Chandra, implying an AGN fraction of 2+-2% for LAE candidate samples. In summary, LAEs at z ~ 3 have rapid star formation, low stellar mass, little dust obscuration and no evidence for a substantial AGN component.

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

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

We present Spitzer 24micron imaging of 1.5 < z < 2.5 Distant Red Galaxies (DRGs) in the 10arcmin by 10arcmin Extended Hubble Deep Field South of the Multiwavelength Survey by Yale-Chile. We detect 65% of the DRGs with K_AB < 23.2 mag at S_24micron > 40uJy, and conclude that the bulk of the DRG population are dusty active galaxies. A mid-infrared (MIR) color analysis with IRAC data suggests that the MIR fluxes are not dominated by buried AGN, and we interpret the high detection rate as evidence for a high average star formation rate of <SFR> = 130+/-30 M/yr. From this, we infer that DRGs are important contributors to the cosmic star formation rate density at z ~ 2, at a level of \~0.02 M/yr/Mpc^3 to our completeness limit of K_AB = 22.9 mag

We present UBVRIz' optical images of the 0.32 square degree Extended Hubble Deep Field South. This is one of four fields comprising the MUSYC survey, which is optimized for the study of galaxies at z=3, AGN demographics, and Galactic structure. We calculate corrected aperture photometry and its uncertainties and find through tests that these provide a significant improvement upon standard techniques. Our photometric catalog of 62968 objects is complete to a total magnitude of R_AB=25. We select z=3 Lyman break galaxy (LBG) candidates from their UVR colors and find a sky surface density of 1.4/sq. arcmin and an angular correlation function w(theta) = 2.3+-1.0 theta^{-0.8}. (Abridged)

We study the observational signatures of the lensing signal produced by dark matter halos with embedded misaligned disks. This issue is of particular interest at the present time since most of the observed multiple lens systems have magnification ratios and image geometries that are not well-fit by standard mass models. The presence of substructure exterior to the lens has been invoked by several authors in the context of Cold Dark Matter (CDM) models in order to explain the anamolous magnification ratios. We emphasize that the anomalous magnification ratios may be an artifact of the simple one-component mass models currently in use; the inclusion of a misaligned disk may be able to mimic the effect of substructure. These slight spatial offsets between the dark matter halo and the disk, which are likely to occur during or as a consequence of interactions or mergers, lead to complex image configurations and non-standard magnification ratios. We investigate two illustrative lens models: a spiral disk embedded within a dark matter halo, and a compact disk-like component within an elliptical galaxy. The expected fraction of galaxies with a misaligned disk is estimated to be of the order of 10%. In such cases we find that the resultant lensing geometries are unusual, with high image multiplicities. The caustic structures - both radial and tangential - are drastically modified and the magnification ratios differ compared to expectations from standard lens models. The additional parameters required to specify the relative alignment of the mass components in the primary lens introduce yet another source of uncertainty in the mass modeling of gravitational lens systems.