We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0<z<1.4, drawn from the Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) survey. Based on more than 500 hours of Gemini/GMOS spectroscopy, and deep multi-band photometry taken with a range of observatories, we probe the SMFs down to a stellar mass limit of 10^9.7 Msun (10^9.5 Msun for star-forming galaxies). At this early epoch, the fraction of quiescent galaxies is already highly elevated in the clusters compared to the field at the same redshift. The quenched fraction excess (QFE) represents the fraction of galaxies that would be star-forming in the field, but are quenched due to their environment. The QFE is strongly mass dependent, and increases from ~30% at Mstar=10^9.7 Msun, to ~80% at Mstar=10^11.0 Msun. Nonetheless, the shapes of the SMFs of the two individual galaxy types, star-forming and quiescent galaxies, are identical between the clusters and the field - to high statistical precision. Yet, along with the different quiescent fractions is the total galaxy SMF environmentally dependent, with a relative deficit of low-mass galaxies in the clusters. These results are in stark contrast with findings in the local Universe, and thus require a substantially different quenching mode to operate at early times. We discuss these results in the light of several popular quenching models.

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields.

Minor mergers have been proposed as the driving mechanism for the size growth of quiescent galaxies with decreasing redshift. The process whereby large star-forming galaxies quench and join the quiescent population at the large size end has also been suggested as an explanation for this size growth. Given the clear association of quenching with clusters, we explore this mechanism by studying the structural properties of 23 spectroscopically identified recently quenched (or "poststarburst" (PSB)) cluster galaxies at $z\sim1$. Despite clear PSB spectral signatures implying rapid and violent quenching, 87\% of these galaxies have symmetric, undisturbed morphologies in the stellar continuum. Remarkably, they follow a mass-size relation lying midway between the star-forming and quiescent field relations, with sizes $0.1$ dex smaller than $z\sim1$ star-forming galaxies at log$(M_{*}/M_{\odot})=10.5$. This implies a rapid change in the light profile without directly effecting the stellar distribution, suggesting changes in the mass-to-light ratio gradients across the galaxy are responsible. We develop fading toy models to explore how star-forming galaxies move across the mass-size plane as their stellar populations fade to match those of the PSBs. "Outside-in" fading has the potential to reproduce the contraction in size and increase in bulge-dominance observed between star-forming and PSB cluster galaxies. Since cluster PSBs lie on the large size end of the quiescent mass-size relation, and our previous work shows cluster galaxies are smaller than field galaxies, the sizes of quiescent galaxies must grow both from the quenching of star-forming galaxies and dry minor mergers.

We present spectroscopic observations obtained at the {\it Large Binocular Telescope} in the field of the cluster XLSSJ0223-0436 at $z=1.22$. We confirm 12 spheroids cluster members and determine stellar velocity dispersion for 7 of them. We combine these data with those in the literature for clusters RXJ0848+4453 at $z=1.27$ (8 galaxies) and XMMJ2235-2557 at $z=1.39$ (7 galaxies) to determine the Fundamental Plane of cluster spheroids. We find that the FP at $z\sim1.3$ is offset and { rotated ($\sim3\sigma$)} with respect to the local FP. The offset corresponds to a mean evolution $\Delta$\rm{log}(M$_{dyn}$/L$_B$)=(-0.5$\pm$0.1)$z$. High-redshift galaxies follow a steeper mass-dependent M$_{dyn}$/L$_B$-M$_{dyn}$ relation than local ones. Assuming $\Delta$ log$(M_{dyn}/L_B)$=$\Delta$ log$(M^*/L_B)$, higher-mass galaxies (log(M$_{dyn}$/M$_\odot$)$\geq$11.5) have a higher-formation redshift ($z_f\geq$6.5) than lower-mass ones ($z_f\leq$2 for log(M$_{dyn}$/M$_\odot$$\leq$10)), with a median $z_f\simeq2.5$ for the whole sample. Also, galaxies with higher stellar mass density host stellar populations formed earlier than those in lower density galaxies. At fixed IMF, M$_{dyn}$/M$^*$ varies systematically with mass and mass density. It follows that the evolution of the stellar populations (M$^*/L_B$) accounts for the observed evolution of M$_{dyn}/L_B$ for M$_{dyn}$$>10^{11}$ M$_\odot$ galaxies, while accounts for $\sim$85\% of the evolution at M$_{dyn}$$<10^{11}$ M$_\odot$. We find no evidence in favour of structural evolution of individual galaxies, while we find evidences that spheroids later added to the population account for the observed discrepancy at masses $<10^{11}$ M$_\odot$. [Abridged]

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

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

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

We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 \mu$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we derive a shallow faint end slope of $ \alpha \sim -0.35^{+0.15}_{-0.15} $ and $ M_{H}^{*} \sim -23.52^{+0.15}_{-0.17} $, suggesting that there is a deficit of faint red sequence galaxies in clusters at high redshift. By comparing the stacked red sequence LF of our sample with a sample of clusters at z~0.6, we find an evolution in the faint end of the red sequence over the ~2.6 Gyr between the two samples, with the mean faint end red sequence luminosity growing by more than a factor of two. The faint-to-luminous ratio of our sample ($0.78^{+0.19}_{-0.15}$) is consistent with the trend of decreasing ratio with increasing redshift as proposed in previous studies. A comparison with the field shows that the faint-to-luminous ratios in clusters are consistent with the field at z~1.15 and exhibit a stronger redshift dependence. Our results support the picture that the build up of the faint red sequence galaxies occurs gradually over time and suggest that faint cluster galaxies, similar to bright cluster galaxies, experience the quenching effect induced by environment already at z~1.15.

The bulk of the star-formation rate density peak at cosmic noon was obscured by dust. How accurately we can assess the role of dust obscured star-formation is affected by inherent biases in our empirical methods -- both those that rely on direct dust emission and those that rely on the inferred dust attenuation of starlight. We use a library of hydrodynamic simulations with radiative transfer to explore these biases. We find that for IR luminous galaxies that are in rapidly quenching systems (e.g. post-coalescence) standard luminosity-to-SFR relations can strongly overestimate the true SFRs. We propose using the $L_{IR}/L_{1.6}$ color to both help identify such systems and provide more accurate SFRs. Conversely, we find that the diagnostic UVJ plot misidentifies a subset of dusty star-forming galaxies. This is due to variability in the effective attenuation curves including being much grayer in the optical-to-near-IR regime than the Calzetti starburst law. This is in agreement with recent observations of IR-selected galaxies at cosmic noon. Our results support the view that we need a panchromatic approach from the rest-frame UV through the IR and SED modeling that includes realistic SFHs and allows for variable attenuation curves if we want to fully account for dust obscured star-formation across the epochs of greatest galaxy build-up.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present the analysis, based on spectra collected at the Large Binocular Telescope, of the stellar populations in seven spheroidal galaxies in the cluster XLSSJ0223 at $z$$\sim$1.22. The aim is to constrain the epoch of their formation and their star formation history. Using absorption line strenghts and full spectral fitting, we derive for the stellar populations of the seven spheroids a median age <Age>=2.4$\pm$0.6 Gyr, corresponding to a median formation redshift $<$$z_f$$>$$\sim2.6_{-0.5}^{+0.7}$ (lookback time = 11$_{-1.0}^{+0.6}$ Gyr). We find a significant scatter in age, showing that massive spheroids, at least in our targeted cluster, are not coeval. The median metallicity is [Z/H]=0.09$\pm$0.16, as for early-types in clusters at 0$<$$z$<0.9. This lack of evolution of [Z/H] over the range 0$<$$z$$<$1.3, corresponding to the last 9 billions years, suggests that no significant additional star formation and chemical enrichment are required for cluster spheroids to reach the present-day population. We do not detect significant correlation between age and velocity dispersion $\sigma_e$, or dynamical mass M$_{dyn}$, or effective stellar mass density $\Sigma_e$. On the contrary, the metallicity [Z/H] of the seven spheroids is correlated to their dynamical mass M$_{dyn}$, according to a relation similar to the one for local spheroids. [Z/H] is also anticorrelated to stellar mass density $\Sigma_e$ because of the anticorrelation between M$_{dyn}$ and $\Sigma_e$. Therefore, the basic trends observed in the local universe were already established at $z\sim1.3$, i.e. more massive spheroids are more metal rich, have lower stellar mass density and tend to be older than lower-mass galaxies.

