Source author record

Pieter van Dokkum

Pieter van Dokkum 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

102works
6topics
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

102 published item(s)

preprint2025arXiv

Kinematic Confirmation of a Remarkable Linear Trail of Galaxies in the NGC 1052 Field, Consistent with Formation in a High-Speed Bullet Dwarf Collision

A unique linear trail of diffuse galaxies was recently identified in the NGC 1052 field. This trail includes the remarkable, ultra-diffuse galaxies DF2 and DF4 which lack dark matter and host unusually luminous globular clusters. It has been proposed that the trail formed via a high-speed collision between two gas-rich dwarf galaxies. This scenario predicts that the trail galaxies are kinematically connected and follow a specific trend in radial velocity as a function of position, based on the known velocities and positions of DF2 and DF4. To test this hypothesis, we measured radial velocities for seven additional galaxies on the trail. While the galaxies' low surface brightnesses presented observational challenges, we employ several methods to obtain measurements for galaxies with effective surface brightnesses up to 28.6 mag arcsec$^{-2}$, including a narrow slit placed over globular clusters and a novel wide slit mode on Keck/LRIS, as well as a 'light bucket' mode on Keck/KCWI. We find that five of our seven targets follow the precise velocity trend predicted by DF2 and DF4, to a degree with just a 2% chance of randomly occurring. Moreover, the trail galaxies' radial velocities are significantly higher than those of the NGC 1052 group, setting it apart as a separate, kinematically connected system. Our findings support the theory that this trail of galaxies, including DF2 and DF4, formed together in a single event. A 'bullet dwarf' collision remains the only known explanation for all the unusual properties of DF2, DF4, and the associated trail of galaxies.

preprint2025arXiv

Quantitative Morphology of Galactic Cirrus in Deep Optical Imaging

Imaging of optical Galactic cirrus, the spatially resolved form of diffuse Galactic light, provides important insights into the properties of the diffuse interstellar medium (ISM) in the Milky Way. While previous investigations have focused mainly on the intensity characteristics of optical cirrus, their morphological properties remain largely unexplored. In this study, we employ several complementary statistical approaches -- local intensity statistics, angular power spectrum / $Δ$-variance analysis, and wavelet scattering transform analysis -- to characterize the morphology of cirrus in deep optical imaging data. We place our investigation of optical cirrus into a multi-wavelength context by comparing the morphology of cirrus seen with the Dragonfly Telephoto Array to that seen with space-based facilities working at longer wavelengths (Herschel 250 $μm$, WISE 12 $μm$, and Planck radiance), as well as with structures seen in the DHIGLS HI column density map. Our statistical methods quantify the similarities and the differences of cirrus morphology in all these datasets. The morphology of cirrus at visible wavelengths resembles that of far-infrared cirrus more closely than that of mid-infrared cirrus; on small scales, anisotropies in the cosmic infrared background and systematics may lead to differences. Across all dust tracers, cirrus morphology can be well described by a power spectrum with a common power-law index $γ\sim-2.9$. We demonstrate quantitatively that optical cirrus exhibits filamentary, coherent structures across a broad range of angular scales. Our results offer promising avenues for linking the analysis of coherent structures in optical cirrus to the underlying physical processes in the ISM that shape them. Furthermore, we demonstrate that these morphological signatures can be leveraged to distinguish and disentangle cirrus from extragalactic light.

preprint2023arXiv

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

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

preprint2023arXiv

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

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

preprint2022arXiv

A giant shell of ionized gas discovered near M82 with the Dragonfly Spectral Line Mapper pathfinder

We present the discovery of a giant cloud of ionized gas in the field of the starbursting galaxy M82. Emission from the cloud is seen in H$α$ and [NII]$λ$6583 in data obtained though a small pathfinder instrument used to test the key ideas that will be implemented in the Dragonfly Spectral Line Mapper, an upcoming ultranarrow-bandpass imaging version of the Dragonfly Telephoto Array. The discovered cloud has a shell-like morphology with a linear extent of $0.8^{\circ}$ and is positioned $0.6^{\circ}$ northwest of M82. At the heliocentric distance of the M81 group, the cloud's extent corresponds to 55 kpc and its projected distance from the nucleus of M82 is 40 kpc. The cloud has an average H$α$ surface brightness of $2\times10^{-18}$ $\mathrm{erg} ~\mathrm{cm}^{-2} ~\mathrm{s}^{-1} ~\mathrm{arcsec}^{-2}$. The [NII]$λ$6583/H$α$ line ratio varies from [NII]/H$α$ $\sim0.2$ to [NII]/H$α$ $\sim1.0$ across the cloud, with higher values found in its eastern end. Follow-up spectra obtained with Keck-LRIS confirm the existence of the cloud and yield line ratios of [NII]$λ$6583/H$α$ = 0.340 $\pm$ 0.003 and [SII]$λ$6716,6731/H$α$ = 0.64 $\pm$ 0.03 in the cloud. This giant cloud of material could be lifted from M82 by tidal interactions or by its powerful starburst. Alternatively, it may be gas infalling from the cosmic web, potentially precipitated by the superwinds of M82. Deeper data are needed to test these ideas further. The upcoming Dragonfly Spectral Line Mapper will have 120 lenses, $40\times$ more than in the pathfinder instrument used to obtain the data presented here.

preprint2022arXiv

A trail of dark matter-free galaxies from a bullet dwarf collision

The ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group share several unusual properties: they both have large sizes, rich populations of overluminous and large globular clusters, and very low velocity dispersions indicating little or no dark matter. It has been suggested that these galaxies were formed in the aftermath of high velocity encounters of gas rich galaxies, events that resemble the collision that created the bullet cluster but on much smaller scales. The gas separates from the dark matter in the collision and subsequent star formation leads to the formation of one or more dark matter-free galaxies. Here we show that the present-day line-of-sight distances and radial velocities of DF2 and DF4 are consistent with their joint formation in the aftermath of a single bullet-dwarf collision, around eight billion years ago. Moreover, we find that DF2 and DF4 are part of an apparent linear substructure of 7-11 large, low-luminosity objects. We propose that these all originated in the same event, forming a trail of dark matter-free galaxies that is more than 2 Mpc long and angled 7 +- 2 degrees from the line of sight. We also tentatively identify the highly dark matter-dominated remnants of the two progenitor galaxies that are expected at the leading edges of the trail.

preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

preprint2022arXiv

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

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

preprint2022arXiv

Light from the Darkness: Detecting Ultra-Diffuse Galaxies in the Perseus Cluster through Over-densities of Globular Clusters with a Log-Gaussian Cox Process

We introduce a new method for detecting ultra-diffuse galaxies by searching for over-densities in intergalactic globular cluster populations. Our approach is based on an application of the log-Gaussian Cox process, which is a commonly used model in the spatial statistics literature but rarely used in astronomy. This method is applied to the globular cluster data obtained from the PIPER survey, a \textit{Hubble Space Telescope} imaging program targeting the Perseus cluster. We successfully detect all confirmed ultra-diffuse galaxies with known globular cluster populations in the survey. We also identify a potential galaxy that has no detected diffuse stellar content. Preliminary analysis shows that it is unlikely to be merely an accidental clump of globular clusters or other objects. If confirmed, this system would be the first of its kind. Simulations are used to assess how the physical parameters of the globular cluster systems within ultra-diffuse galaxies affect their detectability using our method. We quantify the correlation of the detection probability with the total number of globular clusters in the galaxy and the anti-correlation with increasing half-number radius of the globular cluster system. The Sérsic index of the globular cluster distribution has little impact on detectability.

preprint2022arXiv

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

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

preprint2022arXiv

Stellar halos from The Dragonfly Edge-on Galaxies Survey

We present the primary results from the Dragonfly Edge-on Galaxies Survey (DEGS), an exploration of the stellar halos of twelve nearby ($d < 25$ Mpc) edge-on disc galaxies with the Dragonfly Telephoto Array. The edge-on orientation of these galaxies allows their stellar halos to be explored with minimal obscuration by or confusion with the much brighter disc light. Galaxies in the sample span a range of stellar masses from $10^{9.68} - 10^{10.88} M_\odot$. We confirm that the wide range of stellar halo mass fractions previously seen for Milky Way-mass galaxies is also found among less massive spiral galaxies. The scatter in stellar halo mass fraction is large but we do find a significant positive correlation between stellar halo mass fraction and total stellar mass when the former is measured beyond five half-mass radii. Reasonably good agreement is found with predictions from cosmological hydrodynamical simulations, although observed stellar halo fractions appear to be somewhat lower than expected from these simulations.

preprint2022arXiv

The Star Formation Burstiness and Ionizing Efficiency of Low-mass Galaxies

We investigate the burstiness of star formation and the ionizing efficiency of a large sample of galaxies at $0.7 < z < 1.5$ using HST grism spectroscopy and deep ultraviolet (UV) imaging in the GOODS-N and GOODS-S fields. The star formation history (SFH) in these strong emission line low-mass galaxies indicates an elevated star formation rate (SFR) based on the H$α$ emission line at a given stellar mass when compared to the standard main sequence. Moreover, when comparing the H$α$ and UV SFR indicators, we find that an excess in SFR(H$α$) compared to SFR(UV) is preferentially observed in lower-mass galaxies below $10^{9}$ M$\odot$, which are also the highest-EW galaxies. These findings suggest that the burstiness parameters of these strong emission line galaxies may differ from those inferred from hydrodynamical simulations and previous observations. For instance, a larger burstiness duty cycle would explain the observed SFR(H$α$) excess. We also estimate the ionizing photon production efficiency $ξ_{ion}$, finding a median value of Log($ξ_{ion}$/erg$^{-1}$ Hz)$=24.80 \pm 0.26$ when adopting a Galactic dust correction for H$α$ and an SMC one for the stellar component. We observe an increase of $ξ_{ion}$ with redshift, further confirming similar results at higher redshifts. We also find that $ξ_{ion}$ is strongly correlated with EW(H$α$), which provides an approach for deriving $ξ_{ion}$ in early galaxies. Lower-mass, lower-luminosity galaxies have a higher $ξ_{ion}$. Overall, these results provide further support for faint galaxies playing a major role in the reionization of the Universe.

preprint2022arXiv

Tidal Distortions in NGC1052-DF2 and NGC1052-DF4: Independent Evidence for a Lack of Dark Matter

Two ultra diffuse galaxies in the same group, NGC1052-DF2 and NGC1052-DF4, have been found to have little or no dark matter and to host unusually luminous globular cluster populations. Such low mass diffuse objects in a group environment are easily disrupted and are expected to show evidence of tidal distortions. In this work we present deep new imaging of the NGC1052 group, obtained with the Dragonfly Telephoto Array, to test this hypothesis. We find that both galaxies show strong position angle twists and are significantly more elongated at their outskirts than in their interiors. The group's central massive elliptical NGC1052 is the most likely source of these tidal disturbances. The observed distortions imply that the galaxies have a low total mass or are very close to NGC1052. Considering constraints on the galaxies' relative distances, we infer that the dark matter halo masses of these galaxies cannot be much greater than their stellar masses. Calculating pericenters from the distortions, we find that the galaxies are on highly elliptical orbits, with a ratio of pericenter to present-day radius Rperi/R0~0.1 if the galaxies are dark matter-free and Rperi/R0~0.01 if they have a normal dark halo. Our findings provide strong evidence, independent of kinematic constraints, that both galaxies are dark matter deficient. Furthermore, the similarity of the tidal features in NGC1052-DF2 and NGC1052-DF4 strongly suggests that they arose at comparable distances from NGC1052. In Appendix A, we describe sbcontrast, a robust method to determine the surface brightness limit of images.

preprint2021arXiv

A Method To Characterize the Wide-Angle Point Spread Function of Astronomical Images

Uncertainty in the wide-angle Point Spread Function (PSF) at large angles (tens of arcseconds and beyond) is one of the dominant sources of error in a number of important quantities in observational astronomy. Examples include the stellar mass and shape of galactic halos and the maximum extent of starlight in the disks of nearby galaxies. However, modeling the wide-angle PSF has long been a challenge in astronomical imaging. In this paper, we present a self-consistent method to model the wide-angle PSF in images. Scattered light from multiple bright stars is fitted simultaneously with a background model to characterize the extended wing of the PSF using a Bayesian framework operating on pixel-by-pixel level. The method is demonstrated using our software elderflower and is applied to data from the Dragonfly Telephoto Array to model its PSF out to 20-25 arcminutes. We compare the wide-angle PSF of Dragonfly to that of a number of other telescopes, including the SDSS PSF, and show that on scales of arcminutes the scattered light in the Dragonfly PSF is markedly lower than that of other wide-field imaging telescopes. The energy in the wings of the Dragonfly point-spread function is sufficiently low that optical cleanliness plays an important role in defining the PSF. This component of the PSF can be modelled accurately, highlighting the power of our self-contained approach.

