Researcher profile

Pieter G. van Dokkum

Pieter G. van Dokkum contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
78works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

78 published item(s)

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.

preprint2020arXiv

Brackett-$γ$ as a Gold-standard Test of Star Formation Rates Derived from SED Fitting

Using a local reference sample of 21 galaxies, we compare observations of the $λ$2.16 $μ$m Brackett-$γ$ (Br$γ$) hydrogen recombination line with predictions from the Prospector Bayesian inference framework, which was used to fit the broadband photometry of these systems. This is a clean test of the spectral-energy-distribution-derived star formation rates (SFRs), as dust is expected to be optically thin at this wavelength in nearly all galaxies; thus, the internal conversion of SFR to predicted line luminosity does not depend strongly on the adopted dust model and posterior dust parameters, as is the case for shorter-wavelength lines such as H$α$. We find that Prospector predicts Br$γ$ luminosities and equivalent widths with small offsets ($\sim$0.05 dex), and scatter ($\sim$0.2 dex), consistent with measurement uncertainties, though we caution that the derived offset is dependent on the choice of stellar isochrones. We demonstrate that even when the Prospector-derived dust attenuation does not well describe, e.g., H$α$ line properties or observed reddening between H$α$ and Br$γ$, the underlying SFRs are accurate, as verified by the dust-free Br$γ$ comparison. Finally, we discuss in what ways Br$γ$ might be able to help constrain model parameters when treated as an input to the model, and comment on its potential as an accurate monochromatic SFR indicator in the era of JWST multiobject near-IR spectroscopy.

preprint2016arXiv

A massive, quiescent, population II galaxy at a redshift of 2.1

Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.

preprint2016arXiv

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

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

preprint2016arXiv

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

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

ZFOURGE catalogue of AGN candidates: an enhancement of 160μm-derived star-formation rates in active galaxies to $z$ = 3.2

We investigate active galactic nuclei (AGN) candidates within the FourStar Galaxy Evolution Survey (ZFOURGE) to determine the impact they have on star-formation in their host galaxies. We first identify a population of radio, X-ray, and infrared-selected AGN by cross-matching the deep $K_{s}$-band imaging of ZFOURGE with overlapping multi-wavelength data. From this, we construct a mass-complete (log(M$_{*}$/M$_{\odot}$) $\ge$ 9.75), AGN luminosity limited sample of 235 AGN hosts over z = 0.2 - 3.2. We compare the rest-frame U - V versus V - J (UVJ) colours and specific star-formation rates (sSFRs) of the AGN hosts to a mass-matched control sample of inactive (non-AGN) galaxies. UVJ diagnostics reveal AGN tend to be hosted in a lower fraction of quiescent galaxies and a higher fraction of dusty galaxies than the control sample. Using 160μm Herschel PACS data, we find the mean specific star-formation rate of AGN hosts to be elevated by 0.34$\pm$0.07 dex with respect to the control sample across all redshifts. This offset is primarily driven by infrared-selected AGN, where the mean sSFR is found to be elevated by as much as a factor of ~5. The remaining population, comprised predominantly of X-ray AGN hosts, is found mostly consistent with inactive galaxies, exhibiting only a marginal elevation. We discuss scenarios that may explain these findings and postulate that AGN are less likely to be a dominant mechanism for moderating galaxy growth via quenching than has previously been suggested.

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 &#34;normal&#34; 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 &#34;star-formation sequence&#34;, 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

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

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

preprint2015arXiv

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

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

preprint2015arXiv

The 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

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&#39; 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

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

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

On the importance of using appropriate spectral models to derive physical properties of galaxies at 0.7<z<2.8

Interpreting observations of distant galaxies in terms of constraints on physical parameters - such as stellar mass, star-formation rate (SFR) and dust optical depth - requires spectral synthesis modelling. We analyse the reliability of these physical parameters as determined under commonly adopted `classical&#39; assumptions: star-formation histories assumed to be exponentially declining functions of time, a simple dust law and no emission-line contribution. Improved modelling techniques and data quality now allow us to use a more sophisticated approach, including realistic star-formation histories, combined with modern prescriptions for dust attenuation and nebular emission (Pacifici et al. 2012). We present a Bayesian analysis of the spectra and multi-wavelength photometry of 1048 galaxies from the 3D-HST survey in the redshift range 0.7<z<2.8 and in the stellar mass range 9<log(M/Mo)<12. We find that, using the classical spectral library, stellar masses are systematically overestimated (~0.1 dex) and SFRs are systematically underestimated (~0.6 dex) relative to our more sophisticated approach. We also find that the simultaneous fit of photometric fluxes and emission-line equivalent widths helps break a degeneracy between SFR and optical depth of the dust, reducing the uncertainties on these parameters. Finally, we show how the biases of classical approaches can affect the correlation between stellar mass and SFR for star-forming galaxies (the `Star-Formation Main Sequence&#39;). We conclude that the normalization, slope and scatter of this relation strongly depend on the adopted approach and demonstrate that the classical, oversimplified approach cannot recover the true distribution of stellar mass and SFR.

preprint2014arXiv

The Assembly Histories of Quiescent Galaxies Since z=0.7 from Absorption Line Spectroscopy

We present results from modeling the optical spectra of a large sample of quiescent galaxies between 0.1 < z < 0.7 from the Sloan Digital Sky Survey (SDSS) and the AGN and Galaxy Evolution Survey (AGES). We examine how the stellar ages and abundance patterns of galaxies evolve over time as a function of stellar mass from 10^{9.6}-10^{11.8} Msun. Galaxy spectra are stacked in bins of mass and redshift, and modeled over a wavelength range from 4000 A to 5500 A. Full spectrum stellar population synthesis modeling provides estimates of the age and the abundances of the elements Fe, Mg, C, N, and Ca. We find negligible evolution in elemental abundances at fixed stellar mass over roughly 7 Gyr of cosmic time. In addition, the increase in stellar ages with time for massive galaxies is consistent with passive evolution since z = 0.7. Taken together, these results favor a scenario in which the inner ~ 0.3-3 R_e of massive quiescent galaxies have been passively evolving over the last half of cosmic time. Interestingly, the derived stellar ages are considerably younger than the age of the universe at all epochs, consistent with an equivalent single-burst star formation epoch of z < 1.5. These young stellar population ages coupled with the existence of massive quiescent galaxies at z > 1 indicate the inhomogeneous nature of the z < 0.7 quiescent population. The data also permit the addition of newly-quenched galaxies at masses below ~10^{10.5} Msun at z < 0.7. Additionally, we analyze very deep Keck DEIMOS spectra of the two brightest quiescent galaxies in a cluster at z = 0.83. There is tentative evidence that these galaxies are older than their counterparts in low-density environments. In the Appendix, we demonstrate that our full spectrum modeling technique allows for accurate and reliable modeling of galaxy spectra to low S/N (~20 A^{-1}) and/or low spectral resolution (R ~ 500).