Minor mergers are thought to be responsible for the size growth of quiescent field galaxies with decreasing redshift. We test this hypothesis using the cluster environment as a laboratory. Satellite galaxies in clusters move at high velocities, making mergers between them rare. The stellar mass-size relation in ten clusters and in the field is measured and compared at $z~\mathtt{\sim}~1$. Our cluster sample contains 344 spectroscopically-confirmed cluster members with Gemini/GMOS and 182 confirmed with HST WFC3 G141 grism spectroscopy. On average, quiescent and star-forming cluster galaxies are smaller than their field counterparts by ($0.08\pm0.04$) dex and ($0.07\pm0.01$) dex respectively. These size offsets are consistent with the average sizes of quiescent and star-forming field galaxies between $1.2\leqslant z\leqslant1.5$, implying the cluster environment has inhibited size growth between this period and $z~\mathtt{\sim}~1$. The negligible differences measured between the $z~\mathtt{\sim}~0$ field and cluster quiescent mass-size relations in other works imply that the average size of quiescent cluster galaxies must rise with decreasing redshift. Using a toy model, we show that the disappearance of the compact cluster galaxies might be explained if, on average, $\mathtt{\sim}40\%$ of them merge with their brightest cluster galaxies (BCGs) and $\mathtt{\sim}60\%$ are tidally destroyed into the intra-cluster light (ICL) between $0\leqslant z\leqslant1$. This is in agreement with the observed stellar mass growth of BCGs between $0\leqslant z\leqslant1$ and the observed ICL stellar mass fraction at $z~\mathtt{\sim}~0$. Our results support minor mergers as the cause for the size growth in quiescent field galaxies, with cluster-specific processes responsible for the similarity between the field and cluster quiescent mass-size relations at low redshift.

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

We present a study of the influence of galaxy mergers on star formation at 0.3<z<2.5. Major mergers are selected from the CANDELS/3D-HST catalog using a peak-finding algorithm. Mergers have projected galaxy nuclei separation of their members between 3-15 kpc. We compare the star formation activity in merging and non-merging galaxies and find no significant differences. We find that only 12% of the galaxies in major mergers (in which both galaxies have log(M/Msun)>10) are star-bursting (i.e., with SFR above the main sequence of star-forming galaxies by >0.5 dex). Merging galaxies which include galaxies with lower masses show a higher fraction of star-bursting galaxies (20%). The low fraction of star-bursting merging galaxies in this sample suggests that at galaxy nuclei separations of 3-15 kpc merging galaxies are still in a early stage and are yet to reach the maximum level of star formation activity. Furthermore, the level of star formation enhancement and its duration could be arguably reduced compared to local mergers, as shown by simulations of high-z mergers, and might also depend on the physical properties (such as stellar mass and gas fraction) of the merging galaxies. Finally, we compare the specific SFR between merging galaxies. Our results suggest that, as the mass of the merging galaxies increases, the star formation activity in the less massive member in the merger suffers a more dramatic impact than its companion galaxy.

We present spatially resolved emission diagnostics for eight $z\sim0.9$ galaxies that demonstrate extended low ionization emission-line regions (LIERs) over kpc scales. Eight candidates are selected based on their spatial extent and emission line fluxes from slitless spectroscopic observations with the HST/WFC3 G141 and G800L grisms in the well-studied GOODS survey fields. Five of the candidates (62.5%) are matched to X-ray counterparts in the \textit{Chandra X-Ray Observatory} Deep Fields. We modify the traditional Baldwin-Philips-Terlevich (BPT) emission line diagnostic diagram to use [SII]/(H$\alpha$+[NII]) instead of [NII]/H$\alpha$ to overcome the blending of [NII] and H$\alpha$+[NII] in the low resolution slitless grism spectra. We construct emission line ratio maps and place the individual pixels in the modified BPT. The extended LINER-like emission present in all of our candidates, coupled with X-Ray properties consistent with star-forming galaxies and weak [OIII]$\lambda$5007\AA\ detections, is inconsistent with purely nuclear sources (LINERs) driven by active galactic nuclei. While recent ground-based integral field unit spectroscopic surveys have revealed significant evidence for diffuse LINER-like emission in galaxies within the local universe $(z\sim0.04)$, this work provides the first evidence for the non-AGN origin of LINER-like emission out to high redshifts.

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

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

We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1\sigma$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal.

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

Using observed stellar mass functions out to $z=5$, we measure the main progenitor stellar mass growth of descendant galaxies with masses of $\log{M_{*}/M_{\odot}}=11.5,11.0,10.5,10.0$ at $z\sim0.1$ using an evolving cumulative number density selection. From these mass growth histories, we are able to measure the time at which half the total stellar mass of the descendant galaxy was assembled, $t_{a}$, which, in order of decreasing mass corresponds to redshifts of $z_{a}=1.28, 0.92, 0.60$ and $0.51$. We compare this to the median light-weighted stellar age $t_{*}$ ($z_{*} = 2.08, 1.49, 0.82$ and $0.37$) of a sample of low redshift SDSS galaxies (from the literature) and find the timescales are consistent with more massive galaxies forming a higher fraction of their stars ex-situ compared to lower mass descendants. We find that both $t_{*}$ and $t_{a}$ strongly correlate with mass which is in contrast to what is found in the EAGLE hydrodynamical simulation which shows a flat relationship between $t_{a}$ and $M_{*}$. However, the semi-analytic model of \citet{henriques2015} is consistent with the observations in both $t_{a}$ and $t_{*}$ with $M_{*}$, showing the most recent semi-analytic models are better able to decouple the evolution of the baryons from the dark matter in lower-mass galaxies.

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

