Source author record

Charlie Conroy

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.IM astro-ph.HE astro-ph.EP astro-ph hep-ph

Catalog footprint

What is connected

117works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2023arXiv

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

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 massive AGB donor in Scutum X-1: Identification of the first Mira variable in an X-ray binary

The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observations, ii) the first broadband infrared spectrum ($\approx1-5\,μ$m; $R\approx 2000$) obtained with SpeX on the NASA Infrared Telescope Facility and iii) $J$-band light curve from the Palomar Gattini-IR survey. The infrared spectrum is characterized by i) deep water absorption features (H$_2$O index $\approx 40$%), ii) strong TiO, VO and CO features, and iii) weak/absent CN lines. We show that these features are inconsistent with known red supergiants, but suggest a M8-9 III type O-rich Mira donor star. We report the discovery of large amplitude ($ΔJ\approx3.5$ mag) periodic photometric variability suggesting a pulsation period of $621\pm36\,{\rm(systematic)}\pm8\,{\rm(statistical)}$ days, which we use to constrain the donor to be a relatively luminous Mira ($M_K=-8.6\pm0.3$ mag) at a distance of $3.6^{+0.8}_{-0.7}$ kpc. Comparing these characteristics to recent models, we find the donor to be consistent with a $\approx 3-5$ M$_\odot$ star at an age of $\approx 0.1-0.3$ Gyr. Together, we show that Sct X-1 was previously mis-classified as an evolved High Mass X-ray Binary; instead it is an intermediate mass system with the first confirmed Mira donor in an X-ray binary. We discuss the implications of Mira donors in symbiotic X-ray binaries, and highlight the potential of wide-field infrared time domain surveys and broadband infrared spectroscopy to unveil their demographics.

preprint2022arXiv

A Tilt in the Dark Matter Halo of the Galaxy

Recent observations of the stellar halo have uncovered the debris of an ancient merger, Gaia-Sausage-Enceladus, estimated to have occurred ~8 Gyr ago. Follow-up studies have associated GSE with a large-scale tilt in the stellar halo that links two well-known stellar over-densities in diagonally opposing octants of the Galaxy (the Hercules-Aquila Cloud and Virgo Overdensity; HAC and VOD). In this paper, we study the plausibility of such unmixed merger debris persisting over several Gyr in the Galactic halo. We employ the simulated stellar halo from Naidu et al. (2021), which reproduces several key properties of the merger remnant, including the large-scale tilt. By integrating the orbits of these simulated stellar halo particles, we show that adoption of a spherical halo potential results in rapid phase mixing of the asymmetry. However, adopting a tilted halo potential preserves the initial asymmetry in the stellar halo for many Gyr. The asymmetry is preserved even when a realistic growing disk is added to the potential. These results suggest that HAC and VOD are long-lived structures that are associated with GSE and that the dark matter halo of the Galaxy is tilted with respect to the disk and aligned in the direction of HAC-VOD. Such halo-disk misalignment is common in modern cosmological simulations. Lastly, we study the relationship between the local and global stellar halo in light of a tilted global halo comprised of highly radial orbits. We find that the local halo offers a dynamically biased view of the global halo due to its displacement from the Galactic Center.

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

Deriving ages and horizontal branch properties of integrated stellar populations

A major source of uncertainty in the age determination of old ($\sim10$ Gyr) integrated stellar populations is the presence of hot horizontal branch (HB) stars. Here, we describe a simple approach to tackle this problem, and show the performance of this technique that simultaneously models the age, abundances and HB properties of integrated stellar populations. For this we compare the results found during the fits of the integrated spectra of a sample of stellar population benchmarks, against the values obtained from the analysis of their resolved CMDs. We find that the ages derived from our spectral fits for most (26/32) of our targets are within 0.1 dex to their CMDs values. Similarly, for the majority of the targets in our sample we are able to recover successfully the flux contribution from hot HB stars (within $\sim0.15 {\rm ~dex}$ for 18/24 targets) and their mean temperature (14/24 targets within $\sim30 \%$). Finally, we present a diagnostic that can be used to detect spurious solutions in age, that will help identify the few cases when this method fails. These results open a new window for the detailed study of globular clusters beyond the Local Group.

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

How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

One of the most common methods for inferring galaxy attenuation curves is via spectral energy distribution (SED) modeling, where the dust attenuation properties are modeled simultaneously with other galaxy physical properties. In this paper, we assess the ability of SED modeling to infer these dust attenuation curves from broadband photometry, and suggest a new flexible model that greatly improves the accuracy of attenuation curve derivations. To do this, we fit mock SEDs generated from the Simba cosmological simulation with the Prospector SED fitting code. We consider the impact of the commonly-assumed uniform screen model and introduce a new non-uniform screen model parameterized by the fraction of unobscured stellar light. This non-uniform screen model allows for a non-zero fraction of stellar light to remain unattenuated, resulting in a more flexible attenuation curve shape by decoupling the shape of the UV attenuation curve from the optical attenuation curve. The ability to constrain the dust attenuation curve is significantly improved with the use of a non-uniform screen model, with the median offset in UV attenuation decreasing from $-0.30$ dex with a uniform screen model to $-0.17$ dex with the non-uniform screen model. With this increase in dust attenuation modeling accuracy, we also improve the star formation rates (SFRs) inferred with the non-uniform screen model, decreasing the SFR offset on average by $0.12$ dex. We discuss the efficacy of this new model, focusing on caveats with modeling star-dust geometries and the constraining power of available SED observations.

preprint2022arXiv

Live Fast, Die $α$-Enhanced: The Mass-Metallicity-$α$ Relation of the Milky Way's Disrupted Dwarf Galaxies

The Milky Way's satellite galaxies ("surviving dwarfs") have been studied for decades as unique probes of chemical evolution in the low-mass regime. Here we extend such studies to the "disrupted dwarfs", whose debris constitutes the stellar halo. We present abundances ([Fe/H], [$α$/Fe]) and stellar masses for nine disrupted dwarfs with $M_{\star}\approx10^{6}-10^{9}M_{\odot}$ from the H3 Survey (Sagittarius, $Gaia$-Sausage-Enceladus, Helmi Streams, Sequoia, Wukong/LMS-1, Cetus, Thamnos, I'itoi, Orphan/Chenab). The surviving and disrupted dwarfs are chemically distinct: at fixed mass, the disrupted dwarfs are systematically metal-poor and $α$-enhanced. The disrupted dwarfs define a mass-metallicity relation (MZR) with a similar slope as the $z=0$ MZR followed by the surviving dwarfs, but offset to lower metallicities by $Δ$[Fe/H]$\approx0.3-0.4$ dex. Dwarfs with larger offsets from the $z=0$ MZR are more $α$-enhanced. In simulations as well as observations, galaxies with higher $Δ$[Fe/H] formed at higher redshifts -- exploiting this, we infer the disrupted dwarfs have typical star-formation truncation redshifts of $z_{\rm{trunc}}{\sim}1-2$. We compare the chemically inferred $z_{\rm{trunc}}$ with dynamically inferred accretion redshifts and find almost all dwarfs are quenched only after accretion. The differences between disrupted and surviving dwarfs are likely because the disrupted dwarfs assembled their mass rapidly, at higher redshifts, and within denser dark matter halos that formed closer to the Galaxy. Our results place novel archaeological constraints on low-mass galaxies inaccessible to direct high-$z$ studies: (i) the redshift evolution of the MZR along parallel tracks but offset to lower metallicities extends to $M_{\star}\approx10^{6}-10^{9}M_{\odot}$; (ii) galaxies at $z\approx2-3$ are $α$-enhanced with [$α$/Fe]$\approx0.4$.

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

SDSS-IV MaNGA: How the stellar populations of passive central galaxies depend on stellar and halo mass

We analyze spatially resolved and co-added SDSS-IV MaNGA spectra with signal-to-noise ~100 from 2200 passive central galaxies (z~0.05) to understand how central galaxy assembly depends on stellar mass (M*) and halo mass (Mh). We control for systematic errors in Mh by employing a new group catalog from Tinker (2020a,b) and the widely-used Yang et al. (2007) catalog. At fixed M*, the strength of several stellar absorption features varies systematically with Mh. Completely model-free, this is one of the first indications that the stellar populations of centrals with identical M* are affected by the properties of their host halos. To interpret these variations, we applied full spectral fitting with the code alf. At fixed M*, centrals in more massive halos are older, show lower [Fe/H], and have higher [Mg/Fe] with 3.5 sigma confidence. We conclude that halos not only dictate how much M* galaxies assemble, but also modulate their chemical enrichment histories. Turning to our analysis at fixed Mh, high-M* centrals are older, show lower [Fe/H], and have higher [Mg/Fe] for Mh>10^{12}Msun/h with confidence > 4 sigma. While massive passive galaxies are thought to form early and rapidly, our results are among the first to distinguish these trends at fixed Mh. They suggest that high-M* centrals experienced unique early formation histories, either through enhanced collapse and gas fueling, or because their halos were early-forming and highly concentrated, a possible signal of galaxy assembly bias.

preprint2022arXiv

SRGA J181414.6-225604: A new Galactic symbiotic X-ray binary outburst triggered by an intense mass loss episode of a heavily obscured Mira variable

We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and luminous ($M_K\approx -9.9 \pm 0.2$) O-rich Mira donor star located at a distance of $\approx 14.6^{+2.9}_{-2.3}$ kpc. Combining multi-color photometric data over the last $\approx 25$ years, we show that the IR counterpart underwent a recent (starting $\approx 800$ days before the X-ray flare) enhanced mass loss (reaching $\approx 2.1 \times 10^{-5}$ M$_\odot$ yr$^{-1}$) episode resulting in an expanding dust shell obscuring the underlying star. Multi-epoch follow-up from Swift, NICER and NuSTAR reveal a $\approx 200$ day long X-ray outburst reaching a peak luminosity of $L_X \approx 2.5 \times 10^{36}$ erg s$^{-1}$, characterized by a heavily absorbed ($N_{\rm H} \approx 6\times 10^{22}$ cm$^{-2}$) X-ray spectrum consistent with an optically thick Comptonized plasma. The X-ray spectral and timing behavior suggest the presence of clumpy wind accretion together with a dense ionized nebula overabundant in silicate material surrounding the compact object. Together, we show that SRGA J181414.6-225604 is a new symbiotic X-ray binary in outburst, triggered by an intense dust formation episode of a highly evolved donor. Our results offer the first direct confirmation for the speculated connection between enhanced late-stage donor mass loss and active lifetimes of the symbiotic X-ray binaries.

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

Unicorns and Giraffes in the binary zoo: stripped giants with subgiant companions

We analyze two binary systems containing giant stars, V723 Mon ("the Unicorn") and 2M04123153+6738486 ("the Giraffe"). Both giants orbit more massive but less luminous companions, previously proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modeling of the spectra, light curves, and spectral energy distributions robustly constrains the masses, temperatures, and radii of both components: the primaries are luminous, cool giants ($T_{\rm eff,\,giant} = 3,800\,\rm K$ and $4,000\,\rm K$, $R_{\rm giant}= 22.5\,R_{\odot}$ and $25\,R_{\odot}$) with exceptionally low masses ($M_{\rm giant} \approx 0.4\,M_{\odot}$) that likely fill their Roche lobes. The secondaries are only slightly warmer subgiants ($T_{\rm eff,\,2} = 5,800\,\rm K$ and $5,150\,\rm K$, $R_2= 8.3\,R_{\odot}$ and $9\,R_{\odot}$) and thus are consistent with observed UV limits that would rule out main-sequence stars with similar masses ($M_2 \approx 2.8\,M_{\odot}$ and $\approx 1.8\,M_{\odot}$). In the Unicorn, rapid rotation blurs the spectral lines of the subgiant, making it challenging to detect even at wavelengths where it dominates the total light. Both giants have surface abundances indicative of CNO processing and subsequent envelope stripping. The properties of both systems can be reproduced by binary evolution models in which a $1-2\,M_{\odot}$ primary is stripped by a companion as it ascends the giant branch. The fact that the companions are also evolved implies either that the initial mass ratio was very near unity, or that the companions are temporarily inflated due to rapid accretion. The Unicorn and Giraffe offer a window into into a rarely-observed phase of binary evolution preceding the formation of wide-orbit helium white dwarfs, and eventually, compact binaries containing two helium white dwarfs.

