Source author record

Katherine E. Whitaker

Katherine E. Whitaker 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

67works
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

67 published item(s)

preprint2025arXiv

Picture an Astronomer: Best Practices for Retaining Talent in Astrophysics

Women are consistently underrepresented in astrophysics yet are simultaneously subject to disproportionate attrition at every career stage. This disparity between demonstrated efficacy in job performance and ultimate career outcome was the primary motivation for the Picture an Astronomer series, which included both targeted public outreach to increase representation of women in astrophysics and high-level, solution-oriented discussions among professional astronomers. In March 2025, more than 200 astronomers came together in a hybrid-format symposium focused on the state of the field for female scientists, combining scientific exchange with discussions of policies and practices to strengthen retention of talent in the field. This white paper is the result of those discussions, offering a wide range of recommendations developed in the context of gendered attrition in astrophysics but which ultimately support a healthier climate for all scientists alike.

preprint2023arXiv

Early JWST imaging reveals strong optical and NIR color gradients in galaxies at $z\sim2$ driven mostly by dust

Recent studies have shown that galaxies at cosmic noon are redder in the center and bluer in the outskirts, mirroring results in the local universe. These color gradients could be caused by either gradients in the stellar age or dust opacity; however, distinguishing between these two causes is impossible with rest-frame optical photometry alone. Here we investigate the underlying causes of the gradients from spatially-resolved rest-frame $U-V$ vs. $V-J$ color-color diagrams, measured from early observations with the James Webb Space Telescope. We use $1\, μm - 4\, μm$ NIRCam photometry from the CEERS survey of a sample of 54 galaxies with $M_* / M_\odot>10$ at redshifts $1.7<z<2.3$ selected from the 3D-HST catalog. We model the light profiles in the F115W, F200W and F356W NIRCam bands using \texttt{imcascade}, a Bayesian implementation of the Multi-Gaussian expansion (MGE) technique which flexibly represents galaxy profiles using a series of Gaussians. We construct resolved rest-frame $U-V$ and $V-J$ color profiles. The majority of star-forming galaxies have negative gradients (i.e. redder in the center, bluer in the outskirts) in both $U-V$ and $V-J$ colors consistent with radially decreasing dust attenuation. A smaller population (roughly 15\%) of star-forming galaxies have positive $U-V$ but negative $V-J$ gradients implying centrally concentrated star-formation. For quiescent galaxies we find a diversity of UVJ color profiles, with roughly one-third showing star-formation in their center. This study showcases the potential of JWST to study the resolved stellar populations of galaxies at cosmic noon.

preprint2023arXiv

UNCOVER: Illuminating the Early Universe -- JWST/NIRSpec Confirmation of $z > 12$ Galaxies

Observations of high-redshift galaxies provide a critical direct test to the theories of early galaxy formation, yet to date, only three have been spectroscopically confirmed at $z>12$. Due to strong gravitational lensing over a wide area, the galaxy cluster field A2744 is ideal for searching for the earliest galaxies. Here we present JWST/NIRSpec observations of two galaxies: a robust detection at $z_{\rm spec} = 12.393^{+0.004}_{-0.001}$, and a plausible candidate at $z_{\rm spec} = 13.079^{+0.013}_{-0.001}$. The galaxies are discovered in JWST/NIRCam imaging and their distances are inferred with JWST/NIRSpec spectroscopy, all from the JWST Cycle 1 UNCOVER Treasury survey. Detailed stellar population modeling using JWST NIRCam and NIRSpec data corroborates the primeval characteristics of these galaxies: low mass ($\sim 10^8~{\rm M_\odot}$), young, rapidly-assembling, metal-poor, and star-forming. Interestingly, both galaxies are spatially resolved, having lensing-corrected rest-UV effective radii on the order of 300-400 pc, which are notably larger than other spectroscopically confirmed systems at similar redshifts. The observed dynamic range of $z \gtrsim 10$ sizes spans over 1 order of magnitude, implying a significant scatter in the size-mass relation at early times. Deep into the epoch of reionization, these discoveries elucidate the emergence of the first galaxies.

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

ALMA measures molecular gas reservoirs comparable to field galaxies in a low-mass galaxy cluster at z=1.3

We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas measurements, bridging the observational gap between galaxies in the more extreme, well-studied clusters (Log$_{10}$~M$_{\rm halo}/M_{\odot}\gtrsim14$) and those in group or field environments at cosmic noon. The CO sources are concentrated on the sky (within ~1-arcmin diameter) and phase space analysis indicates the gas resides in galaxies already within the cluster environment. We find that CO sources sit in similar phase space as CO-rich galaxies in more massive clusters at similar redshifts (have similar accretion histories) while maintaining field-like molecular gas reservoirs, compared to scaling relations. This work presents the deepest CO survey to date in a galaxy cluster at z>1, uncovering gas reservoirs down to M$_{\rm H_{2}}>1.6\times10^{10}$M$_{\odot}$ (5$σ$ at 50% primary beam). Our deep limits rule out the presence of gas content in excess of the field scaling relations; however, combined with literature CO detections, cluster gas fractions in general appear systematically high, on the upper envelope or above the field. This study is the first demonstration that low mass clusters at z~1-2 can host overdensities of CO emitters with surviving gas reservoirs, in line with the prediction that quenching is delayed after first infall while galaxies consume the gas bound to the disk.

preprint2022arXiv

Quenching and the UVJ diagram in the SIMBA cosmological simulation

Over the past decade, rest-frame color-color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color-color selection tool -- the rest-frame $U-V$ vs. $V-J$ diagram -- by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties in UVJ space and the color-color evolution of massive galaxies as they quench at redshifts $z\sim 1$--$2$. We find that our models broadly reproduce the observed UVJ diagram at $z=1$--$2$, including (for the first time in a cosmological simulation) reproducing the population of extremely dust-reddened galaxies in the top right of the UVJ diagram. However, our models primarily populate this region with low-mass galaxies and do not produce as clear a bimodality between star-forming and quiescent galaxies as is seen in observations. The former issue is due to an excess of dust in low-mass galaxies and relatively gray attenuation curves in high-mass galaxies, while the latter is due to the overpopulation of the green valley in SIMBA. When investigating the time evolution of galaxies on the UVJ diagram, we find that the quenching pathway on the UVJ diagram is independent of the quenching timescale, and instead dependent primarily on the average specific star formation rate in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies.

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

Schrodinger's Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

$JWST$'s first glimpse of the $z>10$ Universe has yielded a surprising abundance of luminous galaxy candidates. Here we present the most extreme of these systems: CEERS-1749. Based on $0.6-5μ$m photometry, this strikingly luminous ($\approx$26 mag) galaxy appears to lie at $z\approx17$. This would make it an $M_{\rm{UV}}\approx-22$, $M_{\rm{\star}}\approx5\times10^{9}M_{\rm{\odot}}$ system that formed a mere $\sim220$ Myrs after the Big Bang. The implied number density of this galaxy and its analogues challenges virtually every early galaxy evolution model that assumes $Λ$CDM cosmology. However, there is strong environmental evidence supporting a secondary redshift solution of $z\approx5$: all three of the galaxy's nearest neighbors at $<2.5$" have photometric redshifts of $z\approx5$. Further, we show that CEERS-1749 may lie in a $z\approx5$ protocluster that is $\gtrsim5\times$ overdense compared to the field. Intense line emission at $z\approx5$ from a quiescent galaxy harboring ionized gas, or from a dusty starburst, may provide satisfactory explanations for CEERS-1749's photometry. The emission lines at $z\approx5$ conspire to boost the $>2μ$m photometry, producing an apparent blue slope as well as a strong break in the SED. Such a perfectly disguised contaminant is possible only in a narrow redshift window ($Δz\lesssim0.1$), implying that the permitted volume for such interlopers may not be a major concern for $z>10$ searches, particularly when medium-bands are deployed. If CEERS-1749 is confirmed to lie at $z\approx5$, it will be the highest-redshift quiescent galaxy, or one of the lowest mass dusty galaxies of the early Universe detected to-date. Both redshift solutions of this intriguing galaxy hold the potential to challenge existing models of early galaxy evolution, making spectroscopic follow-up of this source critical.