We present the analysis and results of a spectroscopic follow-up program of a mass-selected sample of six galaxies at 3 < z < 4 using data from Keck-NIRSPEC and VLT-Xshooter. We confirm the z > 3 redshifts for half of the sample through the detection of strong nebular emission lines, and improve the zphot accuracy for the remainder of the sample through the combination of photometry and spectra. The modeling of the emission-line-corrected spectral energy distributions (SEDs) adopting improved redshifts confirms the very large stellar masses of the sample (M_* ~ 1.5-4 x 10^11 Msun) in the first 2 Gyrs of cosmic history, with a diverse range in stellar ages, star formation rates and dust content. From the analysis of emission line luminosities and widths, and far-infrared (FIR) fluxes we confirm that >80% of the sample are hosts to luminous hidden active galactic nuclei (AGNs), with bolometric luminosities of ~10^(44-46) erg/s. We find that the MIPS 24um photometry is largely contaminated by AGN continuum, rendering the SFRs derived using only 24um photometry to be severely overestimated. By including the emission from the AGN in the modeling of the UV-to-FIR SEDs, we confirm that the presence of the AGN does not bias considerably the stellar masses (< 0.3 dex at 1\sigma). We show evidence for a rapid increase of the AGN fraction from ~30% to ~60-100% over the 1 Gyr between z~2 and z~3. Although we cannot exclude some enhancement of the AGN fraction for our sample due to selection effects, the small measured [OIII] contamination to the observed K-band fluxes suggests that our sample is not significantly biased toward massive galaxies hosting AGNs.

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

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

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

The James Webb Space Telescope's Medium Resolution Spectrometer (MRS), will offer nearly 2 orders of magnitude improvement in sensitivity and >3X improvement in spectral resolution over our previous space-based mid-IR spectrometer, the Spitzer IRS. In this paper, we make predictions for spectroscopic pointed observations and serendipitous detections with the MRS. Specifically, pointed observations of Herschel sources require only a few minutes on source integration for detections of several star-forming and active galactic nucleus lines, out to z$=$3 and beyond. But the same data will also include tens of serendipitous 0$\lesssim$z$\lesssim$4 galaxies per field with infrared luminosities ranging $\sim10^6-10^{13}$L$_{\odot}$. In particular, for the first time and for free we will be able to explore the $L_{IR}<10^{9}L_{\odot}$ regime out to $z\sim3$. We estimate that with $\sim$100 such fields, statistics of these detections will be sufficient to constrain the evolution of the low-$L$ end of the infrared luminosity function, and hence the star formation rate function. The above conclusions hold for a wide range in potential low-$L$ end of the IR luminosity function, and accounting for the PAH deficit in low-$L$, low-metallicity galaxies.

We compared the properties of 56 elliptical galaxies selected from three clusters at $1.2<z<1.4$ with those of field galaxies in the GOODS-S (~30), COSMOS (~180) and CANDELS (~220) fields. We studied the relationships among effective radius, surface brightness, stellar mass, stellar mass density $\Sigma_{Re}$ and central mass density $\Sigma_{1kpc}$ within 1 kpc radius. We find that cluster ellipticals do not differ from field ellipticals: they share the same structural parameters at fixed mass and the same scaling relations. On the other hand, the population of field ellipticals at $z\sim1.3$ shows a significant lack of massive ($M_*> 2\times 10^{11}$ M$_\odot$) and large (R$_e > 4-5$ kpc) ellipticals with respect to the cluster. Nonetheless, at $M*<2\times 10^{11}$ M$_\odot$, the two populations are similar. The size-mass relation of ellipticals at z~1.3 defines two different regimes, above and below a transition mass $m_t\sim 2-3\times10^{10}$ M$_\odot$: at lower masses the relation is nearly flat (R$_e\propto M_*^{-0.1\pm 0.2}$), the mean radius is constant at ~1 kpc and $\Sigma_{Re}\sim \Sigma_{1kpc}$ while, at larger masses, the relation is R$_e\propto M*^{0.64\pm0.09}$. The transition mass marks the mass at which galaxies reach the maximum $\Sigma_{Re}$. Also the $\Sigma_{1kpc}$-mass relation follows two different regimes, $\Sigma_{1kpc}\propto M*^{0.64\ >m_t}_{1.07\ <m_t}$, defining a transition mass density $\Sigma_{1kpc}\sim 2-3\times10^3$ M$_\odot$ pc$^{-2}$. The mass density $\Sigma_{Re}$ does not correlate with mass, dense/compact galaxies can be assembled over a wide mass regime, independently of the environment. The central mass density, $\Sigma_{1kpc}$, besides to be correlated with the mass, is correlated to the age of the stellar population: the higher the central stellar mass density, the higher the mass, the older the age of the stellar population. [Abridged]

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

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

Using a sample of 98 galaxy clusters recently imaged in the near infra-red with the ESO NTT, WIYN and WHT telescopes, supplemented with 33 clusters from the ESO archive, we measure how the stellar mass of the most massive galaxies in the universe, namely Brightest Cluster Galaxies (BCG), increases with time. Most of the BCGs in this new sample lie in the redshift range $0.2<z<0.6$, which has been noted in recent works to mark an epoch over which the growth in the stellar mass of BCGs stalls. From this sample of 132 clusters, we create a subsample of 102 systems that includes only those clusters that have estimates of the cluster mass. We combine the BCGs in this subsample with BCGs from the literature, and find that the growth in stellar mass of BCGs from 10 billion years ago to the present epoch is broadly consistent with recent semi-analytic and semi-empirical models. As in other recent studies, tentative evidence indicates that the stellar mass growth rate of BCGs may be slowing in the past 3.5 billion years. Further work in collecting larger samples, and in better comparing observations with theory using mock images is required if a more detailed comparison between the models and the data is to be made.

The growth of galaxies is a key problem in understanding the structure and evolution of the universe. Galaxies grow their stellar mass by a combination of star formation and mergers, with a relative importance that is redshift dependent. Theoretical models predict quantitatively different contributions from the two channels; measuring these from the data is a crucial constraint. Exploiting the UltraVISTA catalog and a unique sample of progenitors of local ultra massive galaxies selected with an abundance matching approach, we quantify the role of the two mechanisms from z=2 to 0. We also compare our results to two independent incarnations of semi-analytic models. At all redshifts, progenitors are found in a variety of environments, ranging from being isolated to having 5-10 companions with mass ratio at least 1:10 within a projected radius of 500 kpc. In models, progenitors have a systematically larger number of companions, entailing a larger mass growth for mergers than in observations, at all redshifts. Generally, in both observations and models, the inferred and the expected mass growth roughly agree, within the uncertainties. Overall, our analysis confirms the model predictions, showing how the growth history of massive galaxies is dominated by in situ star formation at z~2, both star-formation and mergers at 1<z<2, and by mergers alone at z<1. Nonetheless, detailed comparisons still point out to tensions between the expected mass growth and our results, which might be due to either an incorrect progenitors-descendants selection, uncertainties on star formation rate and mass estimates, or the adopted assumptions on merger rates.

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