preprint2021arXiv

Blue Rest-Frame UV-Optical Colors in z~8 Galaxies from GREATS: Very Young Stellar Populations at ~650 Myr of Cosmic Time

Deep rest-optical observations are required to accurately constrain the stellar populations of $z\sim8$ galaxies. Due to significant limitations in the availability of such data for statistically complete samples, observational results have been limited to modest numbers of bright or lensed sources. To revolutionize the present characterization of $z\sim8$ galaxies, we exploit the ultradeep ($\sim27$ mag, $3σ$) Spitzer/IRAC $3.6μ$m and $4.5μ$m data, probing the rest-frame optical at $z\sim8$, over $\sim200$ arcmin$^2$ of the GOODS fields from the recently completed GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) program, combined with observations in the CANDELS UDS and COSMOS fields. We stacked $\gtrsim100$ $z\sim8$ Lyman-Break galaxies in four bins of UV luminosity ($M_\mathrm{UV}\sim -20.7$ to $-18.4$) and study their $H_\mathrm{160}-[3.6]$ and $[3.6]-[4.5]$ colors. We find young ages ($\lesssim100$ Myr) for the three faintest stacks, inferred from their blue $H_\mathrm{160}-[3.6]\sim 0$ mag colors, consistent with a negative Balmer break. Meanwhile, the redder $H_\mathrm{160}-[3.6]$ color seen in the brightest stack is suggestive of slightly older ages. We explored the existence of a correlation between the UV luminosity and age, and find either no trend or fainter galaxies being younger. The stacked SEDs also exhibit very red $[3.6]-[4.5]\sim0.5$ mag colors, indicative of intense [OIII]+H$β$ nebular emission and SFR. The correspondingly high specific star-formation rates, sSFR$\gtrsim10$Gyr$^{-1}$, are consistent with recent determinations at similar redshifts and higher luminosities, and support the co-evolution between the sSFR and the specific halo mass accretion rate.

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

Measuring distances to low-luminosity galaxies using surface brightness fluctuations

We present an in-depth study of surface brightness fluctuations (SBFs) in low-luminosity stellar systems. Using the MIST models, we compute theoretical predictions for absolute SBF magnitudes in the LSST, HST ACS/WFC, and proposed Roman Space Telescope filter systems. We compare our calculations to observed SBF-color relations of systems that span a wide range of age and metallicity. Consistent with previous studies, we find that single-age population models show excellent agreement with observations of low-mass galaxies with $0.5 \lesssim g - i \lesssim 0.9$. For bluer galaxies, the observed relation is better fit by models with composite stellar populations. To study SBF recovery from low-luminosity systems, we perform detailed image simulations in which we inject fully populated model galaxies into deep ground-based images from real observations. Our simulations show that LSST will provide data of sufficient quality and depth to measure SBF magnitudes with precisions of ${\sim}0.2$-0.5 mag in ultra-faint $\left(\mathrm{10^4 \leq M_\star/M_\odot \leq 10^5}\right)$ and low-mass classical (M$_\star\leq10^7$ M$_\odot$) dwarf galaxies out to ${\sim}4$ Mpc and ${\sim}25$ Mpc, respectively, within the first few years of its deep-wide-fast survey. Many significant practical challenges and systematic uncertainties remain, including an irreducible "sampling scatter" in the SBFs of ultra-faint dwarfs due to their undersampled stellar mass functions. We nonetheless conclude that SBFs in the new generation of wide-field imaging surveys have the potential to play a critical role in the efficient confirmation and characterization of dwarf galaxies in the nearby universe.

preprint2021arXiv

NGC5846-UDG1: A galaxy formed mostly by star formation in massive, extremely dense clumps of gas

It has been shown that ultra-diffuse galaxies (UDGs) have higher specific frequencies of globular clusters on average than other dwarf galaxies with similar luminosities. The UDG NGC5846-UDG1 is among the most extreme examples of globular cluster-rich galaxies found so far. Here we present new Hubble Space Telescope (HST) observations and analysis of this galaxy and its globular cluster system. We find that NGC5846-UDG1 hosts $54 \pm 9$ globular clusters, three to four times more than any previously known galaxy with a similar luminosity, and higher than reported in previous studies. With a galaxy luminosity of $L_{V,\mathrm{gal}} \approx 6 \times 10^7\,{\rm L}_{\odot}$ ($M_\star \approx 1.2 \times 10^8\,{\rm M}_\odot$) and a total globular cluster luminosity of $L_{V,\mathrm{GCs}} \approx 7.6 \times 10^6\,{\rm L}_{\odot}$, we find that the clusters currently comprise $\sim 13 \%$ of the total light. Taking into account the effects of mass loss from clusters during their formation and throughout their lifetime, we infer that most of the stars in the galaxy likely formed in globular clusters, and very little to no "normal" low-density star formation occurred. This result implies that the most extreme conditions during early galaxy formation promoted star formation in massive and dense clumps, in contrast to the dispersed star formation observed in galaxies today.

preprint2021arXiv

Spatially Resolved Stellar Spectroscopy of the Ultra-diffuse Galaxy Dragonfly 44. III. Evidence for an Unexpected Star-Formation History

We use the Keck Cosmic Web Imager integral-field unit spectrograph to: 1) measure the global stellar population parameters for the ultra-diffuse galaxy (UDG) Dragonfly 44 (DF44) to much higher precision than previously possible for any UDG, and 2) for the first time measure spatially-resolved stellar population parameters of a UDG. We find that DF44 falls below the mass--metallicity relation established by canonical dwarf galaxies both in and beyond the Local Group. We measure a flat radial age gradient ($m_{\rm age} \sim +0.01_{-0.08}^{+0.07}$ log Gyr kpc$^{-1}$) and a flat-to-positive metallicity gradient ($m_{\rm [Fe/H]} \sim +0.08_{-0.11}^{+0.11}$ dex kpc$^{-1}$), which are inconsistent with the gradients measured in similarly pressure-supported dwarf galaxies. We also measure a flat-to-negative [Mg/Fe] gradient ($m_{\rm [Mg/Fe]} \sim -0.18_{-0.17}^{+0.17}$ dex kpc$^{-1}$) such that the central $1.5$ kpc of DF44 has stellar population parameters comparable to metal-poor globular clusters. Overall, DF44 does not have internal properties similar to other dwarf galaxies and is inconsistent with it having been puffed up through a prolonged, bursty star-formation history, as suggested by some simulations. Rather, the evidence indicates that DF44 experienced an intense epoch of "inside-out" star formation and then quenched early and catastrophically, such that star-formation was cut off more quickly than in canonical dwarf galaxies.

preprint2021arXiv

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

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

preprint2020arXiv

A missing outskirts problem? Comparisons between stellar halos in the Dragonfly Nearby Galaxies Survey and the TNG100 simulation

Low surface brightness galactic stellar halos provide a challenging but promising path towards unraveling the past assembly histories of individual galaxies. Here, we present detailed comparisons between the stellar halos of Milky Way-mass disk galaxies observed as part of the Dragonfly Nearby Galaxies Survey (DNGS) and stellar mass-matched galaxies in the TNG100 run of the IllustrisTNG project. We produce stellar mass maps as well as mock $g$ and $r$-band images for randomly oriented simulated galaxies, convolving the latter with the Dragonfly PSF and taking care to match the background noise, surface brightness limits and spatial resolution of DNGS. We measure azimuthally averaged stellar mass density and surface brightness profiles, and find that the DNGS galaxies generally have less stellar mass (or light) at large radii (>20 kpc) compared to their mass-matched TNG100 counterparts, and that simulated galaxies with similar surface density profiles tend to have low accreted mass fractions for their stellar mass. We explore potential solutions to this apparent "missing outskirts problem" by implementing several ad-hoc adjustments within TNG100 at the stellar particle level. Although we are unable to identify any single adjustment that fully reconciles the differences between the observed and simulated galaxy outskirts, we find that artificially delaying the disruption of satellite galaxies and reducing the spatial extent of in-situ stellar populations result in improved matches between the outer profile shapes and stellar halo masses, respectively. Further insight can be achieved with higher resolution simulations that are able to better resolve satellite accretion, and with larger samples of observed galaxies.

preprint2020arXiv

A New Census of the 0.2 < z < 3.0 Universe, Part I: The Stellar Mass Function

There has been a long-standing factor-of-two tension between the observed star formation rate density and the observed stellar mass buildup after $z\sim2$. Recently we have proposed that sophisticated panchromatic SED models can resolve this tension, as these methods infer systematically higher masses and lower star formation rates than standard approaches. In a series of papers we now extend this analysis and present a complete, self-consistent census of galaxy formation over $0.2 < z < 3$ inferred with the \texttt{Prospector} galaxy SED-fitting code. In this work, Paper I, we present the evolution of the galaxy stellar mass function using new mass measurements of $\sim$10$^5$ galaxies in the 3D-HST and COSMOS-2015 surveys. We employ a new methodology to infer the mass function from the observed stellar masses: instead of fitting independent mass functions in a series of fixed redshift intervals, we construct a continuity model that directly fits for the redshift evolution of the mass function. This approach ensures a smooth picture of galaxy assembly and makes use of the full, non-Gaussian uncertainty contours in our stellar mass inferences. The resulting mass function has higher number densities at a fixed stellar mass than almost any other measurement in the literature, largely owing to the older stellar ages inferred by \texttt{Prospector}. The stellar mass density is $\sim$50% higher than previous measurements, with the offset peaking at $z\sim1$. The next two papers in this series will present the new measurements of star-forming main sequence and the cosmic star formation rate density, respectively.

preprint2020arXiv

A Tip of the Red Giant Branch Distance to the Dark Matter Deficient Galaxy NGC1052-DF4 from Deep Hubble Space Telescope Data

Previous studies have shown that the large, diffuse galaxies NGC1052-DF2 and NGC1052-DF4 both have populations of unusually luminous globular clusters as well as a very low dark matter content. Here we present newly-obtained deep Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) imaging of one of these galaxies, NGC1052-DF4. We use these data to measure the distance of the galaxy from the location of the tip of the red giant branch (TRGB). We find a rapid increase in the number of detected stars fainter than $m_{F814W} \sim 27.3$, which we identify as the onset of the red giant branch. Using a forward modeling approach that takes the photometric uncertainties into account, we find a TRGB magnitude of $m_{F814W,\rm TRGB}=27.47 \pm 0.16$. The inferred distance, including the uncertainty in the absolute calibration, is $D_{\rm TRGB}=20.0 \pm 1.6$ Mpc. The TRGB distance of NGC1052-DF4 is consistent with the previously-determined surface brightness fluctuation distance of $D_{\rm SBF}=18.7\pm 1.7$ Mpc to NGC1052-DF2 and is consistent with the distance of the bright elliptical galaxy NGC1052. We conclude that the unusual properties of these galaxies cannot be explained by distance errors.

preprint2020arXiv

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

Minor mergers have been proposed as the driving mechanism for the size growth of quiescent galaxies with decreasing redshift. The process whereby large star-forming galaxies quench and join the quiescent population at the large size end has also been suggested as an explanation for this size growth. Given the clear association of quenching with clusters, we explore this mechanism by studying the structural properties of 23 spectroscopically identified recently quenched (or "poststarburst" (PSB)) cluster galaxies at $z\sim1$. Despite clear PSB spectral signatures implying rapid and violent quenching, 87\% of these galaxies have symmetric, undisturbed morphologies in the stellar continuum. Remarkably, they follow a mass-size relation lying midway between the star-forming and quiescent field relations, with sizes $0.1$ dex smaller than $z\sim1$ star-forming galaxies at log$(M_{*}/M_{\odot})=10.5$. This implies a rapid change in the light profile without directly effecting the stellar distribution, suggesting changes in the mass-to-light ratio gradients across the galaxy are responsible. We develop fading toy models to explore how star-forming galaxies move across the mass-size plane as their stellar populations fade to match those of the PSBs. "Outside-in" fading has the potential to reproduce the contraction in size and increase in bulge-dominance observed between star-forming and PSB cluster galaxies. Since cluster PSBs lie on the large size end of the quiescent mass-size relation, and our previous work shows cluster galaxies are smaller than field galaxies, the sizes of quiescent galaxies must grow both from the quenching of star-forming galaxies and dry minor mergers.