preprint2021arXiv

CLIMBER: Galaxy-Halo Connection Constraints from Next-Generation Surveys

In the coming decade, a new generation of massively multiplexed spectroscopic surveys, such as PFS, WAVES, and MOONS, will probe galaxies in the distant universe in vastly greater numbers than was previously possible. In this work, we generate mock catalogs for each of these three planned surveys to help quantify and optimize their scientific output. To assign photometry into the UniverseMachine empirical model, we develop the Calibrating Light: Illuminating Mocks By Empirical Relations (CLIMBER) procedure using UltraVISTA photometry. Using the published empirical selection functions for each aforementioned survey, we quantify the mass completeness of each survey. We compare different targeting strategies by varying the area and targeting completeness, and quantify how these survey parameters affect the uncertainty of the two-point correlation function. We demonstrate that the PFS and MOONS measurements will be primarily dominated by cosmic variance, not shot noise, motivating the need for increasingly large survey areas. On the other hand, the WAVES survey, which covers a much larger area, will strike a good balance between cosmic variance and shot noise. For a fixed number of targets, a 5% increased survey area (and $\sim$5% decreased completeness) would decrease the uncertainty of the correlation function at intermediate scales by 0.15%, 1.2%, and 1.1% for our WAVES, PFS, and MOONS samples, respectively. Meanwhile, for a fixed survey area, 5% increased targeting completeness improves the same constraints by 0.7%, 0.25%, and 0.1%. All of the utilities used to construct our mock catalogs and many of the catalogs themselves are publicly available.

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

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

ALMA measures rapidly depleted molecular gas reservoirs in massive quiescent galaxies at z~1.5

We present ALMA CO(2-1) spectroscopy of 6 massive (log$_{10}$M$_{\rm{*}}/\rm{M}_\odot>$11.3) quiescent galaxies at $z\sim1.5$. These data represent the largest sample using CO emission to trace molecular gas in quiescent galaxies above $z>1$, achieving an average 3$σ$ sensitivity of M$_{\rm{H_{2}}}\sim10^{10}\rm{M}_\odot$. We detect one galaxy at 4$σ$ significance and place upper limits on the molecular gas reservoirs of the other 5, finding molecular gas mass fractions M$_{\rm{H_{2}}}$/M$_{\rm{*}}$=f$_{\rm{H_{2}}}<2-6$% (3$σ$ upper limits). This is 1-2 orders of magnitude lower than coeval star-forming galaxies at similar stellar mass, and comparable to galaxies at $z=0$ with similarly low sSFR. This indicates that their molecular gas reservoirs were rapidly and efficiently used up or destroyed, and that gas fractions are uniformly low ($<$6%) despite the structural diversity of our sample. The implied rapid depletion time of molecular gas (t$_{\rm{dep}}<0.6$ Gyr) disagrees with extrapolations of empirical scaling relations to low sSFR. We find that our low gas fractions are instead in agreement with predictions from both the recent SIMBA cosmological simulation, and from analytical "bathtub" models for gas accretion onto galaxies in massive dark matter halos (log$_{10}M_{\rm{halo}}/\rm{M}_\odot\sim14$ at $z=0$). Such high mass halos reach a critical mass of log$_{10}M_{\rm{halo}}/\rm{M}_\odot>12$ by $z\sim4$ that halt the accretion of baryons early in the Universe. Our data is consistent with a simple picture where galaxies truncate accretion and then consume the existing gas at or faster than typical main sequence rates. Alternatively, we cannot rule out that these galaxies reside in lower mass halos, and low gas fractions may instead reflect either stronger feedback, or more efficient gas consumption.

preprint2020arXiv

Recent Star-formation in a Massive Slowly-Quenched Lensed Quiescent Galaxy at z=1.88

In this letter, we reconstruct the formation pathway of MRG-S0851, a massive, $\log M_*/M_\odot=11.02\pm0.04$, strongly lensed, red, galaxy at $z=1.883\pm0.001$. While the global photometry and spatially-resolved outskirts of MRG-S0851 imply an early-formation scenario with a slowly decreasing or constant star-formation history, a joint fit of 2D grism spectroscopy and photometry reveals a more complex scenario: MRG-S0851 is likely to be experiencing a centrally-concentrated rejuvenation in the inner $\sim$1 kpc in the last $\sim$100 Myr of evolution. We estimate $0.5\pm0.1\%$ of the total stellar mass is formed in this phase. Rejuvenation episodes are suggested to be infrequent for massive galaxies at $z\sim2$, but as our analyses indicate, more examples of complex star-formation histories may yet be hidden within existing data. By adding a FUV color criterion to the standard U-V/V-J diagnostic, thereby heightening our sensitivity to recent star formation, we show that we can select populations of galaxies with similar spectral energy distributions to that of MRG-S0851, but note that deep follow-up spectroscopic observations and/or spatially resolved analyses are necessary to robustly confirm the rejuvenation of these candidates. Using our criteria with MRG-S0851 as a prototype, we estimate that $\sim$1\% of massive quiescent galaxies at $1<z<2$ are potentially rejuvenating.

preprint2020arXiv

Three Dusty Star Forming Galaxies at $z\sim1.5$: Mergers and Disks on the Main Sequence

The main sequence of galaxies, a correlation between the star formation rates and stellar masses of galaxies, has been observed out to $z\sim4$. Galaxies within the scatter of the correlation are typically interpreted to be secularly evolving while galaxies with star formation rates elevated above the main sequence are interpreted to be undergoing interactions or to be Toomre-unstable disks with starbursting clumps. In this paper we investigate the recent merger histories of three dusty star forming galaxies, identified by their bright submillimeter emission at $z\sim1.5$. We analyze rest-frame optical and UV imaging, rest-frame optical emission line kinematics using slit spectra obtained with MOSFIRE on Keck I, and calculate Gini and M$_{20}$ statistics for each galaxy and conclude two are merger-driven while the third is an isolated disk galaxy. The disk galaxy lies $\sim$4$\times$ above the main sequence, one merger lies within the scatter of the main sequence, and one merger lies $\sim$4$\times$ below the main sequence. This hints that the location of a galaxy with respect to the main sequence may not be a useful discriminator of the recent star formation history of high-M$_{\star}$ galaxies at $z\sim1$.

preprint2019arXiv

Strong Lens Models for 37 Clusters of Galaxies from the SDSS Giant Arcs Survey

We present strong gravitational lensing models for 37 galaxy clusters from the SDSS Giant Arcs Survey. We combine data from multi-band Hubble Space Telescope WFC3imaging, with ground-based imaging and spectroscopy from Magellan, Gemini, APO, and MMT, in order to detect and spectroscopically confirm new multiply-lensed background sources behind the clusters. We report spectroscopic or photometric redshifts of sources in these fields, including cluster galaxies and background sources. Based on all available lensing evidence, we construct and present strong lensing mass models for these galaxy clusters.

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

Ages of massive galaxies at $0.5 < z < 2.0$ from 3D-HST rest-frame optical spectroscopy

We present low-resolution near-infrared stacked spectra from the 3D-HST survey up to $z=2.0$ and fit them with commonly used stellar population synthesis models: BC03 (Bruzual & Charlot, 2003), FSPS10 (Flexible Stellar Population Synthesis, Conroy & Gunn 2010), and FSPS-C3K (Conroy, Kurucz, Cargile, Castelli, in prep). The accuracy of the grism redshifts allows the unambiguous detection of many emission and absorption features, and thus a first systematic exploration of the rest-frame optical spectra of galaxies up to $z=2$. We select massive galaxies ($\rm log(M_{*} / M_{\odot}) > 10.8$), we divide them into quiescent and star-forming via a rest-frame color-color technique, and we median-stack the samples in 3 redshift bins between $z=0.5$ and $z=2.0$. We find that stellar population models fit the observations well at wavelengths below $\rm 6500 Å$ rest-frame, but show systematic residuals at redder wavelengths. The FSPS-C3K model generally provides the best fits (evaluated with a $χ^2_{red}$ statistics) for quiescent galaxies, while BC03 performs the best for star-forming galaxies. The stellar ages of quiescent galaxies implied by the models, assuming solar metallicity, vary from 4 Gyr at $z \sim 0.75$ to 1.5 Gyr at $z \sim 1.75$, with an uncertainty of a factor of 2 caused by the unknown metallicity. On average the stellar ages are half the age of the Universe at these redshifts. We show that the inferred evolution of ages of quiescent galaxies is in agreement with fundamental plane measurements, assuming an 8 Gyr age for local galaxies. For star-forming galaxies the inferred ages depend strongly on the stellar population model and the shape of the assumed star-formation history.

