Source author record

Mariska Kriek

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

87works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

87 published item(s)

preprint2022arXiv

CO Emission, Molecular Gas, and Metallicity in Main-Sequence Star-Forming Galaxies at $z\sim2.3$

We present observations of CO(3-2) in 13 main-sequence $z=2.0-2.5$ star-forming galaxies at $\log(M_*/M_{\odot})=10.2-10.6$ that span a wide range in metallicity (O/H) based on rest-optical spectroscopy. We find that CO(3-2)/SFR decreases with decreasing metallicity, implying that the CO luminosity per unit gas mass is lower in low-metallicity galaxies at $z\sim2$. We constrain the CO-to-H$_2$ conversion factor ($α_{\text{CO}}$) and find that $α_{\text{CO}}$ inversely correlates with metallicity at $z\sim2$. We derive molecular gas masses ($M_{\text{mol}}$) and characterize the relations among $M_*$, SFR, $M_{\text{mol}}$, and metallicity. At $z\sim2$, $M_{\text{mol}}$ increases and molecular gas fraction ($M_{\text{mol}}$/$M_*$) decrease with increasing $M_*$, with a significant secondary dependence on SFR. Galaxies at $z\sim2$ lie on a near-linear molecular KS law that is well-described by a constant depletion time of 700 Myr. We find that the scatter about the mean SFR-$M_*$, O/H-$M_*$, and $M_{\text{mol}}$-$M_*$ relations is correlated such that, at fixed $M_*$, $z\sim2$ galaxies with larger $M_{\text{mol}}$ have higher SFR and lower O/H. We thus confirm the existence of a fundamental metallicity relation at $z\sim2$ where O/H is inversely correlated with both SFR and $M_{\text{mol}}$ at fixed $M_*$. These results suggest that the scatter of the $z\sim2$ star-forming main sequence, mass-metallicity relation, and $M_{\text{mol}}$-$M_*$ relation are primarily driven by stochastic variations in gas inflow rates. We place constraints on the mass loading of galactic outflows and perform a metal budget analysis, finding that massive $z\sim2$ star-forming galaxies retain only 30% of metals produced, implying that a large mass of metals resides in the circumgalactic medium.

preprint2022arXiv

IR SED and Dust Masses of Sub-solar Metallicity Galaxies at z~2.3

We present results from ALMA 1.2mm continuum observations of a sample of 27 star-forming galaxies at z=2.1-2.5 from the MOSFIRE Deep Evolution Field (MOSDEF) survey. These galaxies have gas-phase metallicity and star-formation rate measurements from Hb, [OIII], Ha, and [NII]. Using stacks of Spitzer, Herschel, and ALMA photometry (rest-frame ~ 8-400$μ$m), we examine the IR SED of high-redshift subsolar metallicity (~0.5 $Z_{\odot}$) LIRGs. We find that the data agree well with an average SED template of higher luminosity local low-metallicity dwarf galaxies (reduced $χ^2$ of 1.8). When compared with the commonly used templates for solar-metallicity local galaxies or high-redshift LIRGs and ULIRGs, even in the most favorable case (with reduced $χ^2$ of 2.8), the templates are rejected at >98% confidence level. The broader and hotter IR SED of both the local dwarfs and high-redshift subsolar metallicity galaxies may result from different grain properties, a clumpy dust geometry, or a harder/more intense ionizing radiation field that heats the dust to higher temperatures. The obscured SFR indicated by the FIR emission of the subsolar metallicity galaxies is only ~ 60% of the total SFR, which is considerably lower than that of the local LIRGs with ~ 96-97% obscured fractions. Due to the evolving IR SED shape, the local LIRG templates fit to mid-IR data can overestimate the Rayleigh-Jeans tail measurements at z~2 by a factor of 2-20, and these templates underestimate IR luminosities if fit to the observed ALMA fluxes by >0.4dex. At a given stellar mass or metallicity, dust masses at z~2.3 are an order of magnitude higher than those at z~0. Given the predicted molecular gas mass fractions, the observed z~2.3 dust-to-stellar mass ratios suggest lower dust-to-molecular gas masses than in local galaxies at the same metallicity.

preprint2022arXiv

Quenching and the UVJ diagram in the SIMBA cosmological simulation

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

preprint2022arXiv

Reconciling the Results of the z~2 MOSDEF and KBSS-MOSFIRE Surveys

The combination of the MOSDEF and KBSS-MOSFIRE surveys represents the largest joint investment of Keck/MOSFIRE time to date, with ~3000 galaxies at 1.4<=z<=3.8, roughly half of which are at z~2. MOSDEF is photometric- and spectroscopic-redshift selected with a rest-optical magnitude limit, while KBSS-MOSFIRE is primarily selected based on rest-UV colors and a rest-UV magnitude limit. Analyzing both surveys in a uniform manner with consistent spectral-energy-distribution (SED) models, we find that the MOSDEF z~2 targeted sample has a higher median M_* and redder rest U-V color than the KBSS-MOSFIRE z~2 targeted sample, and a smaller median SED-based SFR and sSFR (SFR(SED) and sSFR(SED)). Specifically, MOSDEF targeted a larger population of red galaxies with U-V and V-J >=1.25, while KBSS-MOSFIRE contains more young galaxies with intense star formation. Despite these differences in the z~2 targeted samples, the subsets of the surveys with multiple emission lines detected and analyzed in previously published work are much more similar. All median host-galaxy properties with the exception of stellar population age -- i.e., M_*, SFR(SED), sSFR(SED), A_V, and UVJ colors -- agree within the uncertainties. Additionally, when uniform emission-line fitting and stellar Balmer absorption correction techniques are applied, there is no significant offset between the two samples in the [OIII]$λ$5008/H$β$ vs. [NII]$λ$6585/H$α$ diagnostic diagram, in contrast to previously-reported discrepancies. We can now combine the MOSDEF and KBSS-MOSFIRE surveys to form the largest z~2 sample with moderate-resolution rest-optical spectra and construct the fundamental scaling relations of star-forming galaxies during this important epoch.

preprint2022arXiv

Recovering the star formation histories of recently-quenched galaxies: the impact of model and prior choices

Accurate models of the star formation histories (SFHs) of recently-quenched galaxies can provide constraints on when and how galaxies shut down their star formation. The recent development of "non-parametric" SFH models promises the flexibility required to make these measurements. However, model and prior choices significantly affect derived SFHs, particularly for post-starburst galaxies (PSBs) which have sharp changes in their recent SFH. In this paper, we create mock PSBs, then use the Prospector SED fitting software to test how well four different SFH models recover key properties. We find that a two-component parametric model performs well for our simple mock galaxies, but is sensitive to model mismatches. The fixed- and flexible-bin non-parametric models included in Prospector are able to rapidly quench a major burst of star formation, but systematically underestimate the post-burst age by up to 200 Myr. We develop a custom SFH model that allows for additional flexibility in the recent SFH. Our flexible non-parametric model is able to constrain post-burst ages with no significant offset and just ~90 Myr of scatter. Our results suggest that while standard non-parametric models are able to recover first-order quantities of the SFH (mass, SFR, average age), accurately recovering higher-order quantities (burst fraction, quenching time) requires careful consideration of model flexibility. These mock recovery tests are a critical part of future SFH studies. Finally, we show that our new, public SFH model is able to accurately recover the properties of mock star-forming and quiescent galaxies and is suitable for broader use in the SED fitting community.