preprint2020arXiv

Measuring Star Formation Histories, Distances, and Metallicities with Pixel Color-Magnitude Diagrams II: Applications to Nearby Elliptical Galaxies

We present spatially-resolved measurements of star formation histories (SFHs), metallicities, and distances in three nearby elliptical galaxies and the bulge of M31 derived using the pixel color-magnitude diagram (pCMD) technique. We compute pCMDs from archival $\textit{HST}$ photometry of M87, M49, NGC 3377 and M31, and fit the data using the new code $\texttt{PCMDPy}$. We measure distances to each system that are accurate to $\sim 10\%$. The recovered non-parametric SFHs place reasonable ($\pm 1$ dex) constraints on the recent (< 2 Gyr) star formation in M31 and NGC 3377, both of which show evidence of inside-out growth. The SFHs in M87 and M49 are constrained only at the oldest ages. The pCMD technique is a promising new avenue for studying the evolutionary history of the nearby universe, and is highly complementary to existing stellar population modeling techniques.

preprint2020arXiv

Multi-resolution filtering: an empirical method for isolating faint, extended emission in Dragonfly data and other low resolution images

We describe an empirical, self-contained method to isolate faint, large-scale emission in imaging data of low spatial resolution. Multi-resolution filtering (MRF) uses independent data of superior spatial resolution and point source depth to create a model for all compact and high surface brightness objects in the field. This model is convolved with an appropriate kernel and subtracted from the low resolution image. The halos of bright stars are removed in a separate step and artifacts are masked. The resulting image only contains extended emission fainter than a pre-defined surface brightness limit. The method was developed for the Dragonfly Telephoto Array, which produces images that have excellent low surface brightness sensitivity but poor spatial resolution. We demonstrate the MRF technique using Dragonfly images of a satellite of the spiral galaxy M101, the tidal debris surrounding M51, two ultra-diffuse galaxies in the Coma cluster, and the galaxy NGC5907. As part of the analysis we present a newly-identified very faint galaxy in the filtered Dragonfly image of the M101 field. We also discuss variations of the technique for cases when no low resolution data are available (self-MRF and cross-MRF). The method is implemented in mrf, an open-source MIT licensed Python package.

preprint2020arXiv

The Dragonfly Edge-on Galaxies Survey: Shaping the outer disc of NGC 4565 via accretion

We present extremely deep imaging of the well-known edge-on galaxy NGC 4565 (the "Needle Galaxy"), observed as part of the Dragonfly Edge-on Galaxies Survey. Our imaging traces starlight over scales of ~100 kpc to surface brightnesses of ~32 mag/arcsec^2 in azimuthally averaged radial profiles. In narrow slice profiles along the major axis (with bin sizes ranging from 1.7 x 0.5 kpc to 1.7 x 7.8 kpc), we achieve a depth of ~29-30 mag/arcsec^2. We confirm the previously observed asymmetric disc truncation. In addition, the sharp northwest truncation turns over to a shallower component that coincides with a fan-like feature seen to wrap around the northwest disc limb. We propose that the fan may be a tidal ribbon, and qualitatively replicate the fan with simple simulations. Alternative explanations of the fan and the disc's asymmetry are considered as well. We conclude that we are most likely seeing evidence for accretion-based outer disk growth in NGC 4565.

preprint2020arXiv

The Dragonfly Wide Field Survey. I. Telescope, Survey Design and Data Characterization

We present a description of the Dragonfly Wide Field Survey (DWFS), a deep photometric survey of a wide area of sky. The DWFS covers 330 $\mathrm{deg}^2$ in the equatorial GAMA fields and the Stripe 82 fields in the SDSS $g$ and $r$ bands. It is carried out with the 48-lens Dragonfly Telephoto Array, a telescope that is optimized for the detection of low surface brightness emission. The main goal of the survey is to study the dwarf galaxy population beyond the Local Group. In this paper, we describe the survey design and show early results. We reach $1σ$ depths of $μ_g\approx 31$ mag arcsec$^{-2}$ on arcminute scales and show that Milky Way satellites such as Sextans, Bootes, and Ursa Major should be detectable out to $D\gtrsim 10$ Mpc. We also provide an overview of the elements and operation of the 48-lens Dragonfly telescope and a detailed description of its data reduction pipeline. The pipeline is fully automated, with individual frames subjected to a rigorous series of quality tests. The sky subtraction is performed in two stages, ensuring that emission features with spatial scales up to $\sim 0.^{\circ}9 \times 0.^{\circ}6$ are preserved. The DWFS provides unparalleled sensitivity to low surface brightness features on arcminute scales.

preprint2020arXiv

Timing the Early Assembly of the Milky Way with the H3 Survey

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$α$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$α$ and high-$α$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$α$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$α$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$α$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$α$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

preprint2019arXiv

Anomalously Narrow Linewidths of Compact Massive Star-Forming Galaxies at z~2.3: A Possible Inclination Bias in the Size-Mass Plane

Compact, massive star forming galaxies at $z\sim2.5$ are thought to be building the central regions of giant elliptical galaxies today. However, a significant fraction of these objects were previously shown to have much smaller H$α$ line widths than expected. A possible interpretation is that H$α$ emission from their central regions, where the highest velocities are expected, is typically obscured by dust. Here we present ALMA observations of the CO(3-2) emission line of three compact, massive galaxies with H$α$ line widths of FWHM(H$α$)$\sim$125-260 km s$^{-1}$ to test this hypothesis. Surprisingly, in all three galaxies, the CO line width is similar to the H$α$ line width: we find FWHM(CO)$\sim$165 km s$^{-1}$ for all three galaxies whereas FWHM(CO)$\sim$450-700 km s$^{-1}$ was expected from a simple virial estimator. These results show that the narrow H$α$ linewidths of many compact massive star-forming galaxies are not due to preferential obscuration of the highest velocity gas. An alternative explanation for the narrow line widths is that the galaxies are disks that are viewed nearly face-on. We suggest that there may be an inclination bias in the size-mass plane, such that the apparent rest-frame optical sizes of face-on galaxies are smaller than those of edge-on galaxies. Although not conclusive, this hypothesis is supported by an observed anti-correlation between the size and axis ratio of massive galaxies.

preprint2019arXiv

Spatially Resolved Stellar Kinematics of the Ultra-Diffuse Galaxy Dragonfly 44: II. Constraints on Fuzzy Dark Matter

Given the absence of directly detected dark matter (DM) as weakly interacting massive particles, there is strong interest in the possibility that DM is an ultra-light scalar field, here denoted as fuzzy DM. Ultra-diffuse galaxies, with the sizes of giant galaxies and the luminosities of dwarf galaxies, have a wide range of DM halo masses, thus providing new opportunities for exploring the connections between galaxies and their DM halos. Following up on new integral field unit spectroscopic observations and dynamics modeling of the DM-dominated ultra-diffuse galaxy Dragonfly 44 in the outskirts of the Coma Cluster, we present models of fuzzy DM constrained by the stellar dynamics of this galaxy. We infer a scalar field mass of $\sim 3 \times 10^{-22}$ eV, consistent with other constraints from galaxy dynamics but in tension with constraints from Ly$α$ forest power spectrum modeling. While we are unable to statistically distinguish between fuzzy DM and normal cold DM models, we find that the inferred properties of the fuzzy DM halo satisfy a number of predictions for halos in a fuzzy DM cosmology. In particular, we find good agreement with the predicted core size--halo mass relation and the predicted transition radius between the quantum pressure-dominated inner region and the outer halo region.

preprint2018arXiv

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

We present a targeted follow-up Hubble Space Telescope WFC3 F160W imaging study of very massive galaxies $(\log(M_{\rm{star}}/M_{\odot})> 11.2)$ selected from a combination of ground-based near-infrared galaxy surveys (UltraVISTA, NMBS-II, UKIDSS UDS) at $1.5<z<3$. We find that these galaxies are diverse in their structures, with $\sim1/3$ of the targets being composed of close pairs, and span a wide range in sizes. At $1.5<z<2.5$, the sizes of both star-forming and quiescent galaxies are consistent with the extrapolation of the stellar mass-size relations determined at lower stellar masses. At $2.5<z<3.0$, however, we find evidence that quiescent galaxies are systematically larger than expected based on the extrapolation of the relation derived using lower stellar mass galaxies. We used the observed light profiles of the blended systems to decompose their stellar masses and investigate the effect of the close pairs on the measured number densities of very massive galaxies in the early universe. We estimate correction factors to account for close-pair blends and apply them to the observed stellar mass functions measured using ground-based surveys. Given the large uncertainties associated with this extreme population of galaxies, there is currently little tension between the (blending-corrected) number density estimates and predictions from theoretical models. Although we currently lack the statistics to robustly correct for close-pair blends, we show that this is a systematic effect which can reduce the observed number density of very massive galaxies by up to a factor of $\sim1.5$, and should be accounted for in future studies of stellar mass functions.

preprint2016arXiv

A High Stellar Velocity Dispersion and ~100 Globular Clusters for the Ultra Diffuse Galaxy Dragonfly 44

Recently a population of large, very low surface brightness, spheroidal galaxies was identified in the Coma cluster. The apparent survival of these Ultra Diffuse Galaxies (UDGs) in a rich cluster suggests that they have very high masses. Here we present the stellar kinematics of Dragonfly 44, one of the largest Coma UDGs, using a 33.5 hr integration with DEIMOS on the Keck II telescope. We find a velocity dispersion of 47 km/s, which implies a dynamical mass of M_dyn=0.7x10^10 M_sun within its deprojected half-light radius of r_1/2=4.6 kpc. The mass-to-light ratio is M/L=48 M_sun/L_sun, and the dark matter fraction is 98 percent within the half-light radius. The high mass of Dragonfly 44 is accompanied by a large globular cluster population. From deep Gemini imaging taken in 0.4" seeing we infer that Dragonfly 44 has 94 globular clusters, similar to the counts for other galaxies in this mass range. Our results add to other recent evidence that many UDGs are "failed" galaxies, with the sizes, dark matter content, and globular cluster systems of much more luminous objects. We estimate the total dark halo mass of Dragonfly 44 by comparing the amount of dark matter within r=4.6 kpc to enclosed mass profiles of NFW halos. The enclosed mass suggests a total mass of ~10^12 M_sun, similar to the mass of the Milky Way. The existence of nearly-dark objects with this mass is unexpected, as galaxy formation is thought to be maximally-efficient in this regime.

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

Pixel Color Magnitude Diagrams for Semi-Resolved Stellar Populations: The Star Formation History of Regions within the Disk and Bulge of M31

The analysis of stellar populations has, by and large, been developed for two limiting cases: spatially-resolved stellar populations in the color-magnitude diagram, and integrated light observations of distant systems. In between these two extremes lies the semi-resolved regime, which encompasses a rich and relatively unexplored realm of observational phenomena. Here we develop the concept of pixel color magnitude diagrams (pCMDs) as a powerful technique for analyzing stellar populations in the semi-resolved regime. pCMDs show the distribution of imaging data in the plane of pixel luminosity vs. pixel color. A key feature of pCMDs is that they are sensitive to all stars, including both the evolved giants and the unevolved main sequence stars. An important variable in this regime is the mean number of stars per pixel, $N_{\rm pix}$. Simulated pCMDs demonstrate a strong sensitivity to the star formation history (SFH) and allow one to break degeneracies between age, metallicity and dust based on two filter data for values of $N_{\rm pix}$ up to at least $10^4$. We extract pCMDs from Hubble Space Telescope (HST) optical imaging of M31 and derive non-parametric SFHs from $10^6$ yr to $10^{10}$ yr for both the crowded disk and bulge regions (where $N_{\rm pix}\approx30-10^3$). From analyzing a small region of the disk we find a non-parametric SFH that is smooth and consistent with an exponential decay timescale of 4 Gyr. The bulge SFH is also smooth and consistent with a 2 Gyr decay timescale. pCMDs will likely play an important role in maximizing the science returns from next generation ground and space-based facilities.

preprint2016arXiv