preprint2014arXiv

The Fundamental Plane of massive quiescent galaxies out to z~2

The Fundamental Plane (FP) of early-type galaxies, relating the effective radius, velocity dispersion, and surface brightness, has long been recognized as a unique tool for analyzing galaxy structure and evolution. With the discovery of distant quiescent galaxies and the introduction of high sensitivity near-infrared spectrographs, it is now possible to explore the FP out to z~2. In this Letter we study the evolution of the FP out to z~2 using kinematic measurements of massive quiescent galaxies ($M_{*}>10^{11} M_{\odot}$). We find preliminary evidence for the existence of an FP out to z~2. The scatter of the FP, however, increases from z~0 to z~2, even when taking into account the larger measurement uncertainties at higher redshifts. We find a strong evolution of the zero point from z~2 to z~0: $Δ\log_{10}M/L_g\propto(-0.49\pm0.03)~z$. In order to assess whether our spectroscopic sample is representative of the early-type galaxy population at all redshifts, we compare their rest-frame g-z colors with those from a larger mass complete sample of quiescent galaxies. At z>1 we find that the spectroscopic sample is bluer. We use the color offsets to estimate a mass-to-light ratio (M/L) correction. The implied FP zero point evolution after correction is significantly smaller: $Δ\log_{10}M/L_g\propto(-0.39\pm 0.02)~z$. This is consistent with an apparent formation redshift of $z_{\rm{form}}=6.62^{+3.19}_{-1.44}$ for the underlying population, ignoring the effects of progenitor bias. A more complete spectroscopic sample is required at z~2 to properly measure the M/L evolution from the FP evolution.

preprint2014arXiv

The Relation between Dynamical Mass-to-Light Ratio and Color for Massive Quiescent Galaxies out to z~2 and Comparison with Stellar Population Synthesis Models

We explore the relation between the dynamical mass-to-light ratio ($M/L$) and rest-frame color of massive quiescent galaxies out to z~2. We use a galaxy sample with measured stellar velocity dispersions in combination with Hubble Space Telescope and ground-based multi-band photometry. Our sample spans a large range in $\log M_{dyn}/L_{g}$ (of 1.6~dex) and $\log M_{dyn}/L_{K}$ (of 1.3~dex). There is a strong, approximately linear correlation between the $M/L$ for different wavebands and rest-frame color. The root-mean-scatter scatter in $\log~M_{dyn}/L$ residuals implies that it is possible to estimate the $M/L$ with an accuracy of ~0.25 dex from a single rest-frame optical color. Stellar population synthesis (SPS) models with a Salpeter stellar initial mass function (IMF) can not simultaneously match $M_{dyn}/L_{g}$ vs. $(g-z)_{rest-frame}$ and $M_{dyn}/L_{K}$ vs. $(g-K)_{rest-frame}$. By changing the slope of the IMF we are still unable to explain the M/L of the bluest and reddest galaxies. We find that an IMF with a slope between $α=2.35$ and $α=1.35$ provides the best match. We also explore a broken IMF with a Salpeter slope at $M<1M_{\odot}$ and $M>4M_{\odot}$ and a slope $α$ in the intermediate region. The data favor a slope of $α=1.35$ over $α=2.35$. Nonetheless, our results show that variations between different SPS models are comparable to the IMF variations. In our analysis we assume that the variation in $M/L$ and color is driven by differences in age, and that other contributions (e.g., metallicity evolution, dark matter) are small. These assumptions may be an important source of uncertainty as galaxies evolve in more complex ways.

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 &#39;main sequence of star formation&#39; 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

A Tentative Detection of an Emission Line at 1.6 μm for the z~12 Candidate UDFj-39546284

We present deep WFC3 grism observations of the candidate z~12 galaxy UDFj-39546284 in the HST Ultra Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total exposure time is 40.5 ks and the spectrum covers 1.10 < λ < 1.65μm. We search for faint emission lines by cross-correlating the 2D G141 spectrum with the observed H160 morphology, a technique that is unique to slitless spectroscopy at HST resolution. We find a 2.7σ detection of an emission line at 1.599 μm---just redward of the JH140 filter---with flux 3.5 +/- 1.3x10-18 erg/s/cm2. Assuming the line is real, it contributes 110 +/- 40% of the observed H160 flux and has an observed equivalent width > 7300 Å. If the line is confirmed, it could be Ly-α at z = 12.12. However, a more plausible interpretation, given current results, could be a lower redshift feature such as [O III]λ4959,5007 at z = 2.19. We find two other 3D-HST [O III] emitters within 1000 km/s of that redshift in the GOODS-South field. Additional support for this interpretation comes from the discovery of a bright &#34;[O III] blob&#34; with a secure G141 grism redshift of z = 1.605. This object has a strikingly large observed equivalent width of nearly 9000 Å that results in similar &#34;dropout&#34; colors as UDFj-39546284.

preprint2013arXiv

Dynamical Versus Stellar Masses in Compact Early-Type Galaxies: Further Evidence for Systematic Variation in the Stellar Initial Mass Function

