Researcher profile

Charlie Conroy

Charlie Conroy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
39works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

39 published item(s)

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&#39;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 &#34;late bloomers&#34; are uncommon. In summary, our analysis supports the &#34;grow & quench&#34; 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&#39;s Disrupted Dwarf Galaxies

The Milky Way&#39;s satellite galaxies (&#34;surviving dwarfs&#34;) have been studied for decades as unique probes of chemical evolution in the low-mass regime. Here we extend such studies to the &#34;disrupted dwarfs&#34;, 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&#39;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&#39;s Galaxy Candidate: Puzzlingly Luminous at $z\approx17$, or Dusty/Quenched at $z\approx5$?

$JWST$&#39;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&#39;s nearest neighbors at $<2.5$&#34; 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&#39;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 (&#34;the Unicorn&#34;) and 2M04123153+6738486 (&#34;the Giraffe&#34;). 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&#39;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&#39; 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&#39;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 &#34;sampling scatter&#34; 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 &#34;normal&#34; 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 &#34;inside-out&#34; 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&#39;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&#39;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&#39;) 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 &#34;Needle Galaxy&#34;), 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&#39;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 &#34;in-situ&#34; 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)