Source author record

Adam Muzzin

Adam Muzzin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

78works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

78 published item(s)

preprint2026arXiv

JWST Reveals a Candidate Jellyfish Galaxy at z=1.156

We report the discovery of COSMOS2020-635829 as a candidate jellyfish galaxy undergoing ram pressure stripping in a (proto)cluster at $z > 1$. High-resolution imaging from the James Webb Space Telescope reveals a symmetric stellar disk coupled to a unilateral tail of star-forming knots to the south. Using Gemini GMOS IFU observations, we show that these extra-planar continuum sources are embedded within an ionized gas tail that is kinematically connected to the disk of COSMOS2020-635829. If confirmed, this represents the highest-redshift discovery of a ram pressure stripped ionized gas tail. The tail sources are characterized by extremely young stellar populations ($\lesssim 100\,\mathrm{Myr}$), have stellar masses of ${\sim}10^8\,\mathrm{M_\odot}$, and star formation rates of $0.1\text{--}1\,\mathrm{M_\odot\,yr^{-1}}$. This work shows that ram pressure stripping can potentially perturb group and cluster galaxies at $z > 1$ and may contribute to environmental quenching even near Cosmic Noon.

preprint2024arXiv

Modelling and Subtracting Diffuse Cluster Light in JWST Images: A Relation between the Spatial Distribution of Globular Clusters, Dwarf Galaxies, and Intracluster Light in the Lensing Cluster SMACS 0723

We present a methodology for modeling and removing light from cluster galaxies and intracluster light (ICL) from $James\ Webb\ Space\ Telescope$ ($JWST$) images of gravitational lensing clusters. We apply our method to Webb's First Deep Field the SMACS 0723 Early Release Observations and use the ICL subtracted images to select a sample of globular clusters (GCs) and dwarf galaxies within the cluster. We compare the spatial distributions of these two samples with our models of the galaxy and ICL light, finding significant similarity. Specifically we find that GCs trace the diffuse ICL, while dwarf galaxies are centrally concentrated near the cluster center We quantify the relationship between the surface density of compact sources and total cluster light, demonstrating a significant, tight correlation. We repeat our methodology and compare distributions of GCs with dark matter surface density and find a comparable result. Our findings suggest a common origin for GCs and diffuse ICL, with stripping from massive galaxies as they merge with the cluster being a plausible scenario.

preprint2022arXiv

3D-DASH: The Widest Near-Infrared Hubble Space Telescope Survey

The 3D-Drift And SHift (3D-DASH) program is a \textit{Hubble Space Telescope} WFC3 F160W imaging and G141 grism survey of the equatorial COSMOS field. 3D-DASH extends the legacy of HST near-infrared imaging and spectroscopy to degree-scale swaths of the sky, enabling the identification and study of distant galaxies ($z>2$) that are rare or in short-lived phases of galaxy evolution at rest-frame optical wavelengths. Furthermore, when combined with existing ACS/F814W imaging, the program facilitates spatially-resolved studies of the stellar populations and dust content of intermediate-redshift ($0.5<z<2$) galaxies. Here we present the reduced F160W imaging mosaic available to the community. Observed with the efficient DASH technique, the mosaic comprises 1256 individual WFC3 pointings, corresponding to an area of 1.35 deg$^2$ (1.43 deg$^2$ in 1912 when including archival data). The median $5σ$ point-source limit in $H_{160}$ is 24.74 mag. We also provide tools to determine the local point spread function (PSF), create cutouts, and explore the image at any location within the 3D-DASH footprint. 3D-DASH is the widest \textit{HST}/WFC3 imaging survey in the F160W filter to date, increasing the existing extragalactic survey area in the near-infrared at HST resolution by an order of magnitude.

preprint2022arXiv

An Assessment of the In-Situ Growth of the Intracluster Light in the High Redshift Galaxy Cluster SpARCS1049+56

The formation of the stellar mass within galaxy cluster cores is a poorly understood process. It features the complicated physics of cooling flows, AGN feedback, star formation and more. Here, we study the growth of the stellar mass in the vicinity of the Brightest Cluster Galaxy (BCG) in a z = 1.7 cluster, SpARCS1049+56. We synthesize a reanalysis of existing HST imaging, a previously published measurement of the star formation rate, and the results of new radio molecular gas spectroscopy. These analyses represent the past, present and future star formation respectively within this system. We show that a large amount of stellar mass -- between $(2.2 \pm 0.5) \times 10^{10} \: M_\odot$ and $(6.6 \pm 1.2) \times 10^{10}\: M_\odot$ depending on the data processing -- exists in a long and clumpy tail-like structure that lies roughly 12 kpc off the BCG. Spatially coincident with this stellar mass is a similarly massive reservoir ($(1.0 \pm 0.7) \times 10^{11} \: M_\odot$) of molecular gas that we suggest is the fuel for the immense star formation rate of $860 \pm 130 \: M_\odot$/yr, as measured by infrared observations. Hlavacek-Larrondo et al. 2021 surmised that massive, runaway cooling of the hot intracluster X-ray gas was feeding this star formation, a process that had not been observed before at high-redshift. We conclude, based on the amount of fuel and current stars, that this event may be rare in the lifetime of a cluster, producing roughly 15 to 21% of the Intracluster Light (ICL) mass in one go, though perhaps a common event for all galaxy clusters.

