Showing 1–50 of 135 results for author: Price-Whelan, A

Stellar Streams in the Gaia Era

Authors: Ana Bonaca, Adrian M. Price-Whelan

Abstract: The hierarchical model of galaxy formation predicts that the Milky Way halo is populated by tidal debris of dwarf galaxies and globular clusters. Due to long dynamical times, debris from the lowest mass objects remains coherent as thin and dynamically cold stellar streams for billions of years. The Gaia mission, providing astrometry and spectrophotometry for billions of stars, has brought three fu… ▽ More The hierarchical model of galaxy formation predicts that the Milky Way halo is populated by tidal debris of dwarf galaxies and globular clusters. Due to long dynamical times, debris from the lowest mass objects remains coherent as thin and dynamically cold stellar streams for billions of years. The Gaia mission, providing astrometry and spectrophotometry for billions of stars, has brought three fundamental changes to our view of stellar streams in the Milky Way. First, more than a hundred stellar streams have been discovered and characterized using Gaia data. This is an order of magnitude increase in the number of known streams, thanks to Gaia's capacity for identifying comoving groups of stars among the field Milky Way population. Second, Gaia data have revealed that density variations both along and across stellar streams are common. Dark-matter subhalos, as well as baryonic structures were theoretically predicted to form such features, but observational evidence for density variations was uncertain before Gaia. Third, stream kinematics are now widely available and have constrained the streams' orbits and origins. Gaia has not only provided proper motions directly, but also enabled efficient spectroscopic follow-up of the proper-motion selected targets. These discoveries have established stellar streams as a dense web of sensitive gravitational tracers in the Milky Way halo. We expect the coming decade to bring a full mapping of the Galactic population of stellar streams, as well as develop numerical models that accurately capture their evolution within the Milky Way for a variety of cosmological models. Perhaps most excitingly, the comparison between the two will be able to reveal the presence of dark-matter subhalos below the threshold for galaxy formation (~10^6 Msun), and provide the most stringent test of the cold dark matter paradigm on small scales. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: Review article for New Astronomy Reviews (as part of a Special Issue: "Gaia, the first crop of discoveries"). Comments are welcome!

[X/Fe] Marks the Spot: Mapping Chemical Azimuthal Variations in the Galactic Disk with APOGEE

Authors: Zoe Hackshaw, Keith Hawkins, Carrie Filion, Danny Horta, Chervin F. P. Laporte, Chris Carr, Adrian M. Price-Whelan

Abstract: Chemical cartography of the Galactic disk provides insights to its structure and assembly history over cosmic time. In this work, we use chemical cartography to explore chemical gradients and azimuthal substructure in the Milky Way disk with giant stars from APOGEE DR17. We confirm the existence of a radial metallicity gradient in the disk of $Δ$[Fe/H]/$Δ$R $\sim -0.066 \pm 0.0004$ dex/kpc and a v… ▽ More Chemical cartography of the Galactic disk provides insights to its structure and assembly history over cosmic time. In this work, we use chemical cartography to explore chemical gradients and azimuthal substructure in the Milky Way disk with giant stars from APOGEE DR17. We confirm the existence of a radial metallicity gradient in the disk of $Δ$[Fe/H]/$Δ$R $\sim -0.066 \pm 0.0004$ dex/kpc and a vertical metallicity gradient of $Δ$[Fe/H]/$Δ$Z $\sim -0.164 \pm 0.001$ dex/kpc. We find azimuthal variations ($\pm0.1$ dex) on top of the radial metallicity gradient that have been previously established with other surveys. The APOGEE giants show strong correlations with stellar age and the intensity of azimuthal variations in iron; older stellar populations show the largest deviations from the radial metallicity gradient. Beyond iron, we show that other elements (e.g., Mg, O) display azimuthal variations at the $\pm0.05$ dex-level across the Galactic disk. We illustrate that moving into the orbit-space could help constrain the mechanisms producing these azimuthal metallicity variations. These results suggest that the spiral arms of the Galaxy are not solely responsible for azimuthal metallicity variations and other Galactic processes are at play. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 19 pages, 11 figures

Gaia DR3 detectability of unresolved binary systems

Authors: Alfred Castro-Ginard, Zephyr Penoyre, Andrew R. Casey, Anthony G. A. Brown, Vasily Belokurov, Tristan Cantat-Gaudin, Ronald Drimmel, Morgan Fouesneau, Shourya Khanna, Evgeny P. Kurbatov, Adrian M. Price-Whelan, Hans-Walter Rix, Richard L. Smart

Abstract: Gaia can not individually resolve very close binary systems, however, the collected data can still be used to identify them. A powerful indicator of stellar multiplicity is the sources reported Renormalized Unit Weight Error (ruwe), which effectively captures the astrometric deviations from single-source solutions. We aim to characterise the imprints left on ruwe caused by binarity. By flagging po… ▽ More Gaia can not individually resolve very close binary systems, however, the collected data can still be used to identify them. A powerful indicator of stellar multiplicity is the sources reported Renormalized Unit Weight Error (ruwe), which effectively captures the astrometric deviations from single-source solutions. We aim to characterise the imprints left on ruwe caused by binarity. By flagging potential binary systems based on ruwe, we aim to characterise which of their properties will contribute the most to their detectability. We develop a model to estimate ruwe values for observations of Gaia sources, based on the biases to the single-source astrometric track arising from the presence of an unseen companion. Then, using the recipes from previous GaiaUnlimited selection functions, we estimate the selection probability of sources with high ruwe, and discuss what binary properties contribute to increasing the sources ruwe. We compute the maximum ruwe value which is compatible with single-source solutions as a function of their location on-sky. We see that binary systems selected as sources with a ruwe higher than this sky-varying threshold have a strong detectability window in their orbital period distribution, which peaks at periods equal to the Gaia observation time baseline. We demonstrate how our sky-varying ruwe threshold provides a more complete sample of binary systems when compared to single sky-averaged values by studying the unresolved binary population in the Gaia Catalogue of Nearby Stars. We provide the code and tools used in this study, as well as the sky-varying ruwe threshold through the GaiaUnlimited Python package △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 11 pages, 10 figures. Accepted for publication in A&A

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Authors: Gaia Collaboration, P. Panuzzo, T. Mazeh, F. Arenou, B. Holl, E. Caffau, A. Jorissen, C. Babusiaux, P. Gavras, J. Sahlmann, U. Bastian, Ł. Wyrzykowski, L. Eyer, N. Leclerc, N. Bauchet, A. Bombrun, N. Mowlavi, G. M. Seabroke, D. Teyssier, E. Balbinot, A. Helmi, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne , et al. (390 additional authors not shown)

Abstract: Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is exp… ▽ More Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way. △ Less

Submitted 19 April, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

zoomies: A tool to infer stellar age from vertical action in Gaia data

Authors: Sheila Sagear, Adrian M. Price-Whelan, Sarah Ballard, Yuxi, Lu, Ruth Angus, David W. Hogg

Abstract: Stellar age measurements are fundamental to understanding a wide range of astronomical processes, including galactic dynamics, stellar evolution, and planetary system formation. However, extracting age information from Main Sequence stars is complicated, with techniques often relying on age proxies in the absence of direct measurements. The Gaia data releases have enabled detailed studies of the d… ▽ More Stellar age measurements are fundamental to understanding a wide range of astronomical processes, including galactic dynamics, stellar evolution, and planetary system formation. However, extracting age information from Main Sequence stars is complicated, with techniques often relying on age proxies in the absence of direct measurements. The Gaia data releases have enabled detailed studies of the dynamical properties of stars within the Milky Way, offering new opportunities to understand the relationship between stellar age and dynamics. In this study, we leverage high-precision astrometric data from Gaia DR3 to construct a stellar age prediction model based only on stellar dynamical properties; namely, the vertical action. We calibrate two distinct, hierarchical stellar age--vertical action relations, first employing asteroseismic ages for red giant branch stars, then isochrone ages for main-sequence turn-off stars. We describe a framework called "zoomies" based on this calibration, by which we can infer ages for any star given its vertical action. This tool is open-source and intended for community use. We compare dynamical age estimates from "zoomies" with ages derived from other techniques for a sample of open clusters and main-sequence stars with asteroseismic age measurements. We also compare dynamical age estimates for stellar samples from the Kepler, K2, and TESS exoplanet transit surveys. While dynamical age relations are associated with large uncertainty, they are generally mass-independent and depend on homogeneously measured astrometric data. These age predictions are uniquely useful for large-scale demographic investigations, especially in disentangling the relationship between planet occurrence, metallicity, and age for low-mass stars. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 23 pages, 12 figures, 2 tables. Submitted to AAS Journals

Stellar Mergers or Truly Young? Intermediate-Age Stars on Highly-Radial Orbits in the Milky Way's Stellar Halo

Authors: Danny Horta, Yuxi, Lu, Melissa K. Ness, Mariangela Lisanti, Adrian M. Price-Whelan

Abstract: Reconstructing the mass assembly history of the Milky Way relies on obtaining detailed measurements of the properties of many stars in the Galaxy, especially in the stellar halo. One of the most constraining quantities is stellar age, as it can shed light on the accretion time and quenching of star formation in merging satellites. However, obtaining reliable age estimates for large samples of halo… ▽ More Reconstructing the mass assembly history of the Milky Way relies on obtaining detailed measurements of the properties of many stars in the Galaxy, especially in the stellar halo. One of the most constraining quantities is stellar age, as it can shed light on the accretion time and quenching of star formation in merging satellites. However, obtaining reliable age estimates for large samples of halo stars is difficult. We report published ages of 120 subgiant halo stars with highly-radial orbits that likely belong to the debris of the \textsl{Gaia-Enceladus/Sausage}~(GES) galaxy. The majority of these halo stars are old, with an age distribution characterized by a median of 11.6~Gyr and 16$^{\rm th}$(84$^{\rm th}$) percentile of 10.5~(12.7)~Gyr. However, the distribution is skewed, with a tail of younger stars that span ages down to $\sim6$--$9$ Gyr. All highly-radial halo stars have chemical and kinematic/orbital quantities that associate them with the GES debris. Initial results suggest that these intermediate-age stars are not a product of mass transfer and/or stellar mergers, which can bias their age determination low. If this conclusion is upheld by upcoming spectro-photometric studies, then the presence of these stars will pose an important challenge for constraining the properties of the GES merger and the accretion history of the Galaxy. △ Less

