# Remnants of Galactic Subhalos and Their Impact on Indirect Dark-Matter Searches

Martin Stref, Thomas Lacroix, Julien Lavalle

► To cite this version:

Martin Stref, Thomas Lacroix, Julien Lavalle. Remnants of Galactic Subhalos and Their Impact on Indirect Dark-Matter Searches. Galaxies, 2019, 7 (2), pp.65. 10.3390/galaxies7020065. hal-02123884

HAL Id: hal-02123884

https://hal.science/hal-02123884

Submitted on 5 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Remnants of Galactic Subhalos and Their Impact on Indirect Dark-Matter Searches

Martin Stref <fnref n="1" />,* , Thomas Lacroix <fnref n="2" /> and Julien Lavalle <fnref n="1" />

Received: 6 May 2019; Accepted: 29 May 2019; Published: 4 June 2019

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** Dark-matter subhalos, predicted in large numbers in the cold-dark-matter scenario, should have an impact on dark-matter-particle searches. Recent results show that tidal disruption of these objects in computer simulations is overefficient due to numerical artifacts and resolution effects. Accounting for these results, we re-estimated the subhalo abundance in the Milky Way using semianalytical techniques. In particular, we showed that the boost factor for gamma rays and cosmic-ray antiprotons is increased by roughly a factor of two.

**Keywords:** particle dark matter; subhalos; indirect searches

## 1. Introduction

There is overwhelming evidence that most of the matter in the universe is nonbaryonic [1]. An exciting possibility to account for these puzzling observations is that the universe is filled with exotic particles that interact only very weakly with ordinary matter [2,3]. One of the most elegant and popular dark-matter (DM) particle candidates is the Weakly Interacting Massive Particle (WIMP). These hypothetical particles are being looked for in particle colliders [4–6], in direct detection experiments [7–9], and in cosmic radiation [10,11], so far without success. Although one of the motivations for WIMPs is related to the fact that they emerge naturally in particle theories addressing a hierarchy problem [12–14], WIMPs are also attractive stemming from their very simple thermal-production mechanism in the early universe. Moreover, a large fraction of available parameter space is still unconstrained and currently actively explored [15]. For completeness, it is worth recalling that many alternatives to WIMPs exist that we do not discuss here, like axions [16], sterile neutrinos [17], primordial black holes [18], and extended dark sectors [19]. The cosmological paradigm best supported by current probes is that DM is cold, i.e., collisionless and nonrelativistic. This implies a structuring of matter on scales smaller than typical galaxies, with a model-dependent cutoff [20,21]. Interestingly, subgalactic scales are those where there could be departures from the predictions of the cold DM paradigm because of some observational issues [22]. This might sign new specific properties of the dark matter (e.g., [23]), or it could be due to baryonic effects (e.g., [24]). This motivates a detailed inspection of the impact of DM properties on the smallest scales, irrespective of the underlying scenario.

The small-scale structuring of DM, as treated, for instance, in the WIMP scenario, translates into a large population of subhalos within galactic halos [25–27]. Modeling these subhalos is crucial if one is to make accurate predictions for direct and indirect DM searches. This is a difficult task, as numerical simulations are far from resolving the smallest structures predicted by the cold DM paradigm. To incorporate the smallest structures, one can extrapolate the results of simulations over orders of magnitude in scales (see, e.g., [28]) but this represents a leap of faith. On the other hand,

<footnote n="1">¹ Laboratoire Univers et Particules de Montpellier (LUPM), Université de Montpellier & CNRS, Place Eugène Bataillon, 34095 Montpellier CEDEX 05, France; lavalle@in2p3.fr</footnote>

<footnote n="2">² Instituto de Fisica Teórica, C/ Nicolás Cabrera 13-15, Campus de Cantoblanco UAM, 28049 Madrid, Spain; thomas.lacroix@uam.es</footnote>

<footnote>* Correspondence: martin.stref@umontpellier.fr</footnote>
one can employ semianalytical models (see, e.g., [29–32]). The difficulty with the latter is accounting for the tidal effects experienced by subhalos within the host galaxy. These models can be calibrated on cosmological simulations, which are supposed to consistently describe the tidal stripping of subhalos in their host halo. However, it was recently pointed out by van den Bosch and collaborators [33,34] that simulations are plagued with numerical artifacts that lead to a significant overestimate of the tidal stripping efficiency, and therefore to an underestimate of the actual subhalo population even within the numerical resolution limit. An alternative and complementary way to study the tidal stripping of subhalos is to rely on analytical or semianalytical methods, which are based on first principles and allow to deal with subhalo mass scales, down to the free-streaming scale. Here, we review the semianalytical model developed by Stref and Lavalle [35] (SL17 hereafter), which incorporates a realistic and kinematically constrained Milky Way mass model (including baryons) and predicts the galactic subhalo abundance.<fnref n="1" /> This model accounts for different sources of tidal effects, and can easily accommodate to different prescriptions for tidal disruption efficiency.

This paper is structured as follows. In Section 2, we briefly review the SL17 model and discuss the resilience of subhalos to tidal effects in light of recent analyses of simulation results [33,34]. In Section 3, we compute the DM mass density within subhalos, as well as the number density of these objects in the Milky Way. Finally, in Section 4, we look at the impact of our results on indirect searches for annihilating DM, focusing on gamma rays and cosmic-ray antiprotons.

## 2. Semianalytical Model of Galactic Subhalos

In this section, we review the SL17 Galactic subhalo population model and discuss the tidal effects experienced by subhalos. We then propose a way of incorporating the recent results of van den Bosch and collaborators in the model in a consistent calibration procedure.

### 2.1. Review of the Stref and Lavalle Model

SL17 is a semianalytical model of galactic subhalos that is built upon dynamical constraints and cosmological considerations. The main input of the model is the initial subhalo phase-space density

$$ \frac{dN}{dV dm dc}(\vec{r}, m, c) \propto \frac{dP_v}{dV}(\vec{r}) \times \frac{dP_m}{dm}(m) \times \frac{dP_c}{dc}(c, m), \quad (1) $$

where phase space refers to the position–mass–concentration space. Functions $dP_v/dV$, $dP_m/dm$ and $dP_c/dc$ are the spatial, mass, and concentration distributions, respectively. It is assumed that, should subhalos behave as hard spheres (as is the case for single DM “particles” in a cosmological simulation), they would be spatially distributed as $dP_v/dV \propto \rho_{DM}$ where $\rho_{DM}$ is the total DM density profile of the galaxy. This sets our initial conditions before tidal disruption. The smooth DM mass density is computed through

$$ \rho_{\text{sm}}(\vec{r}) = \rho_{\text{DM}}(\vec{r}) - \langle \rho_{\text{cl}} \rangle (\vec{r}), \quad (2) $$

where $\langle \rho_{cl} \rangle$ is the average DM mass density inside clumps (this quantity is explicitly computed in Section 3). In the following, we use the galactic mass models constrained by McMillan [36] on pre-Gaia data for the DM and baryonic mass distributions. In this framework, Equation (2) ensures the compatibility of our subhalo model with the constrained DM profile $\rho_{DM}$. Although this work is devoted to the study of Milky Way subhalos, SL17 can in principle be used to study the substructure population in any virialized DM system. One only needs a mass model for the system in question and a proper calibration of the subhalo mass fraction through the procedure outlined in Section 2.3.

<footnote n="1">¹ We refer to this model as semianalytical because it involves integrals that must be computed numerically. The model does not rely on numerical simulations except at the level of a calibration described in Section 2.3.</footnote>
Mass *m* and concentration *c* refer to cosmological mass *m*<sub>200</sub> and concentration *c*<sub>200</sub> (defined with respect to the critical density), where we dropped the 200 index for convenience. The subhalo mass function measured in simulations is consistent with a power law [26,27]

$$ \frac{d\mathcal{P}_m}{dm}(m) \propto m^{-\alpha_m} \Theta(m - m_{\text{min}}) \Theta(m_{\text{max}} - m), \qquad (3) $$

