Showing 1–50 of 136 results for author: Sanderson, R

The imprint of dark matter on the Galactic acceleration field

Authors: Arpit Arora, Robyn E. Sanderson, Sukanya Chakrabarti, Andrew Wetzel, Thomas Donlon II, Danny Horta, Sarah R. Loebman, Lina Necib, Micah Oeur

Abstract: Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, from pulsar timing, and from eclipsing binary stars in the Solar Neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density, and from this, by accoun… ▽ More Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, from pulsar timing, and from eclipsing binary stars in the Solar Neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density, and from this, by accounting for the mass in stars, gas, and dust, we can infer the amount of dark matter. Leveraging the FIRE-2 simulations of Milky Way-mass galaxies, we compare vertical acceleration profiles between cold dark matter (CDM) and self-interacting dark matter (SIDM) with constant cross-section of 1 cm$^2$ g$^{-1}$ across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos consistently exhibit higher local stellar and dark matter densities and steeper vertical acceleration gradients, up to 30% steeper near the Solar Neighborhood. SIDM halos also manifest a more oblate halo shape in the Solar Neighborhood. Furthermore, enhanced precision in acceleration measurements and larger datasets promise to provide better constraints on the local dark matter density, complementing our understanding from kinematic analysis of their distribution within galaxies. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 14 pages, 5 figures, and 2 tables. Submitted to APJ

Early and elongated epochs of planetesimal dynamo generation

Authors: Hannah R. Sanderson, James F. J. Bryson, Claire I. O. Nichols

Abstract: Accreting in the first few Ma after Solar System formation, planetesimals record conditions in the protoplanetary disc and are the remnants of planetary formation processes. The meteorite paleomagnetic record carries key insights into the thermal history of planetesimals and their extent of differentiation. The current paradigm splits the paleomagnetic record into three magnetic field generation e… ▽ More Accreting in the first few Ma after Solar System formation, planetesimals record conditions in the protoplanetary disc and are the remnants of planetary formation processes. The meteorite paleomagnetic record carries key insights into the thermal history of planetesimals and their extent of differentiation. The current paradigm splits the paleomagnetic record into three magnetic field generation epochs: an early nebula field (<5Ma after CAI formation), followed by thermal dynamos (5-34 Ma after CAI formation), then a gap in dynamo generation, before the onset of core solidification and compositional dynamos. The split between these epochs has been defined using thermal evolution and dynamo generation models of planetesimals. Here we demonstrate these epochs are not as distinct as previously thought based on our refined thermal evolution model that includes more realistic parametrisations for mantle convection, non-eutectic core solidification and radiogenic $^{60}Fe$ in the core. Inclusion of $^{60}$ in the core brings forward the onset of dynamo generation to 1-2 Ma after CAI formation, which overlaps with the existence of the nebula field. The second epoch of dynamo generation begins prior to the onset of core solidification, suggesting this epoch is not purely compositionally driven. Planetesimal radius is the dominant control on dynamo generation, and the choice of reference viscosity can widen the gap between epochs of dynamo generation from 0-200 Ma. Overall, timings of different planetesimal magnetic field generation mechanisms are more variable. This alters the information we can glean from the meteorite paleomagnetic record about the early Solar System. Evidence for the nebula field requires more careful interpretation and young paleomagnetic remanences, for example in the pallasites, may not be evidence for planetesimal core solidification. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 29 pages, 8 figures. Supplementary materials are a pdf in the TeX source file. This is a preliminary release in order to share the manuscript at a conference. This paper will be submitted to Earth and Planetary Science Letters within the next month

Recommendations for Early Definition Science with the Nancy Grace Roman Space Telescope

Authors: Robyn E. Sanderson, Ryan Hickox, Christopher M. Hirata, Matthew J. Holman, Jessica R. Lu, Ashley Villar

Abstract: The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recom… ▽ More The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship observatory, has significant mission time to be spent on surveys for general astrophysics in addition to its three core community surveys. We considered what types of observations outside the core surveys would most benefit from early definition, given 700 hours of mission time in the first two years of Roman's operation. We recommend that a survey of the Galactic plane be defined early, based on the broad range of stakeholders for such a survey, the added scientific value of a first pass to obtain a baseline for proper motions complementary to Gaia's, and the significant potential synergies with ground-based surveys, notably the Legacy Survey of Space and Time (LSST) on Rubin. We also found strong motivation to follow a community definition process for ultra-deep observations with Roman. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Final report of the committee convened by NASA to assess community proposals for early definition science with Roman

Generating synthetic star catalogs from simulated data for next-gen observatories with py-ananke

Authors: Adrien C. R. Thob, Robyn E. Sanderson, Andrew P. Eden, Farnik Nikakhtar, Nondh Panithanpaisal, Nicolás Garavito-Camargo, Sanjib Sharma

Abstract: We find ourselves on the brink of an exciting era in observational astrophysics, driven by groundbreaking facilities like JWST, Euclid, Rubin, Roman, SKA, or ELT. Simultaneously, computational astrophysics has shown significant strides, yielding highly realistic galaxy formation simulations, thanks to both hardware and software enhancements. Bridging the gap between simulations and observations ha… ▽ More We find ourselves on the brink of an exciting era in observational astrophysics, driven by groundbreaking facilities like JWST, Euclid, Rubin, Roman, SKA, or ELT. Simultaneously, computational astrophysics has shown significant strides, yielding highly realistic galaxy formation simulations, thanks to both hardware and software enhancements. Bridging the gap between simulations and observations has become paramount for meaningful comparisons. We introduce py-ananke, a Python pipeline designed to generate synthetic resolved stellar surveys from cosmological simulations, adaptable to various instruments. Building upon its predecessor, ananke by Sanderson et al. 2020 (arXiv:1806.10564), which produced Gaia DR2 mock star surveys, the py-ananke package offers a user-friendly "plug & play" experience. The pipeline employs cutting-edge phase-space density estimation and initial mass function sampling to convert particle data into synthetic stars, while interpolating pre-computed stellar isochrone tracks for photometry. Additionally, it includes modules for estimating interstellar reddening, dust-induced extinctions, and for quantifying errors through dedicated modeling approaches. py-ananke promises to serve as a vital bridge between computational astrophysics and observational astronomy, facilitating preparations and making scientific predictions for the next generation of telescopes. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: Submitted to the Journal of Open Source Software. 12 pages (6 pages of bibliography) and 1 figure. Software repository at https://github.com/athob/py-ananke

On the co-rotation of Milky Way satellites: LMC-mass satellites induce apparent motions in outer halo tracers

Authors: Nicolas Garavito-Camargo, Adrian M. Price-Whelan, Jenna Samuel, Emily C. Cunningham, Ekta Patel, Andrew Wetzel, Kathryn V. Johnston, Arpit Arora, Robyn E. Sanderson, Lehman Garrison, Danny Horta

Abstract: Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing th… ▽ More Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing the inner regions of the halo with respect to the outer halo. Using the Latte suite of FIRE-2 cosmological simulations of MW-mass galaxies, we confirm that the apparent motion of the outer halo induced by the infall of a massive satellite changes the observed distribution of orbital poles of outer-halo tracers, including satellites. We quantify the changes in the distribution of orbital poles using the two-point angular correlation function and find that all satellites induce changes. However, the most massive satellites with pericentric passages between 30-100kpc induce the largest changes. The best LMC-like satellite analog shows the largest change in orbital pole distribution. The dispersion of orbital poles decreases by 20° during the first two pericentric passages. Even when excluding the satellites brought in with the LMC-like satellite, there is clustering of orbital poles. These results suggest that in the MW, the recent pericentric passage of the LMC should have changed the observed distribution of orbital poles of all other satellites. Therefore, studies of kinematically-coherent planes of satellites that seek to place the MW in a cosmological context should account for the existence of a massive satellite like the LMC. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 20 pages, 10 figures. ApJ submitted, Comments are welcome

The Debris of the "Last Major Merger" is Dynamically Young

Authors: Thomas Donlon II, Heidi Jo Newberg, Robyn Sanderson, Emily Bregou, Danny Horta, Arpit Arora, Nondh Panithanpaisal