Submitted 10 June, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: Accepted for publication to ApJ after minor edits to the manuscript. 18 pages, 20 figures, 1 Table

Radial phase spirals in the Solar neighbourhood

Authors: Jason A. S. Hunt, Adrian M. Price-Whelan, Kathryn V. Johnston, Rachel L. McClure, Carrie Filion, Ben Cassese, Danny Horta

Abstract: The second data release of ESA's $Gaia$ mission revealed numerous signatures of disequilibrium in the Milky Way's disc. These signatures are seen in the planar kinematics of stars, which manifest as ridges and ripples in $R-v_φ$, and in vertical kinematics, where a prominent spiral is seen in the $z-v_z$ phase space. In this work, we show an equivalent $ΔR-v_{\mathrm{R}}$ phase spiral forms follow… ▽ More The second data release of ESA's $Gaia$ mission revealed numerous signatures of disequilibrium in the Milky Way's disc. These signatures are seen in the planar kinematics of stars, which manifest as ridges and ripples in $R-v_φ$, and in vertical kinematics, where a prominent spiral is seen in the $z-v_z$ phase space. In this work, we show an equivalent $ΔR-v_{\mathrm{R}}$ phase spiral forms following a perturbation to the disc. We demonstrate the behaviour of the $ΔR-v_{\mathrm{R}}$ phase spirals in both a toy model and a high resolution $N$-body simulation of a satellite interaction. We then confront these models with the data, where we find partial $ΔR-v_{\mathrm{R}}$ phase spirals in the Solar neighborhood using the most recent data from $Gaia$ DR3. This structure indicates ongoing radial phase mixing in the Galactic disc, suggesting a history of recent perturbations, either through internal or external (e.g., satellite) processes. Future work modelling the $z-v_z$ and $ΔR-v_{\mathrm{R}}$ phase spirals in tandem may help break degeneracy's between possible origins of the perturbation. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 13 pages, 10 figures, published in MNRAS

Journal ref: MNRAS, 2024, 527, 11393

Data-driven Dynamics with Orbital Torus Imaging: A Flexible Model of the Vertical Phase Space of the Galaxy

Authors: Adrian M. Price-Whelan, Jason A. S. Hunt, Danny Horta, Micah Oeur, David W. Hogg, Kathryn V. Johnston, Lawrence Widrow

Abstract: The vertical kinematics of stars near the Sun can be used to measure the total mass distribution near the Galactic disk and to study out-of-equilibrium dynamics. With contemporary stellar surveys, the tracers of vertical dynamics are so numerous and so well measured that the shapes of underlying orbits are almost directly visible in the data through element abundances or even stellar density. Thes… ▽ More The vertical kinematics of stars near the Sun can be used to measure the total mass distribution near the Galactic disk and to study out-of-equilibrium dynamics. With contemporary stellar surveys, the tracers of vertical dynamics are so numerous and so well measured that the shapes of underlying orbits are almost directly visible in the data through element abundances or even stellar density. These orbits can be used to infer a mass model for the Milky Way, enabling constraints on the dark matter distribution in the inner galaxy. Here we present a flexible model for foliating the vertical position-velocity phase space with orbits, for use in data-driven studies of dynamics. The vertical acceleration profile in the vicinity of the disk, along with the orbital actions, angles, and frequencies for individual stars, can all be derived from that orbit foliation. We show that this framework - "Orbital Torus Imaging" (OTI) - is rigorously justified in the context of dynamical theory, and does a good job of fitting orbits to simulated stellar abundance data with varying degrees of realism. OTI (1) does not require a global model for the Milky Way mass distribution, and (2) does not require detailed modeling of the selection function of the input survey data. We discuss the approximations and limitations of the OTI framework, which currently trades dynamical interpretability for flexibility in representing the data in some regimes, and which also presently separates the vertical and radial dynamics. We release an open-source tool, torusimaging, to accompany this article. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 32 pages, 10 figures. See arXiv:2312.07664 for an application of this method with data from APOGEE and Gaia

Spectacular nucleosynthesis from early massive stars

Authors: Alexander P. Ji, Sanjana Curtis, Nicholas Storm, Vedant Chandra, Kevin C. Schlaufman, Keivan G. Stassun, Alexander Heger, Marco Pignatari, Adrian M. Price-Whelan, Maria Bergemann, Guy S. Stringfellow, Carla Frohlich, Henrique Reggiani, Erika M. Holmbeck, Jamie Tayar, Shivani P. Shah, Emily J. Griffith, Chervin F. P. Laporte, Andrew R. Casey, Keith Hawkins, Danny Horta, William Cerny, Pierre Thibodeaux, Sam A. Usman, Joao A. S. Amarante , et al. (17 additional authors not shown)

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

Submitted 4 January, 2024; originally announced January 2024.

Comments: 11 pages + 22 page appendix, accepted to ApJL

Orbital Torus Imaging: Acceleration, density, and dark matter in the Galactic disk measured with element abundance gradients

Authors: Danny Horta, Adrian M. Price-Whelan, David W. Hogg, Kathryn V. Johnston, Lawrence Widrow, Julianne J. Dalcanton, Melissa K. Ness, Jason A. S. Hunt

Abstract: Under the assumption of a simple and time-invariant gravitational potential, many Galactic dynamics techniques infer the Milky Way's mass and dark matter distribution from stellar kinematic observations. These methods typically rely on parameterized potential models of the Galaxy and must take into account non-trivial survey selection effects, because they make use of the density of stars in phase… ▽ More Under the assumption of a simple and time-invariant gravitational potential, many Galactic dynamics techniques infer the Milky Way's mass and dark matter distribution from stellar kinematic observations. These methods typically rely on parameterized potential models of the Galaxy and must take into account non-trivial survey selection effects, because they make use of the density of stars in phase space. Large-scale spectroscopic surveys now supply information beyond kinematics in the form of precise stellar label measurements (especially element abundances). These element abundances are known to correlate with orbital actions or other dynamical invariants. Here, we use the Orbital Torus Imaging (OTI) framework that uses abundance gradients in phase space to map orbits. In many cases these gradients can be measured without detailed knowledge of the selection function. We use stellar surface abundances from the APOGEE survey combined with kinematic data from the Gaia mission. Our method reveals the vertical ($z$-direction) orbit structure in the Galaxy and enables empirical measurements of the vertical acceleration field and orbital frequencies in the disk. From these measurements, we infer the total surface mass density, $Σ$, and midplane volume density, $ρ_0$, as a function of Galactocentric radius and height. Around the Sun, we find $Σ_{\odot}(z=1.1$ kpc)$=72^{+6}_{-9}$M$_{\odot}$pc$^{-2}$ and $ρ_{\odot}(z=0)=0.081^{+0.015}_{-0.009}$ M$_{\odot}$pc$^{-3}$ using the most constraining abundance ratio, [Mg/Fe]. This corresponds to a dark matter contribution in surface density of $Σ_{\odot,\mathrm{DM}}(z=1.1$ kpc)$=24\pm4$ M$_{\odot}$pc$^{-2}$, and in total volume mass density of $ρ_{\odot,\mathrm{DM}}(z=0)=0.011\pm0.002$ M$_{\odot}$pc$^{-3}$. Moreover, using these mass density values we estimate the scale length of the low-$α$ disc to be $h_R=2.24\pm0.06$kpc. △ Less

Submitted 18 December, 2023; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: Accepted for publication in ApJ. 19 pages, 11 figures, 3 Tables

Stream Members Only: Data-Driven Characterization of Stellar Streams with Mixture Density Networks

Authors: Nathaniel Starkman, Jacob Nibauer, Jo Bovy, Jeremy J. Webb, Kiyan Tavangar, Adrian Price-Whelan, Ana Bonaca

Abstract: Stellar streams are sensitive probes of the Milky Way's gravitational potential. The mean track of a stream constrains global properties of the potential, while its fine-grained surface density constrains galactic substructure. A precise characterization of streams from potentially noisy data marks a crucial step in inferring galactic structure, including the dark matter, across orders of magnitud… ▽ More Stellar streams are sensitive probes of the Milky Way's gravitational potential. The mean track of a stream constrains global properties of the potential, while its fine-grained surface density constrains galactic substructure. A precise characterization of streams from potentially noisy data marks a crucial step in inferring galactic structure, including the dark matter, across orders of magnitude in mass scales. Here we present a new method for constructing a smooth probability density model of stellar streams using all of the available astrometric and photometric data. To characterize a stream's morphology and kinematics, we utilize mixture density networks to represent its on-sky track, width, stellar number density, and kinematic distribution. We model the photometry for each stream as a single-stellar population, with a distance track that is simultaneously estimated from the stream's inferred distance modulus (using photometry) and parallax distribution (using astrometry). We use normalizing flows to characterize the distribution of background stars. We apply the method to the stream GD-1, and the tidal tails of Palomar 5. For both streams we obtain a catalog of stellar membership probabilities that are made publicly available. Importantly, our model is capable of handling data with incomplete phase-space observations, making our method applicable to the growing census of Milky Way stellar streams. When applied to a population of streams, the resulting membership probabilities from our model form the required input to infer the Milky Way's dark matter distribution from the scale of the stellar halo down to subhalos. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 35 pages, 13 figures, 4 tables, fully open-source and reproducible using ShowYourWork

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 Future of Astronomical Data Infrastructure: Meeting Report

