Showing 1–50 of 184 results for author: Bovy, J

Euclid: Early Release Observations -- Unveiling the morphology of two Milky Way globular clusters out to their periphery

Authors: D. Massari, E. Dalessandro, D. Erkal, E. Balbinot, J. Bovy, I. McDonald, A. M. N. Ferguson, S. S. Larsen, A. Lançon, F. Annibali, B. Goldman, P. B. Kuzma, K. Voggel, T. Saifollahi, J. -C. Cuillandre, M. Schirmer, M. Kluge, B. Altieri, A. Amara, S. Andreon, N. Auricchio, M. Baldi, A. Balestra, S. Bardelli, A. Basset , et al. (136 additional authors not shown)

Abstract: As part of the Euclid Early Release Observations (ERO) programme, we analyse deep, wide-field imaging from the VIS and NISP instruments of two Milky Way globular clusters (GCs), namely NGC 6254 (M10) and NGC 6397, to look for observational evidence of their dynamical interaction with the Milky Way. We search for such an interaction in the form of structural and morphological features in the cluste… ▽ More As part of the Euclid Early Release Observations (ERO) programme, we analyse deep, wide-field imaging from the VIS and NISP instruments of two Milky Way globular clusters (GCs), namely NGC 6254 (M10) and NGC 6397, to look for observational evidence of their dynamical interaction with the Milky Way. We search for such an interaction in the form of structural and morphological features in the clusters' outermost regions, which are suggestive of the development of tidal tails on scales larger than those sampled by the ERO programme. Our multi-band photometric analysis results in deep and well-behaved colour-magnitude diagrams that, in turn, enable an accurate membership selection. The surface brightness profiles built from these samples of member stars are the deepest ever obtained for these two Milky Way GCs, reaching down to $\sim30.0$ mag~arcsec$^{-2}$, which is about $1.5$ mag arcsec$^{-2}$ below the current limit. The investigation of the two-dimensional density map of NGC 6254 reveals an elongated morphology of the cluster peripheries in the direction and with the amplitude predicted by $N$-body simulations of the cluster's dynamical evolution, at high statistical significance. We interpret this as strong evidence for the first detection of tidally induced morphological distortion around this cluster. The density map of NGC 6397 reveals a slightly elliptical morphology, in agreement with previous studies, which requires further investigation on larger scales to be properly interpreted. This ERO project thus demonstrates the power of Euclid in studying the outer regions of GCs at an unprecedented level of detail, thanks to the combination of large field of view, high spatial resolution, and depth enabled by the telescope. Our results highlight the future Euclid survey as the ideal data set to investigate GC tidal tails and stellar streams. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 15 pages, 18 figures. Paper accepted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Programme overview and pipeline for compact- and diffuse-emission photometry

Authors: J. -C. Cuillandre, E. Bertin, M. Bolzonella, H. Bouy, S. Gwyn, S. Isani, M. Kluge, O. Lai, A. Lançon, D. A. Lang, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, Abdurro'uf, N. Aghanim, B. Altieri, F. Annibali, H. Atek, P. Awad, M. Baes, E. Bañados, D. Barrado, S. Belladitta, V. Belokurov , et al. (240 additional authors not shown)

Abstract: The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline t… ▽ More The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Distribution functions for the modelling of accretion remnants in Milky Way-like galaxies: insights from IllustrisTNG

Authors: James M. M. Lane, Jo Bovy

Abstract: We study accretion remnants around Milky Way analogs in the IllustrisTNG simulations to determine how well commonly used distribution functions (DFs) describe their phase-space distributions. We identify 30 Milky Way analogs and 116 remnants from mergers with stellar mass ratios greater than 1:20. Two-power density profiles, as well as rotating constant-anisotropy and Osipkov-Merritt DFs are fit t… ▽ More We study accretion remnants around Milky Way analogs in the IllustrisTNG simulations to determine how well commonly used distribution functions (DFs) describe their phase-space distributions. We identify 30 Milky Way analogs and 116 remnants from mergers with stellar mass ratios greater than 1:20. Two-power density profiles, as well as rotating constant-anisotropy and Osipkov-Merritt DFs are fit to the remnants. We determine that the remnants are suitable for equilibrium modelling by assessing them in the context of the Jeans equation. Each of the models we consider are reasonably able to fit the stellar remnant energy and angular momentum distribution, as well as the magnitude and shape of velocity dispersion profiles. Case studies matched to two well-known merger remnants in the stellar halo-Gaia-Sausage/Enceladus (GS/E) and Sequoia-are explored in more depth. We find good evidence that remnants with high anisotropy $β$, such as GS/E, are better modelled with a superposition of two Osipkov-Merritt DFs than either a constant-anisotropy model or a single Osipkov-Merritt DF. We estimate an Osipkov-Merritt profile with scale radius between 2-4 kpc would be a good first-order representation of GS/E, and comment on existing observational evidence for this as well as studies which could demonstrate it. Overall, we find that DF-based models work well for describing the kinematics of large merger remnants. Our results will be an important reference for future studies which seek to constrain both the spatial and kinematic properties of merger remnants in the Milky Way stellar halo. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 17 pages, 11 figures. Submitted to MNRAS. Comments welcome

Scaling Laws for Galaxy Images

Authors: Mike Walmsley, Micah Bowles, Anna M. M. Scaife, Jason Shingirai Makechemu, Alexander J. Gordon, Annette M. N. Ferguson, Robert G. Mann, James Pearson, Jürgen J. Popp, Jo Bovy, Josh Speagle, Hugh Dickinson, Lucy Fortson, Tobias Géron, Sandor Kruk, Chris J. Lintott, Kameswara Mantha, Devina Mohan, David O'Ryan, Inigo V. Slijepevic

Abstract: We present the first systematic investigation of supervised scaling laws outside of an ImageNet-like context - on images of galaxies. We use 840k galaxy images and over 100M annotations by Galaxy Zoo volunteers, comparable in scale to Imagenet-1K. We find that adding annotated galaxy images provides a power law improvement in performance across all architectures and all tasks, while adding trainab… ▽ More We present the first systematic investigation of supervised scaling laws outside of an ImageNet-like context - on images of galaxies. We use 840k galaxy images and over 100M annotations by Galaxy Zoo volunteers, comparable in scale to Imagenet-1K. We find that adding annotated galaxy images provides a power law improvement in performance across all architectures and all tasks, while adding trainable parameters is effective only for some (typically more subjectively challenging) tasks. We then compare the downstream performance of finetuned models pretrained on either ImageNet-12k alone vs. additionally pretrained on our galaxy images. We achieve an average relative error rate reduction of 31% across 5 downstream tasks of scientific interest. Our finetuned models are more label-efficient and, unlike their ImageNet-12k-pretrained equivalents, often achieve linear transfer performance equal to that of end-to-end finetuning. We find relatively modest additional downstream benefits from scaling model size, implying that scaling alone is not sufficient to address our domain gap, and suggest that practitioners with qualitatively different images might benefit more from in-domain adaption followed by targeted downstream labelling. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 10+6 pages, 12 figures. Appendix C2 based on arxiv:2206.11927. Code, demos, documentation at https://github.com/mwalmsley/zoobot

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

Towards an astronomical foundation model for stars with a Transformer-based model

Authors: Henry W. Leung, Jo Bovy

Abstract: Rapid strides are currently being made in the field of artificial intelligence using Transformer-based models like Large Language Models (LLMs). The potential of these methods for creating a single, large, versatile model in astronomy has not yet been explored. In this work, we propose a framework for data-driven astronomy that uses the same core techniques and architecture as used by LLMs. Using… ▽ More Rapid strides are currently being made in the field of artificial intelligence using Transformer-based models like Large Language Models (LLMs). The potential of these methods for creating a single, large, versatile model in astronomy has not yet been explored. In this work, we propose a framework for data-driven astronomy that uses the same core techniques and architecture as used by LLMs. Using a variety of observations and labels of stars as an example, we build a Transformer-based model and train it in a self-supervised manner with cross-survey data sets to perform a variety of inference tasks. In particular, we demonstrate that a $\textit{single}$ model can perform both discriminative and generative tasks even if the model was not trained or fine-tuned to do any specific task. For example, on the discriminative task of deriving stellar parameters from Gaia XP spectra, we achieve an accuracy of 47 K in $T_\mathrm{eff}$, 0.11 dex in $\log{g}$, and 0.07 dex in $[\mathrm{M/H}]$, outperforming an expert $\texttt{XGBoost}$ model in the same setting. But the same model can also generate XP spectra from stellar parameters, inpaint unobserved spectral regions, extract empirical stellar loci, and even determine the interstellar extinction curve. Our framework demonstrates that building and training a $\textit{single}$ foundation model without fine-tuning using data and parameters from multiple surveys to predict unmeasured observations and parameters is well within reach. Such "Large Astronomy Models" trained on large quantities of observational data will play a large role in the analysis of current and future large surveys. △ Less

Submitted 2 November, 2023; v1 submitted 21 August, 2023; originally announced August 2023.

Decoding the age-chemical structure of the Milky Way disk: An application of Copulas and Elicitable Maps

Authors: Aarya A. Patil, Jo Bovy, Sebastian Jaimungal, Neige Frankel, Henry W. Leung