preprint2022arXiv

Gemini North Adaptive Optics (GNAO) facility overview and status updates

The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument operating with an f/32 beam can be deployed using GNAO. The GNAO project includes the development of a new laser guide star facility which will consist of four side-launched laser beams supporting the two primary AO modes of GNAO: a wide-field mode providing an improved image quality over natural seeing for a 2-arcminute circular field-of-view and a narrow-field mode providing near diffraction-limited performance over a 20x20 arcsecond square field-of-view. The GNAO wide field mode will enable GIRMOS's multi-IFU configuration in which the science beam to each individual IFU will be additionally corrected using multi-object AO within GIRMOS. The GNAO narrow field mode will feed the GIRMOS tiled IFU configuration in which all IFUs are combined into a "super"-IFU in the center of the field. GNAO also includes the development of a new Real Time Controller, a new GNAO Facility System Controller and finally the development of a new AO Bench. We present in this paper an overview of the GNAO facility and provide a status update of each product.

preprint2022arXiv

LEGA-C: analysis of dynamical masses from ionized gas and stellar kinematics at z~0.8

We compare dynamical mass estimates based on spatially extended stellar and ionized gas kinematics ($\mathrm{M_{dyn,*}}$ and $\mathrm{M_{dyn,eml}}$, respectively) of 157 star forming galaxies at $0.6\leq z<1$. Compared to $z\sim0$, these galaxies have enhanced star formation rates, with stellar feedback likely affecting the dynamics of the gas. We use LEGA-C DR3, the highest redshift dataset providing sufficiently deep measurements of a $K_s-$band limited sample. For $\mathrm{M_{dyn,*}}$ we use Jeans Anisotropic Multi-Gaussian Expansion models. For $\mathrm{M_{dyn,eml}}$ we first fit a custom model of a rotating exponential disk with uniform dispersion, whose light is projected through a slit and corrected for beam smearing. We then apply an asymmetric drift correction based on assumptions common in the literature to the fitted kinematic components to obtain the circular velocity, assuming hydrostatic equilibrium. Within the half-light radius, $\mathrm{M_{dyn,eml}}$ is on average lower than $\mathrm{M_{dyn,*}}$, with a mean offset of $-0.15\pm0.016$ dex and galaxy-to-galaxy scatter of $0.19$ dex, reflecting the combined random uncertainty. While data of higher spatial resolution are needed to understand this small offset, it supports the assumption that the galaxy-wide ionized gas kinematics do not predominantly originate from disruptive events such as star formation driven outflows. However, a similar agreement can be obtained without modeling from the integrated emission line dispersions for axis ratios $q<0.8$. This suggests that our current understanding of gas kinematics is not sufficient to efficiently apply asymmetric drift corrections to improve dynamical mass estimates compared to observations lacking the $S/N$ required for spatially extended dynamics.

preprint2022arXiv

Resolved Stellar Mass Maps of Galaxies in the Hubble Frontier Fields: Evidence for Mass Dependency in Environmental Quenching

One of the challenges in understanding the quenching processes for galaxies is connecting progenitor star-forming populations to their descendant quiescent populations over cosmic time. Here we attempt a novel approach to this challenge by assuming that the underlying stellar mass distribution of galaxies is not significantly altered during environmental quenching processes that solely affect the gas content of cluster galaxies, such as strangulation and ram-pressure stripping. Using the deep, high-resolution photometry of the Hubble Frontier Fields, we create resolved stellar mass maps for both cluster and field galaxies, from which we determine 2D Sérsic profiles, and obtain Sérsic indices and half-mass radii. We classify the quiescent cluster galaxies into disk-like and bulge-like populations based on their Sérsic indices, and find that bulge-like quiescent galaxies dominate the quiescent population at higher masses ($M_\star > 10^{9.5}M_\odot$), whereas disk-like quiescent galaxies dominate at lower masses ($10^{8.5}M_\odot< M_\star < 10^{9.5}M_\odot$). Using both the Sérsic indices and half-mass radii, we identify a population of quiescent galaxies in clusters that are "morphological analogues" of field star-forming galaxies. These analogues are interpreted to be star-forming galaxies that had been environmentally quenched. We use these morphological analogues to compute the environmental-quenching efficiency, and we find that the efficiency decreases with increasing stellar mass. This demonstrates that environmental quenching is more effective on less massive galaxies and that the effect of environment on quenching galaxies is not completely separable from the effect of mass on quenching galaxies.

preprint2022arXiv

The GOGREEN Survey: Constraining the Satellite Quenching Timescale in Massive Clusters at $\boldsymbol{z} \gtrsim 1$