Authors: Michael R. Blanton, Janet D. Evans, Dara Norman, William O'Mullane, Adrian Price-Whelan, Luca Rizzi, Alberto Accomazzi, Megan Ansdell, Stephen Bailey, Paul Barrett, Steven Berukoff, Adam Bolton, Julian Borrill, Kelle Cruz, Julianne Dalcanton, Vandana Desai, Gregory P. Dubois-Felsmann, Frossie Economou, Henry Ferguson, Bryan Field, Dan Foreman-Mackey, Jaime Forero-Romero, Niall Gaffney, Kim Gillies, Matthew J. Graham , et al. (47 additional authors not shown)

Abstract: The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and productio… ▽ More The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and production of both astronomical data and software. This report is the culmination of a workshop held in February 2023 on the Future of Astronomical Data Infrastructure. Attended by 70 scientists and software professionals from ground-based and space-based missions and archives spanning the entire spectrum of astronomical research, the group deliberated on the prevailing state of software and data infrastructure in astronomy, identified pressing issues, and explored potential solutions. In this report, we describe the ecosystem of astronomical data, its existing flaws, and the many gaps, duplication, inconsistencies, barriers to access, drags on productivity, missed opportunities, and risks to the long-term integrity of essential data sets. We also highlight the successes and failures in a set of deep dives into several different illustrative components of the ecosystem, included as an appendix. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 59 pages; please send comments and/or questions to foadi@googlegroups.com

Chemical Cartography of the Sagittarius Stream with Gaia

Authors: Emily C. Cunningham, Jason A. S. Hunt, Adrian M. Price-Whelan, Kathryn V. Johnston, Melissa K. Ness, Yuxi Lu, Ivanna Escala, Ioana A. Stelea

Abstract: The stellar stream connected to the Sagittarius (Sgr) dwarf galaxy is the most massive tidal stream that has been mapped in the Galaxy, and is the dominant contributor to the outer stellar halo of the Milky Way. We present metallicity maps of the Sgr stream, using 34,240 red giant branch stars with inferred metallicities from Gaia BP/RP spectra. This sample is larger than previous samples of Sgr s… ▽ More The stellar stream connected to the Sagittarius (Sgr) dwarf galaxy is the most massive tidal stream that has been mapped in the Galaxy, and is the dominant contributor to the outer stellar halo of the Milky Way. We present metallicity maps of the Sgr stream, using 34,240 red giant branch stars with inferred metallicities from Gaia BP/RP spectra. This sample is larger than previous samples of Sgr stream members with chemical abundances by an order of magnitude. We measure metallicity gradients with respect to Sgr stream coordinates $(Λ, B)$, and highlight the gradient in metallicity with respect to stream latitude coordinate $B$, which has not been observed before. We find $\nabla \mathrm{[M/H]} = -2.48 \pm 0.08 \times 10^{-2}$ dex/deg above the stream track ($B>B_0$ where $B_0=1.5$ deg is the latitude of the Sgr remnant) and $\nabla \mathrm{[M/H]} =- 2.02 \pm 0.08 \times 10^{-2}$ dex/deg below the stream track ($B<B_0$). By painting metallicity gradients onto a tailored N-body simulation of the Sgr stream, we find that the observed metallicities in the stream are consistent with an initial radial metallicity gradient in the Sgr dwarf galaxy of $\sim -0.1$ to $-0.2$ dex/kpc, well within the range of observed metallicity gradients in Local Group dwarf galaxies. Our results provide novel observational constraints for the internal structure of the dwarf galaxy progenitor of the Sgr stream. Leveraging new large datasets in conjunction with tailored simulations, we can connect the present day properties of disrupted dwarfs in the Milky Way to their initial conditions. △ Less

Submitted 17 July, 2023; originally announced July 2023.

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

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

Roman CCS White Paper: Adding Fields Hosting Globular Clusters To The Galactic Bulge Time Domain Survey

Authors: Samuel K. Grunblatt, Robert F. Wilson, Andrew Winter, B. Scott Gaudi, Daniel Huber, Daniel A. Yahalomi, Andrea Bellini, Zachary R. Claytor, Jorge Martinez Palomera, Thomas Barclay, Guangwei Fu, Adrian Price-Whelan

Abstract: Despite multiple previous searches, no transiting planets have yet been identified within a globular cluster. This is believed to be due to a combination of factors: the low metallicities of most globular clusters suggests that there is significantly less planet-forming material per star in most globular clusters relative to the solar neighborhood, the high likelihood of dynamical interactions can… ▽ More Despite multiple previous searches, no transiting planets have yet been identified within a globular cluster. This is believed to be due to a combination of factors: the low metallicities of most globular clusters suggests that there is significantly less planet-forming material per star in most globular clusters relative to the solar neighborhood, the high likelihood of dynamical interactions can also disrupt planetary orbits, and the data available for globular clusters is limited. However, transiting planets have been identified in open clusters, indicating that there may be planets in more massive clusters that have simply gone undetected, or that more massive clusters inhibit planet formation. Less than two degrees away from the nominal Galactic Bulge Time Domain Survey footprint, two globular clusters, NGC 6522 and NGC 6528, can be simultaneously observed by the Roman telescope during the Galactic Bulge Time Domain Survey. These clusters are comparable in mass (1-2 x 10$^5$ solar masses) and age (12 Gyr), but feature drastically different average metallicities: NGC 6522 has an average [Fe/H] $\sim$ -1.3, while NGC 6528 has an average [Fe/H] $\sim$ -0.1. If no transiting planets are detected in one season of time domain observations of these clusters, this would indicate a difference in planet occurrence among field stars and globular clusters at >3-$σ$ significance even after accounting for metallicity, which could be enhanced to >5-$σ$ significance with similar observations of another nearby field hosting a metal-rich globular cluster. Additionally, time domain observations of NGC 6522 and NGC 6528 will detect variable stars in both clusters, testing the connection between stellar variability and binary fraction to metallicity and cluster environment, as well as testing the dependence of exoplanet yields on stellar density and distance from the Galactic midplane. △ Less

Submitted 18 June, 2023; originally announced June 2023.

Comments: 7 pages, 3 figures

Prospects for Detecting Gaps in Globular Cluster Stellar Streams in External Galaxies with the Nancy Grace Roman Space Telescope

Authors: Christian Aganze, Sarah Pearson, Tjitske Starkenburg, Gabriella Contardo, Kathryn V. Johnston, Kiyan Tavangar, Adrian M. Price-Whelan, Adam J. Burgasser

Abstract: Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and, therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter m… ▽ More Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and, therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter models. However, current samples are limited to the Milky Way. With its large field of view, deep imaging sensitivity, and high angular resolution, the upcoming Nancy Grace Roman Space Telescope (Roman) presents a unique opportunity to increase the number of observed streams and gaps significantly. This paper presents a first exploration of the prospects for detecting gaps in streams in M31 and other nearby galaxies with resolved stars. We simulate the formation of gaps in a Palomar-5-like stream and generate mock observations of these gaps with background stars in M31 and the foreground Milky Way stellar fields. We assess Roman's ability to detect gaps out to 10 Mpc through visual inspection and with the gap-finding tool ${\texttt{FindTheGap}}$. We conclude that gaps of $\approx 1.5$ kpc in streams that are created from subhalos of masses $\geq5 \times 10^6$M$_{\odot}$ are detectable within a 2-3 Mpc volume in exposures of 1000s to 1 hour. This volume contains $\approx 150$ galaxies, including $\approx 8$ galaxies with luminosities $>10^{9}~$L$_{\odot}$. Large samples of stream gaps in external galaxies will open up a new era of statistical analyses of gap characteristics in stellar streams and help constrain dark matter models. △ Less

Submitted 15 February, 2024; v1 submitted 19 May, 2023; originally announced May 2023.

Comments: ApJ version

Estimating the selection function of Gaia DR3 sub-samples

Authors: A. Castro-Ginard, A. G. A. Brown, Z. Kostrzewa-Rutkowska, T. Cantat-Gaudin, R. Drimmel, S. Oh, V. Belokurov, A. R. Casey, M. Fouesneau, S. Khanna, A. M. Price-Whelan, H. W. Rix

Abstract: Understanding which sources are present in an astronomical catalogue and which are not is crucial for the accurate interpretation of astronomical data. In particular, for the multidimensional Gaia data, filters and cuts on different parameters or measurements introduces a selection function that may unintentionally alter scientific conclusions in subtle ways. We aim to develop a methodology to est… ▽ More Understanding which sources are present in an astronomical catalogue and which are not is crucial for the accurate interpretation of astronomical data. In particular, for the multidimensional Gaia data, filters and cuts on different parameters or measurements introduces a selection function that may unintentionally alter scientific conclusions in subtle ways. We aim to develop a methodology to estimate the selection function for different sub-samples of stars in the Gaia catalogue. Comparing the number of stars in a given sub-sample to those in the overall Gaia catalogue, provides an estimate of the sub-sample membership probability, as a function of sky position, magnitude and colour. This estimate must differentiate the stochastic absence of sub-sample stars from selection effects. When multiplied with the overall Gaia catalogue selection function this provides the total selection function of the sub-sample. We present the method by estimating the selection function of the sources in Gaia DR3 with heliocentric radial velocity measurements. We also compute the selection function for the stars in the Gaia-Sausage/Enceladus sample, confirming that the apparent asymmetry of its debris across the sky is merely caused by selection effects. The developed method estimates the selection function of the stars present in a sub-sample of Gaia data, given that the sub-sample is completely contained in the Gaia parent catalogue (for which the selection function is known). This tool is made available in a GaiaUnlimited Python package. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 13 pages, 12 figures. Submitted to A&A