preprint2022arXiv

Star Formation Suppresion by Tidal Removal of Cold Molecular Gas from an Intermediate-Redshift Massive Post-Starburst Galaxy

Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all galaxies. We present the discovery of vast molecular tidal features extending up to 64 kpc outside of a massive z=0.646 post-starburst galaxy that recently concluded its primary star-forming episode. The tidal tails contain (1.2 +/- 0.1)x10^10 Msun of molecular gas, 47 +/- 5 % of the total cold gas reservoir of the system. Both the scale and magnitude of the molecular tidal features are unprecedented compared to all known nearby or high-redshift merging systems. We infer that the cold gas was stripped from the host galaxies during the merger, which is most likely responsible for triggering the initial burst phase and the subsequent suppression of star formation. While only a single example, this result shows that galaxy mergers can regulate the cold gas contents in distant galaxies by directly removing a large fraction of the molecular gas fuel, and plausibly suppress star formation directly, a qualitatively different physical mechanism than feedback-driven outflows.

preprint2022arXiv

The Compact Structures of Massive $z\sim0.7$ Post-Starburst Galaxies in the SQuIGG$\vec{L}$E Sample

We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift.

preprint2022arXiv

The MOSDEF Survey: Implications of the Lack of Evolution in the Dust Attenuation-Mass Relation to z~2

We investigate the relationship between dust attenuation and stellar mass ($M_*$) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey at $z\sim2.3$ and the Sloan Digital Sky Survey (SDSS) at $z\sim0$, augmenting the latter optical dataset with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H$α$ and H$β$ emission lines, and photometric measurements of the rest-UV stellar continuum. The H$α$/H$β$ ratio is used to determine the magnitude of attenuation at the wavelength of H$α$, $A_{{\rm H}α}$. Rest-UV colors and spectral-energy-distribution fitting are used to estimate $A_{1600}$, the magnitude of attenuation at a rest wavelength of 1600Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and $M_*$ over the redshift range $z\sim0$ to $z\sim2.3$. Folding in the latest estimates of the evolution of $M_{\rm dust}$, $({M_{\rm dust}}/{M_{\rm gas}})$, and gas surface density at fixed $M_*$, we find that the expected $M_{\rm dust}$ and dust mass surface density are both significantly higher at $z\sim2.3$ than at $z\sim0$. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation vs. $M_*$, it is essential to determine the relationship between metallicity and $({M_{\rm dust}}/{M_{\rm gas}})$, the dust mass absorption coefficient, and dust geometry, and the evolution of these relations and quantities from $z\sim0$ to $z\sim2.3$.

preprint2021arXiv

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

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

preprint2021arXiv

Now you see it, now you don't: Star formation truncation precedes the loss of molecular gas by ~100 Myr in massive post-starburst galaxies at z~0.6

We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of the SQuIGGLE program. Early results showed that two post-starburst galaxies host large H$_2$ reservoirs despite their low inferred star formation rates. Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with $M_{H_2}\gtrsim10^9 M_{\odot}$. Given their high stellar masses, this mass limit corresponds to an average gas fraction of $\langle f_{H_2} \equiv M_{H_2}/M_{\star} \rangle \sim7\%$, or ${\sim}14\%$ using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with the $D_n4000$ spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched $\lesssim 150$ Myr prior to observation host detectable CO(2-1) emission, while older post-starburst galaxies are undetected. The large $\mathrm{H_2}$ reservoirs and low star formation rates in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100-200 Myrs, the SQuIGGLE galaxies require the additional and efficient heating or removal of cold gas to bring their low star formation rates in line with standard $\mathrm{H_2}$ scaling relations.

preprint2021arXiv

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process.

preprint2021arXiv

The MOSDEF Survey: A Comprehensive Analysis of the Rest-optical Emission-line Properties of $z\sim 2.3$ Star-forming Galaxies

We analyze the rest-optical emission-line spectra of $z\sim2.3$ star-forming galaxies in the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. In investigating the origin of the well-known offset between the sequences of high-redshift and local galaxies in the [O III]5008/H$β$ vs. [N II]6585/H$α$ ("[N II] BPT") diagram, we define two populations of $z\sim2.3$ MOSDEF galaxies. These include the "high" population that is offset towards higher [O III]5008/H$β$ and/or [N II]6585/H$α$ with respect to the local SDSS sequence and the "low" population that overlaps the SDSS sequence. These two groups are also segregated within the [O III]5008/H$β$ vs. [S II]6718,6733/H$α$ and the [O III]4960,5008/[O II]3727,3730 (O$_{32}$) vs. ([O III]4960,5008+[O II]3727,3730)/H$β$ (R$_{23}$) diagram, which suggests qualitatively that star-forming regions in the more offset galaxies are characterized by harder ionizing spectra at fixed nebular oxygen abundance. We also investigate many galaxy properties of the split sample and find that the "high" sample is on average smaller in size and less massive, but has higher specific star-formation rate and star-formation-rate surface density values and is slightly younger compared to the "low" population. From Cloudy+BPASS photoionization models, we estimate that the "high" population has a lower stellar metallicity (i.e., harder ionizing spectrum) but slightly higher nebular metallicity and higher ionization parameter compared to the "low" population. While the "high" population is more $α$-enhanced (i.e., higher $α$/Fe) than the "low" population, both samples are significantly more $α$-enhanced compared to local star-forming galaxies with similar rest-optical line ratios. These differences must be accounted for in all high-redshift star-forming galaxies -- not only those "offset" from local excitation sequences.

preprint2021arXiv

The MOSDEF Survey: Environmental dependence of the gas-phase metallicity of galaxies at $1.4 \leq z \leq 2.6$

Using the near-IR spectroscopy of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the role of local environment in the gas-phase metallicity of galaxies. The local environment measurements are derived from accurate and uniformly calculated photometric redshifts with well-calibrated probability distributions. Based on rest-frame optical emission lines, [NII]$\lambda6584$ and H$α$, we measure gas-phase oxygen abundance of 167 galaxies at $1.37\leq z\leq1.7$ and 303 galaxies at $2.09\leq z\leq2.61$, located in diverse environments. We find that at $z\sim1.5$, the average metallicity of galaxies in overdensities with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.8}M_\odot$ is higher relative to their field counterparts by $0.094\pm0.051$, $0.068\pm0.028$ and $0.052\pm0.043$ dex, respectively. However, this metallicity enhancement does not exist at higher redshift, $z\sim2.3$, where, compared to the field galaxies, we find $0.056\pm0.043$, $0.056\pm0.028$ and $0.096\pm 0.034$ dex lower metallicity for galaxies in overdense environments with $M_*\sim10^{9.8}M_\odot, 10^{10.2}M_\odot$ and $10^{10.7}M_\odot$, respectively. Our results suggest that, at $1.37\leq z\leq2.61$, the variation of mass-metallicity relation with local environment is small ($<0.1$dex), and reverses at $z\sim2$. Our results support the hypothesis that, at the early stages of cluster formation, owing to efficient gas cooling, galaxies residing in overdensities host a higher fraction of pristine gas with prominent primordial gas accretion, which lowers their gas-phase metallicity compared to their coeval field galaxies. However, as the Universe evolves to lower redshifts ($z\lesssim2$), shock-heated gas in overdensities cannot cool down efficiently, and galaxies become metal-rich rapidly due to the suppression of pristine gas inflow and re-accretion of metal-enriched outflows in overdensities.