where Θ is the Heaviside step function, and the power-law index is $\alpha_m = 1.9$ or $\alpha_m = 2$. These values of $\alpha_m$ encompass the Press and Schechter [37] mass function and the Sheth and Tormen [38] mass function, as illustrated in Figure 1. These functions can be computed directly from the matter power spectrum in the framework of excursion set theory [39], for spherical collapse (Press–Schechter) and ellipsoidal collapse (Sheth–Tormen). Thus, the two power-law indices we considered bracket the theoretical uncertainties on the small-scale mass function. If the DM is made of WIMPs, mass cutoff $m_{\text{min}}$ can be related to the kinetic decoupling of the DM particle and is found to lie between $10^{-4} M_{\odot}$ and $10^{-10} M_{\odot}$ [29,40–45]. Maximal mass $m_{\text{max}}$ is set to $0.01 \times M_{\text{DM}}$, where $M_{\text{DM}}$ is the total DM mass in the Milky Way. Concentration distribution $d\mathcal{P}_c/dc$ is generically found to exhibit a log-normal distribution for field halos [46,47], which defines our initial concentration distribution (before tidal stripping). We adopt the peak value and variance fit in Sánchez-Conde and Prada [48], which was shown to provide a good description of cosmological simulations run independently by several groups. Subhalos are assumed to have a Navarro–Frenk–White (NFW) profile [49] with parameters set by *m* and *c* (the impact of choosing an Einasto profile [50] instead of NFW was investigated in [35]).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Mass function $dn/dm$ multiplied by $m^2$. We show the prediction of Press and Schechter [37] (red line), Sheth and Tormen [38] (green line) as well as the power-law mass functions with index $\alpha_m = 1.9$ (magenta line) and $\alpha_m = 2$ (blue line). The Press–Schechter and Sheth–Tormen mass functions were computed for the cosmology of Planck 2018 [1] using the transfer function of Eisenstein and Hu [51] and a sharp-k filter. All mass functions are normalized to unity with $M_{\text{min}} = 10^{-10} M_{\odot}$.</figcaption></figure>

The subhalo population is strongly affected by tidal interactions with the potential of the host galaxy [52]. This is accounted for in the model through the calculation of a tidal radius $r_t$ for each subhalo. The tidal radius should be interpreted as the physical extension of a subhalo, which is in general smaller than the extension it would have on a flat background. The physical mass of a subhalo is then

$$ m_t(\vec{r}, m, c) = 4\pi \int_{0}^{r_t(\vec{r}, m, c)} dx \, x^2 \rho_{\text{sub}}(x) \le m, \qquad (4) $$
where $\rho_{sub}$ is the subhalo inner mass density profile. In our modeling, tidal stripping only removes the outer layers of subhalos while leaving the inner parts unchanged. In reality, DM should rearrange itself into a new equilibrium state. Central density however, should be left essentially unchanged [53]. Since we are interested in indirect searches for annihilating DM, and the central density gives the dominant contribution to the annihilation rate, our modeling should lead to a reasonable approximation. Two important contributions are accounted for in SL17: the effect of the smooth galactic potential (including both DM and baryons), and the gravitational shocking induced by the baryonic disk [54,55]. The latter effect turns out to be very efficient at stripping subhalos in the inner 20 kpc of the galaxy, a result also found in numerical studies [56–58]. The strength of SL17 over simulations is that it accounts for the constrained potential of the MW, with a detailed description of baryons.

## 2.2. Subhalo Disruption?

Whether a subhalo can be completely disrupted by tidal effects is an open question. A number of numerical studies found that a subhalo is completely disrupted when the total energy gained through tidal-stripping or disk-shocking effects is comparable to the binding energy [56,59]. On the other hand, some studies [60–62] found that cuspy subhalos almost always survive mass loss, leaving a small bound remnant behind even after gaining an energy far greater than their binding energy. These contradictory results may have been reconciled in a recent series of papers by van den Bosch and collaborators [33,34,63]. In these studies, it is shown that subhalo disruption in N-body simulations can actually be entirely explained by numerical artifacts. In particular, disruption is shown to be highly sensitive to the value of the force-softening length. If this length is taken sufficiently small, the authors showed that subhalos survive tidal mass loss in the form of a small bound remnant. We aim at quantifying the impact of these results on the whole subhalo population. Tidal disruption is modeled in a very simple way in SL17: given a subhalo with scale radius $r_s$ and tidal radius $r_t$, we assume

$$ \frac{r_t(\vec{r}, m, c)}{r_s(m, c)} < \epsilon_t \Leftrightarrow \text{subhalo is disrupted} \qquad (5) $$

In Equation (5), $\epsilon_t$ is a dimensionless free parameter assumed universal, i.e., independent of the subhalo's mass, concentration, or position. In SL17, the value of the disruption parameter was set to $\epsilon_t = 1$ in agreement with numerical results (see, e.g., [59]). The results of van den Bosch and collaborators point toward a much lower value for $\epsilon_t$. In this work, we consider two extreme values: $\epsilon_t = 1$ and $\epsilon_t = 0.01$. The latter means a subhalo is disrupted when it has lost around 99.99% of its mass. In the following, we refer to these two configurations as “fragile subhalos” ($\epsilon_t = 1$) and “resilient subhalos” ($\epsilon_t = 0.01$). The final subhalo phase-space density can now be written:

$$ \frac{dN}{dV dm dc}(\vec{r}, m, c) = \frac{N_{\text{tot}}}{K_{\text{tot}}} \frac{dP_v}{dV}(\vec{r}) \times \frac{dP_m}{dm}(m) \times \frac{dP_c}{dc}(c, m) \Theta\left(\frac{r_t(\vec{r}, m, c)}{r_s(m, c)} - \epsilon_t\right), \qquad (6) $$

where $N_{tot}$ is the total number of substructures within the virial radius of the Milky Way, and $K_{tot}$ is a normalization factor:

$$ K_{\text{tot}} = \int dV \frac{dP_v}{dV}(\vec{r}) \int dm \frac{dP_m}{dm}(m) \int dc \frac{dP_c}{dc}(c, m) \Theta\left(\frac{r_t(\vec{r}, m, c)}{r_s(m, c)} - \epsilon_t\right). \qquad (7) $$

## 2.3. Calibration Procedure

In its current version, the SL17 model requires a calibration of the subhalo abundance in a given mass range (this will change in future versions). To be consistent with results from the highest-resolution simulations available, calibration is done by demanding that the subhalo mass fraction is similar to what is found in the dark matter-only Via Lactea II simulation [26]. This amounts to 11% of the total dark halo mass in the form of subhalos in the virial mass range $[m_1, m_2] = [2.2 \times 10^{-6} M_{\text{DM}}, 8.8 \times 10^{-4} M_{\text{DM}}]$ where $M_{\text{DM}}$ is the total DM mass of the galaxy. We stress
that these numbers are expressed in terms of virial masses, not tidal masses (see [35] for further details). To reproduce the (likely overestimated) tidal disruption efficiency in simulations, we have to set $\epsilon_t = 1$ at the calibration stage, and we also neglect the impact of baryons. Disruption efficiency parameter $\epsilon_t$ can safely be changed after the calibration has been completed. It is much safer to perform this calibration on dark-matter-only simulations because tidal stripping induced by baryons strongly depends on the details of the stellar distribution, which is acutely constrained in the Milky Way. More formally, the normalization procedure reads:

$$f_{\text{sub}}(m_1, m_2) = \frac{1}{M_{\text{DM}}} \int dV \int_{m_1}^{m_2} dm \int dc \times m \times \left. \frac{dN}{dV dm dc} \right|_{\text{DMO}, \epsilon_t=1} . \qquad (8)$$

Fixing $f_{\text{sub}}(m_1, m_2) = 0.11$ leads to the total number of clumps $N_{\text{DMO}, \epsilon_t=1}$ in the simulation-like configuration. Note that this value assumes that $m$ is really $m_{200}$ in the equation above, not the tidal mass.