The Dragonfly Nearby Galaxies Survey. I. Substantial variation in the diffuse stellar halos around spiral galaxies

Galaxies are thought to grow through accretion; as less massive galaxies are disrupted and merge over time, their debris results in diffuse, clumpy stellar halos enveloping the central galaxy. Here we present a study of the variation in the stellar halos of galaxies, using data from the Dragonfly Nearby Galaxies Survey (DNGS). The survey consists of wide field, deep ($μ_{g} > 31$ mag arcsec$^{-2}$) optical imaging of nearby galaxies using the Dragonfly Telephoto Array. Our sample includes eight spiral galaxies with stellar masses similar to that of the Milky Way, inclinations of $16-90$ degrees and distances between $7-18$ Mpc. We construct stellar mass surface density profiles from the observed $g$-band surface brightness in combination with the $g-r$ color as a function of radius, and compute the halo fractions from the excess stellar mass (relative to a disk$+$bulge fit) beyond $5$ half-mass radii. We find a mean halo fraction of $0.009 \pm 0.005$ and a large RMS scatter of $1.01^{+0.9}_{-0.26}$ dex. The peak-to-peak scatter is a factor of $>100$ -- while some galaxies feature strongly structured halos resembling that of M31, three of the eight have halos that are completely undetected in our data. We conclude that spiral galaxies as a class exhibit a rich variety in stellar halo properties, implying that their assembly histories have been highly non-uniform. We find no convincing evidence for an environmental or stellar mass dependence of the halo fraction in the sample.

preprint2016arXiv

The Dragonfly Nearby Galaxies Survey. II. Ultra diffuse galaxies near the elliptical galaxy NGC 5485

We present the unexpected discovery of four ultra diffuse galaxies (UDGs) in a group environment. We recently identified seven extremely low surface brightness galaxies in the vicinity of the spiral galaxy M101, using data from the Dragonfly Telephoto Array. The galaxies have effective radii of $10"-38"$ and central surface brightnesses of $25.6-27.7$ mag arcsec$^{-2}$ in g-band. We subsequently obtained follow-up observations with $HST$ to constrain the distances to these galaxies. Four remain persistently unresolved even with the spatial resolution of $HST$/ACS, which implies distances of $D > 17.5$ Mpc. We show that the galaxies are most likely associated with a background group at $\sim 27$ Mpc containing the massive ellipticals NGC 5485 and NGC 5473. At this distance, the galaxies have sizes of $2.6-4.9$ kpc, and are classified as UDGs, similar to the populations that have been revealed in clusters such as Coma, Virgo and Fornax, yet even more diffuse. The discovery of four UDGs in a galaxy group demonstrates that the UDG phenomenon is not exclusive to cluster environments. Furthermore, their morphologies seem less regular than those of the cluster populations, which may suggest a different formation mechanism or be indicative of a threshold in surface density below which UDGs are unable to maintain stability.

preprint2016arXiv

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

Using a sample of 98 galaxy clusters recently imaged in the near infra-red with the ESO NTT, WIYN and WHT telescopes, supplemented with 33 clusters from the ESO archive, we measure how the stellar mass of the most massive galaxies in the universe, namely Brightest Cluster Galaxies (BCG), increases with time. Most of the BCGs in this new sample lie in the redshift range $0.2<z<0.6$, which has been noted in recent works to mark an epoch over which the growth in the stellar mass of BCGs stalls. From this sample of 132 clusters, we create a subsample of 102 systems that includes only those clusters that have estimates of the cluster mass. We combine the BCGs in this subsample with BCGs from the literature, and find that the growth in stellar mass of BCGs from 10 billion years ago to the present epoch is broadly consistent with recent semi-analytic and semi-empirical models. As in other recent studies, tentative evidence indicates that the stellar mass growth rate of BCGs may be slowing in the past 3.5 billion years. Further work in collecting larger samples, and in better comparing observations with theory using mock images is required if a more detailed comparison between the models and the data is to be made.

preprint2016arXiv

The evolution of metallicity and metallicity gradients from z=2.7-0.6 with KMOS3D

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

preprint2016arXiv

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

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

preprint2015arXiv

Episodic molecular outflow in the very young protostellar cluster Serpens South