Several independent lines of evidence suggest that the stellar initial mass function (IMF) in early-type galaxies becomes increasingly `bottom-heavy&#39; with increasing galaxy mass and/or velocity dispersion, sigma. Here we consider evidence for IMF variation in a sample of relatively compact early-type galaxies drawn from the Sloan Digital Sky Survey. These galaxies are of sufficiently high stellar density that a dark halo likely makes a minor contribution to the total dynamical mass, Mdyn, within one effective radius. We fit our detailed stellar population synthesis models to the stacked absorption line spectra of these galaxies in bins of sigma and find evidence from IMF-sensitive spectral features for a bottom-heavy IMF at high sigma. We also apply simple `mass-follows-light&#39; dynamical models to the same data and find that Mdyn is significantly higher than what would be expected if these galaxies were stellar dominated and had a universal Milky Way IMF. Adopting Mdyn M_* therefore implies that the IMF is `heavier&#39; at high sigma. Most importantly, the quantitative amount of inferred IMF variation is very similar between the two techniques, agreeing to within 0.1 dex in mass. The agreement between two independent techniques, when applied to the same data, provides compelling evidence for systematic variation in the IMF as a function of early-type galaxy velocity dispersion. Any alternative explanations must reproduce both the results from dynamical and stellar population-based techniques.

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&#39;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

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

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

preprint2013arXiv

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

Stellar Kinematics of z~2 Galaxies and the Inside-Out Growth of Quiescent Galaxies

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

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

Tight Correlations Between Massive Galaxy Structural Properties and Dynamics: The Mass Fundamental Plane Was in Place by z~2

The Fundamental Plane (FP) is an empirical relation between the size, surface brightness, and velocity dispersion of early-type galaxies. This relation has been studied extensively for early-type galaxies in the local universe to constrain galaxy formation mechanisms. The evolution of the zeropoint of this plane has been extended to high redshifts to study the luminosity evolution of massive galaxies, under the assumption of structural homology. In this work, we assess this assumption by replacing surface brightness with stellar mass density and present the evolution of the &#34;mass FP&#34; for massive, quiescent galaxies since z~2. By accounting for stellar populations, we thereby isolate and trace structural and dynamical evolution. Despite the observed dramatic evolution in the sizes and morphologies of massive galaxies since z~3, we find that quiescent galaxies lie on the mass FP out to z~2. In contrast with ~1.4 dex evolution in the luminosity FP, average residuals from the z~0 mass FP are less than ~0.15 dex since z~2. Assuming the Hyde & Bernardi (2009) mass FP slope, we find that this minimal offset scales as (1+z)^{-0.095+/-0.043}. This result lends credence to previous studies that derived luminosity evolution from the FP. Therefore, despite their compact sizes and suggestions that massive galaxies are more disk-like at z~2, the relationship between their dynamics and structural properties are consistent with local early-type galaxies. Finally, we find no strong evidence for a tilt of the mass FP relative to the Virial plane, but emphasize the need for full models including selection biases to fully investigate this issue.

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

Mass growth and mergers: direct observations of the luminosity function of LRG satellite galaxies out to z=0.7 from SDSS and BOSS images

We present a statistical study of the luminosity functions of galaxies surrounding luminous red galaxies (LRGs) at average redshifts <z>=0.34 and <z>=0.65. The luminosity functions are derived by extracting source photometry around more than 40,000 LRGs and subtracting foreground and background contamination using randomly selected control fields. We show that at both studied redshifts the average luminosity functions of the LRGs and their satellite galaxies are poorly fitted by a Schechter function due to a luminosity gap between the centrals and their most luminous satellites. We utilize a two-component fit of a Schechter function plus a log-normal distribution to demonstrate that LRGs are typically brighter than their most luminous satellite by roughly 1.3 magnitudes. This luminosity gap implies that interactions within LRG environments are typically restricted to minor mergers with mass ratios of 1:4 or lower. The luminosity functions further imply that roughly 35% of the mass in the environment is locked in the LRG itself, supporting the idea that mass growth through major mergers within the environment is unlikely. Lastly, we show that the luminosity gap may be at least partially explained by the selection of LRGs as the gap can be reproduced by sparsely sampling a Schechter function. In that case LRGs may represent only a small fraction of central galaxies in similar mass halos.

preprint2012arXiv

Observations of dark and luminous matter: the radial distribution of satellite galaxies around massive red galaxies

We study the projected radial distribution of satellite galaxies around more than 28,000 Luminous Red Galaxies (LRGs) at 0.28<z<0.40 and trace the gravitational potential of LRG groups in the range 15<r/kpc<700. We show that at large radii the satellite number density profile is well fitted by a projected NFW profile with r_s~270 kpc and that at small radii this model underestimates the number of satellite galaxies. Utilizing the previously measured stellar light distribution of LRGs from deep imaging stacks we demonstrate that this small scale excess is consistent with a non-negligible baryonic mass contribution to the gravitational potential of massive groups and clusters. The combined NFW+scaled stellar profile provides an excellent fit to the satellite number density profile all the way from 15 kpc to 700 kpc. Dark matter dominates the total mass profile of LRG halos at r>25 kpc whereas baryons account for more than 50% of the mass at smaller radii. We calculate the total dark-to-baryonic mass ratio and show that it is consistent with measurements from weak lensing for environments dominated by massive early type galaxies. Finally, we divide the satellite galaxies in our sample into three luminosity bins and show that the satellite light profiles of all brightness levels are consistent with each other outside of roughly 25 kpc. At smaller radii we find evidence for a mild mass segregation with an increasing fraction of bright satellites close to the central LRG.

preprint2012arXiv

Revealing velocity dispersion as the best indicator of a galaxy&#39;s color, compared to stellar mass, surface mass density or morphology

Using data of nearby galaxies from the Sloan Digital Sky Survey we investigate whether stellar mass, central velocity dispersion, surface mass density, or the Sersic n parameter is best correlated with a galaxy&#39;s rest-frame color. Specifically, we determine how the mean color of galaxies varies with one parameter when another is fixed. When the stellar mass is fixed we see that strong trends remain with all other parameters, whereas residual trends are weaker when surface mass density, n, or velocity dispersion are fixed. Overall velocity dispersion is the best indicator of a galaxy&#39;s typical color, showing the largest residual color dependence when any of the other three parameters are fixed, and stellar mass is the poorest. Other studies have indicated that both the halo and black hole properties are better correlated with velocity dispersion than with stellar mass, surface mass density or Sersic n. Therefore, our results are consistent with a picture where a galaxy&#39;s star formation history and present star formation rate are determined to some significant degree by the current properties and assembly history of its dark matter halo and/or the feedback from its central super massive black hole.