preprint2016arXiv

Evidence for non-stellar rest-frame near-IR emission associated with increased star formation in galaxies at $z \sim 1$

We explore the presence of non-stellar rest-frame near-IR ($2-5 \ μ\mathrm{m}$) emission in galaxies at $z \sim 1$. Previous studies identified this excess in relatively small samples and suggested that such non-stellar emission, which could be linked to the $3.3 \ μ\mathrm{m}$ polycyclic aromatic hydrocarbons feature or hot dust emission, is associated with an increased star formation rate (SFR). In this Letter, we confirm and quantify the presence of an IR excess in a significant fraction of galaxies in the 3D-HST GOODS catalogs. By constructing a matched sample of galaxies with and without strong non-stellar near-IR emission, we find that galaxies with such emission are predominantly star-forming galaxies. Moreover, star-forming galaxies with an excess show increased mid- and far-IR and H$α$ emission compared to other star-forming galaxies without. While galaxies with a near-IR excess show a larger fraction of individually detected X-ray active galactic nuclei (AGNs), an X-ray stacking analysis, together with the IR-colors and H$α$ profiles, shows that AGNs are unlikely to be the dominant source of the excess in the majority of galaxies. Our results suggest that non-stellar near-IR emission is linked to increased SFRs and is ubiquitous among star-forming galaxies. As such, the near-IR emission might be a powerful tool to measure SFRs in the era of the James Webb Space Telescope.

preprint2016arXiv

Low Gas Fractions Connect Compact Star-Forming Galaxies to their z ~ 2 Quiescent Descendants

Early quiescent galaxies at z~2 are known to be remarkably compact compared to their nearby counterparts. Possible progenitors of these systems include galaxies that are structurally similar, but are still rapidly forming stars. Here, we present Karl G. Jansky Very Large Array (VLA) observations of the CO(1-0) line towards three such compact, star-forming galaxies at z~2.3, significantly detecting one. The VLA observations indicate baryonic gas fractions >~5 times lower and gas depletion times >~10 times shorter than normal, extended massive star-forming galaxies at these redshifts. At their current star formation rates, all three objects will deplete their gas reservoirs within 100Myr. These objects are among the most gas-poor objects observed at z>2, and are outliers from standard gas scaling relations, a result which remains true regardless of assumptions about the CO-H2 conversion factor. Our observations are consistent with the idea that compact, star-forming galaxies are in a rapid state of transition to quiescence in tandem with the build-up of the z~2 quenched population. In the detected compact galaxy, we see no evidence of rotation or that the CO-emitting gas is spatially extended relative to the stellar light. This casts doubt on recent suggestions that the gas in these compact galaxies is rotating and significantly extended compared to the stars. Instead, we suggest that, at least for this object, the gas is centrally concentrated, and only traces a small fraction of the total galaxy dynamical mass.

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

The FourStar Galaxy Evolution Survey (ZFOURGE): ultraviolet to far-infrared catalogs, medium-bandwidth photometric redshifts with improved accuracy, stellar masses, and confirmation of quiescent galaxies to z~3.5

The FourStar galaxy evolution survey (ZFOURGE) is a 45 night legacy program with the FourStar near-infrared camera on Magellan and one of the most sensitive surveys to date. ZFOURGE covers a total of $400\ \mathrm{arcmin}^2$ in cosmic fields CDFS, COSMOS and UDS, overlapping CANDELS. We present photometric catalogs comprising $>70,000$ galaxies, selected from ultradeep $K_s$-band detection images ($25.5-26.5$ AB mag, $5σ$, total), and $>80\%$ complete to $K_s<25.3-25.9$ AB. We use 5 near-IR medium-bandwidth filters ($J_1,J_2,J_3,H_s,H_l$) as well as broad-band $K_s$ at $1.05\ - 2.16\ μm$ to $25-26$ AB at a seeing of $\sim0.5$". Each field has ancillary imaging in $26-40$ filters at $0.3-8\ μm$. We derive photometric redshifts and stellar population properties. Comparing with spectroscopic redshifts indicates a photometric redshift uncertainty $σ_z={0.010,0.009}$, and 0.011 in CDFS, COSMOS, and UDS. As spectroscopic samples are often biased towards bright and blue sources, we also inspect the photometric redshift differences between close pairs of galaxies, finding $σ_{z,pairs}= 0.01-0.02$ at $1<z<2.5$. We quantify how $σ_{z,pairs}$ depends on redshift, magnitude, SED type, and the inclusion of FourStar medium bands. $σ_{z,pairs}$ is smallest for bright, blue star-forming samples, while red star-forming galaxies have the worst $σ_{z,pairs}$. Including FourStar medium bands reduces $σ_{z,pairs}$ by 50\% at $1.5<z<2.5$. We calculate SFRs based on ultraviolet and ultradeep far-IR $Spitzer$/MIPS and Herschel/PACS data. We derive rest-frame $U-V$ and $V-J$ colors, and illustrate how these correlate with specific SFR and dust emission to $z=3.5$. We confirm the existence of quiescent galaxies at $z\sim3$, demonstrating their SFRs are suppressed by $>\times15$.

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

First results from the VIRIAL survey: the stellar content of $UVJ$-selected quiescent galaxies at $1.5 < z < 2$ from KMOS

We investigate the stellar populations of 25 massive, galaxies ($\log[M_\ast/M_\odot] \geq 10.9$) at $1.5 < z < 2$ using data obtained with the K-band Multi-Object Spectrograph (KMOS) on the ESO VLT. Targets were selected to be quiescent based on their broadband colors and redshifts using data from the 3D-HST grism survey. The mean redshift of our sample is $\bar{z} = 1.75$, where KMOS YJ-band data probe age- and metallicity-sensitive absorption features in the rest-frame optical, including the $G$ band, Fe I, and high-order Balmer lines. Fitting simple stellar population models to a stack of our KMOS spectra, we derive a mean age of $1.03^{+0.13}_{-0.08}$ Gyr. We confirm previous results suggesting a correlation between color and age for quiescent galaxies, finding mean ages of $1.22^{+0.56}_{-0.19}$ Gyr and $0.85^{+0.08}_{-0.05}$ Gyr for the reddest and bluest galaxies in our sample. Combining our KMOS measurements with those obtained from previous studies at $0.2 < z < 2$ we find evidence for a $2-3$ Gyr spread in the formation epoch of massive galaxies. At $z < 1$ the measured stellar ages are consistent with passive evolution, while at $1 < z \lesssim2$ they appear to saturate at $\sim$1 Gyr, which likely reflects changing demographics of the (mean) progenitor population. By comparing to star-formation histories inferred for "normal" star-forming galaxies, we show that the timescales required to form massive galaxies at $z \gtrsim 1.5$ are consistent with the enhanced $α$-element abundances found in massive local early-type galaxies.

preprint2015arXiv

Forming Compact Massive Galaxies