Journal ref: A&A 677, A37 (2023)

$\texttt{ESCARGOT}$: Mapping Vertical Phase Spiral Characteristics Throughout the Real and Simulated Milky Way

Authors: Elise Darragh-Ford, Jason A. S. Hunt, Adrian M. Price-Whelan, Kathryn V. Johnston

Abstract: The recent discovery of a spiral pattern in the vertical kinematic structure in the solar neighborhood provides a prime opportunity to study non-equilibrium dynamics in the Milky Way from local stellar kinematics. Furthermore, results from simulations indicate that even in a limited volume, differences in stellar orbital histories allow us to trace variations in the initial perturbation across lar… ▽ More The recent discovery of a spiral pattern in the vertical kinematic structure in the solar neighborhood provides a prime opportunity to study non-equilibrium dynamics in the Milky Way from local stellar kinematics. Furthermore, results from simulations indicate that even in a limited volume, differences in stellar orbital histories allow us to trace variations in the initial perturbation across large regions of the disk. We present $\texttt{ESCARGOT}$, a novel algorithm for studying these variations in both simulated and observed data sets. $\texttt{ESCARGOT}$ automatically extracts key quantities from the structure of a given phase spiral, including the time since perturbation and the perturbation mode. We test $\texttt{ESCARGOT}$ on simulated data and show that it is capable of accurately recovering information about the time since the perturbation occurred as well as subtle differences in phase spiral morphology due to stellar location in the disk at the time of perturbation. We apply $\texttt{ESCARGOT}$ to kinematic data from data release 3 of the ${\it Gaia}$ mission in bins of guiding radius. We show that similar structural differences in morphology occur in the ${\it Gaia}$ phase spirals as a function of stellar orbital history. These results indicate that the phase spirals are the product of a complex dynamical response in the disk with large-scale coupling between different regions of phase space. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Comments: 20 pages, 12 figures

Stream Fanning and Bifurcations: Observable Signatures of Resonances in Stellar Stream Morphology

Authors: Tomer D. Yavetz, Kathryn V. Johnston, Sarah Pearson, Adrian M. Price-Whelan, Chris Hamilton

Abstract: Recent observations have revealed a trove of unexpected morphological features in many of the Milky Way's stellar streams. Explanations for such features include time-dependent deformations of the Galactic gravitational potential, local disruptions induced by dark matter substructure, and special configurations of the streams' progenitors. In this paper, we study how these morphologies can also ar… ▽ More Recent observations have revealed a trove of unexpected morphological features in many of the Milky Way's stellar streams. Explanations for such features include time-dependent deformations of the Galactic gravitational potential, local disruptions induced by dark matter substructure, and special configurations of the streams' progenitors. In this paper, we study how these morphologies can also arise in certain static, non-spherical gravitational potentials that host a subset of resonantly-trapped orbit families. The transitions, or separatrices, between these orbit families mark abrupt discontinuities in the orbital structure of the potential. We develop a novel numerical approach for measuring the libration frequencies of resonant and near-resonant orbits, and apply it to study the evolution of stellar streams on these orbits. We reveal two distinct morphological features that arise in streams on near-resonant orbits: fans, that come about due to a large spread in the libration frequencies near a separatrix; and bifurcations, that arise when a separatrix splits the orbital distribution of the stellar stream between two (or more) distinct orbit families. We demonstrate that these effects can arise in some Milky Way streams for certain choices of the dark matter halo potential, and discuss how this might be used to probe and constrain the global shape of the Milky Way's gravitational potential. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 24 pages, 8 figures. Submitted to ApJ; comments welcome!

The Poor Old Heart of the Milky Way

Authors: Hans-Walter Rix, Vedant Chandra, René Andrae, Adrian M. Price-Whelan, David H. Weinberg, Charlie Conroy, Morgan Fouesneau, David W. Hogg, Francesca De Angeli, Rohan P. Naidu, Maosheng Xiang, Daniela Ruz-Mieres

Abstract: Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet pr… ▽ More Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet provide a comprehensive picture of such a metal-poor "heart" of the Milky Way. We use information from Gaia DR3, especially the XP spectra, to construct a sample of 2 million bright (BP $<15.5$ mag) giant stars within $30^\circ$ of the Galactic Center with robust [M/H] estimates, $δ$ [M/H] $\lesssim 0.1$. For most sample members we can calculate orbits based on Gaia RVS velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes $\sim 18,000$ stars with $-2.7<$ [M/H] $<-1.5$, representing a stellar mass of $\gtrsim 5\times 10^7$ M$_\odot$. The spatial distribution of these [M/H] $<-1.5$ stars has a Gaussian extent of only $σ_{\mathrm{R_{GC}}} \sim 2.7$ kpc around the Galactic center, with most of these orbits being confined to the inner Galaxy. At high orbital eccentricities, there is clear evidence for accreted halo stars in their pericentral orbit phase. Stars with [M/H] $< -2$ show no net rotation, whereas those with [M/H] $\sim -1$ are rotation dominated. Most of the tightly bound stars show $[α/\text{Fe}]$-enhancement and [Al/Fe]-[Mn/Fe] abundance patterns expected for an origin in the more massive portions of the proto-Galaxy. These central, metal-poor stars most likely predate the oldest part of the disk ($τ_{\text{age}}\approx 12.5$ Gyrs), which implies that they formed at $z\gtrsim 5$, forging the proto-Milky Way. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 21 pages, 9 figures, submitted to ApJ

From Data to Software to Science with the Rubin Observatory LSST

Authors: Katelyn Breivik, Andrew J. Connolly, K. E. Saavik Ford, Mario Jurić, Rachel Mandelbaum, Adam A. Miller, Dara Norman, Knut Olsen, William O'Mullane, Adrian Price-Whelan, Timothy Sacco, J. L. Sokoloski, Ashley Villar, Viviana Acquaviva, Tomas Ahumada, Yusra AlSayyad, Catarina S. Alves, Igor Andreoni, Timo Anguita, Henry J. Best, Federica B. Bianco, Rosaria Bonito, Andrew Bradshaw, Colin J. Burke, Andresa Rodrigues de Campos , et al. (75 additional authors not shown)

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the po… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

White dwarf binaries Across the H-R Diagram

Authors: Borja Anguiano, Steven R. Majewski, Keivan G. Stassun, Carles Badenes, Christine Mazzola Daher, Don Dixon, Carlos Allende Prieto, Donald P. Schneider, Adrian M. Price-Whelan, Rachael L. Beaton

Abstract: We created the APOGEE-GALEX-\emph{Gaia} catalog to study white dwarfs binaries. This database aims to create a minimally biased sample of WD binary systems identified from a combination of GALEX, {\it Gaia}, and APOGEE data to increase the number of WD binaries with orbital parameters and chemical compositions. We identify 3,414 sources as WD binary candidates, with nondegenerate companions of spe… ▽ More We created the APOGEE-GALEX-\emph{Gaia} catalog to study white dwarfs binaries. This database aims to create a minimally biased sample of WD binary systems identified from a combination of GALEX, {\it Gaia}, and APOGEE data to increase the number of WD binaries with orbital parameters and chemical compositions. We identify 3,414 sources as WD binary candidates, with nondegenerate companions of spectral types between F and M, including main sequence, main sequence binaries, subgiants, sub-subgiants, red giants, and red clump stars. Among our findings are (a) a total of 1,806 systems having inferred WD radii $R < 25$ R$_{\Earth}$, which constitute a more reliable group of WD binary candidates within the main sample; (b) a difference in the metallicity distribution function between WD binary candidates and the control sample of most luminous giants ($M_H < -3.0$); (c) the existence of a population of sub-subgiants with WD companions; (d) evidence for shorter periods in binaries that contain WDs compared to those that do not, as shown by the cumulative distributions of APOGEE radial velocity shifts; (e) evidence for systemic orbital evolution in a sample of 252 WD binaries with orbital periods, based on differences in the period distribution between systems with red clump, main sequence binary, and sub-subgiant companions and systems with main sequence or red giant companions; and (f) evidence for chemical enrichment during common envelope (CE) evolution, shown by lower metallicities in wide WD binary candidates ($P > 100$ days) compared to post-CE ($P < 100$ days) WD binary candidates. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: 22 pages, 18 figures. Accepted for publication in AJ

The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package

Authors: The Astropy Collaboration, Adrian M. Price-Whelan, Pey Lian Lim, Nicholas Earl, Nathaniel Starkman, Larry Bradley, David L. Shupe, Aarya A. Patil, Lia Corrales, C. E. Brasseur, Maximilian Nöthe, Axel Donath, Erik Tollerud, Brett M. Morris, Adam Ginsburg, Eero Vaher, Benjamin A. Weaver, James Tocknell, William Jamieson, Marten H. van Kerkwijk, Thomas P. Robitaille, Bruce Merry, Matteo Bachetti, H. Moritz Günther, Thomas L. Aldcroft , et al. (111 additional authors not shown)

Abstract: The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as… ▽ More The Astropy Project supports and fosters the development of open-source and openly-developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 43 pages, 5 figures. To appear in ApJ. The author list has two parts: the authors that made significant contributions to the writing and/or coordination of the paper, followed by maintainers of and contributors to the Astropy Project. The position in the author list does not correspond to contributions to the Astropy Project as a whole

The Open Cluster Chemical Abundances and Mapping Survey: VI. Galactic Chemical Gradient Analysis from APOGEE DR17