Abstract: In the Milky Way, the distribution of stars in the $[α/\mathrm{Fe}]$ vs. $[\mathrm{Fe/H}]$ and $[\mathrm{Fe/H}]$ vs. age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disk. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the A… ▽ More In the Milky Way, the distribution of stars in the $[α/\mathrm{Fe}]$ vs. $[\mathrm{Fe/H}]$ and $[\mathrm{Fe/H}]$ vs. age planes holds essential information about the history of star formation, accretion, and dynamical evolution of the Galactic disk. We investigate these planes by applying novel statistical methods called copulas and elicitable maps to the ages and abundances of red giants in the APOGEE survey. We find that the low- and high-$α$ disk stars have a clean separation in copula space and use this to provide an automated separation of the $α$ sequences using a purely statistical approach. This separation reveals that the high-$α$ disk ends at the same [$α$/Fe] and age at high $[\mathrm{Fe/H}]$ as the low-$[\mathrm{Fe/H}]$ start of the low-$α$ disk, thus supporting a sequential formation scenario for the high- and low-$α$ disks. We then combine copulas with elicitable maps to precisely obtain the correlation between stellar age $τ$ and metallicity $[\mathrm{Fe/H}]$ conditional on Galactocentric radius $R$ and height $z$ in the range $0 < R < 20$ kpc and $|z| < 2$ kpc. The resulting trends in the age-metallicity correlation with radius, height, and [$α$/Fe] demonstrate a $\approx 0$ correlation wherever kinematically-cold orbits dominate, while the naively-expected negative correlation is present where kinematically-hot orbits dominate. This is consistent with the effects of spiral-driven radial migration, which must be strong enough to completely flatten the age-metallicity structure of the low-$α$ disk. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 18 pages, 16 figures, Submitted to MNRAS

The stellar mass of the Gaia-Sausage/Enceladus accretion remnant

Authors: James Lane, Jo Bovy, Ted Mackereth

Abstract: The \textit{Gaia}-Sausage/Enceladus (GS/E) structure is an accretion remnant which comprises a large fraction of the Milky Way's stellar halo. We study GS/E using high-purity samples of kinematically selected stars from APOGEE DR16 and \textit{Gaia}. Employing a novel framework to account for kinematic selection biases using distribution functions, we fit density profiles to these GS/E samples and… ▽ More The \textit{Gaia}-Sausage/Enceladus (GS/E) structure is an accretion remnant which comprises a large fraction of the Milky Way's stellar halo. We study GS/E using high-purity samples of kinematically selected stars from APOGEE DR16 and \textit{Gaia}. Employing a novel framework to account for kinematic selection biases using distribution functions, we fit density profiles to these GS/E samples and measure their masses. We find that GS/E has a shallow density profile in the inner Galaxy, with a break between 15--25~kpc beyond which the profile steepens. We also find that GS/E is triaxial, with axis ratios 1:0.55:0.45 (nearly prolate), and the major axis is oriented about 80~degrees from the Sun--Galactic centre line and 16 degrees above the plane. We measure a stellar mass for GS/E of $1.45\,^{+0.92}_{-0.51}\,\mathrm{(stat.)}\,^{+0.13}_{-0.37} \mathrm{(sys.)}\ \times10^{8}$~\Msun. Our mass estimate is lower than others in the literature, a finding we attribute to the excellent purity of the samples we work with. We also fit a density profile to the entire Milky Way stellar halo, finding a mass in the range of $6.7-8.4 \times 10^{8}$~\Msun, and implying that GS/E could make up as little as 15-25~per~cent of the mass of the Milky Way stellar halo. Our lower stellar mass combined with standard stellar-mass-to-halo mass relations implies that GS/E constituted a minor 1:8 mass-ratio merger at the time of its accretion. △ Less

Submitted 23 September, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted in MNRAS

Disentangling Stellar Age Estimates from Galactic Chemodynamical Evolution

Authors: Jeff Shen, Joshua S. Speagle, J. Ted Mackereth, Yuan-Sen Ting, Jo Bovy

Abstract: Stellar ages are key for determining the formation history of the Milky Way, but are difficult to measure precisely. Furthermore, methods that use chemical abundances to infer ages may entangle the intrinsic evolution of stars with the chemodynamical evolution of the Galaxy. In this paper, we present a framework for making probabilistic predictions of stellar ages, and then quantify the contributi… ▽ More Stellar ages are key for determining the formation history of the Milky Way, but are difficult to measure precisely. Furthermore, methods that use chemical abundances to infer ages may entangle the intrinsic evolution of stars with the chemodynamical evolution of the Galaxy. In this paper, we present a framework for making probabilistic predictions of stellar ages, and then quantify the contribution of both stellar evolution and Galactic chemical evolution to those predictions using SHAP values. We apply this interpretable prediction framework to both a simulated Milky Way sample containing stars in a variety of evolutionary stages and an APOGEE-mocked sample of red clump stars. We find that in the former case, stellar evolution is the dominant driver for age estimates, while in the latter case, the more restricted evolutionary information causes the model to proxy ages through the chemical evolution model. We show that as a result of the use of non-intrinsic Galactic chemical information, trends estimated with the predicted ages, such as the age-metallicity relation, can deviate from the truth. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 18 pages, 17 figures, submitted to ApJ

Constraints on the Galactic Centre environment from \textit{Gaia} hypervelocity stars III: Insights on a possible companion to Sgr A*

Authors: Fraser A. Evans, Alexander Rasskazov, Amber Remmelzwaal, Tommaso Marchetti, Alfred Castro-Ginard, Elena Maria Rossi, Jo Bovy