In this paper we study a key phase in the formation of massive galaxies: the transition of star forming galaxies into massive (M_stars~10^11 Msun), compact (r_e~1 kpc) quiescent galaxies, which takes place from z~3 to z~1.5. We use HST grism redshifts and extensive photometry in all five 3D-HST/CANDELS fields, more than doubling the area used previously for such studies, and combine these data with Keck MOSFIRE and NIRSPEC spectroscopy. We first confirm that a population of massive, compact, star forming galaxies exists at z~2, using K-band spectroscopy of 25 of these objects at 2.0<z<2.5. They have a median NII/Halpha ratio of 0.6, are highly obscured with SFR(tot)/SFR(Halpha)~10, and have a large range of observed line widths. We infer from the kinematics and spatial distribution of Halpha that the galaxies have rotating disks of ionized gas that are a factor of ~2 more extended than the stellar distribution. By combining measurements of individual galaxies, we find that the kinematics are consistent with a nearly Keplerian fall-off from V_rot~500 km/s at 1 kpc to V_rot~250 km/s at 7 kpc, and that the total mass out to this radius is dominated by the dense stellar component. Next, we study the size and mass evolution of the progenitors of compact massive galaxies. Even though individual galaxies may have had complex histories with periods of compaction and mergers, we show that the population of progenitors likely followed a simple inside-out growth track in the size-mass plane of d(log r_e) ~ 0.3 d(log M_stars). This mode of growth gradually increases the stellar mass within a fixed physical radius, and galaxies quench when they reach a stellar density or velocity dispersion threshold. As shown in other studies, the mode of growth changes after quenching, as dry mergers take the galaxies on a relatively steep track in the size-mass plane.

preprint2015arXiv

Galaxy Structure as a Driver of the Star Formation Sequence Slope and Scatter

It is well established that (1) star-forming galaxies follow a relation between their star formation rate (SFR) and stellar mass (M$_{\star}$), the "star-formation sequence", and (2) the SFRs of galaxies correlate with their structure, where star-forming galaxies are less concentrated than quiescent galaxies at fixed mass. Here, we consider whether the scatter and slope of the star-formation sequence is correlated with systematic variations in the Sersic indices, $n$, of galaxies across the SFR-M$_{\star}$ plane. We use a mass-complete sample of 23,848 galaxies at $0.5<z<2.5$ selected from the 3D-HST photometric catalogs. Galaxy light profiles parameterized by $n$ are based on Hubble Space Telescope CANDELS near-infrared imaging. We use a single SFR indicator empirically-calibrated from stacks of Spitzer/MIPS 24$μ$m imaging, adding the unobscured and obscured star formation. We find that the scatter of the star-formation sequence is related in part to galaxy structure; the scatter due to variations in $n$ at fixed mass for star-forming galaxies ranges from 0.14$\pm$0.02 dex at $z\sim2$ to 0.30$\pm$0.04 dex at $z<1$. While the slope of the log(SFR)-log(M$_{\star}$) relation is of order unity for disk-like galaxies, galaxies with $n>2$ (implying more dominant bulges) have significantly lower SFR/M$_{\star}$ than the main ridgeline of the star-formation sequence. These results suggest that bulges in massive $z\sim2$ galaxies are actively building up, where the stars in the central concentration are relatively young. At $z<1$, the presence of older bulges within star-forming galaxies lowers global SFR/M$_{\star}$, decreasing the slope and contributing significantly to the scatter of the star-formation sequence.

preprint2015arXiv

Leveraging 3D-HST Grism Redshifts to Quantify Photometric Redshift Performance

We present a study of photometric redshift accuracy in the 3D-HST photometric catalogs, using 3D-HST grism redshifts to quantify and dissect trends in redshift accuracy for galaxies brighter than $H_{F140W}<24$ with an unprecedented and representative high-redshift galaxy sample. We find an average scatter of $0.0197\pm0.0003(1+z)$ in the Skelton et al. (2014) photometric redshifts. Photometric redshift accuracy decreases with magnitude and redshift, but does not vary monotonically with color or stellar mass. The 1-$σ$ scatter lies between $0.01-0.03$(1+z) for galaxies of all masses and colors below $z<2.5$ (for $H_{F140W}{<}24$), with the exception of a population of very red ($U-V > 2$), dusty star-forming galaxies for which the scatter increases to $\sim0.1(1+z)$. Although the overall photometric redshift accuracy for quiescent galaxies is better than for star-forming galaxies, scatter depends more strongly on magnitude and redshift than on galaxy type. We verify these trends using the redshift distributions of close pairs and extend the analysis to fainter objects, where photometric redshift errors further increase to $\sim0.046(1+z)$ at $H_{F160W}=26$. We demonstrate that photometric redshift accuracy is strongly filter-dependent and quantify the contribution of multiple filter combinations. We evaluate the widths of redshift probability distribution functions and find that error estimates are underestimated by a factor of $\sim1.1-1.6$, but that uniformly broadening the distribution does not adequately account for fitting outliers. Finally, we suggest possible applications of these data in planning for current and future surveys and simulate photometric redshift performance in the LSST, DES, and combined DES and VHS surveys.

preprint2015arXiv

Spatially-resolved dust maps from Balmer decrements in galaxies at z~1.4

We derive average radial gradients in the dust attenuation towards HII regions in 609 galaxies at z~1.4, using measurements of the Balmer decrement out to r~3kpc. The Balmer decrements are derived from spatially resolved maps of Halpha and Hbeta emission from the 3D-HST survey. We find that with increasing stellar mass (M) both the normalization and strength of the gradient in dust attenuation increases. Galaxies with a mean mass of <log(M)> = 9.2Msun have little dust attenuation at all radii, whereas galaxies with <log(M)>= 10.2Msun have dust attenuation toward Halpha A(Halpha)~2mag in their central regions. We parameterize this as A(Halpha) = b + c log(r), with b = 0.9 + 1.0 log(M10), c = -1.9 - 2.2 log(M10), r in kpc, and M10 the stellar mass in units of 10^10Msun. This expression can be used to correct spatially resolved measurements of Halpha to radial distributions of star formation. When applied to our data, we find that the star formation rates in the central r<1kpc of galaxies in the highest mass bin are ~ 6 Msun/yr, six times higher than before correction and approximately half of the total star formation rate of these galaxies. If this high central star formation rate is maintained for several Gyr, a large fraction of the stars in present-day bulges likely formed in-situ.

preprint2015arXiv

The 3D-HST Survey: Hubble Space Telescope WFC3/G141 grism spectra, redshifts, and emission line measurements for $\sim 100,000$ galaxies

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

preprint2015arXiv

The MOSFIRE Deep Evolution Field (MOSDEF) Survey: Rest-Frame Optical Spectroscopy for ~1500 H-Selected Galaxies at 1.37 < z < 3.8

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

preprint2015arXiv

The SFR-M* Relation and Empirical Star-Formation Histories from ZFOURGE at 0.5 < z < 4

We explore star-formation histories (SFHs) of galaxies based on the evolution of the star-formation rate stellar mass relation (SFR-M*). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M* relation at 0.5 < z < 4. Similar to recent works we find that the average infrared SEDs of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that the SFR-M* relation is not consistent with a single power-law of the form SFR ~ M*^a at any redshift; it has a power-law slope of a~1 at low masses, and becomes shallower above a turnover mass (M_0) that ranges from 10^9.5 - 10^10.8 Msol, with evidence that M_0 increases with redshift. We compare our measurements to results from state-of-the-art cosmological simulations, and find general agreement in the slope of the SFR-M* relation albeit with systematic offsets. We use the evolving SFR-M* sequence to generate SFHs, finding that typical SFRs of individual galaxies rise at early times and decline after reaching a peak. This peak occurs earlier for more massive galaxies. We integrate these SFHs to generate mass-growth histories and compare to the implied mass-growth from the evolution of the stellar mass function. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass-growth rates that are as much as 10x higher than inferred from the stellar mass function. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers.

preprint2015arXiv

The sizes of massive quiescent and star forming galaxies at z~4 with ZFOURGE and CANDELS