Authors: Natalie Myers, John Donor, Taylor Spoo, Peter M. Frinchaboy, Katia Cunha, Adrian M. Price-Whelan, Steven R. Majewski, Rachael L. Beaton, Gail Zasowski, Julia O'Connell, Amy E. Ray, Dmitry Bizyaev, Cristina Chiappini, D. A. García-Hernández, Doug Geisler, Henrik Jönsson, Richard R. Lane, Penélope Longa-Peña, Ivan Minchev, Dante Minniti, Christian Nitschelm, A. Roman-Lopes

Abstract: The goal of the Open Cluster Chemical Abundances and Mapping (OCCAM) survey is to constrain key Galactic dynamic and chemical evolution parameters by the construction and analysis of a large, comprehensive, uniform data set of infrared spectra for stars in hundreds of open clusters. This sixth contribution from the OCCAM survey presents analysis of SDSS/APOGEE Data Release 17 (DR17) results for a… ▽ More The goal of the Open Cluster Chemical Abundances and Mapping (OCCAM) survey is to constrain key Galactic dynamic and chemical evolution parameters by the construction and analysis of a large, comprehensive, uniform data set of infrared spectra for stars in hundreds of open clusters. This sixth contribution from the OCCAM survey presents analysis of SDSS/APOGEE Data Release 17 (DR17) results for a sample of stars in 150 open clusters, 94 of which we designate to be "high quality'' based on the appearance of their color-magnitude diagram. We find the APOGEE DR17-derived [Fe/H] values to be in good agreement with those from previous high resolution spectroscopic open cluster abundance studies. Using a subset of the high quality sample, the Galactic abundance gradients were measured for 16 chemical elements, including [Fe/H], for both Galactocentric radius ($R_{GC}$) and guiding center radius ($R_{Guide}$). We find an overall Galactic [Fe/H] vs $R_{GC}$ gradient of $-0.073 \pm 0.002$ dex/kpc over the range of $6 < R_{GC} < 11.5$ kpc, and a similar gradient is found for [Fe/H] versus $R_{Guide}$. Significant Galactic abundance gradients are also noted for O, Mg, S, Ca, Mn, Na, Al, K and Ce. Our large sample additionally allows us to explore the evolution of the gradients in four age bins for the remaining 15 elements. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 21 pages, 14 Figures, Astronomical Journal Accepted

Constraining the Planet Occurrence Rate around Halo Stars of Potentially Extragalactic Origin

Authors: Stephanie Yoshida, Samuel Grunblatt, Adrian Price-Whelan

Abstract: The search for planets orbiting other stars has recently expanded to include stars from galaxies outside the Milky Way. With the TESS and Gaia surveys, photometric and kinematic information can be combined to identify transiting planet candidates of extragalactic origin. Here, 1,080 low-luminosity red giant branch stars observed by Gaia and TESS with kinematics suggesting a high likelihood of extr… ▽ More The search for planets orbiting other stars has recently expanded to include stars from galaxies outside the Milky Way. With the TESS and Gaia surveys, photometric and kinematic information can be combined to identify transiting planet candidates of extragalactic origin. Here, 1,080 low-luminosity red giant branch stars observed by Gaia and TESS with kinematics suggesting a high likelihood of extragalactic origin were searched for planet transits. Transit injection-recovery tests were performed to measure the sensitivity of the TESS data and completeness of the transit search. Injected signals of planets larger than Jupiter with orbital periods of 10 days or less were recovered in $\approx$44% of cases. Although no planet transits were detected in this sample, we find an upper limit on planet occurrence of 0.52% for hot Jupiters, consistent with previous studies of planet occurrence around similar host stars. As stars in the halo tend to be lower metallicity, and short period giant planet occurrence tends to be strongly correlated with stellar metallicity, we predict that relative to the Galactic disk population, a smaller fraction of halo stars will host planets detectable by transit surveys. Thus, applying the known planet occurrence trends to potential planet detection around halo stars, we predict $\gtrsim$7,000 stars must be searched with similar cadence and precision as the stars studied here before a detection of a planet of extragalactic origin is likely. This may be possible with future data releases from the TESS and Gaia missions. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 9 pages, 4 figures, submitted to AAS Journals, additional comments welcome

paired: A Statistical Framework for Determining Stellar Binarity with Gaia RVs. I. Planet Hosting Binaries

Authors: Quadry Chance, Daniel Foreman-Mackey, Sarah Ballard, Andrew Casey, Trevor David, Adrian Price-Whelan

Abstract: The effect of stellar multiplicity on the formation and evolution of planetary systems is complex. At a demographic level, campaigns with both high-resolution imaging and radial velocity observations indicate that planet formation is strongly disrupted by close binaries, while being relatively unaffected by wide companions. However, the magnitude and distance-limited nature of those tools mean tha… ▽ More The effect of stellar multiplicity on the formation and evolution of planetary systems is complex. At a demographic level, campaigns with both high-resolution imaging and radial velocity observations indicate that planet formation is strongly disrupted by close binaries, while being relatively unaffected by wide companions. However, the magnitude and distance-limited nature of those tools mean that large ranges of mass ratios and separations remain largely unexplored. The Early Data Release 3 (EDR3) from the Gaia Mission includes radial velocity measurements of over 6.5 million targets, which we employ to explore the effect of binary companions within a statistical framework called paired. These companions present as a source of excess radial velocity noise in the Gaia catalog, when compared to the typical noise for stars of similar spectral type and magnitude. Within this framework, we examine the evidence for stellar multiplicity among the stars surveyed by NASA's Kepler and TESS missions. We use radial velocity errors published in Gaia EDR3 to estimate the probability of an unresolved stellar companion for a large subset of the Kepler and TESS Input Catalog stars, where possible benchmarking our inferred radial velocity semi-amplitudes against the those from ground-based radial velocity surveys. We determine that we are typically sensitive out to several AU and mass ratios $>0.1$, dependent upon the stellar magnitude. We aim for paired to be a useful community tool for the exploration of the effects of binarity on planets at a population level, and for efficient identification of false-positive transit candidates. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 29 pages, 7 figures, Submitted to ApJ. The catalog can be found at https://doi.org/10.5281/zenodo.6677444

Carbon-Enhanced Metal-Poor star candidates from BP/RP Spectra in $Gaia$ DR3

Authors: Madeline Lucey, Nariman Al Kharusi, Keith Hawkins, Yuan-Sen Ting, Nesar Ramachandra, Adrian M. Price-Whelan, Timothy C. Beers, Young Sun Lee, Jinmi Yoon

Abstract: Carbon-enhanced metal-poor (CEMP) stars comprise almost a third of stars with [Fe/H] < --2, although their origins are still poorly understood. It is highly likely that one sub-class (CEMP-$s$ stars) is tied to mass-transfer events in binary stars, while another sub-class (CEMP-no stars) are enriched by the nucleosynthetic yields of the first generations of stars. Previous studies of CEMP stars ha… ▽ More Carbon-enhanced metal-poor (CEMP) stars comprise almost a third of stars with [Fe/H] < --2, although their origins are still poorly understood. It is highly likely that one sub-class (CEMP-$s$ stars) is tied to mass-transfer events in binary stars, while another sub-class (CEMP-no stars) are enriched by the nucleosynthetic yields of the first generations of stars. Previous studies of CEMP stars have primarily concentrated on the Galactic halo, but more recently they have also been detected in the thick disk and bulge components of the Milky Way. $Gaia$ DR3 has provided an unprecedented sample of over 200 million low-resolution ($R\approx$ 50) spectra from the BP and RP photometers. Training on the CEMP catalog from the SDSS/SEGUE database, we use XGBoost to identify the largest all-sky sample of CEMP candidate stars to date. In total, we find 58,872 CEMP star candidates, with an estimated contamination rate of 12%. When comparing to literature high-resolution catalogs, we positively identify 60-68% of the CEMP stars in the data, validating our results and indicating a high completeness rate. Our final catalog of CEMP candidates spans from the inner to outer Milky Way, with distances as close as $r \sim$ 0.8 kpc from the Galactic center, and as far as $r >$ 30 kpc. Future higher-resolution spectroscopic follow-up of these candidates will provide validations of their classification and enable investigations of the frequency of CEMP-$s$ and CEMP-no stars throughout the Galaxy, to further constrain the nature of their progenitors. △ Less

Submitted 1 June, 2023; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 19 pages, 14 figures, accepted to MNRAS

Multiple phase spirals suggest multiple origins in Gaia DR3

Authors: Jason A. S. Hunt, Adrian M. Price-Whelan, Kathryn V. Johnston, Elise Darragh-Ford