preprint2012arXiv

Sizes and surface brightness profiles of quiescent galaxies at z ~ 2

We use deep Hubble Space Telescope Wide Field Camera 3 near-infrared imaging obtained of the GOODS-South field as part of the CANDELS survey to investigate a stellar mass-limited sample of quiescent galaxies at 1.5 < z < 2.5. We measure surface brightness profiles for these galaxies using a method that properly measures low surface brightness flux at large radii. We find that quiescent galaxies at z ~ 2 very closely follow Sersic profiles, with n_{median} = 3.7, and have no excess flux at large radii. Their effective radii are a factor ~ 4 smaller than those of low-redshift quiescent galaxies of similar mass. However, there is significant spread in sizes (sigma_{log r_e} = 0.24), with the largest z ~ 2 galaxies lying close to the z = 0 mass-size relation. We compare the stellar mass surface density profiles with those of massive elliptical galaxies in the Virgo cluster and confirm that most of the mass-growth which occurs between z ~ 2 and z = 0 must be due to accretion of material onto the outer regions of the galaxies. Additionally, we investigate the evolution in the size distribution of massive quiescent galaxies. We find that the minimum size growth required for z ~ 2 quiescent galaxies to fall within the z = 0 size distribution is a factor ~ 2 smaller than the total median size growth between z ~ 2 and z = 0.

preprint2012arXiv

Spatially resolved Halpha maps and sizes of 57 strongly star-forming galaxies at z~1 from 3D-HST: evidence for rapid inside-out assembly of disk galaxies

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

preprint2012arXiv

The Star Formation Rate Function for Redshift z~4-7 Galaxies: Evidence for a Uniform Build-Up of Star-Forming Galaxies During the First 3 Gyr of Cosmic Time

We combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-established z~4-7 UV LFs, UV-continuum slope trends with redshift and luminosity, and IRX-beta relations. They are well-described by Schechter relations. We extend the comparison baseline for SFR functions to z~2 by considering recent determinations of the Hα and mid-IR luminosity functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, Δα\sim+0.13, and show no clear evolution with cosmic time (z~0-7). In addition, we find that the characteristic value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth as a function of redshift from ~5 M_sun/yr at z~8, 650 Myr after the Big Bang, to ~100 M_sun/yr at z~2, ~2.5 Gyr later. Recent results at z~10, close to the onset of galaxy formation, are consistent with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the redshift range z~2-8, providing validation for our dust corrections. These results provide strong evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time.

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 &#34;normal&#34; 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 &#34;normal&#34; 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.

preprint2012arXiv

The stellar mass structure of massive galaxies from z=0 to z=2.5; surface density profiles and half-mass radii

We present stellar mass surface density profiles of a mass-selected sample of 177 galaxies at 0.5 < z < 2.5, obtained using very deep HST optical and near-infrared data over the GOODS-South field, including recent CANDELS data. Accurate stellar mass surface density profiles have been measured for the first time for a complete sample of high-redshift galaxies more massive than 10^10.7 M_sun. The key advantage of this study compared to previous work is that the surface brightness profiles are deconvolved for PSF smoothing, allowing accurate measurements of the structure of the galaxies. The surface brightness profiles account for contributions from complex galaxy structures such as rings and faint outer disks. Mass profiles are derived using radial rest-frame u-g color profiles and a well-established empirical relation between these colors and the stellar mass-to-light ratio. We derive stellar half-mass radii from the mass profiles, and find that these are on average ~25% smaller than rest-frame g band half-light radii. This average size difference of 25% is the same at all redshifts, and does not correlate with stellar mass, specific star formation rate, effective surface density, Sersic index, or galaxy size. Although on average the difference between half-mass size and half-light size is modest, for approximately 10% of massive galaxies this difference is more than a factor two. These extreme galaxies are mostly extended, disk-like systems with large central bulges. These results are robust, but could be impacted if the central dust extinction becomes high. ALMA observations can be used to explore this possibility. These results provide added support for galaxy growth scenarios wherein massive galaxies at these epochs grow by accretion onto their outer regions.

preprint2012arXiv

Which galaxy property is the best indicator of its host dark matter halo properties?