We model satellite quenching at $z \sim 1$ by combining $14$ massive ($10^{13.8} < M_{\mathrm{halo}}/\mathrm{M}_{\odot} < 10^{15}$) clusters at $0.8 < z < 1.3$ from the GOGREEN and GCLASS surveys with accretion histories of $56$ redshift-matched analogs from the IllustrisTNG simulation. Our fiducial model, which is parameterized by the satellite quenching timescale ($τ_{\rm quench}$), accounts for quenching in our simulated satellite population both at the time of infall by using the observed coeval field quenched fraction and after infall by tuning $τ_{\rm quench}$ to reproduce the observed satellite quenched fraction versus stellar mass trend. This model successfully reproduces the observed satellite quenched fraction as a function of stellar mass (by construction), projected cluster-centric radius, and redshift and is consistent with the observed field and cluster stellar mass functions at $z \sim 1$. We find that the satellite quenching timescale is mass dependent, in conflict with some previous studies at low and intermediate redshift. Over the stellar mass range probed ($M_{\star}> 10^{10}~\mathrm{M}_{\odot}$), we find that the satellite quenching timescale decreases with increasing satellite stellar mass from $\sim1.6~{\rm Gyr}$ at $10^{10}~\mathrm{M}_{\odot}$ to $\sim 0.6 - 1~{\rm Gyr}$ at $10^{11}~\mathrm{M}_{\odot}$ and is roughly consistent with the total cold gas (H{\scriptsize I}+H$_{2}$) depletion timescales at intermediate $z$, suggesting that starvation may be the dominant driver of environmental quenching at $z < 2$. Finally, while environmental mechanisms are relatively efficient at quenching massive satellites, we find that the majority ($\sim65-80\%$) of ultra-massive satellites ($M_{\star} > 10^{11}~\mathrm{M}_{\odot}$) are quenched prior to infall.

preprint2022arXiv

The LEGA-C and SAMI Galaxy Surveys: Quiescent Stellar Populations and the Mass-Size Plane across 6 Gyr

We investigate the change in mean stellar population age and metallicity ([Z/H]) scaling relations for quiescent galaxies from intermediate redshift ($0.60\leq z\leq0.76$) using the LEGA-C Survey, to low redshift ($0.014\leq z\leq0.10$) using the SAMI Galaxy Survey. We find that, similarly to their low-redshift counterparts, the stellar metallicity of quiescent galaxies at $0.60\leq z\leq 0.76$ closely correlates with $M_*/R_\mathrm{e}$ (a proxy for the gravitational potential or escape velocity), in that galaxies with deeper potential wells are more metal-rich. This supports the hypothesis that the relation arises due to the gravitational potential regulating the retention of metals, by determining the escape velocity required by metal-rich stellar and supernova ejecta to escape the system and avoid being recycled into later stellar generations. On the other hand, we find no correlation between stellar age and $M_*/R_\mathrm{e}^2$ (stellar mass surface density $Σ$) in the LEGA-C sample, despite this being a strong relation at low redshift. We consider this change in the age--$Σ$ relation in the context of the redshift evolution of the star-forming and quiescent populations in the mass--size plane, and find our results can be explained as a consequence of galaxies forming more compactly at higher redshifts, and remaining compact throughout their evolution. Furthermore, galaxies appear to quench at a characteristic surface density that decreases with decreasing redshift. The $z\sim 0$ age--$Σ$ relation is therefore a result of building up the quiescent and star-forming populations with galaxies that formed at a range of redshifts and so a range of surface densities.

preprint2022arXiv

The Mass Scale of High-Redshift Galaxies: Virial Mass Estimates Calibrated with Stellar Dynamical Models from LEGA-C

Dynamical models for $673$ galaxies at $z=0.6-1.0$ with spatially resolved (long-slit) stellar kinematic data from LEGA-C are used to calibrate virial mass estimates defined as $M_{\rm{vir}}=K σ'^2_{\star,\rm{int}} R$, with $K$ a scaling factor, $σ'_{\star,\rm{int}}$ the spatially-integrated stellar velocity second moment from the LEGA-C survey and $R$ the effective radius measured from a Sérsic profile fit to HST imaging. The sample is representative for $M_{\star}>3\times10^{10}~M_{\odot}$ and includes all types of galaxies, irrespective of morphology and color. We demonstrate that using $R=R_{\rm{sma}}$~(the semi-major axis length of the ellipse that encloses 50\% of the light) in combination with an inclination correction on $σ'_{\star,\rm{int}}$~produces an unbiased $M_{\rm{vir}}$. We confirm the importance of projection effects on $σ'_{\star,\rm{int}}$ by showing the existence of a similar residual trend between virial mass estimates and inclination for the nearby early-type galaxies in the ATLAS$^{\rm{3D}}$~survey. Also, as previously shown, when using a Sérsic profile-based $R$ estimate, then a Sérsic index-dependent correction to account for non-homology in the radial profiles is required. With respect to analogous dynamical models for low-redshift galaxies from the ATLAS$^{\rm{3D}}$~survey we find a systematic offset of 0.1 dex in the calibrated virial constant for LEGA-C, which may be due to physical differences between the galaxy samples or an unknown systematic error. Either way, with our work we establish a common mass scale for galaxies across 8 Gyr of cosmic time with a systematic uncertainty of at most 0.1 dex.

preprint2022arXiv

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- II. Wide Field Slitless Spectroscopy

We present the wide field slitless spectroscopy mode of the NIRISS instrument on the James Webb Space Telescope. This mode employs two orthogonal low-resolution (resolving power $\approx 150$) grisms in combination with a set of six blocking filters in the wavelength range 0.8 to $2.3\,μ$m to provide a spectrum of almost every source across the field-of-view. When combined with the low background, high sensitivity and high spatial resolution afforded by the telescope, this mode will enable unprecedented studies of the structure and evolution of distant galaxies. We describe the performance of the as-built hardware relevant to this mode and expected imaging and spectroscopic sensitivity. We discuss operational and calibration procedures to obtain the highest quality data. As examples of the observing mode usage, we present details of two planned Guaranteed Time Observations programs: The Canadian NIRISS Unbiased Cluster Survey (CANUCS) and The NIRISS Survey for Young Brown Dwarfs and Rogue Planets.

