We use deep, spatially resolved spectroscopy from the LEGA-C Survey to study radial variations in the stellar population of 17 spectroscopically-selected post-starburst (PSB) galaxies. We use spectral fitting to measure two Lick indices, $H{\delta}_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $H{\delta}_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $H{\delta}_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ \r{A}/$R_e$), mainly due to seeing convolution. A two-SSP model suggests intrinsic gradients are as strong as observed in local PSB galaxies ($\approx -0.8$ \r{A}$/R_e$). We interpret these results in terms of inside-out growth (for the bulk of the quiescent population) vs star formation occurring last in the centre (for PSB galaxies). At $z\approx0.8$, central starbursts are often the result of gas-rich mergers, as evidenced by the high fraction of PSB galaxies with disturbed morphologies and tidal features (40%). Our results provide additional evidence for multiple paths to quiescence: a standard path, associated with inside-out disc formation and with gradually decreasing star-formation activity, without fundamental structural transformation, and a fast path, associated with centrally-concentrated starbursts, leaving an inverse age gradient and smaller half-light radius.

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 results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs), we find that cluster galaxies have suppressed SFRs at fixed stellar mass in comparison to their field counterparts by a factor of 1.4 $\pm$ 0.1 ($\sim3.3\sigma$) across the stellar mass range: $9.0 < \log(M_{*} /M_{\odot}) < 11.2$. We also find that this modest suppression in the cluster galaxy star-forming main sequence is mass and redshift dependent: the difference between cluster and field increases towards lower stellar masses and lower redshift. When comparing the distribution of cluster and field galaxy SFRs to the star-forming main sequence, we find an overall shift towards lower SFRs in the cluster population, and note the absence of a tail of high SFR galaxies as seen in the field. Given this observed suppression in the cluster galaxy star-forming main sequence, we explore the implications for several scenarios such as formation time differences between cluster and field galaxies, and environmentally-induced star formation quenching and associated timescales.

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.

In the hierarchical structure formation model of the universe, galaxy clusters are assembled through a series of mergers. Accordingly, it is expected that galaxy clusters in the early universe are actively forming and dynamically young. Located at a high redshift of z=1.71, SpARCS1049+56 offers a unique look into the galaxy cluster formation process. This cluster has been shown to be rich in cluster galaxies and to have intense star formation. Its high redshift pushes a weak-lensing analysis beyond the regime of the optical spectrum into that of the infrared. Equipped with deep Hubble Space Telescope Wide Field Camera 3 UVIS and IR observations, we present a weak-lensing characterization of SpARCS1049+56. As few IR weak-lensing studies have been performed, we discuss the details of PSF modeling and galaxy shape measurement for an IR weak-lensing procedure and the systematics that come with the territory. It will be critical to understand these systematics in future weak-lensing studies in the IR with the next generation space telescopes such as JWST, Euclid, and WFIRST. Through a careful analysis, the mass distribution of this young galaxy cluster is mapped and the convergence peak is detected at a 3.3 sigma level. The weak-lensing mass of the cluster is estimated to be $3.5\pm1.2\times10^{14}\ \text{M}_\odot$ and is consistent with the mass derived from a mass-richness scaling relation. This mass is extreme for a cluster at such a high redshift and suggests that SpARCS1049+56 is rare in the standard $\Lambda$CDM universe.

In this letter, we investigate the impact of environment on integrated and spatially-resolved stellar kinematics of a sample of massive, quiescent galaxies at intermediate redshift ($0.6<z<1.0$). For this analysis, we combine photometric and spectroscopic parameters from the UltraVISTA and Large Early Galaxy Astrophysics Census (LEGA-C) surveys in the COSMOS field and environmental measurements. We analyze the trends with overdensity (1+$\delta$) on the rotational support of quiescent galaxies and find no universal trends at either fixed mass or fixed stellar velocity dispersion. This is consistent with previous studies of the local Universe; rotational support of massive galaxies depends primarily on stellar mass. We highlight two populations of massive galaxies ($\log M_\star/M_\odot\geq11$) that deviate from the average mass relation. First, the most massive galaxies in the most under-dense regions ($(1+\delta)\leq1$) exhibit elevated rotational support. Similarly, at the highest masses ($\log M_\star/M_\odot\geq11.25$) the range in rotational support is significant in all but the densest regions. This corresponds to an increasing slow-rotator fraction such that only galaxies in the densest environments ($(1+\delta)\geq3.5$) are primarily (90$\pm$10\%) slow-rotators.This effect is not seen at fixed velocity dispersion, suggesting minor merging as the driving mechanism: only in the densest regions have the most massive galaxies experienced significant minor merging, building stellar mass and diminishing rotation without significantly affecting the central stellar velocity dispersion. In the local Universe, most massive galaxies are slow-rotators, regardless of environment, suggesting minor merging occurs at later cosmic times $(z\lesssim0.6)$ in all but the most dense environments.

We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(H\delta) \geq 4$\AA) at $z\sim0.8$ drawn from DR2 of the LEGA-C survey to test the hypothesis that these galaxies experienced compact, central starbursts before quenching. These recently quenched galaxies, usually referred to as post-starburst galaxies, span a wide range of colors and we find a clear correlation between color and half-light radius, such that bluer galaxies are smaller. We build simple toy models to explain this correlation: a normal star-forming disk plus a central, compact starburst component. Bursts with exponential decay timescale of $\sim$~100 Myr that produce $\sim10\%$ to more than 100\% of the pre-existing masses can reproduce the observed correlation. More significant bursts also produce bluer and smaller descendants. Our findings imply that when galaxies shut down star formation rapidly, they generally had experienced compact, starburst events and that the large, observed spread in sizes and colors mostly reflects a variety of burst strengths. Recently quenched galaxies should have younger stellar ages in the centers; multi-wavelength data with high spatial resolution are required to reveal the age gradient. Highly dissipative processes should be responsible for this type of formation history. While determining the mechanisms for individual galaxies is challenging, some recently quenched galaxies show signs of gravitational interactions, suggesting that mergers are likely an important mechanism in triggering the rapid shut-down of star-formation activities at $z\sim0.8$.

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.