The loss of mass from protostars, in the form of a jet or outflow, is a necessary counterpart to protostellar mass accretion. Outflow ejection events probably vary in their velocity and/or in the rate of mass loss. Such `episodic' ejection events have been observed during the Class 0 protostellar phase (the early accretion stage), and continue during the subsequent class I phase that marks the first one million years of star formation. Previously observed episodic-ejection sources were relatively isolated; however, the most common sites of star formation are clusters. Outflows link protostars with their environment and provide a viable source of turbulence that is necessary for regulating star formation in clusters, but it is not known how an accretion-driven jet or outflow in a clustered environment manifests itself in its earliest stage. This early stage is important in establishing the initial conditions for momentum and energy transfer to the environment as the protostar and cluster evolve. Here we report that an outflow from a very young class 0 protostar, at the hub of the very active and filamentary Serpens South protostellar cluster, shows unambiguous episodic events. The $^{12}$CO (J=2-1) emission from the protostar reveals 22 distinct features of outflow ejecta, the most recent having the highest velocity. The outflow forms bipolar lobes --- one of the first detectable signs of star formation --- which originate from the peak of 1-mm continuum emission. Emission from the surrounding C$^{18}$O envelope shows kinematics consistent with rotation and an infall of material onto the protostar. The data suggest that episodic accretion-driven outflow begins in the earliest phase of protostellar evolution, and that the outflow remains intact in a very clustered environment, probably providing efficient momentum transfer for driving turbulence.

preprint2015arXiv

Forty-Seven Milky Way-Sized, Extremely Diffuse Galaxies in the Coma Cluster

We report the discovery of 47 low surface brightness objects in deep images of a 3 x 3 degree field centered on the Coma cluster, obtained with the Dragonfly Telephoto Array. The objects have central surface brightness mu(g,0) ranging from 24 - 26 mag/arcsec^2 and effective radii r_e = 3"-10", as measured from archival Canada France Hawaii Telescope images. From their spatial distribution we infer that most or all of the objects are galaxies in the Coma cluster. This relatively large distance is surprising as it implies that the galaxies are very large: with r_e = 1.5 - 4.6 kpc their sizes are similar to those of L* galaxies even though their median stellar mass is only ~6 x 10^7 Solar masses. The galaxies are relatively red and round, with <g-i> = 0.8 and <b/a> = 0.74. One of the 47 galaxies is fortuitously covered by a deep Hubble Space Telescope ACS observation. The ACS imaging shows a large spheroidal object with a central surface brightness mu(g,0) = 25.8 mag/arcsec^2, a Sersic index n=0.6, and an effective radius of 7", corresponding to 3.4 kpc at the distance of Coma. The galaxy is not resolved into stars, consistent with expectations for a Coma cluster object. We speculate that UDGs may have lost their gas supply at early times, possibly resulting in very high dark matter fractions.

preprint2015arXiv

Preventing Star Formation in Early-Type Galaxies with Late-Time Stellar Heating

We revisit previous suggestions that the heating provided by the winds of dying low-mass stars plays an important role in preventing star formation in quiescent galaxies. At the end of their asymptotic giant branch phase, intermediate and low-mass stars eject their envelopes rapidly in a super-wind phase, usually giving rise to planetary nebulae. In spheroidal galaxies with high stellar velocity dispersions, the interaction of these ejected envelopes with the ambient diffuse gas can lead to significant, isotropic and steady-state heating that scales as $\dot{M}_\astσ_\ast^2$. We show that cooling of the central regions of the hot diffuse halo gas can be delayed for a Hubble time for halos more massive than $\sim10^{12.5}M_{\odot}$ at $0<z<2$, although stellar heating alone is unlikely to forestall cooling in the most massive clusters at $z=0$. This mechanism provides a natural explanation for the strong trend of galaxy quiescence with stellar surface density and velocity dispersion. In addition, since the ejected material will thermalize to $kT\simσ_\ast^2$, this mechanism provides an explanation for the observed similarity between the central temperature of the hot diffuse gas and $σ_\ast^2$, a result which is not trivial in light of the short inferred cooling times of the hot gas. The main uncertainty in this analysis is the ultimate fate of the stellar ejecta. Preventing accumulation of the ejecta in the central regions may require energy input from another source, such as Type Ia supernovae. Detailed simulations of the interaction of the stellar wind with the ambient gas are required to better quantify the net effect of AGB heating.

preprint2015arXiv

Satellite Quenching and Galactic Conformity at 0.3 < z < 2.5

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3σ$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

preprint2015arXiv

Spectroscopic Confirmation of the Existence of Large, Diffuse Galaxies in the Coma Cluster

We recently identified a population of low surface brightness objects in the field of the z=0.023 Coma cluster, using the Dragonfly Telephoto Array. Here we present Keck spectroscopy of one of the largest of these "ultra-diffuse galaxies" (UDGs), confirming that it is a member of the cluster. The galaxy has prominent absorption features, including the Ca II H+K lines and the G-band, and no detected emission lines. Its radial velocity of cz=6280 +- 120 km/s is within the 1 sigma velocity dispersion of the Coma cluster. The galaxy has an effective radius of 4.3 +- 0.3 kpc and a Sersic index of 0.89 +- 0.06, as measured from Keck imaging. We find no indications of tidal tails or other distortions, at least out to a radius of ~2 r_e. We show that UDGs are located in a previously sparsely populated region of the size - magnitude plane of quiescent stellar systems, as they are ~6 magnitudes fainter than normal early-type galaxies of the same size. It appears that the luminosity distribution of large quiescent galaxies is not continuous, although this could largely be due to selection effects. Dynamical measurements are needed to determine whether the dark matter halos of UDGs are similar to those of galaxies with the same luminosity or to those of galaxies with the same size.

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 sizes of massive quiescent and star forming galaxies at z~4 with ZFOURGE and CANDELS

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

preprint2014arXiv

A Consistent Study of Metallicity Evolution at 0.8 < z < 2.6

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

preprint2014arXiv

A Massive Galaxy in its Core Formation Phase Three Billion Years After the Big Bang

Most massive galaxies are thought to have formed their dense stellar cores at early cosmic epochs. However, cores in their formation phase have not yet been observed. Previous studies have found galaxies with high gas velocity dispersions or small apparent sizes but so far no objects have been identified with both the stellar structure and the gas dynamics of a forming core. Here we present a candidate core in formation 11 billion years ago, at z=2.3. GOODS-N-774 has a stellar mass of 1.0x10^11 Msun, a half-light radius of 1.0 kpc, and a star formation rate of 90[+45-20]Msun/yr. The star forming gas has a velocity dispersion 317+-30 km/s, amongst the highest ever measured. It is similar to the stellar velocity dispersions of the putative descendants of GOODS-N-774, compact quiescent galaxies at z~2 and giant elliptical galaxies in the nearby Universe. Galaxies such as GOODS-N-774 appear to be rare; however, from the star formation rate and size of the galaxy we infer that many star forming cores may be heavily obscured, and could be missed in optical and near-infrared surveys.

preprint2014arXiv

A substantial population of massive quiescent galaxies at z~4 from ZFOURGE

We report the likely identification of a substantial population of massive M~10^11M_Sun galaxies at z~4 with suppressed star formation rates (SFRs), selected on rest-frame optical to near-IR colors from the FourStar Galaxy Evolution Survey. The observed spectral energy distributions show pronounced breaks, sampled by a set of near-IR medium-bandwidth filters, resulting in tightly constrained photometric redshifts. Fitting stellar population models suggests large Balmer/4000AA breaks, relatively old stellar populations, large stellar masses and low SFRs, with a median specific SFR of 2.9+/-1.8 x 10^-11/yr. Ultradeep Herschel/PACS 100micron, 160micron and Spitzer/MIPS 24micron data reveal no dust-obscured SFR activity for 15/19 (79%) galaxies. Two far-IR detected galaxies are obscured QSOs. Stacking the far-IR undetected galaxies yields no detection, consistent with the SED fit, indicating independently that the average specific SFR is at least 10x smaller than of typical star-forming galaxies at z~4. Assuming all far-IR undetected galaxies are indeed quiescent, the volume density is 1.8+/-0.7 x 10^-5Mpc^-3 to a limit of log10M/M_Sun>10.6, which is 10x and 80x lower than at z = 2 and z = 0.1. They comprise a remarkably high fraction (~35%) of z~4 massive galaxies, suggesting that suppression of star formation was efficient even at very high redshift. Given the average stellar age of 0.8Gyr and stellar mass of 0.8x10^11M_Sun, the galaxies likely started forming stars before z =5, with SFRs well in excess of 100M_Sun/yr, far exceeding that of similarly abundant UV-bright galaxies at z>4. This suggests that most of the star-formation in the progenitors of quiescent z~4 galaxies was obscured by dust.

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

Exploring the z=3-4 massive galaxy population with ZFOURGE: the prevalence of dusty and quiescent galaxies

Our understanding of the redshift $z>3$ galaxy population relies largely on samples selected using the popular "dropout" technique, typically consisting of UV-bright galaxies with blue colors and prominent Lyman breaks. As it is currently unknown if these galaxies are representative of the massive galaxy population, we here use the FourStar Galaxy Evolution (ZFOURGE) Survey to create a stellar mass-limited sample at $z=3-4$. Uniquely, ZFOURGE uses deep near-infrared medium-bandwidth filters to derive accurate photometric redshifts and stellar population properties. The mass-complete sample consists of 57 galaxies with log M $>10.6$, reaching below $M^{\star}$ at $z=3-4$. On average, the massive $z=3-4$ galaxies are extremely faint in the observed optical with median $R_{tot}^{AB}=27.48\pm0.41$ (restframe $M_{1700}=-18.05\pm0.37$). They lie far below the UV luminosity-stellar mass relation for Lyman break galaxies and are about $\sim100\times$ fainter at the same mass. The massive galaxies are red ($R-Ks_{AB}=3.9\pm0.2$; restframe UV-slope $β=-0.2\pm0.3$) likely from dust or old stellar ages. We classify the galaxy SEDs by their restframe $U-V$ and $V-J$ colors and find a diverse population: $46^{+6+10}_{-6-17}$% of the massive galaxies are quiescent, $54^{+8+17}_{-8-10}$% are dusty star-forming galaxies, and only $14^{+3+10}_{-3-4}$% resemble luminous blue star forming Lyman break galaxies. This study clearly demonstrates an inherent diversity among massive galaxies at higher redshift than previously known. Furthermore,we uncover a reservoir of dusty star-forming galaxies with $4\times$ lower specific star-formation rates compared to submillimeter-selected starbursts at $z>3$. With $5\times$ higher numbers, the dusty galaxies may represent a more typical mode of star formation compared to submillimeter-bright starbursts.

preprint2014arXiv

First results from the Dragonfly Telephoto Array: the apparent lack of a stellar halo in the massive spiral galaxy M101

We use a new telescope concept, the Dragonfly Telephoto Array, to study the low surface brightness outskirts of the spiral galaxy M101. The radial surface brightness profile is measured down to mu_g ~ 32 mag/arcsec^2, a depth that approaches the sensitivity of star count studies in the Local Group. We convert surface brightness to surface mass density using the radial g-r color profile. The mass density profile shows no significant upturn at large radius and is well-approximated by a simple bulge + disk model out to R = 70 kpc, corresponding to 18 disk scale lengths. Fitting a bulge + disk + halo model we find that the best-fitting halo mass M_halo ~ 1.7 x 10^8 M_sun. The total stellar mass of M101 is M_tot,* ~ 5.3 x 10^10 Msun, and we infer that the halo mass fraction f_halo = M_halo / M_tot,* ~ 0.003. This mass fraction is lower than that of the Milky Way (f_halo ~ 0.02) and M31 (f_halo ~ 0.04). All three galaxies fall below the f_halo - M_tot,* relation predicted by recent cosmological simulations that trace the light of disrupted satellites, with M101's halo mass a factor of ~10 below the median expectation. However, the predicted scatter in this relation is large, and more galaxies are needed to better quantify this possible tension with galaxy formation models. Dragonfly is well suited for this project: as integrated-light surface brightness is independent of distance, large numbers of galaxies can be studied in a uniform way.

preprint2014arXiv

Fluctuation Spectroscopy: A New Probe of Old Stellar Populations

We introduce a new method to determine the relative contributions of different types of stars to the integrated light of nearby early-type galaxies. As is well known, the surface brightness of these galaxies shows pixel-to-pixel fluctuations due to Poisson variations in the number of giant stars. Differential spectroscopy of pixels as a function of fluctuation strength ("fluctuation spectroscopy") effectively measures the spectral variation of stars as a function of their luminosity, information that is otherwise difficult to obtain for individual stars outside of the Local Group. We apply this technique to the elliptical galaxy NGC 4472, using HST/ACS imaging in six narrow-band ramp filters tuned to spectral features in the range 0.8-1.0 micron. Pixels with +- 5% broad-band variations show differential color variations of 0.1% - 1.0% in the narrow-band filters. These variations are primarily due to the systematic increase in TiO absorption strength with increasing luminosity on the upper giant branch. The data are very well reproduced by the same Conroy & van Dokkum (2012) stellar population synthesis model that is the best fit to the integrated light, with residuals in the range 0.03% - 0.09%. Models with ages or metallicities that are significantly different from the integrated-light values do not yield good fits. We can also rule out several modifications to the underlying model, including the presence of a significant (>3% of the light) population of late M giants. The current observations constitute a powerful test of the expected luminosities and temperatures of metal-rich giants in massive early-type galaxies. Studies of pixels with much larger (negative) fluctuations will provide unique information on main sequence stars and the stellar initial mass function.

preprint2014arXiv

Galaxy Stellar Mass Functions from ZFOURGE/CANDELS: An Excess of Low-Mass Galaxies Since z=2 and the Rapid Buildup of Quiescent Galaxies

Using observations from the FourStar Galaxy Evolution Survey (ZFOURGE), we obtain the deepest measurements to date of the galaxy stellar mass function at 0.5 < z < 2.5. ZFOURGE provides well-constrained photometric redshifts made possible through deep medium-bandwidth imaging at 1-2um . We combine this with HST imaging from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS), allowing for the efficient selection of both blue and red galaxies down to stellar masses ~10^9.5 Msol at z ~ 2.5. The total surveyed area is 316 arcmin^2 distributed over three independent fields. We supplement these data with the wider and shallower NEWFIRM Medium-Band Survey (NMBS) to provide stronger constraints at high masses. Several studies at z<=1 have revealed a steepening of the slope at the low-mass end of the stellar mass function (SMF), leading to an upturn at masses <10^10 Msol that is not well-described by a standard single-Schechter function. We find evidence that this feature extends to at least z ~ 2, and that it can be found in both the star-forming and quiescent populations individually. The characteristic mass (M*) and slope at the lowest masses (alpha) of a double-Schechter function fit to the SMF stay roughly constant at Log(M/Msol) ~ 10.65 and ~-1.5 respectively. The SMF of star-forming galaxies has evolved primarily in normalization, while the change in shape is relatively minor. This is not the case for quiescent galaxies: the depth of our imaging allows us to show for the first time significantly more evolution at Log(M/Msol) < 10.5 than at higher masses. We find that the total mass density (down to 10^9 Msol) in star-forming galaxies has increased by a factor of ~2.2 since z ~ 2.5, whereas in quiescent galaxies it has increased by a factor of ~12 .

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

One Plane for All: Massive Star-Forming and Quiescent Galaxies Lie on the Same Mass Fundamental Plane at z~0 and z~0.7

Scaling relations between galaxy structures and dynamics have been studied extensively for early and late-type galaxies, both in the local universe and at high redshifts. The abundant differences between the properties of disky and elliptical, or star-forming and quiescent, galaxies seem to be characteristic of the local Universe; such clear distinctions begin to disintegrate as observations of massive galaxies probe higher redshifts. In this Paper, we investigate the existence the mass fundamental plane of all massive galaxies ($σ\gtrsim$ 100 km/s). This work includes local galaxies (0.05<z<0.07) from the SDSS, in addition to 31 star-forming and 72 quiescent massive galaxies at intermediate redshift (z~0.7) with absorption line kinematics from deep Keck-DEIMOS spectra and structural parameters from HST imaging. In two parameter scaling relations, star-forming and quiescent galaxies differ structurally and dynamically. However, we show that massive star-forming and quiescent galaxies lie on nearly the same mass fundamental plane, or the relationship between stellar mass surface density, stellar velocity dispersion, and effective radius. The scatter in this relation (measured about $\logσ$) is low: 0.072 dex (0.055 dex intrinsic) at z~0 and 0.10 dex (0.08 dex intrinsic) at z~0.7. This three dimensional surface is not unique: virial relations, with or without a dependence on luminosity profile shapes, can connect galaxy structures and stellar dynamics with similar scatter. This result builds on the recent finding that mass fundamental plane has been stable for early-type galaxies since z~2 (Bezanson et al. 2013). As we now find this also holds for star-forming galaxies to z~0.7, this implies that these scaling relations of galaxies will be minimally susceptible to progenitor biases due to the evolving stellar populations, structures, and dynamics of galaxies through cosmic time.

preprint2014arXiv

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

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

preprint2014arXiv

The discovery of seven extremely low surface brightness galaxies in the field of the nearby spiral galaxy M101

Dwarf satellite galaxies are a key probe of dark matter and of galaxy formation on small scales and of the dark matter halo masses of their central galaxies. They have very low surface brightness, which makes it difficult to identify and study them outside of the Local Group. We used a low surface brightness-optimized telescope, the Dragonfly Telephoto Array, to search for dwarf galaxies in the field of the massive spiral galaxy M101. We identify seven large, low surface brightness objects in this field, with effective radii of \(10 - 30\) arcseconds and central surface brightnesses of \(μ_{g} \sim 25.5 - 27.5\) mag arcsec\(^{-2}\). Given their large apparent sizes and low surface brightnesses, these objects would likely be missed by standard galaxy searches in deep fields. Assuming the galaxies are dwarf satellites of M101, their absolute magnitudes are in the range \(-11.6 \lesssim M_{V} \lesssim -9.3\) and their effective radii are \(350\) pc \(-\) \(1.3\) kpc. Their radial surface brightness profiles are well fit by Sersic profiles with a very low Sersic index (\(n \sim 0.3 - 0.7\)). The properties of the sample are similar to those of well-studied dwarf galaxies in the Local Group, such as Sextans I and Phoenix. Distance measurements are required to determine whether these galaxies are in fact associated with M101 or are in its foreground or background.

preprint2014arXiv

The distribution of satellites around massive galaxies at 1<z<3 in ZFOURGE/CANDELS: dependence on star formation activity

We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$σ$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

preprint2014arXiv

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

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

preprint2014arXiv

Ultra Low Surface Brightness Imaging with the Dragonfly Telephoto Array

We describe the Dragonfly Telephoto Array, a robotic imaging system optimized for the detection of extended ultra low surface brightness structures. The array consists of eight Canon 400mm f/2.8 telephoto lenses coupled to eight science-grade commercial CCD cameras. The lenses are mounted on a common framework and are co-aligned to image simultaneously the same position on the sky. The system provides an imaging capability equivalent to a 0.4m aperture f/1.0 refractor with a 2.6 deg X 1.9 deg field of view. The system has no obstructions in the light path, optimized baffling, and internal optical surfaces coated with a new generation of anti-reflection coatings based on sub-wavelength nanostructures. As a result, the array's point spread function has a factor of ~10 less scattered light at large radii than well-baffled reflecting telescopes. The Dragonfly Telephoto Array is capable of imaging extended structures to surface brightness levels below 30 mag/arcsec^2 in 10h integrations (without binning or foreground star removal). This is considerably deeper than the surface brightness limit of any existing wide-field telescope. At present no systematic errors limiting the usefulness of much longer integration times has been identified. With longer integrations (50-100h), foreground star removal and modest binning the Dragonfly Telephoto Array is capable of probing structures with surface brightnesses below 32 mag/arcsec^2. Detection of structures at these surface brightness levels may hold the key to solving the "missing substructure" and "missing satellite" problems of conventional hierarchical galaxy formation models. The Dragonfly Telephoto Array is therefore executing a fully-automated multi-year imaging survey of a complete sample of nearby galaxies in order to undertake the first census of ultra-faint substructures in the nearby Universe.

preprint2013arXiv

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

We present a catalog covering 1.62 deg^2 of the COSMOS/UltraVISTA field with PSF-matched photometry in 30 photometric bands. The catalog covers the wavelength range 0.15um - 24um including the available GALEX, Subaru, CFHT, VISTA and Spitzer data. Catalog sources have been selected from the DR1 UltraVISTA Ks band imaging that reaches a depth of K_{s,tot} = 23.4 AB (90% completeness). The PSF-matched catalog is generated using position-dependent PSFs ensuring accurate colors across the entire field. Also included is a catalog of photometric redshifts (z_phot) for all galaxies computed with the EAZY code. Comparison with spectroscopy from the zCOSMOS 10k bright sample shows that up to z ~ 1.5 the z_phot are accurate to dz/(1 + z) = 0.013, with a catastrophic outlier fraction of only 1.6%. The z_phot also show good agreement with the z_phot from the NEWFIRM Medium Band Survey (NMBS) out to z ~ 3. A catalog of stellar masses and stellar population parameters for galaxies determined using the FAST spectral energy distribution fitting code is provided for all galaxies. Also included are rest-frame U-V and V-J colors, L_2800 and L_IR. The UVJ color-color diagram confirms that the galaxy bi-modality is well-established out to z ~ 2. Star-forming galaxies also obey a star forming "main sequence" out to z ~ 2.5, and this sequence evolves in a manner consistent with previous measurements. The COSMOS/UltraVISTA Ks-selected catalog covers a unique parameter space in both depth, area, and multi-wavelength coverage and promises to be a useful tool for studying the growth of the galaxy population out to z ~ 3 - 4.

preprint2013arXiv

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

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

preprint2013arXiv

Early-Type Galaxy Archeology: Ages, Abundance Ratios, and Effective Temperatures from Full-Spectrum Fitting

The stellar populations of galaxies hold vital clues to their formation histories. In this paper we present results based on modeling stacked spectra of early-type galaxies drawn from the Sloan Digital Sky Survey (SDSS) as a function of velocity dispersion, sigma, from 90 km/s to 300 km/s. The spectra are of extremely high quality, with typical S/N of 1000/A, and a wavelength coverage of 4000A-8800A. Our population synthesis model includes variation in 16 elements from C to Ba, the shift in effective temperature, Delta(Teff), of the stars with respect to a solar metallicity isochrone, amongst other parameters. In our approach we fit the full optical spectra rather than a select number of spectral indices and are able to, for the first time, measure the abundances of the elements V, Cr, Mn, Co, and Ni from the integrated light of distant galaxies. Our main results are as follows: 1) light-weighted stellar ages range from 6-12 Gyr from low to high sigma; 2) [Fe/H] varies by less than 0.1 dex across the entire sample; 3) Mg closely tracks O, and both increase from ~0.0 at low sigma to ~0.25 at high sigma; Si and Ti show a shallower rise with sigma, and Ca tracks Fe rather than O; 4) the iron peak elements V, Cr, Mn, and Ni track Fe, while Co tracks O, suggesting that Co forms primarily in massive stars; 5) C and N track O over the full sample and [C/Fe] and [N/Fe] exceed 0.2 at high sigma; and 6) the variation in Delta(Teff) with total metallicity follows theoretical predictions based on stellar evolution theory. Our derived [Mg/Fe] and [O/Fe] abundance ratios are 0.05-0.1 dex lower than most previous determinations. Under the conventional interpretation that the variation in these ratios is due to star formation timescale variations, our results suggest longer star formation timescales for massive early-type galaxies than previous studies. (ABRIDGED)

preprint2013arXiv

Strontium and Barium In Early-Type Galaxies

The detailed abundance patterns of the stars within galaxies provide a unique window into the history of star formation (SF) at early times. Two widely used `chronometers' include the alpha and iron-peak elements, which are created on short and long timescales, respectively. These two clocks have been widely used to estimate SF timescales from moderate-resolution spectra of early-type galaxies. Elements formed via s-process neutron captures (e.g., Sr and Ba) comprise a third type of chronometer, as the site of the main s-process is believed to be intermediate and low-mass asymptotic giant branch stars. The [alpha/Ba] ratio in particular should provide a powerful new constraint on the SF histories of galaxies, in part because it is insensitive to the uncertain distribution of type Ia SNe detonation times and the overall Ia rate. Here we present new measurements of the abundance of Sr and Ba in nearby early-type galaxies by applying stellar population synthesis tools to high S/N optical spectra. We find a strong anti-correlation between [Mg/Fe] and [Ba/Fe], and a strong positive correlation between [Mg/Ba] and galaxy velocity dispersion. These trends are consistent with the idea that more massive galaxies formed their stars on shorter timescales compared to less massive galaxies, and rule out several other proposed explanations for the observed super-solar [Mg/Fe] values in massive galaxies. In contrast, [Sr/Fe]~0, with no strong variation across the sample. It is difficult to interpret the Sr trends without detailed chemical evolution models owing to the multiplicity of proposed nucleosynthetic sites for Sr.