Abstract: We consider a scenario in which Sgr A* is in a massive black hole binary (MBHB) with an as-of-yet undetected supermassive or intermediate-mass black hole companion. Dynamical encounters between this MBHB and single stars in its immediate vicinity would eject hypervelocity stars (HVSs) with velocities beyond the Galactic escape velocity. In this work, we use existing HVS observations to constrain f… ▽ More We consider a scenario in which Sgr A* is in a massive black hole binary (MBHB) with an as-of-yet undetected supermassive or intermediate-mass black hole companion. Dynamical encounters between this MBHB and single stars in its immediate vicinity would eject hypervelocity stars (HVSs) with velocities beyond the Galactic escape velocity. In this work, we use existing HVS observations to constrain for the first time the existence of a companion to Sgr A*. We simulate the ejection of HVSs via the `MBHB slingshot' scenario and show that the population of HVSs detectable today depends strongly on the companion mass and the separation of the MBHB. We demonstrate that the lack of uncontroversial HVS candidates in \textit{Gaia} Data Release 3 places a firm upper limit on the mass of a possible Sgr A* companion. Within one milliparsec of Sgr A*, our results exclude a companion more massive than $1000 \, \mathrm{M_\odot}$. If Sgr A* recently merged with a companion black hole, our findings indicate that unless this companion was less massive than $500 \, \mathrm{M_\odot}$, this merger must have occurred at least $10$ Myr ago. These results complement and improve upon existing independent constraints on a companion to Sgr A* and show that large regions of its parameter space can now be ruled out. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: Submitted to MNRAS. 16 pages, 9 figures. Comments welcome

A variational encoder-decoder approach to precise spectroscopic age estimation for large Galactic surveys

Authors: Henry W. Leung, Jo Bovy, J. Ted Mackereth, Andrea Miglio

Abstract: Constraints on the formation and evolution of the Milky Way Galaxy require multi-dimensional measurements of kinematics, abundances, and ages for a large population of stars. Ages for luminous giants, which can be seen to large distances, are an essential component of studies of the Milky Way, but they are traditionally very difficult to estimate precisely for a large dataset and often require car… ▽ More Constraints on the formation and evolution of the Milky Way Galaxy require multi-dimensional measurements of kinematics, abundances, and ages for a large population of stars. Ages for luminous giants, which can be seen to large distances, are an essential component of studies of the Milky Way, but they are traditionally very difficult to estimate precisely for a large dataset and often require careful analysis on a star-by-star basis in asteroseismology. Because spectra are easier to obtain for large samples, being able to determine precise ages from spectra allows for large age samples to be constructed, but spectroscopic ages are often imprecise and contaminated by abundance correlations. Here we present an application of a variational encoder-decoder on cross-domain astronomical data to solve these issues. The model is trained on pairs of observations from APOGEE and Kepler of the same star in order to reduce the dimensionality of the APOGEE spectra in a latent space while removing abundance information. The low dimensional latent representation of these spectra can then be trained to predict age with just $\sim$ 1,000 precise seismic ages. We demonstrate that this model produces more precise spectroscopic ages ($\sim$ 22% overall, $\sim$ 11% for red-clump stars) than previous data-driven spectroscopic ages while being less contaminated by abundance information (in particular, our ages do not depend on [$α$/M]). We create a public age catalog for the APOGEE DR17 data set and use it to map the age distribution and the age-[Fe/H]-[$α$/M] distribution across the radial range of the Galactic disk. △ Less

Submitted 26 April, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

The Eighteenth Data Release of the Sloan Digital Sky Surveys: Targeting and First Spectra from SDSS-V

Authors: Andrés Almeida, Scott F. Anderson, Maria Argudo-Fernández, Carles Badenes, Kat Barger, Jorge K. Barrera-Ballesteros, Chad F. Bender, Erika Benitez, Felipe Besser, Dmitry Bizyaev, Michael R. Blanton, John Bochanski, Jo Bovy, William Nielsen Brandt, Joel R. Brownstein, Johannes Buchner, Esra Bulbul, Joseph N. Burchett, Mariana Cano Díaz, Joleen K. Carlberg, Andrew R. Casey, Vedant Chandra, Brian Cherinka, Cristina Chiappini, Abigail A. Coker , et al. (129 additional authors not shown)

Abstract: The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM),… ▽ More The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM), including input catalogs and selection functions for their numerous scientific objectives. We describe the production of the targeting databases and their calibration- and scientifically-focused components. DR18 also includes ~25,000 new SDSS spectra and supplemental information for X-ray sources identified by eROSITA in its eFEDS field. We present updates to some of the SDSS software pipelines and preview changes anticipated for DR19. We also describe three value-added catalogs (VACs) based on SDSS-IV data that have been published since DR17, and one VAC based on the SDSS-V data in the eFEDS field. △ Less

Submitted 6 July, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: Accepted to ApJS

Vertical motion in the Galactic disc: unwinding the Snail

Authors: Neige Frankel, Jo Bovy, Scott Tremaine, David W. Hogg

Abstract: The distribution of stars in the Milky Way disc shows a spiral structure--the Snail--in the space of velocity and position normal to the Galactic mid-plane. The Snail appears as straight lines in the vertical frequency--vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the Sn… ▽ More The distribution of stars in the Milky Way disc shows a spiral structure--the Snail--in the space of velocity and position normal to the Galactic mid-plane. The Snail appears as straight lines in the vertical frequency--vertical phase plane when effects from sample selection are removed. Their slope has the dimension of inverse time, with the simplest interpretation being the inverse age of the Snail. Here, we devise and fit a simple model in which the spiral starts as a lopsided perturbation from steady state, that winds up into the present-day morphology. The winding occurs because the vertical frequency decreases with vertical action. We use data from stars in Gaia EDR3 that have measured radial velocities, pruned by simple distance and photometric selection functions. We divide the data into boxels of dynamical invariants (radial action, angular momentum); our model fits the data well in many of the boxels. The model parameters have physical interpretations: one, $A$, is a perturbation amplitude, and one, $t$, is interpretable in the simplest models as the time since the event that caused the Snail. We find trends relating the strength and age to angular momentum: (i) the amplitude $A$ is small at low angular momentum ($<1\,600\mathrm{\,kpc\,km\,s}^{-1}$ or guiding-centre radius $< 7.3\,$kpc), and over a factor of three larger, with strong variations, in the outer disc; (ii) there is no single well-defined perturbation time, with $t$ varying between 0.2 and 0.6 Gyr. Residuals between the data and the model display systematic trends, implying that the data call for more complex models. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: 9 pages, 8 figures. Submitted to MNRAS (constructive comments are welcome)

The origin and fate of the Gaia phase-space snail

Authors: Scott Tremaine, Neige Frankel, Jo Bovy

Abstract: The Gaia snail is a spiral feature in the distribution of solar-neighbourhood stars in position and velocity normal to the Galactic midplane. The snail probably arises from phase mixing of gravitational disturbances that perturbed the disc in the distant past. The most common hypothesis is that the strongest disturbance resulted from a passage of the Sagittarius dwarf galaxy close to the solar nei… ▽ More The Gaia snail is a spiral feature in the distribution of solar-neighbourhood stars in position and velocity normal to the Galactic midplane. The snail probably arises from phase mixing of gravitational disturbances that perturbed the disc in the distant past. The most common hypothesis is that the strongest disturbance resulted from a passage of the Sagittarius dwarf galaxy close to the solar neighbourhood. In this paper we investigate the alternative hypothesis that the snail is created by many small disturbances rather than one large one, that is, by Gaussian noise in the gravitational potential. Probably most of this noise is due to substructures in the dark-matter halo. We show that this hypothesis naturally reproduces most of the properties of the snail. In particular it predicts correctly, with no free parameters, that the apparent age of the snail will be $ \sim 0.5$ Gyr. An important ingredient of this model is that any snail-like feature in the solar neighbourhood, whatever its cause, is erased by scattering from giant molecular clouds or other small-scale structure on a time-scale $\lesssim 1$ Gyr. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: 9 pages, 6 figures. Submitted to MNRAS

Made-to-Measure Modelling of Globular Clusters

Authors: Jeremy J. Webb, Jason A. S. Hunt, Jo Bovy

Abstract: We present the first application of the made-to-measure method for modelling dynamical systems to globular clusters. Through the made-to-measure algorithm, the masses of individual particles within a model cluster are adjusted while the system evolves forward in time via a gravitational $N$-body code until the model cluster is able to reproduce select properties of an observed cluster. The method… ▽ More We present the first application of the made-to-measure method for modelling dynamical systems to globular clusters. Through the made-to-measure algorithm, the masses of individual particles within a model cluster are adjusted while the system evolves forward in time via a gravitational $N$-body code until the model cluster is able to reproduce select properties of an observed cluster. The method is first applied to observations of mock isotropic and anisotropic clusters while fitting against the cluster's three dimensional or projected density profile, density weighted mean-squared velocity profile, or its density profile with individual mean-squared velocity profiles. We find that a cluster's three-dimensional density profile can easily be reproduced by the made-to-measure method, with minor discrepancies in the outer regions if fitting against a cluster's projected surface density or projected kinematic properties. If an observed cluster is anisotropic, only fitting against the cluster's density profile and individual mean-squared velocity profiles will fully recover the full degree of anisotropy. Partial anisotropy can be recovered as long as two kinematic properties are included in the fit. We further apply the method to observations of the Galactic globular cluster M4 and generate a complete six-dimensional representation of the cluster that reproduces observations of its surface density profile, mean-squared proper motion velocity profile, and mean-squared line of sight velocity profile. The M2M method predicts M4 is primarily isotropic with a mass of $9.2 \pm 0.4 \times 10^4\, M_{\odot}$ and a half-mass radius of $3.7 \pm 0.1$ pc. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 11 pages, 10 figures, submitted to MNRAS for publication

On the Fast Track: Rapid construction of stellar stream paths

Authors: Nathaniel Starkman, Jo Bovy, Jeremy J. Webb, Daniela Calvetti, Erkki Somersalo

Abstract: Stellar streams are sensitive probes of the Galactic potential. The likelihood of a stream model given stream data is often assessed using simulations. However, comparing to simulations is challenging when even the stream paths can be hard to quantify. Here we present a novel application of Self-Organizing Maps and first-order Kalman Filters to reconstruct a stream's path, propagating measurement… ▽ More Stellar streams are sensitive probes of the Galactic potential. The likelihood of a stream model given stream data is often assessed using simulations. However, comparing to simulations is challenging when even the stream paths can be hard to quantify. Here we present a novel application of Self-Organizing Maps and first-order Kalman Filters to reconstruct a stream's path, propagating measurement errors and data sparsity into the stream path uncertainty. The technique is Galactic-model independent, non-parametric, and works on phase-wrapped streams. With this technique, we can uniformly analyze and compare data with simulations, enabling both comparison of simulation techniques and ensemble analysis with stream tracks of many stellar streams. Our method is implemented in the public Python package TrackStream, available at https://github.com/nstarman/trackstream. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 16 pages, 11 figures, preprint

Constraining resonant dark matter self-interactions with strong gravitational lenses

Authors: Daniel Gilman, Yi-Ming Zhong, Jo Bovy

Abstract: We devise a method to constrain self-interacting dark matter (SIDM) from observations of quadruply-imaged quasars, and apply it to five self-interaction potentials with a long-range dark force. We consider several SIDM models with an attractive potential that allows for the formation of quasi-bound states, giving rise to resonant features in the cross section localized at particular velocities bel… ▽ More We devise a method to constrain self-interacting dark matter (SIDM) from observations of quadruply-imaged quasars, and apply it to five self-interaction potentials with a long-range dark force. We consider several SIDM models with an attractive potential that allows for the formation of quasi-bound states, giving rise to resonant features in the cross section localized at particular velocities below $50 \ \rm{km} \ \rm{s^{-1}}$. We propose these resonances, which amplify or suppress the cross section amplitude by over an order of magnitude, accelerate or delay the onset of core collapse in low-mass dark matter halos, and derive constraints on the timescale for core collapse for the five interaction potentials we consider. Our data strongly disfavors scenarios in which a majority of halos core collapse, with the strongest constraints obtained for cross section strengths exceeding $100 \ \rm{cm^2} \rm{g}^{-1}$ at relative velocities below $30 \ \rm{km} \ \rm{s^{-1}}$. This work opens a new avenue to explore the vast landscape of possible SIDM theories. △ Less

Submitted 8 May, 2023; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: accepted for publication in PRD

Quantum fluctuations masquerade as halos: Bounds on ultra-light dark matter from quadruply-imaged quasars

Authors: Alexander Laroche, Daniel Gilman, Xinyu Li, Jo Bovy, Xiaolong Du

Abstract: Ultra-light dark matter (ULDM) refers to a class of theories, including ultra-light axions, in which particles with mass $m_ψ < 10^{-20}\, \rm{eV}$ comprise a significant fraction of the dark matter. A galactic scale de Broglie wavelength distinguishes these theories from cold dark matter (CDM), suppressing the overall abundance of structure on sub-galactic scales, and producing wave-like interfer… ▽ More Ultra-light dark matter (ULDM) refers to a class of theories, including ultra-light axions, in which particles with mass $m_ψ < 10^{-20}\, \rm{eV}$ comprise a significant fraction of the dark matter. A galactic scale de Broglie wavelength distinguishes these theories from cold dark matter (CDM), suppressing the overall abundance of structure on sub-galactic scales, and producing wave-like interference phenomena in the density profiles of halos. With the aim of constraining the particle mass, we analyze the flux ratios in a sample of eleven quadruple-image strong gravitational lenses. We account for the suppression of the halo mass function and concentration-mass relation predicted by ULDM theories, and the wave-like fluctuations in the host halo density profile, calibrating the model for the wave interference against numerical simulations of galactic-scale halos. We show that the granular structure of halo density profiles, in particular, the amplitude of the fluctuations, significantly impacts image flux ratios, and therefore inferences on the particle mass derived from these data. We infer relative likelihoods of CDM to ULDM of 8:1, 7:1, 6:1, and 4:1 for particle masses $\log_{10}(m_ψ/\rm{eV})\in[-22.5,-22.25], [-22.25,-22.0],[-22.0,-21.75], [-21.75,-21.5]$, respectively. Repeating the analysis and omitting fluctuations associated with the wave interference effects, we obtain relative likelihoods of CDM to ULDM with a particle mass in the same ranges of 98:1, 48:1, 26:1 and 18:1, highlighting the significant perturbation to image flux ratios associated with the fluctuations. Nevertheless, our results disfavor the lightest particle masses with $m_ψ < 10^{-21.5}\,\rm{eV}$, adding to mounting pressure on ultra-light axions as a viable dark matter candidate. △ Less

Submitted 15 September, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

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

A direct measurement of the distance to the Galactic center using the kinematics of bar stars

Authors: Henry W. Leung, Jo Bovy, J. Ted Mackereth, Jason A. S. Hunt, Richard R. Lane, John C. Wilson

Abstract: The distance to the Galactic center $R_0$ is a fundamental parameter for understanding the Milky Way, because all observations of our Galaxy are made from our heliocentric reference point. The uncertainty in $R_0$ limits our knowledge of many aspects of the Milky Way, including its total mass and the relative mass of its major components, and any orbital parameters of stars employed in chemo-dynam… ▽ More The distance to the Galactic center $R_0$ is a fundamental parameter for understanding the Milky Way, because all observations of our Galaxy are made from our heliocentric reference point. The uncertainty in $R_0$ limits our knowledge of many aspects of the Milky Way, including its total mass and the relative mass of its major components, and any orbital parameters of stars employed in chemo-dynamical analyses. While measurements of $R_0$ have been improving over a century, measurements in the past few years from a variety of methods still find a wide range of $R_0$ being somewhere within $8.0$ to $8.5\,\mathrm{kpc}$. The most precise measurements to date have to assume that Sgr A$^*$ is at rest at the Galactic center, which may not be the case. In this paper, we use maps of the kinematics of stars in the Galactic bar derived from APOGEE DR17 and Gaia EDR3 data augmented with spectro-photometric distances from the \texttt{astroNN} neural-network method. These maps clearly display the minimum in the rotational velocity $v_T$ and the quadrupolar signature in radial velocity $v_R$ expected for stars orbiting in a bar. From the minimum in $v_T$, we measure $R_0 = 8.23 \pm 0.12\,\mathrm{kpc}$. We validate our measurement using realistic $N$-body simulations of the Milky Way. We further measure the pattern speed of the bar to be $Ω_\mathrm{bar} = 40.08\pm1.78\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}$. Because the bar forms out of the disk, its center is manifestly the barycenter of the bar+disc system and our measurement is therefore the most robust and accurate measurement of $R_0$ to date. △ Less

Submitted 26 April, 2022; originally announced April 2022.

The structure of accreted stellar streams

Authors: Yansong Qian, Yumna Arshad, Jo Bovy

Abstract: Many of the Milky Way's globular clusters are likely accreted from satellite galaxies that have long since merged with the Milky Way. When these globular clusters are susceptible to tidal disruption, this process likely starts already inside the parent satellite leading to an early stellar stream within the satellite. When the parent satellite merges with the Milky Way, the globular cluster and it… ▽ More Many of the Milky Way's globular clusters are likely accreted from satellite galaxies that have long since merged with the Milky Way. When these globular clusters are susceptible to tidal disruption, this process likely starts already inside the parent satellite leading to an early stellar stream within the satellite. When the parent satellite merges with the Milky Way, the globular cluster and its pre-merger stellar stream are accreted in a somewhat chaotic process. Here, we investigate the properties of the accreted stream after the merger as we would see it today using a suite of simulations of accretion events. We find that the accretion process leads to a wide range of behaviors, but generally scatters the accreted stream over a wide, two-dimensional area of the sky. The behavior ranges from a set of a few or more well-defined "sub-streams" extending out from the post-merger thin stream by tens of degrees, to more widely dispersed debris over much of the sky, depending on how close to the center of the MilkyWay the merger happened. Using mock Gaia-like observations of the simulated streams, we demonstrate that an accreted-stream component can explain the off-track features observed in the GD-1 stream. Sub-streams can appear like thin tidal streams themselves that are seemingly unassociated with the post-merger stream, raising the possibility that some of the progenitor-less streams observed in the Milky Way are part of a single or a few accreted streams created in an ancient merger event. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: MNRAS, in press

Simulated Bars May Be Shorter But Are Not Slower Than Observed: TNG50 vs. MaNGA

Authors: Neige Frankel, Annalisa Pillepich, Hans-Walter Rix, Vicente Rodriguez-Gomez, Jason Sanders, Jo Bovy, Juna Kollmeier, Norm Murray, Ted Mackereth

Abstract: Galactic bars are prominent dynamical structures within disk galaxies whose size, formation time, strength, and pattern speed influence the dynamical evolution of their hosts galaxies. Yet, their formation and evolution in a cosmological context is not well understood, as cosmological simulation studies have been limited by the classic trade off between simulation volume and resolution. Here we an… ▽ More Galactic bars are prominent dynamical structures within disk galaxies whose size, formation time, strength, and pattern speed influence the dynamical evolution of their hosts galaxies. Yet, their formation and evolution in a cosmological context is not well understood, as cosmological simulation studies have been limited by the classic trade off between simulation volume and resolution. Here we analyze barred disk galaxies in the cosmological magneto-hydrodynamical simulation TNG50 and quantitatively compare the distributions of bar size and pattern speed to those from MaNGA observations at $z=0$. TNG50 galaxies are selected to match the stellar mass and size distributions of observed galaxies, to account for observational selection effects. We find that the high-resolution of TNG50 yields bars with a wide range of pattern speeds (including those with $\geq 40~\mathrm{km}\,\mathrm{s}^{-1}$\,$\mathrm{kpc}^{-1}$) and a mean value of $\sim36~\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$ consistent with observations within $6\,\mathrm{km}\,\mathrm{s}^{-1}$\,$\mathrm{kpc}^{-1}$, in contrast with previous lower-resolution cosmological simulations that produced bars that were too slow. We find, however, that bars in TNG50 are on average $\sim 35\%$ shorter than observed, although this discrepancy may partly reflect remaining inconsistencies in the simulation-data comparison. This leads to higher values of $\mathcal{R} = R_\mathrm{corot}/R_\mathrm{bar}$ in TNG50, but points to simulated bars being `too short' rather than `too slow'. After repeating the analysis on the lower-resolution run of the same simulation (with the same physical model), we qualitatively reproduce the results obtained in previous studies: this implies that, along with physical model variations, numerical resolution effects may explain the previously found `slowness' of simulated bars. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 14 pages, key figures are Fig. 5 and 6. Submitted to ApJ, comments are welcome

The primordial matter power spectrum on sub-galactic scales

Authors: Daniel Gilman, Andrew Benson, Jo Bovy, Simon Birrer, Tommaso Treu, Anna Nierenberg

Abstract: The primordial matter power spectrum quantifies fluctuations in the distribution of dark matter immediately following inflation. Over cosmic time, over-dense regions of the primordial density field grow and collapse into dark matter halos, whose abundance and density profiles retain memory of the initial conditions. By analyzing the image magnifications in eleven strongly-lensed and quadruply-imag… ▽ More The primordial matter power spectrum quantifies fluctuations in the distribution of dark matter immediately following inflation. Over cosmic time, over-dense regions of the primordial density field grow and collapse into dark matter halos, whose abundance and density profiles retain memory of the initial conditions. By analyzing the image magnifications in eleven strongly-lensed and quadruply-imaged quasars, we infer the abundance and concentrations of low-mass halos, and cast the measurement in terms of the amplitude of the primordial matter power spectrum. We anchor the power spectrum on large scales, isolating the effect of small-scale deviations from the $Λ$CDM prediction. Assuming an analytic model for the power spectrum and accounting for several sources of potential systematic uncertainty, including three different models for the halo mass function, we obtain correlated inferences of $\log_{10}\left(P / P_{Λ\rm{CDM}}\right)$, the power spectrum amplitude relative to the predictions of the concordance cosmological model, of $0.0_{-0.4}^{+0.5}$, $0.1_{-0.6}^{+0.7}$, and $0.2_{-0.9}^{+1.0}$ at k = 10, 25 and 50 $\rm{Mpc^{-1}}$ at $68 \%$ confidence, consistent with cold dark matter and single-field slow-roll inflation. △ Less

Submitted 9 March, 2022; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: accepted to MNRAS; includes expanded discussion of results, tests with simulated data, and additional references

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)

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

Authors: Aarya A. Patil, Jo Bovy, Gwendolyn Eadie, Sebastian Jaimungal

Abstract: High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archaeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the compo… ▽ More High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archaeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of $18,933$ giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ${\approx}$10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: $\mathrm{Fe} \lesssim 0.02$ dex; $\mathrm{C} \lesssim 0.03$ dex; $\mathrm{O}, \mathrm{Mg}, \mathrm{Si}, \mathrm{Ni} \lesssim 0.04$ dex; $\mathrm{Ca} \lesssim 0.05$ dex; $\mathrm{N}, \mathrm{Al} \lesssim 0.07$ dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way. △ Less

Submitted 7 January, 2022; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: 30 pages, 16 figures, 5 tables, Accepted to ApJ, Code available at https://github.com/aaryapatil/specdims

Journal ref: The Astrophysical Journal 926 (2022) 51

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

Variation in the Stellar Mass Function Along Stellar Streams

Authors: Jeremy J. Webb, Jo Bovy

Abstract: Stellar streams are the inevitable end product of star cluster evolution, with the properties of a given stream being related to its progenitor. We consider how the dynamical history of a progenitor cluster, as traced by the evolution of its stellar mass function, is reflected in the resultant stream. We generate model streams by evolving star clusters with a range of initial half-mass relaxation… ▽ More Stellar streams are the inevitable end product of star cluster evolution, with the properties of a given stream being related to its progenitor. We consider how the dynamical history of a progenitor cluster, as traced by the evolution of its stellar mass function, is reflected in the resultant stream. We generate model streams by evolving star clusters with a range of initial half-mass relaxation times and dissolution times via direct N-body simulations. Stellar streams that dissolve quickly show no variation in the stellar mass function along the stream. Variation is, however, observed along streams with progenitor clusters that dissolve after several relaxation times. The mass function at the edges of a stream is approximately primordial as it is populated by the first stars to escape the cluster before segregation occurs. Moving inwards the mass function steepens as the intermediate parts of the stream consist of mostly low-mass stars that escaped the cluster after some segregation has occurred. The centre of the stream is then marked by a flatter mass function, as the region is dominated by high-mass stars that quickly segregated to the progenitor cluster's centre and were the last stars to become unbound. We further find that the maximum slope of the mass function along the stream and the rate at which it decreases with distance from the dissolved progenitor serve as proxies for the dynamical state reached by the progenitor cluster before dissolution; this may be able to be applied to observed streams with near-future observations. △ Less

Submitted 24 November, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 11 pages, 11 figures, accepted by MNRAS for publication

Exploring the Sgr-Milky-Way-disc interaction using high resolution N-body simulations

Authors: Morgan Bennett, Jo Bovy, Jason A. S. Hunt

Abstract: The ongoing merger of the Sagittarius (Sgr) dwarf galaxy with the Milky Way is believed to strongly affect the dynamics of the Milky Way's disc. We present a suite of 13 $N$-body simulations, with 500 million to 1 billion particles, modelling the interaction between the Sagittarius dwarf galaxy (Sgr) and the Galactic disc. To quantify the perturbation to the disc's structure and dynamics in the si… ▽ More The ongoing merger of the Sagittarius (Sgr) dwarf galaxy with the Milky Way is believed to strongly affect the dynamics of the Milky Way's disc. We present a suite of 13 $N$-body simulations, with 500 million to 1 billion particles, modelling the interaction between the Sagittarius dwarf galaxy (Sgr) and the Galactic disc. To quantify the perturbation to the disc's structure and dynamics in the simulation, we compute the number count asymmetry and the mean vertical velocity in a solar-neighbourhood-like volume. We find that overall the trends in the simulations match those seen in a simple one-dimensional model of the interaction. We explore the effects of changing the mass model of Sgr, the orbital kinematics of Sgr, and the mass of the Milky Way halo. We find that none of the simulations match the observations of the vertical perturbation using Gaia Data Release 2. In the simulation which is the most similar, we find that the final mass of Sgr far exceeds the observed mass of the Sgr remnant, the asymmetry wavelength is too large, and the shape of the asymmetry doesn't match past $z \approx 0.7$ kpc. We therefore conclude that our simulations support the conclusion that Sgr alone could not have caused the observed perturbation to the solar neighbourhood. △ Less

Submitted 11 March, 2022; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: ApJ, in press

The Kinematic Properties of Milky Way Stellar Halo Populations

Authors: James M. M. Lane, Jo Bovy, J. Ted Mackereth

Abstract: In the $Gaia$ era stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern datasets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo.… ▽ More In the $Gaia$ era stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern datasets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ``action diamond'') is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially-anisotropic $Gaia$-Enceladus merger remnant. Assuming a 1:1 ratio of $Gaia$-Enceladus stars to more isotropic halo, we find our selection achieves a sample purity of 82 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic datasets can noticeably impact the $E-L_{z}$ distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo. △ Less

Submitted 7 October, 2022; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: MNRAS in press

Journal ref: Monthly Notices of the Royal Astronomical Society (2022), Volume 510, Issue 4, pp.5119-5141

No globular cluster progenitors in Milky Way satellite galaxies

Authors: Pierre Boldrini, Jo Bovy

Abstract: In order to find the possible progenitors of Milky Way globular clusters, we perform orbit integrations to track the orbits of 170 Galactic globular clusters and the eleven classical Milky Way satellite galaxies backwards in time for 11 Gyr in a Milky-Way-plus-satellites potential including the response of the MW to the infall of the Large Magellanic Cloud and the effect of dynamical friction on t… ▽ More In order to find the possible progenitors of Milky Way globular clusters, we perform orbit integrations to track the orbits of 170 Galactic globular clusters and the eleven classical Milky Way satellite galaxies backwards in time for 11 Gyr in a Milky-Way-plus-satellites potential including the response of the MW to the infall of the Large Magellanic Cloud and the effect of dynamical friction on the satellites. To evaluate possible past associations, we devise a globular-cluster--satellite binding criterion based on the satellite's tidal radius and escape velocity and we test it on globular clusters associated with the Sagittarius dwarf and on dwarf galaxies associated with the Large Magellanic Cloud. For these, we successfully recover the dynamical associations highlighted by previous studies and we derive their time of accretion by the Galaxy by using Gaia EDR3 data. Assuming that Milky Way globular clusters are and have been free of dark matter and thus consist of stars alone, we demonstrate that none of the globular clusters show any clear association with the eight classical dwarf spheroidal galaxies even though a large fraction of them are believed to be accreted. This means that accreted globular clusters either came in as part of now-disrupted satellite galaxies or that globular clusters may have had dark matter halos in the past -- as suggested by the similar metallicity between globular clusters and dwarf galaxies. △ Less

Submitted 13 September, 2022; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: 9 pages, 3 figures, accepted by MNRAS

Strong lensing signatures of self-interacting dark matter in low-mass halos

Authors: Daniel Gilman, Jo Bovy, Tommaso Treu, Anna Nierenberg, Simon Birrer, Andrew Benson, Omid Sameie

Abstract: Core formation and runaway core collapse in models with self-interacting dark matter (SIDM) significantly alter the central density profiles of collapsed halos. Using a forward modeling inference framework with simulated datasets, we demonstrate that flux ratios in quadruple image strong gravitational lenses can detect the unique structural properties of SIDM halos, and statistically constrain the… ▽ More Core formation and runaway core collapse in models with self-interacting dark matter (SIDM) significantly alter the central density profiles of collapsed halos. Using a forward modeling inference framework with simulated datasets, we demonstrate that flux ratios in quadruple image strong gravitational lenses can detect the unique structural properties of SIDM halos, and statistically constrain the amplitude and velocity dependence of the interaction cross section in halos with masses between $10^6 - 10^{10} M_{\odot}$. Measurements on these scales probe self-interactions at velocities below $30 \ \rm{km} \ \rm{s^{-1}}$, a relatively unexplored regime of parameter space, complimenting constraints at higher velocities from galaxies and clusters. We cast constraints on the amplitude and velocity dependence of the interaction cross section in terms of $σ_{20}$, the cross section amplitude at $20 \ \rm{km} \ \rm{s^{-1}}$. With 50 lenses, a sample size available in the near future, and flux ratios measured from spatially compact mid-IR emission around the background quasar, we forecast $σ_{20} < 11-23 \ \rm{cm^2} \rm{g^{-1}}$ at $95 \%$ CI, depending on the amplitude of the subhalo mass function, and assuming cold dark matter (CDM). Alternatively, if $σ_{20} = 19.2 \ \rm{cm^2}\rm{g^{-1}}$ we can rule out CDM with a likelihood ratio of 20:1, assuming an amplitude of the subhalo mass function that results from doubly-efficient tidal disruption in the Milky Way relative to massive elliptical galaxies. These results demonstrate that strong lensing of compact, unresolved sources can constrain SIDM structure on sub-galactic scales across cosmological distances, and the evolution of SIDM density profiles over several Gyr of cosmic time. △ Less

Submitted 5 September, 2021; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: matches version accepted to MNRAS

On $N$-body simulations of globular cluster stream

Authors: Nilanjan Banik, Jo Bovy

Abstract: Stellar tidal streams are sensitive tracers of the properties of the gravitational potential in which they orbit and detailed observations of their density structure can be used to place stringent constraints on fluctuations in the potential caused by, e.g., the expected populations of dark matter subhalos in the standard cold dark matter paradigm (CDM). Simulations of the evolution of stellar str… ▽ More Stellar tidal streams are sensitive tracers of the properties of the gravitational potential in which they orbit and detailed observations of their density structure can be used to place stringent constraints on fluctuations in the potential caused by, e.g., the expected populations of dark matter subhalos in the standard cold dark matter paradigm (CDM). Simulations of the evolution of stellar streams in live $N$-body halos without low-mass dark-matter subhalos, however, indicate that streams exhibit significant perturbations on small scales even in the absence of substructure. Here we demonstrate, using high-resolution $N$-body simulations combined with sophisticated semi-analytic and simple analytic models, that the mass resolutions of $10^4$--$10^5\,\rm{M}_{\odot}$ commonly used to perform such simulations cause spurious stream density variations with a similar magnitude on large scales as those expected from a CDM-like subhalo population and an order of magnitude larger on small, yet observable, scales. We estimate that mass resolutions of $\approx100\,\rm{M}_{\odot}$ ($\approx1\,\rm{M}_{\odot}$) are necessary for spurious, numerical density variations to be well below the CDM subhalo expectation on large (small) scales. That streams are sensitive to a simulation's particle mass down to such small masses indicates that streams are sensitive to dark matter clustering down to these low masses if a significant fraction of the dark matter is clustered or concentrated in this way, for example, in MACHO models with masses of $10$--$100\,\rm{M}_{\odot}$. △ Less

Submitted 28 January, 2021; originally announced January 2021.

Comments: 6 pages, 6 figures

Using LSST Microlensing to Constrain Dark Compact Objects in Spherical and Disk Configurations

Authors: Harrison Winch, Jack Setford, Jo Bovy, David Curtin

Abstract: The Legacy Survey of Space and Time (LSST) with the Vera Rubin Observatory will provide strong microlensing constraints on dark compact objects (DCOs) in our Galaxy. However, most current forecasts limit their analysis to Primordial Black Holes (PBH). It is unclear how well LSST microlensing will be able to constrain alternative models of DCOs with different Galactic spatial profile distributions… ▽ More The Legacy Survey of Space and Time (LSST) with the Vera Rubin Observatory will provide strong microlensing constraints on dark compact objects (DCOs) in our Galaxy. However, most current forecasts limit their analysis to Primordial Black Holes (PBH). It is unclear how well LSST microlensing will be able to constrain alternative models of DCOs with different Galactic spatial profile distributions at a subdominant DM fraction. In this work, we investigate how well LSST microlensing will constrain spherical or disk-like Galactic spatial distributions of DCOs, taking into account extended observing times, baryonic microlensing background, and sky distribution of LSST sources. These extensions represent significant improvements over existing microlensing forecasts in terms of both accuracy and versatility. We demonstrate this power by deriving new LSST sensitivity projections for DCOs in spherical and disk-like distributions. We forecast that LSST will be able to constrain one solar mass PBHs to have a DM fraction under $4.1\times10^{-4}$. One-solar-mass objects in a dark disk distribution with the same dimensions as the Galactic disk will be constrained below $3.1\times10^{-4}$, while those with $m = 10^5M_{\odot}$ will be constrained to below $3.4\times10^{-5}$. We find that compressed dark disks can be constrained up to a factor of $\sim10$ better than ones with identical dimensions to the baryonic disk. We also find that dark disks become less tightly constrained when they are tilted with respect to our own disk. This forecasting software is a versatile tool, capable of constraining any model of DCOs in the Milky Way with microlensing, and is made publically available at {https://github.com/HarrisonWinch96/DarkDisk_Microlensing}. △ Less

Submitted 13 June, 2022; v1 submitted 13 December, 2020; originally announced December 2020.

Comments: 16 pages, 10 figures

Journal ref: ApJ 933 177 (2022)

A purely acceleration-based measurement of the fundamental Galactic parameters

Authors: Jo Bovy

Abstract: Klioner et al. have used the Gaia EDR3 data to directly measure the solar system's acceleration within the Milky Way using the apparent proper motions of quasars. Here I show that this single absolute acceleration measurement in combination with relative accelerations obtained from pulsar orbital decay allows one to determine all of the parameters describing the dynamics of our local Galactic envi… ▽ More Klioner et al. have used the Gaia EDR3 data to directly measure the solar system's acceleration within the Milky Way using the apparent proper motions of quasars. Here I show that this single absolute acceleration measurement in combination with relative accelerations obtained from pulsar orbital decay allows one to determine all of the parameters describing the dynamics of our local Galactic environment, including the circular velocity at the Sun $V_0 = 244 \pm 8$ km/s and its derivative $V'_0 = 2 \pm 9$ km/s/kpc, the local angular frequency, the Oort constants, and the Sun's motion with respect to the LSR. This is the first determination of these parameters that only uses the general theory of relativity without the need for additional assumptions. △ Less

Submitted 27 July, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Don't cross the streams: caustics from Fuzzy Dark Matter

Authors: Neal Dalal, Jo Bovy, Lam Hui, Xinyu Li

Abstract: We study how tidal streams from globular clusters may be used to constrain the mass of ultra-light dark matter particles, called `fuzzy' dark matter (FDM). A general feature of FDM models is the presence of ubiquitous density fluctuations in bound, virialized dark matter structures, on the scale of the de Broglie wavelength, arising from wave interference in the evolving dark matter distribution.… ▽ More We study how tidal streams from globular clusters may be used to constrain the mass of ultra-light dark matter particles, called `fuzzy' dark matter (FDM). A general feature of FDM models is the presence of ubiquitous density fluctuations in bound, virialized dark matter structures, on the scale of the de Broglie wavelength, arising from wave interference in the evolving dark matter distribution. These time-varying fluctuations can disturb the motions of stars, leading to potentially observable signatures in cold thin tidal streams in our own Galaxy. The study of this effect has been hindered by the difficulty in simulating the FDM wavefunction in Milky Way-sized systems. We present a simple method to evolve realistic wavefunctions in nearly static potentials, that should provide an accurate estimate of this granulation effect. We quantify the impact of FDM perturbations on tidal streams, and show that initially, while stream perturbations are small in amplitude, their power spectra exhibit a sharp cutoff corresponding to the de Broglie wavelength of the FDM potential fluctuations. Eventually, when stream perturbations become nonlinear, fold caustics generically arise that lead to density fluctuations with universal behavior. This erases the signature of the de Broglie wavelength in the stream density power spectrum, but we show that it will still be possible to determine the FDM mass in this regime, by considering the fluctuations in quantities like angular momenta or actions. △ Less