In this work we investigate the link between galaxy velocity dispersion, mass and other properties (color, morphology) with the properties of dark matter halos by comparing the clustering of galaxies at both fixed mass and velocity dispersion. We use the Sloan Digital Sky Survey to define a volume limited sample of massive galaxies complete in both stellar mass (>6e10 Msun) and velocity dispersion (>75 km/s). Using this sample we show that at fixed velocity dispersion there is no dependence of the clustering amplitude on stellar or dynamical mass. Conversely when stellar or dynamical mass are fixed there is a clear dependence of the clustering amplitude on velocity dispersion with higher dispersion galaxies showing a higher clustering amplitude. We also show that whilst when stellar or dynamical mass are fixed there remains a dependence of clustering amplitude on morphology, there is no such dependency when dispersion is fixed. However, we do see a dependence of the clustering amplitude on color when both mass and dispersion are fixed. Despite this, even when we restrict our samples to only elliptical or red galaxies the relationship between dispersion and clustering amplitude at fixed mass remains. It seems likely that the residual correlation with color is driven by satellite galaxies in massive halos being redder at fixed dispersion. The lack of a similar morphology dependence implies that the mechanism turning satellites red is not changing their morphology. Our central result is that velocity dispersion is more closely related to the clustering amplitude of galaxies than either stellar or dynamical mass. This implies that velocity dispersion is more tightly correlated with the halo properties that determine clustering, either halo mass or age, and supports the notion that the star formation history of a galaxy is more closely related to its halo properties than its overall mass.

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 (&#34;post-starburst&#34;) 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

Morphological evolution of galaxies from ultradeep HST WFC3 imaging: the Hubble sequence at z~2

We use ultradeep HST WFC3/IR imaging of the HUDF to investigate the rest-frame optical morphologies of a mass-selected sample of galaxies at z~2. We find a large variety of galaxy morphologies, ranging from large, blue, disk-like galaxies to compact, red, early-type galaxies. We derive rest-frame u-g color profiles for these galaxies and show that most z~2 galaxies in our sample have negative color gradients such that their cores are red. Although these color gradients may partly be caused by radial variations in dust content, they point to the existence of older stellar populations in the centers of z~2 galaxies. This result is consistent with an &#34;inside-out&#34; scenario of galaxy growth. We find that the median color gradient is fairly constant with redshift: [Delta(u-g_{rest})/Delta(log r)]_{median} = -0.46, -0.44 and -0.49 for z~2, z~1 and z=0, respectively. Using structural parameters derived from surface brightness profiles we confirm that at z~2 galaxy morphology correlates well with specific star formation rate. At the same mass, star forming galaxies have larger effective radii, bluer rest-frame u-g colors and lower Sersic indices than quiescent galaxies. These correlations are very similar to those at lower redshift, suggesting that the relations that give rise to the Hubble sequence at z=0 are already in place for massive galaxies at this early epoch.

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

Confirmation of the compactness of a z=1.91 quiescent galaxy with Hubble Space Telescope&#39;s Wide Field Camera 3

We present very deep Wide Field Camera 3 (WFC3) photometry of a massive, compact galaxy located in the Hubble Ultra Deep Field. This quiescent galaxy has a spectroscopic redshift z=1.91 and has been identified as an extremely compact galaxy by Daddi et al. 2005. We use new H-F160W imaging data obtained with Hubble Space Telescope/WFC3 to measure the deconvolved surface brightness profile to H = 28 mag arcsec**-2. We find that the surface brightness profile is well approximated by an n=3.7 Sersic profile. Our deconvolved profile is constructed by a new technique which corrects the best-fit Sersic profile with the residual of the fit to the observed image. This allows for galaxy profiles which deviate from a Sersic profile. We determine the effective radius of this galaxy: r_e=0.42 +- 0.14 kpc in the observed H-F160W-band. We show that this result is robust to deviations from the Sersic model used in the fit. We test the sensitivity of our analysis to faint &#34;wings&#34; in the profile using simulated galaxy images consisting of a bright compact component and a faint extended component. We find that due to the combination of the WFC3 imaging depth and our method&#39;s sensitivity to extended faint emission we can accurately trace the intrinsic surface brightness profile, and that we can therefore confidently rule out the existence of a faint extended envelope around the observed galaxy down to our surface brightness limit. These results confirm that the galaxy lies a factor of 10 off from the local 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&#39;. 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

Hubble Space Telescope WFC3 Grism Spectroscopy and Imaging of a Growing Compact Galaxy at z=1.9

We present HST/WFC3 grism spectroscopy of the brightest galaxy at z>1.5 in the GOODS-South WFC3 Early Release Science grism pointing, covering the wavelength range 0.9-1.7 micron. The spectrum is of remarkable quality and shows the redshifted Balmer lines Hbeta, Hgamma, and Hdelta in absorption at z=1.902, correcting previous erroneous redshift measurements from the rest-frame UV. The average rest-frame equivalent width of the Balmer lines is 8+-1 Angstrom, which can be produced by a post-starburst stellar population with a luminosity-weighted age of ~0.5 Gyr. The M/L ratio inferred from the spectrum implies a stellar mass of ~4x10^11 Msun. We determine the morphology of the galaxy from a deep WFC3 F160W image. Similar to other massive galaxies at z~2 the galaxy is compact, with an effective radius of 2.1+-0.3 kpc. Although most of the light is in a compact core, the galaxy has two red, smooth spiral arms that appear to be tidally-induced. The spatially-resolved spectroscopy demonstrates that the center of the galaxy is quiescent and the surrounding disk is forming stars, as it shows Hbeta in emission. The galaxy is interacting with a companion at a projected distance of 18 kpc, which also shows prominent tidal features. The companion has a slightly redder spectrum than the primary galaxy but is a factor of ~10 fainter and may have a lower metallicity. It is tempting to interpret these observations as &#34;smoking gun&#34; evidence for the growth of compact, quiescent high redshift galaxies through minor mergers, which has been proposed by several recent observational and theoretical studies. Interestingly both objects host luminous AGNs, as indicated by their X-ray luminosities, which implies that these mergers can be accompanied by significant black hole growth. This study illustrates the power of moderate dispersion, low background near-IR spectroscopy at HST resolution, which is now available with the WFC3 grism.

preprint2010arXiv

On Sizes, Kinematics, M/L Gradients, and Light Profiles of Massive Compact Galaxies at z~2

We present a detailed analysis of the structure and resolved stellar populations of simulated merger remnants, and compare them to observations of compact quiescent galaxies at z ~ 2. We find that major merging is a viable mechanism to produce systems of ~ 10^11 Msun and ~ 1 kpc size, provided the gas fraction at the time of final coalescence is high (~ 40%), and provided that the progenitors are compact star-forming galaxies, as expected at high redshift. Their integrated spectral energy distributions and velocity dispersions are in good agreement with the observations, and their position in the (v_{maj}/sigma, ellipticity) diagram traces the upper envelope of the distribution of lower redshift early-type galaxies. The simulated merger remnants show time- and sightline-dependent M/L ratio gradients that result from a superposition of radially dependent stellar age, stellar metallicity, and extinction. The median ratio of effective radius in rest-frame V-band light to that in mass surface density is ~ 2 during the quiescent remnant phase. This is typically expressed by a negative color gradient (i.e., red core), which we expect to correlate with the integrated color of the system. Finally, the simulations differ from the observations in their surface brightness profile shape. The simulated remnants are typically best fit by high (n >> 4) Sersic indices, whereas observed quiescent galaxies at z ~ 2 tend to be less cuspy (median n ~ 2.3). Limiting early star formation in the progenitors may be required to prevent the simulated merger remnants from having extended wings.

preprint2010arXiv

On the Masses of Galaxies in the Local Universe

We compare estimates of stellar mass, Mstar, and dynamical mass,Mdyn,for a sample of galaxies from the Sloan Digital Sky Survey (SDSS). We assume dynamical homology (i.e., Mdyn = dispersion**2 * Reff, and we find a tight but strongly non-linear relation: the best fit relation is Mstar = Mdyn**0.73, with an observed scatter of 0.15 dex. We also find that, at fixed Mstar, the ratio Mstar/Mdyn depends strongly on galaxy structure, as parameterized by Sersic index, n. The size of the differential effect is on the order of 0.6 dex across 2 < n < 10. The apparent n-dependence of Mstar/Mdyn is similar to expectations from simple models, indicating that assuming homology gives the wrong dynamical mass. We have also derived dynamical mass estimates that explicitly account for differences in galaxies&#39; profiles. Using this `structure-corrected&#39; dynamical mass estimator, M(dyn,n), the best fit relation is Mstar = M(dyn,n)**(0.92 +- 0.08) with an observed scatter of 0.13 dex. While the data are thus consistent with a linear relation, they do prefer a slightly shallower slope. Further, we see only a small residual trend in Mstar/M(dyn,n) with n. We find no statistically significant systematic trends in Mstar/M(dyn,n) as a function of observed quantities (e.g, apparent magnitude, redshift), or as a function of tracers of stellar populations. The net differential bias in Mstar/M(dyn,n) across a wide range of stellar populations and star formation activities is <= 0.12 dex. The very good agreement between stellar mass and structure-corrected dynamical mass strongly suggests that: 1.) galaxy non-homology has a major impact on dynamical mass estimates, and 2. there are not strong systematic biases in the stellar mass-to-light ratios derived from broadband optical SEDs. Further, these results suggest that that the central dark-to-luminous mass ratio has a relatively weak mass dependence.

preprint2010arXiv

Searching for a Gravitational Heating Signature in Nearby Luminous Ellipticals

We present a new deep optical study of a luminosity limited sample of nearby elliptical galaxies, attempting to observe the effects of gravitational interactions on the ISM of these objects. This study is motivated by recent observations of M86, a nearby elliptical galaxy that shows possible evidence for gas heating through a recent gravitational interaction. The complete sample includes luminous ellipticals in clusters, groups and the field. For each of the galaxies we objectively derive a tidal parameter which measures the deviation of the stellar body from a smooth, relaxed model and find that 73% of them show tidal disturbance signatures in their stellar bodies. This is the first time that such an analysis is done on a statistically complete sample and it confirms that elliptical galaxies continue to grow and evolve through gravitational interactions even in the local Universe. Our study of ellipticals in a wide range of interaction stages, along with available ISM data will attempt to shed light on this possibly alternative mechanism for maintaining the observed ISM temperatures of elliptical galaxies.

preprint2010arXiv

The Hubble Sequence beyond z=2 for Massive Galaxies: Contrasting Large Star-Forming and Compact Quiescent Galaxies

We present Hubble Space Telescope NIC2 morphologies of a spectroscopic sample of massive galaxies at z~2.3, by extending our sample of 9 compact quiescent galaxies (r_e~0.9 kpc) with 10 massive emission-line galaxies. The emission-line galaxies are classified by the nature of their ionized emission; there are six star-forming galaxies and four galaxies hosting an active galactic nucleus (AGN). The star-forming galaxies are the largest among the emission-line galaxies, with a median size of re = 2.8 kpc. The three galaxies with the highest star formation rates (> 100M \odot/yr) have irregular and clumpy morphologies. The AGN host galaxies are more similar to the compact quiescent galaxies in terms of their structures (re~1.1 kpc) and spectral energy distributions. The total sample clearly separates into two classes in a color-mass diagram: the large star-forming galaxies that form the blue cloud, and the compact quiescent galaxies on the red sequence. However, it is unclear how or even if the two classes are evolutionary related. Three out of six massive star-forming galaxies have dense cores and thus may passively evolve into compact galaxies due to fading of outer star-forming regions. For these galaxies a reverse scenario, in which compact galaxies grow inside-out by star formation is also plausible. We do caution though that the sample is small. Nonetheless, it is evident that a Hubble sequence of massive galaxies with strongly correlated galaxy properties is already in place at z > 2.

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&#39;s optical, with median r=26.1, and rest-frame U-V=1.6. About 60% of the sample have optical colors satisfying either the U- or the B-dropout color criteria, although ~50% of these galaxies have r>25.5. About 30% of the sample has SFRs from SED modeling consistent with zero. However, >80% of the sample is detected at 24 micron, with total infrared luminosities in the range (0.5-4.0)x10^13 Lsun. This implies the presence of either dust-enshrouded starburst activity (with SFRs of 600-4300 Msun/yr) and/or highly-obscured active galactic nuclei (AGN). The contribution of galaxies with M_star>2.5x10^11 Msun to the total stellar mass budget at z=3.5 is ~8%. We find an evolution by a factor of 2-7 and 3-22 from z~5 and z~6, respectively, to z=3.5. The previously found disagreement at the high-mass end between observed and model-predicted SMFs is now significant at the 3sigma level. However, systematic uncertainties dominate the total error budget, with errors up to a factor of ~8 in the densities, bringing the observed SMF in marginal agreement with the predicted SMF. Additional systematic uncertainties on the high-mass end could be introduced by either 1) the intense star-formation and/or the very common AGN activities as inferred from the MIPS 24 micron detections, and/or 2) contamination by a significant population of massive, old, and dusty galaxies at z~2.6.

preprint2010arXiv

The Multiwavelength Survey by Yale-Chile (MUSYC): Deep Medium-Band optical imaging and high quality 32-band photometric redshifts in the ECDF-S

We present deep optical 18-medium-band photometry from the Subaru telescope over the ~30&#39; x 30&#39; Extended Chandra Deep Field-South (ECDF-S), as part of the Multiwavelength Survey by Yale-Chile (MUSYC). This field has a wealth of ground- and space-based ancillary data, and contains the GOODS-South field and the Hubble Ultra Deep Field. We combine the Subaru imaging with existing UBVRIzJHK and Spitzer IRAC images to create a uniform catalog. Detecting sources in the MUSYC BVR image we find ~40,000 galaxies with R_AB<25.3, the median 5 sigma limit of the 18 medium bands. Photometric redshifts are determined using the EAZY code and compared to ~2000 spectroscopic redshifts in this field. The medium band filters provide very accurate redshifts for the (bright) subset of galaxies with spectroscopic redshifts, particularly at 0.1 < z < 1.2 and at z > 3.5. For 0.1 < z < 1.2, we find a 1 sigma scatter in Δz/(1+z) of 0.007, similar to results obtained with a similar filter set in the COSMOS field. As a demonstration of the data quality, we show that the red sequence and blue cloud can be cleanly identified in rest-frame color-magnitude diagrams at 0.1 < z < 1.2. We find that ~20% of the red-sequence-galaxies show evidence of dust-emission at longer rest-frame wavelengths. The reduced images, photometric catalog, and photometric redshifts are provided through the public MUSYC website.

preprint2009arXiv

A high stellar velocity dispersion for a compact massive galaxy at z=2.2

Recent studies have found that the oldest and most luminous galaxies in the early Universe are surprisingly compact, having stellar masses similar to present-day elliptical galaxies but much smaller sizes. This finding has attracted considerable attention as it suggests that massive galaxies have grown by a factor of ~five in size over the past ten billion years. A key test of these results is a determination of the stellar kinematics of one of the compact galaxies: if the sizes of these objects are as extreme as has been claimed, their stars are expected to have much higher velocities than those in present-day galaxies of the same mass. Here we report a measurement of the stellar velocity dispersion of a massive compact galaxy at redshift z=2.186, corresponding to a look-back time of 10.7 billion years. The velocity dispersion is very high at 510 (+165, -95) km/s, consistent with the mass and compactness of the galaxy inferred from photometric data and indicating significant recent structural and dynamical evolution of massive galaxies. The uncertainty in the dispersion was determined from simulations which include the effects of noise and template mismatch. However, we caution that some subtle systematic effect may influence the analysis given the low signal-to-noise ratio of our spectrum.

preprint2009arXiv

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

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

preprint2009arXiv

Optical Spectroscopy of Distant Red Galaxies

We present optical spectroscopic follow-up of a sample of Distant Red Galaxies (DRGs) with K < 22.5 (Vega), selected by J-K > 2.3, in the Hubble Deep Field South, the MS 1054-03 field, and the Chandra Deep Field South. Spectroscopic redshifts were obtained for 15 DRGs. Only 2 out of 15 DRGs are located at z < 2, suggesting a high efficiency to select high-redshift sources. From other spectroscopic surveys in the CDFS targeting intermediate to high redshift populations selected with different criteria, we find spectroscopic redshifts for a further 30 DRGs. We use the sample of spectroscopically confirmed DRGs to establish the high quality (scatter in Δz/(1+z) of ~ 0.05) of their photometric redshifts in the considered deep fields, as derived with EAZY (Brammer et al. 2008). Combining the spectroscopic and photometric redshifts, we find that 74% of DRGs with K < 22.5 lie at z > 2. The combined spectroscopic and photometric sample is used to analyze the distinct intrinsic and observed properties of DRGs at z < 2 and z > 2. In our photometric sample to K < 22.5, low-redshift DRGs are brighter in K than high-redshift DRGs by 0.7 mag, and more extincted by 1.2 mag in Av. Our analysis shows that the DRG criterion selects galaxies with different properties at different redshifts. Such biases can be largely avoided by selecting galaxies based on their rest-frame properties, which requires very good multi-band photometry and high quality photometric redshifts.

preprint2009arXiv

Recovering Stellar Population Properties and Redshifts from Broad-Band Photometry of Simulated Galaxies: Lessons for SED Modeling

We present a detailed analysis of our ability to determine stellar masses, ages, reddening and extinction values, and star formation rates of high-redshift galaxies by modeling broad-band SEDs with stellar population synthesis. In order to do so, we computed synthetic optical-to-NIR SEDs for model galaxies taken from hydrodynamical merger simulations placed at redshifts 1.5 < z < 3. Viewed under different angles and during different evolutionary phases, the simulations represent a wide variety of galaxy types (disks, mergers, spheroids). We show that simulated galaxies span a wide range in SEDs and color, comparable to these of observed galaxies. In all star-forming phases, dust attenuation has a large effect on colors, SEDs, and fluxes. The broad-band SEDs were then fed to a standard SED modeling procedure and resulting stellar population parameters were compared to their true values. Disk galaxies generally show a decent median correspondence between the true and estimated mass and age, but suffer from large uncertainties. During the merger itself, we find larger offsets (e.g., log M_recovered - log M_true = -0.13^{+0.10}_{-0.14}). E(B-V) values are generally recovered well, but the estimated total visual absorption Av is consistently too low, increasingly so for larger optical depths. Since the largest optical depths occur during the phases of most intense star formation, it is for the highest SFRs that we find the largest underestimates. The masses, ages, E(B-V), Av, and SFR of merger remnants (spheroids) are very well reproduced. We discuss possible biases in SED modeling results caused by mismatch between the true and template star formation history, dust distribution, metallicity variations and AGN contribution.

preprint2009arXiv

Star formation in the CDFS: observations confront simulations

We investigate the star formation history of the universe using FIREWORKS, a multiwavelength survey of the CDFS. We study the evolution of the specific star formation rate (sSFR) with redshift in different mass bins from z = 0 to z ~ 3. We find that the sSFR increases with redshift for all masses. The logarithmic increase of the sSFR with redshift is nearly independent of mass, but this cannot yet be verified at the lowest-mass bins at z > 0.8, due to incompleteness. We convert the sSFRs to a dimensionless growth rate to facilitate a comparison with a semi-analytic galaxy formation model that was implemented on the Millennium Simulation. The model predicts that the growth rates and sSFRs increase similarly with redshift for all masses, consistent with the observations. However, we find that for all masses, the inferred observed growth rates increase more rapidly with redshift than the model predictions. We discuss several possible causes for this discrepancy, ranging from field-to-field variance, conversions to SFR, and shape of the IMF. We find that none of these can solve the discrepancy completely. We conclude that the models need to be adapted to produce the steep increase in growth rate between redshift z=0 and z=1.

preprint2009arXiv

The Evolution of the Stellar Mass Function of Galaxies from z=4.0 and the First Comprehensive Analysis of its Uncertainties: Evidence for Mass-dependent Evolution

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

preprint2009arXiv

The Frequency of Tidal Features Associated with Nearby Luminous Elliptical Galaxies from a Statistically Complete Sample

We present a deep broadband optical imaging study of a complete sample of luminous elliptical galaxies (M_B<-20) at distances 15 Mpc - 50 Mpc, selected from the Tully catalog of nearby galaxies. The images are flat to ~0.35% across the 20&#39; field and reach a V band depth of 27.7 mag arcsec^-2. We derive an objective tidal interaction parameter for all galaxies and find that 73% of them show tidal disturbance signatures in their stellar bodies. This is the first time that such an analysis is done on a statistically complete sample and it confirms that tidal features in ellipticals are common even in the local Universe. From the dynamical time of the sample galaxies at the innermost radius where tidal features are detected we estimate the mass assembly rate of nearby ellipticals to be dM/M 0.2 per Gyr with large uncertainty. We explore the relation between gravitational interaction signatures and the galaxy environment and find that galaxies in clusters are less disturbed than group and field galaxies. We also study how these interactions affect the broadband colors of ellipticals and find a moderate correlation, suggesting that the mergers are not accompanied by significant star-formation. Lastly, we find no correlation between AGN activity, as measured by 6cm radio emission, and large scale tidal distortions. This implies that gravitational interactions are not the only, and perhaps not the most important, trigger of nuclear activity. In summary, we find that elliptical galaxies in groups and low density environments continue to grow at the present day through mostly &#34;dry&#34; mergers involving little star formation.

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 &#34;spectra&#34; 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 &#34;Y&#34; dwarfs could be identified using only two of the filters.

preprint2009arXiv

The Relation Between Compact, Quiescent High Redshift Galaxies and Massive Nearby Elliptical Galaxies: Evidence for Hierarchical, Inside-Out Growth

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

preprint2006arXiv

Comparing Dynamical and Photometric Mass Estimates of Low- and High-Redshift Galaxies: Random and Systematic Uncertainties

We determine the importance of redshift-dependent systematic effects in the determination of stellar masses from broad band spectral energy distributions (SEDs), using high quality kinematic and photometric data of early-type galaxies at z~1 and z~0. We find that photometric masses of z~1 galaxies can be systematically different, by up to a factor of 2, from photometric masses of z~0 galaxies with the same dynamical mass. The magnitude of this bias depends on the choice of stellar population synthesis model and the rest-frame wavelength range used in the fits. The best result, i.e., without significant bias, is obtained when rest-frame optical SEDs are fitted with models from Bruzual&Charlot (2003). When the SEDs are extended to the rest-frame near-IR, a bias is introduced: photometric masses of the z~1 galaxies increase by a factor of 2 relative to the photometric masses of the z~0 galaxies. When we use the Maraston (2005) models, the photometric masses of the z~1 galaxies are low relative to the photometric masses of the z~0 galaxies by a factor of ~1.8. This offset occurs both for fits based on rest-frame optical SEDs, and fits based on rest-frame optical+near-IR SEDs. The results indicate that model uncertainties produce uncertainties as high as a factor of 2.5 in mass estimates from rest-frame near-IR photometry, independent of uncertainties due to unknown star formation histories.

preprint2002arXiv

The Fundamental Plane at z=1.27: First Calibration of the Mass Scale of Red Galaxies at Redshifts z>1

We present results on the Fundamental Plane (FP) of early-type galaxies in the cluster RDCS J0848+4453 at z=1.27. Internal velocity dispersions of three K-selected early-type galaxies are determined from deep Keck spectra. Structural parameters are determined from HST NICMOS images. The galaxies show substantial offsets from the FP of the nearby Coma cluster, as expected from passive evolution of their stellar populations. The offsets from the FP can be expressed as offsets in M/L ratio. The M/L ratios of the two most massive galaxies are consistent with an extrapolation of results obtained at z=0.02-0.83. The evolution of early-type galaxies with masses >10^11 M_sun is well described by ln M/L(B) = (-1.06 +- 0.09) z, corresponding to passive evolution of -1.50 +- 0.13 mag at z=1.3. Ignoring selection effects, the best fitting stellar formation redshift is z*=2.6, corresponding to a luminosity weighted age at the epoch of observation of ~2 Gyr. The M/L ratios of these two galaxies are also in excellent agreement with predictions from models that include progenitor bias. The third galaxy is a factor ~10 less massive than the other two, shows strong Balmer absorption lines in its spectrum, and is offset from the Coma Fundamental Plane by 2.9 mag in rest-frame B. Despite their large range in M/L ratios, all three galaxies fall in the ``Extremely Red Object&#39;&#39; (ERO) class with I-H>3 and R-K>5, and our results show that it is hazardous to use simple models for converting luminosity to mass for these objects. Measurements of M/L ratios at high redshift can be considered first steps to empirically disentangle luminosity and mass evolution at the high mass end of the galaxy population, lifting an important degeneracy in the interpretation of evolution of the luminosity function. [SHORTENED]

preprint2001arXiv

Cosmic-Ray Rejection by Laplacian Edge Detection

Conventional algorithms for rejecting cosmic-rays in single CCD exposures rely on the contrast between cosmic-rays and their surroundings, and may produce erroneous results if the Point Spread Function (PSF) is smaller than the largest cosmic-rays. This paper describes a robust algorithm for cosmic-ray rejection, based on a variation of Laplacian edge detection. The algorithm identifies cosmic-rays of arbitrary shapes and sizes by the sharpness of their edges, and reliably discriminates between poorly sampled point sources and cosmic-rays. Examples of its performance are given for spectroscopic and imaging data, including HST WFPC2 images.