preprint2022arXiv

Wide binaries from the H3 survey: the thick disk and halo have similar wide binary fractions

Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-confidence (a contamination rate of 4%) wide binaries and two resolved triples, with binary separations mostly between $10^3$-$10^5$ AU and a lowest [Fe/H] of $-2.7$. Based on their Galactic kinematics, 33 of them are halo wide binaries, and most of those are associated with the accreted Gaia-Sausage-Enceladus galaxy. The wide binary fraction in the thick disk decreases toward the low metallicity end, consistent with the previous findings for the thin disk. Our key finding is that the halo wide binary fraction is consistent with the thick-disk stars at a fixed [Fe/H]. There is no significant dependence of the wide binary fraction on the $α$-captured abundance. Therefore, the wide binary fraction is mainly determined by the iron abundance, not their disk or halo origin nor the $α$-captured abundance. Our results suggest that the formation environments play a major role for the wide binary fraction, instead of other processes like radial migration that only apply to disk stars.

preprint2021arXiv

Evidence from Disrupted Halo Dwarfs that $r$-process Enrichment via Neutron Star Mergers is Delayed by $\gtrsim500$ Myrs

The astrophysical origins of $r$-process elements remain elusive. Neutron star mergers (NSMs) and special classes of core-collapse supernovae (rCCSNe) are leading candidates. Due to these channels' distinct characteristic timescales (rCCSNe: prompt, NSMs: delayed), measuring $r$-process enrichment in galaxies of similar mass, but differing star-formation durations might prove informative. Two recently discovered disrupted dwarfs in the Milky Way's stellar halo, Kraken and \textit{Gaia}-Sausage Enceladus (GSE), afford precisely this opportunity: both have $M_{\star}\approx10^{8}M_{\rm{\odot}}$, but differing star-formation durations of ${\approx}2$ Gyrs and ${\approx}3.6$ Gyrs. Here we present $R\approx50,000$ Magellan/MIKE spectroscopy for 31 stars from these systems, detecting the $r$-process element Eu in all stars. Stars from both systems have similar [Mg/H]$\approx-1$, but Kraken has a median [Eu/Mg]$\approx-0.1$ while GSE has an elevated [Eu/Mg]$\approx0.2$. With simple models we argue NSM enrichment must be delayed by $500-1000$ Myrs to produce this difference. rCCSNe must also contribute, especially at early epochs, otherwise stars formed during the delay period would be Eu-free. In this picture, rCCSNe account for $\approx50\%$ of the Eu in Kraken, $\approx25\%$ in GSE, and $\approx15\%$ in dwarfs with extended star-formation durations like Sagittarius. The inferred delay time for NSM enrichment is $10-100\times$ longer than merger delay times from stellar population synthesis -- this is not necessarily surprising because the enrichment delay includes time taken for NSM ejecta to be incorporated into subsequent generations of stars. For example, this may be due to natal kicks that result in $r$-enriched material deposited far from star-forming gas, which then takes $\approx10^{8}-10^{9}$ years to cool in these galaxies.

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

Stellar labels for hot stars from low-resolution spectra - I. the HotPayne method and results for 330,000 stars from LAMOST DR6

We set out to determine stellar labels from low-resolution survey spectra of hot, OBA stars with effective temperature (Teff) higher than 7500K. This fills a gap in the scientific analysis of large spectroscopic stellar surveys such as LAMOST, which offers spectra for millions of stars at R=1800. We first explore the theoretical information content of such spectra for determining stellar labels, via the Cramér-Rao bound. We show that in the limit of perfect model spectra and observed spectra with S/N of 100, precise estimates are possible for a wide range of stellar labels: not only the effective temperature Teff, surface gravity logg, and projected rotation velocity vsini, but also the micro-turbulence velocity, Helium abundance and the elemental abundances [C/H], [N/H], [O/H], [Si/H], [S/H], and [Fe/H]. Our analysis illustrates that the temperature regime of around 9500K is challenging, as the dominant Balmer and Paschen line strength vary little with Teff. We implement the simultaneous fitting of these 11 stellar labels to LAMOST hot-star spectra using the Payne approach, drawing on Kurucz's ATLAS12/SYNTHE LTE spectra as the underlying models. We then obtain stellar parameter estimates for a sample of about 330,000 hot stars with LAMOST spectra, an increase by about two orders of magnitude in sample size. Among them, about 260,000 have good Gaia parallaxes (S/N>5), and more than 95 percent of them are luminous stars, mostly on the main sequence; the rest reflects lower luminosity evolved stars, such as hot subdwarfs and white dwarfs. We show that the fidelity of the abundance estimates is limited by the systematics of the underlying models, as they do not account for NLTE effects. Finally, we show the detailed distribution of vsini of stars with 8000-15,000K, illustrating that it extends to a sharp cut-off at the critical rotation velocity, across a wide range of temperatures.

preprint2021arXiv

The Mass of the Milky Way from the H3 Survey

The mass of the Milky Way is a critical quantity which, despite decades of research, remains uncertain within a factor of two. Until recently, most studies have used dynamical tracers in the inner regions of the halo, relying on extrapolations to estimate the mass of the Milky Way. In this paper, we extend the hierarchical Bayesian model applied in Eadie & Jurić (2019) to study the mass distribution of the Milky Way halo; the new model allows for the use of all available 6D phase-space measurements. We use kinematic data of halo stars out to $142~{\rm kpc}$, obtained from the H3 Survey and $\textit{Gaia}$ EDR3, to infer the mass of the Galaxy. Inference is carried out with the No-U-Turn sampler, a fast and scalable extension of Hamiltonian Monte Carlo. We report a median mass enclosed within $100~{\rm kpc}$ of $\rm M(<100 \; kpc) = 0.69_{-0.04}^{+0.05} \times 10^{12} \; M_\odot$ (68% Bayesian credible interval), or a virial mass of $\rm M_{200} = M(<216.2_{-7.5}^{+7.5} \; kpc) = 1.08_{-0.11}^{+0.12} \times 10^{12} \; M_\odot$, in good agreement with other recent estimates. We analyze our results using posterior predictive checks and find limitations in the model's ability to describe the data. In particular, we find sensitivity with respect to substructure in the halo, which limits the precision of our mass estimates to $\sim 15\%$.

preprint2020arXiv

A Diffuse Metal-Poor Component of the Sagittarius Stream Revealed by the H3 Survey

The tidal disruption of the Sagittarius dwarf galaxy has generated a spectacular stream of stars wrapping around the entire Galaxy. We use data from $Gaia$ and the H3 Stellar Spectroscopic Survey to identify 823 high-quality Sagittarius members based on their angular momenta. The H3 Survey is largely unbiased in metallicity, and so our sample of Sagittarius members is similarly unbiased. Stream stars span a wide range in [Fe/H] from $-0.2$ to $\approx -3.0$, with a mean overall metallicity of $\langle$[Fe/H]$\rangle=-0.99$. We identify a strong metallicity-dependence to the kinematics of the stream members. At [Fe/H]$\gt -0.8$ nearly all members belong to the well-known cold ($σ_v \lt 20$ km/s) leading and trailing arms. At intermediate metallicities ($-1.9 \lt$[Fe/H]$\lt -0.8$) a significant population (24$\%$) emerges of stars that are kinematically offset from the cold arms. These stars also appear to have hotter kinematics. At the lowest metallicities ([Fe/H]$\lesssim-2$), the majority of stars (69$\%$) belong to this kinematically-offset diffuse population. Comparison to simulations suggests that the diffuse component was stripped from the Sagittarius progenitor at earlier epochs, and therefore resided at larger radius on average, compared to the colder metal-rich component. We speculate that this kinematically diffuse, low metallicity, population is the stellar halo of the Sagittarius progenitor system.

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

Ancient Very Metal-Poor Stars Associated With the Galactic Disk in the H3 Survey

Ancient, very metal-poor stars offer a window into the earliest epochs of galaxy formation and assembly. We combine data from the H3 Spectroscopic Survey and Gaia to measure metallicities, abundances of $α$ elements, stellar ages, and orbital properties of a sample of 482 very metal-poor (VMP; [Fe/H]$<-2$) stars in order to constrain their origins. This sample is confined to $1\lesssim |Z| \lesssim3$ kpc from the Galactic plane. We find that >70% of VMP stars near the disk are on prograde orbits and this fraction increases toward lower metallicities. This result unexpected if metal-poor stars are predominantly accreted from many small systems with no preferred orientation, as such a scenario would imply a mostly isotropic distribution. Furthermore, we find there is some evidence for higher fractions of prograde orbits amongst stars with lower [$α$/Fe]. Isochrone-based ages for main sequence turn-off stars reveal that these VMP stars are uniformly old ($\approx12$ Gyr) irrespective of the $α$ abundance and metallicity, suggesting that the metal-poor population was not born from the same well-mixed gas disk. We speculate that the VMP population has a heterogeneous origin, including both in-situ formation in the ancient disk and accretion from a satellite with the same direction of rotation as the ancient disk at early times. Our precisely measured ages for these VMP stars on prograde orbits show that the Galaxy has had a relatively quiescent merging history over most of cosmic time, and implies the angular momentum alignment of the Galaxy has been in place for at least 12 Gyr.

preprint2020arXiv

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

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

preprint2020arXiv