Submitted 26 November, 2020; originally announced November 2020.

Comments: Comments welcome. To be submitted to JCAP

Testing the chemical homogeneity of chemically tagged dissolved birth clusters

Authors: Chloe M. Cheng, Natalie Price-Jones, Jo Bovy

Abstract: Chemically tagging stars back to common formation sites in the Milky Way and establishing a high level of chemical homogeneity in these chemically tagged birth clusters is crucial for understanding the chemical and dynamical history of the Galactic disc. We constrain the intrinsic abundance scatter in 17 newly chemically tagged dissolved birth clusters found in the APOGEE survey by modelling APOGE… ▽ More Chemically tagging stars back to common formation sites in the Milky Way and establishing a high level of chemical homogeneity in these chemically tagged birth clusters is crucial for understanding the chemical and dynamical history of the Galactic disc. We constrain the intrinsic abundance scatter in 17 newly chemically tagged dissolved birth clusters found in the APOGEE survey by modelling APOGEE spectra as a one-dimensional function of initial stellar mass, performing forward modelling of the observed stellar spectra, then comparing the data and simulations using Approximate Bayesian Computation. We exemplify this method with the well-known open clusters M67, NGC 6819, and NGC 6791. We study 15 elements measured by APOGEE and find that, in general, we are able to obtain very strong constraints on the intrinsic abundance scatter of most elements in the chemically tagged birth clusters, with upper limits of $\lesssim 0.02$ dex for C, $\lesssim 0.03$ dex for O, Mn, and Fe, $\lesssim 0.04$ dex for Si and Ni, and $\lesssim 0.05$ dex for N, Mg, and Ca. While we find some evidence for a small amount of chemical inhomogeneity in the remaining elements (i.e. Na, Al, S, K, Ti, and V), we are still able to obtain similar or stronger limits compared to those found for open clusters, consistent with previous findings. By strongly constraining their level of chemical homogeneity, we can strengthen the statement that these groups of stars represent birth clusters, with promising implications for future chemical tagging studies. △ Less

