skip to main content

Title: Sizing from the smallest scales: the mass of the Milky Way

As the Milky Way and its satellite system become more entrenched in near field cosmology efforts, the need for an accurate mass estimate of the Milky Way’s dark matter halo is increasingly critical. With the second and early third data releases of stellar proper motions from Gaia, several groups calculated full 6D phase-space information for the population of Milky Way satellite galaxies. Utilizing these data in comparison to subhalo properties drawn from the Phat ELVIS simulations, we constrain the Milky Way dark matter halo mass to be ∼1–1.2 × 1012 M⊙. We find that the kinematics of subhaloes drawn from more- or less-massive hosts (i.e. >1.2 × 1012 M⊙ or <1012 M⊙) are inconsistent, at the 3σ confidence level, with the observed velocities of the Milky Way satellites. The preferred host halo mass for the Milky Way is largely insensitive to the exclusion of systems associated with the Large Magellanic Cloud, changes in galaxy formation thresholds, and variations in observational completeness. As more Milky Way satellites are discovered, their velocities (radial, tangential, and total) plus Galactocentric distances will provide further insight into the mass of the Milky Way dark matter halo.

more » « less
Award ID(s):
1945310 1910346 1752913 2138144
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
p. 4968-4982
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We measure the enclosed Milky Way mass profile to Galactocentric distances of ∼70 and ∼50 kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration (SDSS/SEGUE) K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and Gaia DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to ∼70 and ∼50 kpc, the Jeans method yields total enclosed masses of 4.3 ± 0.95 (random) ±0.6 (systematic) × 1011 M⊙ and 4.1 ± 1.2 (random) ±0.6 (systematic) × 1011 M⊙ for the KG and BHB stars, respectively. For the KG and BHB samples, we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) ±0.083 (systematic) × 1012 M⊙ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) ±0.15 (systematic) × 1012 M⊙, respectively.

    more » « less
  2. Abstract

    Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. BelowM∼ 109M, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way (MW). Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume (D≲ 12 Mpc) to date. We use 250 confirmed and 71 candidate dwarf satellites around 27 MW-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semianalyticalSatGenmodel for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMFs), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low-mass end are populated by a relation of the formMMpeakα, with a moderate slope ofαconst=2.10±0.01and a low scatter, constant as a function of the peak halo mass, ofσconst=0.060.05+0.07. A model with a steeper slope (αgrow= 2.39 ± 0.06) and a scatter that grows with decreasingMpeakis also consistent with the observed SSMF but is not required. Our new model for the satellite–subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only MW satellites.

    more » « less
  3. Abstract

    The shape and orientation of dark matter (DM) halos are sensitive to the microphysics of the DM particles, yet in many mass models, the symmetry axes of the Milky Way’s DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well motivated for the inner DM halo, but not for the outer halo. We use zoomed-in cosmological baryonic simulations from the Latte suite of FIRE-2 Milky Way–mass galaxies to explore the evolution of the DM halo’s orientation with radius and time, with or without a major merger with a Large Magellanic Cloud analog, and when varying the DM model. In three of the four cold DM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20° beyond about 30 galactocentric kpc, reaching a maximum of 30°–90°, depending on the individual halo’s formation history. In identical simulations using a model of self-interacting DM withσ= 1 cm2g−1, the halo remains aligned with the stellar disk out to ∼200–400 kpc. Interactions with massive satellites (M≳ 4 × 1010Mat pericenter;M≳ 3.3 × 1010Mat infall) affect the orientation of the halo significantly, aligning the halo’s major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of self-interacting DM, if the effects of massive satellites can be accounted for.

    more » « less
  4. ABSTRACT The satellite populations of the Milky Way, and Milky Way mass galaxies in the local Universe, have been extensively studied to constrain dark matter and galaxy evolution physics. Recently, there has been a shift to studying satellites of hosts with stellar masses between that of the Large Magellanic Cloud and the Milky Way, since they can provide further insight on hierarchical structure formation, environmental effects on satellites, and the nature of dark matter. Most work is focused on the Local Volume, and little is still known about low-mass host galaxies at higher redshift. To improve our understanding of the evolution of satellite populations of low-mass hosts, we study satellite galaxy populations as a function of host stellar mass 9.5 < log (M*/M⊙) < 10.5 and redshifts 0.1 < $z$ < 0.8 in the COSMOS survey, making this the first study of satellite systems of low-mass hosts across half the age of the universe. We find that the satellite populations of low-mass host galaxies, which we measure down to satellite masses equivalent to the Fornax dwarf spheroidal satellite of the Milky Way, remain mostly unchanged through time. We observe a weak dependence between host stellar mass and number of satellites per host, which suggests that the stellar masses of the hosts are in the power-law regime of the stellar mass to halo mass relation (M*–Mhalo) for low-mass galaxies. Finally, we test the constraining power of our measured cumulative luminosity function to calculate the low-mass end slope of the M*–Mhalo relation. These new satellite luminosity function measurements are consistent with Lamda cold dark matter predictions. 
    more » « less

    In this paper, we construct the circular velocity curve of the Milky Way out to ∼30 kpc, providing an updated model of the dark matter density profile. We derive precise parallaxes for 120 309 stars with a data-driven model, using APOGEE DR17 spectra combined with GaiaDR3, 2MASS, and WISE photometry. At outer galactic radii up to 30 kpc, we find a significantly faster decline in the circular velocity curve compared to the inner parts. This decline is better fit with a cored Einasto profile with a slope parameter $0.91^{+0.04}_{-0.05}$ than a generalized Navarro–Frenk–White (NFW) profile. The virial mass of the best-fitting dark matter halo profile is only $1.81^{+0.06}_{-0.05}\times 10^{11}$ M⊙, significantly lower than what a generalized NFW profile delivers. We present a study of the potential systematics, affecting mainly large radii. Such a low mass for the Galaxy is driven by the functional forms tested, given that it probes beyond our measurements. It is found to be in tension with mass measurements from globular clusters, dwarf satellites, and streams. Our best-fitting profile also lowers the expected dark matter annihilation signal flux from the galactic centre by more than an order of magnitude, compared to an NFW profile-fit. In future work, we will explore profiles with more flexible functional forms to more fully leverage the circular velocity curve and observationally constrain the properties of the Milky Way’s dark matter halo.

    more » « less