High-resolution spectroscopy of the GD-1 stellar stream localizes the perturber near the orbital plane of Sagittarius

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $σ_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [α/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

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

MINESweeper: Spectrophotometric Modeling of Stars in the Gaia Era

We present MINESweeper, a tool to measure stellar parameters by jointly fitting observed spectra and broadband photometry to model isochrones and spectral libraries. This approach enables the measurement of spectrophotometric distances, in addition to stellar parameters such as Teff, log(g), [Fe/H], [a/Fe], and radial velocity. MINESweeper employs a Bayesian framework and can easily incorporate a variety of priors, including Gaia parallaxes. Mock data are fit in order to demonstrate how the precision of derived parameters depends on evolutionary phase and SNR. We then fit a selection of data in order to validate the model outputs. Fits to a variety of benchmark stars including Procyon, Arcturus, and the Sun result in derived stellar parameters that are in good agreement with the literature. We then fit combined spectra and photometry of stars in the open and globular clusters M92, M13, M3, M107, M71, and M67. Derived distances, [Fe/H], [a/Fe], and log(g)-Teff, relations are in overall good agreement with literature values, although there are trends between metallicity and log(g), within clusters that point to systematic uncertainties at the ~0.1 dex level. Finally, we fit a large sample of stars from the H3 Spectroscopic Survey in which high quality Gaia parallaxes are also available. These stars are fit without the Gaia parallaxes so that the geometric parallaxes can serve as an independent test of the spectrophotometric distances. Comparison between the two reveals good agreement within their formal uncertainties after accounting for the Gaia zero point uncertainties.

preprint2020arXiv

Powderday: Dust Radiative Transfer for Galaxy Simulations

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

preprint2020arXiv

Rapid Reionization by the Oligarchs: The Case for Massive, UV-Bright, Star-Forming Galaxies with High Escape Fractions

The protagonists of cosmic reionization remain elusive. Faint star-forming galaxies are leading candidates because they are numerous and may have significant ionizing photon escape fractions ($f_{esc}$). Here we update this picture via an empirical model that successfully predicts latest observations (e.g., the drop in star-formation density at z>8). We generate an ionizing spectrum for each galaxy in our model and constrain $f_{esc}$ using latest measurements of the reionization timeline (e.g., Ly$α$ damping of quasars and galaxies at z>7). Assuming a constant $f_{esc}$, we find $M_{UV}$<-13.5 galaxies need $f_{esc}=0.21^{+0.06}_{-0.04}$ to complete reionization. The inferred IGM neutral fraction is [0.9, 0.5, 0.1] at z=[8.2, 6.8, 6.2]$\pm$0.2, i.e., the bulk of reionization transpires in 300 Myrs. Inspired by the emergent sample of Lyman Continuum (LyC) leakers that overwhelmingly displays higher-than-average star-formation surface density ($Σ$), we propose a model relating $f_{esc}$ to $Σ$ and find $f_{esc}\proptoΣ^{0.4\pm0.1}$. Since $Σ$ falls by ~2.5 dex between z=8 and z=0, our model explains the humble upper limits on $f_{esc}$ at lower redshifts and its required evolution to ~0.2 at z>6. Within this model, strikingly, <5% of galaxies with $M_{UV}$<-18 (the `oligarchs') account for >80% of the reionization budget. In fact, faint sources ($M_{UV}$>-16) must be relegated to a limited role to ensure high neutral fractions at z=7-8. Shallow faint-end slopes of the UV luminosity function ($α$>-2) and/or $f_{esc}$ distributions skewed toward bright galaxies produce the required late and rapid reionization. We predict LyC leakers like COLA1 (z=6.6, $f_{esc}$~30%, $M_{UV}$=-21.5) become increasingly common towards z~6 and that the drivers of reionization do not lie hidden across the faint-end of the luminosity function, but are already known to us. (abridged)

preprint2020arXiv

SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry

We present SPECULATOR - a fast, accurate, and flexible framework for emulating stellar population synthesis (SPS) models for predicting galaxy spectra and photometry. For emulating spectra, we use principal component analysis to construct a set of basis functions, and neural networks to learn the basis coefficients as a function of the SPS model parameters. For photometry, we parameterize the magnitudes (for the filters of interest) as a function of SPS parameters by a neural network. The resulting emulators are able to predict spectra and photometry under both simple and complicated SPS model parameterizations to percent-level accuracy, giving a factor of $10^3$-$10^4$ speed up over direct SPS computation. They have readily-computable derivatives, making them amenable to gradient-based inference and optimization methods. The emulators are also straightforward to call from a GPU, giving an additional order-of-magnitude speed-up. Rapid SPS computations delivered by emulation offers a massive reduction in the computational resources required to infer the physical properties of galaxies from observed spectra or photometry and simulate galaxy populations under SPS models, whilst maintaining the accuracy required for a range of applications.

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

A Lower Limit on the Mass of Our Galaxy from the H3 Survey

The timing argument provides a lower limit on the mass of the Milky Way. We find, using a sample of 32 stars at $R > 60$ kpc drawn from the H3 Spectroscopic Survey and mock catalogs created from published numerical simulations, that M$_{200} > 0.91\times 10^{12}$ M$_\odot$ with 90% confidence. We recommend using this limit to refine the allowed prior mass range in more complex and sophisticated statistical treatments of Milky Way dynamics. The use of such a prior would have significantly reduced many previously published uncertainty ranges. Our analysis suggests that the most likely value of M$_{200}$ is $\sim 1.4 \times 10^{12}$ M$_\odot$, but establishing this as the Milky Way mass requires a larger sample of outer halo stars and a more complete analysis of the inner halo stars in H3. The imminent growth in the sample of outer halo stars due to ongoing and planned surveys will make this possible.

preprint2019arXiv

Combined Effects of Rotation and Age Spreads on Extended Main Sequence Turn Offs

The extended main sequence turn offs (eMSTOs) of several young to intermediate age clusters are examined in the Magellanic Clouds and the Milky Way. We explore the effects of extended star formation (eSF) and a range of stellar rotation rates on the behavior of the color-magnitude diagram (CMD), paying particular attention to the MSTO. We create synthetic stellar populations based on MESA stellar models to simulate observed Hubble Space Telescope and Gaia star cluster data. We model the effect of rotation as a non-parametric distribution, allowing for maximum flexibility. In our models the slow rotators comprise the blueward, and fast rotators the redward portion of the eMSTO. We simulate data under three scenarios: non-rotating eSF, a range of rotation rates with a single age, and a combination of age and rotation effects. We find that two of the five clusters (the youngest and oldest) favor an age spread, but these also achieve the overall worst fits of all clusters. The other three clusters show comparable statistical evidence between rotation and an age spread. In all five cases, a rotation rate distribution alone is capable of qualitatively matching the observed eMSTO structure. In future work, we aim to compare our predicted Vsin(i) with observations in order to better constrain the physics related to stellar rotation.

preprint2019arXiv

Resolving the Metallicity Distribution of the Stellar Halo with the H3 Survey

The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars, and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically-selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, [Fe/H]=-1.2, and there is no discernable metallicity gradient over the range $6<R_{\rm gal}<100$ kpc. However, the halo metallicity distribution is highly structured including distinct metal-rich and metal-poor components at $R_{\rm gal}<10$ kpc and $R_{\rm gal}>30$ kpc, respectively. Metal-poor stars with [Fe/H]$<-2$ are a small population of the halo at all distances and orbital categories. We associate the "in-situ" stellar halo with stars displaying thick-disk chemistry on halo-like orbits; such stars are confined to $|z|<10$ kpc. The majority of the stellar halo is resolved into discrete features in orbital-chemical space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. (ABRIDGED)

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

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

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

preprint2016arXiv

Accelerated Fitting of Stellar Spectra

Stellar spectra are often modeled and fit by interpolating within a rectilinear grid of synthetic spectra to derive the stars' labels: stellar parameters and elemental abundances. However, the number of synthetic spectra needed for a rectilinear grid grows exponentially with the label space dimensions, precluding the simultaneous and self-consistent fitting of more than a few elemental abundances. Shortcuts such as fitting subsets of labels separately can introduce unknown systematics and do not produce correct error covariances in the derived labels. In this paper we present a new approach -- CHAT (Convex Hull Adaptive Tessellation) -- which includes several new ideas for inexpensively generating a sufficient stellar synthetic library, using linear algebra and the concept of an adaptive, data-driven grid. A convex hull approximates the region where the data lie in the label space. A variety of tests with mock datasets demonstrate that CHAT can reduce the number of required synthetic model calculations by three orders of magnitude in an 8D label space. The reduction will be even larger for higher-dimensional label spaces. In CHAT the computational effort increases only linearly with the number of labels that are fit simultaneously. Around each of these grid points in label space an approximate synthetic spectrum can be generated through linear expansion using a set of "gradient spectra" that represent flux derivatives at every wavelength point with respect to all labels. These techniques provide new opportunities to fit the full stellar spectra from large surveys with 15-30 labels simultaneously.

preprint2016arXiv

Constructing Polynomial Spectral Models for Stars

Stellar spectra depend on the stellar parameters and on dozens of photospheric elemental abundances. Simultaneous fitting of these $\mathcal{N}\sim 10-40$ model labels to observed spectra has been deemed unfeasible, because the number of ab initio spectral model grid calculations scales exponentially with $\mathcal{N}$. We suggest instead the construction of a polynomial spectral model (PSM) of order $\mathcal{O}$ for the model flux at each wavelength. Building this approximation requires a minimum of only ${\mathcal{N}+\mathcal{O}\choose\mathcal{O}}$ calculations: e.g. a quadratic spectral model ($\mathcal{O}=2$) to fit $\mathcal{N}=20$ labels simultaneously, can be constructed from as few as $231$ ab initio spectral model calculations; in practice, a somewhat larger number ($\sim 300-1000$) of randomly chosen models lead to a better performing PSM. Such a PSM can be a good approximation only over a portion of label space, which will vary case by case. Yet, taking the APOGEE survey as an example, a single quadratic PSM provides a remarkably good approximation to the exact ab initio spectral models across much of this survey: for random labels within that survey the PSM approximates the flux to within $10^{-3}$, and recovers the abundances to within $\sim 0.02$ dex rms of the exact models. This enormous speed-up enables the simultaneous many-label fitting of spectra with computationally expensive ab initio models for stellar spectra, such as non-LTE models. A PSM also enables the simultaneous fitting of observational parameters, such as the spectrum's continuum or line-spread function.

preprint2016arXiv

Exploring Systematic Effects in the Relation Between Stellar Mass, Gas Phase Metallicity, and Star Formation Rate

There is evidence that the well-established mass-metallicity relation in galaxies is correlated with a third parameter: star formation rate (SFR). The strength of this correlation may be used to disentangle the relative importance of different physical processes (e.g., infall of pristine gas, metal-enriched outflows) in governing chemical evolution. However, all three parameters are susceptible to biases that might affect the observed strength of the relation between them. We analyze possible sources of systematic error, including sample bias, application of signal-to-noise ratio cuts on emission lines, choice of metallicity calibration, uncertainty in stellar mass determination, aperture effects, and dust. We present the first analysis of the relation between stellar mass, gas phase metallicity, and SFR using strong line abundance diagnostics from Dopita et al. (2013) for ~130,000 star-forming galaxies in the Sloan Digital Sky Survey and provide a detailed comparison of these diagnostics in an appendix. Using these abundance diagnostics yields a 30-55% weaker anti-correlation between metallicity and SFR at fixed stellar mass than that reported by Mannucci et al. (2010). We find that, for all abundance diagnostics, the anti-correlation with SFR is stronger for the relatively few galaxies whose current SFRs are elevated above their past average SFRs. This is also true for the new abundance diagnostic of Dopita et al. (2016), which gives anti-correlation between metallicity and SFR only in the high specific star formation rate (sSFR) regime, in contrast to the recent results of Kashino et al. (2016). The poorly constrained strength of the relation between stellar mass, metallicity, and SFR must be carefully accounted for in theoretical studies of chemical evolution.

preprint2016arXiv

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

This is the first of a series of papers presenting the Modules for Experiments in Stellar Astrophysics (MESA) Isochrones and Stellar Tracks (MIST) project, a new comprehensive set of stellar evolutionary tracks and isochrones computed using MESA, a state-of-the-art open-source 1D stellar evolution package. In this work, we present models with solar-scaled abundance ratios covering a wide range of ages ($5 \leq \rm \log(Age)\;[yr] \leq 10.3$), masses ($0.1 \leq M/M_{\odot} \leq 300$), and metallicities ($-2.0 \leq \rm [Z/H] \leq 0.5$). The models are self-consistently and continuously evolved from the pre-main sequence to the end of hydrogen burning, the white dwarf cooling sequence, or the end of carbon burning, depending on the initial mass. We also provide a grid of models evolved from the pre-main sequence to the end of core helium burning for $-4.0 \leq \rm [Z/H] < -2.0$. We showcase extensive comparisons with observational constraints as well as with some of the most widely used existing models in the literature. The evolutionary tracks and isochrones can be downloaded from the project website at http://waps.cfa.harvard.edu/MIST/.

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

Redshift Measurement and Spectral Classification for eBOSS Galaxies with the Redmonster Software

We describe the redmonster automated redshift measurement and spectral classification software designed for the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey IV (SDSS-IV). We describe the algorithms, the template standard and requirements, and the newly developed galaxy templates to be used on eBOSS spectra. We present results from testing on early data from eBOSS, where we have found a 90.5% automated redshift and spectral classification success rate for the luminous red galaxy sample (redshifts 0.6$\lesssim z \lesssim$1.0). The redmonster performance meets the eBOSS cosmology requirements for redshift classification and catastrophic failures, and represents a significant improvement over the previous pipeline. We describe the empirical processes used to determine the optimum number of additive polynomial terms in our models and an acceptable $Δχ_r^2$ threshold for declaring statistical confidence. Statistical errors on redshift measurement due to photon shot noise are assessed, and we find typical values of a few tens of km s$^{-1}$. An investigation of redshift differences in repeat observations scaled by error estimates yields a distribution with a Gaussian mean and standard deviation of $μ\sim$0.01 and $σ\sim$0.65, respectively, suggesting the reported statistical redshift uncertainties are over-estimated by $\sim$54%. We assess the effects of object magnitude, signal-to-noise ratio, fiber number, and fiber head location on the pipeline's redshift success rate. Finally, we describe directions of ongoing development.

preprint2016arXiv

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

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

preprint2015arXiv

APOGEE chemical tagging constraint on the maximum star cluster mass in the $α$-enhanced Galactic disk

Stars born from the same molecular cloud should be nearly homogeneous in their element abundances. The concept of chemical tagging is to identify members of disrupted clusters by their clustering in element abundance space. Chemical tagging requires large samples of stars with precise abundances for many individual elements. With uncertainties of $σ_{[X/{\rm Fe}]}$ and $σ_{\rm [Fe/H]} \simeq 0.05$ for 10 elements measured for $>10^4$ stars, the APOGEE DR12 spectra may be the first well-suited data set to put this idea into practice. We find that even APOGEE data offer only $\sim 500$ independent volume elements in the 10-dimensional abundance space, when we focus on the $α$-enhanced Galactic disk. We develop and apply a new algorithm to search for chemically homogeneous sets of stars against a dominant background. By injecting star clusters into the APOGEE data set we show that chemically homogeneous clusters with masses $\gtrsim 3 \times 10^7 \, {\rm M}_\odot$ would be easily detectable and yet no such signal is seen in the data. By generalizing this approach, we put a first abundance-based constraint on the cluster mass function for the old disk stars in the Milky Way.

preprint2015arXiv

Beacons In the Dark: Using Novae and Supernovae to Detect Dwarf Galaxies in the Local Universe

We propose that luminous transients, including novae and supernovae, can be used to detect the faintest galaxies in the universe. Beyond a few Mpc, dwarf galaxies with stellar masses $<10^6 M_{\odot}$ will likely be too faint and/or too low in surface brightness to be directly detected in upcoming large area ground-based photometric surveys. However, single epoch LSST photometry will be able to detect novae to distances of $\sim30$ Mpc and SNe to Gpc-scale distances. Depending on the form of the stellar mass-halo mass relation and the underlying star formation histories of low mass dwarfs, the expected nova rates will be a few to $\sim100$ yr$^{-1}$ and the expected SN rates (including both type Ia and core-collapse) will be $\sim10^2-10^4$ within the observable ($4π$ sr) volume. The transient rate associated with intrahalo stars will be comparably large, but these transients will be located close to bright galaxies, in contrast to the dwarfs, which should trace the underlying large scale structure of the cosmic web. Aggressive follow-up of hostless transients has the potential to uncover the predicted enormous population of low mass field dwarf galaxies.

preprint2015arXiv

Circumstellar Dust Around AGB Stars and Implications for Infrared Emission from Galaxies

Stellar population synthesis (SPS) models are used to infer many galactic properties including star formation histories, metallicities, and stellar and dust masses. However, most SPS models neglect the effect of circumstellar dust shells around evolved stars and it is unclear to what extent they impact the analysis of SEDs. To overcome this shortcoming we have created a new set of circumstellar dust models, using the radiative transfer code DUSTY Ivezic et al. 1999, for asymptotic giant branch (AGB) stars and incorporated them into the Flexible Stellar Population Synthesis code. The circumstellar dust models provide a good fit to individual AGB stars as well as the IR color-magnitude diagrams of the Large and Small Magellanic Clouds. IR luminosity functions from the Large and Small Magellanic Clouds are not well-fit by the 2008 Padova isochrones when coupled to our circumstellar dust models, and so we adjusted the lifetimes of AGB stars in the models to provide a match to the data. We show, in agreement with previous work, that circumstellar dust from AGB stars can make a significant contribution to the IR ($\gtrsim4μm$) emission from galaxies that contain relatively little diffuse dust, including low-metallicity and/or non-star forming galaxies. Our models provide a good fit to the mid-IR spectra of early-type galaxies. Circumstellar dust around AGB stars appears to have a small effect on the IR SEDs of metal-rich star-forming galaxies (i.e., when A$_{\rm V}$ $\gtrsim$~0.1). Stellar population models that include circumstellar dust will be needed to accurately interpret data from the James Webb Space Telescope (JWST) and other IR facilities.

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

From Cosmic Birth to Living Earths: The Future of UVOIR Space Astronomy

For the first time in history, humans have reached the point where it is possible to construct a revolutionary space-based observatory that has the capability to find dozens of Earth-like worlds, and possibly some with signs of life. This same telescope, designed as a long-lived facility, would also produce transformational scientific advances in every area of astronomy and astrophysics from black hole physics to galaxy formation, from star and planet formation to the origins of the Solar System. The Association of Universities for Research in Astronomy (AURA) commissioned a study on a next-generation UVOIR space observatory with the highest possible scientific impact in the era following JWST. This community-based study focuses on the future space-based options for UV and optical astronomy that significantly advance our understanding of the origin and evolution of the cosmos and the life within it. The committee concludes that a space telescope equipped with a 12-meter class primary mirror can find and characterize dozens of Earth-like planets and make fundamental advances across nearly all fields of astrophysics. The concept is called the High Definition Space Telescope (HDST). The telescope would be located at the Sun-Earth L2 point and would cover a spectral range that, at a minimum, runs from 0.1 to 2 microns. Unlike JWST, HDST will not need to operate at cryogenic temperatures. HDST can be made to be serviceable on orbit but does not require servicing to complete its primary scientific objectives. We present the scientific and technical requirements for HDST and show that it could allow us to determine whether or not life is common outside the Solar System. We do not propose a specific design for such a telescope, but show that designing, building and funding such a facility is feasible beginning in the next decade - if the necessary strategic investments in technology begin now.

preprint2015arXiv

Hiding in plain sight: record-breaking compact stellar systems in the Sloan Digital Sky Survey

Motivated by the recent, serendipitous discovery of the densest known galaxy, M60-UCD1, we present two initial findings from a follow-up search, using the Sloan Digital Sky Survey, Subaru/Suprime-Cam and Hubble Space Telescope imaging, and SOuthern Astrophysical Research (SOAR)/Goodman spectroscopy. The first object discovered, M59-UCD3, has a similar size to M60-UCD1 (half-light radius of r_h ~ 20 pc) but is 40% more luminous (M_V ~ -14.6), making it the new densest-known galaxy. The second, M85-HCC1, has a size like a typical globular cluster (GC; r_h ~ 1.8 pc) but is much more luminous (M_V ~ -12.5). This hypercompact cluster is by far the densest confirmed free-floating stellar system, and is equivalent to the densest known nuclear star clusters. From spectroscopy, we find that both objects are relatively young (~9 Gyr and ~3 Gyr, respectively), with metal-abundances that resemble those of galaxy centers. Their host galaxies show clear signs of large-scale disturbances, and we conclude that these dense objects are the remnant nuclei of recently accreted galaxies. M59-UCD3 is an ideal target for follow-up with high-resolution imaging and spectroscopy to search for an overweight central supermassive black hole as was discovered in M60-UCD1. These findings also emphasize the potential value of ultra-compact dwarfs and massive GCs as tracers of the assembly histories of galaxies.

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

Prospects for chemically tagging stars in the Galaxy

It is now well-established that the elemental abundance patterns of stars holds key clues not only to their formation but also to the assembly histories of galaxies. One of the most exciting possibilities is the use of stellar abundance patterns as "chemical tags" to identify stars that were born in the same molecular cloud. In this paper we assess the prospects of chemical tagging as a function of several key underlying parameters. We show that in the fiducial case of $10^4$ distinct cells in chemical space and $10^5-10^6$ stars in the survey, one can expect to detect $\sim 10^2-10^3$ groups that are $\ge 5σ$ overdensities in the chemical space. However, we find that even very large overdensities in chemical space do not guarantee that the overdensity is due to a single set of stars from a common birth cloud. In fact, for our fiducial model parameters, the typical $5σ$ overdensity is comprised of stars from a wide range of clusters with the most dominant cluster contributing only 25% of the stars. The most important factors limiting the identification of disrupted clusters via chemical tagging are the number of chemical cells in the chemical space and the survey sampling rate of the underlying stellar population. Both of these factors can be improved through strategic observational plans. While recovering individual clusters through chemical tagging may prove challenging, we show, in agreement with previous work, that different CMFs imprint different degrees of clumpiness in chemical space. These differences provide the opportunity to statistically reconstruct the slope and high mass cutoff of CMF and its evolution through cosmic time.

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

Stellar Populations of Barred Quiescent Galaxies

Selecting centrally quiescent galaxies from the Sloan Digital Sky Survey (SDSS) to create high signal-to-noise (>100) stacked spectra with minimal emission line contamination, we accurately and precisely model the central stellar populations of barred and unbarred quiescent disk galaxies. By splitting our sample by redshift, we can use the fixed size of the SDSS fiber to model the stellar populations at different radii within galaxies. At 0.02<z<0.04, the SDSS fiber radius corresponds to ~1 kpc, which is the typical half-light radii of both classical bulges and disky pseudobulges. Assuming that the SDSS fiber primarily covers the bulges at these redshifts, our analysis shows that there are no significant differences in the stellar populations, i.e., stellar age, [Fe/H], [Mg/Fe], and [N/Fe], of the bulges of barred vs. unbarred quiescent disk galaxies. Modeling the stellar populations at different redshift intervals from z=0.020 to z=0.085 at fixed stellar masses produces an estimate of the stellar population gradients out to about half the typical effective radius of our sample, assuming null evolution over this ~1 Gyr epoch. We find that there are no noticeable differences in the slopes of the azimuthally averaged gradients of barred vs. unbarred quiescent disk galaxies. These results suggest that bars are not a strong influence on the chemical evolution of quiescent disk galaxies.

preprint2015arXiv

The IMF-sensitive 1.14-micron Na I doublet in early-type galaxies

We present J-band spectroscopy of passive galaxies focusing on the Na I doublet at 1.14 μm. Like the Na I 0.82 μm doublet, this feature is strong in low-mass stars and hence may provide a useful probe of the initial mass function (IMF). From high signal-to-noise composite spectra, we find that Na I 1.14 μm increases steeply with increasing velocity dispersion, σ, and for the most massive galaxies (σ> 300 km/s) is much stronger than predicted from synthetic spectra with Milky-Way-like IMFs and solar abundances. Reproducing Na I 1.14 μm at high σ likely requires either a very high [Na/H], or a bottom-heavy IMF, or a combination of both. Using the Na D line to break the degeneracy between IMF and abundance, we infer [Na/H] $\approx$ +0.5 and a steep IMF (single-slope-equivalent x $\approx$ 3.2, where x = 2.35 for Salpeter), for the high-σgalaxies. At lower mass (σ = 50-100 km/s), the line strengths are compatible with MW-like IMFs and near-solar [Na/H]. We highlight two galaxies in our sample where strong gravitational lensing masses favour MW-like IMFs. Like the high-σ sample on average, these galaxies have strong Na I 1.14 μm; taken in isolation their sodium indices imply bottom-heavy IMFs which are hard to reconcile with the lensing masses. An alternative full-spectrum-fitting approach, applied to the high-σsample, recovers an IMF less heavy than Salpeter, but under-predicts the Na I 1.14 μm line at the 5σ level. We conclude that current models struggle to reproduce this feature in the most massive galaxies without breaking other constraints, and caution against over-reliance on the sodium lines in spectroscopic IMF studies.

preprint2015arXiv

The Local Group as a time machine: studying the high-redshift Universe with nearby galaxies

We infer the UV luminosities of Local Group galaxies at early cosmic times ($z \sim 2$ and $z \sim 7$) by combining stellar population synthesis modeling with star formation histories derived from deep color-magnitude diagrams constructed from Hubble Space Telescope (HST) observations. Our analysis provides a basis for understanding high-$z$ galaxies - including those that may be unobservable even with the James Webb Space Telescope (JWST) - in the context of familiar, well-studied objects in the very low-$z$ Universe. We find that, at the epoch of reionization, all Local Group dwarfs were less luminous than the faintest galaxies detectable in deep HST observations of blank fields. We predict that JWST will observe $z \sim 7$ progenitors of galaxies similar to the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may already be observing such galaxies, highlighting the power of gravitational lensing. Consensus reionization models require an extrapolation of the observed blank-field luminosity function at $z \approx 7$ by at least two orders of magnitude in order to maintain reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf spheroidal galaxies to be contributors to the ionizing background at $z \sim 7$. Combined with numerical simulations, our results argue for a break in the UV luminosity function from a faint-end slope of $α\sim -2$ at $M_{\rm UV} < -13$ to $α\sim -1.2$ at lower luminosities. Applied to photometric samples at lower redshifts, our analysis suggests that HST observations in lensing fields at $z \sim 2$ are capable of probing galaxies with luminosities comparable to the expected progenitor of Fornax.

preprint2015arXiv

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

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

preprint2015arXiv

The SINFONI Nearby Elliptical Lens Locator Survey: Discovery of two new low-redshift strong lenses and implications for the initial mass function in giant early-type galaxies

We present results from a blind survey to identify strong gravitational lenses among the population of low-redshift early-type galaxies. The SINFONI Nearby Elliptical Lens Locator Survey (SNELLS) uses integral-field infrared spectroscopy to search for lensed emission line sources behind massive lens candidates at $z$<0.055. From 27 galaxies observed, we have recovered one previously-known lens (ESO325-G004) at $z$=0.034, and discovered two new systems, at $z$=0.031 and $z$=0.052. All three lens galaxies have high velocity dispersions (σ>300 km/s) and α-element abundances ([Mg/Fe]>0.3). From the lensing configurations we derive total J-band mass-to-light ratios of 1.8$\pm$0.1, 2.1$\pm$0.1 and 1.9$\pm$0.2 within the $\sim$2 kpc Einstein radius. Correcting for estimated dark-matter contributions, and comparing to stellar population models with a Milky Way (Kroupa) initial mass function (IMF), we determine the "mass excess factor", α. Assuming the lens galaxies have "old" stellar populations (10$\pm$1 Gyr), the average IMF mass factor is $\langleα\rangle$=1.10$\pm$0.08$\pm$0.10, where the first error is random and the second is systematic. If we instead fit the stellar populations from 6dF optical survey spectra, all three galaxies are consistent with being old, but the age errors are 3-4 Gyr, due to limited signal-to-noise. The IMF constraints are therefore looser in this case, with $\langleα\rangle$ = $1.23^{+0.16}_{-0.13}\pm{0.10}$. Our results are thus consistent with a Kroupa IMF (α=1.00) on average, and strongly reject very heavy IMFs with α>2. A Salpeter IMF (α=1.55) is inconsistent at the 3.5$σ$ level if the galaxies are old, but cannot be excluded using age constraints derived from the currently-available optical spectra.

preprint2014arXiv

Comparing Simple Quasar Demographics Models

This paper explores several simple model variations for the connections among quasars, galaxies, and dark matter halos for redshifts 1 < z < 6. A key component of these models is that we enforce a self-consistent black hole (BH) history by tracking both BH mass and BH growth rate at all redshifts. We connect objects across redshift with a simple constant-number-density procedure, and choose a fiducial model with a relationship between BH and galaxy growth rates that is linear and evolves in a simple way with redshift. Within this fiducial model, we find the quasar luminosity function (QLF) by calculating an "intrinsic" luminosity based on either the BH mass or BH growth rate, and then choosing a model of quasar variability with either a lognormal or truncated power-law distribution of instantaneous luminosities. This gives four model variations, which we fit to the observed QLF at each redshift. With the best-fit models in hand, we undertake a detailed comparison of the four fiducial models, and explore changes to our fiducial model of the BH-galaxy relationship. Each model variation can successfully fit the observed QLF, the shape of which is generally set by the "intrinsic" luminosity at the faint end and by the scatter due to variability at the bright end. We focus on accounting for the reasons that physically different models can make such similar predictions, and on identifying what observational data or physical arguments are most essential in breaking the degeneracies among models.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

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 evolution near groups and clusters: ejected satellites and the spatial extent of environmental quenching

Galaxies that are several virial radii beyond groups/clusters show preferentially quiescent star formation rates. Using a galaxy group/cluster catalog from the Sloan Digital Sky Survey, together with a cosmological N-body simulation, we examine the origin of this environmental quenching beyond the virial radius. Accounting for the clustering of groups/clusters, we show that central galaxies show enhanced SFR quenching out to 2.5 virial radii beyond groups/clusters, and we demonstrate that this extended environmental enhancement can be explained simply by 'ejected' satellite galaxies that orbit beyond their host halo's virial radius. We show that ejected satellites typically orbit for several Gyr beyond the virial radius before falling back in, and thus they compose up to 40% of all central galaxies near groups/clusters. We show that a model in which ejected satellites experience the same SFR quenching as satellites within a host halo can explain essentially all environmental dependence of galaxy quenching. Furthermore, ejected satellites (continue to) lose significant halo mass, an effect that is potentially observable via gravitational lensing. The SFRs/colors and stellar-to-halo masses of ejected satellites highlight the importance of environmental history and present challenges to models of galaxy occupation that ignore such history.

preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

preprint2014arXiv

Parametrising Star Formation Histories

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

preprint2014arXiv

Reversal of Fortune: Increased Star Formation Efficiencies in the Early Histories of Dwarf Galaxies?

On dwarf galaxy scales, the different shapes of the galaxy stellar mass function and the dark halo mass function require a star-formation efficiency (SFE) in these systems that is currently more than 1 dex lower than that of Milky Way-size halos. Here, we argue that this trend may actually be reversed at high redshift. Specifically, by combining the resolved star-formation histories of nearby isolated dwarfs with the simulated mass-growth rates of dark matter halos, we show that the assembly of these systems occurs in two phases: (1) an early, fast halo accretion phase with a rapidly deepening potential well, characterized by a high SFE; and (2) a late slow halo accretion phase where, perhaps as a consequence of reionization, the SFE is low. Nearby dwarfs have more old stars than predicted by assuming a constant or decreasing SFE with redshift, a behavior that appears to deviate qualitatively from the trends seen amongst more massive systems. Taken at face value, the data suggest that, at sufficiently early epochs, dwarf galaxy halos above the atomic cooling mass limit can be among the most efficient sites of star formation in the universe.

preprint2014arXiv

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

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

preprint2014arXiv

Stellar Models of Multiple Populations in Globular Clusters. I. The Main Sequence of NGC 6752

We present stellar atmosphere and evolution models of main sequence stars in two stellar populations of the Galactic globular cluster NGC 6752. These populations represent the two extremes of light-element abundance variations in the cluster. NGC 6752 is a benchmark cluster in the study of multiple stellar populations because of the rich array of spectroscopic abundances and panchromatic Hubble Space Telescope photometry. The spectroscopic abundances are used to compute stellar atmosphere and evolution models. The synthetic spectra for the two populations show significant differences in the ultraviolet and, for the coolest temperatures, in the near-infrared. The stellar evolution models exhibit insignificant differences in the H-R diagram except on the lower main sequence. The appearance of multiple sequences in the colour-magnitude diagrams (CMDs) of NGC 6752 is almost exclusively due to spectral effects caused by the abundance variations. The models reproduce the observed splitting and/or broadening of sequences in a range of CMDs. The ultraviolet CMDs are sensitive to variations in carbon, nitrogen, and oxygen but the models are not reliable enough to directly estimate abundance variations from photometry. On the other hand, the widening of the lower main sequence in the near-infrared CMD, driven by oxygen-variation via the water molecule, is well-described by the models and can be used to estimate the range of oxygen present in a cluster from photometry. We confirm that it is possible to use multiband photometry to estimate helium variations among the different populations, with the caveat that the estimated amount of helium-enhancement is model-dependent.

preprint2014arXiv

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

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

preprint2014arXiv

The Baryon Cycle of Dwarf Galaxies: Dark, Bursty, Gas-Rich Polluters

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

preprint2014arXiv

The SLUGGS survey: Globular cluster stellar population trends from weak absorption lines in stacked spectra

As part of the SLUGGS survey, we stack 1137 Keck DEIMOS spectra of globular clusters from 10 galaxies to study their stellar populations in detail. The stacked spectra have median signal to noise ratios of $\sim 90$ Å$^{-1}$. Besides the calcium triplet, we study weaker sodium, magnesium, titanium and iron lines as well as the H$α$ and higher order Paschen hydrogen lines. In general, the stacked spectra are consistent with old ages and a Milky Way-like initial mass function. However, we see different metal line index strengths at fixed colour and magnitude, and differences in the calcium triplet--colour relation from galaxy to galaxy. We interpret this as strong evidence for variations in the globular cluster colour--metallicity relation between galaxies. Two possible explanations for the colour--metallicity relation variations are that the average ages of globular clusters vary from galaxy to galaxy or that the average abundances of light elements (i.e. He, C, N and O) differ between galaxies. Stacking spectra by magnitude, we see that the colours become redder and metal line indices stronger with brighter magnitudes. These trends are consistent with the previously reported `blue tilts' being mass--metallicity relations.

preprint2014arXiv

The Very Faint End of the UV Luminosity Function over Cosmic Time: Constraints from the Local Group Fossil Record

We present a new technique to estimate the evolution of the very faint end of the UV luminosity function (LF) out to $z\sim5$. Measured star formation histories (SFHs) from the fossil record of Local Group galaxies are used to reconstruct the LF down to M$_{UV}\sim-5$ at $z\sim5$ and M$_{UV}\sim-1.5$ at $z<1$. Such faint limits are well beyond the current observational limits and are likely to remain beyond the limits of next generation facilities. The reconstructed LFs, when combined with direct measurements of the LFs at higher luminosity, are well-fit by a standard Schechter function with no evidence of a break to the faintest limits probed by this technique. The derived faint end slope, $α$, steepens from $\approx-1.2$ at $z<1$ to $\approx-1.6$ at $4<z<5$. We test the effects of burstiness in the SFHs and find the recovered LFs to be only modestly affected. Incompleteness corrections for the faintest Local Group galaxies and the (unlikely) possibility of significant luminosity-dependent destruction of dwarf galaxies between high redshift and the present epoch are important uncertainties. These and other uncertainties can be mitigated with more detailed modeling and future observations. The reconstructed faint end LF from the fossil record can therefore be a powerful and complementary probe of the high redshift faint galaxies believed to play a key role in the reionization of the Universe.

preprint2013arXiv

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

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

preprint2013arXiv

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

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

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

preprint2013arXiv

Galaxy evolution in groups and clusters: satellite star formation histories and quenching timescales in a hierarchical Universe

Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFR) and lie on the red-sequence than central (field) galaxies. Using galaxy group/cluster catalogs from SDSS DR7, together with a cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching timescales of satellites of M_star > 5 x 10^9 M_sun at z=0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z~0.5, or ~5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs and quiescent fractions. With this, we constrain the importance and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at M_star < 10^10 M_sun. We then constrain satellite star formation histories, finding a 'delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2-4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of < 0.8 Gyr. These quenching timescales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in satellite quiescent fraction with halo mass arises simply because of satellites quenching in a lower mass group prior to infall (group preprocessing), which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact provides key physical insight into the subhalo abundance matching method.

preprint2013arXiv

Modeling the Panchromatic Spectral Energy Distributions of Galaxies

The spectral energy distributions (SEDs) of galaxies are shaped by nearly every physical property of the system, including the star formation history, metal content, abundance pattern, dust mass, grain size distribution, star-dust geometry, and interstellar radiation field. The principal goal of stellar population synthesis (SPS) is to extract these variables from observed SEDs. In this review I provide an overview of the SPS technique and discuss what can be reliably measured from galaxy SEDs. Topics include stellar masses, star formation rates and histories, metallicities and abundance patterns, dust properties, and the stellar initial mass function.

preprint2013arXiv

Star Clusters in M31. V. Evidence for Self-Enrichment in Old M31 Clusters from Integrated Spectroscopy

In the past decade, the notion that globular clusters (GCs) are composed of coeval stars with homogeneous initial chemical compositions has been challenged by growing evidence that they host an intricate stellar population mix, likely indicative of a complex history of star formation and chemical enrichment. Several models have been proposed to explain the existence of multiple stellar populations in GCs, but no single model provides a fully satisfactory match to existing data. Correlations between chemistry and global parameters such as cluster mass or luminosity are fundamental clues to the physics of GC formation. In this Letter, we present an analysis of the mean abundances of Fe, Mg, C, N, and Ca for 72 old GCs from the Andromeda galaxy. We show for the first time that there is a correlation between the masses of GCs and the mean stellar abundances of nitrogen, spanning almost two decades in mass. This result sheds new light on the formation of GCs, providing important constraints on their internal chemical evolution and mass loss history.

preprint2013arXiv

Stellar mass-gap as a probe of halo assembly history and concentration: youth hidden among old fossils

We investigate the use of the halo mass-gap statistic --- defined as the logarithmic difference in mass between the host halo and its most massive satellite subhalo --- as a probe of halo age and concentration. A cosmological N-body simulation is used to study N ~25, 000 group/cluster sized halos in the mass range 10^12.5 < M_halo/M_sol < 10^14.5. In agreement with previous work, we find that halo mass-gap is related to halo formation time and concentration. On average, older and more highly concentrated halos have larger halo mass-gaps, and this trend is stronger than the mass-concentration relation over a similar dynamic range. However, there is a large amount of scatter owing to the transitory nature of the satellite subhalo population, which limits the use of the halo mass-gap statistic on an object-by-object basis. For example, we find that 20% of very large halo mass-gap systems (akin to "fossil groups") are young, and have likely experienced a recent merger between a massive satellite subhalo and the central subhalo. We relate halo mass-gap to the observable stellar mass-gap via abundance matching. Using a galaxy group catalog constructed from the Sloan Digital Sky Survey Data Release 7, we find that the star formation and structural properties of galaxies at fixed mass show no trend with stellar mass-gap. This is despite a variation in halo-age of ~2.5 Gyr over ~1.2 dex in stellar mass-gap. Thus, we find no evidence to suggest that the halo formation history significantly affects galaxy properties.

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 Assembly of Milky Way-like Galaxies Since z~2.5

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

preprint2013arXiv

The Average Star Formation Histories of Galaxies in Dark Matter Halos from z=0-8

We present a robust method to constrain average galaxy star formation rates, star formation histories, and the intracluster light as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates, and cosmic star formation rates from z=0 to z=8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, star formation rates, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z=8. We also provide new compilations of cosmic and specific star formation rates; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 10^12 Msun are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for star formation histories that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the star formation history of galaxies into a self-consistent framework based on the modern understanding of structure formation in LCDM. Constraints on the stellar mass-halo mass relationship and star formation rates are available for download at http://www.peterbehroozi.com/data.html .

preprint2013arXiv

The Densest Galaxy

We report the discovery of a remarkable ultra-compact dwarf galaxy around the massive Virgo elliptical galaxy NGC 4649 (M60), which we term M60-UCD1. With a dynamical mass of 2.0 x 10^8 M_sun but a half-light radius of only ~ 24 pc, M60-UCD1 is more massive than any ultra-compact dwarfs of comparable size, and is arguably the densest galaxy known in the local universe. It has a two-component structure well-fit by a sum of Sersic functions, with an elliptical, compact (r_h=14 pc; n ~ 3.3) inner component and a round, exponential, extended (r_h=49 pc) outer component. Chandra data reveal a variable central X-ray source with L_X ~ 10^38 erg/s that could be an active galactic nucleus associated with a massive black hole or a low-mass X-ray binary. Analysis of optical spectroscopy shows the object to be old (~> 10 Gyr) and of solar metallicity, with elevated [Mg/Fe] and strongly enhanced [N/Fe] that indicates light element self-enrichment; such self-enrichment may be generically present in dense stellar systems. The velocity dispersion (~ 70 km/s) and resulting dynamical mass-to-light ratio (M/L_V=4.9 +/- 0.7) are consistent with---but slightly higher than---expectations for an old, metal-rich stellar population with a Kroupa initial mass function. The presence of a massive black hole or a mild increase in low-mass stars or stellar remnants is therefore also consistent with this M/L_V. The stellar density of the galaxy is so high that no dynamical signature of dark matter is expected. However, the properties of M60-UCD1 suggest an origin in the tidal stripping of a nucleated galaxy with M_B ~ -18 to -19.

preprint2013arXiv

The Dust Attenuation Law in Distant Galaxies: Evidence for Variation with Spectral Type

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

preprint2013arXiv

The stellar IMF determined in early-type galaxies from a non-degenerate set of optical line indices

We investigate the optical spectral region of spectra of 1000 stars searching for IMF-sensitive features to constrain the low-mass end of the initial mass function (IMF) slope in elliptical galaxies. We use the MILES stellar library in the wavelength range [3500-7500] Å to select indices that are sensitive to cool dwarf stars and that only weakly depend on age and metallicity. In this wavelength range, the response of a change in the effective temperature of the cool red giant population is similar to the response of a change in the number of dwarf stars in the galaxy. We investigate the degeneracy between IMF variation and Teff,RGB and show that it is possible to break it with the new IMF indicators defined here. In particular, we define a CaH1 index around λ6380 Å that arises purely from cool dwarfs, does not strongly depend on age and is anti-correlated with [a/Fe]. This index allows the determination of the the IMF slope when combined with TiO lines and age- and metallicity-dependent features. The use of several indicators is crucial to break degeneracies between IMF, age, abundance pattern and Teff of the RGB population. We measure line-index strengths of the optical IMF indicators in the SSP models and compare these with index strengths of the same features in stacked SSDSS early-type galaxy spectra with varying velocity dispersions. Using 10 different indicators, we find a clear trend of a steepening IMF with increasing velocity dispersion described by the equation x=(2.3+-0.1)*logσ200+(2.13+-0.15), where x is the IMF slope and σ200 is the central stellar velocity dispersion measured in units of 200km/s. We compute the 'IMF mismatch' parameter as the ratio of stellar mass-to-light ratio predicted from the x-σ200 relation to that inferred from SSP models assuming a Salpeter IMF and find good agreement with independent published results.

preprint2013arXiv

The Strikingly Similar Relation between Satellite and Central Galaxies and Their Dark Matter Halos Since z=2

Satellite galaxies in rich clusters are subject to numerous physical processes that can significantly influence their evolution. However, the typical L* satellite galaxy resides in much lower mass galaxy groups, where the processes capable of altering their evolution are generally weaker and have had less time to operate. To investigate the extent to which satellite and central galaxy evolution differs, we separately model the stellar mass - halo mass (M* -Mh) relation for these two populations over the redshift interval 0 < z < 1. This relation for central galaxies is constrained by the galaxy stellar mass function while the relation for satellite galaxies is constrained against recent measurements of the galaxy two-point correlation function (2PCF). At z ~ 0 the satellites, on average, have ~10% larger stellar masses at fixed peak subhalo mass compared to central galaxies of the same halo mass. This is required in order to reproduce the observed stellar mass-dependent 2PCF and satellite fractions. At low masses our model slightly under-predicts the correlation function at ~1 Mpc scales. At z ~ 1 the satellite and central galaxy M*-Mh relations are consistent within the errors, and the model provides an excellent fit to the clustering data. At present, the errors on the clustering data at z ~ 2 are too large to constrain the satellite model. A simple model in which satellite and central galaxies share the same M*-Mh relation is able to reproduce the extant z ~ 2 clustering data. We speculate that the striking similarity between the satellite and central galaxy M*-Mh relations since z ~ 2 arises because the central galaxy relation evolves very weakly with time and because the stellar mass of the typical satellite galaxy has not changed significantly since it was accreted. [Abridged]

preprint2012arXiv

A Simple Model for Quasar Demographics

We present a simple model for the relationship between quasars, galaxies, and dark matter halos from 0.5<z<6. In the model, black hole (BH) mass is linearly related to galaxy mass, and galaxies are connected to dark matter halos via empirically constrained relations. A simple "scattered" light bulb model for quasars is adopted, wherein BHs shine at a fixed fraction of the Eddington luminosity during accretion episodes, and Eddington ratios are drawn from a lognormal distribution that is redshift-independent. This model has two free, physically meaningful parameters at each redshift: the normalization of the Mbh-Mgal relation and the quasar duty cycle; these parameters are fit to the observed quasar luminosity function (LF) over the interval 0.5<z<6. This simple model provides an excellent fit to the LF at all epochs, and also successfully predicts the observed projected two-point correlation of quasars from 0.5<z<2.5. It is significant that a single quasar duty cycle at each redshift is capable of reproducing the extant observations. The data are therefore consistent with a scenario wherein quasars are equally likely to exist in galaxies, and therefore dark matter halos, over a wide range in masses. The knee in the quasar LF is a reflection of the knee in the stellar mass-halo mass relation. Future constraints on the quasar LF and quasar clustering at high redshift will provide strong constraints on the model. In the model, the autocorrelation function of quasars becomes a strong function of luminosity only at the very highest luminosities, and will be difficult to observe because such quasars are so rare. Cross-correlation techniques may provide useful constraints on the bias of such rare objects.

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

NIMBUS: The Near-Infrared Multi-Band Ultraprecise Spectroimager for SOFIA

We present a new and innovative near-infrared multi-band ultraprecise spectroimager (NIMBUS) for SOFIA. This design is capable of characterizing a large sample of extrasolar planet atmospheres by measuring elemental and molecular abundances during primary transit and occultation. This wide-field spectroimager would also provide new insights into Trans-Neptunian Objects (TNO), Solar System occultations, brown dwarf atmospheres, carbon chemistry in globular clusters, chemical gradients in nearby galaxies, and galaxy photometric redshifts. NIMBUS would be the premier ultraprecise spectroimager by taking advantage of the SOFIA observatory and state of the art infrared technologies. This optical design splits the beam into eight separate spectral bandpasses, centered around key molecular bands from 1 to 4 microns. Each spectral channel has a wide field of view for simultaneous observations of a reference star that can decorrelate time-variable atmospheric and optical assembly effects, allowing the instrument to achieve ultraprecise calibration for imaging and photometry for a wide variety of astrophysical sources. NIMBUS produces the same data products as a low-resolution integral field spectrograph over a large spectral bandpass, but this design obviates many of the problems that preclude high-precision measurements with traditional slit and integral field spectrographs. This instrument concept is currently not funded for development.

preprint2012arXiv

On The Birth Masses of the Ancient Globular Clusters

All globular clusters (GCs) studied to date show evidence for internal variation in their light element abundances. These variations have been interpreted as evidence for multiple star formation episodes within GCs, with secondary episodes fueled, at least in part, by the ejecta of AGB stars from a first generation of stars. A major puzzle emerging from this otherwise plausible scenario is that the fraction of stars associated with the second episode of star formation is observed to be much larger than expected for a standard IMF. The present work investigates this tension by modeling the observed anti-correlation between [Na/Fe] and [O/Fe] for 20 Galactic GCs. If the abundance pattern of the retained AGB ejecta does not depend on GC mass at fixed [Fe/H], then a strong correlation is found between the fraction of current GC stellar mass comprised of pure AGB ejecta, f_p, and GC mass. This fraction varies from 0.20 at low masses (10^4.5 Msun) to 0.45 at high masses (10^6.5 Msun). The fraction of mass associated with pure AGB ejecta is directly related to the total mass of the cluster at birth; the ratio between the initial and present mass in stars can therefore be derived. Assuming a star formation efficiency of 50%, the observed Na-O anti-correlations imply that GCs were at least 10-20 times more massive at birth. These factors are lower limits because any mass-loss mechanism that removes first and second generation stars equally will leave f_p unchanged. The mass-dependence of f_p probably arises because lower mass GCs are unable to retain all of the AGB ejecta from the first stellar generation. Recent observations of elemental abundances in intermediate-age LMC clusters are re-interpreted and shown to be consistent with this basic scenario. A convincing explanation of these trends is currently lacking.

preprint2012arXiv

On the Lack of Evolution in Galaxy Star Formation Efficiency

Using reconstructed galaxy star formation histories, we calculate the instantaneous efficiency of galaxy star formation (i.e., the star formation rate divided by the baryon accretion rate) from $z=8$ to the present day. This efficiency exhibits a clear peak near a characteristic halo mass of 10^11.7 Msun, which coincides with longstanding theoretical predictions for the mass scale relevant to virial shock heating of accreted gas. Above the characteristic halo mass, the efficiency falls off as the mass to the minus four-thirds power; below the characteristic mass, the efficiency falls off at an average scaling of mass to the two-thirds power. By comparison, the shape and normalization of the efficiency change very little since z=4. We show that a time-independent star formation efficiency simply explains the shape of the cosmic star formation rate since z=4 in terms of dark matter accretion rates. The rise in the cosmic star formation from early times until z=2 is especially sensitive to galaxy formation efficiency. The mass dependence of the efficiency strongly limits where most star formation occurs, with the result that two-thirds of all star formation has occurred inside halos within a factor of three of the characteristic mass, a range that includes the mass of the Milky Way.

preprint2012arXiv

Runaway Stars and the Escape of Ionizing Radiation from High-Redshift Galaxies

Approximately 30% of all massive stars in the Galaxy are runaways with velocities exceeding 30 km/s. Their high speeds allow them to travel ~0.1-1 kpc away from their birth place before they explode at the end of their several Myr lifetimes. At high redshift, when galaxies were much smaller than in the local universe, runaways could venture far from the dense inner regions of their host galaxies. From these large radii, and therefore low column densities, much of their ionizing radiation is able to escape into the intergalactic medium. Runaways may therefore significantly enhance the overall escape fraction of ionizing radiation, fesc, from small galaxies at high redshift. We present simple models of the high-redshift runaway population and its impact on fesc as a function of halo mass, size, and redshift. We find that the inclusion of runaways enhances fesc by factors of ~1.1-8, depending on halo mass, galaxy geometry, and the mechanism of runaway production, implying that runaways may contribute 50-90% of the total ionizing radiation escaping from high-redshift galaxies. Runaways may therefore play an important role in reionizing the universe.

preprint2012arXiv

The DEEP2 Galaxy Redshift Survey: Design, Observations, Data Reduction, and Redshifts

We describe the design and data sample from the DEEP2 Galaxy Redshift Survey, the densest and largest precision-redshift survey of galaxies at z ~ 1 completed to date. The survey has conducted a comprehensive census of massive galaxies, their properties, environments, and large-scale structure down to absolute magnitude M_B = -20 at z ~ 1 via ~90 nights of observation on the DEIMOS spectrograph at Keck Observatory. DEEP2 covers an area of 2.8 deg^2 divided into four separate fields, observed to a limiting apparent magnitude of R_AB=24.1. Objects with z < 0.7 are rejected based on BRI photometry in three of the four DEEP2 fields, allowing galaxies with z > 0.7 to be targeted ~2.5 times more efficiently than in a purely magnitude-limited sample. Approximately sixty percent of eligible targets are chosen for spectroscopy, yielding nearly 53,000 spectra and more than 38,000 reliable redshift measurements. Most of the targets which fail to yield secure redshifts are blue objects that lie beyond z ~ 1.45. The DEIMOS 1200-line/mm grating used for the survey delivers high spectral resolution (R~6000), accurate and secure redshifts, and unique internal kinematic information. Extensive ancillary data are available in the DEEP2 fields, particularly in the Extended Groth Strip, which has evolved into one of the richest multiwavelength regions on the sky. DEEP2 surpasses other deep precision-redshift surveys at z ~ 1 in terms of galaxy numbers, redshift accuracy, sample number density, and amount of spectral information. We also provide an overview of the scientific highlights of the DEEP2 survey thus far. This paper is intended as a handbook for users of the DEEP2 Data Release 4, which includes all DEEP2 spectra and redshifts, as well as for the publicly-available DEEP2 DEIMOS data reduction pipelines. [Abridged]

preprint2012arXiv

The progenitors of present-day massive red galaxies up to z ~ 0.7 - finding passive galaxies using SDSS-I/II and SDSS-III

We present a comprehensive study of 250,000 galaxies targeted by the Baryon Oscillation Spectroscopic Survey (BOSS) up to z ~ 0.7 with the specific goal of identifying and characterising a population of galaxies that has evolved without significant merging. We compute a likelihood that each BOSS galaxy is a progenitor of the Luminous Red Galaxies (LRGs) sample, targeted by SDSS-I/II up z ~ 0.5, by using the fossil record of LRGs and their inferred star-formation histories, metallicity histories and dust content. We determine merger rates, luminosity growth rates and the evolution of the large-scale clustering between the two surveys, and we investigate the effect of using different stellar population synthesis models in our conclusions. We demonstrate that our sample is slowly evolving (of the order of 2 +/- 1.5% per Gyr by merging). Our conclusions refer to the bright and massive end of the galaxy population, with Mi0.55 < -22, and M* > 1E11.2 Msolar, corresponding roughly to 95% and 40% of the LRGs and BOSS galaxy populations, respectively. Our analysis further shows that any possible excess of flux in BOSS galaxies, when compared to LRGs, from potentially unresolved targets at z ~ 0.55 must be less than 1% in the r0.55-band (approximately equivalent to the g-band in the rest-frame of galaxies at z=0.55). When weighting the BOSS galaxies based on the predicted properties of the LRGs, and restricting the analysis to the reddest BOSS galaxies, we find an evolution of the large-scale clustering that is consistent with dynamical passive evolution, assuming a standard cosmology. We conclude that our likelihoods give a weighted sample that is as clean and as close to passive evolution (in dynamical terms, i.e. no or negligible merging) as possible, and that is optimal for cosmological studies.

preprint2012arXiv

The SLUGGS Survey: NGC 3115, A Critical Test Case for Metallicity Bimodality in Globular Cluster Systems

Due to its proximity (9 Mpc) and the strongly bimodal color distribution of its spectroscopically well-sampled globular cluster (GC) system, the early-type galaxy NGC 3115 provides one of the best available tests of whether the color bimodality widely observed in GC systems generally reflects a true metallicity bimodality. Color bimodality has alternatively been attributed to a strongly nonlinear color--metallicity relation reflecting the influence of hot horizontal branch stars. Here we couple Subaru Suprime-Cam gi photometry with Keck/DEIMOS spectroscopy to accurately measure GC colors and a CaT index that measures the CaII triplet. We find the NGC 3115 GC system to be unambiguously bimodal in both color and the CaT index. Using simple stellar population models, we show that the CaT index is essentially unaffected by variations in horizontal branch morphology over the range of metallicities relevant to GC systems (and is thus a robust indicator of metallicity) and confirm bimodality in the metallicity distribution. We assess the existing evidence for and against multiple metallicity subpopulations in early and late-type galaxies and conclude that metallicity bi/multimodality is common. We briefly discuss how this fundamental characteristic links directly to the star formation and assembly histories of galaxies.

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

Are Halo and Galaxy Formation Histories Correlated?

The properties of dark matter halos, including mass growth, correlate with larger scale environment at fixed mass, an effect known as assembly bias. However, whether this environmental dependence manifests itself in galaxy properties remains unclear. We apply a group-finding algorithm to DR7 of the SDSS to estimate the halo mass of each galaxy and to decompose galaxies into those that exist at the centers of distinct halos and those that orbit as satellites within larger halos. Using the 4000-A break as a measure of star formation history, we examine the correlation between the quenched fraction of galaxies, f_q, and large-scale environment, rho. At all galaxy magnitudes, there is a positive, monotonic relationship between f_q and rho. We use the group catalog to decompose this correlation into the contribution from central and satellite galaxies as a function of halo mass. Because satellites are more likely to be quenched than central galaxies, the observed f_q-rho correlation is primarily due to variations of the halo mass function with environment, which causes a larger fraction of satellite galaxies at high rho. For low-mass central galaxies (Mgal <~ 10^10.0 Msol/h^2), there is no correlation between f_q and rho. These results are inconsistent with the strong assembly bias of dark matter halos seen in this mass regime if recent galaxy growth at all correlates with recent halo growth, as we demonstrate through a high resolution N-body simulation. We also find that the mean stellar age of quenched central galaxies is independent of rho at fixed Mgal, while the formation times of low mass halos vary significantly. We conclude that the processes that halt the star formation of low mass central galaxies are not correlated to the formation histories of their host halos, and old galaxies do not reside preferentially in old halos. (Abridged)

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

Dark halo response and the stellar initial mass function in early-type and late-type galaxies

We investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

preprint2011arXiv

Empirical Constraints on the Coevolution of Supermassive Black Holes and their Host Spheroids

We investigate the evolution of the MBH-σ relation by examining the relationship between the intrinsic scatter in the MBH-σ relation and galaxy bolometric nuclear luminosity, the latter being a probe of the accretion rate of the black hole (BH). Our sample is composed of galaxies with classical bulges when possible, of which 38 have dynamically measured BHs masses, and 17 have BHs masses measured by reverberation mapping. In order to obtain the bolometric nuclear luminosity for galaxies with low nuclear luminosity, we convert the X-ray nuclear luminosity measured by Chandra to bolometric luminosity. We find that the scatter in the MBH-σ relation is uncorrelated with nuclear luminosity over seven orders of magnitude in luminosity, with the high luminosity end approaching the Eddington luminosity. This suggests that at the present epoch galaxies evolve along the MBH-σ relation. This conclusion is consistent with the standard paradigm that BHs grow contemporaneously with their host stellar spheroids.

preprint2011arXiv

Evidence Against Dark Matter Halos Surrounding the Globular Clusters MGC1 and NGC 2419

The conjecture that the ancient globular clusters (GCs) formed at the center of their own dark matter halos was first proposed by Peebles (1984), and has recently been revived to explain the puzzling abundance patterns observed within many GCs. In this paper we demonstrate that the outer stellar density profile of isolated GCs is very sensitive to the presence of an extended dark halo. The GCs NGC 2419, located at 90 kpc from the center of our Galaxy, and MGC1, located at ~200 kpc from the center of M31, are ideal laboratories for testing the scenario that GCs formed at the centers of massive dark halos. Comparing analytic models to observations of these GCs, we conclude that these GCs cannot be embedded within dark halos with a virial mass greater than 10^6 Msun, or, equivalently, the dark matter halo mass-to-stellar mass ratio must be Mdm/M_*<1. If these GCs have indeed orbited within weak tidal fields throughout their lifetimes, then these limits imply that these GCs did not form within their own dark halos. Recent observations of an extended stellar halo in the GC NGC 1851 are also interpreted in the context of our analytic models. Implications of these results for the formation of GCs are briefly discussed.

preprint2011arXiv

Improved Constraints on Type Ia Supernova Host Galaxy Properties using Multi-Wavelength Photometry and their Correlations with Supernova Properties

We improve estimates of stellar mass and mass-weighted average age of Type Ia supernova (SN Ia) host galaxies by combining UV and near-IR photometry with optical photometry in our analysis. Using 206 SNe Ia drawn from the full three-year SDSS-II Supernova Survey (median redshift of z {\approx} 0.2) and multi-wavelength host-galaxy photometry from SDSS, GALEX, and UKIDSS, we present evidence of a correlation (1.9σ confidence level) between the residuals of SNe Ia about the best-fit Hubble relation and the mass-weighted average age of their host galaxies. The trend is such that older galaxies host SNe Ia that are brighter than average after standard light-curve corrections are made. We also confirm, at the 3.0σ level, the trend seen by previous studies that more massive galaxies often host brighter SNe Ia after light-curve correction.

preprint2011arXiv

The integrated stellar content of dark matter halos

Measurements of the total amount of stars locked up in galaxies as a function of host halo mass contain key clues about the efficiency of processes that regulate star formation. We derive the total stellar mass fraction f_star as a function of halo mass M500c from z=0.2 to z=1 using two complementary methods. First, we derive f_star using a statistical Halo Occupation Distribution model jointly constrained by data from lensing, clustering, and the stellar mass function. This method enables us to probe f_star over a much wider halo mass range than with group or cluster catalogs. Second, we derive f_star at group scales using a COSMOS X-ray group catalog and we show that the two methods agree to within 30%. We quantify the systematic uncertainty on f_star using abundance matching methods and we show that the statistical uncertainty on f_star (~10%) is dwarfed by systematic uncertainties associated with stellar mass measurements (~45% excluding IMF uncertainties). Assuming a Chabrier IMF, we find 0.012<f_star<0.025 at M500c=10^13 Msun and 0.0057<f_star<0.015 at M500c=10^14 Msun. These values are significantly lower than previously published estimates. We investigate the cause of this difference and find that previous work has overestimated f_star due to a combination of inaccurate stellar mass estimators and/or because they have assumed that all galaxies in groups are early type galaxies with a constant mass-to-light ratio. Contrary to previous claims, our results suggest that the mean value of f_star is always significantly lower than f_gas for halos above 10^13 Msun. Combining our results with recently published gas mas fractions, we find a shortfall in f_star+f_gas at R500c compared to the cosmic mean. This shortfall varies with halo mass and becomes larger towards lower halos masses.

preprint2010arXiv

A Comprehensive Analysis of Uncertainties Affecting the Stellar Mass - Halo Mass Relation for 0<z<4

We conduct a comprehensive analysis of the relationship between central galaxies and their host dark matter halos, as characterized by the stellar mass-halo mass (SM-HM) relation, with rigorous consideration of uncertainties. Our analysis focuses on results from the abundance matching technique, which assumes that every dark matter halo or subhalo above a specific mass threshold hosts one galaxy. We discuss the quantitative effects of uncertainties in observed galaxy stellar mass functions (GSMFs) (including stellar mass estimates and counting uncertainties), halo mass functions (including cosmology and uncertainties from substructure), and the abundance matching technique used to link galaxies to halos (including scatter in this connection). Our analysis results in a robust estimate of the SM-HM relation and its evolution from z=0 to z=4. The shape and evolution are well constrained for z < 1. The largest uncertainties at these redshifts are due to stellar mass estimates; however, failure to account for scatter in stellar masses at fixed halo mass can lead to errors of similar magnitude in the SM-HM relation for central galaxies in massive halos. We also investigate the SM-HM relation to z=4, although the shape of the relation at higher redshifts remains fairly unconstrained when uncertainties are taken into account. These results will provide a powerful tool to inform galaxy evolution models. [Abridged]

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

Dust attenuation in disk-dominated galaxies: evidence for the 2175A dust feature

The attenuation of starlight by interstellar dust is investigated in a sample of low redshift, disk-dominated star-forming galaxies using photometry from GALEX and SDSS. By considering broadband colors as a function of galaxy inclination we are able to confidently separate trends arising from increasing dust opacity from possible differences in stellar populations, since stellar populations do not correlate with inclination. All commonly employed dust attenuation curves (such as the Calzetti curve for starbursts, or a power-law curve) provide poor fits to the ultraviolet colors for moderately and highly inclined galaxies. This conclusion rests on the fact that the average FUV-NUV color varies little from face-on to edge-on galaxies, while other colors such as NUV-u and u-r vary strongly with inclination. After considering a number of model variations, we are led to speculate that the presence of the strong dust extinction feature at 2175A seen in the Milky Way (MW) extinction curve is responsible for the observed trends. Independent of our interpretation, these results imply that the modeling of dust attenuation in the ultraviolet is significantly more complicated than traditionally assumed. These results also imply a very weak dependence of the FUV-NUV color on total FUV attenuation, and we conclude from this that it is extremely difficult to use only the observed UV spectral slope to infer the total UV dust attenuation, as is commonly done. We propose several simple tests that might finally identify the grain population responsible for the 2175A feature.

preprint2010arXiv

On the Formation of Multiple Stellar Populations in Globular Clusters

Nearly all globular clusters (GCs) studied to date show evidence for multiple stellar populations, in stark contrast to the conventional view that GCs are a mono-metallic, coeval population of stars. Building on earlier work, we propose a simple physical model for the early evolution (several 10^8 yr) of GCs. We consider the effects of stellar mass-loss, type II and prompt type Ia supernovae, ram pressure, and accretion from the ambient ISM on the development of a young GC's own gas reservoir. In our model, type II SNe from a first generation of star formation clears the GC of its initial gas reservoir. Over the next several 10^8 yr, mass lost from AGB stars and matter accreted from the ambient ISM collect at the center of the GC. This material must remain quite cool (T~10^2K), but does not catastrophically cool on a crossing time because of the high Lyman-Werner flux density in young GCs. The collection of gas within the GC must compete with ram pressure from the ambient ISM. After several 10^8 yr, the Lyman-Werner photon flux density drops by more than three orders of magnitude, allowing molecular hydrogen and then stars to form. After this second generation of star formation, type II SNe from the second generation and then prompt type Ia SNe associated with the first generation maintain a gas-free GC, thereby ending the cycle of star formation events. Our model makes clear predictions for the presence or absence of multiple stellar populations within GCs as a function of GC mass and formation environment. Analyzing intermediate-age LMC clusters, we find evidence for a mass threshold of ~10^4 Msun below which LMC clusters appear to be truly coeval. This threshold mass is consistent with our predictions for the mass at which ram pressure is capable of clearing gas from clusters in the LMC at the present epoch. (ABRIDGED)

preprint2010arXiv

The Kinematic Connection Between Galaxies and Dark Matter Haloes

Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt \ne V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.

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.

preprint2009arXiv

Connecting Galaxies, Halos, and Star Formation Rates Across Cosmic Time

A simple, observationally-motivated model is presented for understanding how halo masses, galaxy stellar masses, and star formation rates are related, and how these relations evolve with time. The relation between halo mass and galaxy stellar mass is determined by matching the observed spatial abundance of galaxies to the expected spatial abundance of halos at multiple epochs -- i.e. more massive galaxies are assigned to more massive halos at each epoch. Halos at different epochs are connected by halo mass accretion histories estimated from N-body simulations. The halo--galaxy connection at fixed epochs in conjunction with the connection between halos across time provides a connection between observed galaxies across time. With approximations for the impact of merging and accretion on the growth of galaxies, one can then directly infer the star formation histories of galaxies as a function of stellar and halo mass. This model is tuned to match both the observed evolution of the stellar mass function and the normalization of the observed star formation rate -- stellar mass relation to z~1. The data demands, for example, that the star formation rate density is dominated by galaxies with Mstar~10^{10.0-10.5} Msun from 0<z<1, and that such galaxies over these epochs reside in halos with M~10^{11.5-12.5} Msun. The star formation rate -- halo mass relation is approximately Gaussian over the range 0<z<1 with a mildly evolving mean and normalization. This model is then used to shed light on a number of issues, including 1) a clarification of ``downsizing'', 2) the lack of a sharp characteristic halo mass at which star formation is truncated, and 3) the dominance of star formation over merging to the stellar build-up of galaxies with Mstar<10^11 Msun at z<1.