Now that the model is properly calibrated, we incorporate all the effects that are not included in the calibration, i.e., the tidal effects due to the baryons and possibly $\epsilon_t < 1$. This is done by assuming that subhalos in the outskirts of the galaxy are not affected by baryonic tides or the value of $\epsilon_t$. This is motivated by the observation that tidal effects are inefficient far from the center of the galaxy, and subhalos almost behave like isolated halos. The DM mass within clumps per unit of volume can be expressed as

$$\langle \rho_{\text{cl}} \rangle (\vec{r}) = \int_{m_{\text{min}}}^{m_{\text{max}}} dm \int_{1}^{\infty} dc \frac{dN}{dV dm dc} m_t(\vec{r}, m, c), \qquad (9)$$

where $m_t$ is the tidal subhalo mass introduced in Equation (4). Equating $\langle \rho_{\text{cl}} \rangle (r_{200})$, where $r_{200}$ is the virial radius of the galaxy, in the DM-only + $\epsilon_t = 1$ configuration, to the same quantity in the realistic configuration (including baryons and $\epsilon_t \le 1$) leads to the simple relation

$$\frac{N_{\text{DMO}, \epsilon_t=1}}{K_{\text{DMO}, \epsilon_t=1}} = \frac{N_{\text{tot}}}{K_{\text{tot}}}. \qquad (10)$$

The two normalization factors $K$ can be computed using Equation (7), and we obtain the value of $N_{\text{tot}}$. The number of subhalos within the solar radius $r_{\odot} = 8.21$ kpc is shown in Table 1. This number is highly sensitive to the parameters of mass function $\alpha_m$ and $m_{\text{min}}$, as already shown in Reference [35]. Furthermore, it is quite sensitive to $\epsilon_t$: going from $\epsilon_t = 1$ to $\epsilon_t = 0.01$, the number of subhalos increases by at least an order of magnitude. The impact of $\epsilon_t$ on subhalo mass and number density is investigated in the next section.

<table><caption>**Table 1.** Top panel: number of subhalos within $r_{\odot} = 8.21$ kpc, for different values of mass function parameters $\alpha_m$ and $m_{\text{min}}$, for $\epsilon_t = 1$. Bottom panel: same as top panel, for $\epsilon_t = 0.01$.</caption><thead><tr><th>ε<sub>t</sub> = 1</th><th>m<sub>min</sub> = 10<sup>-4</sup> M<sub>☉</sub></th><th>m<sub>min</sub> = 10<sup>-10</sup> M<sub>☉</sub></th></tr></thead><tbody><tr><td>α<sub>m</sub> = 1.9</td><td>1.90 × 10<sup>10</sup></td><td>1.55 × 10<sup>16</sup></td></tr><tr><td>α<sub>m</sub> = 2</td><td>2.64 × 10<sup>11</sup></td><td>8.40 × 10<sup>17</sup></td></tr><tr><th>ε<sub>t</sub> = 0.01</th><th>m<sub>min</sub> = 10<sup>-4</sup> M<sub>☉</sub></th><th>m<sub>min</sub> = 10<sup>-10</sup> M<sub>☉</sub></th></tr><tr><td>α<sub>m</sub> = 1.9</td><td>6.64 × 10<sup>11</sup></td><td>1.68 × 10<sup>17</sup></td></tr><tr><td>α<sub>m</sub> = 2</td><td>9.06 × 10<sup>12</sup></td><td>9.10 × 10<sup>18</sup></td></tr></tbody></table>

### 3. Mass and Number Densities of Subhalos

In this section, we computed the mass density within subhalos, as well as the subhalo number density. Subhalo mass density is defined in Equation (9). Once subhalo density is known, it is used to
determine the amount of DM smoothly distributed across the galaxy through Equation (2). The DM mass inside subhalos was compared with the total DM density $\rho_{DM}$ in Figure 2. Mass density in the form of subhalos is predicted to be much higher, by orders of magnitude, for resilient subhalos than for fragile subhalos. The former case is also more theoretically justified, although the latter one allows us to compare with very conservative assumptions. At the position of the Solar System, the impact is around one order of magnitude. Although these are large differences, we note the subhalo mass density is still far below the total DM density. This means that most of the DM mass within the orbit of the Sun is smoothly distributed rather than clumpy, irrespective of the efficiency of the tidal disruption set by $\epsilon_t$.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Left panel: Dark matter (DM) mass density inside subhalos $\langle\rho_{cl}\rangle$ for $m_{min} = 10^{-4} M_{\odot}$. The mass function index is $\alpha_m = 2$ (red) or $\alpha_m = 1.9$ (blue). We show the result for $\epsilon_t = 1$ (dashed) and $\epsilon_t = 0.01$ (solid). Total DM density is shown as a black solid curve for comparison. Right panel: same as left panel, for $m_{min} = 10^{-10} M_{\odot}$.</figcaption></figure>

The subhalo number density in SL17 can be formally written:

$$ \frac{dN}{dV}(\vec{r}) = \int dm \int dc \frac{dN}{dV dm dc} \qquad (11) $$

$$ = \frac{N_{\text{tot}}}{K_{\text{tot}}} \frac{dP_v}{dV}(\vec{r}) \int_{m_{\text{min}}}^{m_{\text{max}}} dm \frac{dP_m}{dm}(m) \int_{1}^{\infty} dc \frac{dP_c}{dc}(c, m) \Theta\left(\frac{r_t(\vec{r}, m, c)}{r_s} - \epsilon_t\right). \qquad (12) $$

The obtained results are shown in Figure 3. Number density, just like mass density, is highly sensitive to $\alpha_m$ and $m_{min}$, as well as the disruption parameter. Interestingly, the values we get in the Solar neighborhood are comparable to the local number density of stars $n_* \sim 1 \text{ pc}^{-3}$. For a low value of minimal mass $m_{min}$, the subhalo number density can even be much higher, possibly going as high as $10^5 \text{ pc}^{-3}$. This could have a number of interesting implications for the interactions between subhalos and stars. The tidal heating of subhalos by stars has been investigated in a number of studies [29,60,64–68], with different conclusions. In the next section, we look at the impact of our results on indirect DM searches.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Left panel: subhalo number density for $m_{min} = 10^{-4}$ M$_{\odot}$. Mass function index is $\alpha_m = 2$ (red) or $\alpha_m = 1.9$ (blue). We show the result for $\epsilon_t = 1$ (dashed) and $\epsilon_t = 0.01$ (solid). Right panel: same as left panel, for $m_{min} = 10^{-10}$ M$_{\odot}$.</figcaption></figure>

## 4. Impact on Indirect Searches for Annihilating Dark Matter

In this section, we quantify the impact of galactic clumps on indirect searches for self-annihilating DM. Inhomogeneities are known to enhance the DM annihilation rate in galactic halos [69]. We computed the local DM self-annihilation rate and evaluated the enhancement due to the survival of clump remnants, referred to as the boost factor. Two complementary channels were then investigated: gamma rays and antiproton cosmic rays.

### 4.1. Annihilation Profiles and Local Boost Factors

The number of self-annihilation of DM particles at position $\vec{r}$ is proportional to $\rho^2(\vec{r})$, where $\rho$ is DM mass density. If subhalos are discarded, the galactic annihilation profile is

$$ \mathcal{L}_0(\vec{r}) = \rho_{\text{DM}}^2(\vec{r}). \qquad (13) $$

Let us now consistently include the contribution of subhalos. The luminosity of a single clump is

$$ L_t(\vec{r}, m, c) = \int_{V_t} d^3\vec{r} \rho_{\text{sub}}^2(\vec{r}), \qquad (14) $$

where $V_t(\vec{r}, m, c)$ is the volume of the clump within its tidal radius. The annihilation of the full subhalo population is simply obtained by integrating the luminosity of a single object over the subhalo phase-space number density:

$$ \mathcal{L}_{\text{cl}}(\vec{r}) = \int_{m_{\text{min}}}^{m_{\text{max}}} dm \int_{1}^{+\infty} dc \frac{dN}{dV dm dc} L_t(\vec{r}, m, c). \qquad (15) $$

The full annihilation profile must also incorporate the annihilation in the smooth halo (different from $\mathcal{L}_0$, which is the density assuming all the DM is smoothly distributed), as well as the annihilation of subhalo particles onto smooth halo particles. The first contribution can be written as:

$$ \mathcal{L}_{\text{sm}}(\vec{r}) = \rho_{\text{sm}}^2(\vec{r}), \qquad (16) $$
where $\rho_{\text{sm}}(\vec{r}) = \rho_{\text{DM}}(\vec{r}) - \langle \rho_{\text{cl}} \rangle(\vec{r})$ is the smooth DM density. The clump-smooth contribution is:

$$ \mathcal{L}_{\text{cs}}(\vec{r}) = 2 \rho_{\text{sm}}(\vec{r}) \langle \rho_{\text{cl}} \rangle(\vec{r}) \qquad (17) $$

$$ = 2 \rho_{\text{sm}}(\vec{r}) \int_{m_{\text{min}}}^{m_{\text{max}}} dm \int_{1}^{+\infty} dc \frac{dN}{dV dm dc} m_t(\vec{r}, m, c). \qquad (18) $$

The total annihilation profile is simply the sum of all contributions:

$$ \mathcal{L}(\vec{r}) = \mathcal{L}_{\text{cl}}(\vec{r}) + \mathcal{L}_{\text{sm}}(\vec{r}) + \mathcal{L}_{\text{cs}}(\vec{r}). \qquad (19) $$

This should be compared to $\mathcal{L}_0(\vec{r})$ to evaluate the impact of clustering on the annihilation rate. This is usually done in terms of a boost factor, which we define as

$$ 1 + \mathcal{B}(\vec{r}) = \frac{\mathcal{L}(\vec{r})}{\mathcal{L}_0(\vec{r})}. \qquad (20) $$

This is not quite the boost factor used in indirect searches, which is defined through a ratio of fluxes (see Equations (23) and 29)). The boost in Equation (20) is rather the local increase in the annihilation rate due to clustering. According to this definition, the boost is zero if $\mathcal{L} = \mathcal{L}_0$, i.e., substructures are not included.<fnref n="2" />

The annihilation profiles are shown on the top panels in Figure 4, and the associated boost factors are shown on the bottom panels. As already shown a long time ago, see e.g., References [29,70], the boost is an increasing function of the galactocentric radius $r = |\vec{r}|$. This is due to the morphology of the annihilation profiles that is modified by the inclusion of clumps: we have $\mathcal{L}_{\text{cl}} \propto \rho_{\text{DM}}$ while $\mathcal{L}_0 \propto \rho_{\text{DM}}^2$. The high sensitivity of the annihilation profile to the mass function index $\alpha_m$ is also noticeable, which is by far the largest source of uncertainty on clump contribution. The value of disruption parameter $\epsilon_t$ has almost no impact on the boost above 20 kpc due to the ineffectiveness of tidal effects far from the center of the galaxy. Below 20 kpc, the boost is strongly sensitive to the disruption parameter. In the inner few kiloparsecs, fixing $\epsilon_t = 0.01$ leads to a boost orders of magnitude larger than in the $\epsilon_t = 1$ configuration. We, however, have $\mathcal{B}(\vec{r}) \ll 1$ below 3 kpc regardless of the value of $\epsilon_t$, meaning $\mathcal{L} \simeq \mathcal{L}_0$ and subhalos do not have any impact on the annihilation rate in that region. The region where the impact of $\epsilon_t$ on the annihilation profile is the most important is located between 3 and 10 kpc. This region coincidentally includes the Solar System, located at $r \simeq 8$ kpc. This motivates a more detailed investigation of two standard annihilation channels: gamma rays and cosmic-ray antiprotons, which are sensitive to different annihilation regions.

<footnote n="2">² This differs by one unit from the definition used in Reference [35].</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Total luminosity density profiles as defined in Equation (19), for a mass function index $\alpha_m = 2$ (red) and $\alpha_m = 1.9$ (blue). We show the results for efficient tidal disruption ($\epsilon_t = 1$, dashed) and very resilient clumps ($\epsilon_t = 0.01$, solid). (Left panel) Results for $m_{min} = 10^{-4}$ M$_\odot$, (right panel) results for $m_{min} = 10^{-10}$ M$_\odot$. The total luminosity density without clumps is displayed as a solid black line on each panel.</figcaption></figure>

## 4.2. Application to Gamma Rays

The energy-differential flux of gamma rays originating from DM self-annihilation is, on a given line of sight,

$$ \frac{d\Phi_{\gamma}}{dE d\Omega} = \frac{1}{4\pi 2 m_{\chi}^2} \frac{\langle \sigma v \rangle}{dE} \int ds \rho^2, \qquad (21) $$

where $\langle \sigma v \rangle$ is the thermally averaged annihilation cross-section, $m_\chi$ is the DM mass, $dN_\gamma/dE$ is the gamma-ray spectrum at annihilation, and $s$ the distance coordinate along the line of sight.<fnref n="3" /> If the annihilation cross-section is velocity-independent, astrophysical ingredients only enter through the J-factor, defined as

$$ J(\psi) = \int ds \rho^2(s, \psi), \qquad (22) $$

where $\psi$ is the angle between the direction of the galactic center and the line of sight (spherical symmetry of the dark halo is assumed). The impact of small-scale clustering on this J-factor has been considered in a number of studies [32,48,70–78]. We define the gamma-ray boost factor as

$$ 1 + \mathcal{B}_{\gamma}(\psi) = \frac{J(\psi)}{J_0(\psi)}, \qquad (23) $$

where $J_0 = \int ds \mathcal{L}_0$ is the J-factor without subhalos. Unlike local boost $\mathcal{B}$ in Equation (20), gamma-ray boost $\mathcal{B}_\gamma$ depends on line of sight rather than the position in the galaxy [70]. The boost is shown as a function of $\psi$ in Figure 5. The growth of local boost $\mathcal{B}(\vec{r})$ as a function of $r$ translates into a growth of $\mathcal{B}_\gamma(\psi)$ as a function of $\psi$, i.e., the maximal gamma-ray boost is reached at the anticenter. This maximal boost ranges from 0.5 to 9, depending on the values of $\alpha_m$ and $m_{min}$. The survival of clumps noticeably increases the boost at all latitudes. The gain is greater at small latitudes where substructures are more

<footnote n="3">³ The expression of $d\Phi_\gamma/dE$ should be multiplied by 1/2 if the DM particle is not its own antiparticle.</footnote>
impacted by tidal effects. Below $\psi \simeq 40$ deg, the boost is increased by a factor of at least two in all configurations. This should have important consequences for indirect searches using gamma rays, especially at high latitudes. Interestingly, high latitudes have been shown to be a very sensitive probe of DM annihilation even without the inclusion of clumps [79].

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Left panel: gamma-ray boost factor as defined in Equation (23), as a function of angle $\psi$ between the direction of the galactic center and the line of sight, for a minimal subhalo mass of $m_{min} = 10^{-4} M_{\odot}$. We show the results for efficient tidal disruption ($\epsilon_t = 1$, dashed) and very resilient clumps ($\epsilon_t = 0.01$, solid). Right panel: same as left panel, for $m_{min} = 10^{-10} M_{\odot}$.</figcaption></figure>

### 4.3. Application to Cosmic-Ray Antiprotons

Charged cosmic rays constitute an indirect detection channel, complementary to gamma rays [10,80]. Since their original proposal as a probe of DM annihilation [81], cosmic-ray antiprotons have been shown to be especially sensitive, see e.g., References [82–88]. Antiprotons have been the subject of much scrutiny since the latest measurement of the antiproton flux performed by the AMS-2 collaboration [89]. A number of studies [90–95] have found a discrepancy between the measured flux and a purely secondary origin of antiprotons. This discrepancy could be interpreted as evidence for annihilating DM, although the significance of the excess is debated as it depends on the propagation model used and the modeling of systematic uncertainties. In this context, it is worth evaluating systematic uncertainties coming from small-scale structuring.