Abstract: The Milky Way's (MW) inner stellar halo contains an [Fe/H]-rich component with highly eccentric orbits, often referred to as the "last major merger." Hypotheses for the origin of this component include Gaia-Sausage/Enceladus (GSE), where the progenitor collided with the MW proto-disk 8-11 Gyr ago, and the Virgo Radial Merger (VRM), where the progenitor collided with the MW disk within the last 3 G… ▽ More The Milky Way's (MW) inner stellar halo contains an [Fe/H]-rich component with highly eccentric orbits, often referred to as the "last major merger." Hypotheses for the origin of this component include Gaia-Sausage/Enceladus (GSE), where the progenitor collided with the MW proto-disk 8-11 Gyr ago, and the Virgo Radial Merger (VRM), where the progenitor collided with the MW disk within the last 3 Gyr. These two scenarios make different predictions about observable structure in local phase space, because the morphology of debris depends on how long it has had to phase mix. The recently-identified phase-space folds in Gaia DR3 have positive caustic velocities, making them fundamentally different than the phase-mixed chevrons found in simulations at late times. Roughly 20\% of the stars in the prograde local stellar halo are associated with the observed caustics. Based on a simple phase-mixing model, the observed number of caustics are consistent with a merger that occurred 1--2 Gyr ago. We also compare the observed phase-space distribution to FIRE-2 Latte simulations of GSE-like mergers, using a quantitative measurement of phase mixing (2D causticality). The observed local phase-space distribution best matches the simulated data 1--2 Gyr after collision, and certainly not later than 3 Gyr. This is further evidence that the progenitor of the "last major merger" did not collide with the MW proto-disk at early times, as is thought for the GSE, but instead collided with the MW disk within the last few Gyr, consistent with the body of work surrounding the VRM. △ Less

Submitted 8 May, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Bar formation and destruction in the FIRE-2 simulations

Authors: Sioree Ansar, Sarah Pearson, Robyn E. Sanderson, Arpit Arora, Philip F. Hopkins, Andrew Wetzel, Emily C. Cunningham, Jenny Quinn

Abstract: The physical mechanisms responsible for bar formation and destruction in galaxies remain a subject of debate. While we have gained valuable insight into how bars form and evolve from isolated idealized simulations, in the cosmological domain, galactic bars evolve in complex environments with mergers, gas accretion events, in presence of turbulent Inter Stellar Medium (ISM) with multiple star forma… ▽ More The physical mechanisms responsible for bar formation and destruction in galaxies remain a subject of debate. While we have gained valuable insight into how bars form and evolve from isolated idealized simulations, in the cosmological domain, galactic bars evolve in complex environments with mergers, gas accretion events, in presence of turbulent Inter Stellar Medium (ISM) with multiple star formation episodes, in addition to coupling to their host galaxies' dark matter halos. We investigate bar formation in 13 Milky Way-mass galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations. 8 of the 13 simulated galaxies form bars at some point during their history: three from tidal interactions and five from internal evolution of the disk. The bars in FIRE-2 are generally shorter than the corotation radius (mean bar radius $\sim 1.53$ kpc), have a wide range of pattern speeds (36--97 km s$^{-1}$kpc$^{-1}$), and live for a wide range of dynamical times (2--160 bar rotations). We find that bar formation in FIRE-2 galaxies is influenced by satellite interactions and the stellar-to-dark matter mass ratio in the inner galaxy, but neither is a sufficient condition for bar formation. Bar formation is more likely to occur, and the bars formed are stronger and longer-lived, if the disks are kinematically cold; galaxies with high central gas fractions and/or vigorous star formation, on the other hand, tend to form weaker bars. In the case of the FIRE-2 galaxies these properties combine to produce ellipsoidal bars with strengths $A_2/A_0 \sim$ 0.1--0.2. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: 36 pages (including Appendix), 22 figures, to be submitted to the Astrophysical Journal

Analysis of the Usability of Automatically Enriched Cultural Heritage Data

Authors: Julien Antoine Raemy, Robert Sanderson

Abstract: This chapter presents the potential of interoperability and standardised data publication for cultural heritage resources, with a focus on community-driven approaches and web standards for usability. The Linked Open Usable Data (LOUD) design principles, which rely on JSON-LD as lingua franca, serve as the foundation. We begin by exploring the significant advances made by the International Image… ▽ More This chapter presents the potential of interoperability and standardised data publication for cultural heritage resources, with a focus on community-driven approaches and web standards for usability. The Linked Open Usable Data (LOUD) design principles, which rely on JSON-LD as lingua franca, serve as the foundation. We begin by exploring the significant advances made by the International Image Interoperability Framework (IIIF) in promoting interoperability for image-based resources. The principles and practices of IIIF have paved the way for Linked Art, which expands the use of linked data by demonstrating how it can easily facilitate the integration and sharing of semantic cultural heritage data across portals and institutions. To provide a practical demonstration of the concepts discussed, the chapter highlights the implementation of LUX, the Yale Collections Discovery platform. LUX serves as a compelling case study for the use of linked data at scale, demonstrating the real-world application of automated enrichment in the cultural heritage domain. Rooted in empirical study, the analysis presented in this chapter delves into the broader context of community practices and semantic interoperability. By examining the collaborative efforts and integration of diverse cultural heritage resources, the research sheds light on the potential benefits and challenges associated with LOUD. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: This is the preprint version of a chapter submitted to be included in the book "Decoding Cultural Heritage: a critical dissection and taxonomy of human creativity through digital tools", to be published by Springer Nature. The chapter is currently undergoing peer review for potential inclusion in the book

ACM Class: D.2.12; H.1.0; J.5

LMC-driven anisotropic boosts in stream--subhalo interactions

Authors: Arpit Arora, Nicolás Garavito-Camargo, Robyn E. Sanderson, Emily C. Cunningham, Andrew Wetzel, Nondh Panithanpaisal, Megan Barry

Abstract: Dark Matter (DM) subhalos are predicted to perturb stellar streams; stream morphologies and dynamics can constrain the mass distribution of subhalos. Using FIRE-2 simulations of Milky Way-mass galaxies, we show that presence of a Large Magellanic Cloud (LMC)--analog significantly changes stream-subhalo encounter rates. Three key factors drive these changes. First, the LMC--analog brings in many su… ▽ More Dark Matter (DM) subhalos are predicted to perturb stellar streams; stream morphologies and dynamics can constrain the mass distribution of subhalos. Using FIRE-2 simulations of Milky Way-mass galaxies, we show that presence of a Large Magellanic Cloud (LMC)--analog significantly changes stream-subhalo encounter rates. Three key factors drive these changes. First, the LMC--analog brings in many subhalos, increasing encounter rates for streams near the massive satellite by up to 20--40%. Second, the LMC--analog displaces the host from its center-of-mass (inducing reflex motion), causing a north-south asymmetry in the density and radial velocity distribution of subhalos. This asymmetry results in encounter rates varying by 50--70% across the sky at the same distance. Finally, the LMC--mass satellite induces a density wake in the host's DM halo, further boosting the encounter rates near the LMC--analog. We also explore the influence of stream orbital properties, finding a 50% increase in encounters for streams moving retrograde to the LMC--analog's orbit in the opposite hemisphere. The dependence of encounter rates on stream location and orbit has important implications for where to search for new streams with spurs and gaps in the Milky Way. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 26 pages, 15 figures, submitted to APJ

Testing the near-far connection with FIRE simulations: inferring the stellar mass function of the proto-Local Group at z > 6 using the fossil record of present-day galaxies

Authors: Pratik J. Gandhi, Andrew Wetzel, Michael Boylan-Kolchin, Robyn E. Sanderson, Alessandro Savino, Daniel R. Weisz, Erik J. Tollerud, Guochao Sun, Claude-Andre Faucher-Giguere

Abstract: The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil r… ▽ More The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful `near-far' technique for studying the reionisation-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass (Mstar < 10^9 Msun) galaxies across 11 Milky Way- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterise their progenitors at z ~ 6 - 9, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with Mstar ~ 10^5 Msun (10^9 Msun) had ~1 (~30) progenitors at z ~ 7, and its main progenitor comprised ~100% (~50%) of the total stellar mass of all its progenitors at z ~ 7. We show that although only ~ 15% of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at z ~ 6 - 9. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just the low-mass galaxies at z = 0 and not the MW/M31 hosts correctly recovers the slope of the SMF down to Mstar ~ 10^4.5 Msun at z > 6. Thus, we validate the `near-far' approach as an unbiased tool for probing low-mass reionisation-era galaxies. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 6 figures, 3 appendix figures, and 1 table. Submitted to MNRAS