We study the rest-frame ultra-violet sizes of massive (~0.8 x 10^11 M_Sun) galaxies at 3.4<z<4.2, selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles to HST/WFC3/F160W images from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). Massive quiescent galaxies are very compact, with a median circularized half-light radius r_e = 0.63 +/- 0.18 kpc. Removing 5/16 (31%) sources with signs of AGN activity does not change the result. Star-forming galaxies have r_e = 2.0 +/- 0.60 kpc, 3.2 +/- 1.3 x larger than quiescent galaxies. Quiescent galaxies at z~4 are on average 6.0 +\- 0.17 x smaller than at z~0 and 1.9 +/- 0.7 x smaller than at z~2. Star-forming galaxies of the same stellar mass are 2.4 +/- 0.7 x smaller than at z~0. Overall, the size evolution at 0<z<4 is well described by a powerlaw, with r_e = 5.08 +/- 0.28 (1+z)^(-1.44+/-0.08) kpc for quiescent and r_e = 6.02 +/- 0.28 (1+z)^(-0.72+/-0.05) kpc for star-forming galaxies. Compact star-forming galaxies are rare in our sample: we find only 1/14 (7%) with r_e / (M / 10^11 M_Sun)^0.75 < 1.5, whereas 13/16 (81%) of the quiescent galaxies is compact. The number density of compact quiescent galaxies at z~4 is 1.8 +/- 0.8 x 10^-5 Mpc^-3 and increases rapidly, by >5 x, between 2<z<4. The paucity of compact star-forming galaxies at z~4 and their large rest-frame ultra-violet median sizes suggest that the formation phase of compact cores is very short and/or highly dust obscured.

preprint2015arXiv

Where stars form: inside-out growth and coherent star formation from HST Halpha maps of 2676 galaxies across the main sequence at z~1

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

preprint2014arXiv

3D-HST WFC3-selected Photometric Catalogs in the Five CANDELS/3D-HST Fields: Photometry, Photometric Redshifts and Stellar Masses

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

preprint2014arXiv

Bulge Growth and Quenching since z = 2.5 in CANDELS/3D-HST

Exploiting the deep high-resolution imaging of all 5 CANDELS fields, and accurate redshift information provided by 3D-HST, we investigate the relation between structure and stellar populations for a mass-selected sample of 6764 galaxies above 10^10 Msun, spanning the redshift range 0.5 < z < 2.5. For the first time, we fit 2-dimensional models comprising a single Sersic fit and two-component (i.e., bulge + disk) decompositions not only to the H-band light distributions, but also to the stellar mass maps reconstructed from resolved stellar population modeling. We confirm that the increased bulge prominence among quiescent galaxies, as reported previously based on rest-optical observations, remains in place when considering the distributions of stellar mass. Moreover, we observe an increase of the typical Sersic index and bulge-to-total ratio (with median B/T reaching 40-50%) among star-forming galaxies above 10^11 Msun. Given that quenching for these most massive systems is likely to be imminent, our findings suggest that significant bulge growth precedes a departure from the star-forming main sequence. We demonstrate that the bulge mass (and ideally knowledge of the bulge and total mass) is a more reliable predictor of the star-forming versus quiescent state of a galaxy than the total stellar mass. The same trends are predicted by the state-of-the-art semi-analytic model by Somerville et al. In the latter, bulges and black holes grow hand in hand through merging and/or disk instabilities, and AGN-feedback shuts off star formation. Further observations will be required to pin down star formation quenching mechanisms, but our results imply they must be internal to the galaxies and closely associated with bulge growth.

preprint2014arXiv

Constraining the Low-Mass Slope of the Star Formation Sequence at 0.5<z<2.5

We constrain the slope of the star formation rate ($\logΨ$) to stellar mass ($\log\mathrm{M_{\star}}$) relation down to $\log(\mathrm{M_{\star}/M_{\odot}})=8.4$ ($\log(\mathrm{M_{\star}/M_{\odot}})=9.2$) at $z=0.5$ ($z=2.5$) with a mass-complete sample of 39,106 star-forming galaxies selected from the 3D-HST photometric catalogs, using deep photometry in the CANDELS fields. For the first time, we find that the slope is dependent on stellar mass, such that it is steeper at low masses ($\log\mathrmΨ\propto\log\mathrm{M_{\star}}$) than at high masses ($\log\mathrmΨ\propto(0.3-0.6)\log\mathrm{M_{\star}}$). These steeper low mass slopes are found for three different star formation indicators: the combination of the ultraviolet (UV) and infrared (IR), calibrated from a stacking analysis of Spitzer/MIPS 24$μ$m imaging; $β$-corrected UV SFRs; and H$α$ SFRs. The normalization of the sequence evolves differently in distinct mass regimes as well: for galaxies less massive than $\log(\mathrm{M_{\star}/M_{\odot}})<10$ the specific SFR ($Ψ/\mathrm{M_{\star}}$) is observed to be roughly self-similar with $Ψ/\mathrm{M_{\star}}\propto(1+z)^{1.9}$, whereas more massive galaxies show a stronger evolution with $Ψ/\mathrm{M_{\star}}\propto(1+z)^{2.2-3.5}$ for $\log(\mathrm{M_{\star}/M_{\odot}})=10.2-11.2$. The fact that we find a steep slope of the star formation sequence for the lower mass galaxies will help reconcile theoretical galaxy formation models with the observations. The results of this study support the analytical conclusions of Leja et al. (2014).

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

Direct measurements of dust attenuation in z~1.5 star-forming galaxies from 3D-HST: Implications for dust geometry and star formation rates

The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct comparison between dust around star-forming regions ($A_{V,\mathrm{HII}}$) and the integrated dust content ($A_{V,\mathrm{star}}$). We select a sample of 163 galaxies between $1.36\le{}z\le1.5$ with H$α$ signal-to-noise ratio $\ge5$ and measure Balmer decrements from stacked spectra to calculate $A_{V,\mathrm{HII}}$. First, we stack spectra in bins of $A_{V,\mathrm{star}}$, and find that $A_{V,\mathrm{HII}}=1.86\,A_{V,\mathrm{star}}$, with a significance of $σ=1.7$. Our result is consistent with the two-component dust model, in which galaxies contain both diffuse and stellar birth cloud dust. Next, we stack spectra in bins of specific star formation rate ($\log\,\mathrm{SSFR}$), star formation rate ($\log\,\mathrm{SFR}$), and stellar mass ($\log{}M_*$). We find that on average $A_{V,\mathrm{HII}}$ increases with SFR and mass, but decreases with increasing SSFR. Interestingly, the data hint that the amount of extra attenuation decreases with increasing SSFR. This trend is expected from the two-component model, as the extra attenuation will increase once older stars outside the star-forming regions become more dominant in the galaxy spectrum. Finally, using Balmer decrements we derive dust-corrected H$α$ SFRs, and find that stellar population modeling produces incorrect SFRs if rapidly declining star formation histories are included in the explored parameter space.

preprint2014arXiv

How dead are dead galaxies? Mid-Infrared fluxes of quiescent galaxies at redshift 0.3 < z < 2.5: implications for star formation rates and dust heating

We investigate the star formation rates of quiescent galaxies at high redshift (0.3 < z < 2.5) using 3D-HST WFC3 grism spectroscopy and Spitzer mid-infrared data. We select quiescent galaxies on the basis of the widely used UVJ color-color criteria. Spectral energy distribution fitting (rest frame optical and near-IR) indicates very low star formation rates for quiescent galaxies (sSFR ~ 10^-12 yr^-1). However, SED fitting can miss star formation if it is hidden behind high dust obscuration and ionizing radiation is re-emitted in the mid-infrared. It is therefore fundamental to measure the dust-obscured SFRs with a mid-IR indicator. We stack the MIPS-24um images of quiescent objects in five redshift bins centered on z = 0.5, 0.9, 1.2, 1.7, 2.2 and perform aperture photometry. Including direct 24um detections, we find sSFR ~ 10^-11.9 * (1+z)^4 yr^-1. These values are higher than those indicated by SED fitting, but at each redshift they are 20-40 times lower than those of typical star forming galaxies. The true SFRs of quiescent galaxies might be even lower, as we show that the mid-IR fluxes can be due to processes unrelated to ongoing star formation, such as cirrus dust heated by old stellar populations and circumstellar dust. Our measurements show that star formation quenching is very efficient at every redshift. The measured SFR values are at z > 1.5 marginally consistent with the ones expected from gas recycling (assuming that mass loss from evolved stars refuels star formation) and well above that at lower redshifts.

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