Abstract: {\it Gaia} Data Release 2 (DR2) revealed that the Milky Way contains significant indications of departures from equilibrium in the form of asymmetric features in the phase space density of stars in the Solar neighborhood. One such feature is the $z$--$v_z$ phase spiral, interpreted as the response of the disk to the influence of a perturbation perpendicular to the disk plane, which could be extern… ▽ More {\it Gaia} Data Release 2 (DR2) revealed that the Milky Way contains significant indications of departures from equilibrium in the form of asymmetric features in the phase space density of stars in the Solar neighborhood. One such feature is the $z$--$v_z$ phase spiral, interpreted as the response of the disk to the influence of a perturbation perpendicular to the disk plane, which could be external (e.g., a satellite) or internal (e.g., the bar or spiral arms). In this work we use {\it Gaia} DR3 to dissect the phase spiral by dividing the local data set into groups with similar azimuthal actions, $J_φ$, and conjugate angles, $θ_φ$, which selects stars on similar orbits and at similar orbital phases, thus having experienced similar perturbations in the past. These divisions allow us to explore areas of the Galactic disk larger than the surveyed region. The separation improves the clarity of the $z$--$v_z$ phase spiral and exposes changes to its morphology across the different action-angle groups. In particular, we discover a transition to two armed `breathing spirals' in the inner Milky Way. We conclude that the local data contains signatures of not one, but multiple perturbations with the prospect to use their distinct properties to infer the properties of the interactions that caused them. △ Less

Submitted 27 July, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 5 pages, 3 figures, accepted in MNRAS letters

Dynamically constraining the length of the Milky Way bar

Authors: Madeline Lucey, Sarah Pearson, Jason A. S. Hunt, Keith Hawkins, Melissa Ness, Michael S. Petersen, Adrian M. Price-Whelan, Martin D. Weinberg

Abstract: We present a novel method for constraining the length of the Galactic bar using 6D phase space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named $R_{Freq}$, based on the maximal extent of trapped bar orbits. We find the $R_{Freq}$ measured from orbits is consistent with the $R_{Freq}$ of the assumed potential only when the length of the bar and pattern… ▽ More We present a novel method for constraining the length of the Galactic bar using 6D phase space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named $R_{Freq}$, based on the maximal extent of trapped bar orbits. We find the $R_{Freq}$ measured from orbits is consistent with the $R_{Freq}$ of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model's or the Milky Way's bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured $R_{Freq}$. When we apply this method to $\approx$210,000 stars in APOGEE DR17 and $Gaia$ eDR3 data, we find a consistent result only for potential models with a dynamical bar length of $\approx$3.5 kpc. We find the Milky Way's trapped bar orbits extend out to only $\approx$3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential. △ Less

Submitted 2 February, 2023; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 15 pages, 8 figures, 2 tables, accepted to MNRAS, comments welcome

On the Hunt for the Origins of the Orphan--Chenab Stream: Detailed Element Abundances with APOGEE and Gaia

Authors: Keith Hawkins, Adrian M. Price-Whelan, Allyson A. Sheffield, Aidan Z. Subrahimovic, Rachael L. Beaton, Vasily Belokurov, Denis Erkal, Sergey E. Koposov, Richard R. Lane, Chervin F. P. Laporte, Christian Nitschelm

Abstract: Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the… ▽ More Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA's Gaia mission (Gaia EDR3) and radial velocity information from the SDSS-IV APOGEE survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 13 stars) of the OC stream in the $α$ (O, Mg, Ca, Si, Ti, S), odd-Z (Al, K, V), Fe-peak (Fe, Ni, Mn, Co, Cr) and neutron capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a mono-metallic population and have a median metallicity of --1.92~dex with a dispersion of 0.28 dex. Our results also indicate that the $α$-elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second generation stars from globular clusters. The detailed chemical pattern of these stars indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ~10$^6$ M$_\odot$. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 13 Pages; 4 Figures. Submitted to AAS Journals, comments welcome

The impact of a massive Sagittarius dSph on GD-1-like streams

Authors: Adam M. Dillamore, Vasily Belokurov, N. Wyn Evans, Adrian M. Price-Whelan

Abstract: We investigate the effects of a massive ($\gtrsim4\times10^{10}M_\odot$) Sagittarius dwarf spheroidal galaxy (Sgr) on stellar streams using test particle simulations in a realistic Milky Way potential. We find that Sgr can easily disrupt streams formed more than $\sim3$ Gyr ago, while stars stripped more recently are generally unaffected. In certain realizations, Sgr is able to produce asymmetry b… ▽ More We investigate the effects of a massive ($\gtrsim4\times10^{10}M_\odot$) Sagittarius dwarf spheroidal galaxy (Sgr) on stellar streams using test particle simulations in a realistic Milky Way potential. We find that Sgr can easily disrupt streams formed more than $\sim3$ Gyr ago, while stars stripped more recently are generally unaffected. In certain realizations, Sgr is able to produce asymmetry between the leading and trailing tails of Pal 5, qualitatively similar to observations. Using data from the Gaia space telescope and elsewhere, we fit models to the GD-1 stream in the presence of a Sgr with various initial masses. While the best-fitting models do show perturbations resulting from interactions with Sgr, we find that the level of disruption is not significantly greater than in the observed stream. To investigate the general effects of Sgr on a population of streams, we generate 1000 mock streams on GD-1-like orbits with randomized orientations. Some streams show clear evidence of disruption, becoming folded on the sky or developing asymmetry betweeen their two tails. However, many survive unaffected and the peak surface brightness of stars is decreased by no more than $\sim0.3$ mag/arcsec$^2$ on average. We conclude that Sgr having an initial mass of $\gtrsim4\times10^{10}M_\odot$ is compatible with the survival and detection of streams formed more than 3 Gyr ago. △ Less

Submitted 1 September, 2022; v1 submitted 26 May, 2022; originally announced May 2022.

Comments: 19 pages, 12 figures, published in MNRAS. Animated figures can be found at https://www.youtube.com/playlist?list=PLEleLLhXAwEMx6GcSror-iF-QsskHrXYM

Mapping Dark Matter with Extragalactic Stellar Streams: the Case of Centaurus A

Authors: Sarah Pearson, Adrian M. Price-Whelan, David W. Hogg, Anil C. Seth, David J. Sand, Jason A. S. Hunt, Denija Crnojevic

Abstract: In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A's (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned… ▽ More In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A's (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned about the Cen A dark-matter halo. We develop a novel external galaxy stream-fitting technique and generate model stellar streams that reproduce the stream morphology visible in the imaging. We find that there are many viable stream models that fit the data well, with reasonable parameters, provided that Cen A has a halo mass larger than M$_{200}$ $>4.70\times 10^{12}$ M$_{\odot}$. There is a second stream in Cen A's halo that is also reproduced within the context of this same dynamical model. However, stream morphology in the imaging alone does not uniquely determine the mass or mass distribution for the Cen A halo. In particular, the stream models with high likelihood show covariances between the inferred Cen A mass distribution, the inferred Dw3 progenitor mass, the Dw3 velocity, and the Dw3 line-of-sight position. We show that these degeneracies can be broken with radial-velocity measurements along the stream, and that a single radial velocity measurement puts a substantial lower limit on the halo mass. These results suggest that targeted radial-velocity measurements will be critical if we want to learn about dark matter from extragalactic stellar streams. △ Less

Submitted 20 October, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 11 figures, 23 pages, Accepted to ApJ

The 3D Galactocentric velocities of Kepler stars: marginalizing over missing RVs

Authors: Ruth Angus, Adrian M. Price-Whelan, Joel C. Zinn, Megan Bedell, Yuxi, Lu, Daniel Foreman-Mackey

Abstract: Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to dire… ▽ More Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities, however one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three-quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: Accepted for publication in AAS Journals

Implications of the Milky Way travel velocity for dynamical mass estimates of the Local Group

Authors: Katie Chamberlain, Adrian M. Price-Whelan, Gurtina Besla, Emily C. Cunningham, Nicolás Garavito-Camargo, Jorge Peñarrubia, Michael S. Petersen

Abstract: The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of… ▽ More The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history including the recent pericentric passage of the LMC, and recent work has found that the MW disk is moving with a lower bound "travel velocity" of $\sim 32~{\rm km}~{\rm s}^{-1}$ with respect to the outer stellar halo. Previous Timing Argument measurements attempt to account for this non-equilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10--12\% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume. △ Less

Submitted 30 December, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 14 pages, 3 figures; Published in ApJ; updated to match published version

Journal ref: ApJ 942 18 (2023)

Snowmass2021 Cosmic Frontier White Paper: Dark Matter Physics from Halo Measurements

Authors: Keith Bechtol, Simon Birrer, Francis-Yan Cyr-Racine, Katelin Schutz, Susmita Adhikari, Mustafa Amin, Arka Banerjee, Simeon Bird, Nikita Blinov, Kimberly K. Boddy, Celine Boehm, Kevin Bundy, Malte Buschmann, Sukanya Chakrabarti, David Curtin, Liang Dai, Alex Drlica-Wagner, Cora Dvorkin, Adrienne L. Erickcek, Daniel Gilman, Saniya Heeba, Stacy Kim, Vid Iršič, Alexie Leauthaud, Mark Lovell , et al. (19 additional authors not shown)

Abstract: The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determin… ▽ More The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics. △ Less

Submitted 24 April, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: White paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). 88 pages, 9 figures. Comments welcome

Astrophysical and Cosmological Probes of Dark Matter

Authors: Kimberly K. Boddy, Mariangela Lisanti, Samuel D. McDermott, Nicholas L. Rodd, Christoph Weniger, Yacine Ali-Haïmoud, Malte Buschmann, Ilias Cholis, Djuna Croon, Adrienne L. Erickcek, Vera Gluscevic, Rebecca K. Leane, Siddharth Mishra-Sharma, Julian B. Muñoz, Ethan O. Nadler, Priyamvada Natarajan, Adrian Price-Whelan, Simona Vegetti, Samuel J. Witte

Abstract: While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for t… ▽ More While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade. △ Less

Submitted 12 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021, 24 pages, 4 figures, comments welcome

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

Authors: Kareem El-Badry, Rhys Seeburger, Tharindu Jayasinghe, Hans-Walter Rix, Silvia Almada, Charlie Conroy, Adrian M. Price-Whelan, Kevin Burdge

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

Submitted 21 March, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 20 pages, 16 figures. Accepted to MNRAS

Stellar Abundance Maps of the Milky Way Disk

Authors: Anna-Christina Eilers, David W. Hogg, Hans-Walter Rix, Melissa K. Ness, Adrian M. Price-Whelan, Szabolcs Meszaros, Christian Nitschelm

Abstract: To understand the formation of the Milky Way's prominent bar it is important to know whether stars in the bar differ in the chemical element composition of their birth material as compared to disk stars. This requires stellar abundance measurements for large samples across the Milky Way's body. Such samples, e.g. luminous red giant stars observed by SDSS's Apogee survey, will inevitably span a ran… ▽ More To understand the formation of the Milky Way's prominent bar it is important to know whether stars in the bar differ in the chemical element composition of their birth material as compared to disk stars. This requires stellar abundance measurements for large samples across the Milky Way's body. Such samples, e.g. luminous red giant stars observed by SDSS's Apogee survey, will inevitably span a range of stellar parameters; as a consequence, both modelling imperfections and stellar evolution may preclude consistent and precise estimates of their chemical composition at a level of purported bar signatures, which has left current analyses of a chemically distinct bar inconclusive. Here, we develop a new self-calibration approach to eliminate both modelling and astrophysical abundance systematics among red giant branch (RGB) stars of different luminosities (and hence surface gravity $\log g$). We apply our method to $48,853$ luminous Apogee DR16 RGB stars to construct spatial abundance maps of $20$ chemical elements near the Milky Way's mid-plane, covering Galactocentric radii of $0\,{\rm kpc}<R_{\rm GC}<20\,\rm kpc$. Our results indicate that there are no abundance variations whose geometry matches that of the bar, and that the mean abundance gradients vary smoothly and monotonically with Galactocentric radius. We confirm that the high-$α$ disk is chemically homogeneous, without spatial gradients. Furthermore, we present the most precise [Fe/H] vs. $R_{\rm GC}$ gradient to date with a slope of $-0.057\pm0.001\rm~dex\,kpc^{-1}$ out to approximately $15$ kpc. △ Less

Submitted 18 February, 2022; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: ApJ accepted

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)

TOI-2076 and TOI-1807: Two young, comoving planetary systems within 50 pc identified by TESS that are ideal candidates for further follow-up

Authors: Christina Hedges, Alex Hughes, Steven Giacalone, George Zhou, Trevor J. David, Juliette Becker, Andrew Vanderburg, Joseph E. Rodriguez, Shaun Atherton, Samueln. Quinn, Courtney D. Dressing, Allyson Bieryla, Tara Fetherolf, Adrian Price-whelan, Megan Bedell, David W. Latham, Georger. Ricker, Roland K. Vanderspek, Sara Seager, Joshua N. Winn, Jon M. Jenkins, Rene Tronsgaard, Lars A. Buchhave, Karen A. Collins, Tianjun Gan , et al. (20 additional authors not shown)

Abstract: We report the discovery of two planetary systems around comoving stars; TOI-2076 (TIC 27491137) and TOI-1807 (TIC 180695581). TOI-2076 is a nearby (41.9 pc) multi-planetary system orbiting a young (204$\pm$50 Myr), bright (K = 7.115 in TIC v8.1). TOI-1807 hosts a single transiting planet, and is similarly nearby (42.58pc), similarly young (180$\pm$40 Myr), and bright. Both targets exhibit signific… ▽ More We report the discovery of two planetary systems around comoving stars; TOI-2076 (TIC 27491137) and TOI-1807 (TIC 180695581). TOI-2076 is a nearby (41.9 pc) multi-planetary system orbiting a young (204$\pm$50 Myr), bright (K = 7.115 in TIC v8.1). TOI-1807 hosts a single transiting planet, and is similarly nearby (42.58pc), similarly young (180$\pm$40 Myr), and bright. Both targets exhibit significant, periodic variability due to star spots, characteristic of their young ages. Using photometric data collected by TESS we identify three transiting planets around TOI-2076 with radii of R$_b$=3.3$\pm$0.04$R_\oplus$, R$_c$=4.4$\pm$0.05$R_\oplus$, and R$_d$=4.1$\pm$0.07$R_\oplus$. Planet TOI-2076b has a period of P$_b$=10.356 d. For both TOI 2076c and d, TESS observed only two transits, separated by a 2-year interval in which no data were collected, preventing a unique period determination. A range of long periods (>17d) are consistent with the data. We identify a short-period planet around TOI-1807 with a radius of R$_b$=1.8$\pm$0.04$R_\oplus$ and a period of P$_b$=0.549 d. Their close proximity, and bright, cool host stars, and young ages, make these planets excellent candidates for follow-up. TOI-1807b is one of the best known small ($R<2R_\oplus$) planets for characterization via eclipse spectroscopy and phase curves with JWST. TOI-1807b is the youngest ultra-short period planet discovered to date, providing valuable constraints on formation time-scales of short period planets. Given the rarity of young planets, particularly in multiple planet systems, these planets present an unprecedented opportunity to study and compare exoplanet formation, and young planet atmospheres, at a crucial transition age for formation theory. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 27 pages, 13 figures, 6 tables

Journal ref: Astrophysical Journal, 162, 54 (2021)

The homogeneity of the star forming environment of the Milky Way disk over time

Authors: Melissa K. Ness, Adam J. Wheeler, Kevin McKinnon, Danny Horta, Andrew R. Casey, Emily C. Cunningham, Adrian M. Price-Whelan

Abstract: Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$α$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using meas… ▽ More Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$α$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using measurements from $\sim$27,000 APOGEE stars (R=22,500, SNR$>$200), we build simple local linear models to predict a sample of elements (X = Si, O, Ca, Ti, Ni, Al, Mn, Cr) using (Fe, Mg) abundances alone, as fiducial tracers of supernovae production channels. Given [Fe/H] and [Mg/H], we predict these elements, [X/H], to about double the uncertainty of their measurements. The intrinsic dispersion, after subtracting measurement errors in quadrature is $\approx 0.015-0.04$~dex. The residuals of the prediction (measurement $-$ model) for each element demonstrate that each element has an individual link to birth properties at fixed (Fe, Mg). Residuals from primarily massive-star supernovae (i.e. Si, O, Al) partially correlate with guiding radius. Residuals from primarily supernovae Ia (i.e. Mn, Ni) partially correlate with age. A fraction of the intrinsic scatter that persists at fixed (Fe, Mg), however, after accounting for correlations, does not appear to further discriminate between birth properties that can be traced with present-day measurements. Presumably, this is because the residuals are also, in part, a measure of the typical (in)-homogeneity of the disk's stellar birth environments, previously inferred only using open-cluster systems. Our study implies at fixed birth radius and time, there is a median scatter of $\approx 0.01-0.015$ dex in elements generated in supernovae sources. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: submitted to ApJ, 16 pages, 12 figures

APOGEE Chemical Abundance Patterns of the Massive Milky Way Satellites

Authors: Sten Hasselquist, Christian R. Hayes, Jianhui Lian, David H. Weinberg, Gail Zasowski, Danny Horta, Rachael Beaton, Diane K. Feuillet, Elisa R. Garro, Carme Gallart, Verne V. Smith, Jon A. Holtzman, Dante Minniti, Ivan Lacerna, Matthew Shetrone, Henrik Jönsson, Maria-Rosa L. Cioni, Sean P. Fillingham, Katia Cunha, Robert OĆonnell, José G. Fernández-Trincado, Ricardo R. Muñoz, Ricardo Schiavon, Andres Almeida, Borja Anguiano , et al. (20 additional authors not shown)

Abstract: The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf (Sgr), Fornax (Fnx), and the now fully disrupted \emph{Gaia} Sausage/Enceladus (GSE) system. We present… ▽ More The SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey has obtained high-resolution spectra for thousands of red giant stars distributed among the massive satellite galaxies of the Milky Way (MW): the Large and Small Magellanic Clouds (LMC/SMC), the Sagittarius Dwarf (Sgr), Fornax (Fnx), and the now fully disrupted \emph{Gaia} Sausage/Enceladus (GSE) system. We present and analyze the APOGEE chemical abundance patterns of each galaxy to draw robust conclusions about their star formation histories, by quantifying the relative abundance trends of multiple elements (C, N, O, Mg, Al, Si, Ca, Fe, Ni, and Ce), as well as by fitting chemical evolution models to the [$α$/Fe]-[Fe/H] abundance plane for each galaxy. Results show that the chemical signatures of the starburst in the MCs observed by Nidever et al. in the $α$-element abundances extend to C+N, Al, and Ni, with the major burst in the SMC occurring some 3-4 Gyr before the burst in the LMC. We find that Sgr and Fnx also exhibit chemical abundance patterns suggestive of secondary star formation epochs, but these events were weaker and earlier ($\sim$~5-7 Gyr ago) than those observed in the MCs. There is no chemical evidence of a second starburst in GSE, but this galaxy shows the strongest initial star formation as compared to the other four galaxies. All dwarf galaxies had greater relative contributions of AGB stars to their enrichment than the MW. Comparing and contrasting these chemical patterns highlight the importance of galaxy environment on its chemical evolution. △ Less

Submitted 30 September, 2021; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 37 pages, 21 figures, accepted for publication in ApJ

Final Targeting Strategy for the SDSS-IV APOGEE-2S Survey

Authors: Felipe A. Santana, Rachael L. Beaton, Kevin R. Covey, Julia E. O'Connell, Penélope Longa-Peña, Roger Cohen, José G. Fernández-Trincado, Christian R. Hayes, Gail Zasowski, Jennifer S. Sobeck, Steven R. Majewski, S. D. Chojnowski, Nathan De Lee, Ryan J. Oelkers, Guy S. Stringfellow, Andrés Almeida, Borja Anguiano, John Donor, Peter M. Frinchaboy, Sten Hasselquist, Jennifer A. Johnson, Juna A. Kollmeier, David L. Nidever, Adrian. M. Price-Whelan, Alvaro Rojas-Arriagada , et al. (21 additional authors not shown)

Abstract: APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, A… ▽ More APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, APOGEE is uniquely able to probe the dust-hidden inner regions of the Milky Way that are best accessed from the Southern Hemisphere. In this paper we present the targeting strategy of APOGEE-2S, with special attention to documenting modifications to the original, previously published plan. The motivation for these changes is explained as well as an assessment of their effectiveness in achieving their intended scientific objective. In anticipation of this being the last paper detailing APOGEE targeting, we present an accounting of all such information complete through the end of the APOGEE-2S project; this includes several main survey programs dedicated to exploration of major stellar populations and regions of the Milky Way, as well as a full list of programs contributing to the APOGEE database through allocations of observing time by the Chilean National Time Allocation Committee (CNTAC) and the Carnegie Institution for Science (CIS). This work was presented along with a companion article, R. Beaton et al. (submitted; AAS29028), presenting the final target selection strategy adopted for APOGEE-2 in the Northern Hemisphere. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: 43 pages; 5 figures; 6 Tables; 1 Appendix; Submitted to Journal and Under Review; Posting to accompany papers using the SDSS-IV/APOGEE-2 Data Release 17 scheduled for December 2021

Final Targeting Strategy for the SDSS-IV APOGEE-2N Survey

Authors: Rachael L. Beaton, Ryan J. Oelkers, Christian R. Hayes, Kevin R. Covey, S. D. Chojnowski, Nathan De Lee, Jennifer S. Sobeck, Steven R. Majewski, Roger Cohen, Jose Fernandez-Trincado, Penelope Longa-Pena, Julia E. O'Connell, Felipe A. Santana, Guy S. Stringfellow, Gail Zasowski, Conny Aerts, Borja Anguiano, Chad Bender, Caleb I. Canas, Katia Cunha, John Donor Scott W. Fleming, Peter M. Frinchaboy, Diane Feuillet, Paul Harding, Sten Hasselquist , et al. (35 additional authors not shown)

Abstract: APOGEE-2 is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemo-dynamical mapping of the Milky Way Galaxy. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation… ▽ More APOGEE-2 is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemo-dynamical mapping of the Milky Way Galaxy. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs" competitively awarded to SDSS-IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5-year expansion of the survey, known as the Bright Time Extension, made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The Bright Time Extension permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new datasets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, Santana et al. (submitted), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: 59 pages; 11 Figures; 7 Tables; 2 Appendices; Submitted to Journal and Under Review; Posting to accompany papers using the SDSS-IV/APOGEE-2 Data Release 17 scheduled for December 2021

The Clustering of Orbital Poles Induced by the LMC: Hints for the Origin of Planes of Satellites

Authors: Nicolas Garavito-Camargo, Ekta Patel, Gurtina Besla, Adrian M. Price-Whelan, Facundo A. Gomez, Chervin F. P Laporte, Kathryn V. Johnston

Abstract: A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the… ▽ More A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter (DM) particles from the present-day random distribution. Instead, the orbital poles of particles beyond $R\approx 50$kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor $δρ_{max}$=30\%(50\%) with respect to underdense regions, and $δρ_{iso}$=15\%(30\%) relative to the isolated MW simulation (no LMC) between 50-150 kpc (150-300 kpc). The clustering appears after the LMC's pericenter ($\approx$ 50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: 1) the LMC shifts the velocity and position of the central density of the MW's halo and disk; 2) the DM dynamical friction wake and collective response induced by the LMC changes the kinematics of particles; 3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least $\approx$ 1:10 mass ratio), causing the clustering of orbital poles of halo tracers. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ. 23 pages and 16 figures. Comments are welcome

Double-lined spectroscopic binaries in the APOGEE DR16 and DR17 data

Authors: Marina Kounkel, Kevin R. Covey, Keivan G. Stassun, Adrian M. Price-Whelan, Jon Holtzman, Drew Chojnowski, Penélope Longa-Peña, Carlos G. Román-Zúñiga, Jesus Hernandez, Javier Serna, Carles Badenes, Nathan De Lee, Steven Majewski, Guy S. Stringfellow, Kaitlin M. Kratter, Maxwell Moe, Peter M. Frinchaboy, Rachael L. Beaton, José G. Fernández-Trincado, Suvrath Mahadevan, Dante Minniti, Timothy C. Beers, Donald P. Schneider, Rodolfo H. Barbá, Joel R. Brownstein , et al. (3 additional authors not shown)

Abstract: APOGEE spectra offer $\lesssim$1 km s$^{-1}$ precision in the measurement of stellar radial velocities (RVs). This holds even when multiple stars are captured in the same spectrum, as happens most commonly with double-lined spectroscopic binaries (SB2s), although random line of sight alignments of unrelated stars can also occur. We develop a code that autonomously identifies SB2s and higher order… ▽ More APOGEE spectra offer $\lesssim$1 km s$^{-1}$ precision in the measurement of stellar radial velocities (RVs). This holds even when multiple stars are captured in the same spectrum, as happens most commonly with double-lined spectroscopic binaries (SB2s), although random line of sight alignments of unrelated stars can also occur. We develop a code that autonomously identifies SB2s and higher order multiples in the APOGEE spectra, resulting in 7273 candidate SB2s, 813 SB3s, and 19 SB4s. We estimate the mass ratios of binaries, and for a subset of these systems with sufficient number of measurements we perform a complete orbital fit, confirming that most systems with period $<$10 days have circularized. Overall, we find a SB2 fraction ($F_{SB2}$) $\sim$3\% among main sequence dwarfs, and that there is not a significant trend in $F_{SB2}$ with temperature of a star. We are also able to recover a higher $F_{SB2}$ in sources with lower metallicity, however there are some observational biases. We also examine light curves from TESS to determine which of these spectroscopic binaries are also eclipsing. Such systems, particularly those that are also pre- and post-main sequence, are good candidates for a follow-up analysis to determine their masses and temperatures. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted to AJ. 23 pages, 17 figures

Snails Across Scales: Local and Global Phase-Mixing Structures as Probes of the Past and Future Milky Way

Authors: Suroor S. Gandhi, Kathryn V. Johnston, Jason A. S. Hunt, Adrian M. Price-Whelan, Chervin F. P. Laporte, David W. Hogg

Abstract: Signatures of vertical disequilibrium have been observed across the Milky Way's disk. These signatures manifest locally as unmixed phase-spirals in $z$--$v_z$ space ("snails-in-phase") and globally as nonzero mean $z$ and $v_z$ which wraps around as a physical spiral across the $x$--$y$ plane ("snails-in-space"). We explore the connection between these local and global spirals through the example… ▽ More Signatures of vertical disequilibrium have been observed across the Milky Way's disk. These signatures manifest locally as unmixed phase-spirals in $z$--$v_z$ space ("snails-in-phase") and globally as nonzero mean $z$ and $v_z$ which wraps around as a physical spiral across the $x$--$y$ plane ("snails-in-space"). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle Milky Way (MW)-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using a $z$--$v_z$ asymmetry metric we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal action $J_φ$, excited by a single event (in our case, a satellite disk-crossing); (b) these distinct phase-spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase-spirals offer a global view of the MW's perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)-like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars within $R \lesssim 10$ kpc. Hence such perturbations have widespread influence which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW-Sgr interaction. While Sgr has not yet crossed the disk (currently, $z_{Sgr} \approx -6$ kpc, $v_{z,Sgr} \approx 210$ km/s), we demonstrate that the peak of the impact has already passed. Sgr's pull over the past 150 Myr creates a global $v_z$ signature with amplitude $\propto M_{Sgr}$, which might be detectable in future spectroscopic surveys. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Journal ref: ApJ 928 80 (2022)

Selection Functions in Astronomical Data Modeling, with the Space Density of White Dwarfs as Worked Example

Authors: Hans-Walter Rix, David W. Hogg, Douglas Boubert, Anthony G. A. Brown, Andrew Casey, Ronald Drimmel, Andrew Everall, Morgan Fouesneau, Adrian M. Price-Whelan

Abstract: Statistical studies of astronomical data sets, in particular of cataloged properties for discrete objects, are central to astrophysics. One cannot model those objects' population properties or incidences without a quantitative understanding of the conditions under which these objects ended up in a catalog or sample, the sample's selection function. As systematic and didactic introductions to this… ▽ More Statistical studies of astronomical data sets, in particular of cataloged properties for discrete objects, are central to astrophysics. One cannot model those objects' population properties or incidences without a quantitative understanding of the conditions under which these objects ended up in a catalog or sample, the sample's selection function. As systematic and didactic introductions to this topic are scarce in the astrophysical literature, we aim to provide one, addressing generically the following questions: What is a selection function? What arguments $\vec{q}$ should a selection function depend on? Over what domain must a selection function be defined? What approximations and simplifications can be made? And, how is a selection function used in `modelling'? We argue that volume-complete samples, with the volume drastically curtailed by the faintest objects, reflect a highly sub-optimal selection function that needlessly reduces the number of bright and usually rare objects in the sample. We illustrate these points by a worked example, deriving the space density of white dwarfs (WD) in the Galactic neighbourhood as a function of their luminosity and Gaia color, $Φ_0(M_G,B-R)$ in [mag$^{-2}$pc$^{-3}$]. We construct a sample of $10^5$ presumed WDs through straightforward selection cuts on the Gaia EDR3 catalog, in magnitude, color, parallax, and astrometric fidelity $\vec{q}=(m_G,B-R,\varpi,p_{af})$. We then combine a simple model for $Φ_0$ with the effective survey volume derived from this selection function $S_C(\vec{q})$ to derive a detailed and robust estimate of $Φ_0(M_G,B-R)$. This resulting white dwarf luminosity-color function $Φ_0(M_G,B-R)$ differs dramatically from the initial number density distribution in the luminosity-color plane: by orders of magnitude in density and by four magnitudes in density peak location. △ Less

Submitted 14 June, 2021; originally announced June 2021.

Journal ref: ApJ 2021