We present the first spectroscopic confirmation of an ultra-massive galaxy at redshift z>3 using data from Keck-NIRSPEC, VLT-Xshooter, and GTC-Osiris. We detect strong [OIII] and Ly$\alpha$ emission, and weak [OII], CIV, and HeII, placing C1-23152 at a spectroscopic redshift of $z_{spec}$=3.351. The modeling of the emission-line corrected spectral energy distribution results in a best-fit stellar mass of $M_{*}=3.1^{+0.6}_{-0.7}\times10^{11} M_{\odot}$, a star-formation rate of <7 $M_{\odot} yr^{-1}$, and negligible dust extinction. The stars appear to have formed in a short intense burst ~300-500 Myr prior to the observation epoch, setting the formation redshift of this galaxy at z~4.1. From the analysis of the line ratios and widths, and the observed flux at 24$\mu$m, we confirm the presence of a luminous hidden active galactic nucleus (AGN), with bolometric luminosity of ~$10^{46}erg s^{-1}$. Potential contamination to the observed SED from the AGN continuum is constrained, placing a lower limit on the stellar mass of $2\times10^{11} M_{\odot}$. HST/WFC3 $H_{160}$ and ACS $I_{814}$ images are modeled, resulting in an effective radius of $r_{e}$~1 kpc in the $H_{160}$ band and a Sersic index n~4.4. This object may be a prototype of the progenitors of local most massive elliptical galaxies in the first 2 Gyr of cosmic history, having formed most of its stars at z>4 in a highly dissipative, intense, and short burst of star formation. C1-23152 is completing its transition to a post-starburst phase while hosting a powerful AGN, potentially responsible for the quenching of the star formation activity.

We build a Spitzer IRAC complete catalog of objects, obtained by complementing the $K_\mathrm{s}$-band selected UltraVISTA catalog with objects detected in IRAC only. With the aim of identifying massive (i.e., $\log(M_*/M_\odot)>11$) galaxies at $4<z<7$, we consider the systematic effects on the measured photometric redshifts from the introduction of an old and dusty SED template and from the introduction of a bayesian prior taking into account the brightness of the objects, as well as the systematic effects from different star formation histories (SFHs) and from nebular emission lines in the recovery of stellar population parameters. We show that our results are most affected by the bayesian luminosity prior, while nebular emission lines and SFHs only introduce a small dispersion in the measurements. Specifically, the number of $4<z<7$ galaxies ranges from 52 to 382 depending on the adopted configuration. Using these results we investigate, for the first time, the evolution of the massive end of the stellar mass functions (SMFs) at $4<z<7$. Given the rarity of very massive galaxies in the early universe, major contributions to the total error budget come from cosmic variance and poisson noise. The SMF obtained without the introduction of the bayesian luminosity prior does not show any evolution from $z\sim6.5$ to $z\sim 3.5$, implying that massive galaxies could already be present when the Universe was $\sim0.9$~Gyr old. However, the introduction of the bayesian luminosity prior reduces the number of $z>4$ galaxies with best fit masses $\log(M_*/M_\odot)>11$ by 83%, implying a rapid growth of very massive galaxies in the first 1.5 Gyr of cosmic history. From the stellar-mass complete sample, we identify one candidate of a very massive ($\log(M_*/M_\odot)\sim11.5$), quiescent galaxy at $z\sim5.4$, with MIPS $24\mu$m detection suggesting the presence of a powerful obscured AGN.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

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.

Using the UltraVISTA catalogs, we investigate the evolution in the 11.4~Gyr since $z=3$ of the progenitors of local ultra-massive galaxies ($\log{(M_{\rm star}/M_{\odot})}\approx11.8$; UMGs), providing a complete and consistent picture of how the most massive galaxies at $z=0$ have assembled. By selecting the progenitors with a semi-empirical approach using abundance matching, we infer a growth in stellar mass of 0.56$^{+0.35}_{-0.25}$ dex, 0.45$^{+0.16}_{-0.20}$~dex, and 0.27$^{+0.08}_{-0.12}$ dex from $z=3$, $z=2$, and $z=1$, respectively, to $z=0$. At $z<1$, the progenitors of UMGs constitute a homogeneous population of only quiescent galaxies with old stellar populations. At $z>1$, the contribution from star-forming galaxies progressively increases, with the progenitors at $2<z<3$ being dominated by massive ($M_{\rm star} \approx 2 \times 10^{11}$M$_{\odot}$), dusty ($A_{\rm V}\sim$1--2.2 mag), star-forming (SFR$\sim$100--400~M$_{\odot}$ yr$^{-1}$) galaxies with a large range in stellar ages. At $z=2.75$, $\sim$15\% of the progenitors are quiescent, with properties typical of post-starburst galaxies with little dust extinction and strong Balmer break, and showing a large scatter in color. Our findings indicate that at least half of the stellar content of local UMGs was assembled at $z>1$, whereas the remaining was assembled via merging from $z\sim 1$ to the present. Most of the quenching of the star-forming progenitors happened between $z=2.75$ and $z=1.25$, in good agreement with the typical formation redshift and scatter in age of $z=0$ UMGs as derived from their fossil records. The progenitors of local UMGs, including the star-forming ones, never lived on the blue cloud since $z=3$. We propose an alternative path for the formation of local UMGs that refines previously proposed pictures and that is fully consistent with our findings.

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

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

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

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.

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.

Comparing galaxies across redshifts at fixed cumulative number density is a popular way to estimate the evolution of specific galaxy populations. This method ignores scatter in mass accretion histories and galaxy-galaxy mergers, which can lead to errors when comparing galaxies over large redshift ranges (Delta z > 1). We use abundance matching in the LCDM paradigm to estimate the median change in number density with redshift and provide a simple fit (+0.16 dex per unit Delta z) for progenitors of z = 0 galaxies. We find that galaxy descendants do not evolve in the same way as galaxy progenitors, largely due to scatter in mass accretion histories. We also provide estimates for the 1-sigma range of number densities corresponding to galaxy progenitors and descendants. Finally, we discuss some limits on number density comparisons, which arise due to difficulties measuring physical quantities (e.g., stellar mass) consistently across redshifts. A public tool to calculate number density evolution for galaxies, as well as approximate halo masses, is available online.

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

Quiescent massive galaxies at z~2 are thought to be the progenitors of present-day massive ellipticals. Observations revealed them to be extraordinarily compact. The determination of stellar ages, star formation rates and dust properties via spectroscopic measurements has up to now only been feasible for the most luminous and massive specimens (~3x M*). Here we present a spectroscopic study of two near-infrared selected galaxies which are close to the characteristic stellar mass M* (~0.9x M* and ~1.3x M*) and whose observed brightness has been boosted by the gravitational lensing effect. We measure the redshifts of the two galaxies to be z=1.71\pm0.02 and z=2.15\pm0.01. By fitting stellar population synthesis models to their spectro-photometric SEDs we determine their ages to be 2.4^{+0.8}_{-0.6} Gyr and 1.7\pm0.3 Gyr, respectively, which implies that the two galaxies have higher mass-to-light ratios than most quiescent z~2 galaxies in other studies. We find no direct evidence for active star-formation or AGN activity in either of the two galaxies, based on the non-detection of emission lines. Based on the derived redshifts and stellar ages we estimate the formation redshifts to be z=4.3^{+3.4}_{-1.2} and z=4.3^{+1.0}_{-0.6}, respectively. We use the increased spatial resolution due to the gravitational lensing to derive constraints on the morphology. Fitting Sersic profiles to the de-lensed images of the two galaxies confirms their compactness, with one of them being spheroid-like, and the other providing the first confirmation of a passive lenticular galaxy at a spectroscopically derived redshift z~2.