preprint2020arXiv

Color gradients along the quiescent galaxy sequence: clues to quenching and structural growth

This Letter examines how the sizes, structures, and color gradients of galaxies change along the quiescent sequence. Our sample consists of ~400 quiescent galaxies at $1.0\le z\le2.5$ and $10.1 \le \log{M_*/M_\odot}\le11.6$ in three CANDELS fields. We exploit deep multi-band HST imaging to derive accurate mass profiles and color gradients, then use an empirical calibration from rest-frame UVJ colors to estimate galaxy ages. We find that -- contrary to previous results -- the youngest quiescent galaxies are not significantly smaller than older quiescent galaxies at fixed stellar mass. These `post-starburst' galaxies only appear smaller in half-light radii because they have systematically flatter color gradients. The strength of color gradients in quiescent galaxies is a clear function of age, with older galaxies exhibiting stronger negative color gradients (i.e., redder centers). Furthermore, we find that the central mass surface density $Σ_1$ is independent of age at fixed stellar mass, and only weakly depends on redshift. This finding implies that the central mass profiles of quiescent galaxies do not significantly change with age; however, we find that older quiescent galaxies have additional mass at large radii. Our results support the idea that building a massive core is a necessary requirement for quenching beyond $z=1$, and indicate that post-starburst galaxies are the result of a rapid quenching process that requires structural change. Furthermore, our observed color gradient and mass profile evolution supports a scenario where quiescent galaxies grow inside-out via minor mergers.

preprint2020arXiv

The MOSDEF Survey: Calibrating the relationship between H$α$ star-formation rate and radio continuum luminosity at $1.4 < z < 2.6$

The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$α$ SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 < z < 2.6$ with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis we compare the observed radio continuum luminosities with those predicted from the dust-corrected H$α$ SFR assuming a range of $z\sim0$ relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H$α$ SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to $z\sim 2$.

preprint2020arXiv

The MOSDEF Survey: Kinematic and Structural Evolution of Star-Forming Galaxies at $1.4\leq z\leq 3.8$

We present ionized gas kinematics for 681 galaxies at $z\sim 1.4-3.8$ from the MOSFIRE Deep Evolution Field survey, measured using models which account for random galaxy-slit misalignments together with structural parameters derived from CANDELS Hubble Space Telescope (HST) imaging. Kinematics and sizes are used to derive dynamical masses. Baryonic masses are estimated from stellar masses and inferred gas masses from dust-corrected star formation rates (SFRs) and the Kennicutt-Schmidt relation. We measure resolved rotation for 105 galaxies. For the remaining 576 galaxies we use models based on HST imaging structural parameters together with integrated velocity dispersions and baryonic masses to statistically constrain the median ratio of intrinsic ordered to disordered motion, $V/σ_{V,0}$. We find that $V/σ_{V,0}$ increases with increasing stellar mass and decreasing specific SFR (sSFR). These trends may reflect marginal disk stability, where systems with higher gas fractions have thicker disks. For galaxies with detected rotation we assess trends between their kinematics and mass, sSFR, and baryon surface density ($Σ_{\mathrm{bar},e}$). Intrinsic dispersion correlates most with $Σ_{\mathrm{bar},e}$ and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at $z\sim 1.4-3.8$ are baryon dominated within their effective radii ($R_E$), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within $R_E$, $f_{\mathrm{bar}}$, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, $f_{\mathrm{bar}}$ does not appear to vary with stellar mass but increases with decreasing $R_E$ and increasing $Σ_{\mathrm{bar},e}$. For galaxies at $z\geq2$, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

preprint2020arXiv

The MOSDEF Survey: the Variation of the Dust Attenuation Curve with Metallicity

We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$α$, H$β$, and [NII]$λ6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$α$-to-H$β$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3σ$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies.

preprint2020arXiv

The MOSDEF-LRIS Survey: The Interplay Between Massive Stars and Ionized Gas in High-Redshift Star-Forming Galaxies

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.

preprint2020arXiv

The Role of Active Galactic Nuclei in the Quenching of Massive Galaxies in the SQuiGGLE Survey

We study the incidence of nuclear activity in a large sample of massive post-starburst galaxies at z~0.7 selected from the Sloan Digital Sky Survey, and identify active galactic nuclei based on radio continuum and optical emission lines. Over our mass range of 10^10.6-10^11.5 Msun, the incidence of radio activity is weakly dependent on stellar mass and independent of stellar age, while radio luminosity depends strongly on stellar mass. Optical nuclear activity incidence depends most strongly on the Dn4000 line index, a proxy for stellar age, with an active fraction that is ~ten times higher in the youngest versus oldest post-starburst galaxies. Since a similar trend is seen between age and molecular gas fractions, we argue that, like in local galaxies, the age trend reflects a peak in available fueling rather than feedback from the central black hole on the surrounding galaxy.

preprint2019arXiv

Half-mass radii of quiescent and star-forming galaxies evolve slowly from 0 < z < 2.5: implications for galaxy assembly histories

We use high-resolution, multi-band imaging of ~16,500 galaxies in the CANDELS fields at 0 < z < 2.5 to study the evolution of color gradients and half-mass radii over cosmic time. We find that galaxy color gradients at fixed mass evolve rapidly between z~2.5 and z~1, but remain roughly constant below z~1. This result implies that the sizes of both star-forming and quiescent galaxies increase much more slowly than previous studies found using half-light radii. The half-mass radius evolution of quiescent galaxies is fully consistent with a model which uses observed minor merger rates to predict the increase in sizes due to the accretion of small galaxies. Progenitor bias may still contribute to the growth of quiescent galaxies, particularly if we assume a slower timescale for the minor merger growth model. The slower half-mass radius evolution of star-forming galaxies is in tension with cosmological simulations and semi-analytic galaxy models. Further detailed, consistent comparisons with simulations are required to place these results in context.

preprint2019arXiv

The MOSDEF survey: a census of AGN-driven ionized outflows at $z = 1.4-3.8$

Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] emission line spectra. We detect outflows in $17\%$ of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of $\sim 400-3500 ~\mathrm{km~s}^{-1}$ and spatial extents of $0.3-11.0$ kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically $0.1-1$ and increase with AGN luminosity, capable of exceeding unity at $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~erg~s}^{-1}$. In these more luminous sources the ionized outflow alone is likely sufficient to regulate star formation, and when combined with outflowing neutral and molecular gas may be able to quench star formation in their host galaxies.

preprint2019arXiv

The MOSDEF Survey: The Metallicity Dependence of X-ray Binary Populations at $z\sim2$

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.

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

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

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

preprint2016arXiv

Science Impacts of the SPHEREx All-Sky Optical to Near-Infrared Spectral Survey: Report of a Community Workshop Examining Extragalactic, Galactic, Stellar and Planetary Science

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $μ$m [with R$\sim$41.4] and 4.18 and 5.00 $μ$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

preprint2016arXiv

The MOSDEF Survey: Detection of [OIII]$λ$4363 and the direct-method oxygen abundance of a star-forming galaxy at z=3.08

We present measurements of the electron-temperature based oxygen abundance for a highly star-forming galaxy at z=3.08, COSMOS-1908. This is the highest redshift at which [OIII]$λ$4363 has been detected, and the first time that this line has been measured at z>2. We estimate an oxygen abundance of 12+log(O/H)$=8.00^{+0.13}_{-0.14}$. This galaxy is a low-mass ($10^{9.3}$ M$_{\odot}$), highly star-forming ($\sim50$ M$_{\odot}$ yr$^{-1}$) system that hosts a young stellar population ($\sim160$ Myr). We investigate the physical conditions of the ionized gas in COSMOS-1908 and find that this galaxy has a high ionization parameter, little nebular reddening ($E(B-V)_{\rm gas}<0.14$), and a high electron density ($n_e\sim500$ cm$^{-3}$). We compare the ratios of strong oxygen, neon, and hydrogen lines to the direct-method oxygen abundance for COSMOS-1908 and additional star-forming galaxies at z=0-1.8 with [OIII]$λ$4363 measurements, and show that galaxies at z$\sim$1-3 follow the same strong-line correlations as galaxies in the local universe. This agreement suggests that the relationship between ionization parameter and O/H is similar for z$\sim$0 and high-redshift galaxies. These results imply that metallicity calibrations based on lines of oxygen, neon, and hydrogen do not strongly evolve with redshift and can reliably estimate abundances out to z$\sim$3, paving the way for robust measurements of the evolution of the mass-metallicity relation to high redshift.

preprint2016arXiv

The MOSDEF Survey: Dynamical and Baryonic Masses and Kinematic Structures of Star-Forming Galaxies at $1.4 \leq z \leq 2.6$

We present H$α$ gas kinematics for 178 star-forming galaxies at z~2 from the MOSFIRE Deep Evolution Field survey. We have developed models to interpret the kinematic measurements from fixed-angle multi-object spectroscopy, using structural parameters derived from CANDELS HST/F160W imaging. For 35 galaxies we measure resolved rotation with a median $(V/σ_{V,0})_{R_E}=2.1$. We derive dynamical masses from the kinematics and sizes and compare them to baryonic masses, with gas masses estimated from dust-corrected H$α$ star formation rates (SFRs) and the Kennicutt-Schmidt relation. When assuming that galaxies with and without observed rotation have the same median $(V/σ_{V,0})_{R_E}$, we find good agreement between the dynamical and baryonic masses, with a scatter of $σ_{RMS}=0.34$dex and a median offset of $Δ\log_{10}M=0.04$dex. This comparison implies a low dark matter fraction (8% within an effective radius) for a Chabrier initial mass function (IMF), and disfavors a Salpeter IMF. Moreover, the requirement that $M_{dyn}/M_{baryon}$ should be independent of inclination yields a median value of $(V/σ_{V,0})_{R_E}=2.1$ for galaxies without observed rotation. If instead we treat the galaxies without detected rotation as early-type galaxies, the masses are also in reasonable agreement ($Δ\log_{10}M=-0.07$dex, $σ_{RMS}=0.37$dex). The inclusion of gas masses is critical in this comparison; if gas masses are excluded there is an increasing trend of $M_{dyn}/M_{*}$ with higher specific SFR (SSFR). Furthermore, we find indications that $V/σ$ decreases with increasing H$α$ SSFR for our full sample, which may reflect disk settling. We also study the Tully-Fisher relation and find that at fixed stellar mass $S_{0.5}=(0.5V_{2.2}^2+σ_{V,0}^2)^{1/2}$ was higher at earlier times. At fixed baryonic mass, we observe the opposite trend. [abridged]

preprint2016arXiv

The MOSDEF Survey: The strong agreement between H-alpha and UV-to-FIR star formation rates for z~2 star-forming galaxies

We present the first direct comparison between Balmer line and panchromatic SED-based SFRs for z~2 galaxies. For this comparison we used 17 star-forming galaxies selected from the MOSFIRE Deep Evolution Field (MOSDEF) survey, with $3σ$ detections for H$α$ and at least two IR bands (Spitzer/MIPS 24$μ$m and Herschel/PACS 100 and 160$μ$m, and in some cases Herschel/SPIRE 250, 350, and 500$μ$m). The galaxies have total IR (8-1000$μ$m) luminosities of $\sim10^{11.4}-10^{12.4}\,\textrm{L}_\odot$ and star-formation rates (SFRs) of $\sim30-250\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We fit the UV-to-far-IR SEDs with flexible stellar population synthesis (FSPS) models - which include both stellar and dust emission - and compare the inferred SFRs with the SFR(H$α$,H$β$) values corrected for dust attenuation using Balmer decrements. The two SFRs agree with a scatter of 0.17 dex. Our results imply that the Balmer decrement accurately predicts the obscuration of the nebular lines and can be used to robustly calculate SFRs for star-forming galaxies at z~2 with SFRs up to $\sim200\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We also use our data to assess SFR indicators based on modeling the UV-to-mid-IR SEDs or by adding SFR(UV) and SFR(IR), for which the latter is based on the mid-IR only or on the full IR SED. All these SFRs show a poorer agreement with SFR(H$α$,H$β$) and in some cases large systematic biases are observed. Finally, we show that the SFR and dust attenuation derived from the UV-to-near-IR SED alone are unbiased when assuming a delayed exponentially declining star-formation history.

preprint2016arXiv

The Relation between Galaxy Structure and Spectral Type: Implications for the Buildup of the Quiescent Galaxy Population at 0.5<z<2.0

We present the relation between galaxy structure and spectral type, using a K-selected galaxy sample at 0.5<z<2.0. Based on similarities between the UV-to-NIR spectral energy distributions, we classify galaxies into 32 spectral types. The different types span a wide range in evolutionary phases, and thus -- in combination with available CANDELS/F160W imaging -- are ideal to study the structural evolution of galaxies. Effective radii (R_e) and Sersic parameters (n) have been measured for 572 individual galaxies, and for each type, we determine R_e at fixed stellar mass by correcting for the mass-size relation. We use the rest-frame U-V vs. V-J diagram to investigate evolutionary trends. When moving into the direction perpendicular to the star-forming sequence, in which we see the Halpha equivalent width and the specific star formation rate (sSFR) decrease, we find a decrease in R_e and an increase in n. On the quiescent sequence we find an opposite trend, with older redder galaxies being larger. When splitting the sample into redshift bins, we find that young post-starburst galaxies are most prevalent at z>1.5 and significantly smaller than all other galaxy types at the same redshift. This result suggests that the suppression of star formation may be associated with significant structural evolution at z>1.5. At z<1, galaxy types with intermediate sSFRs (10^{-11.5}-10^{-10.5} yr^-1) do not have post-starburst SED shapes. These galaxies have similar sizes as older quiescent galaxies, implying that they can passively evolve onto the quiescent sequence, without increasing the average size of the quiescent galaxy population.

preprint2016arXiv

UV to IR Luminosities and Dust Attenuation Determined from ~4000 K-Selected Galaxies at 1<z<3 in the ZFOURGE Survey

We build a set of composite galaxy SEDs by de-redshifting and scaling multi-wavelength photometry from galaxies in the ZFOURGE survey, covering the CDFS, COSMOS, and UDS fields. From a sample of ~4000 K_s-band selected galaxies, we define 38 composite galaxy SEDs that yield continuous low-resolution spectra (R~45) over the rest-frame range 0.1-4 um. Additionally, we include far infrared photometry from the Spitzer Space Telescope and the Herschel Space Observatory to characterize the infrared properties of our diverse set of composite SEDs. From these composite SEDs we analyze the rest-frame UVJ colors, as well as the ratio of IR to UV light (IRX) and the UV slope ($β$) in the IRX$-β$ dust relation at 1<z<3. Blue star-forming composite SEDs show IRX and $β$ values consistent with local relations; dusty star-forming galaxies have considerable scatter, as found for local IR bright sources, but on average appear bluer than expected for their IR fluxes. We measure a tight linear relation between rest-frame UVJ colors and dust attenuation for star-forming composites, providing a direct method for estimating dust content from either (U-V) or (V-J) rest-frame colors for star-forming galaxies at intermediate redshifts.

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

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 Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

preprint2015arXiv

The MOSDEF Survey: Dissecting the star-formation rate vs. stellar mass relation using H$α$ and H$β$ emission lines at z ~ 2

We present results on the star-formation rate (SFR) versus stellar mass ($M_*$) relation (i.e., the "main sequence") among star-forming galaxies at 1.37$\leq$$z$$\leq$2.61 using the MOSFIRE Deep Evolution Field (MOSDEF) survey. Based on a sample of 261 galaxies with H$α$ and H$β$ spectroscopy, we have estimated robust dust-corrected instantaneous SFRs over a large range in $M_*$ ($\sim10^{9.5}-10^{11.5}M_\odot$). We find a correlation between log(SFR(H$α$)) and log($M_*$) with a slope of 0.65$\pm$0.08 (0.58$\pm$0.10) at 1.4<z<2.6 (2.1<z<2.6). We find that different assumptions for the dust correction, such as using the color-excess of the stellar continuum to correct the nebular lines, sample selection biases against red star-forming galaxies, and not accounting for Balmer absorption can yield steeper slopes of the log(SFR)-log($M_*$) relation. Our sample is immune from these biases as it is rest-frame optically selected, H$α$ and H$β$ are corrected for Balmer absorption, and the H$α$ luminosity is dust-corrected using the nebular color-excess computed from the Balmer decrement. The scatter of the log(SFR(H$α$))-log($M_*$) relation, after accounting for the measurement uncertainties, is 0.31 dex at 2.1<z<2.6, which is 0.05 dex larger than the scatter in log(SFR(UV))-log($M_*$). Based on comparisons to a simulated SFR-$M_*$ relation with some intrinsic scatter, we argue that in the absence of direct measurements of galaxy-to-galaxy variations in the attenuation/extinction curves and the IMF, one cannot use the difference in the scatter of the SFR(H$α$)- and SFR(UV)-$M_*$ relations to constrain the stochasticity of star formation in high-redshift galaxies.

preprint2015arXiv

The MOSDEF Survey: Electron Density and Ionization Parameter at $z\sim2.3$

Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the physical conditions of star-forming regions in $z\sim2.3$ galaxies, specifically the electron density and ionization state. From measurements of the [O II]$λλ$3726,3729 and [S II]$λλ$6716,6731 doublets, we find a median electron density of $\sim250$ cm$^{-3}$ at $z\sim2.3$, an increase of an order of magnitude compared to measurements of galaxies at $z\sim0$. While $z\sim2.3$ galaxies are offset towards significantly higher O$_{32}$ values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O$_{32}$ vs. R$_{23}$ and O3N2 diagrams. Based on these results, we propose that $z\sim2.3$ star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and $z\sim2.3$ galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from $z\sim0$ to $z\sim2.3$. We find that $z\sim2.3$ galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving [N II]$λ$6584. We conclude that the offset of $z\sim2.3$ galaxies from the local star-forming sequence in the [N II] BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at $z\sim0$ and $z\sim2$.

preprint2015arXiv

The MOSDEF Survey: Excitation Properties of $z\sim 2.3$ Star-forming Galaxies

We present results on the excitation properties of z~2.3 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) Survey. With its coverage of the full suite of strong rest-frame optical emission lines, MOSDEF provides an unprecedented view of the rest-frame optical spectra of a representative sample of distant star-forming galaxies. We investigate the locations of z~2.3 MOSDEF galaxies in multiple emission-line diagnostic diagrams. These include the [OIII]/Hb vs. [NII]/Ha and [OIII]/Hb vs. [SII]/Ha "BPT" diagrams, as well as the O_32 vs. R_23 excitation diagram. We recover the well-known offset in the star-forming sequence of high-redshift galaxies in the [OIII]/Hb vs. [NII]/Ha BPT diagram relative to SDSS star-forming galaxies. However, the shift for our rest-frame optically selected sample is less significant than for rest-frame-UV selected and emission-line selected galaxies at z~2. Furthermore, we find that the offset is mass-dependent, only appearing within the low-mass half of the z~2.3 MOSDEF sample, where galaxies are shifted towards higher [NII]/Ha at fixed [OIII]/Hb. Within the [OIII]/Hb vs. [SII]/Ha and O_32 vs. R_23 diagrams, we find that z~2.3 galaxies are distributed like local ones, and therefore attribute the shift in the [OIII]/Hb vs. [NII]/Ha BPT diagram to elevated N/O abundance ratios among lower-mass (M_*<10^10 M_sun) high-redshift galaxies. The variation in N/O ratios calls into question the use at high redshift of oxygen abundance indicators based on nitrogen lines, but the apparent invariance with redshift of the excitation sequence in the O_32 vs. R_23 diagram paves the way for using the combination of O_32 and R_23 as an unbiased metallicity indicator over a wide range in redshift. This indicator will allow for an accurate characterization of the shape and normalization of the mass-metallicity relationship over more than 10 Gyr.

preprint2015arXiv

The MOSDEF Survey: Measurements of Balmer Decrements and the Dust Attenuation Curve at Redshifts z~1.4-2.6

We present results on the dust attenuation curve of z~2 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies with nebular spectroscopic redshifts in the range z= 1.36-2.59 and high S/N measurements of, or upper limits on, the H-alpha and H-beta emission lines obtained with Keck/MOSFIRE. We construct composite SEDs of galaxies in bins of specific SFR and Balmer optical depth in order to directly constrain the dust attenuation curve from the UV through near-IR for typical star-forming galaxies at high redshift. Our results imply an attenuation curve that is very similar to the SMC extinction curve at wavelengths redward of 2500 Angstroms. At shorter wavelengths, the shape of the curve is identical to that of the Calzetti relation, but with a lower normalization (R_V). Hence, the new attenuation curve results in SFRs that are ~20% lower, and log stellar masses that are 0.16 dex lower, than those obtained with the Calzetti attenuation curve. Moreover, we find that the difference in the reddening---and the total attenuation---of the ionized gas and stellar continuum correlates strongly with SFR, such that for dust-corrected SFRs larger than 20 Msun/yr assuming a Chabrier IMF, the nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple model that can account for these trends is one in which the UV through optical stellar continuum is dominated by a population of less reddened stars, while the nebular line and bolometric luminosities become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based SFRs may underestimate the total SFR at even modest levels of ~20 Msun/yr. [Abridged]

preprint2015arXiv

The MOSDEF Survey: Optical AGN Diagnostics at z~2.3

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on rest-frame optical AGN identification and completeness at z~2.3. With our sample of 50 galaxies and 10 X-ray and IR-selected AGN with measured H-beta, [OIII], H-alpha, and [NII] emission lines, we investigate the location of AGN in the BPT, MEx (mass-excitation), and CEx (color-excitation) diagrams. We find that the BPT diagram works well to identify AGN at z~2.3 and that the z~0 AGN/star-forming galaxy classifications do not need to shift substantially at z~2.3 to robustly separate these populations. However, the MEx diagram fails to identify all of the AGN identified in the BPT diagram, and the CEx diagram is substantially contaminated at high redshift. We further show that AGN samples selected using the BPT diagram have selection biases in terms of both host stellar mass and stellar population, in that AGN in low mass and/or high specific star formation rate galaxies are difficult to identify using the BPT diagram. These selection biases become increasingly severe at high redshift, such that optically-selected AGN samples at high redshift will necessarily be incomplete. We also find that the gas in the narrow-line region appears to be more enriched than gas in the host galaxy for at least some MOSDEF AGN. However, AGN at z~2 are generally less enriched than local AGN with the same host stellar mass.

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

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

Simultaneous Modeling of the Stellar and Dust Emission in Distant Galaxies: Implications for Star Formation Rate Measurements

We have used near-ultraviolet (NUV) to mid-infrared (MIR) composite spectral energy distributions (SEDs) to simultaneously model the attenuated stellar and dust emission of 0.5 < z < 2.0 galaxies. These composite SEDs were previously constructed from the photometric catalogs of the NEWFIRM Medium-Band Survey, by stacking the observed photometry of galaxies that have similar rest-frame NUV-to-NIR SEDs. In this work, we include a stacked MIPS 24 micron measurement for each SED type to extend the SEDs to rest-frame MIR wavelengths. Consistent with previous studies, the observed MIR emission for most SED types is higher than expected from only the attenuated stellar emission. We fit the NUV-to-MIR composite SEDs by the Flexible Stellar Population Synthesis (SPS) models, which include both stellar and dust emission. We compare the best-fit star formation rates (SFRs) to the SFRs based on simple UV+IR estimators. Interestingly, the UV and IR luminosities overestimate SFRs - compared to the model SFRs - by more than ~ 1 dex for quiescent galaxies, while for the highest star-forming galaxies in our sample the two SFRs are broadly consistent. The difference in specific SFRs also shows a gradually increasing trend with declining specific SFR, implying that quiescent galaxies have even lower specific SFRs than previously found. Contributions from evolved stellar populations to both the UV and the MIR SEDs most likely explain the discrepancy. Based on this work, we conclude that SFRs should be determined from modeling the attenuated stellar and dust emission simultaneously, instead of employing simple UV+IR-based SFR estimators.

preprint2014arXiv

Synthetic Galaxy Images and Spectra from the Illustris Simulation

We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift z = 0 galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.

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 MOSDEF Survey: Mass, Metallicity, and Star-formation Rate at z~2.3

We present results on the z~2.3 mass-metallicity relation (MZR) using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We use an initial sample of 87 star-forming galaxies with spectroscopic coverage of Hβ, [OIII]λ5007, Hα, and [NII]λ6584 rest-frame optical emission lines, and estimate the gas-phase oxygen abundance based on the N2 and O3N2 strong-line indicators. We find a positive correlation between stellar mass and metallicity among individual z~2.3 galaxies using both the N2 and O3N2 indicators. We also measure the emission-line ratios and corresponding oxygen abundances for composite spectra in bins of stellar mass. Among composite spectra, we find a monotonic increase in metallicity with increasing stellar mass, offset ~0.15-0.3 dex below the local MZR. When the sample is divided at the median star-formation rate (SFR), we do not observe significant SFR dependence of the z~2.3 MZR among either individual galaxies or composite spectra. We furthermore find that z~2.3 galaxies have metallicities ~0.1 dex lower at a given stellar mass and SFR than is observed locally. This offset suggests that high-redshift galaxies do not fall on the local "fundamental metallicity relation" among stellar mass, metallicity, and SFR, and may provide evidence of a phase of galaxy growth in which the gas reservoir is built up due to inflow rates that are higher than star-formation and outflow rates. However, robust conclusions regarding the gas-phase oxygen abundances of high-redshift galaxies await a systematic reappraisal of the application of locally calibrated metallicity indicators at high redshift.

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

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

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

preprint2013arXiv

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 Dust Attenuation Law in Distant Galaxies: Evidence for Variation with Spectral Type

This letter utilizes composite spectral energy distributions (SEDs) constructed from NEWFIRM Medium-Band Survey photometry to constrain the dust attenuation curve in 0.5<z<2.0 galaxies. Based on similarities between the full SED shapes (0.3-8 micron), we have divided galaxies in 32 different spectral classes and stacked their photometry. As each class contains galaxies over a range in redshift, the resulting rest-frame SEDs are well-sampled in wavelength and show various spectral features including Halpha and the UV dust bump at 2175 Angstrom. We fit all composite SEDs with flexible stellar population synthesis models, while exploring attenuation curves with varying slopes and UV bump strengths. The Milky Way and Calzetti law provide poor fits at UV wavelengths for nearly all SEDs. Consistent with previous studies, we find that the best-fit attenuation law varies with spectral type. There is a strong correlation between the best-fit dust slope and UV bump strength, with steeper laws having stronger bumps. Moreover, the attenuation curve correlates with specific star formation rate (SFR), with more active galaxies having shallower dust curves and weaker bumps. There is also a weak correlation with inclination. The observed trends can be explained by differences in the dust-to-star geometry, a varying grain size distribution, or a combination of both. Our results have several implications for galaxy evolution studies. First, the assumption of a universal dust model leads to biases in derived galaxy properties. Second, the presence of a dust bump may result in underestimated values for the UV slope, used to correct SFRs of distant galaxies.

preprint2013arXiv

The Radial Distribution of Star Formation in Galaxies at z~1 from the 3D-HST Survey

The assembly of galaxies can be described by the distribution of their star formation as a function of cosmic time. Thanks to the WFC3 grism on HST it is now possible to measure this beyond the local Universe. Here we present the spatial distribution of Halpha emission for a sample of 54 strongly star-forming galaxies at z~1 in the 3D-HST Treasury survey. By stacking the Halpha emission we find that star formation occurred in approximately exponential distributions at z~1, with median Sersic index of n=1.0+-0.2. The stacks are elongated with median axis ratios of b/a=0.58+-0.09 in Halpha, consistent with (possibly thick) disks at random orientation angles. Keck spectra obtained for a subset of eight of the galaxies show clear evidence for rotation, with inclination-corrected velocities of 90 to 330 km/s. The most straightforward interpretation of our results is that star formation in strongly star-forming galaxies at z~1 generally occurred in disks. The disks appear to be "scaled-up" versions of nearby spiral galaxies: they have EW(Halpha)~100 Angstroms out to the solar orbit and they have star formation surface densities above the threshold for driving galactic scale winds.

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 "mass FP" 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: A wide-field grism spectroscopic survey with the Hubble Space Telescope

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

preprint2012arXiv

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

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

preprint2012arXiv

Large-Scale Star Formation-Driven Outflows at 1<z<2 in the 3D-HST Survey

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

preprint2012arXiv

Massive and Newly Dead: Discovery of a Significant Population of Galaxies with High Velocity Dispersions and Strong Balmer Lines at z~1.5 from Deep Keck Spectra and HST/WFC3 Imaging

We present deep Keck/LRIS spectroscopy and HST/WFC3 imaging in the rest-frame optical for a sample of eight galaxies at z~1.5 with high photometrically-determined stellar masses. The data are combined with VLT/XShooter spectra of five galaxies from van de Sande et al. (2011, 2012 to be submitted). We find that these thirteen galaxies have high velocity dispersions, with a median of sigma=301 km s^{-1}. This high value is consistent with their relatively high stellar masses and compact sizes. We study their stellar populations using the strength of Balmer absorption lines, which are not sensitive to dust absorption. We find a large range in Balmer absorption strength, with many galaxies showing very strong lines indicating young ages. The median Hdelta_A equivalent width, determined directly or inferred from the H10 line, is 5.4 Angstroms, indicating a luminosity-weighted age of ~1 Gyr. Although this value may be biased towards higher values because of selection effects,high-dispersion galaxies with such young ages are extremely rare in the local Universe. Interestingly we do not find a simple correlation with rest-frame U-V color: some of the reddest galaxies have very strong Balmer absorption lines. These results demonstrate that many high-dispersion galaxies at z~1.5 were quenched recently. This implies that there must be a population of star-forming progenitors at z~2 with high velocity dispersions or linewidths, which are notoriously absent from CO/Halpha selected surveys.

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.

preprint2011arXiv

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

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

preprint2011arXiv

Evolution of Galaxy Stellar Mass Functions, Mass Densities, and Mass to Light Ratios from z 7 to z 4

We derive stellar masses from SED fitting to rest-frame optical and UV fluxes for 401 star-forming galaxies at z 4, 5, and 6 from Hubble-WFC3/IR observations of the ERS combined with the deep GOODS-S Spitzer/IRAC data (and include a previously-published z 7 sample). A mass-luminosity relation with strongly luminosity-dependent M/Luv ratios is found for the largest sample (299 galaxies) at z 4. The relation M \propto L_{UV,1500}^(1.7+/-0.2) has a well-determined intrinsic sample variance of 0.5 dex. This relation is also consistent with the more limited samples at z 5-7. This z 4 mass-luminosity relation, and the well-established faint UV luminosity functions at z 4-7, are used to derive galaxy mass functions (MF) to masses M~10^8 at z 4-7. A bootstap approach is used to derive the MFs to account for the large scatter in the M--Luv relation and the luminosity function uncertainties, along with an analytical crosscheck. The MFs are also corrected for the effects of incompleteness. The incompleteness-corrected MFs are steeper than previously found, with slopes α_M-1.4 to -1.6 at low masses. These slopes are, however, still substantially flatter than the MFs obtained from recent hydrodynamical simulations. We use these MFs to estimate the stellar mass density (SMD) of the universe to a fixed M_{UV,AB}<-18 as a function of redshift and find a SMD growth \propto(1+z)^{-3.4 +/-0.8} from z 7 to z 4. We also derive the SMD from the completeness-corrected MFs to a mass limit M~10^{8} Msun. Such completeness-corrected MFs and the derived SMDs will be particularly important for model comparisons as future MFs reach to lower masses.

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

Redshift Evolution of the Galaxy Velocity Dispersion Function

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

The evolution of the mass-size relation to z=3.5 for UV-bright galaxies and sub-mm galaxies in the GOODS-NORTH field

We study the evolution of the size - stellar mass relation for a large spectroscopic sample of galaxies in the GOODs North field up to $z \sim 3.5$. The sizes of the galaxies are measured from $\textit{K}_{s}$-band images (corresponding to rest-frame optical/NIR) from the Subaru 8m telescope. We reproduce earlier results based on photometric redshifts that the sizes of galaxies at a given mass evolve with redshift. Specifically, we compare sizes of UV-bright galaxies at a range of redshifts: Lyman break galaxies (LBGs) selected through the U-drop technique ($z \sim 2.5-3.5$), BM/BX galaxies at $z \sim 1.5-2.5$, and GALEX LBGs at low redshift ($z \sim 0.6-1.5$). The median sizes of these UV-bright galaxies evolve as $(1+z)^{-1.11\pm0.13}$ between $z \sim 0.5-3.5$. The UV-bright galaxies are significantly larger than quiescent galaxies at the same mass and redshift by $0.45\pm0.09$ dex. We also verify the correlation between color and stellar mass density of galaxies to high redshifts. The sizes of sub-mm galaxies in the same field are measured and compared with BM/BX galaxies. We find that median half-light radii of SMGs is $2.90 \pm 0.45$ kpc and there is little difference in their size distribution to the UV-bright star forming galaxies.

preprint2010arXiv

The evolving relations between size, mass, surface density, and star formation in 3x10^4 galaxies since z=2

The presence of massive, compact, quiescent galaxies at z>2 presents a major challenge for theoretical models of galaxy formation and evolution. Using one of the deepest large public near-IR surveys to date, we investigate in detail the correlations between star formation and galaxy structural parameters (size, stellar mass, and surface density) from z=2 to the present. At all redshifts, massive quiescent galaxies (i.e. those with little or no star formation) occupy the extreme high end of the surface density distribution and follow a tight mass-size correlation, while star-forming galaxies show a broad range of both densities and sizes. Conversely, galaxies with the highest surface densities comprise a nearly-homogeneous population with little or no ongoing star formation, while less dense galaxies exhibit high star-formation rates and varying levels of dust obscuration. Both the sizes and surface densities of quiescent galaxies evolve strongly from z=2-0; we parameterize this evolution for both populations with simple power law functions and present best-fit parameters for comparison to future theoretical models. Higher-mass quiescent galaxies undergo faster structural evolution, consistent with previous results. Interestingly, star-forming galaxies' sizes and densities evolve at rates similar to those of quiescent galaxies. It is therefore possible that the same physical processes drive the structural evolution of both populations, suggesting that "dry mergers" may not be the sole culprit in this size evolution.

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

preprint2010arXiv

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

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

preprint2010arXiv

The Stellar Mass Density and Specific Star Formation Rate of the Universe at z~7

We use a robust sample of 11 z~7 galaxies (z-dropouts) to estimate the stellar mass density of the universe when it was only ~750 Myr old. We combine the very deep optical to near-Infrared photometry from the HST ACS and NICMOS cameras with mid-Infrared Spitzer IRAC imaging available through the GOODS program. After carefully removing the flux from contaminating foreground sources we have obtained reliable photometry in the 3.6 and 4.5 micron IRAC channels. The spectral shapes of these sources, including their rest frame optical colors, strongly support their being at z~7 with a mean photometric redshift of <z>=7.2+/-0.5. We use Bruzual & Charlot (2003) synthetic stellar population models to constrain their stellar masses and star formation histories. We find stellar masses that range over 0.1 -12x10^9 M_sol and average ages from 20 Myr to up to 425 Myr with a mean of ~300 Myr, suggesting that in some of these galaxies most of the stars were formed at z>8 (and probably at z>~10). The best fits to the observed SEDs are consistent with little or no dust extinction, in agreement with recent results at z~4-8. The star formation rates (SFR) are in the range from 5-20 M_sol/yr. From this sample we measure a stellar mass density of 6.6_{-3.3}^{+5.4}x10^5 M_sol/Mpc^3 to a limit of M_{UV,AB}<-20 (or 0.4 L*(z=3)). Combined with a fiducial lower limit for their ages (80 Myr) this implies a maximum SFR density of 0.008 M_sol/yr/Mpc^3. This is well below the critical level needed to reionize the universe at z~8 using standard assumptions. However, this result is based on luminous sources (>L*) and does not include the dominant contribution of the fainter galaxies. Strikingly, we find that the specific SFR is constant from z~7 to z~2 but drops substantially at more recent times.

preprint2010arXiv

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

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

preprint2009arXiv

A 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

A Spatially Resolved Map of the Kinematics, Star-Formation and Stellar Mass Assembly in a Star-Forming Galaxy at z=4.9

We present a detailed study of the spatially resolved kinematics, star-formation and stellar mass in a highly amplified galaxy at z=4.92 behind the lensing cluster MS1358+62. We use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the velocity structure of the galaxy from the nebular [OII] emission. Using a detailed gravitational lens model, we account for lensing amplification factor 12.+/-2.0 and find that this intrinsically L* galaxy has a stellar mass of M*=7+/-2x10^8Mo, a dynamical mass of Mdyn=3+/-1x10^9csc^2(i)Mo (within of 2kpc) and a star-formation rate of 42+/-8Mo/yr. The source-plane UV/optical morphology of this galaxy is dominated by five discrete star-forming regions. Exploiting the dynamical information we derive masses for individual star-forming regions of Mcl~10^(8-9)Mo with sizes of ~200pc. We find that, at a fixed size, the star-formation rate density within these HII regions is approximately two orders of magnitude greater than those observed in local spiral/starburst galaxies, but consistent with the most massive HII regions in the local Universe such as 30Dor. Finally, we compare the spatially resolved nebular emission line velocity with the Ly-alpha and UV-ISM lines and find that this galaxy is surrounded by a galactic scale outflow in which the Ly-alpha appears redshifted by ~150km/s and the UV-ISM lines blue-shifted by -200km/s from the (systemic) nebular emission. The velocity structure of the outflow mirrors that of the nebular emission suggesting the outflow is young (~15Myr), and has yet to burst out of the system. Taken together, these results suggest that this young galaxy is undergoing its first major epoch of mass assembly.

preprint2009arXiv

How Massive are Massive Compact Galaxies?

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

preprint2009arXiv

The Growth of Massive Galaxies Since z=2

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

preprint2009arXiv

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

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

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.

preprint2008arXiv

Confirmation of the remarkable compactness of massive quiescent galaxies at z~2.3: early-type galaxies did not form in a simple monolithic collapse

Using deep near-infrared spectroscopy Kriek et al. (2006) found that ~45% of massive galaxies at z~2.3 have evolved stellar populations and little or no ongoing star formation. Here we determine the sizes of these quiescent galaxies using deep, high-resolution images obtained with HST/NIC2 and laser guide star-assisted Keck/AO. Considering that their median stellar mass is 1.7x10^11 Solar masses the galaxies are remarkably small, with a median effective radius of 0.9 kpc. Galaxies of similar mass in the nearby Universe have sizes of ~5 kpc and average stellar densities which are two orders of magnitude lower than the z~2.3 galaxies. These results extend earlier work at z~1.5 and confirm previous studies at z>2 which lacked spectroscopic redshifts and imaging of sufficient resolution to resolve the galaxies. Our findings demonstrate that fully assembled early-type galaxies make up at most ~10% of the population of K-selected quiescent galaxies at z~2.3, effectively ruling out simple monolithic models for their formation. The galaxies must evolve significantly after z~2.3, through dry mergers or other processes, consistent with predictions from hierarchical models.

preprint2006arXiv

Spectroscopic Identification of Massive Galaxies at z~2.3 with Strongly Suppressed Star Formation

We present first results of a spectroscopic survey targeting K-selected galaxies at z=2.0-2.7 using the GNIRS instrument on Gemini-South. We obtained near-infrared spectra with a wavelength coverage of 1.0-2.5 micron for 26 K-bright galaxies (K<19.7) selected from the MUSYC survey using photometric redshifts. We successfully derived spectroscopic redshifts for all 26 galaxies using rest-frame optical emission lines or the redshifted Balmer/4000 Angstrom break. Twenty galaxies have spectroscopic redshifts in the range 2.0<z<2.7, for which bright emission lines like Halpha and [OIII] fall in atmospheric windows. Surprisingly, we detected no emission lines for nine of these 20 galaxies. The median 2 sigma upper limit on the rest-frame equivalent width of Halpha for these nine galaxies is ~10 Angstrom. The stellar continuum emission of these same nine galaxies is best fitted by evolved stellar population models. The best-fit star formation rate (SFR) is zero for five out of nine galaxies, and consistent with zero within 1 sigma for the remaining four. Thus, both the Halpha measurements and the independent stellar continuum modeling imply that 45% of our K-selected galaxies are not forming stars intensely. This high fraction of galaxies without detected line emission and low SFRs may imply that the suppression of star formation in massive galaxies occurs at higher redshift than is predicted by current CDM galaxy formation models. However, obscured star formation may have been missed, and deep mid-infrared imaging is needed to clarify this situation.