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 more »

Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
Award ID(s):
Publication Date:
NSF-PAR ID:
10368238
Journal Name:
The Astrophysical Journal
Volume:
932
Issue:
2
Page Range or eLocation-ID:
Article No. 128
ISSN:
0004-637X
Publisher:
DOI PREFIX: 10.3847
National Science Foundation
##### More Like this
1. Abstract

Signatures of vertical disequilibrium have been observed across the Milky Way’s (MW’s) disk. These signatures manifest locally as unmixed phase spirals inzvzspace (“snails-in-phase”), and globally as nonzero meanzandvz, wrapping around the disk into physical spirals in thexyplane (“snails-in-space”). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle MW-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using azvzasymmetry metric, we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal actionJϕ, excited by a single event (in our case, a satellite disk crossing); (b) these distinct phase spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase spirals offer a global view of the MW’s perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)–like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars withinR≲ 10 kpc. Hence such perturbations have widespread influence, which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW–Sgr interaction. Whilemore »

2. Abstract

We combine Gaia early data release 3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we use a Monte Carlo sample for the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare the dSph radial orbital phase to subhalosmore »

3. ABSTRACT

The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($M_{*}\lesssim 10^7\, \rm {M}_{\odot }$) are mostly quiescent and higher mass satellites ($M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$) are mostly star forming, with intermediate-mass satellites ($M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group ormore »

4. (Ed.)
ABSTRACT We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the Feedback In Realistic Environments-2 simulations. We select the 14 most massive satellites by stellar mass within $d_\mathrm{host}\le 300\mathrm{\, kpc}$ of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1–2 per cent of snapshots during z = 0−0.2, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for <500 Myr. However, if we select hosts with a Large Magellanic Cloud-like satellite near first pericentre, the fraction of snapshots with MW-like planes increases dramatically to 7–16 per cent, with lifetimes of  0.7–1 Gyr, likely because of group accretion of satellites. We find that M31’s satellite distribution is much more common: M31’s satellites lie within ∼1σ of the simulation median for every plane metric we consider. We find no significant difference in average satellite planarity for isolated hosts versus hosts in LG-like pairs. Baryonic and dark matter-only simulations exhibitmore »
5. ABSTRACT Reticulum II (Ret II) is a satellite galaxy of the Milky Way (MW) and presents a prime target to investigate the nature of dark matter (DM) because of its high mass-to-light ratio. We evaluate a dedicated INTEGRAL observation campaign data set to obtain γ-ray fluxes from Ret II and compare those with expectations from DM. Ret II is not detected in the γ-ray band 25–8000 keV, and we derive a flux limit of ${\lesssim}10^{-8}\, \mathrm{erg\, cm^{-2}\, s^{-1}}$. The previously reported 511 keV line is not seen, and we find a flux limit of ${\lesssim}1.7 \times 10^{-4}\, \mathrm{ph\, cm^{-2}\, s^{-1}}$. We construct spectral models for primordial black hole (PBH) evaporation and annihilation/decay of particle DM, and subsequent annihilation of e+s produced in these processes. We exclude that the totality of DM in Ret II is made of a monochromatic distribution of PBHs of masses ${\lesssim}8 \times 10^{15}\, \mathrm{g}$. Our limits on the velocity-averaged DM annihilation cross section into e+e− are $\langle \sigma v \rangle \lesssim 5 \times 10^{-28} \left(m_{\rm DM} / \mathrm{MeV} \right)^{2.5}\, \mathrm{cm^3\, s^{-1}}$. We conclude that analysing isolated targets in the MeV γ-ray band can set strong bounds on DM properties without multi-year data sets of the entire MW, and encourage follow-up observations ofmore »