We use reconstructed star-formation histories (SFHs) of quiescent galaxies at $z=0.6-1$ in the LEGA-C survey to identify secondary star-formation episodes that, after an initial period of quiescence, moved the galaxies back to the star-forming main sequence (blue cloud). $16\pm3$\% of the $z\sim0.8$ quiescent population has experienced such rejuvenation events in the redshift range $0.7<z<1.5$ after reaching quiescence at some earlier time. On average, these galaxies first became quiescent at $z=1.2$, and those that rejuvenated, remained quiescent for $\sim1$Gyr before their secondary SF episode which lasted $\sim0.7$Gyr. The stellar mass attributed to rejuvenation is on average 10\% of the galaxy stellar mass, with rare instances of an increase of more than a factor 2. Overall, rejuvenation events only contribute $\sim2$\% of the total stellar mass in $z\sim0.8$ quiescent galaxies and we conclude that rejuvenation is not an important evolutionary channel when considering the growth of the red sequence. However, our results complicate the interpretation of galaxy demographics in color space: the galaxies with rejuvenation events tend to lie in the so-called `green valley', yet their progenitors were quiescent at $z\sim2$.

We present deep, multiwavelength radio observations of SpARCS104922.6+564032.5, a z = 1.71 galaxy cluster with a starbusting core. Observations were made with the Karl G. Jansky Very Large Array (JVLA) in 3 bands: 1-2 GHz, 4-8 GHz and 8-12 GHz. We detect a radio source coincident with the Brightest Cluster Galaxy (BCG) that has a spectral index of {\alpha}=0.44\pm 0.29 and is indicative of emission from an Active Galactic Nucleus. The radio luminosity is consistent with the average luminosity of the lower redshift BCG sample, but the flux densities are 6{\sigma} below the predicted values of the star-forming Spectral Energy Distribution based on far infrared data. Our new fit fails to simultaneously describe the far infrared and radio fluxes. This, coupled with the fact that no other bright source is detected in the vicinity of the BCG implies that the star formation region, traced by the infrared emission, is extended or clumpy and not located directly within the BCG. Thus, we suggest that the star-forming core might not be driven by a single major wet merger, but rather by several smaller galaxies stripped of their gas or by a displaced cooling flow, although more data are needed to confirm any of those scenarios.

We study the star-forming (SF) population of galaxies within a sample of 209 IR-selected galaxy clusters at 0.3$\,\leq\,z\,\leq\,$1.1 in the ELAIS-N1 and XMM-LSS fields, exploiting the first HSC-SSP data release. The large area and depth of these data allows us to analyze the dependence of the SF fraction, $f_{SF}$, on stellar mass and environment separately. Using $R/R_{200}$ to trace environment, we observe a decrease in $f_{SF}$ from the field towards the cluster core, which strongly depends on stellar mass and redshift. The data show an accelerated growth of the quiescent population within the cluster environment: the $f_{SF}$ vs. stellar mass relation of the cluster core ($R/R_{200}\,\leq\,$0.4) is always below that of the field (4$\,\leq\,R/R_{200}\,<\,$6). Finally, we find that environmental and mass quenching efficiencies depend on galaxy stellar mass and distance to the center of the cluster, demonstrating that the two effects are not separable in the cluster environment. We suggest that the increase of the mass quenching efficiency in the cluster core may emerge from an initial population of galaxies formed ``in situ.'' The dependence of the environmental quenching efficiency on stellar mass favors models in which galaxies exhaust their reservoir of gas through star formation and outflows, after new gas supply is truncated when galaxies enter the cluster.

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

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.

We present a comparison of the observed, spatially integrated stellar and ionized gas velocity dispersions of $\sim1000$ massive ($\log M_{\star}/M_{\odot}\gtrsim\,10.3$) galaxies in the Large Early Galaxy Astrophysics Census (LEGA-C) survey at $0.6\lesssim\,z\lesssim1.0$. The high $S/N\sim20{\rm\AA^{-1}}$ afforded by 20 hour VLT/VIMOS spectra allows for joint modeling of the stellar continuum and emission lines in all galaxies, spanning the full range of galaxy colors and morphologies. These observed integrated velocity dispersions (denoted as $\sigma'_{g, int}$ and $\sigma'_{\star, int}$) are related to the intrinsic velocity dispersions of ionized gas or stars, but also include rotational motions through beam smearing and spectral extraction. We find good average agreement between observed velocity dispersions, with $\langle\log(\sigma'_{g, int}/\sigma'_{\star, int})\rangle=-0.003$. This result does not depend strongly on stellar population, structural properties, or alignment with respect to the slit. However, in all regimes we find significant scatter between $\sigma'_{g, int}$ and $\sigma'_{\star, int}$, with an overall scatter of 0.13 dex of which 0.05 dex is due to observational uncertainties. For an individual galaxy, the scatter between $\sigma'_{g, int}$ and $\sigma'_{\star, int}$ translates to an additional uncertainty of $\sim0.24\rm{dex}$ on dynamical mass derived from $\sigma'_{g, int}$, on top of measurement errors and uncertainties from Virial constant or size estimates. We measure the $z\sim0.8$ stellar mass Faber-Jackson relation and demonstrate that emission line widths can be used to measure scaling relations. However, these relations will exhibit increased scatter and slopes that are artificially steepened by selecting on subsets of galaxies with progressively brighter emission lines.

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 analyze stellar age indicators (D$_n$4000 and EW(H$\delta$)) and sizes of 467 quiescent galaxies with $M_\ast \geq 10^{10} M_\odot$ at $z\sim0.7$ drawn from DR2 of the LEGA-C survey. Interpreting index variations in terms of equivalent single stellar population age, we find that the median stellar population is younger for larger galaxies at fixed stellar mass. The effect is significant, yet small; the ages of the larger and the smaller subsets differ by only $<500$ Myr, much less than the age variation among individual galaxies ($\sim1.5$ Gyr). At the same time, quiescent galaxies with the strongest H$\delta$ absorption --- those experienced recent and rapid quenching events --- tend to be smaller than the average. These co-existing trends unify seemingly contradictory results in the literature; the complex correlations between size and age indicators revealed by our large sample of galaxies with high-quality spectra suggest that there are multiple evolutionary pathways to quiescence. Regardless of the specific physical mechanisms responsible for the cessation of star formation in massive galaxies, the large scatter in D$_n$4000 and EW(H$\delta$) immediately implies that galaxies follow a large variety in evolutionary pathways. On the one hand, we see evidence for a process that slowly shuts off star-formation and transforms star-forming galaxies to quiescent galaxies without necessarily changing their structures. On the other hand, there is likely a mechanism that rapidly quenches galaxies, an event that coincides with dramatic structural changes, producing small post-starburst galaxies.

We present the second data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO 130-night public spectroscopic survey conducted with VIMOS on the Very Large Telescope. We release 1988 spectra with typical continuum S / N ~= 20 /Angstrom of galaxies at 0.6 ~< z ~< 1.0, each observed for ~20 hours and fully reduced with a custom-built pipeline. We also release a catalog with spectroscopic redshifts, emission line fluxes, Lick/IDS indices, and observed stellar and gas velocity dispersions that are spatially integrated quantities including both rotational motions and genuine dispersion. To illustrate the new parameter space in the intermediate redshift regime probed by LEGA-C we explore relationships between dynamical and stellar population properties. The star-forming galaxies typically have observed stellar velocity dispersions of ~150 km/s and strong Hdelta absorption (Hd_A ~ 5 Angstrom), while passive galaxies have higher observed stellar velocity dispersions ~200 km/s and weak Hdelta absortion (Hd_A ~ 0 Angstrom). Strong [O III]5007 / Hbeta ratios tend to occur mostly for galaxies with weak Hd_A or galaxies with higher observed velocity dispersion. Beyond these broad trends, we find a large diversity of possible combinations of rest-frame colors, absorption line strengths and emission line detections, illustrating the utility of spectroscopic measurements to more accurately understand galaxy evolution. By making the spectra and value-added catalogs publicly available we encourage the community to take advantage of this very substantial investment in telescope time provided by ESO.

We present the first spatially-resolved observations of molecular gas in a sample of cluster galaxies beyond z>0.1. Using ALMA, we detect CO (2-1) in 8 z~1.6 cluster galaxies, all within a single 70" primary beam, in under 3 hours of integration time. The cluster, SpARCS-J0225, is replete with gas-rich galaxies in close proximity. It thus affords an efficient multiplexing strategy to build up the first sample of resolved CO in distant galaxy clusters. Mapping out the kinematic structure and morphology of the molecular gas on 3.5 kpc scales reveals rotating gas disks in the majority of the galaxies, as evidenced by smooth velocity gradients. Detailed velocity maps also uncover kinematic peculiarities, including a central gas void, a merger, and a few one-sided gas tails. We compare the extent of the molecular gas component to that of the optical stellar component, measured with rest-frame optical HST imaging. We find that the cluster galaxies, while broadly consistent with a ratio of unity for stellar-to-gas effective radii, have a moderately larger ratio compared to the coeval field; this is consistent with the more pronounced trend in the low-redshift Universe. Thus, at first glance, the z~1.6 cluster galaxies generally look like galaxies infalling from the field, with typical main-sequence star formation rates and massive molecular gas reservoirs situated in rotating disks. However, there are potentially important differences from their field counterparts, including elevated gas fractions, slightly smaller CO disks, and possible asymmetric gas tails. Taken in tandem, these signatures are tentative evidence for gas-stripping in the z~1.6 cluster. However, the current sample size of spatially-resolved molecular gas in galaxies at high redshift is small, and verification of these trends will require much larger samples of both cluster and field galaxies.

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.

The Gemini Infrared Multi-Object Spectrograph (GIRMOS) is a powerful new instrument being built to facility-class standards for the Gemini telescope. It takes advantage of the latest developments in adaptive optics and integral field spectrographs. GIRMOS will carry out simultaneous high-angular-resolution, spatially-resolved infrared ($1-2.4$ $\mu$m) spectroscopy of four objects within a two-arcminute field-of-regard by taking advantage of multi-object adaptive optics. This capability does not currently exist anywhere in the world and therefore offers significant scientific gains over a very broad range of topics in astronomical research. For example, current programs for high redshift galaxies are pushing the limits of what is possible with infrared spectroscopy at $8-10$-meter class facilities by requiring up to several nights of observing time per target. Therefore, the observation of multiple objects simultaneously with adaptive optics is absolutely necessary to make effective use of telescope time and obtain statistically significant samples for high redshift science. With an expected commissioning date of 2023, GIRMOS's capabilities will also make it a key followup instrument for the James Webb Space Telescope when it is launched in 2021, as well as a true scientific and technical pathfinder for future Thirty Meter Telescope (TMT) multi-object spectroscopic instrumentation. In this paper, we will present an overview of this instrument's capabilities and overall architecture. We also highlight how this instrument lays the ground work for a future TMT early-light instrument.

This paper describes new deep 3.6 and 4.5 micron imaging of three UltraVISTA near-infrared survey stripes within the COSMOS field. The observations were carried out with Spitzer's Infrared Array Camera (IRAC) for the Spitzer Matching Survey of the Ultra-VISTA Deep Stripes (SMUVS). In this work we present our data reduction techniques, and document the resulting mosaics, coverage maps, and catalogs in both IRAC passbands for the three easternmost UltraVISTA survey stripes, covering a combined area of about 0.66 square degrees, of which 0.45 square degrees have at least 20 hr integration time. SMUVS reaches point-source sensitivities of about 25.0 AB mag at both 3.6 and 4.5 microns with a significance of 4-sigma accounting for both survey sensitivity and source confusion. To this limit the SMUVS catalogs contain a total of about 350,000 sources, each of which is detected significantly in at least one IRAC band. Because of its uniform and high sensitivity, relatively large area coverage, and the wide array of ancillary data available in COSMOS, the SMUVS survey will be useful for a large number of cosmological investigations. We will make all images and catalogues described herein publicly available via the Spitzer Science Center.

We study the abundance of star-forming and quiescent galaxies in a sample of 21 massive clusters at 0.5<z<0.7, detected with the Planck satellite. We measure the cluster galaxy stellar mass function (SMF), which is a fundamental observable to study and constrain the formation and evolution of galaxies. Our measurements are based on homogeneous and deep multi-band photometry spanning u- to the Ks-band for each cluster and are supported by spectroscopic data from different programs. The galaxy population is separated between quiescent and star-forming galaxies based on their rest-frame U-V and V-J colours. The SMF is compared to that of field galaxies at the same redshifts, using data from the COSMOS/UltraVISTA survey. We find that the shape of the SMF of star-forming galaxies does not depend on environment, while the SMF of quiescent galaxies has a significantly steeper low-mass slope in the clusters compared to the field. We estimate the environmental quenching efficiency (f_EQ), i.e. the probability for a galaxy that would normally be star forming in the field, to be quenched due to its environment. The f_EQ shows no stellar-mass dependence in any environment, but it increases from 40% in the cluster outskirts to ~90% in the cluster centres. The radial signature of f_EQ provides constraints on where the dominant quenching mechanism operates in these clusters and on what timescale. Exploring these using a simple model based on galaxy orbits obtained from an N-body simulation, we find a clear degeneracy between both parameters. For example, the quenching process may either be triggered on a long (~3 Gyr) time scale at large radii (r~8R_500), or happen well within 1 Gyr at r<R_500. The radius where quenching is triggered is at least r_quench> 0.67R_500 (95%CL). The ICM density at this location suggests that ram-pressure stripping of the cold gas is a likely cause of quenching. [Abridged]

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

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

We present a comparison of two samples of z>0.8 galaxy clusters selected using different wavelength-dependent techniques and examine the physical differences between them. We consider 18 clusters from the X-ray selected XMM-LSS distant cluster survey and 92 clusters from the optical-MIR selected SpARCS cluster survey. Both samples are selected from the same approximately 9 square degree sky area and we examine them using common XMM-Newton, Spitzer-SWIRE and CFHT Legacy Survey data. Clusters from each sample are compared employing aperture measures of X-ray and MIR emission. We divide the SpARCS distant cluster sample into three sub-samples: a) X-ray bright, b) X-ray faint, MIR bright, and c) X-ray faint, MIR faint clusters. We determine that X-ray and MIR selected clusters display very similar surface brightness distributions of galaxy MIR light. In addition, the average location and amplitude of the galaxy red sequence as measured from stacked colour histograms is very similar in the X-ray and MIR-selected samples. The sub-sample of X-ray faint, MIR bright clusters displays a distribution of BCG-barycentre position offsets which extends to higher values than all other samples. This observation indicates that such clusters may exist in a more disturbed state compared to the majority of the distant cluster population sampled by XMM-LSS and SpARCS. This conclusion is supported by stacked X-ray images for the X-ray faint, MIR bright cluster sub-sample that display weak, centrally-concentrated X-ray emission, consistent with a population of growing clusters accreting from an extended envelope of material.

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

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.

Using a sample of 4 galaxy clusters at $1.35 < z < 1.65$ and 10 galaxy clusters at $0.85 < z < 1.35$, we measure the environmental quenching timescale, $t_Q$, corresponding to the time required after a galaxy is accreted by a cluster for it to fully cease star formation. Cluster members are selected by a photometric-redshift criterion, and categorized as star-forming, quiescent, or intermediate according to their dust-corrected rest-frame colors and magnitudes. We employ a "delayed-then-rapid" quenching model that relates a simulated cluster mass accretion rate to the observed numbers of each type of galaxy in the cluster to constrain $t_Q$. For galaxies of mass $M_* \gtrsim 10^{10.5}~ \mathrm{M}_\odot$, we find a quenching timescale of $t_Q=$ 1.24 Gyr in the $z\sim1.5$ cluster sample, and $t_Q=$ 1.50 Gyr at $z\sim1$. Using values drawn from the literature, we compare the redshift evolution of $t_Q$ to timescales predicted for different physical quenching mechanisms. We find $t_Q$ to depend on host halo mass such that quenching occurs over faster timescales in clusters relative to groups, suggesting that properties of the host halo are responsible for quenching high-mass galaxies. Between $z=0$ and $z=1.5$, we find that $t_Q$ evolves faster than the molecular gas depletion timescale and slower than an SFR-outflow timescale, but is consistent with the evolution of the dynamical time. This suggests that environmental quenching in these galaxies is driven by the motion of satellites relative to the cluster environment, although due to uncertainties in the atomic gas budget at high redshift, we cannot rule out quenching due to simple gas depletion.

Drawing from the LEGA-C dataset, we present the spectroscopic view of the stellar population across a large volume- and mass-selected sample of galaxies at large lookback time. We measure the 4000\AA\ break (D$_n$4000) and Balmer absorption line strengths (probed by H$\delta$) from 1019 high-quality spectra of $z=0.6 - 1.0$ galaxies with $M_\ast = 2 \times 10^{10} M_\odot - 3 \times 10^{11} M_\odot$. Our analysis serves as a first illustration of the power of high-resolution, high-S/N continuum spectroscopy at intermediate redshifts as a qualitatively new tool to constrain galaxy formation models. The observed D$_n$4000-EW(H$\delta$) distribution of our sample overlaps with the distribution traced by present-day galaxies, but $z\sim 0.8$ galaxies populate that locus in a fundamentally different manner. While old galaxies dominate the present-day population at all stellar masses $> 2\times10^{10} M_\odot$, we see a bimodal D$_n$4000-EW(H$\delta$) distribution at $z\sim0.8$, implying a bimodal light-weighted age distribution. The light-weighted age depends strongly on stellar mass, with the most massive galaxies $>1\times10^{11}M_\odot$ being almost all older than 2 Gyr. At the same time we estimate that galaxies in this high mass range are only $\sim3$ Gyr younger than their $z\sim0.1$ counterparts, at odd with pure passive evolution given a difference in lookback time of $>5$ Gyr; younger galaxies must grow to $>10^{11}M_\odot$ in the meantime, and/or small amounts of young stars must keep the light-weighted ages young. Star-forming galaxies at $z\sim0.8$ have stronger H$\delta$ absorption than present-day galaxies with the same D$_n$4000, implying larger short-term variations in star-formation activity.

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

The recent discovery of thousands of ultra diffuse galaxies (UDGs) in nearby galaxy clusters has opened a new window into the process of galaxy formation and evolution. Several scenarios have been proposed to explain the formation history of UDGs, and their ability to survive in the harsh cluster environments. A key requirement to distinguish between these scenarios is a measurement of their halo masses which, due to their low surface brightnesses, has proven difficult if one relies on stellar tracers of the potential. We exploit weak gravitational lensing, a technique that does not depend on these baryonic tracers, to measure the average subhalo mass of 784 UDGs selected in 18 clusters at $z\leq0.09$. Our sample of UDGs has a median stellar mass $\langle m_\star\rangle=2\times10^8\,\mathrm{M}_\odot$ and a median effective radius $\langle r_\mathrm{eff}\rangle=2.8$ kpc. We constrain the average mass of subhaloes within 30 kpc to $\log m_\mathrm{UDG}(r<30\,\mathrm{kpc})/\mathrm{M}_\odot\leq10.99$ at 95 per cent credibility, implying an effective virial mass $\log m_{200}/\mathrm{M}_\odot\leq11.80$, and a lower limit on the stellar mass fraction within 10 kpc of 1.0 per cent. Such mass is consistent with a simple extrapolation of the subhalo-to-stellar mass relation of typical satellite galaxies in massive clusters. However, our analysis is not sensitive to scatter about this mean mass; the possibility remains that extreme UDGs reside in haloes as massive as the Milky Way.

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.

Measuring and calibrating relations between cluster observables is critical for resource-limited studies. The mass-richness relation of clusters offers an observationally inexpensive way of estimating masses. Its calibration is essential for cluster and cosmological studies, especially for high-redshift clusters. Weak gravitational lensing magnification is a promising and complementary method to shear studies, that can be applied at higher redshifts. We employed the weak lensing magnification method to calibrate the mass-richness relation up to a redshift of 1.4. We used the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS) galaxy cluster candidates ($0.2<z<1.4$) and optical data from the Canada France Hawaii Telescope (CFHT) to test whether magnification can be effectively used to constrain the mass of high-redshift clusters. Lyman-Break Galaxies (LBGs) selected using the $u$-band dropout technique and their colours were used as a background sample of sources. LBG positions were cross-correlated with the centres of the sample of SpARCS clusters to estimate the magnification signal measured for cluster sub-samples, binned in both redshift and richness. We detected a weak lensing magnification signal for all bins at a detection significance of 2.6-5.5$\sigma$. In particular, the significance of the measurement for clusters with $z>1.0$ is 4.1$\sigma$; for the entire cluster sample we obtained an average M$_{200}$ of $1.28^{+0.23}_{-0.21}$ $\times 10^{14} \, \textrm{M}_{\odot}$. Our measurements demonstrated the feasibility of using weak lensing magnification as a viable tool for determining the average halo masses for samples of high redshift galaxy clusters. The results also established the success of using galaxy over-densities to select massive clusters at $z > 1$. Additional studies are necessary for further modelling of the various systematic effects we discussed.

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

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

In recent years, multiple studies have reported substantial populations of large, low-surface-brightness galaxies in local galaxy clusters. Various theories that aim to explain the presence of such ultra-diffuse galaxies (UDGs) have since been proposed. A key question that will help to differentiate between models is whether UDGs have counterparts in lower-mass host haloes, and what their abundance as a function of halo mass is. In this study we extend our previous study of UDGs in galaxy clusters to galaxy groups. We measure the abundance of UDGs in 325 spectroscopically-selected groups from the Galaxy And Mass Assembly (GAMA) survey. We make use of the overlapping imaging from the ESO Kilo-Degree Survey (KiDS), from which we can identify galaxies with mean surface brightnesses within their effective radii down to ~25.5 mag arcsec$^{-2}$ in the r-band. We are able to measure a significant overdensity of UDGs (with sizes r_eff > 1.5 kpc) in galaxy groups down to M200=10^12 Msun, a regime where approximately only 1 in 10 groups contains a UDG that we can detect. We combine measurements of the abundance of UDGs in haloes that cover three orders of magnitude in halo mass, finding that their numbers scale quite steeply with halo mass; N_UDG (R<R200) $\propto$ M200^(1.11+/-0.07). To better interpret this, we also measure the mass-richness relation for brighter galaxies down to $M^*_r$+2.5 in the same GAMA groups, and find a much shallower relation of N_Bright (R<R200) $\propto$ M200^(0.78+/-0.05). This shows that UDGs are relatively more abundant, compared to bright galaxies, in massive clusters than in groups. We discuss implications, but whether this difference is related to a higher destruction rate of UDGs in groups, or whether massive haloes have a positive effect on their formation, is not yet clear.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.

We present ALMA CO (2-1) detections in 11 gas-rich cluster galaxies at z~1.6, constituting the largest sample of molecular gas measurements in z>1.5 clusters to date. The observations span three galaxy clusters, derived from the Spitzer Adaptation of the Red-sequence Cluster Survey. We augment the >5sigma detections of the CO (2-1) fluxes with multi-band photometry, yielding stellar masses and infrared-derived star formation rates, to place some of the first constraints on molecular gas properties in z~1.6 cluster environments. We measure sizable gas reservoirs of 0.5-2x10^11 solar masses in these objects, with high gas fractions and long depletion timescales, averaging 62% and 1.4 Gyr, respectively. We compare our cluster galaxies to the scaling relations of the coeval field, in the context of how gas fractions and depletion timescales vary with respect to the star-forming main sequence. We find that our cluster galaxies lie systematically off the field scaling relations at z=1.6 toward enhanced gas fractions, at a level of ~4sigma, but have consistent depletion timescales. Exploiting CO detections in lower-redshift clusters from the literature, we investigate the evolution of the gas fraction in cluster galaxies, finding it to mimic the strong rise with redshift in the field. We emphasize the utility of detecting abundant gas-rich galaxies in high-redshift clusters, deeming them as crucial laboratories for future statistical studies.

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 report the detection of CO(2-1) emission coincident with the brightest cluster galaxy (BCG) of the high-redshift galaxy cluster SpARCS1049+56, with the Redshift Search Receiver (RSR) on the Large Millimetre Telescope (LMT). We confirm a spectroscopic redshift for the gas of z = 1.7091+/-0.0004, which is consistent with the systemic redshift of the cluster galaxies of z = 1.709. The line is well-fit by a single component Gaussian with a RSR resolution-corrected FWHM of 569+/-63 km/s. We see no evidence for multiple velocity components in the gas, as might be expected from the multiple image components seen in near-infrared imaging with the Hubble Space Telescope. We measure the integrated flux of the line to be 3.6+/-0.3 Jy km/s and, using alpha_CO = 0.8 Msun (K km s^-1 pc^2)^-1 we estimate a total molecular gas mass of 1.1+/-0.1x10^11 Msun and a M_H2/M_star ~ 0.4. This is the largest gas reservoir detected in a BCG above z > 1 to date. Given the infrared-estimated star formation rate of 860+/-130 Msun/yr, this corresponds to a gas depletion timescale of ~0.1Gyr. We discuss several possible mechanisms for depositing such a large gas reservoir to the cluster center -- e.g., a cooling flow, a major galaxy-galaxy merger or the stripping of gas from several galaxies -- but conclude that these LMT data are not sufficient to differentiate between them.

We compile a sample of spectroscopically- and photometrically-selected cluster galaxies from four high-redshift galaxy clusters ($1.59 < z < 1.71$) from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS), and a comparison field sample selected from the UKIDSS Deep Survey. Using near-infrared imaging from the \textit{Hubble Space Telescope} we classify potential mergers involving massive ($M_* \geq 3\times 10^{10}\mathrm{M}_\odot$) cluster members by eye, based on morphological properties such as tidal distortions, double nuclei, and projected near neighbors within 20 kpc. With a catalogue of 23 spectroscopic and 32 photometric massive cluster members across the four clusters and 65 spectroscopic and 26 photometric comparable field galaxies, we find that after taking into account contamination from interlopers, $11.0 ^{+7.0}_{-5.6}\%$ of the cluster members are involved in potential mergers, compared to $24.7^{+5.3}_{-4.6}\%$ of the field galaxies. We see no evidence of merger enhancement in the central cluster environment with respect to the field, suggesting that galaxy-galaxy merging is not a stronger source of galaxy evolution in cluster environments compared to the field at these redshifts.

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 present the results of a Spitzer/Herschel infrared photometric analysis of the largest (716) and highest-redshift (z=1.8) sample of Brightest Cluster Galaxies (BCGs), those from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). Given the tension that exists between model predictions and recent observations of BCGs at z<2, we aim to uncover the dominant physical mechanism(s) guiding the stellar-mass buildup of this special class of galaxies, the most massive in the Universe uniquely residing at the centres of galaxy clusters. Through a comparison of their stacked, broadband, infrared spectral energy distributions (SEDs) to a variety of SED model templates in the literature, we identify the major sources of their infrared energy output, in multiple redshift bins between 0 < z < 1.8. We derive estimates of various BCG physical parameters from the stacked {\nu}L{\nu} SEDs, from which we infer a star-forming, as opposed to a 'red and dead' population of galaxies, producing tens to hundreds of solar masses per year down to z=0.5. This discovery challenges the accepted belief that BCGs should only passively evolve through a series of gas-poor, minor mergers since z~4 (De Lucia & Blaizot 2007), but agrees with the improved semi-analytic model of hierarchical structure formation of Tonini et al. (2012), which predicts star-forming BCGs throughout the epoch considered. We attribute the star formation inferred from the stacked infrared SEDs to both major and minor 'wet' (gas-rich) mergers, based on a lack of key signatures (to date) of the cluster cooling flows to which BCG star formation is typically attributed, as well as a number of observational and simulation-based studies that support this scenario.

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

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

We analyse the evolution of environmental quenching efficiency, the fraction of quenched cluster galaxies that would be star-forming if they were in the field, as a function of redshift in 14 spectroscopically confirmed galaxy clusters with 0.87 < z < 1.63 from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). The clusters are the richest in the survey at each redshift. Passive fractions rise from $42_{-13}^{+10}$\% at z ~ 1.6 to $80_{-9}^{+12}$\% at z ~ 1.3 and $88_{-3}^{+4}$\% at z < 1.1, outpacing the change in passive fraction in the field. Environmental quenching efficiency rises dramatically from $16_{-19}^{+15}$ at z ~ 1.6 to $62_{-15}^{+21}\% at z ~ 1.3 and $73_{-7}^{+8}$\% at z $\lesssim$ 1.1. This work is the first to show direct observational evidence for a rapid increase in the strength of environmental quenching in galaxy clusters at z ~ 1.5, where simulations show cluster-mass halos undergo non-linear collapse and virialisation.

We constrain the internal dynamics of a stack of 10 clusters from the GCLASS survey at 0.87<z<1.34. We determine the stack cluster mass profile M(r) using the MAMPOSSt algorithm of Mamon et al., the velocity anisotropy profile beta(r) from the inversion of the Jeans equation, and the pseudo-phase-space density profiles Q(r) and Qr(r), obtained from the ratio between the mass density profile and the third power of the (total and, respectively, radial) velocity dispersion profiles of cluster galaxies. Several M(r) models are statistically acceptable for the stack cluster (Burkert, Einasto, Hernquist, NFW). The total mass distribution has a concentration c=r200/r-2=4.0-0.6+1.0, in agreement with theoretical expectations, and is less concentrated than the cluster stellar-mass distribution. The stack cluster beta(r) is similar for passive and star-forming galaxies and indicates isotropic galaxy orbits near the cluster center and increasingly radially elongated with increasing cluster-centric distance. Q(r) and Qr(r) are almost power-law relations with slopes similar to those predicted from numerical simulations of dark matter halos. Combined with results obtained for lower-z clusters we determine the dynamical evolution of galaxy clusters, and compare it with theoretical predictions. We discuss possible physical mechanisms responsible for the differential evolution of total and stellar mass concentrations, and of passive and star-forming galaxy orbits [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 the stellar mass functions (SMFs) of passive and star-forming galaxies with a limiting mass of 10$^{10.1}$ M$_{\odot}$ in four spectroscopically confirmed Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) galaxy clusters at 1.37 $<$ z $<$ 1.63. The clusters have 113 spectroscopically confirmed members combined, with 8-45 confirmed members each. We construct $Ks$-band-selected photometric catalogs for each cluster with an average of 11 photometric bands ranging from $u$ to 8 $\mu$m. We compare our cluster galaxies to a field sample derived from a similar $Ks$-band-selected catalog in the UltraVISTA/COSMOS field. The SMFs resemble those of the field, but with signs of environmental quenching. We find that 30 $\pm$ 20\% of galaxies that would normally be forming stars in the field are quenched in the clusters. The environmental quenching efficiency shows little dependence on projected cluster-centric distance out to $\sim$ 4 Mpc, providing tentative evidence of pre-processing and/or galactic conformity in this redshift range. We also compile the available data on environmental quenching efficiencies from the literature, and find that the quenching efficiency in clusters and in groups appears to decline with increasing redshift in a manner consistent with previous results and expectations based on halo mass growth.

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.

Recent observations have highlighted a significant population of faint but large (r_eff>1.5 kpc) galaxies in the Coma cluster. The origin of these Ultra Diffuse Galaxies (UDGs) remains puzzling, as the interpretation of the observational results has been hindered by the subjective selection of UDGs, and the limited study of only the Coma (and some examples in the Virgo-) cluster. We extend the study of UDGs using 8 clusters in the redshift range 0.044<z<0.063 with deep g- and r-band imaging data taken with MegaCam at the CFHT. We describe an automatic selection pipeline for quantitative identification, tested for completeness using image simulations of these galaxies. We find that the abundance of the UDGs we can detect increases with cluster mass, reaching ~200 in typical haloes of M200~10^15 Msun. The cluster UDGs have colours consistent with the cluster red sequence, and have a steep size distribution that declines as n~r_eff^-3.4. Their radial distribution is significantly steeper than NFW in the outskirts, and is significantly shallower in the inner parts. They follow the same radial distribution as the more massive quiescent galaxies in the clusters, except within the core region of r<0.15XR200 (or <300 kpc). Within this region the number density of UDGs drops and is consistent with zero. These diffuse galaxies can only resist tidal forces down to this cluster-centric distance if they are highly centrally dark-matter dominated. The observation that the radial distribution of more compact dwarf galaxies (r_eff<1.0 kpc) with similar luminosities follows the same distribution as the UDGs, but exist down to a smaller distance of 100kpc from the cluster centres, indicates that they may have similarly massive sub-haloes as the UDGs. Although several scenarios can give rise to the UDG population, our results point to differences in the formation history as the most plausible explanation.

The Large Early Galaxy Census (LEGA-C) is a Public Spectroscopic Survey of $\sim3200$ $K$-band selected galaxies at redshifts $z=0.6-1.0$ with stellar masses M_star > 1e10M_sun, conducted with VIMOS on ESO's Very Large Telescope. The survey is embedded in the COSMOS field ($R.A. = 10h00$; $Dec.=+2\deg$). The 20-hour long integrations produce high-$S/N$ continuum spectra that reveal ages, metallicities and velocity dispersions of the stellar populations. LEGA-C's unique combination of sample size and depth will enable us for the first time to map the stellar content at large look-back time, across galaxies of different types and star-formation activity. Observations started in December 2014 and are planned to be completed by mid 2018, with early data releases of the spectra and value-added products. In this paper we present the science case, the observing strategy, an overview of the data reduction process and data products, and a first look at the relationship between galaxy structure and spectral properties, as it existed 7 Gyr ago.

We present a new technique for wide and shallow observations using the near-infrared channel of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). Wide-field near-IR surveys with HST are generally inefficient, as guide star acquisitions make it impractical to observe more than one pointing per orbit. This limitation can be circumvented by guiding with gyros alone, which is possible as long as the telescope has three functional gyros. The method presented here allows us to observe mosaics of eight independent WFC3-IR pointings in a single orbit by utilizing the fact that HST drifts by only a very small amount in the 25 seconds between non-destructive reads of unguided exposures. By shifting the reads and treating them as independent exposures the full resolution of WFC3 can be restored. We use this "drift and shift" (DASH) method in the Cycle 23 COSMOS-DASH program, which will obtain 456 WFC3 $H_{160}$ pointings in 57 orbits, covering an area of 0.6 degree$^2$ in the COSMOS field down to $H_{160} = 25$. When completed, the program will more than triple the area of extra-galactic survey fields covered by near-IR imaging at HST resolution. We demonstrate the viability of the method with the first four orbits (32 pointings) of this program. We show that the resolution of the WFC3 camera is preserved, and that structural parameters of galaxies are consistent with those measured in guided observations.

Measuring stellar velocity dispersions of quiescent galaxies beyond $z\sim2$ is observationally challenging. Such measurements require near-infrared spectra with a continuum detection of at least moderate signal-to-noise, often necessitating long integrations. In this paper, we present deep X-Shooter spectroscopy of one of only two known gravitationally-lensed massive quiescent galaxies at $z>2$. This galaxy is quadruply imaged, with the brightest images magnified by a factor of $\sim5$. The total exposure time of our data is 9.8 hours on-source; however the magnification, and the slit placement encompassing 2 images, provides a total equivalent exposure time of 215 hours. From this deep spectrum we measure a redshift ($z_{\mathrm{spec}}=2.756\pm0.001$), making this one of the highest redshift quiescent galaxies that is spectroscopically confirmed. We simultaneously fit both the spectroscopic and photometric data to determine stellar population parameters and conclude this galaxy is relatively young, intermediate-mass, consistent with low dust content, and has quenched only relatively recently. This recent quenching is confirmed by strong Balmer absorption, particularly $H\delta$. Remarkably, this proves that at least some intermediate-mass galaxies have already quenched as early as $z\sim2.8$. Additionally, we have measured a velocity dispersion ($\sigma=187\pm43~\mathrm{km/s}$), making this the highest-redshift quiescent galaxy with a dispersion measurement. We confirm that this galaxy falls on the same mass fundamental plane (MFP) as galaxies at z=2.2, consistent with little to no evolution in the MFP up to z=2.8. Overall this galaxy is proof of existence of intermediate-mass quenched galaxies in the distant universe, and that lensing is a powerful tool for determining their properties with improved accuracy.

We present an analysis of galaxies in groups and clusters at $0.8<z<1.2$, from the GCLASS and GEEC2 spectroscopic surveys. We compute a "conversion fraction" $f_{\rm convert}$ that represents the fraction of galaxies that were prematurely quenched by their environment. For massive galaxies, $M_{\rm star}>10^{10.3}M_\odot$, we find $f_{\rm convert}\sim 0.4$ in the groups and $\sim 0.6$ in the clusters, similar to comparable measurements at $z=0$. This means the time between first accretion into a more massive halo and final star formation quenching is $t_p\sim 2$ Gyr. This is substantially longer than the estimated time required for a galaxy's star formation rate to become zero once it starts to decline, suggesting there is a long delay time during which little differential evolution occurs. In contrast with local observations we find evidence that this delay timescale may depend on stellar mass, with $t_p$ approaching $t_{\rm Hubble}$ for $M_{\rm star}\sim 10^{9.5}M_\odot$. The result suggests that the delay time must not only be much shorter than it is today, but may also depend on stellar mass in a way that is not consistent with a simple evolution in proportion to the dynamical time. Instead, we find the data are well-matched by a model in which the decline in star formation is due to "overconsumption", the exhaustion of a gas reservoir through star formation and expulsion via modest outflows in the absence of cosmological accretion. Dynamical gas removal processes, which are likely dominant in quenching newly accreted satellites today, may play only a secondary role at $z=1$.

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 a five-band Herschel study (100-500um) of three galaxy clusters at z~1.2 from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). With a sample of 120 spectroscopically-confirmed cluster members, we investigate the role of environment on galaxy properties utilizing the projected cluster phase space (line-of-sight velocity versus clustercentric radius), which probes the time-averaged galaxy density to which a galaxy has been exposed. We divide cluster galaxies into phase-space bins of (r/r200) x (v/sigma_v), tracing a sequence of accretion histories in phase space. Stacking optically star-forming cluster members on the Herschel maps, we measure average infrared star formation rates, and, for the first time in high-redshift galaxy clusters, dust temperatures for dynamically distinct galaxy populations---namely, recent infalls and those that were accreted onto the cluster at an earlier epoch. Proceeding from the infalling to virialized (central) regions of phase space, we find a steady decrease in the specific star formation rate and increase in the stellar age of star-forming cluster galaxies. We perform a probability analysis to investigate all acceptable infrared spectral energy distributions within the full parameter space and measure a ~4 sigma drop in the average dust temperature of cluster galaxies in an intermediate phase-space bin, compared to an otherwise flat trend with phase space. We suggest one plausible quenching mechanism which may be consistent with these trends, invoking ram-pressure stripping of the warmer dust for galaxies within this intermediate accretion phase.

We study the slope, intercept, and scatter of the color-magnitude and color-mass relations for a sample of ten infrared red-sequence-selected clusters at z ~ 1. The quiescent galaxies in these clusters formed the bulk of their stars above z ~ 3 with an age spread {\Delta}t ~ 1 Gyr. We compare UVJ color-color and spectroscopic-based galaxy selection techniques, and find a 15% difference in the galaxy populations classified as quiescent by these methods. We compare the color-magnitude relations from our red-sequence selected sample with X-ray- and photometric- redshift-selected cluster samples of similar mass and redshift. Within uncertainties, we are unable to detect any difference in the ages and star formation histories of quiescent cluster members in clusters selected by different methods, suggesting that the dominant quenching mechanism is insensitive to cluster baryon partitioning at z ~ 1.

We present the results of a MIPS-24um study of the Brightest Cluster Galaxies (BCGs) of 535 high-redshift galaxy clusters. The clusters are drawn from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS), which effectively provides a sample selected on total stellar mass, over 0.2 < z < 1.8 within the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Survey fields. 20%, or 106 clusters have spectroscopically confirmed redshifts, and the rest have redshifts estimated from the color of their red sequence. A comparison with the public SWIRE images detects 125 individual BCGs at 24um > 100uJy, or 23%. The luminosity-limited detection rate of BCGs in similar richness clusters (Ngal> 12) increases rapidly with redshift. Above z ~ 1, an average of ~20\% of the sample have 24um-inferred infrared luminosities of LIR > 10^12 Lsun, while the fraction below z ~ 1 exhibiting such luminosities is < 1 \%. The Spitzer-IRAC colors indicate the bulk of the 24um-detected population is predominantly powered by star formation, with only 7/125 galaxies lying within the color region inhabited by Active Galactic Nuclei (AGN). Simple arguments limit the star-formation activity to several hundred million years and this may therefore be indicative of the timescale for AGN feedback to halt the star formation. Below redshift z ~ 1 there is not enough star formation to significantly contribute to the overall stellar mass of the BCG population, and therefore BCG growth is likely dominated by dry-mergers. Above z~ 1, however, the inferred star formation would double the stellar mass of the BCGs and is comparable to the mass assembly predicted by simulations through dry mergers. We cannot yet constrain the process driving the star formation for the overall sample, though a single object studied in detail is consistent with a gas-rich merger.

We have discovered an optically rich galaxy cluster at z=1.7089 with star formation occurring in close proximity to the central galaxy. The system, SpARCS104922.6+564032.5, was detected within the Spitzer Adaptation of the red-sequence Cluster Survey, (SpARCS), and confirmed through Keck-MOSFIRE spectroscopy. The rest-frame optical richness of Ngal(500kpc) = 30+/-8 implies a total halo mass, within 500kpc, of ~3.8+/-1.2 x 10^14 Msun, comparable to other clusters at or above this redshift. There is a wealth of ancillary data available, including Canada-France-Hawaii Telescope optical, UKIRT-K, Spitzer-IRAC/MIPS, and Herschel-SPIRE. This work adds submillimeter imaging with the SCUBA2 camera on the James Clerk Maxwell Telescope and near-infrared imaging with the Hubble Space Telescope (HST). The mid/far-infrared (M/FIR) data detect an Ultra-luminous Infrared Galaxy spatially coincident with the central galaxy, with LIR = 6.2+/-0.9 x 10^12 Lsun. The detection of polycyclic aromatic hydrocarbons (PAHs) at z=1.7 in a Spitzer-IRS spectrum of the source implies the FIR luminosity is dominated by star formation (an Active Galactic Nucleus contribution of 20%) with a rate of ~860+/-30 Msun/yr. The optical source corresponding to the IR emission is likely a chain of of > 10 individual clumps arranged as "beads on a string" over a linear scale of 66 kpc. Its morphology and proximity to the Brightest Cluster Galaxy imply a gas-rich interaction at the center of the cluster triggered the star formation. This system indicates that wet mergers may be an important process in forming the stellar mass of BCGs at early times.

Finding a sample of the most massive clusters with redshifts $z>0.6$ can provide an interesting consistency check of the $\Lambda$ cold dark matter ($\Lambda$CDM) model. Here we present results from our search for clusters with $0.6\lesssim z\lesssim1.0$ where the initial candidates were selected by cross-correlating the RASS faint and bright source catalogues with red galaxies from the Sloan Digital Sky Survey DR8. Our survey thus covers $\approx10,000\,\rm{deg^2}$, much larger than previous studies of this kind. Deeper follow-up observations in three bands using the William Herschel Telescope and the Large Binocular Telescope were performed to confirm the candidates, resulting in a sample of 44 clusters for which we present richnesses and red sequence redshifts, as well as spectroscopic redshifts for a subset. At least two of the clusters in our sample are comparable in richness to RCS2-$J$232727.7$-$020437, one of the richest systems discovered to date. We also obtained new observations with the Combined Array for Research in Millimeter Astronomy for a subsample of 21 clusters. For 11 of those we detect the Sunyaev-Zel'dovich effect signature. The Sunyaev-Zel'dovich signal allows us to estimate $M_{200}$ and check for tension with the cosmological standard model. We find no tension between our cluster masses and the $\Lambda$CDM model.

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.

We study the radial number density and stellar mass density distributions of satellite galaxies in a sample of 60 massive clusters at 0.04<z<0.26 selected from the Multi-Epoch Nearby Cluster Survey (MENeaCS) and the Canadian Cluster Comparison Project (CCCP). In addition to ~10,000 spectroscopically confirmed member galaxies, we use deep ugri-band imaging to estimate photometric redshifts and stellar masses, and then statistically subtract fore-, and background sources using data from the COSMOS survey. We measure the galaxy number density and stellar mass density distributions in logarithmically spaced bins over 2 orders of magnitude in radial distance from the BCGs. For projected distances in the range 0.1<R/R200<2.0, we find that the stellar mass distribution is well-described by an NFW profile with a concentration of c=2.03+/-0.20. However, at smaller radii we measure a significant excess in the stellar mass in satellite galaxies of about $10^{11}$ Msun per cluster, compared to these NFW profiles. We do obtain good fits to generalized NFW profiles with free inner slopes, and to Einasto profiles. To examine how clusters assemble their stellar mass component over cosmic time, we compare this local sample to the GCLASS cluster sample at z~1, which represents the approximate progenitor sample of the low-z clusters. This allows for a direct comparison, which suggests that the central parts (R<0.4 Mpc) of the stellar mass distributions of satellites in local galaxy clusters are already in place at z~1, and contain sufficient excess material for further BCG growth. Evolving towards z=0, clusters appear to assemble their stellar mass primarily onto the outskirts, making them grow in an inside-out fashion.

Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogs that include new deep near-infrared observations. Nonetheless we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1-b=0.76+-0.05(stat)}+-0.06(syst)}, which does not resolve the tension with the measurements from the primary cosmic microwave background.

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.

In April 2013, a workshop entitled "What Regulates Galaxy Evolution" was held at the Lorentz Center. The aim of the workshop was to bring together the observational and theoretical community working on galaxy evolution, and to discuss in depth of the current problems in the subject, as well as to review the most recent observational constraints. A total of 42 astrophysicists attended the workshop. A significant fraction of the time was devoted to identifying the most interesting "open questions" in the field, and to discuss how progress can be made. This review discusses the four questions (one for each day of the workshop) that, in our opinion, were the focus of the most intense debate. We present each question in its context, and close with a discussion of what future directions should be pursued in order to make progress on these problems.

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.

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 investigate the velocity vs. position phase space of z ~ 1 cluster galaxies using a set of 424 spectroscopic redshifts in 9 clusters drawn from the GCLASS survey. Dividing the galaxy population into three categories: quiescent, star-forming, and poststarburst, we find that these populations have distinct distributions in phase space. Most striking are the poststarburst galaxies, which are commonly found at small clustercentric radii with high clustercentric velocities, and appear to trace a coherent ``ring" in phase space. Using several zoom simulations of clusters we show that the coherent distribution of the poststarbursts can be reasonably well-reproduced using a simple quenching scenario. Specifically, the phase space is best reproduced if satellite quenching occurs on a rapid timescale (0.1 < tau_{Q} < 0.5 Gyr) after galaxies make their first passage of R ~ 0.5R_{200}, a process that takes a total time of ~ 1 Gyr after first infall. We compare this quenching timescale to the timescale implied by the stellar populations of the poststarburst galaxies and find that the poststarburst spectra are well-fit by a rapid quenching (tau_{Q} = 0.4^{+0.3}_{-0.4} Gyr) of a typical star-forming galaxy. The similarity between the quenching timescales derived from these independent indicators is a strong consistency check of the quenching model. Given that the model implies satellite quenching is rapid, and occurs well within R_{200}, this would suggest that ram-pressure stripping of either the hot or cold gas component of galaxies are the most plausible candidates for the physical mechanism. The high cold gas consumption rates at z ~ 1 make it difficult to determine if hot or cold gas stripping is dominant; however, measurements of the redshift evolution of the satellite quenching timescale and location may be capable of distinguishing between the two.

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

We study the stellar mass content of massive haloes in the redshift range 0.86<z<1.34, by measuring: (1) The stellar mass in the central galaxy versus total dynamical halo mass. (2) The total stellar mass (including satellites) versus total halo mass. (3) The radial stellar mass and number density profiles for the ensemble halo. We use a Ks-band selected catalogue for the 10 clusters in the Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS), with photometric redshifts and stellar masses measured from 11-band SED fitting. Combining the photometric catalogues with the deep spectroscopic component of GCLASS, we correct the cluster galaxy sample for interlopers. We also perform a dynamical analysis of the cluster galaxies to estimate the halo mass M200 for each cluster based on a measurement of its velocity dispersion. (1) We find that the central galaxy stellar mass fraction decreases with total halo mass, and that this is in reasonable quantitative agreement with measurements from abundance matching studies at z~1. (2) The total stellar mass fractions of these systems decrease with halo mass, indicating that lower mass systems are more efficient at transforming baryons into stars. We find the total stellar mass to be a good proxy for total halo mass, with a small intrinsic scatter. When we compare these results with literature measurements, we find that the stellar mass fraction at fixed halo mass shows no significant evolution in the range 0<z<1. (3) We measure a relatively high NFW concentration parameter c ~ 7 for the stellar mass distribution in these clusters. A simple model shows that the stellar mass content of GCLASS can evolve in typical distributions observed at lower redshifts if the clusters primarily accrete stellar mass onto the outskirts. [Abridged]

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.

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

The Spitzer-South Pole Telescope Deep Field (SSDF) is a wide-area survey using Spitzer's Infrared Array Camera (IRAC) to cover 94 square degrees of extragalactic sky, making it the largest IRAC survey completed to date outside the Milky Way midplane. The SSDF is centered at 23:30,-55:00, in a region that combines observations spanning a broad wavelength range from numerous facilities. These include millimeter imaging from the South Pole Telescope, far-infrared observations from Herschel/SPIRE, X-ray observations from the XMM XXL survey, near-infrared observations from the VISTA Hemisphere Survey, and radio-wavelength imaging from the Australia Telescope Compact Array, in a panchromatic project designed to address major outstanding questions surrounding galaxy clusters and the baryon budget. Here we describe the Spitzer/IRAC observations of the SSDF, including the survey design, observations, processing, source extraction, and publicly available data products. In particular, we present two band-merged catalogs, one for each of the two warm IRAC selection bands. They contain roughly 5.5 and 3.7 million distinct sources, the vast majority of which are galaxies, down to the SSDF 5-sigma sensitivity limits of 19.0 and 18.2 Vega mag (7.0 and 9.4 microJy) at 3.6 and 4.5 microns, respectively.

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

Recent independent results from numerical simulations and observations have shown that brightest cluster galaxies (BCGs) have increased their stellar mass by a factor of almost two between z~0.9 and z~0.2. The numerical simulations further suggest that more than half this mass is accreted through major mergers. Using a sample of 18 distant galaxy clusters with over 600 spectroscopically confirmed cluster members between them, we search for observational evidence that major mergers do play a significant role. We find a major merger rate of 0.38 +/- 0.14 mergers per Gyr at z~1. While the uncertainties, which stem from the small size of our sample, are relatively large, our rate is consistent with the results that are derived from numerical simulations. If we assume that this rate continues to the present day and that half of the mass of the companion is accreted onto the BCG during these mergers, then we find that this rate can explain the growth in the stellar mass of the BCGs that is observed and predicted by simulations. Major mergers therefore appear to be playing an important role, perhaps even the dominant one, in the build up of stellar mass in these extraordinary galaxies.

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 present an IR study of high-redshift galaxy clusters with the MIPS camera on the Spitzer Space Telescope. Employing a sample of 42 clusters from the RCS-1 over the redshift range 0.3 < z < 1.0 and spanning an approximate range in mass of 10^{14-15} Msun, we show the number of IR-luminous galaxies in clusters above a fixed IR luminosity of 2x10^{11} Msun per unit cluster mass evolves as (1+z)^{5.1+/-1.9}. These results assume a single star forming galaxy template; the presence of AGN, and an evolution in their relative contribution to the mid-IR galaxy emission, will alter the overall number counts per cluster and their rate of evolution. We infer the total SFR per unit cluster mass and find T_SFR/M_c ~ (1+z)^{5.4+/-1.9}. This evolution can be attributed entirely to the change in the in-falling field galaxy population. The T_SFR/M_c (binned over all redshift) decreases with increasing cluster mass with a slope (T_SFR/M_c ~ M_c^{-1.5+/-0.4}) consistent with the dependence of the stellar-to-total mass per unit cluster mass seen locally. The inferred star formation seen here could produce ~5-10% of the total stellar mass in massive clusters at z = 0. Finally, we show a clear decrease in the number of IR-bright galaxies per unit optical galaxy in the cluster cores, confirming star formation continues to avoid the highest density regions of the universe at z ~ 0.75 (the average redshift of the high-redshift clusters). While several previous studies appear to show enhanced star formation in high-redshift clusters relative to the field we note that these papers have not accounted for the overall increase in galaxy or dark matter density at the location of clusters. Once this is done, clusters at z ~ 0.75 have the same or less star formation per unit mass or galaxy as the field.

We present an infrared study of a z=0.872 cluster, SpARCS J161314+564930, with the primary aim of distinguishing the dynamical histories of spectroscopically confirmed star-forming members to assess the role of cluster environment. We utilize deep MIPS imaging and a mass-limited sample of 85 spectroscopic members to identify 16 24um-bright sources within the cluster, and measure their 24um star formation rates (SFRs) down to ~6 Msolar/year. Based on their line-of-sight velocities and stellar ages, MIPS cluster members appear to be an infalling population that was recently accreted from the field with minimal environmental dependency on their star formation. However, we identify a double-sequenced distribution of star-forming galaxies amongst the members, with one branch exhibiting declining specific SFRs with mass. The members along this sub-main sequence contain spectral features suggestive of passive galaxies. Using caustic diagrams, we kinematically identify these galaxies as a virialized and/or backsplash population. Moreover, we find a mix of dynamical histories at all projected radii, indicating that standard definitions of environment (i.e., radius and density) are contaminated with recently accreted interlopers, which could contribute to a lack of environmental trends for star-forming galaxies. A cleaner narrative of their dynamical past begins to unfold when using a proxy for accretion histories through profiles of constant (r/r_200)x(Delta v/sigma_v); galaxies accreted at earlier times possess lower values of (r/r_200)x(Delta v/sigma_v) with minimal contamination from the distinct infalling population. Therefore, adopting a time-averaged definition for density (as traced by accretion histories) rather than an instantaneous density yields a depressed specific SFR within the dynamical cluster core.

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

We present a new two-color algorithm, the "Stellar Bump Sequence" (SBS), that is optimized for robustly identifying candidate high-redshift galaxy clusters in combined wide-field optical and mid-infrared (MIR) data. The SBS algorithm is a fusion of the well-tested cluster red-sequence method of Gladders & Yee (2000) with the MIR 3.6um - 4.5um cluster detection method developed by Papovich (2008). As with the cluster red-sequence method, the SBS identifies candidate overdensities within 3.6um - 4.5um color slices, which are the equivalent of a rest-frame 1.6um stellar bump "red-sequence". In addition to employing the MIR colors of galaxies, the SBS algorithm incorporates an optical/MIR (z' - 3.6um) color cut. This cut effectively eliminates foreground 0.2 < z < 0.4 galaxies which have 3.6um - 4.5um colors that are similarly red as z > 1.0 galaxies and add noise when searching for high-redshift galaxy overdensities. We demonstrate using the z ~ 1 GCLASS cluster sample that similar to the red sequence, the stellar bump sequence appears to be a ubiquitous feature of high-redshift clusters, and that within that sample the color of the stellar bump sequence increases monotonically with redshift and provides photometric redshifts accurate to dz = 0.05. We apply the SBS method in the XMM-LSS SWIRE field and show that it robustly recovers the majority of confirmed optical, MIR, and X-ray-selected clusters at z > 1.0 in that field. Lastly, we present confirmation of SpARCS J022427-032354 at z = 1.63, a new cluster detected with the method and confirmed with 12 high-confidence spectroscopic redshifts obtained using FORS2 on the VLT. We conclude with a discussion of future prospects for using the algorithm.

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.

Using new and published data, we construct a sample of 160 brightest cluster galaxies (BCGs) spanning the redshift interval 0.03 < z < 1.63. We use this sample, which covers 70% of the history of the universe, to measure the growth in the stellar mass of BCGs after correcting for the correlation between the stellar mass of the BCG and the mass of the cluster in which it lives. We find that the stellar mass of BCGs increase by a factor of 1.8 between z=0.9 and z=0.2. Compared to earlier works, our result is closer to the predictions of semi-analytic models. However, BCGs at z=0.9, relative to BCGs at z=0.2, are still a factor of 1.5 more massive than the predictions of these models. Star formation rates in BCGs at z~1 are generally to low to result in significant amounts of mass. Instead, it is likely that most of the mass build up occurs through mainly dry mergers in which perhaps half of the mass is lost to the intra-cluster medium of the cluster.

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 evaluate the effects of environment and stellar mass on galaxy properties at 0.85 < z < 1.20 using a 3.6um-selected spectroscopic sample of 797 cluster and field galaxies drawn from the GCLASS survey. We confirm that for galaxies with LogM* > 9.3 the well-known correlations between environment and properties such as star-forming fraction (f_SF), SFR, SSFR, D(4000), and color are already in place at z ~ 1. We separate the effects of environment and stellar mass on galaxies by comparing the properties of star-forming and quiescent galaxies at fixed environment, and fixed stellar mass. The SSFR of star-forming galaxies at fixed environment is correlated with stellar mass; however, at fixed stellar mass it is independent of environment. The same trend exists for the D(4000) measures of both the star-forming and quiescent galaxies and shows that their properties are determined primarily by their stellar mass, not by their environment. Instead, it appears that environment's primary role is to control the fraction of star-forming galaxies. Using the spectra we identify candidate poststarburst galaxies and find that those with 9.3 < LogM* < 10.7 are 3.1 +/- 1.1 times more common in high-density regions compared to low-density regions. The clear association of poststarbursts with high-density regions as well as the lack of a correlation between the SSFRs and D(4000)s of star-forming galaxies with their environment suggests that at z ~ 1 the environmental-quenching timescale must be rapid. Lastly, we construct a simple quenching model which demonstrates that the lack of a correlation between the D(4000) of quiescent galaxies and their environment results naturally if self quenching dominates over environmental quenching at z > 1, or if the evolution of the self-quenching rate mirrors the evolution of the environmental-quenching rate at z > 1, regardless of which dominates.

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

We 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 introduce a novel method to measure the masses of galaxy clusters at high redshift selected from optical and IR Spitzer data via the red-sequence technique. Lyman-break galaxies are used as a well understood, high-redshift background sample allowing mass measurements of lenses at unprecedented high redshifts using weak lensing magnification. By stacking a significant number of clusters at different redshifts with average masses of ~1-3x10^14M_sun, as estimated from their richness, we can calibrate the normalisation of the mass-richness relation. With the current data set (area: 6 deg^2) we detect a magnification signal at the >3-sigma level. There is good agreement between the masses estimated from the richness of the clusters and the average masses estimated from magnification, albeit with large uncertainties. We perform tests that suggest the absence of strong systematic effects and support the robustness of the measurement. This method - when applied to larger data sets in the future - will yield an accurate calibration of the mass-observable relations at z>~1 which will represent an invaluable input for cosmological studies using the galaxy cluster mass function and astrophysical studies of cluster formation. Furthermore this method will probably be the least expensive way to measure masses of large numbers of z>1 clusters detected in future IR-imaging surveys.

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.

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

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

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a z'-passband imaging survey of the 50 deg^2 Spitzer SWIRE Legacy fields, designed with the primary aim of creating the first large, homogeneously selected sample of massive clusters at z>1. SpARCS uses an infrared adaptation of the two-filter cluster red-sequence technique. In this paper we report Keck/LRIS spectroscopic confirmation of two new exceptionally rich galaxy clusters, SpARCS J161315+564930 at z=0.871+/-0.002, with 14 high-confidence members and a rest-frame velocity dispersion of sigma_v= 1230+/-320 km s^-1, and SpARCS J161641+554513 at z=1.161+/-0.003, with seven high-confidence members (including one AGN) and a rest-frame velocity dispersion of sigma_v=950+/-330 km s^-1. We also report confirmation of a third new system, SpARCS J161037+552417 at z=1.210+/-0.002, with seven high-confidence members and a rest-frame velocity dispersion of sigma_v=410+/-300 km s^-1. These three new spectroscopically confirmed clusters further demonstrate the efficiency and effectiveness of two-filter imaging for detecting bona fide galaxy clusters at high redshift. We conclude by demonstrating that prospects are good for the current generation of surveys aiming to estimate cluster redshifts and masses at z>~1 directly from optical-infrared imaging.

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

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

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a z'-passband imaging survey, consisting of deep (z' ~ 24 AB) observations made from both hemispheres using the CFHT 3.6m and CTIO 4m telescopes. The survey was designed with the primary aim of detecting galaxy clusters at z >~ 1. In tandem with pre-existing 3.6um observations from the Spitzer Space Telescope SWIRE Legacy Survey, SpARCS detects clusters using an infrared adaptation of the two-filter red-sequence cluster technique. The total effective area of the SpARCS cluster survey is 41.9 deg^2. In this paper, we provide an overview of the 13.6 deg^2 Southern CTIO/MOSAICII observations. The 28.3 deg^2 Northern CFHT/MegaCam observations are summarized in a companion paper by Muzzin et al. (2008). In this paper, we also report spectroscopic confirmation of SpARCS J003550-431224, a very rich galaxy cluster at z = 1.335, discovered in the ELAIS-S1 field. To date, this is the highest spectroscopically confirmed redshift for a galaxy cluster discovered using the red-sequence technique. Based on nine confirmed members, SpARCS J003550-431224 has a preliminary velocity dispersion of 1050 +/- 230 km/s. With its proven capability for efficient cluster detection, SpARCS is a demonstration that we have entered an era of large, homogeneously-selected z > 1 cluster surveys.

The Spitzer Adaptation of the Red-sequence Cluster Survey (SpARCS) is a deep z'-band imaging survey covering the Spitzer SWIRE Legacy fields designed to create the first large homogeneously-selected sample of massive clusters at z > 1 using an infrared adaptation of the cluster red-sequence method. We present an overview of the northern component of the survey which has been observed with CFHT/MegaCam and covers 28.3 deg^2. The southern component of the survey was observed with CTIO/MOSAICII, covers 13.6 deg^2, and is summarized in a companion paper by Wilson et al. (2008). We also present spectroscopic confirmation of two rich cluster candidates at z ~ 1.2. Based on Nod-and-Shuffle spectroscopy from GMOS-N on Gemini there are 17 and 28 confirmed cluster members in SpARCS J163435+402151 and SpARCS J163852+403843 which have spectroscopic redshifts of 1.1798 and 1.1963, respectively. The clusters have velocity dispersions of 490 +/- 140 km/s and 650 +/- 160 km/s, respectively which imply masses (M200) of (1.0 +/- 0.9) x 10^{14} M_{solar} and (2.4 +/- 1.8) x 10^{14} M_{solar}. Confirmation of these candidates as bona fide massive clusters demonstrates that two-filter imaging is an effective, yet observationally efficient, method for selecting clusters at z > 1.

We present a catalogue of 99 candidate clusters and groups of galaxies in the redshift range 0.1 < z_{phot} < 1.3 discovered in the Spitzer First Look Survey (FLS). The clusters are selected by their R - 3.6um galaxy color-magnitude relation using the cluster red sequence algorithm. Using this cluster sample we compute the 3.6um, 4.5um, 5.8um, & 8.0um cluster luminosity functions (LFs). Similar to previous studies, we find that for the bands that trace stellar mass at these redshifts (3.6um, 4.5um) the evolution in M* is consistent with a passively evolving population of galaxies with a high formation redshift (z_{f} > 1.5). Using the 3.6um LF as a proxy for stellar luminosity we remove this component from the MIR (5.8um & 8.0um) cluster LFs and measure the LF of dusty star formation/AGN in clusters. We find that at z < 0.4 the bright end of the cluster 8.0um LF is well-described by a composite population of quiescent galaxies and regular star forming galaxies with a mix consistent with typical cluster blue fractions; however, at z > 0.4, an additional population of dusty starburst galaxies is required to properly model the 8.0um LFs. Comparison to field studies at similar redshifts shows a strong differential evolution in the field and cluster 8.0um LFs with redshift. At z ~ 0.65 8.0um-detected galaxies are more abundant in clusters compared to the field, but thereafter the number of 8.0um sources in clusters declines with decreasing redshift and by z ~ 0.15, clusters are underdense relative to the field by a factor of ~5. The rapid differential evolution between the cluster and field LFs is qualitatively consistent with recent field galaxy studies that show the star formation rates of galaxies in high density environments are larger than those in low density environments at higher redshift.

Previous studies have indicated that many of the RR Lyrae variables in the LMC have properties similar to the ones in the Galactic globular cluster M3. Assuming that the M3 RR Lyrae variables follow the same relationships among period, temperature, amplitude and Fourier phase parameter phi31 as their LMC counterparts, we have used the M3 phi31-logP relation to identify the M3-like unevolved first overtone RR Lyrae variables in 16 fields near the LMC bar. The temperatures of these variables were calculated from the M3 logP-logTe relation so that the extinction could be derived for each star separately. Since blended stars have lower amplitudes for a given period, the period amplitude relation should be a useful tool for identifying which stars are affected by crowding. We find that the low amplitude stars are brighter. We remove them from the sample and derive an LMC distance modulus 18.49+/-0.11.

Using K-band imaging for 15 of the Canadian Network for Observational Cosmology (CNOC1) clusters we examine the near-infrared properties of moderate-redshift (0.19 < z < 0.55) galaxy clusters. We find that the number of K-band selected cluster galaxies within R_{500} (the Halo Occupation Number, HON) is well-correlated with the the cluster dynamical mass (M_{500}) and X-ray Temperature (T_{x}); however, the intrinsic scatter in these scaling relations is 37% and 46% respectively. Comparison with clusters in the local universe shows that the HON-M_{500} relation does not evolve significantly between z = 0 and z ~ 0.3. This suggests that if dark matter halos are disrupted or undergo significant tidal-stripping in high-density regions as seen in numerical simulations, the stellar mass within the halos is tightly bound, and not removed during the process. The total K-band cluster light (L_{200,K}) and K-band selected richness (parameterized by B_{gc,K}) are also correlated with both the cluster T_{x} and M_{200}. The total (intrinsic) scatter in the L_{200,K}-M_{200} and B_{gc,K}-M_{200} relations are 43%(31%) and 35%(18%) respectively and indicates that for massive clusters both L_{200,K} and B_{gc,K} can predict M_{200} with similar accuracy as T_{x}, L_{x} or optical richness (B_{gc}). Examination of the mass-to-light ratios of the clusters shows that similar to local clusters, the K-band mass-to-light ratio is an increasing function of halo mass. Using the K-band mass-to-light ratios of the clusters, we apply the Oort technique and find Omega_{m,0} = 0.22 pm 0.02, which agrees well with recent combined concordance cosmology parameters, but, similar to previous cluster studies, is on the low-density end of preferred values.

We present K-band imaging for 15 of the Canadian Network for Observational Cosmology (CNOC1) clusters. The cluster luminosity density and number density profiles can be described by NFW models with concentration parameters of c_{l} = 4.28 \pm 0.70 and c_{g} = 4.13 \pm 0.57 respectively. Comparing these to the dynamical mass analysis of the same clusters shows that the galaxy luminosity and number density profiles are similar to the dark matter profile, and are not less concentrated like in local clusters. The luminosity functions show that the evolution of K^{*} over the redshift range 0.2 < z < 0.5 is consistent with a scenario where the majority of stars in cluster galaxies form at high-redshift (z_{f} > 1.5) and evolve passively thereafter. The best-fit for the faint-end slope of the luminosity function is alpha = -0.84 pm 0.08, which indicates that it does not evolve between z = 0 and z = 0.3. Using Principal Component Analysis of the spectra we classify cluster galaxies as either star-forming/recently-star-forming (EM+BAL) or non-starforming (ELL) and compute their respective luminosity functions. The faint-end slope of the ELL luminosity function is much shallower than for the EM+BAL galaxies at z = 0.3, and suggests the number of faint ELL galaxies in clusters decreases by a factor of ~3 from z = 0 to z = 0.3. The redshift evolution of K^{*} for both EM+BAL and ELL types is consistent with a passively evolving stellar population formed at high-redshift. Passive evolution in both classes, as well as the total cluster luminosity function, demonstrates that the bulk of the stellar population in all bright cluster galaxies is formed at high-redshift and subsequent transformations in morphology/color/spectral-type have little effect on the total stellar mass.

As the densest galaxy environments in the universe, clusters are vital to our understanding of the role that environment plays in galaxy formation and evolution. Unfortunately, the evolution of high-redshift cluster galaxies is poorly understood because of the ``cluster desert'' that exists at 1 < z < 2. The SpARCS collaboration is currently carrying out a 1-passband (z') imaging survey which, when combined with the pre-existing 50 square degree 3.6 micron Spitzer SWIRE Legacy Survey data, will efficiently detect hundreds of clusters in the cluster desert using an infrared application of the well-proven cluster red-sequence technique. We have already tested this 1-color (z' - [3.6]) approach using a 6 square degree ``pilot patch'' and shown it to be extremely successful at detecting clusters at 1 < z < 2. The clusters discovered in this project will be the first large sample of ``nascent'' galaxy clusters which connect the star-forming proto-cluster regions at z > 2 to the quiescent population at z < 1. The existing seven-passband Spitzer data (3.6, 4.5, 5.8, 8.0, 24, 70, 160 micron) will allow us to make the first measurements of the evolution of the cluster red-sequence, IR luminosity function, and the mid-IR dust-obscured star-formation rate for 1 < z < 2 clusters.

We have measured the 3.6 micron cluster luminosity function (LF) using a sample of 123 galaxy clusters selected from the 4 degree^2 Spitzer First Look Survey (FLS). The clusters were selected on the basis of their R - 3.6 micron colors using the cluster red-sequence technique of Gladders & Yee (2000). The binned LFs are well-fit by a Schechter function at all redshifts. However, we note two interesting trends. Firstly, the evolution of Mstar with redshift is consistent with models that form galaxies in a single-burst of star-formation at z > 2.0, and evolve passively thereafter. Secondly, the faint-end slope of the LF appears to become shallower at higher redshift. We conclude that the most massive galaxies were formed in the cluster at high redshift (z > 2), while lower-mass galaxies have subsequently been accreted from the field at lower redshift (z ~ 0.5). These results are consistent with the ``downsizing'' picture seen in previous cluster studies using smaller samples.

We present an infrared adaptation of the Cluster Red-Sequence method. We modify the two filter technique of Gladders & Yee (2000) to identify clusters based on their R-[3.6] color. We apply the technique to the 4 degree^2 Spitzer First Look Survey and detect 123 clusters spanning the redshift range 0.09 < z < 1.4. Our results demonstrate that the Spitzer Space Telescope will play an important role in the discovery of large samples of high redshift galaxy clusters.

We present early results from the spectroscopic follow-up of a sample of candidate obscured AGN selected in the mid-infrared from the Spitzer First Look Survey. Our selection allows a direct comparison of the numbers of obscured and unobscured AGN at a given luminosity for the first time, and shows that the ratio of obscured to unobscured AGN at infrared luminosities corresponding to low luminosity quasars is ~1:1 at z~0.5. Most of our optically-faint candidate obscured AGN have the high-ionization, narrow-line spectra expected from type-2 AGN. A composite spectrum shows evidence for Balmer absorption lines, indicating recent star-formation activity in the host galaxies. There is tentative evidence for a decrease in the obscured AGN fraction with increasing AGN luminosity.

Fourier coefficents have been derived for the $V$ and $R$ light curves of 785 overtone RR Lyrae variables in 16 MACHO fields near the bar of the LMC. The $\phi_{31}$ and $R_{21}$ coefficients have been compared with those of the first overtone RR Lyrae variables in the Galactic globular clusters NGC 6441, M107, M5, M3, M2, $\omega$ Centauri and M68. The results indicate that many of the LMC variables have properties similar to the ones in M2, M3, M5 and the Oosterhoff type I variables in $\omega$ Cen, but they are different from the Oosterhoff type II variables in $\omega$ Cen. Equations derived from hydrodynamic pulsation models have been used to calculate the luminosity and temperature for the 330 bona fide first-overtone variables. The results indicate that they have $\log L$ in the range 1.6 to $1.8\lsun$ and $\log T_{eff}$ between 3.85 and 3.87. Based on these temperatures, a mean color excess $E(V-R) =0.08$ mag, equivalent to $E(B-V)=0.14$ mag, has been estimated for these 330 stars. The 80 M5-like variables (selected according to their location in the $\phi_{31}-\log P$ plot) are used to determine a LMC distance. After correcting for the effects of extinction and crowding, a mean apparent magnitude $<V_0>=18.99 \pm 0.02$ (statistical) $\pm 0.16$ (systematic) has been estimated for these 80 stars. Combining this with a mean absolute magnitude $M_V=0.56\pm 0.06$ for M5-like stars derived from Baade-Wesselink analyses, main sequence fitting, Fourier parameters and the trigonometric parallax of RR Lyrae, we derive an LMC distance modulus $\mu=18.43\pm 0.06$ (statistical) $\pm 0.16$ (systematic) mag. The large systematic error arises from the difficulties of correcting for interstellar extinction and for crowding.

Based on a search of the literature up to May 2001, the number of known variable stars in Galactic globular clusters is approximately 3000. Of these, more than 2200 have known periods and the majority (approximately 1800) are of the RR Lyrae type. In addition to the RR Lyrae population, there are approximately 100 eclipsing binaries, 120 SX Phe variables, 60 Cepheids (including population II Cepheids, anomalous Cepheids and RV Tauri) and 120 SR/red variables. The mean period of the fundamental mode RR Lyrae variables is 0.585, for the overtone variables it is 0.342 (0.349 for the first-overtone pulsators and 0.296 for the second-overtone pulsators) and approximately 30% are overtone pulsators. These numbers indicate that about 65% of RR Lyrae variables in Galactic globular clusters belong to Oosterhoff type I systems. The mean period of the RR Lyrae variables in the Oosterhoff type I clusters seems to be correlated with metal abundance in the sense that the periods are longer in the more metal poor clusters. Such a correlation does not exist for the Oosterhoff type II clusters. Most of the Cepheids are in clusters with blue horizontal branches.