We present measurements of the stellar mass functions (SMFs) of star-forming and quiescent galaxies to z = 4 using a sample of 95 675 galaxies in the COSMOS/UltraVISTA field. Sources have been selected from the DR1 UltraVISTA K_{s}-band imaging which covers a unique combination of a wide area (1.62 deg^2), to a significant depth (K_{s,tot} = 23.4). The SMFs of the combined population are in good agreement with previous measurements and show that the stellar mass density of the universe was only 50%, 10% and 1% of its current value at z ~ 0.75, 2.0, and 3.5, respectively. The quiescent population drives most of the overall growth, with the stellar mass density of these galaxies increasing by 2.71^{+0.93}_{-0.22} dex since z = 3.5. At z > 2.5, star-forming galaxies dominate the total SMF at all stellar masses, although a nonzero population of quiescent galaxies persists to z = 4. Comparisons of the K_{s}-selected star-forming galaxy SMFs to UV-selected SMFs at 2.5 < z < 4 show reasonable agreement and suggests UV-selected samples are representative of the majority of the stellar mass density at z > 3.5. We estimate the average mass growth of individual galaxies by selecting galaxies at fixed cumulative number density. The average galaxy with Log(M_{*}/M_{sun}) = 11.5 at z = 0.3 has grown in mass by only 0.2 dex (0.3 dex) since z = 2.0(3.5), whereas those with Log(M_{*}/M_{sun}) = 10.5 have grown by > 1.0 dex since z = 2. At z < 2, the time derivatives of the mass growth are always larger for lower-mass galaxies, which demonstrates that the mass growth in galaxies since that redshift is mass-dependent and primarily bottom-up. Lastly, we examine potential sources of systematic uncertainties on the SMFs and find that those from photo-z templates, SPS modeling, and the definition of quiescent galaxies dominate the total error budget in the SMFs.

We present improved constraints on the formation time scale of massive cluster galaxies based on rest-frame optical spectra of galaxies in a forming cluster located at z=2.16. The spectra are obtained with MOIRCS on the Subaru telescope with an integration time of ~7 hours. We achieve accurate redshift measurements by fitting SEDs using the spectra and broad-band photometry simultaneously, allowing us to identify probable cluster members. Clusters at low redshifts are dominated by quiescent galaxies, but we find that quiescent galaxies and star forming galaxies co-exist in this z=2 system. Interestingly, the quiescent galaxies form a weak red sequence in the process of forming. By stacking the spectra of star forming galaxies, we observe strong emission lines such as [OII] and [OIII] and we obtain a tentative hint of AGN activities in these galaxies. On the other hand, the stacked spectrum of the quiescent galaxies reveals a clear 4000A break with a possible CaII H+K absorption feature and strong emission lines such as [OII] are absent in the spectrum, confirming the quiescent nature of these galaxies. We then perform detailed spectral analyses of the stacked spectrum, which suggest that these massive quiescent galaxies formed at redshifts between 3 and 4 on a time scale of <~0.5 Gyr. This short formation time scale is not reproduced in recent numerical simulations. We discuss possible mechanisms for how these galaxies form 10^11 Msun stellar mass on a short time scale and become red and quiescent by z=2.

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

We present the stellar mass functions (SMFs) of star-forming and quiescent galaxies from observations of 10 rich clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS) in the redshift range 0.86<z<1.34. We compare our results with field measurements at similar redshifts using data from a Ks-band selected catalogue of the COSMOS/UltraVISTA field. We construct a Ks-band selected multi-colour catalogue for the clusters in 11 photometric bands covering u-8um, and estimate photometric redshifts and stellar masses using SED fitting techniques. To correct for interlopers in our cluster sample, we use the deep spectroscopic component of GCLASS, which contains spectra for 1282 identified cluster and field galaxies taken with Gemini/GMOS. Both the photometric and spectroscopic samples are sufficiently deep that we can probe the SMF down to masses of 10^10 Msun. We distinguish between star-forming and quiescent galaxies using the rest-frame U-V versus V-J diagram, and find that the best-fitting Schechter parameters alpha and M* are similar within the uncertainties for these galaxy types within the different environments. However, there is a significant difference in the shape and normalisation of the total SMF between the clusters and the field sample. This difference in the total SMF is primarily a reflection of the increased fraction of quiescent galaxies in high-density environments. We apply a simple quenching model that includes components of mass- and environment-driven quenching, and find that in this picture 45% of the star-forming galaxies, which normally would be forming stars in the field, are quenched by the cluster. If galaxies in clusters and the field quench their star formation via different mechanisms, these processes have to conspire in such a way that the shapes of the quiescent and star-forming SMF remain similar in these different environments.

We present a catalog covering 1.62 deg^2 of the COSMOS/UltraVISTA field with PSF-matched photometry in 30 photometric bands. The catalog covers the wavelength range 0.15um - 24um including the available GALEX, Subaru, CFHT, VISTA and Spitzer data. Catalog sources have been selected from the DR1 UltraVISTA Ks band imaging that reaches a depth of K_{s,tot} = 23.4 AB (90% completeness). The PSF-matched catalog is generated using position-dependent PSFs ensuring accurate colors across the entire field. Also included is a catalog of photometric redshifts (z_phot) for all galaxies computed with the EAZY code. Comparison with spectroscopy from the zCOSMOS 10k bright sample shows that up to z ~ 1.5 the z_phot are accurate to dz/(1 + z) = 0.013, with a catastrophic outlier fraction of only 1.6%. The z_phot also show good agreement with the z_phot from the NEWFIRM Medium Band Survey (NMBS) out to z ~ 3. A catalog of stellar masses and stellar population parameters for galaxies determined using the FAST spectral energy distribution fitting code is provided for all galaxies. Also included are rest-frame U-V and V-J colors, L_2800 and L_IR. The UVJ color-color diagram confirms that the galaxy bi-modality is well-established out to z ~ 2. Star-forming galaxies also obey a star forming "main sequence" out to z ~ 2.5, and this sequence evolves in a manner consistent with previous measurements. The COSMOS/UltraVISTA Ks-selected catalog covers a unique parameter space in both depth, area, and multi-wavelength coverage and promises to be a useful tool for studying the growth of the galaxy population out to z ~ 3 - 4.

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.

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

We present the rest-frame J- and H-band luminosity function (LF) of field galaxies, based on a deep multi-wavelength composite sample from the MUSYC, FIRES and FIREWORKS survey public catalogues, covering a total area of 450 arcmin^2. The availability of flux measurements in the Spitzer IRAC 3.6, 4.5, 5.8, and 8 um channels allows us to compute absolute magnitudes in the rest-frame J and H bands up to z=3.5 minimizing the dependence on the stellar evolution models. We compute the LF in the four redshift bins 1.5<z<2.0, 2.0<z<2.5, 2.5<z<3.0 and 3.0<z<3.5. Combining our results with those already available at lower redshifts, we find that (1) the faint end slope is consistent with being constant up to z=3.5, with alpha=-1.05+/-0.03 for the rest-frame J band and alpha=-1.15+/-0.02 for the rest-frame H band; (2) the normalization phi* decreases by a factor of 6 between z=0 and z~1.75 and by a factor 3 between z~1.75 and z=3.25; (3) the characteristic magnitude M* shows a brightening from z=0 to z~2 followed by a slow dimming to z=3.25. We finally compute the luminosity density (LD) in both rest-frame J and H bands. The analysis of our results together with those available in the literature shows that the LD is approximately constant up to z~1, and it then decreases by a factor of 6 up to z=3.5.