preprint2021arXiv

Diagnosing DASH: A Catalog of Structural Properties for the COSMOS-DASH Survey

We present the $H_{160}$ morphological catalogs for the COSMOS-DASH survey, the largest area near-IR survey using HST-WFC3 to date. Utilizing the "Drift And SHift" observing technique for HST-WFC3 imaging, the COSMOS-DASH survey imaged approximately 0.5 deg$^2$ of the UltraVISTA deep stripes (0.7 deg$^2$ when combined with archival data). Global structural parameters are measured for 51,586 galaxies within COSMOS-DASH using GALFIT (excluding the CANDELS area) with detection using a deep multi-band HST image. We recover consistent results with those from the deeper 3D-HST morphological catalogs, finding that, in general, sizes and Sérsic indices of typical galaxies are accurate to limiting magnitudes of $H_{160}<23$ and $H_{160}<22$ ABmag, respectively. In size-mass parameter space, galaxies in COSMOS-DASH demonstrate robust morphological measurements out to $z\sim2$ and down to $\log(M_\star/M_\odot)\sim9$. With the advantage of the larger area of COSMOS-DASH, we measure a flattening of the quiescent size-mass relation below $\log(M_\star/M_\odot)\sim10.5$ that persists out to $z\sim2$. We show that environment is not the primary driver of this flattening, at least out to $z=1.2$, whereas internal physical processes may instead govern the structural evolution.

preprint2021arXiv

Finite Resolution Deconvolution of Multi-Wavelength Imaging of 20,000 Galaxies in the COSMOS Field: The Evolution of Clumpy Galaxies Over Cosmic Time

Compact star-forming clumps observed in distant galaxies are often suggested to play a crucial role in galaxy assembly. In this paper, we use a novel approach of applying finite resolution deconvolution on ground-based images of the COSMOS field to resolve 20,185 star-forming galaxies (SFG) at 0.5<z<2 to an angular resolution of 0.3", and study their clumpy fractions. A comparison between the deconvolved and HST images across four different filters shows good agreement and validates the deconvolution. We model spectral energy distributions using the deconvolved 14-band images to provide resolved surface brightness and stellar mass density maps for these galaxies. We find that the fraction of clumpy galaxies decreases with increasing stellar masses, and with increasing redshift: from ~30% at z ~ 0.7 to ~50% at z ~ 1.7. Using abundance matching, we also trace the progenitors for galaxies at z ~ 0.7 and measure the fractional mass contribution of clumps toward their total mass budget. Clumps are observed to have a higher fractional mass contribution toward galaxies at higher redshift: increasing from ~1% at z ~ 0.7 to ~5% at z ~ 1.7. Finally, the majority of clumpy SFGs have higher specific star formation rates (sSFR) compared to the average SFGs at fixed stellar mass. We discuss the implication of this result to in-situ clump formation due to disk instability.

preprint2021arXiv

HST/WFC3 grism observations of $z\sim1$ clusters: Evidence for rapid outside-in environmental quenching from spatially resolved H$α$ maps