Reconciling the observed star-forming sequence with the observed stellar mass function

We examine the connection between the observed star-forming sequence (SFR $\propto$ $M^α$) and the observed evolution of the stellar mass function between $0.2 < z < 2.5$. We find the star-forming sequence cannot have a slope $α$ $\lesssim$ 0.9 at all masses and redshifts, as this would result in a much higher number density at $10 < \log(\mathrm{M/M_{\odot}}) < 11$ by $z=1$ than is observed. We show that a transition in the slope of the star-forming sequence, such that $α=1$ at $\log(\mathrm{M/M_{\odot}})<10.5$ and $α=0.7-0.13z$ ({Whitaker} {et~al.} 2012) at $\log(\mathrm{M/M_{\odot}})>10.5$, greatly improves agreement with the evolution of the stellar mass function. We then derive a star-forming sequence which reproduces the evolution of the mass function by design. This star-forming sequence is also well-described by a broken-power law, with a shallow slope at high masses and a steep slope at low masses. At $z=2$, it is offset by $\sim$0.3 dex from the observed star-forming sequence, consistent with the mild disagreement between the cosmic SFR and recent observations of the growth of the stellar mass density. It is unclear whether this problem stems from errors in stellar mass estimates, errors in SFRs, or other effects. We show that a mass-dependent slope is also seen in other self-consistent models of galaxy evolution, including semi-analytical, hydrodynamical, and abundance-matching models. As part of the analysis, we demonstrate that neither mergers nor hidden low-mass quiescent galaxies are likely to reconcile the evolution of the mass function and the star-forming sequence. These results are supported by observations from {Whitaker} {et~al.} (2014).

preprint2014arXiv

Resolved Star Formation on Sub-galactic Scales in a Merger at z=1.7

We present a detailed analysis of Hubble Space Telescope (HST), Wide Field Camera 3 (WFC3) G141 grism spectroscopy for seven star-forming regions of the highly magnified lensed starburst galaxy RCSGA 032727-132609 at z=1.704. We measure the spatial variations of the extinction in RCS0327 through the observed H$γ$/H$β$ emission line ratios, finding a constant average extinction of $\mathrm{E(B-V)_{gas}}=0.40\pm0.07$. We infer that the star formation is enhanced as a result of an ongoing interaction, with measured star formation rates derived from demagnified, extinction-corrected H$β$ line fluxes for the individual star-forming clumps falling >1-2 dex above the star formation sequence. When combining the HST/WFC3 [OIII]$λ$5007/H$β$ emission line ratio measurements with [NII]/H$α$ line ratios from Wuyts et a. (2014), we find that the majority of the individual star-forming regions fall along the local "normal" abundance sequence. With the first detections of the He I $λ$5876 and He II $λ$4686 recombination lines in a distant galaxy, we probe the massive-star content of the star-forming regions in RCS0327. The majority of the star-forming regions have a He I $λ$5876 to H$β$ ratio consistent with the saturated maximum value, which is only possible if they still contain hot O-stars. Two regions have lower ratios, implying that their last burst of new star formation ended ~5 Myr ago. Together, the He I $λ$5876 and He II $λ$4686 to H$β$ line ratios provide indirect evidence for the order in which star formation is stopping in individual star-forming knots of this high redshift merger. We place the spatial variations of the extinction, star formation rate and ionization conditions in the context of the star formation history of RCS0327.

preprint2014arXiv

The Nature of Extreme Emission Line Galaxies at z=1-2: Kinematics and Metallicities from Near-Infrared Spectroscopy

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

preprint2013arXiv

A CANDELS - 3D-HST Synergy: Resolved Star Formation Patterns at 0.7 < z < 1.5

We analyze the resolved stellar populations of 473 massive star-forming galaxies at 0.7 < z < 1.5, with multi-wavelength broad-band imaging from CANDELS and Halpha surface brightness profiles at the same kiloparsec resolution from 3D-HST. Together, this unique data set sheds light on how the assembled stellar mass is distributed within galaxies, and where new stars are being formed. We find the Halpha morphologies to resemble more closely those observed in the ACS I band than in the WFC3 H band, especially for the larger systems. We next derive a novel prescription for Halpha dust corrections, which accounts for extra extinction towards HII regions. The prescription leads to consistent SFR estimates and reproduces the observed relation between the Halpha/UV luminosity ratio and visual extinction, both on a pixel-by-pixel and on a galaxy-integrated level. We find the surface density of star formation to correlate with the surface density of assembled stellar mass for spatially resolved regions within galaxies, akin to the so-called 'main sequence of star formation' established on a galaxy-integrated level. Deviations from this relation towards lower equivalent widths are found in the inner regions of galaxies. Clumps and spiral features, on the other hand, are associated with enhanced Halpha equivalent widths, bluer colors, and higher specific star formation rates compared to the underlying disk. Their Halpha/UV luminosity ratio is lower than that of the underlying disk, suggesting the ACS clump selection preferentially picks up those regions of elevated star formation activity that are the least obscured by dust. Our analysis emphasizes that monochromatic studies of galaxy structure can be severely limited by mass-to-light ratio variations due to dust and spatially inhomogeneous star formation histories.

preprint2013arXiv

Confirmation of Small Dynamical and Stellar Masses for Extreme Emission Line Galaxies at z~2

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

preprint2013arXiv

Exploring the chemical link between local ellipticals and their high-redshift progenitors

We present Keck/MOSFIRE K-band spectroscopy of the first mass-selected sample of galaxies at $z\sim2.3$. Targets are selected from the 3D-HST Treasury survey. The six detected galaxies have a mean [NII]$λ$6584/H$α$ ratio of $0.27\pm0.01$, with a small standard deviation of 0.05. This mean value is similar to that of UV-selected galaxies of the same mass. The mean gas-phase oxygen abundance inferred from the [NII]/H$α$ ratios depends on the calibration method, and ranges from 12+log(O/H)$_{gas}=8.57$ for the {Pettini} & {Pagel} (2004) calibration to 12+log(O/H)$_{gas}= 8.87$ for the {Maiolino} {et~al.} (2008) calibration. Measurements of the stellar oxygen abundance in nearby quiescent galaxies with the same number density indicate 12+log(O/H)$_{stars}= 8.95$, similar to the gas-phase abundances of the $z\sim2.3$ galaxies if the {Maiolino} {et~al.} (2008) calibration is used. This suggests that these high-redshift star forming galaxies may be progenitors of today's massive early-type galaxies. The main uncertainties are the absolute calibration of the gas-phase oxygen abundance and the incompleteness of the $z\sim2.3$ sample: the galaxies with detected H$α$ tend to be larger and have higher star formation rates than the galaxies without detected H$α$, and we may still be missing the most dust-obscured progenitors.

preprint2013arXiv

Galaxy environments over cosmic time: the non-evolving radial galaxy distributions around massive galaxies since z=1.6

We present a statistical study of the environments of massive galaxies in four redshift bins between z=0.04 and z=1.6, using data from the Sloan Digital Sky Survey (SDSS) and the NEWFIRM Medium Band Survey (NMBS). We measure the projected radial distribution of galaxies in cylinders around a constant number density selected sample of massive galaxies and utilize a statistical subtraction of contaminating sources. Our analysis shows that massive primary galaxies typically live in group halos and are surrounded by 2 to 3 satellites with masses more than one-tenth of the primary galaxy mass. The cumulative stellar mass in these satellites roughly equals the mass of the primary galaxy itself. We further find that the radial number density profile of galaxies around massive primaries has not evolved significantly in either slope or overall normalization in the past 9.5 Gyr. A simplistic interpretation of this result can be taken as evidence for a lack of mergers in the studied groups and as support for a static evolution model of halos containing massive primaries. Alternatively, there exists a tight balance between mergers and accretion of new satellites such that the overall distribution of galaxies in and around the halo is preserved. The latter interpretation is supported by a comparison to a semi-analytic model, which shows a similar constant average satellite distribution over the same redshift range.