preprint2013arXiv

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

We present measurements of the stellar mass functions (SMFs) of star-forming and quiescent galaxies to z = 4 using a sample of 95 675 galaxies in the COSMOS/UltraVISTA field. Sources have been selected from the DR1 UltraVISTA K_{s}-band imaging which covers a unique combination of a wide area (1.62 deg^2), to a significant depth (K_{s,tot} = 23.4). The SMFs of the combined population are in good agreement with previous measurements and show that the stellar mass density of the universe was only 50%, 10% and 1% of its current value at z ~ 0.75, 2.0, and 3.5, respectively. The quiescent population drives most of the overall growth, with the stellar mass density of these galaxies increasing by 2.71^{+0.93}_{-0.22} dex since z = 3.5. At z > 2.5, star-forming galaxies dominate the total SMF at all stellar masses, although a nonzero population of quiescent galaxies persists to z = 4. Comparisons of the K_{s}-selected star-forming galaxy SMFs to UV-selected SMFs at 2.5 < z < 4 show reasonable agreement and suggests UV-selected samples are representative of the majority of the stellar mass density at z > 3.5. We estimate the average mass growth of individual galaxies by selecting galaxies at fixed cumulative number density. The average galaxy with Log(M_{*}/M_{sun}) = 11.5 at z = 0.3 has grown in mass by only 0.2 dex (0.3 dex) since z = 2.0(3.5), whereas those with Log(M_{*}/M_{sun}) = 10.5 have grown by > 1.0 dex since z = 2. At z < 2, the time derivatives of the mass growth are always larger for lower-mass galaxies, which demonstrates that the mass growth in galaxies since that redshift is mass-dependent and primarily bottom-up. Lastly, we examine potential sources of systematic uncertainties on the SMFs and find that those from photo-z templates, SPS modeling, and the definition of quiescent galaxies dominate the total error budget in the SMFs.

preprint2013arXiv

The 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

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

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

preprint2013arXiv

The velocity function of dark matter halos at r=20 kpc: Remarkably little evolution since z ~ 4

We investigate the evolution in the dark matter halo circular velocity function, measured at a fixed physical radius of 20 kpc (v_20), which is likely to be a good proxy for galaxy circular velocity, in the Millennium-II simulation. We find that the v_20 function evolves remarkably little since z ~ 4. We analyze the histories of the main progenitors of halos, and we find that the dark matter distribution within the central 20 kpc of massive halos has been in place since early times. This provides evidence for the inside-out growth of haloes. The constancy of the central circular velocity of halos may offer a natural explanation for the observational finding that the galaxy circular velocity is an excellent predictor of various galaxy properties. Our results also indicate that we can expect a significant number of galaxies with high circular velocities already at z=4 (more than one per 10^6 (Mpc/h)^3 with circular velocities in excess of 450 km/s, and more than one per 10^4.5 (Mpc/h)^3 with circular velocities in excess of 350 km/s). Finally, adding baryonic mass and using a simple model for halo adiabatic contraction, we find remarkable agreement with the velocity dispersion functions inferred observationally by Bezanson et al. (2011) up to z ~ 1 and down to about 220 km/s.

preprint2013arXiv

Tracing Galaxies Through Cosmic Time with Number Density Selection

A central challenge in observational studies of galaxy formation is how to associate progenitor galaxies with their descendants at lower redshifts. One promising approach is to link galaxies at fixed number density, rather than fixed luminosity or mass. This method is effective if stellar mass rank order is broadly conserved through cosmic time. In this paper, we use the Guo et al. (2011) semi-analytical model to analyze under what circumstances this assumption is valid in the context of a cosmological simulation. Specifically, we select progenitor galaxies at a constant number density and compare the stellar mass evolution of their descendants to the evolution at a constant number density. The median stellar mass of the descendants increases by a factor of four (0.6 dex) from z = 3 to z = 0. Constant number density selection reproduces this to within 40% (0.15 dex) over a wide range of number densities. We show that the discrepancy primarily results from scatter in the stellar mass growth rates and merging. After applying simple, observationally-based corrections for these processes, the discrepancy is reduced to 12% (0.05 dex). We conclude that number density selection can be used to predict the median descendant mass of high-redshift progenitor galaxies. The main uncertainty in this study is that semi-analytical models do not reproduce the observed mass evolution of galaxies, which makes the quantitative aggregate effects of star formation, merging, and quenching on the rank order of galaxies somewhat uncertain.

preprint2013arXiv

Tracing the Mass Growth and Star Formation Rate Evolution of Massive Galaxies from z~6 to z~1 in the Hubble Ultra-Deep Field

We present an analysis of $\sim$1500 H160-selected photometric galaxies detected to a limiting magnitude of 27.8 in the HUDF, using imaging from the HST WFC3/IR camera in combination with archival UV, optical, and NIR imaging. We fit photometric redshifts and stellar population estimates for all galaxies with well-determined Spitzer IRAC fluxes, allowing for the determination of the cumulative mass function within the range $1<z<6$. By selecting samples of galaxies at a constant cumulative number density, we explore the co-evolution of stellar masses and star formation rates from z$\sim$6. We find a steady increase in the SFRs of galaxies at constant number density from z$\sim$6 to z$\sim$3. The peak epoch of star formation is found to shift to later times for galaxies with increasing number densities, in agreement with the expectations from cosmic downsizing. The observed SFRs can fully account for the mass growth to z$\sim$2 amongst galaxies with cumulative number densities greater than 10$^{-3.5}$ Mpc$^{-3}$. For galaxies with a lower constant number density we find the observed stellar masses are $\sim$3 times greater than that which may be accounted for by the observed star formation alone at late times, implying that growth from mergers plays an important role at $z<2$. We additionally observe a decreasing sSFR, equivalent to approximately one order of magnitude, from z$\sim$6 to z$\sim$2 amongst galaxies with number densities less than 10$^{-3.5}$ Mpc$^{-3}$ along with significant evidence that at any redshift the sSFR is higher for galaxies at higher number density. The combination of these findings can qualitatively explain the previous findings of a sSFR plateau at high redshift. Tracing the evolution of the fraction of quiescent galaxies for samples matched in cumulative number density over this redshift range, we find no unambiguous examples of quiescent galaxies at $z>4$.

preprint2012arXiv

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

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

preprint2012arXiv

A Nearby Analog of z~2 Compact Quiescent Galaxies with a Rotating Disk

Recent studies have identified a population of compact quiescent galaxies at z\sim2. These galaxies are very rare today and establishing the existence of a nearby analog could allow us to study its structure in greater detail than is possible at high redshift. Here we present such a local analog, NGC 5845, which has a dynamical mass of M_dyn = 4.3\pm0.6\times10^10 M_sun and an effective radius of only r_e = 0.45\pm0.05kpc. We study the structure and kinematics with HST/WFPC2 data and previously published spatially resolved kinematics. We find that NGC 5845 is similar to compact quiescent galaxies at z\sim2 in terms of size versus dynamical mass (r_e-M_dyn), effective velocity dispersion versus size (sigma_e-r_e), and effective velocity dispersion versus dynamical mass (sigma_e-M_dyn). The galaxy has a prominent rotating disk evident in both the photometry and the kinematics: it extends to well beyond \geq1/3 effective radius and contribute to \geq1/4 of the total light of the galaxy. Our results lend support to the idea that a fraction of z\sim2 compact galaxies have prominent disks and positive mass-to-light ratio gradients, although we caution that NGC 5845 may have had a different formation history than the more massive compact quiescent galaxies at z\sim2.

preprint2012arXiv

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