We report the discovery of a massive ultra-compact quiescent galaxy that has been strongly-lensed into multiple images by a foreground galaxy at z = 0.960. This system was serendipitously discovered as a set of extremely Ks-bright high-redshift galaxies with red J - Ks colors using new data from the UltraVISTA YJHKs near-infrared survey. The system was also previously identified as an optically-faint lens/source system using the COSMOS ACS imaging by Faure et al. (2008, 2011). Photometric redshifts for the three brightest images of the source galaxy determined from twenty-seven band photometry place the source at z = 2.4 +/- 0.1. We provide an updated lens model for the system which is a good fit to the positions and morphologies of the galaxies in the ACS image. The lens model implies that the magnification of the three brightest images is a factor of 4 - 5. We use the lens model, combined with the Ks-band image to constrain the size and Sersic profile of the galaxy. The best-fit model is an ultra-compact galaxy (Re = 0.64^{+0.08}_{-0.18} kpc, lensing-corrected), with a Sersic profile that is intermediate between a disk and bulge profile (n = 2.2^{+2.3}_{-0.9}). We present aperture photometry for the source galaxy images which have been corrected for flux contamination from the central lens. The best-fit stellar population model is a massive galaxy (Log(M_{star}/M_{sol}) = 10.8^{+0.1}_{-0.1}, lensing-corrected) with an age of 1.0^{+1.0}_{-0.4} Gyr, moderate dust extinction (Av = 0.8^{+0.5}_{-0.6}), and a low specific star formation rate (Log(SSFR) < -11.0 yr^{-1}). This is typical of massive "red-and-dead" galaxies at this redshift and confirms that this source is the first bona fide strongly-lensed massive ultra-compact quiescent galaxy to be discovered. We conclude with a discussion of the prospects of finding a larger sample of these galaxies.

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

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

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

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

Aims: We present a new method that uses luminosity or stellar mass functions combined with clustering measurements to select samples of galaxies at different redshifts likely to follow a progenitor-to-descendant relationship. As the method uses clustering information, we refer to galaxy samples selected this way as clustering-selected samples. We apply this method to infer the number of mergers during the evolution of MUSYC early-type galaxies (ETGs) from z~1 to the present-day. Methods: The method consists in using clustering information to infer the typical dark-matter halo mass of the hosts of the selected progenitor galaxies. Using LambdaCDM predictions, it is then possible to follow these haloes to a later time where the sample of descendants will be that with the clustering of these descendant haloes. Results: This technique shows that ETGs at a given redshift evolve into brighter galaxies at lower redshifts (considering rest-frame, passively evolved optical luminosities). This indicates that the stellar mass of these galaxies increases with time and that, in principle, a stellar mass selection at different redshifts does not provide samples of galaxies in a progenitor-descendant relationship. Conclusions: The comparison between high redshift ETGs and their likely descendants at z=0 points to a higher number density for the progenitors by a factor 5.5+-4.0, implying the need for mergers to decrease their number density by today. Because the luminosity densities of progenitors and descendants are consistent, our results show no need for significant star-formation in ETGs since z=1, which indicates that the needed mergers are dry, i.e. gas free.

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

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

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

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

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

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

We use ~88 arcmin**2 of deep (>~26.5 mag at 5 sigma) NICMOS data over the two GOODS fields and the HDF South to conduct a search for bright z>~7 galaxy candidates. This search takes advantage of an efficient preselection over 58 arcmin**2 of NICMOS H-band data where only plausible z>~7 candidates are followed up with NICMOS J-band observations. ~248 arcmin**2 of deep ground-based near-infrared data (>~25.5 mag, 5 sigma) is also considered in the search. In total, we report 15 z-dropout candidates over this area -- 7 of which are new to these search fields. Two possible z~9 J-dropout candidates are also found, but seem unlikely to correspond to z~9 galaxies. The present z~9 search is used to set upper limits on the prevalence of such sources. Rigorous testing is undertaken to establish the level of contamination of our selections by photometric scatter, low mass stars, supernovae (SNe), and spurious sources. The estimated contamination rate of our z~7 selection is ~24%. Through careful simulations, the effective volume available to our z>~7 selections is estimated and used to establish constraints on the volume density of luminous (L*(z=3), or -21 mag) galaxies from these searches. We find that the volume density of luminous star-forming galaxies at z~7 is 13_{-5}^{+8}x lower than at z~4 and >25x lower (1 sigma) at z~9 than at z~4. This is the most stringent constraint yet available on the volume density of >~L* galaxies at z~9. The present wide-area, multi-field search limits cosmic variance to <20%. The evolution we find at the bright end of the UV LF is similar to that found from recent Subaru Suprime-Cam, HAWK-I or ERS WFC3/IR searches. The present paper also includes a complete summary of our final z~7 z-dropout sample (18 candidates) identified from all NICMOS observations to date (over the two GOODS fields, the HUDF, galaxy clusters).

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

We present MIPS 24 micron observations of the Hubble Deep Field South taken with the Spitzer Space Telescope. The resulting image is 254 arcmin^2 in size and has a sensitivity ranging between ~12 to ~30 microJy rms, with a median sensitivity of ~20 microJy rms. A total of 495 sources have been cataloged with a signal-to-noise ratio greater than 5 sigma. The source catalog is presented as well as source counts which have been corrected for completeness and flux boosting. The IR sources are then combined with MUSYC optical/NIR and ATHDFS radio observations to obtain redshifts and radio flux densities of the sample. We use the IR/radio flux density ratio (q_24) to explore the IR-radio correlation for this IR sample and find q_24 = 0.71 +- 0.31 for sources detected in both IR and radio. The results are extended by stacking IR sources not detected in the radio observations and we derive an average q_24 for redshift bins between 0 < z < 2.5. We find the high redshift (z > 1) sources have an average q_{24} ratio which is better fit by local LIRG SEDs rather than local ULIRG SEDs, indicating that high redshift ULIRGs differ in their IR/radio properties. So ULIRGs at high redshift have SEDs different from those found locally. Infrared faint radio sources are examined, and while nine radio sources do not have a MIPS detection and are therefore radio-loud AGN, only one radio source has an extreme IRAC 3.6 micron to radio flux density ratio indicating it is a radio-loud AGN at z > 1.

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

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

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

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

We present a new maximum likelihood method for the calculation of galaxy luminosity functions from multi-band photometric surveys without spectroscopic data. The method evaluates the likelihood of a trial luminosity function by directly comparing the predicted distribution of fluxes in a multi-dimensional photometric space to the observations, and thus does not require the intermediate step of calculating photometric redshifts. We apply this algorithm to ~27,000 galaxies with m_R<=25 in the MUSYC-ECDFS field, with a focus on recovering the luminosity function of field galaxies at z<1.2. Our deepest LFs reach M_r=-14 and show that the field galaxy LF deviates from a Schechter function, exhibiting a steep upturn at intermediate magnitudes that is due to galaxies of late spectral types.