Since antiprotons have a random motion due to their diffusion on the inhomogeneities of the magnetic halo, their detection gives little information on their source. This implies that the antiproton boost factor, and the boost for charged cosmic rays in general, is not direction-dependent, unlike for gamma rays. Instead, this boost is energy-dependent [96] and has been shown to be mild, at most a factor of two [31,74]. Although smaller than the gamma-ray boost, this can still be larger than the systematic uncertainties on cosmic-ray propagation. This motivates a new computation of the boost, which we performed here. The antiproton boost factor is defined as:

$$1 + \mathcal{B}_{\bar{p}}(T) = \frac{\Phi_{\bar{p}}(T)}{\Phi_{\bar{p},0}(T)}, \qquad (24)$$
where T is the antiproton kinetic energy, $\Phi_{\bar{p}}$ is the DM-induced antiproton flux including the subhalo contribution, and $\Phi_{\bar{p},0}$ the same flux assuming all the DM is smoothly distributed. To obtain the flux, one must solve the cosmic-ray steady-state propagation equation

$$-K\Delta\Psi + \vec{\nabla} \cdot (\vec{V}_c\Psi) + \partial_E [b\Psi - K_{EE}\partial_E\Psi] + 2h\delta(z)\Gamma_{\text{ann}}\Psi = Q_{\text{DM}}, \quad (25)$$

which accounts for spatial diffusion, convection, energy losses, diffusive re-acceleration, and spallation processes in the disk (taken as infinitely thin). In Equation (25), $\Psi$ is the antiproton number density per unit energy that is related to the flux through $\Phi_{\bar{p}} = v_{\bar{p}}/(4\pi) \times \Psi$ where $v_{\bar{p}}$ is the antiproton speed. Antiprotons are sourced by DM annihilation:

$$Q_{\text{DM}}(E, \vec{r}) = \frac{\langle \sigma v \rangle}{2} \frac{dN_{\bar{p}}}{dE} \left( \frac{\rho(\vec{r})}{m_{\chi}} \right)^2, \quad (26)$$

where $dN_{\bar{p}}/dE$ is the antiproton spectrum at annihilation.<fnref n="4" /> Several unknown propagation parameters enter Equation (25). These can be constrained using the measured boron-to-carbon ratio (B/C) [97]. We used the best-fit model derived by Reinert and Winkler [92], which includes an energy break in the diffusion coefficient. The B/C ratio can only constrain $K_0/L$, where $K_0$ is the normalization of the diffusion coefficient, and $L$ is the half-height of the magnetic halo. As shown in Reference [83], the DM-induced antiproton flux crucially depends on $L$; hence, we considered two extremal values in this work. A lower bound on $L$ can be obtained from low-energy positron data [98], and the authors of Reference [92] found $L = 4.1$ kpc. For the largest value, we took $L = 15$ kpc. According to the analysis of Reference [92], the B/C data are consistent with negligible re-acceleration. Furthermore, we neglected energy losses that are unimportant for high-energy antiprotons. The resulting transport equation can be solved semianalytically using Green's function formalism (see Reference [31] for the solution), and the differential flux can be written as:

$$\frac{d\Phi_{\bar{p}}}{dT d\Omega} = \frac{v_{\bar{p}} \langle \sigma v \rangle}{4\pi 2 m_{\chi}^2} \int dE_s \int d^3\vec{r}_s G(E \leftarrow E_s; \vec{r}_{\odot} \leftarrow \vec{r}_s) \frac{dN_{\bar{p}}}{dE}(E_s) \rho^2(\vec{r}_s). \quad (27)$$

Since all energy-dependent terms have been neglected in Equation (25), the energy part of Green's function is trivial

$$G(E \leftarrow E_s; \vec{r}_{\odot} \leftarrow \vec{r}_s) = \delta(E - E_s) \times \bar{G}(\vec{r}_{\odot} \leftarrow \vec{r}_s), \quad (28)$$

and the boost factor can be simply written:

$$1 + B_{\bar{p}}(T) = \frac{\int d^3\vec{r}_s \bar{G}(\vec{r}_{\odot} \leftarrow \vec{r}_s) \mathcal{L}(\vec{r}_s)}{\int d^3\vec{r}_s \bar{G}(\vec{r}_{\odot} \leftarrow \vec{r}_s) \mathcal{L}_0(\vec{r}_s)}. \quad (29)$$

The boost factor is shown as a function of the antiproton kinetic energy in Figure 6. We first note that the boost is roughly energy-independent. This is because antiprotons probe the entire volume of the magnetic halo during their lifetime, independently of their energy. This is not true at low energies, below a few GeVs, where energy losses become relevant. The half-height of the magnetic halo has a small impact on the boost, with $L = 15$ kpc leading to a slightly larger $B_{\bar{p}}$ than $L = 4.1$ kpc. The main source of uncertainties are coming from subhalo parameters $a_m$, $m_{\text{min}}$, and $\epsilon_t$. Survival parameter $\epsilon_t$ had significant impact on the result, with a small value $\epsilon_t = 0.01$ leading to a boost roughly twice as large as in the $\epsilon_t = 1$ case. As for gamma rays, the most critical parameter is $a_m$. For a low-value $a_m = 1.9$, the boost never exceeded 10%, while it was always higher for $a_m = 2$.

<footnote n="4">⁴ If DM is not its own antiparticle, Q should be divided by 2.</footnote>
Overall, playing with the propagation and subhalo parameters, we found that the antiproton boost can conservatively range from 2% to 140%. These values are in agreement with earlier results. Although it is conservative to ignore small-scale clustering when deriving limits on the annihilation cross-section using data, the boost should be included when interpreting an excess as a signature of DM annihilation. Indeed, a factor of two in the DM contribution would change the inferred mass and cross-section of the hypothetical DM particle.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Top panel: antiproton boost factor as a function of kinetic energy, as defined in Equation (29). We show the result for a half-height of the magnetic halo of $L = 15$ kpc, and a minimal subhalo mass of $m_{\min} = 10^{-4}$ M$_{\odot}$ (left) or $m_{\min} = 10^{-10}$ M$_{\odot}$ (right). Bottom panel: same as top panel, for half-height $L = 4.1$ kpc.</figcaption></figure>

**5. Conclusions**

Subhalos suffer mass loss due to their interaction with the tidal field of the galaxy, which makes their modeling very challenging. Consequently, most subhalo models rely at least partly on numerical simulations to calibrate their predictions. However, it was recently shown that numerical simulations might not properly account for the tidal disruption of subhalos, as artificial effects lead to a serious overestimation of the efficiency of these processes [33,34]. We note that the resistance of subhalos to tidal stripping is further supported by theoretical arguments, like adiabatic invariance that should
prevail in their inner parts [55], as already emphasized in Reference [35]. We derived some of the consequences of these results using the semianalytical galactic subhalo population model of Stref and Lavalle [35], assuming tidal disruption efficiency $\epsilon_t = 0.01$. We predicted the spatial dependence of the subhalo properties due to tides induced both by the global gravitational potential and baryonic disk shocking. We remind the reader that this model is built from constrained mass models for the Milky Way and is therefore consistent with current kinematic constraints, which is usually not the case in extrapolations from “Milky Way-like” simulations. We found that the local mass density is still dominated by the smooth component of the dark halo. The local number density of subhalos is increased by roughly one order of magnitude with respect to estimates based on a tidal disruption efficiency similar to that inferred from simulations ($\epsilon_t = 1$). This makes the subhalo number density comparable, for a broad range of minimal subhalo mass $m_{min}$, to the local star number density. Since our description of the subhalo population relies solely on gravitational principles, it should be valid for a wide range of cold DM candidate such as WIMPs, axions, or primordial black holes. One only needs to modify the mass function, in particular, $m_{min}$, to explore alternatives to the WIMP scenario. The resilience of subhalos increases the local WIMP annihilation rate, which, in turn, affects predictions for indirect searches. For gamma rays, we found that the boost factor was increased by at least a factor of two for $\psi < 40$ deg, and slightly less for higher values of $\psi$. The boost factor for antiprotons is also increased by a rough factor of two if subhalos are resilient to tidal disruption. For a complementary study comparing the SL17 model to simulation results regarding indirect searches in gamma rays and neutrinos, we refer the reader to Reference [99].

