skip to main content


Title: Milky Way Satellite Census. IV. Constraints on Decaying Dark Matter from Observations of Milky Way Satellite Galaxies
Abstract

We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetimeτcomparable to the age of the universe to a lighter DM particle (with mass splittingϵ) and to a dark radiation species. These decays impart a characteristic “kick velocity,”Vkick=ϵc, on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function ofτandVkickusing a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy–halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models withτ< 18 Gyr (29 Gyr) forVkick= 20 kms−1(40 kms−1) at 95% confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble andS8tensions.

 
more » « less
Award ID(s):
2009441
NSF-PAR ID:
10368238
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;   « less
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
932
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 128
Size(s):
["Article No. 128"]
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract

    We develop a linear perturbative formalism to compute the response of an inhomogeneous stellar disk embedded in a nonresponsive dark matter (DM) halo to various perturbations like bars, spiral arms, and encounters with satellite galaxies. Without self-gravity to reinforce it, the response of a Fourier mode phase mixes away due to an intrinsic spread in the vertical (Ωz), radial (Ωr), and azimuthal (Ωϕ) frequencies, triggering local phase-space spirals. The detailed galactic potential dictates the shape of phase spirals: phase mixing occurs more slowly and thus phase spirals are more loosely wound in the outer disk and in the presence of an ambient DM halo. Collisional diffusion due to scattering of stars by structures like giant molecular clouds causes superexponential damping of the phase spiral amplitude. Thezvzphase spiral is one-armed (two-armed) for vertically antisymmetric (symmetric) bending (breathing) modes. Only transient perturbations with timescales (τP) comparable to the vertical oscillation period (τz∼ 1/Ωz) can trigger vertical phase spirals. Each (n,l,m) mode of the response to impulsive (τP<τ= 1/(nΩz+lΩr+mΩϕ)) perturbations is power-law (∼τP/τ) suppressed, but that to adiabatic (τP>τ) perturbations is exponentially weak (expτP/τα) except for resonant (τ→ ∞ ) modes. Slower (τP>τz) perturbations, e.g., distant encounters with satellite galaxies, induce stronger bending modes. Sagittarius (Sgr) dominates the solar neighborhood response of the Milky Way (MW) disk to satellite encounters. Thus, if the Gaia phase spiral was triggered by a MW satellite, Sgr is the leading contender. However, the survival of the phase spiral against collisional damping necessitates an impact ∼0.6–0.7 Gyr ago.

     
    more » « less
  4. Abstract

    To better understand the formation of large, low-surface-brightness galaxies, we measure the correlation function between ultradiffuse galaxy (UDG) candidates and Milky Way analogs (MWAs). We find that: (1) the projected radial distribution of UDG satellites (projected surface density ∝r−0.84±0.06) is consistent with that of normal satellite galaxies; (2) the number of UDG satellites per MWA (SUDG) is ∼0.5 ± 0.1 over projected radii from 20 to 250 kpc and −17 <Mr< −13.5; (3)SUDGis consistent with a linear extrapolation of the relationship between the number of UDGs per halo versus halo mass obtained over galaxy group and cluster scales; (4) red UDG satellites dominate the population of UDG satellites (∼80%); (5) over the range of satellite magnitudes studied, UDG satellites comprise ∼10% of the satellite galaxy population of MWAs; and (6) a significant fraction of these (∼13%) have estimated total masses >1010.9Mor, equivalently, at least half the halo mass of the LMC, and populate a large fraction (∼18%) of the expected subhalos down to these masses. All of these results suggest a close association between the overall low-mass galaxy population and UDGs, which we interpret as favoring models where UDG formation principally occurs within the general context of low-mass galaxy formation over models invoking more exotic physical processes specifically invoked to form UDGs.

     
    more » « less
  5. ABSTRACT

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

     
    more » « less