preprint2009arXiv

Dust attenuation in the restframe ultraviolet: constraints from star-forming galaxies at z~1

A novel technique is employed for estimating attenuation curves in galaxies where only photometry and spectroscopic redshifts are available. This technique provides a powerful measure of particular extinction features such as the UV bump at 2175\A, which has been observed in environments ranging from the Milky Way to high-redshift star-forming galaxies. Knowledge of the typical strength of the UV bump as a function of environment and redshift is crucial for converting restframe UV flux into star formation rates. The UV bump will impart a unique signature as it moves through various filters due to redshifting; its presence can therefore be disentangled from other stellar population effects. The utility of this technique is demonstrated with a large sample of galaxies drawn from the DEEP2 Galaxy Redshift Survey. The observed B-R color of star-forming galaxies at 0.6<z<1.4 disfavors the presence of a UV bump as strong as observed in the Milky Way, and instead favors restframe UV (1800A<lambda<3000A) attenuation curves similar to the Milky Way without a UV bump or a power-law with index delta=-0.7. Stronger constraints on the strength of the UV bump in galaxies can be achieved if independent constraints on the V-band optical depth are available.

preprint2009arXiv

The propagation of uncertainties in stellar population synthesis modeling III: model calibration, comparison, and evaluation

Stellar population synthesis (SPS) provides the link between the stellar and dust content of galaxies and their observed spectral energy distributions. In the present work we perform a comprehensive calibration of our own flexible SPS (FSPS) model against a suite of data. Several public SPS models are intercompared, including the models of Bruzual & Charlot (BC03), Maraston (M05) and FSPS. The relative strengths and weaknesses of these models are evaluated, with the following conclusions: 1) The FSPS and BC03 models compare favorably with MC data at all ages, whereas M05 colors are too red and the age-dependence is incorrect; 2) All models yield similar optical and near-IR colors for old metal-poor systems, and yet they all provide poor fits to the integrated J-K and V-K colors of both MW and M31 star clusters; 4) All models predict ugr colors too red, D4000 strengths too strong and Hdelta strengths too weak compared to massive red sequence galaxies, under the assumption that such galaxies are composed solely of old metal-rich stars; 5) FSPS and, to a lesser extent, BC03 can reproduce the optical and near-IR colors of post-starburst galaxies, while M05 cannot. Reasons for these discrepancies are explored. The failure at predicting the ugr colors, D4000, and Hdelta strengths can be explained by some combination of a minority population of metal-poor stars, young stars, blue straggler and/or blue horizontal branch stars, but not by appealing to inadequacies in either theoretical stellar atmospheres or canonical evolutionary phases (e.g., the main sequence turn-off). We emphasize that due to a lack of calibrating star cluster data in regions of the metallicity-age plane relevant for galaxies, all of these models continue to suffer from serious uncertainties that are difficult to quantify. (ABRIDGED)