Modeling the orbital histories of satellites of Milky Way-mass galaxies: testing static host potentials against cosmological simulations

Authors: Isaiah B. Santistevan, Andrew Wetzel, Erik Tollerud, Robyn E Sanderson, Jorge Moreno, Ekta Patel

Abstract: Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modeling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmolog… ▽ More Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modeling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass halos. When a typical surviving satellite fell in ($3.4-9.7$ Gyr ago), the host halo mass and radius were typically $26-86$ per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $r\lesssim50$ kpc, opposite to dark-matter-only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at $z=0$ and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $1.6-4.7$ Gyr ago. Most orbital properties are minimally biased, $\lesssim10$ per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are $\approx20$ per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are $\approx40-80$ per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at $z=0$. △ Less

Submitted 8 December, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 24 pages, 12 figures, 3 appendices, 2 appendix figures. Accepted for publication in MNRAS

The proto-galaxy of Milky Way-mass haloes in the FIRE simulations

Authors: Danny Horta, Emily C. Cunningham, Robyn Sanderson, Kathryn V. Johnston, Alis Deason, Andrew Wetzel, Fiona McCluskey, Nicolás Garavito-Camargo, Lina Necib, Claude-André Faucher-Giguère, Arpit Arora, Pratik J. Gandhi

Abstract: Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological… ▽ More Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: i) can have a stellar mass range between $1\times10^{8}<\mathrm{M}_{\star}<2\times10^{10}[\mathrm{M}_{\odot}]$, a virial mass range between $3\times10^{10}<\mathrm{M}_{\star}<6\times10^{11}[\mathrm{M}_{\odot}]$, and be as young as $8 \lesssim \mathrm{Age} \lesssim 12.8$ [Gyr] ($1\lesssim z \lesssim 6$); ii) are predominantly centrally concentrated, with $\sim50\%$ of the stars contained within $5-10$ kpc; iii) on average show weak but systematic net rotation in the plane of the host's disc at $z=0$ (i.e., $0.25\lesssim\langleκ/κ_{\mathrm{disc}}\rangle\lesssim0.8$); iv) present [$α$/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host's old disc; v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that ~60% of the proto-Milky Way galaxies are comprised by 1 dominant system ($1/5\lesssim$M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim4/5$) and 4-5 lower mass systems (M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim1/10$); the other ~40% are comprised by 2 dominant systems and 3-4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems. △ Less

Submitted 11 December, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted in MNRAS. 16 pages, 12 figures and 3 tables

Synthetic Gaia DR3 surveys from the FIRE cosmological simulations of Milky-Way-mass galaxies

Authors: Tri Nguyen, Xiaowei Ou, Nondh Panithanpaisal, Nora Shipp, Lina Necib, Robyn Sanderson, Andrew Wetzel

Abstract: The third data release (DR3) of Gaia has provided a five-fold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky-Way-mass galaxies of the Latte suite of th… ▽ More The third data release (DR3) of Gaia has provided a five-fold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky-Way-mass galaxies of the Latte suite of the Fire-2 cosmological simulations. These synthetic surveys match the selection function, radial velocity measurements, and photometry of Gaia DR3, adapting the code base Ananke, previously used to match the Gaia DR2 release in Sanderson et al. 2020. The synthetic surveys are publicly available and can be found at http://ananke.hub.yt/. Similarly to the previous release of Ananke, these surveys are based on cosmological simulations and thus able to model non-equilibrium dynamical effects, making them a useful tool in testing and interpreting Gaia DR3. △ Less

Submitted 17 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 17 pages, 7 tables, 6 figures

The Galactic Center with Roman

Authors: Sean K. Terry, Matthew W. Hosek Jr., Jessica R. Lu, Casey Lam, Natasha Abrams, Arash Bahramian, Richard Barry, Jean-Phillipe Beaulieu, Aparna Bhattacharya, Devin Chu, Anna Ciurlo, Will Clarkson, Tuan Do, Kareem El-Badry, Ryan Felton, Matthew Freeman, Abhimat Gautam, Andrea Ghez, Daniel Huber, Jason Hunt, Macy Huston, Tharindu Jayasinghe, Naoki Koshimoto, Madeline Lucey, Florian Peißker , et al. (9 additional authors not shown)

Abstract: We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational… ▽ More We advocate for a Galactic center (GC) field to be added to the Galactic Bulge Time Domain Survey (GBTDS). The new field would yield high-cadence photometric and astrometric measurements of an unprecedented ${\sim}$3.3 million stars toward the GC. This would enable a wide range of science cases, such as finding star-compact object binaries that may ultimately merge as LISA-detectable gravitational wave sources, constraining the mass function of stars and compact objects in different environments, detecting populations of microlensing and transiting exoplanets, studying stellar flares and variability in young and old stars, and monitoring accretion onto the central supermassive black hole. In addition, high-precision proper motions and parallaxes would open a new window into the large-scale dynamics of stellar populations at the GC, yielding insights into the formation and evolution of galactic nuclei and their co-evolution with the growth of the supermassive black hole. We discuss the possible trade-offs between the notional GBTDS and the addition of a GC field with either an optimal or minimal cadence. Ultimately, the addition of a GC field to the GBTDS would dramatically increase the science return of Roman and provide a legacy dataset to study the mid-plane and innermost regions of our Galaxy. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 19 pages, 3 figures. Submitted to the NASA Roman Core Community Surveys White Paper Call

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

The far side of the Galactic bar/bulge revealed through semi-regular variables

Authors: Daniel R. Hey, Daniel Huber, Benjamin J. Shappee, Joss Bland-Hawthorn, Thor Tepper-García, Robyn Sanderson, Sukanya Chakrabarti, Nicholas Saunders, Jason A. S. Hunt, Timothy R. Bedding, John Tonry

Abstract: The Galactic bulge and bar are critical to our understanding of the Milky Way. However, due to the lack of reliable stellar distances, the structure and kinematics of the bulge/bar beyond the Galactic center have remained largely unexplored. Here, we present a method to measure distances of luminous red giants using a period-amplitude-luminosity relation anchored to the Large Magellanic Cloud, wit… ▽ More The Galactic bulge and bar are critical to our understanding of the Milky Way. However, due to the lack of reliable stellar distances, the structure and kinematics of the bulge/bar beyond the Galactic center have remained largely unexplored. Here, we present a method to measure distances of luminous red giants using a period-amplitude-luminosity relation anchored to the Large Magellanic Cloud, with random uncertainties of 10-15% and systematic errors below 1-2%. We apply this method to data from the Optical Gravitational Lensing Experiment (OGLE) to measure distances to $190,302$ stars in the Galactic bulge and beyond out to 20 kpc. Using this sample we measure a distance to the Galactic center of $R_0$ = $8108\pm106_{\rm stat}\pm93_{\rm sys}$ pc, consistent with astrometric monitoring of stars orbiting Sgr A*. We cross-match our distance catalog with Gaia DR3 and use the subset of $39,566$ overlapping stars to provide the first constraints on the Milky Way's velocity field ($V_R,V_φ,V_z$) beyond the Galactic center. We show that the $V_R$ quadrupole from the bar's near side is reflected with respect to the Galactic center, indicating that the bar is both bi-symmetric and aligned with the inner disk, and therefore dynamically settled along its full extent. We also find that the vertical height $V_Z$ map has no major structure in the region of the Galactic bulge, which is inconsistent with a current episode of bar buckling. Finally, we demonstrate with N-body simulations that distance uncertainty plays a major factor in the alignment of the major and kinematic axes of the bar and distribution of velocities, necessitating caution when interpreting results for distant stars. △ Less

Submitted 21 November, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted to the Astronomical Journal

The dark side of FIRE: predicting the population of dark matter subhaloes around Milky Way-mass galaxies

Authors: Megan Barry, Andrew Wetzel, Sierra Chapman, Jenna Samuel, Robyn Sanderson, Arpit Arora