Submitted 4 August, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: MNRAS, in press

Did Sgr cause the vertical waves in the solar neighbourhood?

Authors: Morgan Bennett, Jo Bovy

Abstract: The vertical distribution of stars in the solar neighbourhood is not in equilibrium but contains a wave signature in both density and velocity space originating from a perturbation. With the discovery of the phase-space spiral in Gaia data release 2, determining the origin of this perturbation has become even more urgent. We develop and test a fast method for calculating the perturbation from a pa… ▽ More The vertical distribution of stars in the solar neighbourhood is not in equilibrium but contains a wave signature in both density and velocity space originating from a perturbation. With the discovery of the phase-space spiral in Gaia data release 2, determining the origin of this perturbation has become even more urgent. We develop and test a fast method for calculating the perturbation from a passing satellite on the vertical component of a part of a disc galaxy. This fast method allows us to test a large variety of possible perturbations to the vertical disc very quickly. We apply our method to the range of possible perturbations to the solar neighbourhood stemming from the recent passage of the Sagittarius dwarf galaxy (Sgr), varying its mass, mass profile, and present-day position within their observational uncertainties, and its orbit within different realistic models for the Milky Way's gravitational potential. We find that we are unable to reproduce the observed asymmetry in the vertical number counts and its concomitant breathing mode in velocity space for any plausible combination of Sgr and Milky-Way properties. In all cases, either the amplitude or the perturbation wavelength of the number-count asymmetry and of the oscillations in the mean vertical velocity produced by the passage of Sgr is in large disagreement with the observations from Gaia DR2. We conclude that Sgr cannot have caused the observed oscillations in the vertical disc or the Gaia phase-space spiral. △ Less