preprint2013arXiv

Quiescent Galaxies in the 3D-HST Survey: Spectroscopic Confirmation of a Large Number of Galaxies with Relatively Old Stellar Populations at z~2

Quiescent galaxies at z~2 have been identified in large numbers based on rest-frame colors, but only a small number of these galaxies have been spectroscopically confirmed to show that their rest-frame optical spectra show either strong Balmer or metal absorption lines. Here, we median stack the rest-frame optical spectra for 171 photometrically-quiescent galaxies at 1.4 < z < 2.2 from the 3D-HST grism survey. In addition to Hbeta (4861A), we unambiguously identify metal absorption lines in the stacked spectrum, including the G-band (4304A), Mg I (5175A), and Na I (5894A). This finding demonstrates that galaxies with relatively old stellar populations already existed when the universe was ~3 Gyr old, and that rest-frame color selection techniques can efficiently select them. We find an average age of 1.3^0.1_0.3 Gyr when fitting a simple stellar population to the entire stack. We confirm our previous result from medium-band photometry that the stellar age varies with the colors of quiescent galaxies: the reddest 80% of galaxies are dominated by metal lines and have a relatively old mean age of 1.6^0.5_0.4 Gyr, whereas the bluest (and brightest) galaxies have strong Balmer lines and a spectroscopic age of 0.9^0.2_0.1 Gyr. Although the spectrum is dominated by an evolved stellar population, we also find [OIII] and Hbeta emission. Interestingly, this emission is more centrally concentrated than the continuum with L_[OIII] = 1.7 +/- 0.3 x 10^40 erg s^-1, indicating residual central star formation or nuclear activity.

preprint2013arXiv

The Assembly of Milky Way-like Galaxies Since z~2.5

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

preprint2013arXiv

The Spatial Extent and Distribution of Star Formation in 3D-HST Mergers at z~1.5

We present an analysis of the spatial distribution of star formation in a sample of 60 visually identified galaxy merger candidates at z>1. Our sample, drawn from the 3D-HST survey, is flux-limited and was selected to have high star formation rates based on fits of their broad-band, low spatial resolution spectral energy distributions. It includes plausible pre-merger (close pairs) and post-merger (single objects with tidal features) systems, with total stellar masses and star formation rates derived from multi-wavelength photometry. Here we use near-infrared slitless spectra from 3D-HST which produce Halpha or [OIII] emission line maps as proxies for star-formation maps. This provides a first comprehensive high-resolution, empirical picture of where star formation occurred in galaxy mergers at the epoch of peak cosmic star formation rate. We find that detectable star formation can occur in one or both galaxy centres, or in tidal tails. The most common case (58%) is that star formation is largely concentrated in a single, compact region, coincident with the centre of (one of) the merger components. No correlations between star formation morphology and redshift, total stellar mass, or star formation rate are found. A restricted set of hydrodynamical merger simulations between similarly massive and gas-rich objects implies that star formation should be detectable in both merger components, when the gas fractions of the individual components are the same. This suggests that z~1.5 mergers typically occur between galaxies whose gas fractions, masses, and/or star formation rates are distinctly different from one another.

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.

preprint2013arXiv

What Are the Progenitors of Compact, Massive, Quiescent Galaxies at z=2.3? The Population of Massive Galaxies at z>3 from NMBS and CANDELS

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

preprint2013arXiv

X-ray properties of K-selected galaxies at 0.5<z<2.0: Investigating trends with stellar mass, redshift and spectral type

We examine how the total X-ray luminosity correlates with stellar mass, stellar population, and redshift for a K-band limited sample of ~3500 galaxies at 0.5<z<2.0 from the NEWFIRM Medium Band Survey in the COSMOS field. The galaxy sample is divided into 32 different galaxy types, based on similarities between the spectral energy distributions. For each galaxy type, we further divide the sample into bins of redshift and stellar mass, and perform an X-ray stacking analysis using the Chandra COSMOS (C-COSMOS) data. We find that full band X-ray luminosity is primarily increasing with stellar mass, and at similar mass and spectral type is higher at larger redshifts. When comparing at the same stellar mass, we find that the X-ray luminosity is slightly higher for younger galaxies (i.e., weaker 4000Åbreaks), but the scatter in this relation is large. We compare the observed X-ray luminosities to those expected from low and high mass X-ray binaries (XRBs). For blue galaxies, XRBs can almost fully account for the observed emission, while for older galaxies with larger 4000Åbreaks, active galactic nuclei (AGN) or hot gas dominate the measured X-ray flux. After correcting for XRBs, the X-ray luminosity is still slightly higher in younger galaxies, although this correlation is not significant. AGN appear to be a larger component of galaxy X-ray luminosity at earlier times, as the hardness ratio increases with redshift. Together with the slight increase in X-ray luminosity this may indicate more obscured AGNs or higher accretion rates at earlier times.

preprint2012arXiv

3D-HST Grism Spectroscopy of a Gravitationally Lensed, Low-metallicity Starburst Galaxy at z=1.847

We present Hubble Space Telescope (HST) imaging and spectroscopy of the gravitational lens SL2SJ02176-0513, a cusp arc at z=1.847. The UV continuum of the lensed galaxy is very blue, which is seemingly at odds with its redder optical colors. The 3D-HST WFC3/G141 near-infrared spectrum of the lens reveals the source of this discrepancy to be extremely strong [OIII]5007 and H-beta emission lines with rest-frame equivalent widths of 2000 +/- 100 and 520 +/- 40 Angstroms, respectively. The source has a stellar mass ~10^8 Msun, sSFR \sim 100/Gyr, and detection of [OIII]4363 yields a metallicity of 12 + log(O/H) = 7.5 +/- 0.2. We identify local blue compact dwarf analogs to SL2SJ02176-0513, which are among the most metal-poor galaxies in the SDSS. The local analogs resemble the lensed galaxy in many ways, including UV/optical SED, spatial morphology and emission line equivalent widths and ratios. Common to SL2SJ02176-0513 and its local counterparts is an upturn at mid-IR wavelengths likely arising from hot dust heated by starbursts. The emission lines of SL2SJ02176-0513 are spatially resolved owing to the combination of the lens and the high spatial resolution of HST. The lensed galaxy is composed of two clumps with combined size r_e \sim 300 pc, and we resolve significant differences in UV color and emission line equivalent width between them. Though it has characteristics occasionally attributed to active galactic nuclei, we conclude that SL2SJ02176-0513 is a low-metallicity star-bursting dwarf galaxy. Such galaxies will be found in significant numbers in the full 3D-HST grism survey.

preprint2012arXiv

Halpha Equivalent Widths from the 3D-HST survey: evolution with redshift and dependence on stellar mass

We investigate the evolution of the Halpha equivalent width, EW(Halpha), with redshift and its dependence on stellar mass, taking advantage of the first data from the 3D-HST survey, a large spectroscopic Treasury program with the Hubble Space Telescope WFC3. Combining our Halpha measurements of 854 galaxies at 0.8<z<1.5 with those of ground based surveys at lower and higher redshift, we can consistently determine the evolution of the EW(Halpha) distribution from z=0 to z=2.2. We find that at all masses the characteristic EW(Halpha) is decreasing towards the present epoch, and that at each redshift the EW(Halpha) is lower for high-mass galaxies. We measure a slope of EW(Halpha) ~ (1+z)^(1.8) with little mass dependence. Qualitatively, this measurement is a model-independent confirmation of the evolution of star forming galaxies with redshift. A quantitative conversion of EW(Halpha) to sSFR is very model dependent, because of differential reddening corrections between the continuum SED and the Balmer lines. The observed EW(Halpha) can be reproduced with a simple model in which the SFR for galaxies rises to the epoch of z~2.5 and then decreases with time to z = 0. The model implies that the EW(Halpha) rises to 400 A at z=8. The sSFR evolves faster than EW(Halpha), as the mass-to-light ratio also evolves with redshift. In this context, we find that the sSFR evolves as (1+z)^(3.2), nearly independent of mass, consistent with previous reddening insensitive estimates. We confirm previous results that the observed slope of the sSFR-z relation is steeper than the one predicted by models, but models and observations agree in finding little mass dependence.