Abstract: A variety of observational campaigns seek to test dark-matter models by measuring dark-matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within LambdaCDM, we examine 11 Milky Way (MW)-mass haloes fro… ▽ More A variety of observational campaigns seek to test dark-matter models by measuring dark-matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within LambdaCDM, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to 10^6 Msun, distances < 50 kpc of the galactic center, across z = 0 - 1 (lookback time 0 - 8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytic models. A typical MW-mass halo contains ~16 subhaloes >10^7 Msun (~1 subhalo >10^8 Msun) within 50 kpc at z = 0. We compare our results with dark-matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2-10x, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ~10x higher at z = 1 than at z = 0. Subhaloes have nearly isotropic orbital velocity distributions at z = 0. Across our simulations, we also identified 4 analogs of Large Magellanic Cloud satellite passages; these analogs enhance subhalo counts by 1.4-2.7 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ~5 per Gyr for a stream like GD-1, sufficient to make subhalo-stream interactions a promising method of measuring dark subhaloes. △ Less

Submitted 13 July, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: 14 pages, Published in MNRAS

Orientations of DM Halos in FIRE-2 Milky Way-mass Galaxies

Authors: Jay Baptista, Robyn Sanderson, Dan Huber, Andrew Wetzel, Omid Sameie, Michael Boylan-Kolchin, Jeremy Bailin, Philip F. Hopkins, Claude-Andre Faucher-Giguere, Sukanya Chakrabarti, Drona Vargya, Nondh Panithanpaisal, Arpit Arora, Emily Cunningham

Abstract: The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FI… ▽ More The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $σ= 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Submitted to ApJ, 19 pages, 12 figures, 3 tables

The observable properties of galaxy accretion events in Milky Way-like galaxies in the FIRE-2 cosmological simulations

Authors: Danny Horta, Emily C. Cunningham, Robyn E. Sanderson, Kathryn V. Johnston, Nondh Panithanpaisal, Arpit Arora, Lina Necib, Andrew Wetzel, Jeremy Bailin, Claude-André Faucher-Giguère

Abstract: In the $Λ$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations al… ▽ More In the $Λ$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M$_{*}$$\gtrsim$10$\times$$^{7}$M$_{\odot}$, enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that, the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [$α$/Fe]-[Fe/H] compositions of the stellar debris of such mergers $at$ $present$ $day$. These parameters also govern the resulting dynamical state of an accreted galaxy at $z=0$, leading to the expectation that the inner regions of the stellar halo (R$_{\mathrm{GC}}$ $\lesssim$30 kpc) should contain fully phase-mixed debris from both lower and higher mass satellites. In addition, we find that a significant fraction of the lower mass satellites accreted at early times deposit debris in the outer halo (R$_{\mathrm{GC}}$ $>$50 kpc) that are $not$ fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen $in$ $situ$ in higher redshift surveys. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: Submitted for publication in ApJ

Building an Acceleration Ladder with Tidal Streams and Pulsar Timing

Authors: Peter Craig, Sukanya Chakrabarti, Robyn E. Sanderson, Farnik Nikakhtar

Abstract: We analyze stellar streams in action-angle coordinates combined with recent local direct acceleration measurements to provide joint constraints on the potential of our Galaxy. Our stream analysis uses the Kullback-Leibler divergence with a likelihood analysis based on the two-point correlation function. We provide joint constraints from pulsar accelerations and stellar streams for local and global… ▽ More We analyze stellar streams in action-angle coordinates combined with recent local direct acceleration measurements to provide joint constraints on the potential of our Galaxy. Our stream analysis uses the Kullback-Leibler divergence with a likelihood analysis based on the two-point correlation function. We provide joint constraints from pulsar accelerations and stellar streams for local and global parameters that describe the potential of the Milky Way (MW). Our goal is to build an ``acceleration ladder", where direct acceleration measurements that are currently limited in dynamic range are combined with indirect techniques that can access a much larger volume of the MW. To constrain the MW potential with stellar streams, we consider the Palomar 5, Orphan, Nyx, Helmi and GD1 streams. Of the potential models that we have considered here, the preferred potential for the streams is a two-component Staeckel potential. We also compare the vertical accelerations from stellar streams and pulsar timing, defining a function $f(z) = α_{1pulsar}z - \frac{\partialΦ}{\partial z}$, where $Φ$ is the MW potential determined from stellar streams, and $α_{1~\rm pulsar}z$ is the vertical acceleration determined from pulsar timing observations. Our analysis indicates that the Oort limit determined from streams is consistently (regardless of the choice of potential) lower than that determined from pulsar timing observations. The calibration we have derived here may be used to correct the estimate of the acceleration from stellar streams. △ Less

Submitted 4 May, 2023; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: 8 pages, 4 figures, 1 table. Accepted by ApJ Letters

Constraining the Tilt of the Milky Way's Dark Matter Halo with the Sagittarius Stream

Authors: Nondh Panithanpaisal, Robyn E. Sanderson, Arpit Arora, Emily C. Cunningham, Jay Baptista

Abstract: Recent studies have suggested that the Milky Way (MW)'s Dark Matter (DM) halo may be significantly tilted with respect to its central stellar disk, a feature that might be linked to its formation history. In this work, we demonstrate a method of constraining the orientation of the minor axis of the DM halo using the angle and frequency variables. This method is complementary to other traditional t… ▽ More Recent studies have suggested that the Milky Way (MW)'s Dark Matter (DM) halo may be significantly tilted with respect to its central stellar disk, a feature that might be linked to its formation history. In this work, we demonstrate a method of constraining the orientation of the minor axis of the DM halo using the angle and frequency variables. This method is complementary to other traditional techniques, such as orbit fitting. We first test the method using a simulated tidal stream evolving in a realistic environment inside an MW-mass host from the FIRE cosmological simulation, showing that the theoretical description of a stream in the action-angle-frequency formalism still holds for a realistic dwarf galaxy stream in a cosmological potential. Utilizing the slopes of the line in angle and frequency space, we show that the correct rotation frame yields a minimal slope difference, allowing us to put a constraint on the minor axis location. Finally, we apply this method to the Sagittarius stream's leading arm. We report that the MW's DM halo is oblate with the flattening parameter in the potential $q\sim0.7-0.9$ and the minor axis pointing toward $(\ell,b) = (42^{o},48^{o})$. Our constraint on the minor axis location is weak and disagrees with the estimates from other works; we argue that the inconsistency can be attributed in part to the observational uncertainties and in part to the influence of the Large Magellanic Cloud. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: 16 pages, 12 figures

Orbital dynamics and histories of satellite galaxies around Milky Way-mass galaxies in the FIRE simulations

Authors: Isaiah B. Santistevan, Andrew Wetzel, Erik Tollerud, Robyn E. Sanderson, Jenna Samuel

Abstract: The orbits of satellite galaxies encode rich information about their histories. We investigate the orbital dynamics and histories of satellite galaxies around Milky Way (MW)-mass host galaxies using the FIRE-2 cosmological simulations, which, as previous works have shown, produce satellite mass functions and spatial distributions that broadly agree with observations. We first examine trends in orb… ▽ More The orbits of satellite galaxies encode rich information about their histories. We investigate the orbital dynamics and histories of satellite galaxies around Milky Way (MW)-mass host galaxies using the FIRE-2 cosmological simulations, which, as previous works have shown, produce satellite mass functions and spatial distributions that broadly agree with observations. We first examine trends in orbital dynamics at z = 0, including total velocity, specific angular momentum, and specific total energy: the time of infall into the MW-mass halo primarily determines these orbital properties. We then examine orbital histories, focusing on the lookback time of first infall into a host halo and pericenter distances, times, and counts. Roughly 37 per cent of galaxies with Mstar < 10^7 Msun were `pre-processed' as a satellite in a lower-mass group, typically ~2.7 Gyr before falling into the MW-mass halo. Half of all satellites at z = 0 experienced multiple pericenters about their MW-mass host. Remarkably, for most (67 per cent) of these satellites, their most recent pericenter was not their minimum pericenter: the minimum typically was ~40 per cent smaller and occurred ~6 Gyr earlier. These satellites with growing pericenters appear to have multiple origins: for about half, their specific angular momentum gradually increased over time, while for the other half, most rapidly increased near their first apocenter, suggesting that a combination of a time-dependent MW-mass halo potential and dynamical perturbations in the outer halo caused these satellites' pericenters to grow. Our results highlight the limitations of idealized, static orbit modeling, especially for pericenter histories. △ Less

Submitted 6 December, 2022; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 21 pages, 10 figues, 2 appendices, 4 appendix figures. Accepted in MNRAS for publication

Streams on FIRE: Populations of Detectable Stellar Streams in the Milky Way and FIRE

Authors: Nora Shipp, Nondh Panithanpaisal, Lina Necib, Robyn Sanderson, Denis Erkal, Ting S. Li, Isaiah B. Santistevan, Andrew Wetzel, Lara R. Cullinane, Alexander P. Ji, Sergey E. Koposov, Kyler Kuehn, Geraint F. Lewis, Andrew B. Pace, Daniel B. Zucker, Joss Bland-Hawthorn, Emily C. Cunningham, Stacy Y. Kim, Sophia Lilleengen, Jorge Moreno, Sanjib Sharma

Abstract: We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we pro… ▽ More We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to > 110 kpc) and surviving only at larger apocenters (> 40 kpc) than those observed in the Milky Way. We find that the population of high-stellar mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as satellites in DES-like observations, since their tidal tails are too low-surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low surface brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 24 pages, 13 figures, 3 tables, submitted to ApJ

On the stability of tidal streams in action space

Authors: Arpit Arora, Robyn E. Sanderson, Nondh Panithanpaisal, Emily C. Cunningham, Andrew Wetzel, Nicolás Garavito-Camargo

Abstract: In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological baryonic simulations… ▽ More In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological baryonic simulations. Using these potentials, we calculate actions for star particles in tidal streams around three galaxies with varying merger histories at each snapshot from 7 Gyr ago to the present day. We determine the action-space coherence preserved by each model using the Kullback-Leibler Divergence to gauge the degree of clustering in actions and the relative stability of the clusters over time. We find that all models produce a clustered action space for simulations with no significant mergers. However, a massive (mass ratio prior to infall more similar than 1:8) interacting galaxy not present in the model will result in mischaracterized orbits for stars most affected by the interaction. The locations of the action space clusters (i.e. the orbits of the stream stars) are only preserved by the time-evolving model, while the time-independent models can lose significant amounts of information as soon as 0.5--1 Gyr ago, even if the system does not undergo a significant merger. Our results imply that reverse-integration of stream orbits in the MW using a fixed potential is likely to give incorrect results if integrated longer than 0.5 Gyr into the past. △ Less

Submitted 26 October, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Journal ref: The Astrophysical Journal, vol. 939, no. 1, 2022, p. 2

The relationship between age, metallicity, and abundances for disk stars in a simulated Milky Way galaxy

Authors: Andreia Carrillo, Melissa K. Ness, Keith Hawkins, Robyn Sanderson, Kaile Wang, Andrew Wetzel, Matthew A. Bellardini

Abstract: Observations of the Milky Way's low-$α$ disk show that at fixed metallicity, [Fe/H], several element abundance, [X/Fe], correlate with age, with unique slopes and small scatters around the age-[X/Fe] relations. In this study, we turn to simulations to explore the age-[X/Fe] relations for the elements C, N, O, Mg, Si, S, and Ca that are traced in a FIRE-2 cosmological zoom-in simulation of a Milky… ▽ More Observations of the Milky Way's low-$α$ disk show that at fixed metallicity, [Fe/H], several element abundance, [X/Fe], correlate with age, with unique slopes and small scatters around the age-[X/Fe] relations. In this study, we turn to simulations to explore the age-[X/Fe] relations for the elements C, N, O, Mg, Si, S, and Ca that are traced in a FIRE-2 cosmological zoom-in simulation of a Milky Way-like galaxy, m12i, and understand what physical conditions give rise to the observed age-[X/Fe] trends. We first explore the distributions of mono-age populations in their birth and current locations, [Fe/H], and [X/Fe], and find evidence for inside-out radial growth for stars with ages < 7 Gyr. We then examine the age-[X/Fe] relations across m12i's disk and find that the direction of the trends agree with observations, apart from C, O, and Ca, with remarkably small intrinsic scatters, $σ_{int}$ (0.01-0.04 dex). This $σ_{int}$ measured in the simulations is also metallicity-dependent, with $σ_{int}$ $\approx$ 0.025 dex at [Fe/H]=-0.25 dex versus $σ_{int}$ $\approx$ 0.015 dex at [Fe/H]=0 dex, and a similar metallicity dependence is seen in the GALAH survey for the elements in common. Additionally, we find that $σ_{int}$ is higher in the inner galaxy, where stars are older and formed in less chemically-homogeneous environments. The age-[X/Fe] relations and the small scatter around them indicate that simulations capture similar chemical enrichment variance as observed in the Milky Way, arising from stars sharing similar element abundances at a given birth place and time. △ Less

Submitted 24 April, 2022; originally announced April 2022.

Comments: 27 pages, 13 figures. Submitted to ApJ

Kinematics and Metallicity of Red Giant Branch Stars in the Northeast Shelf of M31

Authors: Ivanna Escala, Karoline M. Gilbert, Mark Fardal, Puragra Guhathakurta, Robyn E. Sanderson, Jason S. Kalirai, Bahram Mobasher

Abstract: We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge… ▽ More We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge pattern agrees with expectations of a tidal shell formed in a radial merger and provides strong evidence in favor of predictions of Giant Stellar Stream (GSS) formation models in which the NE shelf originates from the second orbital wrap of the tidal debris. The observed concentric wedge patterns of the NE, West (W), and Southeast (SE) shelves corroborate this interpretation independently of the models. We do not detect a kinematical signature in the NE shelf region corresponding to an intact progenitor core, favoring GSS formation models in which the progenitor is completely disrupted. The shelf's photometric metallicity distribution implies that it is dominated by tidal material, as opposed to the phase-mixed stellar halo or the disk. The metallicity distribution ([Fe/H]$_{\rm phot}$ = $-0.42$ $\pm$ $0.01$) also matches the GSS, and consequently the W and SE shelves, further supporting a direct physical association between the tidal features. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 21 pages main text, 15 figures, 5 tables (including Appendix). Submitted to AJ

FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations

Authors: Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

Abstract: Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numeric… ▽ More Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

R2-D2: Roman and Rubin -- From Data to Discovery

Authors: Suvi Gezari, Misty Bentz, Kishalay De, K. Decker French, Aaron Meisner, Michelle Ntampaka, Robert Jedicke, Ekta Patel, Daniel Perley, Robyn Sanderson, Christian Aganze, Igor Andreoni, Eric F. Bell, Edo Berger, Ian Dell'Antonio, Ryan Foley, Henry Hsieh, Mansi Kasliwal, Joel Kastner, Charles D. Kilpatrick, J. Davy Kirkpatrick, Casey Lam, Karen Meech, Dante Minniti, Ethan O. Nadler , et al. (6 additional authors not shown)

Abstract: The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to… ▽ More The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation

Authors: Andrew Wetzel, Christopher C. Hayward, Robyn E. Sanderson, Xiangcheng Ma, Daniel Angles-Alcazar, Robert Feldmann, T. K Chan, Kareem El-Badry, Coral Wheeler, Shea Garrison-Kimmel, Farnik Nikakhtar, Nondh Panithanpaisal, Arpit Arora, Alexander B. Gurvich, Jenna Samuel, Omid Sameie, Viraj Pandya, Zachary Hafen, Cameron Hummels, Sarah Loebman, Michael Boylan-Kolchin, James S. Bullock, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert , et al. (1 additional authors not shown)

Abstract: We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, includ… ▽ More We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. △ Less

Submitted 29 March, 2023; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies… ▽ More This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

Applying the metallicity-dependent binary fraction to double white dwarf formation: Implications for LISA

Authors: Sarah Thiele, Katelyn Breivik, Robyn E. Sanderson, Rodrigo Luger

Abstract: Short-period double white dwarf (DWD) binaries will be the most prolific source of gravitational waves (GWs) for the Laser Interferometer Space Antenna (LISA). DWDs with GW frequencies below $\sim1$ mHz will be the dominant contributor to a stochastic foreground caused by overlapping GW signals. Population modeling of Galactic DWDs typically assumes a binary fraction of 50% and a log-uniform Zero… ▽ More Short-period double white dwarf (DWD) binaries will be the most prolific source of gravitational waves (GWs) for the Laser Interferometer Space Antenna (LISA). DWDs with GW frequencies below $\sim1$ mHz will be the dominant contributor to a stochastic foreground caused by overlapping GW signals. Population modeling of Galactic DWDs typically assumes a binary fraction of 50% and a log-uniform Zero Age Main Sequence (ZAMS) orbital period distribution. However, recent observations have shown that the binary fraction of close, solar-type stars exhibits a strong anti-correlation with metallicity which modulates the ZAMS orbital period distribution below $10^4$ days. In this study we perform the first simulation of the Galactic DWD population observable by LISA which incorporates an empirically-derived metallicity-dependent binary fraction, using the binary population synthesis suite COSMIC and a metallicity-dependent star formation history. We compare two models: one which assumes a metallicity-dependent binary fraction, and one with a binary fraction of 50%. We repeat our analysis for three different assumptions for Roche-lobe overflow interactions. We find that while metallicity impacts the evolution and intrinsic properties of our simulated DWD progenitor binaries, the LISA-resolvable populations of the two models remain roughly indistinguishable. However, the size of the total Galactic DWD population orbiting in the LISA frequency band is reduced by more than half when accounting for a metallicity-dependent binary fraction for two of our four variations, which also lowers the effective foreground. The LISA population remains unchanged in number for two variations, highlighting the sensitivity of the population to binary evolution prescriptions. △ Less

Submitted 16 December, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 24 pages, 11 figures. Accepted version for publication in ApJ

Gaia May Detect Hundreds of Well-characterised Stellar Black Holes

Authors: Chirag Chawla, Sourav Chatterjee, Katelyn Breivik, Chaithanya Krishna Moorthy, Jeff J. Andrews, Robyn E. Sanderson

Abstract: Detection of black holes (BHs) with detached luminous companions (LCs) can be instrumental in connecting the BH properties with their progenitors' since the latter can be inferred from the observable properties of the LC. Past studies showed the promise of Gaia astrometry in detecting BH-LC binaries. We build upon these studies by: 1) initialising the zero-age binary properties based on realistic,… ▽ More Detection of black holes (BHs) with detached luminous companions (LCs) can be instrumental in connecting the BH properties with their progenitors' since the latter can be inferred from the observable properties of the LC. Past studies showed the promise of Gaia astrometry in detecting BH-LC binaries. We build upon these studies by: 1) initialising the zero-age binary properties based on realistic, metallicity-dependent star-formation history in the Milky Way (MW), 2) evolving these binaries to current epoch to generate realistic MW populations of BH-LC binaries, 3) distributing these binaries in the MW preserving the complex age-metallicity-Galactic position correlations, 4) accounting for extinction and reddening using three-dimensional dust maps, 5) examining the extended Gaia mission's ability to resolve BH-LC binaries. We restrict ourselves to detached BH-LC binaries with orbital period <10 yr such that Gaia can observe at least one full orbit. We find: 1) the extended Gaia mission can astrometrically resolve 30-300 detached BH-LC binaries depending on our assumptions of supernova physics and astrometric detection threshold; 2) Gaia's astrometry alone can indicate BH candidates for 10-100 BH-LC binaries by constraining the dark primary mass >3 Msun; 3) distributions of observables including orbital periods, eccentricities, and component masses are sensitive to the adopted binary evolution model, hence can directly inform binary evolution models. Finally, we comment on the potential to further characterise these BH binaries through radial velocity measurements and observation of X-ray counterparts. △ Less

Submitted 17 March, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 26 pages, 18 figures, and 1 table; accepted for publication in The Astrophysical Journal

Reading the CARDs: the Imprint of Accretion History in the Chemical Abundances of the Milky Way's Stellar Halo

Authors: Emily C. Cunningham, Robyn E. Sanderson, Kathryn V. Johnston, Nondh Panithanpaisal, Melissa K. Ness, Andrew Wetzel, Sarah R. Loebman, Ivanna Escala, Danny Horta, Claude-André Faucher-Giguère

Abstract: In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of… ▽ More In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies' assembly histories by modelling the CARDs of the stellar halos of the Latte galaxies as a linear combination of template CARDs from disrupted dwarfs, with different stellar masses $M_{\star}$ and quenching times $t_{100}$. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of $<10\%$. For the fraction of mass accreted as a function of $t_{100}$, we find residuals of $20-30\%$ for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Though the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 32 pages, 20 figures. Submitted to ApJ; comments welcome!

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Authors: Erin Kado-Fong, Robyn E. Sanderson, Jenny E. Greene, Emily C. Cunningham, Coral Wheeler, T. K. Chan, Kareem El-Badry, Philip F. Hopkins, Andrew Wetzel, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Song Huang, Eliot Quataert, Tjitske Starkenburg

Abstract: Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosm… ▽ More Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies. △ Less

Submitted 2 May, 2022; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 28 pages, 15 figures. Accepted to ApJ

Galactic Mass Estimates using Dwarf Galaxies as Kinematic Tracers

Authors: Anika Slizewski, Xander Dufresne, Keslen Murdock, Gwendolyn Eadie, Robyn Sanderson, Andrew Wetzel, Mario Juric

Abstract: New mass estimates and cumulative mass profiles with Bayesian credible regions (c.r.) for the Milky Way (MW) are found using the Galactic Mass Estimator (GME) code and dwarf galaxy (DG) kinematic data from multiple sources. GME takes a hierarchical Bayesian approach to simultaneously estimate the true positions and velocities of the DGs, their velocity anisotropy, and the model parameters for the… ▽ More New mass estimates and cumulative mass profiles with Bayesian credible regions (c.r.) for the Milky Way (MW) are found using the Galactic Mass Estimator (GME) code and dwarf galaxy (DG) kinematic data from multiple sources. GME takes a hierarchical Bayesian approach to simultaneously estimate the true positions and velocities of the DGs, their velocity anisotropy, and the model parameters for the Galaxy's total gravitational potential. In this study, we incorporate meaningful prior information from past studies and simulations. The prior distributions for the physical model are informed by the results of Eadie & Juric (2019), which used globular clusters instead of DGs, as well as by the subhalo distributions of the Ananke Gaia-like surveys from Feedback In Realistic Environments-2 (Fire-2) cosmological simulations (see Sanderson et al. 2020). Using DGs beyond 45 kpc, we report median and 95% c.r estimates for $r_{200}$ = 212.8 (191.12,238.44) kpc, and for the total enclosed mass $M_{200}$ = 1.19 (0.87,1.68)$\times10^{12}M_{\odot}$ (adopting $Δ_c=200$). Median mass estimates at specific radii are also reported (e.g., $M(<50\text{ kpc})=0.52\times10^{12}M_{\odot}$ and $M(100\text{ kpc})=0.78\times10^{12}M_{\odot}$). Estimates are comparable to other recent studies using GAIA DR2 and DGs, but notably different from the estimates of Eadie & Juric (2019). We perform a sensitivity analysis to investigate whether individual DGs and/or a more massive Large Magellanic Cloud (LMC) on the order of $10^{11}M_{\odot}$ may be affecting our mass estimates. We find possible supporting evidence for the idea that some DGs are affected by a massive LMC and are not in equilibrium with the MW. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ, 16 pages, 12 figures, 4 tables

Orbital phase-driven biases in Galactic mass constraints from stellar streams

Authors: Stella Reino, Robyn E. Sanderson, Nondh Panithanpaisal, Elena M. Rossi, Konrad Kuijken

Abstract: One of the most promising tracers of the Galactic potential in the halo region are stellar streams. However, individual stream fits can be limited by systematic biases. To study these individual stream systematics, we fit streams in Milky Way-like galaxies from FIRE cosmological galaxy formation simulations with an analytic gravitational potential by maximizing the clustering of stream stars in ac… ▽ More One of the most promising tracers of the Galactic potential in the halo region are stellar streams. However, individual stream fits can be limited by systematic biases. To study these individual stream systematics, we fit streams in Milky Way-like galaxies from FIRE cosmological galaxy formation simulations with an analytic gravitational potential by maximizing the clustering of stream stars in action space. We show that for coherent streams the quality of the constraints depends on the orbital phase of the observed stream stars, despite the fact that the phase information is discarded in action-clustering methods. Streams on intermediate phases give the most accurate results, whereas pericentre streams can be highly biased. This behaviour is tied to the amount of correlation present between positions and momenta in each stream's data: weak correlation in pericentre streams prohibits efficient differentiation between potentials, while strong correlation in intermediate streams promotes it. Although simultaneous fitting of multiple streams is generally prescribed as the remedy to combat individual stream biases, we find that combining multiple pericentric streams is not enough to yield a bias-free result. We finally show that adopting the two-component Stäckel model does not fundamentally induce a biased mass estimate. With our full data set of two multi-wrap streams, we recovered the true rotation curve of the simulated galaxy within $12\%$ over the entire range of radii covered by our set of stars (10 - 176 kpc) and within $6.5\%$ between the 5 and 95-percentile distance range (23 - 109 kpc). △ Less

Submitted 29 October, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: 17 pages, 13 figures. Accepted for publication in MNRAS

Shapes of Milky-Way-Mass Galaxies with Self-Interacting Dark Matter

Authors: Drona Vargya, Robyn Sanderson, Omid Sameie, Michael Boylan-Kolchin, Philip F. Hopkins, Andrew Wetzel, Andrew Graus

Abstract: Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environment… ▽ More Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section $σ/m$. We compare the shape of the main dark matter (DM) halo at redshift $z=0$ predicted by SIDM simulations (at $σ/m=0.1$, $1$, and $10$ cm$^2$ g$^{-1}$) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on $σ/m$. This transition could potentially be used to set limits on the SIDM cross-section in the MW. △ Less

Submitted 20 July, 2022; v1 submitted 28 April, 2021; originally announced April 2021.

Comments: 17 pages, 14 figures, 2 tables, accepted for publication in MNRAS

The Galaxy Progenitors of Stellar Streams around Milky Way-mass Galaxies in the FIRE Cosmological Simulations

Authors: Nondh Panithanpaisal, Robyn E. Sanderson, Andrew Wetzel, Emily C. Cunningham, Jeremy Bailin, Claude-André Faucher-Giguère

Abstract: Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at $z=0$, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: https://flathub.flatironinstitute.org/sapfire). In total, we identify… ▽ More Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at $z=0$, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: https://flathub.flatironinstitute.org/sapfire). In total, we identify 106 simulated stellar streams, with no significant differences in the number of streams and masses of their progenitors between the isolated and paired environments. We resolve simulated streams with stellar masses ranging from $\sim 5\times10^5$ up to $\sim 10^{9} M_\odot$, similar to the mass range between the Orphan and Sagittarius streams in the MW. We confirm that present-day simulated satellite galaxies are good proxies for stellar stream progenitors, with similar properties including their stellar mass function, velocity dispersion, [Fe/H] and [$α$/H] evolution tracks, and orbital distribution with respect to the galactic disk plane. Each progenitor's lifetime is marked by several important timescales: its infall, star-formation quenching, and stream-formation times. We show that the ordering of these timescales is different between progenitors with stellar masses higher and lower than $\sim 2\times10^6 M_\odot$. Finally, we show that the main factor controlling the rate of phase-mixing, and therefore fading, of tidal streams from satellite galaxies in MW-mass hosts is non-adiabatic evolution of the host potential. Other factors commonly used to predict phase-mixing timescales, such as progenitor mass and orbital circularity, show virtually no correlation with the number of dynamical times required for a stream to become phase-mixed. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 26 pages, 19 figures, public dataset available at https://flathub.flatironinstitute.org/sapfire

New families in our Solar neighborhood: applying Gaussian Mixture models for objective classification of structures in the Milky Way and in simulations

Authors: Farnik Nikakhtar, Robyn E. Sanderson, Andrew Wetzel, Sarah Loebman, Sanjib Sharma, Rachael Beaton, J. Ted Mackereth, Vijith Jacob Poovelil, Gail Zasowski, Ana Bonaca, Sarah Martell, Henrik Jonsson, Claude-Andre Faucher-Giguere

Abstract: The standard picture of galaxy formation motivates the decomposition of the Milky Way into 3--4 stellar populations with distinct kinematic and elemental abundance distributions: the thin disk, thick disk, bulge, and stellar halo. To test this idea, we construct a Gaussian mixture model (GMM) for both simulated and observed stars in the Solar neighborhood, using measured velocities and iron abunda… ▽ More The standard picture of galaxy formation motivates the decomposition of the Milky Way into 3--4 stellar populations with distinct kinematic and elemental abundance distributions: the thin disk, thick disk, bulge, and stellar halo. To test this idea, we construct a Gaussian mixture model (GMM) for both simulated and observed stars in the Solar neighborhood, using measured velocities and iron abundances (i.e., an augmented Toomre diagram) as the distributions to be decomposed. We compare results for the Gaia-APOGEE DR16 crossmatch catalog of the Solar neighborhood with those from a suite of synthetic Gaia-APOGEE crossmatches constructed from FIRE-2 cosmological simulations of Milky Way-mass galaxies. We find that in both the synthetic and real data, the best-fit GMM uses five independent components, some of whose properties resemble the standard populations predicted by galaxy formation theory. Two components can be identified unambiguously as the thin disk and another as the halo. However, instead of a single counterpart to the thick disk, there are three intermediate components with different age and alpha abundance distributions (although these data are not used to construct the model). We use decompositions of the synthetic data to show that the classified components indeed correspond to stars with different origins. By analogy with the simulated data, we show that our mixture model of the real Gaia-APOGEE crossmatch distinguishes the following components: (1) a classic thin disk of young stars on circular orbits (46%), (2) thin disk stars heated by interactions with satellites (22%), (3, 4) two components representing the velocity asymmetry of the alpha-enhanced thick disk (27%), and (5) a stellar halo consistent with early, massive accretion (4%). △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: 20 pages, 14 figures, 1 table

Using Action Space Clustering to Constrain the Accretion History of Milky Way like Galaxies

Authors: Youjia Wu, Monica Valluri, Nondh Panithanpaisal, Robyn E. Sanderson, Katherine Freese, Andrew Wetzel, Sanjib Sharma

Abstract: In the currently favored cosmological paradigm galaxies form hierarchically through the accretion of numerous satellite galaxies. Since the satellites are much less massive than the host halo, they occupy a small fraction of the volume in action space defined by the potential of the host halo. Since actions are conserved when the potential of the host halo changes adiabatically, stars from an accr… ▽ More In the currently favored cosmological paradigm galaxies form hierarchically through the accretion of numerous satellite galaxies. Since the satellites are much less massive than the host halo, they occupy a small fraction of the volume in action space defined by the potential of the host halo. Since actions are conserved when the potential of the host halo changes adiabatically, stars from an accreted satellite are expected to remain clustered in action space as the host halo evolves. In this paper, we identify accreted satellites in three Milky Way like disk galaxies from the cosmological baryonic FIRE-2 simulations by tracking satellite galaxies through simulation snapshots. We then try to recover these satellites by applying the cluster analysis algorithm Enlink to the orbital actions of accreted star particles in the present-day snapshot. We define several metrics to quantify the success of the clustering algorithm and use these metrics to identify well-recovered and poorly-recovered satellites. We plot these satellites in the infall time-progenitor mass (or stellar mass) space, and determine the boundaries between the well-recovered and poorly-recovered satellites in these two spaces with classification tree method. The groups found by Enlink are more likely to correspond to a real satellite if they have high significance, a quantity assigned by Enlink. Since cosmological simulations predict that most stellar halos have a population of insitu stars, we test the ability of Enlink to recover satellites when the sample is contaminated by 10-50% of insitu star particles, and show that most of the satellites well-recovered by Enlink in the absence of insitu stars, stay well-recovered even with 50% contamination. We thus expect that, in the future, cluster analysis in action space will be useful in upcoming data sets (e.g. Gaia) for identifying accreted satellites in the Milky Way. △ Less

Submitted 3 May, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: 16 pages, 9 figures, fixed some typos and added two references

Report number: NORDITA-2020-094, UTTG-14-2020, LCTP-21-10

The central densities of Milky Way-mass galaxies in cold and self-interacting dark matter models

Authors: Omid Sameie, Michael Boylan-Kolchin, Robyn Sanderson, Drona Vargya, Philip Hopkins, Andrew Wetzel, James Bullock, Andrew Graus, Victor Robles

Abstract: We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the ``Feedback In Realistic Environment'' (FIRE) project. The three simulated haloes have virial masses of $\sim 10^{12}\, \text{M}_\odot$ at $z=0$, and we study velocity-independent self-interaction cross sections of 1 and 10 ${\rm cm^2 \, g^{-1}}$. We study star formation rates an… ▽ More We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the ``Feedback In Realistic Environment'' (FIRE) project. The three simulated haloes have virial masses of $\sim 10^{12}\, \text{M}_\odot$ at $z=0$, and we study velocity-independent self-interaction cross sections of 1 and 10 ${\rm cm^2 \, g^{-1}}$. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at $z\leq1$, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts. △ Less

Submitted 8 September, 2021; v1 submitted 24 February, 2021; originally announced February 2021.

Comments: 10 pages, 5 figures. Accepted for publication in MNRAS

The origin of metal-poor stars on prograde disk orbits in FIRE simulations of Milky Way-mass galaxies

Authors: Isaiah B. Santistevan, Andrew Wetzel, Robyn E. Sanderson, Kareem El-Badry, Jenna Samuel, Claude-André Faucher-Giguère

Abstract: In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower-mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disk. However, recent observations of stars in the MW show that metal-poor ([Fe/H] < -2) stars are prefere… ▽ More In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower-mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disk. However, recent observations of stars in the MW show that metal-poor ([Fe/H] < -2) stars are preferentially on prograde orbits with respect to the disk. Using the FIRE-2 suite of cosmological zoom-in simulations of MW/M31-mass galaxies, we investigate the prevalence and origin of prograde metal-poor stars. Almost all (11 of 12) of our simulations have metal-poor stars on preferentially prograde orbits today and throughout most of their history: we thus predict that this is a generic feature of MW/M31-mass galaxies. The typical prograde-to-retrograde ratio is ~2:1, which depends weakly on stellar metallicity at Fe/H] < -1. These trends predicted by our simulations agree well with MW observations. Prograde metal-poor stars originate largely from a single LMC/SMC-mass gas-rich merger 7 - 12.5 Gyr ago, which deposited existing metal-poor stars and significant gas on an orbital vector that sparked the formation of and/or shaped the orientation of a long-lived stellar disk, giving rise to a prograde bias for all low-metallicity stars. We find sub-dominant contributions from in-situ stars formed in the host galaxy before this merger, and in some cases, additional massive mergers. We find few clear correlations between any properties of our MW/M31-mass galaxies at z = 0 and the degree of this prograde bias as a result of diverse merger scenarios. △ Less

Submitted 7 May, 2021; v1 submitted 5 February, 2021; originally announced February 2021.

Comments: 18 pages, 8 figures. Accepted for publication in MNRAS

Galactic potential constraints from clustering in action space of combined stellar stream data

Authors: Stella Reino, Elena M. Rossi, Robyn E. Sanderson, Elena Sellentin, Amina Helmi, Helmer H. Koppelman, Sanjib Sharma

Abstract: Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational poten… ▽ More Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighboring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Pal 5, Orphan and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a Stäckel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1, Pal 5 and Orphan streams whose combined data set yields $M(< 20\ \mathrm{kpc}) = 2.96^{+0.25}_{-0.26} \times 10^{11} \ M_{\odot}$. When including the Helmi stream in the data set, the mass uncertainty increases to $M(< 20\ \mathrm{kpc}) = 3.12^{+3.21}_{-0.46} \times 10^{11} \ M_{\odot}$. △ Less

Submitted 2 February, 2021; v1 submitted 1 July, 2020; originally announced July 2020.

Comments: 25 pages, 15 figures. Accepted for publication in MNRAS

The Milky Way's Shell Structure Reveals the Time of a Radial Collision

Authors: Thomas Donlon II, Heidi Jo Newberg, Robyn Sanderson, Lawrence M. Widrow

Abstract: We identify shell structures in the Milky Way for the first time. We find 2 shells in the Virgo Overdensity (VOD) region and 2 shells in the Hercules Aquila Cloud (HAC) region using Sloan Digital Sky Survey, Gaia, and LAMOST data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenit… ▽ More We identify shell structures in the Milky Way for the first time. We find 2 shells in the Virgo Overdensity (VOD) region and 2 shells in the Hercules Aquila Cloud (HAC) region using Sloan Digital Sky Survey, Gaia, and LAMOST data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenitor dwarf galaxy passed through the Galactic center 2.7 +/- 0.2 Gyr ago. Based on the time of collision, it is also possible that the VRM is related to the phenomenon that created phase-space spirals in the vertical motion of the disk and/or the Splash, and could have caused a burst of star formation in the inner disk. We analyze phase mixing in a collection of radial merger N-body simulations, and find that shell structure similar to that observed in Milky Way data disappears by 5 Gyr after collision with the Galactic center. The method used to calculate the merger time of the VRM was able to reliably recover the correct merger times for these simulations. Previous work supports the idea that the VRM and the Gaia Sausage/Gaia-Enceladus Merger are the same. However, the Gaia Sausage is widely believed to be 8--11 Gyr old. The disparate ages could be reconciled if the larger age is associated with an infall time when the progenitor crossed the virial radius; we do not constrain the time at which the progenitor became bound to the Milky Way. Alternatively, the Gaia Sausage could be younger than previously thought. △ Less

Submitted 9 September, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics

Authors: Emily C. Cunningham, Nicolas Garavito-Camargo, Alis J. Deason, Kathryn V. Johnston, Denis Erkal, Chervin F. P. Laporte, Gurtina Besla, Rodrigo Luger, Robyn E. Sanderson

Abstract: The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter G… ▽ More The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 25 pages, 13 figures. Accepted for publication in ApJ

Beyond Gaia: Asteroseismic Distances of M giants using Ground-Based Transient Surveys

Authors: Connor Auge, Daniel Huber, Aren Heinze, B. J. Shappee, John Tonry, Sukanya Chakrabarti, Robyn E. Sanderson, Larry Denneau, Heather Flewelling, Thomas W. -S. Holoien, C. S. Kochanek, Giuliano Pignata, Amanda Sickafoose, Brian Stalder, K. Z. Stanek, Dennis Stello, Todd A. Thompson

Abstract: Evolved stars near the tip of the red giant branch (TRGB) show solar-like oscillations with periods spanning hours to months and amplitudes ranging from $\sim$1 mmag to $\sim$100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently acce… ▽ More Evolved stars near the tip of the red giant branch (TRGB) show solar-like oscillations with periods spanning hours to months and amplitudes ranging from $\sim$1 mmag to $\sim$100 mmag. The systematic detection of the resulting photometric variations with ground-based telescopes would enable the application of asteroseismology to a much larger and more distant sample of stars than is currently accessible with space-based telescopes such as \textit{Kepler} or the ongoing Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. We present an asteroseismic analysis of 493 M giants using data from two ground-based surveys: the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the All-Sky Automated Survey for Supernovae (ASAS-SN). By comparing the extracted frequencies with constraints from \textit{Kepler}, the Sloan Digital Sky Survey Apache Point Observatory Galaxy Evolution Experiment (APOGEE), and Gaia we demonstrate that ground-based transient surveys allow accurate distance measurements to oscillating M giants with a precision of $\sim$15$\%$. Using stellar population synthesis models we predict that ATLAS and ASAS-SN can provide asteroseismic distances to $\sim$2$\times$10$^{6}$ galactic M giants out to typical distances of $20-50 \; \rm{kpc}$, vastly improving the reach of Gaia and providing critical constraints for Galactic archaeology and galactic dynamics. △ Less

Submitted 8 May, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 13 pages, 8 figures, 2 tables, accepted for publication in The Astronomical Journal

Stars made in outflows may populate the stellar halo of the Milky Way

Authors: Sijie Yu, James S. Bullock, Andrew Wetzel, Robyn E. Sanderson, Andrew S. Graus, Michael Boylan-Kolchin, Anna M. Nierenberg, Michael Y. Grudić, Philip F. Hopkins, Dušan Kereš, Claude-André Faucher-Giguère

Abstract: We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The… ▽ More We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: 19 pages, 20 figures, submitted to MNRAS

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Authors: Lina Necib, Bryan Ostdiek, Mariangela Lisanti, Timothy Cohen, Marat Freytsis, Shea Garrison-Kimmel, Philip F. Hopkins, Andrew Wetzel, Robyn Sanderson

Abstract: Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify… ▽ More Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor. △ Less

Submitted 28 January, 2022; v1 submitted 16 July, 2019; originally announced July 2019.

Comments: 8+5 pages, 4+7 figures. Nyx star sample updated with a modified probability calculation, conclusions unchanged. Stellar IDs can be found here https://linoush.github.io/nyx_stream/

Journal ref: Nature Astronomy volume 4, pages1078-1083 (2020)

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

Authors: Bryan Ostdiek, Lina Necib, Timothy Cohen, Marat Freytsis, Mariangela Lisanti, Shea Garrison-Kimmel, Andrew Wetzel, Robyn E. Sanderson, Philip F. Hopkins

Abstract: The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significan… ▽ More The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution. △ Less

Submitted 22 April, 2020; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: v1: 24 pages and 13 Figures + 4 appendices, v2: Journal version, fixed bug resulting in incorrect final catalog scores, updated plots, conclusions unchanged

Journal ref: A&A 636, A75 (2020)