Charlie Conroy

What is connected

Connect this record

See the researcher in context

Building this map preview

117 published item(s)

Spectacular nucleosynthesis from early massive stars

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

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

A massive AGB donor in Scutum X-1: Identification of the first Mira variable in an X-ray binary

A Tilt in the Dark Matter Halo of the Galaxy

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

Deriving ages and horizontal branch properties of integrated stellar populations

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

How Well Can We Measure Galaxy Dust Attenuation Curves? The Impact of the Assumed Star-Dust Geometry Model in SED Fitting

Live Fast, Die $α$-Enhanced: The Mass-Metallicity-$α$ Relation of the Milky Way's Disrupted Dwarf Galaxies

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

SDSS-IV MaNGA: How the stellar populations of passive central galaxies depend on stellar and halo mass

SRGA J181414.6-225604: A new Galactic symbiotic X-ray binary outburst triggered by an intense mass loss episode of a heavily obscured Mira variable

Stellar halos from The Dragonfly Edge-on Galaxies Survey

Unicorns and Giraffes in the binary zoo: stripped giants with subgiant companions

Wide binaries from the H3 survey: the thick disk and halo have similar wide binary fractions

Evidence from Disrupted Halo Dwarfs that $r$-process Enrichment via Neutron Star Mergers is Delayed by $\gtrsim500$ Myrs