We study the relation between size and star formation activity in a complete sample of 225 massive (M > 5 x 10^10 Msun) galaxies at 1.5<z<2.5, selected from the FIREWORKS UV-IR catalog of the CDFS. Based on stellar population synthesis model fits to the observed restframe UV-NIR SEDs, and independent MIPS 24 micron observations, 65% of galaxies are actively forming stars, while 35% are quiescent. Using sizes derived from 2D surface brightness profile fits to high resolution (FWHM_{PSF}~0.45 arcsec) groundbased ISAAC data, we confirm and improve the significance of the relation between star formation activity and compactness found in previous studies, using a large, complete mass-limited sample. At z~2, massive quiescent galaxies are significantly smaller than massive star forming galaxies, and a median factor of 0.34+/-0.02 smaller than galaxies of similar mass in the local universe. 13% of the quiescent galaxies are unresolved in the ISAAC data, corresponding to sizes <1 kpc, more than 5 times smaller than galaxies of similar mass locally. The quiescent galaxies span a Kormendy relation which, compared to the relation for local early types, is shifted to smaller sizes and brighter surface brightnesses and is incompatible with passive evolution. The progenitors of the quiescent galaxies, were likely dominated by highly concentrated, intense nuclear star bursts at z~3-4, in contrast to star forming galaxies at z~2 which are extended and dominated by distributed star formation.

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.

[Abridged] We present the evolution of the stellar mass function (SMF) of galaxies from z=4.0 to z=1.3 measured from a sample constructed from the deep NIR MUSYC, the FIRES, and the GOODS-CDFS surveys, all having very high-quality optical to mid-infrared data. This sample, unique in that it combines data from surveys with a large range of depths and areas in a self-consistent way, allowed us to 1) minimize the uncertainty due to cosmic variance and empirically quantify its contribution to the total error budget; 2) simultaneously probe the high-mass end and the low-mass end (down to ~0.05 times the characteristic stellar mass) of the SMF with good statistics; and 3) empirically derive the redshift-dependent completeness limits in stellar mass. We provide, for the first time, a comprehensive analysis of random and systematic uncertainties affecting the derived SMFs. We find that the mass density evolves by a factor of ~17(+7,-10) since z=4.0, mostly driven by a change in the normalization Phi*. If only random errors are taken into account, we find evidence for mass-dependent evolution, with the low-mass end evolving more rapidly than the high-mass end. However, we show that this result is no longer robust when systematic uncertainties due to the SED-modeling assumptions are taken into account. Taking our results at face value, we find that they are in conflict with semi-analytic models of galaxy formation. The models predict SMFs that are in general too steep, with too many low-mass galaxies and too few high-mass galaxies. The discrepancy at the high-mass end is susceptible to uncertainties in the models and the data, but the discrepancy at the low-mass end may be more difficult to explain.

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

We present a comparison between the observed color distribution, number and mass density of massive galaxies at 1.5 < z < 3 and a model by Hopkins et al. that relates the quasar and galaxy population on the basis of gas-rich mergers. In order to test the hypothesis that quiescent red galaxies are formed after a gas-rich merger involving quasar activity, we confront photometry of massive (M > 4x10^10 Msun) galaxies extracted from the FIRES, GOODS-South, and MUSYC surveys, together spanning an area of 496 arcmin^2, with synthetic photometry from hydrodynamical merger simulations. As in the Hopkins et al. (2006b) model, we use the observed quasar luminosity function to estimate the merger rate. We find that the synthetic U-V and V-J colors of galaxies that had a quasar phase in their past match the colors of observed galaxies that are best characterized by a quiescent stellar population. At z ~ 2.6, the observed number and mass density of quiescent red galaxies with M > 4x10^10 Msun is consistent with the model in which every quiescent massive galaxy underwent a quasar phase in the past. At z ~ 1.9, 2.8 times less quiescent galaxies are observed than predicted by the model as descendants of higher redshift quasars. The merger model also predicts a large number of galaxies undergoing merger-driven star formation. We find that the predicted number and mass density accounts for 30-50% of the observed massive star-forming galaxies. However, their colors do not match those of observed star-forming galaxies. In particular, the colors of dusty red galaxies are not reproduced by the simulations. Several possible origins of this discrepancy are discussed. The observational constraints on the validity of the model are currently limited by cosmic variance and uncertainties in stellar population synthesis and radiative transfer.

We present the rest-frame optical luminosity function (LF) of red sequence galaxies in 16 clusters at 0.4<z<0.8 drawn from the ESO Distant Cluster Survey (EDisCS). We compare our clusters to an analogous sample from the Sloan Digital Sky Survey (SDSS) and match the EDisCS clusters to their most likely descendants. We measure all LFs down to M M* + (2.5 - 3.5). At z<0.8, the bright end of the LF is consistent with passive evolution but there is a significant build-up of the faint end of the red sequence towards lower redshift. There is a weak dependence of the LF on cluster velocity dispersion for EDisCS but no such dependence for the SDSS clusters. We find tentative evidence that red sequence galaxies brighter than a threshold magnitude are already in place, and that this threshold evolves to fainter magnitudes toward lower redshifts. We compare the EDisCS LFs with the LF of co-eval red sequence galaxies in the field and find that the bright end of the LFs agree. However, relative to the number of bright red galaxies, the field has more faint red galaxies than clusters at 0.6<z<0.8 but fewer at 0.4<z<0.6, implying differential evolution. We compare the total light in the EDisCS cluster red sequences to the total red sequence light in our SDSS cluster sample. Clusters at 0.4<z<0.8 must increase their luminosity on the red sequence (and therefore stellar mass in red galaxies) by a factor of 1-3 by z=0. The necessary processes that add mass to the red sequence in clusters predict local clusters that are over-luminous as compared to those observed in the SDSS. The predicted cluster luminosities can be reconciled with observed local cluster luminosities by combining multiple previously known effects.

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

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

We compare recently measured rest-frame V-band luminosity functions (LFs) of galaxies at redshifts 2.0<z<3.3 to predictions of semianalytic models by De Lucia & Blaizot and Bower et al. and hydrodynamic simulations by Dave et al. The models succeed for some luminosity and redshift ranges and fail for others. A notable success is that the Bower et al. model provides a good match to the observed LF at z~3. However, all models predict an increase with time of the rest-frame V-band luminosity density, whereas the observations show a decrease. The models also have difficulty matching the observed rest-frame colors of galaxies. In all models the luminosity density of red galaxies increases sharply from z~3 to z~2.2, whereas it is approximately constant in the observations. Conversely, in the models the luminosity density of blue galaxies is approximately constant, whereas it decreases in the observations. These discrepancies cannot be entirely remedied by changing the treatment of dust and suggest that current models do not yet provide an adequate description of galaxy formation and evolution since z~3.

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.