We present and publicly release (https://www.gclasshst.com) the first spatially resolved H$α$ maps of star-forming cluster galaxies at $z\sim1$, made possible with the Wide Field Camera 3 (WFC3) G141 grism on the Hubble Space Telescope (HST). Using a similar but updated method to 3D-HST in the field environment, we stack the H$α$ maps in bins of stellar mass, measure the half-light radius of the H$α$ distribution and compare it to the stellar continuum. The ratio of the H$α$ to stellar continuum half-light radius, $R[\mathrm{H}α/\mathrm{C}]=\frac{R_{\mathrm{eff, H}α}}{R_{\mathrm{eff, Cont}}}$, is smaller in the clusters by $(6\pm9)\%$, but statistically consistent within $1σ$ uncertainties. A negligible difference in $R[\mathrm{H}α/\mathrm{C}]$ with environment is surprising, given the higher quenched fractions in the clusters relative to the field. We postulate that the combination of high quenched fractions and no change in $R[\mathrm{H}α/\mathrm{C}]$ with environment can be reconciled if environmental quenching proceeds rapidly. We investigate this hypothesis by performing similar analysis on the spectroscopically-confirmed recently quenched cluster galaxies. 87% have H$α$ detections, with star formation rates $8\pm1$ times lower than star-forming cluster galaxies of similar stellar mass. Importantly, these galaxies have a $R[\mathrm{H}α/\mathrm{C}]$ that is $(81\pm8)\%$ smaller than coeval star-forming field galaxies at fixed stellar mass. This suggests the environmental quenching process occurred outside-in. We conclude that disk truncation due to ram-pressure stripping is occurring in cluster galaxies at $z\sim1$, but more rapidly and/or efficiently than in $z\lesssim0.5$ clusters, such that the effects on $R[\mathrm{H}α/\mathrm{C}]$ become observable just after the cluster galaxy has recently quenched.

preprint2021arXiv

Low Surface Brightness Galaxies in z > 1 Galaxy Clusters: HST approaches the Progenitors of Local Ultra Diffuse Galaxies

Ultra Diffuse Galaxies (UDGs), a type of large Low Surface Brightness (LSB) galaxies with particularly large effective radii (r_eff > 1.5 kpc), are now routinely studied in the local (z<0.1) universe. While they are found to be abundant in clusters, groups, and in the field, their formation mechanisms remain elusive and an active topic of debate. New insights may be found by studying their counterparts at higher redshifts (z>1.0), even though cosmological surface brightness dimming makes them particularly diffcult to detect and study there. This work uses the deepest Hubble Space Telescope (HST) imaging stacks of z > 1 clusters, namely: SPT-CL J2106-5844 and MOO J1014+0038. These two clusters, at z=1.13 and z=1.23, were monitored as part of the HST See-Change program. Compared to the Hubble Extreme Deep Field (XDF) as reference field, we find statistical over-densities of large LSB galaxies in both clusters. Based on stellar population modelling and assuming no size evolution, we find that the faintest sources we can detect are about as bright as expected for the progenitors of the brightest local UDGs. We find that the LSBs we detect in SPT-CL J2106-5844 and MOO J1014-5844 already have old stellar populations that place them on the red sequence. Correcting for incompleteness, and based on an extrapolation of local scaling relations, we estimate that distant UDGs are relatively under-abundant compared to local UDGs by a factor ~ 3. A plausible explanation for the implied increase with time would be a significant size growth of these galaxies in the last ~ 8 Gyr, as also suggested by hydrodynamical simulations.

preprint2021arXiv

SITELLE Hα Imaging Spectroscopy of z~0.25 Clusters: Emission Line Galaxy Detection and Ionized Gas Offset in Abell 2390 & Abell 2465

Environmental effects are crucial to the understanding of the evolution of galaxies in dense environments, such as galaxy clusters. Using the large field-of-view of SITELLE, the unique imaging fourier transform spectrograph at CFHT, we are able to obtain 2D spectral information for a large and complete sample of cluster galaxies out to the infall region. We describe a pipeline developed to identify emission line galaxies (ELGs) from the datacube using cross-correlation techniques. We present results based on the spatial offsets between the emission-line regions and stellar continua in ELGs from two z$\sim$0.25 galaxy clusters, Abell 2390 and Abell 2465. We find a preference in the offsets being pointed away from the cluster center. Combining the two clusters, there is a 3$σ$ excess for high-velocity galaxies within the virial radius having the offsets to be pointed away from the cluster center. Assuming the offset being a proxy for the velocity vector of a galaxy, as expected from ram pressure stripping, this excess indicates that ram pressure stripping occurs most effectively during the first passage of an infalling galaxy, leading to the quenching of its star formation. We also find that, outside the virial region, the continuum-normalized H$α$ line flux for infalling galaxies with large offsets are on average lower than those with small or no measurable offset, further supporting ram pressure as a dominant quenching mechanism during the initial infall stages.

preprint2021arXiv

Spectroscopic Confirmation of a Protocluster at $z=3.37$ with a High Fraction of Quiescent Galaxies

We report the discovery of MAGAZ3NE J095924+022537, a spectroscopically-confirmed protocluster at $z = 3.3665^{+0.0009}_{-0.0012}$ around a spectroscopically-confirmed $UVJ$-quiescent ultra-massive galaxy (UMG; $M_{\star}=2.34^{+0.23}_{-0.34}\times10^{11} {\rm M}_\odot$) in the COSMOS UltraVISTA field. We present a total of 38 protocluster members (14 spectroscopic and 24 photometric), including the UMG. Notably, and in marked contrast to protoclusters previously reported at this epoch which have been found to contain predominantly star-forming members, we measure an elevated fraction of quiescent galaxies relative to the coeval field ($73.3^{+26.7}_{-16.9}\%$ versus $11.6^{+7.1}_{-4.9}\%$ for galaxies with stellar mass $M_{\star} \geq 10^{11} {\rm M}_\odot$). This high quenched fraction provides a striking and important counterexample to the seeming ubiquitousness of star-forming galaxies in protoclusters at $z>2$ and suggests, rather, that protoclusters exist in a diversity of evolutionary states in the early Universe. We discuss the possibility that we might be observing either "early mass quenching" or non-classical "environmental quenching." We also present the discovery of MAGAZ3NE J100028+023349, a second spectroscopically-confirmed protocluster, at a very similar redshift of $z = 3.3801^{+0.0213}_{-0.0281}$. We present a total of 20 protocluster members, 12 of which are photometric and 8 spectroscopic including a post-starburst UMG ($M_{\star}=2.95^{+0.21}_{-0.20}\times10^{11} {\rm M}_\odot$). Protoclusters MAGAZ3NE J0959 and MAGAZ3NE J1000 are separated by 18 arcminutes on the sky (35 comoving Mpc), in good agreement with predictions from simulations for the size of "Coma"-type cluster progenitors at this epoch. It is highly likely that the two UMGs are the progenitors of Brightest Cluster Galaxies (BCGs) seen in massive virialized clusters at lower redshift.

preprint2021arXiv

The LEGA-C of nature and nurture in stellar populations of galaxies at z~0.6-1.0: D4000 and H-delta reveal different assembly histories for quiescent galaxies in different environments

Galaxy evolution is driven by a variety of physical processes which are predicted to proceed at different rates for different dark matter haloes and environments across cosmic times. A record of this evolution is preserved in galaxy stellar populations, which we can access using absorption-line spectroscopy. Here we explore the large LEGA-C survey (DR3) to investigate the role of the environment and stellar mass on stellar populations at z~0.6-1.0 in the COSMOS field. Leveraging the statistical power and depth of LEGA-C, we reveal significant gradients in D4000 and H-delta equivalent widths (EWs) distributions over the stellar mass vs environment 2D spaces for the massive galaxy population (M>10^10 M$_{\odot}$) at z~0.6-1.0. D4000 and H-delta EWs primarily depend on stellar mass, but they also depend on environment at fixed stellar mass. By splitting the sample into centrals and satellites, and in terms of star-forming galaxies and quiescent galaxies, we reveal that the significant environmental trends of D4000 and H-delta EW when controlling for stellar mass are driven by quiescent galaxies. Regardless of being centrals or satellites, star-forming galaxies reveal D4000 and H-delta EWs which depend strongly on their stellar mass and are completely independent of the environment at 0.6<z<1.0. The environmental trends seen for satellite galaxies are fully driven by the trends that hold only for quiescent galaxies, combined with the strong environmental dependency of the quiescent fraction at fixed stellar mass. Our results are consistent with recent predictions from simulations that point towards massive galaxies forming first in over-densities or the most compact dark matter haloes.

preprint2021arXiv

Ubiquitous [OII] Emission in Quiescent Galaxies at z ~ 0.85

Using deep rest-frame optical spectroscopy from the Large Early Galaxy Astrophysical Census (LEGA-C) survey, conducted using VIMOS on the ESO Very Large Telescope, we systematically search for low-ionization [OII] 3726,3729 emission in the spectra of a mass-complete sample of z~0.85 galaxies. Intriguingly, we find that 59 percent of UVJ-quiescent (i.e. non star-forming) galaxies in the sample have ionized gas, as traced by [OII] emission, detected above our completeness limit of 1.5 Angstroms. The median stacked spectrum of the lowest equivalent width quiescent galaxies also shows [OII] emission. The overall fraction of sources with [OII] above our equivalent width limit is comparable to what we find in the low-redshift Universe from GAMA and MASSIVE, except perhaps at the highest stellar masses (log Mstar/Msol > 11.5). However, stacked spectra for the individual low-equivalent width systems uniquely indicates ubiquitous [OII] emission in the higher-z LEGA-C sample, with typical [OII] luminosities per unit stellar mass that are a factor of 3 larger than the lower-z GAMA sample. Star formation in these otherwise quiescent galaxies could play a role in producing the [OII] emission at higher-z, although it is unlikely to provide the bulk of the ionizing photons. More work is required to fully quantify the contributions of evolved stellar populations or active galactic nuclei to the observed spectra.

preprint2020arXiv

Constraining the Mass of the Emerging Galaxy Cluster SpARCS1049+56 at z=1.71 with Infrared Weak Lensing

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 $Λ$CDM universe.

preprint2020arXiv

HST/WFC3 grism observations of $z\sim1$ clusters: evidence for evolution in the mass-size relation of quiescent galaxies from poststarburst galaxies

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.

preprint2020arXiv

Inverse stellar population age gradients of post-starburst galaxies at z=0.8 with LEGA-C

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δ_A$ and $Fe4383$, and find that, on average, PSB galaxies have radially decreasing $Hδ_A$ and increasing $Fe4383$ profiles. In contrast, a control sample of quiescent, non-PSB galaxies in the same mass range shows outwardly increasing $Hδ_A$ and decreasing $Fe4383$. The observed gradients are weak ($\approx-0.2$ Å/$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_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.

preprint2020arXiv

Stellar Kinematics and Environment at z~0.8 in the LEGA-C Survey: Massive, Slow-Rotators are Built First in Overdense Environments

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+$δ$) 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+δ)\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+δ)\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.