Measuring distances to low-luminosity galaxies using surface brightness fluctuations

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

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

Stellar labels for hot stars from low-resolution spectra - I. the HotPayne method and results for 330,000 stars from LAMOST DR6

The Mass of the Milky Way from the H3 Survey

A Diffuse Metal-Poor Component of the Sagittarius Stream Revealed by the H3 Survey

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

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

Ancient Very Metal-Poor Stars Associated With the Galactic Disk in the H3 Survey

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

High-resolution spectroscopy of the GD-1 stellar stream localizes the perturber near the orbital plane of Sagittarius

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

MINESweeper: Spectrophotometric Modeling of Stars in the Gaia Era

Powderday: Dust Radiative Transfer for Galaxy Simulations

Rapid Reionization by the Oligarchs: The Case for Massive, UV-Bright, Star-Forming Galaxies with High Escape Fractions

SPECULATOR: Emulating stellar population synthesis for fast and accurate galaxy spectra and photometry

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

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

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

A Lower Limit on the Mass of Our Galaxy from the H3 Survey

Combined Effects of Rotation and Age Spreads on Extended Main Sequence Turn Offs

Resolving the Metallicity Distribution of the Stellar Halo with the H3 Survey

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

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

Accelerated Fitting of Stellar Spectra

Constructing Polynomial Spectral Models for Stars