Recent galaxy formation models successfully reproduce the local luminosity function (LF) of galaxies by invoking mechanisms to suppress star formation in low- and high-mass galaxies. As these models are optimized to fit the LF at low redshift, a crucial question is how well they predict the LF at earlier times. Here we compare recently measured rest-frame V-band LFs of galaxies at redshifts 2.0<z<3.3 to predictions of semianalytic models by De Lucia & Blaizot and Bower et al. and hydrodynamic simulations by Dave et al.. The models succeed for some luminosity and redshift ranges and fail for others. A notable success is that the Bower et al. model provides a good match to the observed LF at z~3. However, all models predict an increase with time of the rest-frame V-band luminosity density, whereas the observations show a decrease. The models also have difficulty matching the observed rest-frame colors of galaxies. In all models the luminosity density of red galaxies increases sharply from z~3 to z~2.2, whereas it is approximately constant in the observations. Conversely, in the models the luminosity density of blue galaxies is approximately constant, whereas it decreases in the observations. These discrepancies cannot be entirely remedied by changing the treatment of dust and suggest that current models do not yet provide a complete description of galaxy formation and evolution since z~3.

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

We describe the results of an extremely deep, 0.28 deg^2 survey for z = 3.1 Ly-alpha emission-line galaxies in the Extended Chandra Deep Field South. By using a narrow-band 5000 Anstrom filter and complementary broadband photometry from the MUSYC survey, we identify a statistically complete sample of 162 galaxies with monochromatic fluxes brighter than 1.5 x 10^-17 ergs cm^-2 s^-1 and observers frame equivalent widths greater than 80 Angstroms. We show that the equivalent width distribution of these objects follows an exponential with a rest-frame scale length of w_0 = 76 +/- 10 Angstroms. In addition, we show that in the emission line, the luminosity function of Ly-alpha galaxies has a faint-end power-law slope of alpha = -1.49 +/- 0.4, a bright-end cutoff of log L^* = 42.64 +/- 0.2, and a space density above our detection thresholds of 1.46 +/- 0.12 x 10^-3 h70^3 galaxies Mpc^-3. Finally, by comparing the emission-line and continuum properties of the LAEs, we show that the star-formation rates derived from Ly-alpha are ~3 times lower than those inferred from the rest-frame UV continuum. We use this offset to deduce the existence of a small amount of internal extinction within the host galaxies. This extinction, coupled with the lack of extremely-high equivalent width emitters, argues that these galaxies are not primordial Pop III objects, though they are young and relatively chemically unevolved.

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 analyze the near-infrared luminosity profiles and photometric parameters of the host galaxies of 3CR radio sources with z<0.3, to investigate their physical nature. Our sample includes 82 galaxies, of which 22 (27%) are FR Is and 60 (73%) are FR IIs. Using near-infrared data taken both with NICMOS onboard the Hubble Space Telescope and from the ground with the Telescopio Nazionale Galileo, we find that luminosity profiles are very well described by a single Sersic law in 52% of the cases and for the remaining objects (48%) it is necessary to include an exponential profile, which could indicate the presence of a disk. The average bulge to disk luminosity ratio for the galaxies is (b/d) ~ 1.1. The analysis of the photometric parameters of the sub samples indicates that FR Is and FR IIs show rather similar bulges in terms of effective surface magnitude, effective radius, and Sersic index. On the other hand, disks in FR Is and FR IIs hosts show, on average, different properties. Central surface magnitudes are dimmer and scale lengths are greater by a factor of 2 in FR Is when compared to FR IIs. We also estimate the black hole mass associated with each galaxy using two different methods that claim tight correlations between the black hole mass (M_BH) with the infrared bulge luminosity (L_bulge) and with the Sersic index (n). Our data indicate that masses obtained through these two methods show a high dispersion and M_BH obtained through L_bulge are systematically higher (by a factor of ~3) than those obtained using n. This result may reflect the fact that for our sample galaxies we do not find any correlation between L_bulge and n.

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 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 the first near-infrared (K'-band) homogeneous observations of a complete sub-sample of the 3CR radio catalogue comprising all High Excitation Galaxies (HEGs) at z<0.3. After showing that the surface brightness decomposition technique to measure central point-like sources is affected by significant uncertainties for the objects in the studied sample, we present a new, more accurate method based on the R-K' color profile. Via this method we find a substantial nuclear K'-band excess in all but two HEGs -- most likely directly associated to their nuclear emission -- and we measure the corresponding 2.12 $\mu$m nuclear luminosities. Within the frame of the unification scheme for radio-loud active galactic nuclei, it appears that obscuration alone is not able to account for the different nuclear properties of the majority of the HEGs and Broad Line Radio Galaxies (BLRGs), and also scattering of the (optically) hidden nuclear light from a compact region must be invoked. More precisely, for ~70% of the HEGs the observed point-like optical emission is dominated by the scattered component, while in the K'-band both scattered and direct light passing through the torus contribute to the observed nuclear luminosity. The estimated fraction of scattered light ranges from a few tenths to a few percent, while the torus extinction is between 15<A_{V,torus}<50 mag with only a few exceptions with lower obscuration.

We present evidence for the presence of a transition in the accretion properties of radio-loud sources. For a sample of radio galaxies and radio-loud quasars, selected based on their extended radio properties, the accretion rate is estimated from the black hole mass and nuclear luminosity. The inferred distribution is bimodal, with paucity of sources at accretion rates, in Eddington units, of order ~ 10^{-2} - assuming a radiative efficiency of 10 per cent - and possibly spanning one-two orders of magnitude. Selection biases are unlikely to be responsible for such behavior; we discuss possible physical explanations, including a fast transition to low accretion rates, a change in the accretion mode/actual accretion rate/radiative efficiency, the lack of stable disc solutions at intermediate accretion rates or the inefficiency of the jet formation processes in geometrically thin flows. This transition might be analogous to spectral states (and jet) transitions in black hole binary systems.

We present high resolution Halpha rotation curves of 4 late-type dwarf galaxies and 2 low surface brightness galaxies (LSB) for which accurate HI rotation curves are available from the literature. Observations are carried out at Telescopio Nazionale Galileo (TNG). For LSB F583-1 an innovative dispersing element was used, the Volume Phase Holographic (VPH) with a dispersion of about 0.35 A/pxl. We find good agreement between the Halpha data and the HI observations and conclude that the HI data for these galaxies suffer very little from beam smearing. We show that the optical rotation curves of these dark matter dominated galaxies are best fitted by the Burkert profile. In the centers of galaxies, where the N-body simulations predict cuspy cores and fast rising rotation curves, our data seem to be in better agreement with the presence of soft cores.

We present high spatial resolution H\alpha rotation curves of late-type dwarf and LSB galaxies. From our analysis we find good agreement between our H\alpha data and the HI observations taken from the literature, concluding that the HI rotation curves for these galaxies suffer very little from beam smearing. A preliminary analysis of our data rules out the CDM model in the inner regions of these galaxies.

This work is focused on the preliminary results of the observations of H\alpha rotation curves for some of the late type dwarf and LSB galaxies carried out at the TNG telescope. In light of the observational data and of the N-body simulations we have performed recently, we discuss some of the implications on the nature of the DM particles and the formation of the dark haloes.