preprint2020arXiv

The GOGREEN Survey: A deep stellar mass function of cluster galaxies at 1.0<z<1.4 and the complex nature of satellite quenching

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.

preprint2020arXiv

The GOGREEN survey: Post-infall environmental quenching fails to predict the observed age difference between quiescent field and cluster galaxies at z>1

We study the star formation histories (SFHs) and mass-weighted ages of 331 UVJ-selected quiescent galaxies in 11 galaxy clusters and in the field at 1<z<1.5 from the Gemini Observations of Galaxies in Rich Early ENvironments (GOGREEN) survey. We determine the SFHs of individual galaxies by simultaneously fitting rest-frame optical spectroscopy and broadband photometry to stellar population models. We confirm that the SFHs are consistent with more massive galaxies having on average earlier formation times. Comparing galaxies found in massive clusters with those in the field, we find galaxies with $M_\ast<10^{11.3}$ M$_{\odot}$ in the field have more extended SFHs. From the SFHs we calculate the mass-weighted ages, and compare age distributions of galaxies between the two environments, at fixed mass. We constrain the difference in mass-weighted ages between field and cluster galaxies to $0.31_{^{-0.33}}^{_{+0.51}}$ Gyr, in the sense that cluster galaxies are older. We place this result in the context of two simple quenching models and show that neither environmental quenching based on time since infall (without pre-processing) nor a difference in formation times alone can reproduce both the average age difference and relative quenched fractions. This is distinctly different from local clusters, for which the majority of the quenched population is consistent with having been environmentally quenched upon infall. Our results suggest that quenched population in galaxy clusters at z>1 has been driven by different physical processes than those at play at z=0.