For future work, we plan on including a more detailed mass function, directly deriving from the primordial power spectrum, as well as the tidal heating of subhalos due to individual stars, and studying the consequences of having a large population of small objects in the Solar neighborhood.

**Author Contributions:** conceptualization, all authors; software, M.S.; investigation, M.S.; writing—original-draft preparation, M.S.; writing—review and editing, T.L. and J.L.; supervision, J.L.

**Acknowledgments:** We wish to thank M. A. Sánchez-Conde and M. Doro for inviting us to contribute to this topical review.

**Funding:** J.L. and M.S. are partly supported by the Agence Nationale pour la Recherche (ANR) Project No. ANR-18-CE31-0006, the Origines, Constituants, et EVolution de l'Univers (OCEVU) Labex (No. ANR-11-LABX-0060), the CNRS IN2P3-Theory/INSU-PNHE-PNCG project "Galactic Dark Matter," and the European Union's Horizon 2020 Research and Innovation Program under Marie Skłodowska-Curie Grant Agreements No. 690575 and No. 674896, in addition to recurrent funding by the Centre National de la Recherche Scientifique (CNRS) and the University of Montpellier. T.L. is supported by the European Union's Horizon 2020 Research and Innovation Program under the Marie Skłodowska-Curie grant agreement No. 713366. The work of TL was also supported by the Spanish Agencia Estatal de Investigación through grants PGC2018-095161-B-100, IFT Centro de Excelencia Severo Ochoa SEV-2016-0597, and Red Consolider MultiDark FPA2017-90566-REDC.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

## References

1. Aghanim, N.; Akrami, Y.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Ballardini, M.; Banday, A.J.; Barreiro, R.B.; Bartolo, N.; Basak, S.; et al. Planck 2018 results. VI. Cosmological parameters. arXiv **2018**, arXiv:1807.06209.

2. Bertone, G.; Hooper, D.; Silk, J. Particle dark matter: Evidence, candidates and constraints. Phys. Rep. **2005**, *405*, 279-390. [CrossRef]

3. Feng, J.L. Dark Matter Candidates from Particle Physics and Methods of Detection. Ann. Rev. Astron. Astrophys. **2010**, *48*, 495-545. [CrossRef]

4. Fairbairn, M.; Kraan, A.C.; Milstead, D.A.; Sjöstrand, T.; Skands, P.; Sloan, T. Stable massive particles at colliders. Phys. Rep. **2007**, *438*, 1-63. [CrossRef]

5. Goodman, J.; Ibe, M.; Rajaraman, A.; Shepherd, W.; Tait, T.M.P.; Yu, H.B. Constraints on Dark Matter from Colliders. Phys. Rev. **2010**, D82, 116010. [CrossRef]
6. Kahlhoefer, F. Review of LHC Dark Matter Searches. *Int. J. Mod. Phys. **2017***, A32, 1730006. [CrossRef]

7. Cerdeno, D.G.; Green, A.M. Direct detection of WIMPs. In *Particle Dark Matter: Observations, Models and Searches*; Cambridge University Press: Cambridge, UK, 2010; pp. 347–369.

8. Baudis, L. Direct dark matter detection: The next decade. *Phys. Dark Univ. **2012***, 1, 94–108. [CrossRef]

9. Freese, K.; Lisanti, M.; Savage, C. Colloquium: Annual modulation of dark matter. *Rev. Mod. Phys. **2013***, 85, 1561–1581. [CrossRef]

10. Lavalle, J.; Salati, P. Dark matter indirect signatures. *C. R. Phys. **2012***, 13, 740–782. [CrossRef]

11. Gaskins, J.M. A review of indirect searches for particle dark matter. *Contemp. Phys. **2016***, 57, 496–525. [CrossRef]

12. Jungman, G.; Kamionkowski, M.; Griest, K. Supersymmetric dark matter. *Phys. Rep. **1996***, 267, 195–373. [CrossRef]

13. Servant, G.; Tait, T.M.P. Is the lightest Kaluza-Klein particle a viable dark matter candidate? *Nuclear Phys. B **2003***, 650, 391–419. [CrossRef]

14. Agashe, K.; Servant, G. Warped Unification, Proton Stability, and Dark Matter. *Phys. Rev. Lett. **2004***, 93, 231805. [CrossRef]

15. Leane, R.K.; Slatyer, T.R.; Beacom, J.F.; Ng, K.C.Y. GeV-scale thermal WIMPs: Not even slightly ruled out. *Phys. Rev. **2018***, D98, 023016. [CrossRef]

16. Marsh, D.J.E. Axion Cosmology. *Phys. Rep. **2016***, 643, 1–79. [CrossRef]

17. Boyarsky, A.; Ruchayskiy, O.; Shaposhnikov, M. The Role of sterile neutrinos in cosmology and astrophysics. *Ann. Rev. Nucl. Part. Sci. **2009***, 59, 191–214. [CrossRef]

18. Carr, B.; Kuhnel, F.; Sandstad, M. Primordial Black Holes as Dark Matter. *Phys. Rev. **2016***, D94, 083504. [CrossRef]

19. Alexander, J.; Battaglieri, M.; Echenard, B.; Essig, R.; Graham, M.; Izaguirre, G.; Jaros, J.; Krnjaic, G.; Mardon, J.; Morrissey, D.; et al. Dark Sectors 2016 Workshop: Community Report. *arXiv **2016***, arXiv:1608.08632.

20. Blumenthal, G.R.; Faber, S.M.; Primack, J.R.; Rees, M.J. Formation of galaxies and large-scale structure with cold dark matter. *Nature **1984***, 311, 517–525. [CrossRef]

21. Davis, M.; Efstathiou, G.; Frenk, C.S.; White, S.D.M. The evolution of large-scale structure in a universe dominated by cold dark matter. *Astrophys. J. **1985***, 292, 371–394. [CrossRef]

22. Bullock, J.S.; Boylan-Kolchin, M. Small-Scale Challenges to the ΛCDM Paradigm. *Ann. Rev. Astron. Astrophys. **2017***, 55, 343–387. [CrossRef]

23. Tulin, S.; Yu, H.B. Dark Matter Self-interactions and Small Scale Structure. *Phys. Rep. **2018***, 730, 1–57. [CrossRef]

24. Pontzen, A.; Governato, F. How supernova feedback turns dark matter cusps into cores. *Mon. Not. R. Astron. Soc. **2012***, 421, 3464. [CrossRef]

25. Moore, B.; Ghigna, S.; Governato, F.; Lake, G.; Quinn, T.R.; Stadel, J.; Tozzi, P. Dark matter substructure within galactic halos. *Astrophys. J. **1999***, 524, L19–L22. [CrossRef]

26. Diemand, J.; Kuhlen, M.; Madau, P.; Zemp, M.; Moore, B.; Potter, D.; Stadel, J. Clumps and streams in the local dark matter distribution. *Nature **2008***, 454, 735–738. [CrossRef]

27. Springel, V.; Wang, J.; Vogelsberger, M.; Ludlow, A.; Jenkins, A.; Helmi, A.; Navarro, J.F.; Frenk, C.S.; White, S.D.M. The Aquarius Project: The subhalos of galactic halos. *Mon. Not. R. Astron. Soc. **2008***, 391, 1685–1711. [CrossRef]

28. Springel, V.; White, S.D.M.; Frenk, C.S.; Navarro, J.F.; Jenkins, A.; Vogelsberger, M.; Wang, J.; Ludlow, A.; Helmi, A. Prospects for detecting supersymmetric dark matter in the Galactic halo. *Nature **2008***, 456, 73–76. [CrossRef]

29. Berezinsky, V.; Dokuchaev, V.; Eroshenko, Y. Small-scale clumps in the galactic halo and dark matter annihilation. *Phys. Rev. **2003***, D68, 103003. [CrossRef]

30. Penarrubia, J.; Benson, A.J. Effects of dynamical evolution on the distribution of substructures. *Mon. Not. R. Astron. Soc. **2005***, 364, 977–989. [CrossRef]