Submitted 14 June, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

The contribution of N-rich stars to the Galactic stellar halo using APOGEE red giants

Authors: Danny Horta, J. Ted Mackereth, Ricardo P. Schiavon, Sten Hasselquist, Jo Bovy, Carlos Allende Prieto, Timothy C. Beers, Katia Cunha, D. A. García-Hernández, Shobhit S. Kisku, Richard R. Lane, Steven R. Majewski, Andrew C. Mason, David M. Nataf, Alexandre Roman-Lopes, Mathias Schultheis

Abstract: The contribution of dissolved globular clusters (GCs) to the stellar content of the Galactic halo is a key constraint on models for GC formation and destruction, and the mass assembly history of the Milky Way. Earlier results from APOGEE pointed to a large contribution of destroyed GCs to the stellar content of the inner halo, by as much as 25$\%$, which is an order of magnitude larger than previo… ▽ More The contribution of dissolved globular clusters (GCs) to the stellar content of the Galactic halo is a key constraint on models for GC formation and destruction, and the mass assembly history of the Milky Way. Earlier results from APOGEE pointed to a large contribution of destroyed GCs to the stellar content of the inner halo, by as much as 25$\%$, which is an order of magnitude larger than previous estimates for more distant regions of the halo. We set out to measure the ratio between N-rich and normal halo field stars, as a function of distance, by performing density modelling of halo field populations in APOGEE DR16. Our results show that at 1.5 kpc from the Galactic Centre, N-rich stars contribute a much higher 16.8$^{+10.0}_{-7.0}$$\%$ fraction to the total stellar halo mass budget than the 2.7$^{+1.0}_{-0.8}$$\%$ ratio contributed at 10 kpc. Under the assumption that N-rich stars are former GC members that now reside in the stellar halo field, and assuming the ratio between first-and second-population GC stars being 1:2, we estimate a total contribution from disrupted GC stars of the order of 27.5$^{+15.4}_{-11.5}$$\%$ at r = 1.5 kpc and 4.2$^{+1.5}_{-1.3}$$\%$ at r = 10 kpc. Furthermore, since our methodology requires fitting a density model to the stellar halo, we integrate such density within a spherical shell from 1.5-15 kpc in radius, and find a total stellar mass arising from dissolved and/or evaporated GCs of $M_{\mathrm{GC,total}}$ = 9.6$^{+4.0}_{-2.6}$ $\times$ 10$^{7}$ M$\odot$. △ Less

Submitted 12 October, 2022; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: Paper accepted for Publication in MNRAS, correction of Fig 7

The Effects of Dwarf Galaxies on the Orbital Evolution of Galactic Globular Clusters

Authors: Turner Garrow, Jeremy J. Webb, Jo Bovy