preprint2020arXiv

The Number Densities and Stellar Populations of Massive Galaxies at 3 < z < 6: A Diverse, Rapidly Forming Population in the Early Universe

We present the census of massive (log(M$_{*}$/M$_{\odot}$)$\geq 11$) galaxies at $3<z<6$ identified over the COSMOS/UltraVISTA Ultra-Deep field stripes: consisting of $\approx100$ and $\approx20$ high-confidence candidates at $3<z<4$ and $4<z<6$, respectively. The $3<z<4$ population is comprised of post-starburst, UV star-forming and dusty-star forming galaxies in roughly equal fractions, while UV-star-forming galaxies dominate at $4<z<6$ . We account for various sources of biases in SED modelling, finding that the treatment of emission line contamination is essential for understanding the number densities and mass growth histories of massive galaxies at $z>3$. The significant increase in observed number densities at $z\sim4$ ($>\times$ 5 in $\lesssim600$ Myrs) implies that this is the epoch at which log(M$_{*}$/M$_{\odot}$)$\geq 11$ galaxies emerge in significant numbers, with stellar ages ($\approx500-900$ Myrs) indicating rapid formation epochs as early as $z\sim7$. Leveraging ancillary multi-wavelength datasets, we perform panchromatic SED modelling to constrain the total star-formation activity of the sample. The star-formation activity of the sample is generally consistent with being on the star-formation main sequence at the considered redshifts, with $\approx15-25\%$ of the population showing evidence of suppressed star-formation rates, indicating that quenching mechanisms are already at play by $z\sim4$. We stack available HST imaging, confirming their compact nature ($r_{e}\lesssim2.2$ kpc), consistent with expected sizes of high-$z$ star-forming galaxies. Finally, we discuss how our results are in-line with the early formation epochs and short formation timescales inferred from the fossil records of the most massive galaxies in the Universe.

preprint2019arXiv

An Extremely Massive Quiescent Galaxy at $z=3.493$: Evidence of Insufficiently Rapid Quenching Mechanisms in Theoretical Models

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.

preprint2019arXiv

The colors and sizes of recently quenched galaxies: a result of compact starburst before quenching

We analyze the colors and sizes of 32 quiescent (UVJ-selected) galaxies with strong Balmer absorption ($\mbox{EW}(Hδ) \geq 4$Å) 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$.

preprint2018arXiv

HST F160W Imaging of Very Massive Galaxies at $1.5<z<3.0$: Diversity of Structures and the Effect of Close Pairs on Number Density Estimates

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.

preprint2016arXiv

A New Method for Wide-Field Near-IR Imaging with the Hubble Space Telescope

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.

preprint2016arXiv

Stellar mass function of cluster galaxies at z $\sim$ 1.5: evidence for reduced quenching efficiency at high redshift

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

preprint2016arXiv

Stellar Velocity Dispersion for a Strongly-Lensed, Intermediate-Mass Quiescent Galaxy at z=2.8

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δ$. 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 ($σ=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.

preprint2016arXiv

The abundance and spatial distribution of ultra-diffuse galaxies in nearby galaxy clusters

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.

preprint2016arXiv

The evolution in the stellar mass of Brightest Cluster Galaxies over the past 10 billion years

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.

preprint2016arXiv

The Evolution of the Fractions of Quiescent and Star-forming Galaxies as a Function of Stellar Mass Since z=3: Increasing Importance of Massive, Dusty Star-forming Galaxies in the Early Universe

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

preprint2016arXiv

Ultra-deep Ks-band Imaging of the Hubble Frontier Fields

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

preprint2015arXiv

An Extreme Starburst in Close Proximity to the Central Galaxy of a Rich Galaxy Cluster at z=1.7

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.

preprint2015arXiv

Evidence for a change in the dominant satellite galaxy quenching mechanism at z=1

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

preprint2015arXiv

Evidence for the inside-out growth of the stellar mass distribution in galaxy clusters since z~1

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.

preprint2015arXiv

Evidence for the Universality of Properties of Red-Sequence Galaxies in X-ray- and Red-Sequence-Selected Clusters at z ~ 1

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

preprint2015arXiv

Mergers and Star Formation: The environment and Stellar Mass Growth of the Progenitors of Ultra-Massive Galaxies since z = 2

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.

preprint2015arXiv

Spectroscopic Confirmation of an Ultra Massive and Compact Galaxy at z=3.35: A Detailed Look at an Early Progenitor of Local Most Massive Ellipticals

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

preprint2015arXiv

Stellar mass functions of galaxies at 4<z<7 from an IRAC-selected sample in COSMOS/UltraVISTA: limits on the abundance of very massive galaxies

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μ$m detection suggesting the presence of a powerful obscured AGN.

preprint2015arXiv

The Canadian Cluster Comparison Project: detailed study of systematics and updated weak lensing masses

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.

preprint2015arXiv

The Star-Formation History of BCGs to z = 1.8 from the SpARCS/SWIRE Survey: Evidence for significant in-situ star formation at high-redshift

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.

preprint2014arXiv

Cookie-cutter halos: the remarkable constancy of the stellar mass function of satellite galaxies at 0.2<z<1.2

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.

preprint2014arXiv

Dense cores in galaxies out to z=2.5 in SDSS, UltraVISTA, and the five 3D-HST/CANDELS fields

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.

preprint2014arXiv

Observations of environmental quenching in groups in the 11 Gyr since z=2.5: different quenching for central and satellite galaxies

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.

preprint2014arXiv

The Phase Space and Stellar Populations of Cluster Galaxies at z ~ 1: Simultaneous Constraints on the Location and Timescale of Satellite Quenching

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.

preprint2014arXiv

The Progenitors of Local Ultra-massive Galaxies Across Cosmic Time: from Dusty Star-bursting to Quiescent Stellar Populations

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.

preprint2014arXiv

What Regulates Galaxy Evolution? Open Questions in Our Understanding of Galaxy Formation and Evolution

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.

preprint2013arXiv

A Census of Stellar Mass in 10 Massive Haloes at z~1 from the GCLASS Survey

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]