Exploring Systematic Effects in the Relation Between Stellar Mass, Gas Phase Metallicity, and Star Formation Rate

MESA Isochrones and Stellar Tracks (MIST). I: Solar-Scaled Models

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

Redshift Measurement and Spectral Classification for eBOSS Galaxies with the Redmonster Software

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

APOGEE chemical tagging constraint on the maximum star cluster mass in the $α$-enhanced Galactic disk

Beacons In the Dark: Using Novae and Supernovae to Detect Dwarf Galaxies in the Local Universe

Circumstellar Dust Around AGB Stars and Implications for Infrared Emission from Galaxies

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

From Cosmic Birth to Living Earths: The Future of UVOIR Space Astronomy

Hiding in plain sight: record-breaking compact stellar systems in the Sloan Digital Sky Survey

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

Prospects for chemically tagging stars in the Galaxy

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

Stellar Populations of Barred Quiescent Galaxies

The IMF-sensitive 1.14-micron Na I doublet in early-type galaxies

The Local Group as a time machine: studying the high-redshift Universe with nearby galaxies

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

The SINFONI Nearby Elliptical Lens Locator Survey: Discovery of two new low-redshift strong lenses and implications for the initial mass function in giant early-type galaxies

Comparing Simple Quasar Demographics Models

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

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

Fluctuation Spectroscopy: A New Probe of Old Stellar Populations

Galaxy evolution near groups and clusters: ejected satellites and the spatial extent of environmental quenching

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

Parametrising Star Formation Histories

Reversal of Fortune: Increased Star Formation Efficiencies in the Early Histories of Dwarf Galaxies?

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

Stellar Models of Multiple Populations in Globular Clusters. I. The Main Sequence of NGC 6752

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

The Baryon Cycle of Dwarf Galaxies: Dark, Bursty, Gas-Rich Polluters