Abstract: We investigate the effect that dwarf galaxies have on the orbits, tidal histories, and assumed formation environment of Milky Way globular clusters. We determine the orbits of the Milky Way's 150 globular clusters in a gravitational potential both with and without dwarf galaxies. We find that the presence of a small number of satellite galaxies can affect the orbits of many of the globular cluster… ▽ More We investigate the effect that dwarf galaxies have on the orbits, tidal histories, and assumed formation environment of Milky Way globular clusters. We determine the orbits of the Milky Way's 150 globular clusters in a gravitational potential both with and without dwarf galaxies. We find that the presence of a small number of satellite galaxies can affect the orbits of many of the globular clusters. Over 12 Gyr, we find that the semi-major axis and orbital eccentricity of individual clusters fluctuate with dispersions on the order of $\sim 10\%$ and $\sim 4\%$, respectively. Outer clusters are more strongly affected by dwarf galaxies than inner clusters, with their semi-major axis and orbital eccentricities fluctuating by more than $\sim 15\%$ and $\sim 5\%$, respectively. Using detailed $N$-body simulations of select clusters, we find that altering their orbital histories can lead to different mass loss rates and structural evolution. Furthermore, we caution against using kinematics alone to identify whether a Galactic cluster formed in-situ or was accreted during a past merger event as these values are no longer conserved. The presence of dwarf galaxies causes the orbital energies and actions of individual clusters to evolve over time, spanning a wider range than that coming from random uncertainties in a cluster's proper motions and radial velocity. △ Less

Submitted 20 December, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 10 pages, 6 figures, 2 tables, Accepted to MNRAS

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological Implications from two Decades of Spectroscopic Surveys at the Apache Point observatory

Authors: eBOSS Collaboration, Shadab Alam, Marie Aubert, Santiago Avila, Christophe Balland, Julian E. Bautista, Matthew A. Bershady, Dmitry Bizyaev, Michael R. Blanton, Adam S. Bolton, Jo Bovy, Jonathan Brinkmann, Joel R. Brownstein, Etienne Burtin, Solene Chabanier, Michael J. Chapman, Peter Doohyun Choi, Chia-Hsun Chuang, Johan Comparat, Andrei Cuceu, Kyle S. Dawson, Axel de la Macorra, Sylvain de la Torre, Arnaud de Mattia, Victoria de Sainte Agathe , et al. (75 additional authors not shown)

Abstract: We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Ly$α$ forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter dist… ▽ More We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Ly$α$ forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, $r_d$, from eight different samples and six measurements of the growth rate parameter, $fσ_8$, from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, $Λ$CDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization the BAO data provide nearly an order of magnitude improvement on curvature constraints. The RSD measurements indicate a growth rate that is consistent with predictions from Planck primary data and with General Relativity. When combining the results of SDSS BAO and RSD with external data, all multiple-parameter extensions remain consistent with a $Λ$CDM model. Regardless of cosmological model, the precision on $Ω_Λ$, $H_0$, and $σ_8$, remains at roughly 1\%, showing changes of less than 0.6\% in the central values between models. The inverse distance ladder measurement under a o$w_0w_a$CDM yields $H_0= 68.20 \pm 0.81 \, \rm km\, s^{-1} Mpc^{-1}$, remaining in tension with several direct determination methods. (abridged) △ Less

Submitted 9 July, 2024; v1 submitted 17 July, 2020; originally announced July 2020.

Comments: A summary of all SDSS BAO and RSD measurements with accompanying legacy figures can be found here: https://www.sdss.org/science/final-bao-and-rsd-measurements/ . The full cosmological interpretation of these measurements can be found here: https://www.sdss.org/science/cosmology-results-from-eboss/

Journal ref: Phys. Rev. D 103, 083533 (2021)

APOGEE Data and Spectral Analysis from SDSS Data Release 16: Seven Years of Observations Including First Results from APOGEE-South

Authors: Henrik Jönsson, Jon A. Holtzman, Carlos Allende Prieto, Katia Cunha, D. A. García-Hernández, Sten Hasselquist, Thomas Masseron, Yeisson Osorio, Matthew Shetrone, Verne Smith, Guy S. Stringfellow, Dmitry Bizyaev, Bengt Edvardsson, Steven R. Majewski, Szabolcs Mészáros, Diogo Souto, Olga Zamora, Rachael L. Beaton, Jo Bovy, John Donor, Marc H. Pinsonneault, Vijith Jacob Poovelil, Jennifer Sobeck

Abstract: The spectral analysis and data products in Data Release 16 (DR16; December 2019) from the high-resolution near-infrared APOGEE-2/SDSS-IV survey are described. Compared to the previous APOGEE data release (DR14; July 2017), APOGEE DR16 includes about 200000 new stellar spectra, of which 100000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observato… ▽ More The spectral analysis and data products in Data Release 16 (DR16; December 2019) from the high-resolution near-infrared APOGEE-2/SDSS-IV survey are described. Compared to the previous APOGEE data release (DR14; July 2017), APOGEE DR16 includes about 200000 new stellar spectra, of which 100000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observatory in Chile. DR16 includes all data taken up to August 2018, including data released in previous data releases. All of the data have been re-reduced and re-analyzed using the latest pipelines, resulting in a total of 473307 spectra of 437445 stars. Changes to the analysis methods for this release include, but are not limited to, the use of MARCS model atmospheres for calculation of the entire main grid of synthetic spectra used in the analysis, a new method for filling "holes" in the grids due to unconverged model atmospheres, and a new scheme for continuum normalization. Abundances of the neutron capture element Ce are included for the first time. A new scheme for estimating uncertainties of the derived quantities using stars with multiple observations has been applied, and calibrated values of surface gravities for dwarf stars are now supplied. Compared to DR14, the radial velocities derived for this release more closely match those in the Gaia DR2 data base, and a clear improvement in the spectral analysis of the coolest giants can be seen. The reduced spectra as well as the result of the analysis can be downloaded using links provided in the SDSS DR16 web page. △ Less

Submitted 10 July, 2020; originally announced July 2020.

And In The Darkness Unbind Them: High-Resolution Simulations of Dark Matter Subhalo Disruption in a Milky Way-like Tidal Field

Authors: Jeremy J. Webb, Jo Bovy

Abstract: We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a… ▽ More We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a basis for our simulations of subhalos in a dark matter-only tidal field, while an observationally motivated model for the Milky Way is used for the tidal field that is comprised of a dark matter halo, a stellar disk, and a stellar bulge. Our simulations indicate that including stellar components in the tidal field results in the number of subhalos in Milky Way-like galaxies being only $65\%$ of what is predicted by $Λ$ Cold Dark Matter ($Λ$CDM). For subhalos with small pericentres $(r_p \lesssim 25$ kpc), the subhalo abundance is reduced further to $40\%$, with the surviving subhalos being less dense and having a tangentially-anisotropic orbital distribution. Conversely, subhalos with larger pericentres are minimally affected by the inclusion of a stellar component in the tidal field, with the total number of outer subhalos $\approx 75\%$ of the $Λ$CDM prediction. The densities of outer subhalos are comparable to predictions from $Λ$CDM, with the subhalos having an isotropic distribution of orbits. These ratios are higher than those found in previous studies that include the effects baryonic matter, which are affected by spurious disruption caused by low resolution. △ Less

Submitted 16 December, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 14 pages, 14 figures, 1 table, Accepted for publication in MNRAS

Strong chemical tagging with APOGEE: 21 candidate star clusters that have dissolved across the Milky Way disc

Authors: Natalie Price-Jones, Jo Bovy, Jeremy J. Webb, Carlos Allende Prieto, Rachael Beaton, Joel R. Brownstein, Roger E. Cohen, Katia Cunha, John Donor, Peter M. Frinchaboy, D. A. García-Hernández, Richard R. Lane, Steven R. Majewski, David L. Nidever, Alexandre Roman-Lopes

Abstract: Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe]… ▽ More Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups. △ Less

Submitted 26 June, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: 15 pages, 9 figures, accepted by MNRAS

The building blocks of the Milky Way halo using APOGEE and Gaia -- or -- Is the Galaxy a typical galaxy?

Authors: Ricardo P. Schiavon, J. Ted Mackereth, Joel Pfeffer, Rob A. Crain, Jo Bovy

Abstract: We summarise recent results from analysis of APOGEE/Gaia data for stellar populations in the Galactic halo, disk, and bulge, leading to constraints on the contribution of dwarf galaxies and globular clusters to the stellar content of the Milky Way halo. Interpretation of the extant data in light of cosmological numerical simulations suggests that the Milky Way has been subject to an unusually inte… ▽ More We summarise recent results from analysis of APOGEE/Gaia data for stellar populations in the Galactic halo, disk, and bulge, leading to constraints on the contribution of dwarf galaxies and globular clusters to the stellar content of the Milky Way halo. Interpretation of the extant data in light of cosmological numerical simulations suggests that the Milky Way has been subject to an unusually intense accretion history at z >~ 1.5. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: To Appear in "Star Clusters: From the Milky Way to the Early Universe ", Proceedings of IAU Symposium No. 351, 2019, A. Bragaglia, M.B. Davies, A. Sills & E. Vesperini, eds. Four pages, no figures

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

Authors: Romina Ahumada, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Riccardo Arcodia, Eric Armengaud, Marie Aubert, Santiago Avila, Vladimir Avila-Reese, Carles Badenes, Christophe Balland, Kat Barger, Jorge K. Barrera-Ballesteros, Sarbani Basu, Julian Bautista, Rachael L. Beaton, Timothy C. Beers, B. Izamar T. Benavides, Chad F. Bender, Mariangela Bernardi, Matthew Bershady, Florian Beutler , et al. (289 additional authors not shown)

Abstract: This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the… ▽ More This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17). △ Less

Submitted 11 May, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: DR16 release: Monday Dec 9th 2019. This is the alphabetical order SDSS-IV collaboration data release paper. 25 pages, 6 figures, accepted by ApJS on 11th May 2020. Minor changes clarify or improve text and figures relative to v1

Novel constraints on the particle nature of dark matter from stellar streams

Authors: Nilanjan Banik, Jo Bovy, Gianfranco Bertone, Denis Erkal, T. J. L. de Boer