preprint2013arXiv

A Public Ks-selected Catalog in the COSMOS/UltraVISTA Field: Photometry, Photometric Redshifts and Stellar Population Parameters

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.

preprint2013arXiv

Discovery of a Rich Cluster at z = 1.63 using the Rest-Frame 1.6um "Stellar Bump Sequence" Method

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.

preprint2013arXiv

HST/WFC3 Confirmation of the Inside-Out Growth of Massive Galaxies at 0<z<2 and Identification of their Star Forming Progenitors at z~3

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.

preprint2013arXiv

The environmental dependence of the stellar mass function at z~1: Comparing cluster and field between the GCLASS and UltraVISTA surveys

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.

preprint2013arXiv

The Evolution of Dusty Star formation in Galaxy Clusters to z = 1: Spitzer IR Observations of the First Red-Sequence Cluster Survey

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.

preprint2013arXiv

The Evolution of the Stellar Mass Functions of Star-Forming and Quiescent Galaxies to z = 4 from the COSMOS/UltraVISTA Survey

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.

preprint2013arXiv

The Structural Evolution of Milky Way-like Star Forming Galaxies since z~1.3

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.

preprint2012arXiv

3D-HST: A wide-field grism spectroscopic survey with the Hubble Space Telescope

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.

preprint2012arXiv

A Strongly-Lensed Massive Ultra-Compact Quiescent Galaxy at z ~ 2.4 in the COSMOS/UltraVISTA Field

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.

preprint2011arXiv

First Results from the 3D-HST Survey: The Striking Diversity of Massive Galaxies at z>1

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.

preprint2011arXiv

Lensing Magnification: A novel method to weigh high-redshift clusters and its application to SpARCS

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.

preprint2011arXiv

The Gemini Cluster Astrophysics Spectroscopic Survey (GCLASS): The Role of Environment and Self-Regulation in Galaxy Evolution at z ~ 1

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.

preprint2011arXiv

The NEWFIRM Medium-band Survey: Photometric Catalogs, Redshifts and the Bimodal Color Distribution of Galaxies out to z~3

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.

preprint2011arXiv

The number density and mass density of star-forming and quiescent galaxies at 0.4 < z < 2.2

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

preprint2010arXiv

Galaxy clustering in the NEWFIRM Medium Band Survey: the relationship between stellar mass and dark matter halo mass at 1 < z < 2

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.

preprint2010arXiv

Spectroscopic Confirmation of Three Red-Sequence Selected Galaxy Clusters at z=0.87, 1.16 and 1.21 from the SpARCS Survey

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.

preprint2010arXiv

The Most Massive Galaxies at 3.0<z<4.0 in the NEWFIRM Medium-Band Survey: Properties and Improved Constraints on the Stellar Mass Function

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

preprint2010arXiv

The Spectral Energy Distribution of Post-Starburst Galaxies in the NEWFIRM Medium-Band Survey: A Low Contribution from TP-AGB Stars

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.

preprint2010arXiv

Well-Sampled Far-Infrared Spectral Energy Distributions of z ~ 2 Galaxies: Evidence for Scaled up Cool Galaxies

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.

preprint2009arXiv

A Near-Infrared Spectroscopic Survey of K-selected Galaxies at z ~ 2.3: Comparison of Stellar Population Synthesis Codes and Constraints from the Rest-frame NIR

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.

preprint2009arXiv

How Massive are Massive Compact Galaxies?

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.

preprint2009arXiv

The Growth of Massive Galaxies Since z=2

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.

preprint2009arXiv

The NEWFIRM Medium-Band Survey: Filter Definitions and First Results

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.

preprint2008arXiv

Spectroscopic Confirmation of a Massive Red-Sequence-Selected Galaxy Cluster at z = 1.34 in the SpARCS-South Cluster Survey

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.

preprint2008arXiv

Spectroscopic Confirmation of Two Massive Red-Sequence-Selected Galaxy Clusters at z ~ 1.2 in the SpARCS-North Cluster Survey

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.

preprint2006arXiv

Near-Infrared Properties of Moderate-Redshift Galaxy Clusters: Luminosity Functions and Density Profiles

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.