31. Lavalle, J.; Yuan, Q.; Maurin, D.; Bi, X.J. Full calculation of clumpiness boost factors for antimatter cosmic rays in the light of ΛCDM N-body simulation results. Abandoning hope in clumpiness enhancement? *Astron. Astroph. **2008***, 479, 427–452. [CrossRef]
32. Bartels, R.; Ando, S. Boosting the annihilation boost: Tidal effects on dark matter subhalos and consistent luminosity modeling. *Phys. Rev.* **2015**, *D92*, 123508. [CrossRef]

33. van den Bosch, F.C.; Ogiya, G.; Hahn, O.; Burkert, A. Disruption of Dark Matter Substructure: Fact or Fiction? *Mon. Not. R. Astron. Soc.* **2018**, *474*, 3043–3066. [CrossRef]

34. van den Bosch, F.C.; Ogiya, G. Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption. *Mon. Not. R. Astron. Soc.* **2018**, *475*, 4066–4087. [CrossRef]

35. Stref, M.; Lavalle, J. Modeling dark matter subhalos in a constrained galaxy: Global mass and boosted annihilation profiles. *Phys. Rev.* **2017**, *D95*, 063003. [CrossRef]

36. McMillan, P.J. The mass distribution and gravitational potential of the Milky Way. *Mon. Not. R. Astron. Soc.* **2017**, *465*, 76–94. [CrossRef]

37. Press, W.H.; Schechter, P. Formation of Galaxies and Clusters of Galaxies by Self-Similar Gravitational Condensation. *Astrophys. J.* **1974**, *187*, 425–438. [CrossRef]

38. Sheth, R.K.; Tormen, G. Large-scale bias and the peak background split. *Mon. Not. R. Astron. Soc.* **1999**, *308*, 119–126. [CrossRef]

39. Bond, J.R.; Cole, S.; Efstathiou, G.; Kaiser, N. Excursion set mass functions for hierarchical Gaussian fluctuations. *Astrophys. J.* **1991**, *379*, 440–460. [CrossRef]

40. Hofmann, S.; Schwarz, D.J.; Stoecker, H. Damping scales of neutralino cold dark matter. *Phys. Rev.* **2001**, *D64*, 083507. [CrossRef]

41. Loeb, A.; Zaldarriaga, M. The Small-scale power spectrum of cold dark matter. *Phys. Rev.* **2005**, *D71*, 103520. [CrossRef]

42. Green, A.M.; Hofmann, S.; Schwarz, D.J. The First wimpy halos. *J. Cosmol. Astropart. Phys.* **2005**, *0508*, 003. [CrossRef]

43. Bertschinger, E. The Effects of Cold Dark Matter Decoupling and Pair Annihilation on Cosmological Perturbations. *Phys. Rev.* **2006**, *D74*, 063509. [CrossRef]

44. Bringmann, T.; Hofmann, S. Thermal decoupling of WIMPs from first principles. *J. Cosmol. Astropart. Phys.* **2007**, *0704*, 016. [CrossRef]

45. Bringmann, T. Particle Models and the Small-Scale Structure of Dark Matter. *New J. Phys.* **2009**, *11*, 105027. [CrossRef]

46. Bullock, J.S.; Kolatt, T.S.; Sigad, Y.; Somerville, R.S.; Kravtsov, A.V.; Klypin, A.A.; Primack, J.R.; Dekel, A. Profiles of dark haloes: Evolution, scatter and environment. *Mon. Not. R. Astron. Soc.* **2001**, *321*, 559–575. [CrossRef]

47. Macciò, A.V.; Dutton, A.A.; van den Bosch, F.C. Concentration, spin and shape of dark matter haloes as a function of the cosmological model: WMAP1, WMAP3 and WMAP5 results. *Mon. Not. R. Astron. Soc.* **2008**, *391*, 1940–1954. [CrossRef]

48. Sánchez-Conde, M.A.; Prada, F. The flattening of the concentration mass relation towards low halo masses and its implications for the annihilation signal boost. *Mon. Not. R. Astron. Soc.* **2014**, *442*, 2271–2277. [CrossRef]

49. Navarro, J.F.; Frenk, C.S.; White, S.D.M. A Universal Density Profile from Hierarchical Clustering. *Astrophys. J.* **1997**, *490*, 493–508. [CrossRef]

50. Einasto, J. On the Construction of a Composite Model for the Galaxy and on the Determination of the System of Galactic Parameters. *Trudy Astrofizicheskogo Instituta Alma-Ata* **1965**, *5*, 87–100.

51. Eisenstein, D.J.; Hu, W. Baryonic features in the matter transfer function. *Astrophys. J.* **1998**, *496*, 605. [CrossRef]

52. Binney, J.; Tremaine, S. Galactic Dynamics; Princeton University Press: Princeton, NJ, USA, 1987.

53. Drakos, N.E.; Taylor, J.E.; Benson, A.J. The phase-space structure of tidally stripped haloes. *Mon. Not. R. Astron. Soc.* **2017**, *468*, 2345–2358. [CrossRef]

54. Ostriker, J.P.; Spitzer, L., Jr.; Chevalier, R.A. On the Evolution of Globular Clusters. *Astrophys. J. Lett.* **1972**, *176*, L51. [CrossRef]

55. Gnedin, O.Y.; Ostriker, J.P. On the self-consistent response of stellar systems to gravitational shocks. *Astrophys. J.* **1999**, *513*, 626. [CrossRef]

56. D'Onghia, E.; Springel, V.; Hernquist, L.; Keres, D. Substructure Depletion in the Milky Way Halo by the Disk. *Astrophys. J.* **2010**, *709*, 1138–1147. [CrossRef]
57. Errani, R.; Peñarrubia, J.; Laporte, C.F.P.; Gómez, F.A. The effect of a disc on the population of cuspy and cored dark matter substructures in Milky Way-like galaxies. *Mon. Not. R. Astron. Soc.* **2017**, *465*, L59–L63. [CrossRef]

58. Kelley, T.; Bullock, J.S.; Garrison-Kimmel, S.; Boylan-Kolchin, M.; Pawlowski, M.S.; Graus, A.S. Phat ELVIS: The inevitable effect of the Milky Way’s disk on its dark matter subhaloes. *arXiv* **2018**, arXiv:1811.12413.

59. Hayashi, E.; Navarro, J.F.; Taylor, J.E.; Stadel, J.; Quinn, T.R. The Structural evolution of substructure. *Astrophys. J.* **2003**, *584*, 541–558. [CrossRef]

60. Goerdt, T.; Gnedin, O.Y.; Moore, B.; Diemand, J.; Stadel, J. The survival and disruption of CDM micro-haloes: Implications for direct and indirect detection experiments. *Mon. Not. R. Astron. Soc.* **2007**, *375*, 191–198. [CrossRef]

61. Berezinsky, V.; Dokuchaev, V.; Eroshenko, Y. Remnants of dark matter clumps. *Phys. Rev.* **2008**, *D77*, 083519. [CrossRef]

62. Peñarrubia, J.; Benson, A.J.; Walker, M.G.; Gilmore, G.; McConnachie, A.W.; Mayer, L. The impact of dark matter cusps and cores on the satellite galaxy population around spiral galaxies. *Mon. Not. R. Astron. Soc.* **2010**, *406*, 1290–1305. [CrossRef]

63. van den Bosch, F.C. Dissecting the Evolution of Dark Matter Subhaloes in the Bolshoi Simulation. *Mon. Not. R. Astron. Soc.* **2017**, *468*, 885–909. [CrossRef]

64. Berezinsky, V.; Dokuchaev, V.; Eroshenko, Y. Destruction of small-scale dark matter clumps in the hierarchical structures and galaxies. *Phys. Rev.* **2006**, *D73*, 063504. [CrossRef]

65. Zhao, H.; Taylor, J.E.; Silk, J.; Hooper, D. Tidal Disruption of the First Dark Microhalos. *Astrophys. J.* **2007**, *654*, 697–701. [CrossRef]

66. Angus, G.W.; Zhao, H. Analysis of galactic tides and stars on CDM microhalos. *Mon. Not. R. Astron. Soc.* **2007**, *375*, 1146–1156. [CrossRef]