We report the discovery of a massive ultra-compact quiescent galaxy that has been strongly-lensed into multiple images by a foreground galaxy at z = 0.960. This system was serendipitously discovered as a set of extremely Ks-bright high-redshift galaxies with red J - Ks colors using new data from the UltraVISTA YJHKs near-infrared survey. The system was also previously identified as an optically-faint lens/source system using the COSMOS ACS imaging by Faure et al. (2008, 2011). Photometric redshifts for the three brightest images of the source galaxy determined from twenty-seven band photometry place the source at z = 2.4 +/- 0.1. We provide an updated lens model for the system which is a good fit to the positions and morphologies of the galaxies in the ACS image. The lens model implies that the magnification of the three brightest images is a factor of 4 - 5. We use the lens model, combined with the Ks-band image to constrain the size and Sersic profile of the galaxy. The best-fit model is an ultra-compact galaxy (Re = 0.64^{+0.08}_{-0.18} kpc, lensing-corrected), with a Sersic profile that is intermediate between a disk and bulge profile (n = 2.2^{+2.3}_{-0.9}). We present aperture photometry for the source galaxy images which have been corrected for flux contamination from the central lens. The best-fit stellar population model is a massive galaxy (Log(M_{star}/M_{sol}) = 10.8^{+0.1}_{-0.1}, lensing-corrected) with an age of 1.0^{+1.0}_{-0.4} Gyr, moderate dust extinction (Av = 0.8^{+0.5}_{-0.6}), and a low specific star formation rate (Log(SSFR) < -11.0 yr^{-1}). This is typical of massive "red-and-dead" galaxies at this redshift and confirms that this source is the first bona fide strongly-lensed massive ultra-compact quiescent galaxy to be discovered. We conclude with a discussion of the prospects of finding a larger sample of these galaxies.

preprint2012arXiv

Counting Low-Mass Stars in Integrated Light

Low-mass stars (M<0.4Msun) are thought to comprise the bulk of the stellar mass of galaxies but they constitute only of order a percent of the bolometric luminosity of an old stellar population. Directly estimating the number of low-mass stars from integrated flux measurements of old stellar systems is therefore possible but very challenging given the numerous variables that can affect the light at the percent level. Here we present a new population synthesis model created specifically for the purpose of measuring the low-mass initial mass function (IMF) down to ~0.1Msun for metal-rich stellar populations with ages in the range 3-13.5 Gyr. Our fiducial model is based on the synthesis of three separate isochrones and a combination of optical and near-IR empirical stellar libraries in order to produce integrated light spectra over the wavelength interval 0.35mu<lambda<2.4mu at a resolving power of R~2000. New synthetic stellar atmospheres and spectra have been computed in order to model the spectral variations due to changes in individual elemental abundances including C, N, Na, Mg, Si, Ca, Ti, Fe, and generic alpha elements. We demonstrate the power of combining blue spectral features with surface gravity-sensitive near-IR features in order to simultaneously constrain the low-mass IMF, stellar population age, metallicity, and abundance pattern from integrated light measurements. Finally, we show that the shape of the low-mass IMF can also be directly constrained by employing a suite of surface gravity-sensitive spectral features, each of which is most sensitive to a particular mass interval.

preprint2012arXiv

Evolution of Quiescent and Star-Forming Galaxies Since z~1.5 as a Function of Their Velocity Dispersions

We measure stellar masses and structural parameters for 5,500 quiescent and 20,000 star-forming galaxies at 0.3<z\leq1.5 in the Newfirm Medium Band Survey COSMOS and UKIDSS UDS fields. We combine these measurements to infer velocity dispersions and determine how the number density of galaxies at fixed inferred dispersion, or the Velocity Dispersion Function (VDF), evolves with time for each population. We show that the number of galaxies with high velocity dispersions appears to be surprisingly stable with time, regardless of their star formation history. Furthermore, the overall VDF for star-forming galaxies is constant with redshift, extending down to the lowest velocity dispersions probed by this study. The only galaxy population showing strong evolution are quiescent galaxies with low inferred dispersions, whose number density increases by a factor of ~4 since z=1.5. This build-up leads to an evolution in the quiescent fraction of galaxies such that the threshold dispersion above which quiescent galaxies dominate the counts moves to lower velocity dispersion with time. We show that our results are qualitatively consistent with a simple model in which star-forming galaxies quench and are added to the quiescent population. In order to compensate for the migration into the quiescent population, the velocity dispersions of star-forming galaxies must increase, with a rate that increases with dispersion.

preprint2012arXiv

First Results from Z-FOURGE: Discovery of a Candidate Cluster at z = 2.2 in COSMOS

We report the first results from the Z-FOURGE survey: the discovery of a candidate galaxy cluster at z = 2.2 consisting of two compact overdensities with red galaxies detected at >20σ above the mean surface density. The discovery was made possible by a new deep (Ks < 24.8 AB 5σ) Magellan/FourStar near-IR imaging survey with 5 custom medium-bandwidth filters. The filters pinpoint the location of the Balmer/4000Å break in evolved stellar populations at 1.5 < z < 3.5, yielding significantly more accurate photometric redshifts than possible with broadband imaging alone. The overdensities are within 1' of each other in the COSMOS field and appear to be embedded in a larger structure that contains at least one additional overdensity (\sim10σ). Considering the global properties of the overdensities, the z = 2.2 system appears to be the most distant example of a galaxy cluster with a population of red galaxies. A comparison to a large ΛCDM simulation suggests that the system may consist of merging subclusters, with properties in between those of z > 2 protoclusters with more diffuse distributions of blue galaxies and the lower-redshift galaxy clusters with prominent red sequences. The structure is completely absent in public optical catalogs in COSMOS and only weakly visible in a shallower near-IR survey. The discovery showcases the potential of deep near-IR surveys with medium-band filters to advance the understanding of environment and galaxy evolution at z > 1.5.

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

The Stellar Initial Mass Function in Early-Type Galaxies From Absorption Line Spectroscopy. I. Data and Empirical Trends

The strength of gravity-sensitive absorption lines in the integrated light of old stellar populations is one of the few direct probes of the stellar initial mass function (IMF) outside of the Milky Way. Owing to the advent of fully depleted CCDs with little or no fringing it has recently become possible to obtain accurate measurements of these features. Here we present spectra covering the wavelength ranges 0.35 - 0.55 micron and 0.72 - 1.03 micron for the bulge of M31 and 34 early-type galaxies from the SAURON sample, obtained with the Low Resolution Imaging Spectrometer on Keck. The signal-to-noise ratio is >200 per Angstrom out to 1 micron, which is sufficient to measure gravity-sensitive features for individual galaxies and to determine how they depend on other properties of the galaxies. Combining the new data with previously obtained spectra for globular clusters in M31 and the most massive elliptical galaxies in the Virgo cluster we find that the dwarf-sensitive Na I doublet and the FeH Wing-Ford band increase systematically with velocity dispersion, while the giant-sensitive Ca II triplet decreases with dispersion. These trends are consistent with a varying IMF, such that galaxies with deeper potential wells have more dwarf-enriched mass functions. In a companion paper (Conroy & van Dokkum 2012) we use a comprehensive stellar population synthesis model to demonstrate that IMF effects can be separated from age and abundance variations and quantify the IMF variation among early-type galaxies.

preprint2012arXiv

The Stellar Initial Mass Function in Early-Type Galaxies From Absorption Line Spectroscopy. II. Results

The spectral absorption lines in early-type galaxies contain a wealth of information regarding the detailed abundance pattern, star formation history, and stellar initial mass function (IMF) of the underlying stellar population. Using our new population synthesis model that accounts for the effect of variable abundance ratios of 11 elements, we analyze very high quality absorption line spectra of 38 early-type galaxies and the nuclear bulge of M31. These data extend to 1um and they therefore include the IMF-sensitive spectral features NaI, CaII, and FeH at 0.82um, 0.86um and 0.99um, respectively. The models fit the data well, with typical rms residuals ~1%. Strong constraints on the IMF and therefore the stellar mass-to-light ratio, (M/L)_stars, are derived for individual galaxies. We find that the IMF becomes increasingly bottom-heavy with increasing velocity dispersion and [Mg/Fe]. At the lowest dispersions and [Mg/Fe] values the derived IMF is consistent with the Milky Way IMF, while at the highest dispersions and [Mg/Fe] values the derived IMF contains more low-mass stars (is more bottom-heavy) than even a Salpeter IMF. Our best-fit (M/L)_stars values do not exceed dynamically-based M/L values. We also apply our models to stacked spectra of four metal-rich globular clusters in M31 and find an (M/L)_stars that implies fewer low-mass stars than a Milky Way IMF, again agreeing with dynamical constraints. We discuss other possible explanations for the observed trends and conclude that variation in the IMF is the simplest and most plausible.

preprint2011arXiv

Confirmation of Enhanced Dwarf-sensitive Absorption Features in the Spectra of Massive Elliptical Galaxies: Further Evidence for a Non-universal Initial Mass Function

We recently found that massive cluster elliptical galaxies have strong Na I 8183,8195 and FeH 9916 Wing-Ford band absorption, indicating the presence of a very large population of stars with masses <~ 0.3 M_sun. Here we test this result by comparing the elliptical galaxy spectra to those of luminous globular clusters associated with M31. These globular clusters have similar metallicities, abundance ratios and ages as massive elliptical galaxies but their low dynamical mass-to-light ratios rule out steep stellar initial mass functions (IMFs). From high quality Keck spectra we find that the dwarf-sensitive absorption lines in globular clusters are significantly weaker than in elliptical galaxies, and consistent with normal IMFs. The differences in the Na I and Wing-Ford indices are 0.027 +- 0.007 mag and 0.017 +- 0.006 mag respectively. We directly compare the two classes of objects by subtracting the averaged globular cluster spectrum from the averaged elliptical galaxy spectrum. The difference spectrum is well fit by the difference between a stellar population synthesis model with a bottom-heavy IMF and one with a bottom-light IMF. We speculate that the slope of the IMF may vary with velocity dispersion, although it is not yet clear what physical mechanism would be responsible for such a relation.

preprint2011arXiv

The faint stellar halos of massive red galaxies from stacks of more than 42000 SDSS LRG images

We study the properties of massive galaxies at an average redshift of z~0.34 through stacking more than 42000 images of Luminous Red Galaxies from the Sloan Digital Sky Survey. This is the largest dataset ever used for such an analysis and it allows us to explore the outskirts of massive red galaxies at unprecedented physical scales. Our image stacks extend farther than 400 kpc, where the r-band profile surface brightness reaches 30 mag arcsec-2. This analysis confirms that the stellar bodies of luminous red galaxies follow a simple Sersic profile out to 100 kpc. At larger radii the profiles deviate from the best-fit Sersic models and exhibit extra light in the g, r, i and z-band stacks. This excess light can probably be attributed to unresolved intragroup or intracluster light or a change in the light profile itself. We further show that standard analyses of SDSS-depth images typically miss 20% of the total stellar light and underestimate the size of LRGs by 10% compared to our best fit r-band Sersic model of n=5.5 and r_e=13.1 kpc. If the excess light at r>100 kpc is considered to be part of the galaxy, the best fit r-band Sersic parameters are n=5.8 and r_e=13.6 kpc. In addition we study the radially dependent stack ellipticity and find an increase with radius from e=0.25 at r=10 kpc to e=0.3 at r=100 kpc. This provides support that the stellar light that we trace out to at least 100 kpc is physically associated with the galaxies themselves and may confirm that the halos of individual LRGs have higher ellipticities than their central parts. Lastly we show that the broadband color gradients of the stacked images are flat beyond roughly 40 kpc, suggesting that the stellar populations do not vary significantly with radius in the outer parts of massive ellipticals.

preprint2011arXiv

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

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

preprint2010arXiv

A Substantial Population of Low Mass Stars in Luminous Elliptical Galaxies

The stellar initial mass function (IMF) describes the mass distribution of stars at the time of their formation and is of fundamental importance for many areas of astrophysics. The IMF is reasonably well constrained in the disk of the Milky Way but we have very little direct information on the form of the IMF in other galaxies and at earlier cosmic epochs. Here we investigate the stellar mass function in elliptical galaxies by measuring the strength of the Na I doublet and the Wing-Ford molecular FeH band in their spectra. These lines are strong in stars with masses <0.3 Msun and weak or absent in all other types of stars. We unambiguously detect both signatures, consistent with previous studies that were based on data of lower signal-to-noise ratio. The direct detection of the light of low mass stars implies that they are very abundant in elliptical galaxies, making up >80% of the total number of stars and contributing >60% of the total stellar mass. We infer that the IMF in massive star-forming galaxies in the early Universe produced many more low mass stars than the IMF in the Milky Way disk, and was probably slightly steeper than the Salpeter form in the mass range 0.1 - 1 Msun.

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 Spectral Energy Distribution of Post-Starburst Galaxies in the NEWFIRM Medium-Band Survey: A Low Contribution from TP-AGB Stars

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