preprint2012arXiv

The Evolution of the Rest-frame V-band Luminosity Function from z=4: A Constant Faint-end Slope over the Last 12 Gyr of Cosmic History

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

preprint2012arXiv

The Star-formation Mass Sequence out to z=2.5

We study the star formation rate (SFR) - stellar mass (M*) relation in a self-consistent manner from 0 < z < 2.5 with a sample of galaxies selected from the NEWFIRM Medium-Band Survey. We find a significant non-linear slope of the relation, SFR \propto M*^0.6, and a constant observed scatter of 0.34 dex, independent of redshift and M*. However, if we select only blue galaxies we find a linear relation SFR \propto M*, similar to previous results at z = 0 by Peng et al. (2010). This selection excludes red, dusty, star-forming galaxies with higher masses, which brings down the slope. By selecting on L_IR/L_UV (a proxy for dust obscuration) and the rest-frame U-V colors, we show that star-forming galaxies fall in three distinct regions of the log(SFR)-log(M*) plane: 1) actively star-forming galaxies with "normal" dust obscuration and associated colors (54% for log(M*) > 10 at 1 < z < 1.5), 2) red star-forming galaxies with low levels of dust obscuration and low specific SFRs (11%), and 3) dusty, blue star-forming galaxies with high specific SFRs (7%). The remaining 28% comprises quiescent galaxies. Galaxies on the "normal" star formation sequence show strong trends of increasing dust attenuation with stellar mass and a decreasing specific SFR, with an observed scatter of 0.25 dex (0.17 dex intrinsic scatter). The dusty, blue galaxies reside in the upper envelope of the star formation sequence with remarkably similar spectral shapes at all masses, suggesting that the same physical process is dominating the stellar light. The red, low-dust star-forming galaxies may be in the process of shutting off and migrating to the quiescent population.

preprint2011arXiv

A Large Population of Massive Compact Post-Starburst Galaxies at z>1: Implications for the Size Evolution and Quenching Mechanism of Quiescent Galaxies

We study the growth of the red sequence through the number density and structural evolution of a sample of young and old quiescent galaxies at 0<z<2. The galaxies are selected from the NEWFIRM Medium-Band Survey (NMBS) in the Cosmic Evolution Survey (COSMOS) field. We find a large population of massive young recently quenched ("post-starburst") galaxies at z>1 that are almost non-existent at z<1; their number density is 5 x 10^{-5} Mpc^{-3} at z=2, whereas it is a factor of 10 less at z=0.5. The observed number densities of young and old quiescent galaxies at z>1 are consistent with a simple model in which all old quiescent galaxies were once identified as post-starburst galaxies. We find that the overall population of quiescent galaxies have smaller sizes and slightly more elongated shapes at higher redshift, in agreement with other recent studies. Interestingly, the most recently quenched galaxies at 1<z<2 are not larger, and possibly even smaller, than older galaxies at those redshifts. This result is inconsistent with the idea that the evolution of the average size of quiescent galaxies is largely driven by continuous transformations of larger, star-forming galaxies: in that case, the youngest quiescent galaxies would also be the largest. Instead, mergers or other mechanisms appear to be required to explain the size growth of quiescent galaxies from z=2 to the present.

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

Halpha and 4000 Angstrom Break Measurements for ~3500 K-selected Galaxies at 0.5<z<2.0

We measure spectral features of ~3500 K-selected galaxies at 0.5<z<2.0 from high quality medium-band photometry using a new technique. First, we divide the galaxy sample in 32 subsamples based on the similarities between the full spectral energy distributions (SEDs) of the galaxies. For each of these 32 galaxy types we construct a composite SED by de-redshifting and scaling the observed photometry. This approach increases the signal-to-noise ratio and sampling of galaxy SEDs and allows for model-independent stellar population studies. The composite SEDs are of spectroscopic quality, and facilitate -- for the first time -- Halpha measurement for a large magnitude-limited sample of distant galaxies. The linewidths indicate a photometric redshift uncertainty of dz<0.02x(1+z). The composite SEDs also show the Balmer and 4000 Angstrom breaks, MgII absorption at ~2800 Angstrom, the dust absorption feature at 2175 Angstrom, and blended [OIII]+Hbeta emission. We compare the total equivalent width of Halpha, [NII], and [SII] (W_Halpha+) with the strength of the 4000 Angstrom break (D(4000)) and the best-fit specific star formation rate, and find that all these properties are strongly correlated. This is a reassuring result, as currently most distant stellar population studies are based on just continuum emission. Furthermore, the relation between W_Halpha+ and Dn(4000) provides interesting clues to the SFHs of galaxies, as these features are sensitive to different stellar ages. We find that the correlation between W_Halpha+ and D(4000) at 0.5<z<2.0 is similar to z~0, and that the suppression of star formation in galaxies at z<2 is generally not abrupt, but a gradual process.

preprint2011arXiv

Redshift Evolution of the Galaxy Velocity Dispersion Function

We present a study of the evolution of the galaxy Velocity Dispersion Function (VDF) from z=0 to z=1.5 using photometric data from the UKIDSS Ultra Deep Survey (UDS) and Newfirm Medium Band Survey (NMBS) COSMOS surveys. The VDF has been measured locally using direct kinematic measurements from the Sloan Digital Sky Survey, but direct studies of the VDF at high redshift are difficult as they require velocity dispersion measurements of many thousands of galaxies. Taylor et al. (2010) demonstrated that dynamical and stellar mass are linearly related when the structure of the galaxy is accounted for. We show that the stellar mass, size and Sérsic index can reliably predict the velocity dispersions of SDSS galaxies. We apply this relation to galaxies at high redshift and determine the evolution of the inferred VDF. We find that the VDF at z~0.5 is very similar to the VDF at z=0. At higher redshifts, we find that the number density of galaxies with dispersions <~200 km/s is lower, but the number of high dispersion galaxies is constant or even higher. At fixed cumulative number density, the velocity dispersions of galaxies with logN[Mpc^{-3}]<-3.5 increase with time by a factor of ~1.4 from z~1.5-0, whereas the dispersions of galaxies with lower number density are approximately constant or decrease with time. The VDF appears to show less evolution than the stellar mass function, particularly at the lowest number densities. We note that these results are still somewhat uncertain and we suggest several avenues for further calibrating the inferred velocity dispersions.

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

preprint2011arXiv

The stellar velocity dispersion of a compact massive galaxy at z=1.80 using X-Shooter: confirmation of the evolution in the mass-size and mass-dispersion relations

Recent photometric studies have shown that early-type galaxies at fixed stellar mass were smaller and denser at earlier times. In this paper we assess that finding by deriving the dynamical mass of such a compact quiescent galaxy at z=1.8. We have obtained a high-quality spectrum with full UV-NIR wavelength coverage of galaxy NMBS-C7447 using X-Shooter on the VLT. We determined a velocity dispersion of 294 +- 51 km/s. Given this velocity dispersion and the effective radius of 1.64 +- 0.15 kpc (as determined from HST-WFC3 F160W observations) we derive a dynamical mass of 1.7 +- 0.5 x 10^11 Msun. Comparison of the full spectrum with stellar population synthesis models indicates that NMBS-C774 has a relatively young stellar population (0.40 Gyr) with little or no star formation and a stellar mass of ~1.5 x 10^11 Msun. The dynamical and photometric stellar mass are in good agreement. Thus, our study supports the conclusion that the mass densities of quiescent galaxies were indeed higher at earlier times, and this earlier result is not caused by systematic measurement errors. By combining available spectroscopic measurements at different redshifts, we find that the velocity dispersion at fixed dynamical mass was a factor of ~1.8 higher at z=1.8 compared to z=0. Finally, we show that the apparent discrepancies between the few available velocity dispersion measurements at z>1.5 are consistent with the intrinsic scatter of the mass-size relation.

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

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.

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.