67. Green, A.M.; Goodwin, S.P. On mini-halo encounters with stars. *Mon. Not. R. Astron. Soc.* **2007**, *375*, 1111–1120. [CrossRef]

68. Schneider, A.; Krauss, L.; Moore, B. Impact of Dark Matter Microhalos on Signatures for Direct and Indirect Detection. *Phys. Rev.* **2010**, *D82*, 063525. [CrossRef]

69. Silk, J.; Stebbins, A. Clumpy cold dark matter. *Astrophys. J.* **1993**, *411*, 439–449. [CrossRef]

70. Bergstrom, L.; Edsjo, J.; Gondolo, P.; Ullio, P. Clumpy neutralino dark matter. *Phys. Rev.* **1999**, *D59*, 043506. [CrossRef]

71. Diemand, J.; Kuhlen, M.; Madau, P. Dark matter substructure and gamma-ray annihilation in the Milky Way halo. *Astrophys. J.* **2007**, *657*, 262–270. [CrossRef]

72. Strigari, L.E.; Koushiappas, S.M.; Bullock, J.S.; Kaplinghat, M. Precise constraints on the dark matter content of Milky Way dwarf galaxies for gamma-ray experiments. *Phys. Rev.* **2007**, *D75*, 083526. [CrossRef]

73. Pieri, L.; Bertone, G.; Branchini, E. Dark Matter Annihilation in Substructures Revised. *Mon. Not. R. Astron. Soc.* **2008**, *384*, 1627. [CrossRef]

74. Pieri, L.; Lavalle, J.; Bertone, G.; Branchini, E. Implications of high-resolution simulations on indirect dark matter searches. *Phys. Rev. D* **2011**, *83*, 023518. [CrossRef]

75. Blanchet, S.; Lavalle, J. Diffuse gamma-ray constraints on dark matter revisited. I: The impact of subhalos. *J. Cosmol. Astropart. Phys.* **2012**, *1211*, 021. [CrossRef]

76. Ishiyama, T. Hierarchical Formation of Dark Matter Halos and the Free Streaming Scale. *Astrophys. J.* **2014**, *788*, 27. [CrossRef]

77. Moliné, A.; Sánchez-Conde, M.A.; Palomares-Ruiz, S.; Prada, F. Characterization of subhalo structural properties and implications for dark matter annihilation signals. *Mon. Not. R. Astron. Soc.* **2017**, *466*, 4974–4990. [CrossRef]

78. Hiroshima, N.; Ando, S.; Ishiyama, T. Modeling evolution of dark matter substructure and annihilation boost. *Phys. Rev.* **2018**, *D97*, 123002. [CrossRef]

79. Chang, L.J.; Lisanti, M.; Mishra-Sharma, S. Search for dark matter annihilation in the Milky Way halo. *Phys. Rev.* **2018**, *D98*, 123004. [CrossRef]

80. Maurin, D.; Taillet, R.; Donato, F.; Salati, P.; Barrau, A.; Boudoul, G. Galactic cosmic ray nuclei as a tool for astroparticle physics. *arXiv* **2002**, arXiv:astro-ph/0212111.

81. Silk, J.; Srednicki, M. Cosmic-ray antiprotons as a probe of a photino-dominated universe. *Phys. Rev. Lett.* **1984**, *53*, 624–627. [CrossRef]
82. Bergström, L.; Edsjö, J.; Ullio, P. Cosmic Antiprotons as a Probe for Supersymmetric Dark Matter? *Astrophys. J.* **1999**, *526*, 215–235. [CrossRef]

83. Donato, F.; Fornengo, N.; Maurin, D.; Salati, P. Antiprotons in cosmic rays from neutralino annihilation. *Phys. Rev.* **2004**, *D69*, 063501. [CrossRef]

84. Barrau, A.; Salati, P.; Servant, G.; Donato, F.; Grain, J.; Maurin, D.; Taillet, R. Kaluza-Klein dark matter and Galactic antiprotons. *Phys. Rev.* **2005**, *D72*, 063507. [CrossRef]

85. Bringmann, T.; Salati, P. The galactic antiproton spectrum at high energies: Background expectation vs. exotic contributions. *Phys. Rev.* **2007**, *D75*, 083006. [CrossRef]

86. Cirelli, M.; Kadastik, M.; Raidal, M.; Strumia, A. Model-independent implications of the e+, anti-proton cosmic ray spectra on properties of Dark Matter. *Nucl. Phys.* **2009**, *B813*, 1–21. [CrossRef]

87. Boudaud, M.; Cirelli, M.; Giesen, G.; Salati, P. A fussy revisitation of antiprotons as a tool for Dark Matter searches. *J. Cosmol. Astropart. Phys.* **2015**, *5*, 13. [CrossRef]

88. Giesen, G.; Boudaud, M.; Génolini, Y.; Poulin, V.; Cirelli, M.; Salati, P.; Serpico, P.D. AMS-02 antiprotons, at last! Secondary astrophysical component and immediate implications for Dark Matter. *J. Cosmol. Astropart. Phys.* **2015**, *1509*, 023. [CrossRef]

89. Aguilar, M.; Ali Cavasonza, L.; Alpat, B.; Ambrosi, G.; Arruda, L.; Attig, N.; Aupetit, S.; Azzarello, P.; Bachlechner, A.; Barao, F.; et al. Antiproton Flux, Antiproton-to-Proton Flux Ratio, and Properties of Elementary Particle Fluxes in Primary Cosmic Rays Measured with the Alpha Magnetic Spectrometer on the International Space Station. *Phys. Rev. Lett.* **2016**, *117*, 091103. [CrossRef]

90. Cuoco, A.; Krämer, M.; Korsmeier, M. Novel Dark Matter Constraints from Antiprotons in Light of AMS-02. *Phys. Rev. Lett.* **2017**, *118*, 191102. [CrossRef]

91. Cui, M.Y.; Yuan, Q.; Tsai, Y.L.S.; Fan, Y.Z. Possible dark matter annihilation signal in the AMS-02 antiproton data. *Phys. Rev. Lett.* **2017**, *118*, 191101. [CrossRef]

92. Reinert, A.; Winkler, M.W. A Precision Search for WIMPs with Charged Cosmic Rays. *J. Cosmol. Astropart. Phys.* **2018**, *1801*, 055. [CrossRef]

93. Cui, M.Y.; Pan, X.; Yuan, Q.; Fan, Y.Z.; Zong, H.S. Revisit of cosmic ray antiprotons from dark matter annihilation with updated constraints on the background model from AMS-02 and collider data. *J. Cosmol. Astropart. Phys.* **2018**, *1806*, 024. [CrossRef]

94. Cuoco, A.; Heisig, J.; Klamt, L.; Korsmeier, M.; Krämer, M. Scrutinizing the evidence for dark matter in cosmic-ray antiprotons. *arXiv* **2019**, arXiv:1903.01472.

95. Cholis, I.; Linden, T.; Hooper, D. A Robust Excess in the Cosmic-Ray Antiproton Spectrum: Implications for Annihilating Dark Matter. *arXiv* **2019**, arXiv:1903.02549.

96. Lavalle, J.; Pochon, J.; Salati, P.; Taillet, R. Clumpiness of dark matter and positron annihilation signal: Computing the odds of the galactic lottery. *Astron. Astrophys.* **2007**, *462*, 827–848. [CrossRef]

97. Strong, A.W.; Moskalenko, I.V. Propagation of cosmic-ray nucleons in the galaxy. *Astrophys. J.* **1998**, *509*, 212–228. [CrossRef]

98. Lavalle, J.; Maurin, D.; Putze, A. Direct constraints on diffusion models from cosmic-ray positron data: Excluding the minimal model for dark matter searches. *Phys. Rev.* **2014**, *D90*, 081301. [CrossRef]

99. Hütten, M.; Stref, M.; Combet, C.; Lavalle, J.; Maurin, D. γ-ray and ν searches for dark matter subhalos in the Milky Way with a baryonic potential. *arXiv* **2019**, arXiv:1904.10935.

<figure><img src="image_9.png" /></figure>