Abstract: Tidal streams are highly sensitive to perturbations from passing dark matter (DM) subhalos and thus provide a means of measuring their abundance. In a recent paper, we analyzed the distribution of stars along the GD-1 stream with a combination of data from the Gaia satellite and the Pan-STARRS survey, and we demonstrated that the population of DM subhalos predicted by the cold dark matter (CDM) pa… ▽ More Tidal streams are highly sensitive to perturbations from passing dark matter (DM) subhalos and thus provide a means of measuring their abundance. In a recent paper, we analyzed the distribution of stars along the GD-1 stream with a combination of data from the Gaia satellite and the Pan-STARRS survey, and we demonstrated that the population of DM subhalos predicted by the cold dark matter (CDM) paradigm are necessary and sufficient to explain the perturbations observed in the linear density of stars. In this paper, we use the measurements of the subhalo mass function (SHMF) from the GD-1 data combined with a similar analysis of the Pal 5 stream to provide novel constraints on alternative DM scenarios that predict a suppression of the SHMF on scales smaller than the mass of dwarf galaxies, marginalizing over uncertainties in the slope and normalization of the unsuppressed SHMF and the susceptibility of DM subhalos in the inner Milky Way to tidal disruption. In particular, we derive a 95% lower limit on the mass of warm dark matter (WDM) thermal relics $m_{\rm WDM}>3.6\,\mathrm{keV}$ from streams alone that strengthens to $m_{\rm WDM}>6.2\,\mathrm{keV}$ when adding dwarf satellite counts. Similarly, we constrain the axion mass in ultra-light ("fuzzy") dark matter (FDM) models to be $m_{\rm FDM}>1.4\times10^{-21}\,\mathrm{eV}$ from streams alone or $m_{\rm FDM}>2.2\times10^{-21}\,\mathrm{eV}$ when adding dwarf satellite counts. Because we make use of simple approximate forms of the streams' SHMF measurement, our analysis is easy to replicate with other alternative DM models that lead to a suppression of the SHMF. △ Less

Submitted 27 September, 2021; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: JCAP, in press. Code available to fully replicate the analysis at https://github.com/jobovy/dark-matter-constraints-from-stellar-streams

Evidence of a population of dark subhalos from Gaia and Pan-STARRS observations of the GD-1 stream

Authors: Nilanjan Banik, Jo Bovy, Gianfranco Bertone, Denis Erkal, T. J. L. de Boer

Abstract: New data from the $\textit{Gaia}$ satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyze the stellar density perturbations in the GD-1 stream and show that they cannot be due to known baryonic structures like giant molecular clouds, globular clusters, or the Milky Way's bar or spiral arms. A… ▽ More New data from the $\textit{Gaia}$ satellite, when combined with accurate photometry from the Pan-STARRS survey, allow us to accurately estimate the properties of the GD-1 stream. Here, we analyze the stellar density perturbations in the GD-1 stream and show that they cannot be due to known baryonic structures like giant molecular clouds, globular clusters, or the Milky Way's bar or spiral arms. A joint analysis of the GD-1 and Pal 5 streams instead requires a population of dark substructures with masses $\approx 10^{7}$ to $10^9 \ M_{\rm{\odot}}$. We infer a total abundance of dark subhalos normalised to standard cold dark matter $n_{\rm sub}/n_{\rm sub, CDM} = 0.4 ^{+0.3}_{-0.2}$ ($68 \%$), which corresponds to a mass fraction contained in the subhalos $f_{\rm{sub}} = 0.14 ^{+0.11}_{-0.07} \%$, compatible with the predictions of hydrodynamical simulation of cold dark matter with baryons. △ Less

Submitted 21 January, 2021; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: 20 pages, 14 figures. Also see companion paper : "Novel constraints on the particle nature of dark matter from stellar streams". Accepted in MNRAS. Modified main text slightly. Updated Appendix B and added a plot. Main results and conclusions remain unchanged

LRP2020: The cosmic origin and evolution of the elements

Authors: Rodrigo Fernández, Falk Herwig, Samar Safi-Harb, Iris Dillmann, Kim A. Venn, Benoit Côté, Craig O. Heinke, Erik Rosolowsky, Tyrone E. Woods, Daryl Haggard, Luis Lehner, John J. Ruan, Daniel M. Siegel, Jo Bovy, Alan A. Chen, Andrew Cumming, Barry Davids, Maria R. Drout, Reiner Krüecken

Abstract: The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear… ▽ More The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization. △ Less

Submitted 21 October, 2019; originally announced October 2019.

Comments: White paper submitted to the Canadian Long Range Plan 2020. Minor formatting changes relative to submitted version

LRP2020: Astrostatistics in Canada

Authors: Gwendolyn Eadie, Arash Bahramian, Pauline Barmby, Radu Craiu, Derek Bingham, Renée Hložek, JJ Kavelaars, David Stenning, Samantha Benincasa, Guillaume Thomas, Karun Thanjavur, Jo Bovy, Jan Cami, Ray Carlberg, Sam Lawler, Adrian Liu, Henry Ngo, Mubdi Rahman, Michael Rupen

Abstract: (Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summe… ▽ More (Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada. △ Less

Submitted 19 October, 2019; originally announced October 2019.

Comments: White paper E017 submitted to the Canadian Long Range Plan LRP2020

Weighing the stellar constituents of the Galactic halo with APOGEE red giant stars

Authors: J. Ted Mackereth, Jo Bovy

Abstract: The stellar mass in the halo of the Milky Way is notoriously difficult to determine, owing to the paucity of its stars in the solar neighbourhood. With tentative evidence from $\it{Gaia}$ that the nearby stellar halo is dominated by a massive accretion event -- referred to as $\it{Gaia-Enceladus}$ or Sausage -- these constraints are now increasingly urgent. We measure the mass in kinematically sel… ▽ More The stellar mass in the halo of the Milky Way is notoriously difficult to determine, owing to the paucity of its stars in the solar neighbourhood. With tentative evidence from $\it{Gaia}$ that the nearby stellar halo is dominated by a massive accretion event -- referred to as $\it{Gaia-Enceladus}$ or Sausage -- these constraints are now increasingly urgent. We measure the mass in kinematically selected mono-abundance populations (MAPs) of the stellar halo between $-3 <$ [Fe/H] $< -1$ and $0.0 <$ [Mg/Fe] $< 0.4$ using red giant star counts from APOGEE DR14. We find that MAPs are well fit by single power laws on triaxial ellipsoidal surfaces, and we show that that the power law slope $α$ changes such that high $\mathrm{[Mg/Fe]}$ populations have $α\sim 4$, whereas low [Mg/Fe] MAPs are more extended with shallow slopes, $α\sim 2$. We estimate the total stellar mass to be $M_{*,\mathrm{tot}} = 1.3^{+0.3}_{-0.2}\times10^{9}\ \mathrm{M_{\odot}}$, of which we estimate $\sim 0.9^{+0.2}_{-0.1} \times 10^{9} \mathrm{M_{\odot}}$ to be accreted. We estimate that the mass of accreted stars with $e > 0.7$ is $M_{*,\mathrm{accreted},e>0.7} = 3 \pm 1\ \mathrm{(stat.)}\pm 1\ \mathrm{(syst.)}\times10^{8}\ \mathrm{M_\odot}$, or $\sim 30-50\%$ of the accreted halo mass. If the majority of these stars $\it{are}$ the progeny of a massive accreted dwarf, this places an upper limit on its stellar mass, and implies a halo mass for the progenitor of $\sim 10^{10.2\pm0.2}\ \mathrm{M_\odot}$. This constraint not only shows that the $\it{Gaia-Enceladus}$/Sausage progenitor may not be as massive as originally suggested, but that the majority of the Milky Way stellar halo was accreted. These measurements are an important step towards fully reconstructing the assembly history of the Milky Way. △ Less

Submitted 7 January, 2020; v1 submitted 8 October, 2019; originally announced October 2019.

Comments: 17 Pages, 11 Figures, Accepted for publication in MNRAS following minor revisions. Code for reproduction of results is available at https://github.com/jmackereth/halo-mass

LRP2020: Probing Diverse Phenomena through Data-Intensive Astronomy

Authors: Mubdi Rahman, Dustin Lang, Renée Hložek, Jo Bovy, Laurence Perreault-Levasseur

Abstract: The era of data-intensive astronomy is being ushered in with the increasing size and complexity of observational data across wavelength and time domains, the development of algorithms to extract information from this complexity, and the computational power to apply these algorithms to the growing repositories of data. Data-intensive approaches are pushing the boundaries of nearly all fields of ast… ▽ More The era of data-intensive astronomy is being ushered in with the increasing size and complexity of observational data across wavelength and time domains, the development of algorithms to extract information from this complexity, and the computational power to apply these algorithms to the growing repositories of data. Data-intensive approaches are pushing the boundaries of nearly all fields of astronomy, from exoplanet science to cosmology, and they are becoming a critical modality for how we understand the universe. The success of these approaches range from the discovery of rare or unexpected phenomena, to characterizing processes that are now accessible with precision astrophysics and a deep statistical understanding of the datasets, to developing algorithms that maximize the science that can be extracted from any set of observations. In this white paper, we propose a number of initiatives to maximize Canada's ability to compete in this data-intensive era. We propose joining international collaborations and leveraging Canadian facilities for legacy data potential. We propose continuing to build a more agile computing infrastructure that's responsive to the needs of tackling larger and more complex data, as well as enabling quick prototyping and scaling of algorithms. We recognize that developing the fundamental skills of the field will be critical for Canadian astronomers, and discuss avenues through with the appropriate computational and statistical training could occur. Finally, we note that the transition to data-intensive techniques is not limited to astronomy, and we should coordinate with other disciplines to develop and make use of best practises in methods, infrastructure, and education. △ Less

Submitted 4 October, 2019; originally announced October 2019.

Comments: A white paper submitted for consideration for the Canadian Long Range Plan 2020 (E075)