preprint2010arXiv

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

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

preprint2010arXiv

z~7 galaxy candidates from NICMOS observations over the HDF South and the CDF-S and HDF-N GOODS fields

We use ~88 arcmin**2 of deep (>~26.5 mag at 5 sigma) NICMOS data over the two GOODS fields and the HDF South to conduct a search for bright z>~7 galaxy candidates. This search takes advantage of an efficient preselection over 58 arcmin**2 of NICMOS H-band data where only plausible z>~7 candidates are followed up with NICMOS J-band observations. ~248 arcmin**2 of deep ground-based near-infrared data (>~25.5 mag, 5 sigma) is also considered in the search. In total, we report 15 z-dropout candidates over this area -- 7 of which are new to these search fields. Two possible z~9 J-dropout candidates are also found, but seem unlikely to correspond to z~9 galaxies. The present z~9 search is used to set upper limits on the prevalence of such sources. Rigorous testing is undertaken to establish the level of contamination of our selections by photometric scatter, low mass stars, supernovae (SNe), and spurious sources. The estimated contamination rate of our z~7 selection is ~24%. Through careful simulations, the effective volume available to our z>~7 selections is estimated and used to establish constraints on the volume density of luminous (L*(z=3), or -21 mag) galaxies from these searches. We find that the volume density of luminous star-forming galaxies at z~7 is 13_{-5}^{+8}x lower than at z~4 and >25x lower (1 sigma) at z~9 than at z~4. This is the most stringent constraint yet available on the volume density of >~L* galaxies at z~9. The present wide-area, multi-field search limits cosmic variance to <20%. The evolution we find at the bright end of the UV LF is similar to that found from recent Subaru Suprime-Cam, HAWK-I or ERS WFC3/IR searches. The present paper also includes a complete summary of our final z~7 z-dropout sample (18 candidates) identified from all NICMOS observations to date (over the two GOODS fields, the HUDF, galaxy clusters).

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.

preprint2008arXiv

Assessing the Predictive Power of Galaxy Formation Models with the Rest-Frame Optical Luminosity Functions at 2.0<z<3.3

We compare recently measured rest-frame V-band luminosity functions (LFs) of galaxies at redshifts 2.0<z<3.3 to predictions of semianalytic models by De Lucia & Blaizot and Bower et al. and hydrodynamic simulations by Dave et al. The models succeed for some luminosity and redshift ranges and fail for others. A notable success is that the Bower et al. model provides a good match to the observed LF at z~3. However, all models predict an increase with time of the rest-frame V-band luminosity density, whereas the observations show a decrease. The models also have difficulty matching the observed rest-frame colors of galaxies. In all models the luminosity density of red galaxies increases sharply from z~3 to z~2.2, whereas it is approximately constant in the observations. Conversely, in the models the luminosity density of blue galaxies is approximately constant, whereas it decreases in the observations. These discrepancies cannot be entirely remedied by changing the treatment of dust and suggest that current models do not yet provide an adequate description of galaxy formation and evolution since z~3.

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.

preprint2007arXiv

A Keck Spectroscopic Survey of MS1054-03 (z=0.83): Forming the Red Sequence

Using a magnitude-limited, spectroscopic survey of the X-ray luminous galaxy cluster MS1054-03, we isolate 153 cluster galaxies and measure MS1054's redshift and velocity dispersion to be z=0.8307 and sigma_z=1156 km/s. The absorption-line, post-starburst (``E+A''), and emission-line galaxies respectively make up 63%, 15%, and 23% of the cluster population. With photometry from HST/ACS, we find that the absorption-line members define an exceptionally tight red sequence over a span of ~3.5 magnitudes in i775: their intrinsic scatter in (V-i) color is only 0.048, corresponding to a (U-B)_z scatter of 0.041. Their color scatter is comparable to that of the ellipticals (sigma_Vi=0.055), but measurably smaller than that of the combined E+S0 sample (sigma_Vi=0.072). The color scatter of MS1054's absorption-line population is approximately twice that of the ellipticals in Coma; this difference is consistent with passive evolution where most of the absorption-line members (>75%) formed by z~2, and all of them by z~1.2. For red members, we find a trend (>95% confidence) of weakening Hdelta absorption with redder colors that we conclude is due to age: in MS1054, the color scatter on the red sequence is driven by differences in mean stellar age of up to ~1.5 Gyr. We also generate composite spectra and estimate that the average S0 in MS1054 is ~0.5-1 Gyr younger than the average elliptical; this difference in mean stellar age is mainly due to a number of S0s that are blue (18%) and/or are post-starburst systems (21%).

preprint2007arXiv

Evidence of Cosmic Evolution of the Stellar Initial Mass Function

Theoretical arguments and indirect observational evidence suggest that the stellar initial mass function (IMF) may evolve with time, such that it is more weighted toward high mass stars at higher redshift. Here we test this idea by comparing the rate of luminosity evolution of massive early-type galaxies in clusters at 0.02<z<0.83 to the rate of their color evolution. A combined fit to the rest-frame U-V color evolution and the previously measured evolution of the M/L_B ratio gives x~-0.3 for the logarithmic slope of the IMF in the region around 1 Solar mass, flatter than the present-day value in the Milky Way disk of x=1.3. The best-fitting luminosity-weighted formation redshift of the stars in massive cluster galaxies is ~3.7, and a possible interpretation is that the characteristic mass m_c had a value of ~2 Solar masses at z~4 (compared to m_c~0.1 Solar masses today), in qualitative agreement with models in which the characteristic mass is a function of the Jeans mass in molecular clouds. Such a ``bottom-light'' IMF for massive cluster galaxies has significant implications for the interpretation of measurements of galaxy formation and evolution. Applying a simple form of IMF evolution to literature data, we find that the volume-averaged star formation rate at high redshift may have been overestimated (by a factor of 3-4 at z>4), and the cosmic star formation history may have a fairly well-defined peak at z~1.5. The M/L_V ratios of galaxies are less affected than their star formation rates, and future data on the stellar mass density at z>3 will provide further constraints on IMF evolution. The formal errors likely underestimate the uncertainties, and confirmation of thes results requires a larger sample of clusters and the inclusion of redder rest-frame colors in the analysis.

preprint2007arXiv

Spitzer Constraints on the Stellar Populations of Lyman-Alpha Emitting Galaxies at z = 3.1

We investigate the stellar populations of a sample of 162 Lyman-alpha emitting galaxies (LAEs) at z = 3.1 in the Extended Chandra Deep Field South, using deep Spitzer IRAC data available from the GOODS and SIMPLE surveys to derive reliable stellar population estimates. We divide the LAEs according to their rest-frame near-IR luminosities into IRAC-detected and IRAC-undetected samples. About 70% of the LAEs are undetected in 3.6 micron down to [3.6] = 25.2 AB. Stacking analysis reveals that the average stellar population of the IRAC-undetected sample has an age of ~ 200 Myr and a mass of ~ 3x10^8 solar masses, consistent with the expectation that LAEs are mostly young and low-mass galaxies. On the other hand, the IRAC-detected LAEs are on average significantly older and more massive, with an average age > 1 Gyr and mass ~ 10^10 solar masses. Comparing the IRAC colors and magnitudes of the LAEs to z ~ 3 Lyman break galaxies (LBGs) shows that the IRAC-detected LAEs lie at the faint blue end of the LBG color-magnitude distribution, suggesting that IRAC-detected LAEs may be the low mass extension of the LBG population. We also present tentative evidence for a small fraction (~ 5%) of obscured AGN within the LAE sample. Our results suggest that LAEs posses a wide range of ages and masses. Additionally, the presence of evolved stellar populations inside LAEs suggests that the Lyman-alpha luminous phase of galaxies may either be a long-lasting or recurring phenomenon.

preprint2007arXiv

The Warm Spitzer Mission: Prospects for Studies of the Distant Universe

IRAC excels at detecting distant objects. Due to a combination of the shapes of the spectral energy distributions of galaxies and the low background achieved from space, IRAC reaches greater depth in comparable exposure time at 3.6 and 4.5 micron than any ground- or space-based facility currently can at 2.2 micron. Furthermore, the longer wavelengths probed by IRAC enable studies of the rest-frame optical and near-infrared light of galaxies and AGN to much higher redshift than is possible from the ground. This white paper explores the merits of different survey strategies for studying the distant universe during the warm mission. A three-tiered approach serves a wide range of science goals and uses the spacecraft effectively: 1) an ultra-deep survey of ~0.04 square degrees to a depth of ~250 hrs (in conjunction with an HST/WFC3 program), to study the Universe at 7<z<14; 2) a survey of ~2 square degrees to the GOODS depth of 20 hrs, to identify luminous galaxies at z>6 and characterize the relation between the build-up of dark matter halos and their constituent galaxies at 2<z<6, and 3) a 500 square degree survey to the SWIRE depth of 120 s, to systematically study large scale structure at 1<z<2 and characterize high redshift AGN. One or more of these programs could conceivably be implemented by the SSC, following the example of the Hubble Deep Field campaigns. As priorities in this field continuously shift it is also crucial that a fraction of the exposure time remains unassigned, thus enabling science that will reflect the frontiers of 2010 and beyond rather than those of 2007.

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.

preprint2004arXiv

Stellar Populations and Kinematics of Red Galaxies at z>2: Implications for the Formation of Massive Galaxies

We recently identified a substantial population of galaxies at z>2 with red rest-frame optical colors. These distant red galaxies (DRGs) are efficiently selected by the simple observed color criterion J-K>2.3. In this paper we present NIR spectroscopy with Keck/NIRSPEC of six DRGs at 2.4<z<3.2. We detect continuum emission and emission lines of all observed galaxies. Equivalent widths of H alpha are 20-30 Ang, smaller than measured for LBGs and nearby LIRGs, and comparable to normal nearby galaxies. The modest equivalent widths imply that the galaxies either have a decreasing star formation rate, or that they are very dusty. Fitting both the photometry and the H alpha lines, we find continuum extinction A_V=1-2 mag, ages 1-2.5 Gyr, star formation rates 200-400 solar masses/yr, and stellar masses 1-5x10^11 solar masses for models with constant star formation rates. From [NII]/H alpha ratios we infer that the metallicities are high, 1-1.5 x Solar. For four galaxies we can determine line widths from the optical emission lines. The widths are high, ranging from 130-240 km/s, and by combining data for LBGs and DRGs we find significant correlations between linewidth and restframe U-V color, and between linewidth and stellar mass. The latter correlation has a similar slope and offset as the ``baryonic Tully-Fisher relation'' for nearby galaxies. The median dynamical mass is ~2x10^11 solar masses, supporting the high stellar masses inferred from the photometry. We find that the median M/L_V ~ 0.8, a factor of ~5 higher than measured for LBGs. We infer from our small sample that DRGs are dustier, more metal rich, more massive, and have higher ages than z=3 LBGs of the same rest-frame V-band luminosity. Their high M/L ratios imply that they contribute significantly to the stellar mass density at z~2.5. [ABRIDGED]

preprint1996arXiv

The Fundamenal Plane in CL0024 at z=0.4: Implications for the Evolution of the Mass-to-Light Ratio

We present results on the Fundamental Plane of early-type galaxies in the rich cluster CL 0024+16 at z=0.391. The galaxies satisfy a tight Fundamental Plane relation which is similar to that at low redshift. The scatter is 15 % in log re. The data show that massive early-type galaxies existed at z=0.4. The evolution of the M/L ratio is derived from the Fundamental Plane, by a comparison with Coma. The M/L ratio increases by 31 +- 12 % between z=0.391, and z=0.023. The evolution is low when compared to models for stellar populations. The expected evolution depends on the IMF, q0, and the formation redshift of the galaxies. The data are in agreement with high formation redshifts. The most serious bias in the modeling may be the progenitor bias: if the progenitors of some current day early-types are spirals at z=0.4, they would not be included in the sample, and the sample would be biased towards the oldest galaxies. More data are needed to measure the evolution of the Fundamental Plane more precisely, and its scatter. There is a hint that the form of the Fundamental Plane changes with redshift, and this needs to be determined better.