Reducing theoretical uncertainties in Galactic dark matter (DM) searches is an important challenge as several experiments are now delving into the parameter space relevant to popular (particle or not) candidates. Since many DM signal predictions rely on the knowledge of the DM velocity distribution---direct searches, capture by stars, $p$-wave-suppressed or Sommerfeld-enhanced annihilation rate, microlensing of primordial black holes, etc.---it is necessary to assess the accuracy of our current theoretical handle. Beyond Maxwellian approximations or ad-hoc extrapolations of fits on cosmological simulations, approaches have been proposed to self-consistently derive the DM phase-space distribution only from the detailed mass content of the Galaxy and some symmetry assumptions (e.g. the Eddington inversion and its anisotropic extensions). Although theoretically sound, these methods are still based on simplifying assumptions and their relevance to real galaxies can be questioned. In this paper, we use zoomed-in cosmological simulations to quantify the associated uncertainties. Assuming isotropy, we predict the speed distribution and its moments from the DM and baryonic content measured in simulations, and compare them with the true ones. We reach a predictivity down to ~ 10% for some observables, significantly better than some Maxwellian models. This moderate theoretical error is particularly encouraging at a time when stellar surveys like the Gaia mission should allow us to improve constraints on Galactic mass models.

Searches for "dark" subhaloes in gamma-ray point-like source catalogues are among promising strategies for indirect dark matter detection. Such a search is nevertheless affected by uncertainties related, on the one hand, to the modelling of the dark matter subhalo distribution in Milky-Way-like galaxies, and, on the other hand, to the sensitivity of gamma-ray instruments to the dark matter subhalo signals. In the present work, we assess the detectability of dark matter subhaloes in Fermi-LAT catalogues, taking into accounts uncertainties associated with the modelling of the Galactic subhalo population. We use four different halo models bracketing a large set of uncertainties. For each model, adopting an accurate detection threshold of the LAT to dark matter subhalo signals and comparing model predictions with the number of unassociated point-sources in Fermi-LAT catalogues, we derive upper limits on the annihilation cross section as a function of dark matter mass. Our results show that, even in the best-case scenario (i.e. DMonly subhalo model), which does not include tidal disruption from baryons, the limits on the dark matter parameter space are less stringent than current gamma-ray limits from dwarf spheroidal galaxies. Comparing the results obtained with the different subhalo models, we find that baryonic effects on the subhalo population are significant and lead to dark matter constraints that are less stringent by a factor of about 2 to 5. This uncertainty comes from the unknown resilience of dark matter subhaloes to tidal disruption.

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

The distribution of dark matter (DM) subhalos in our Galaxy remains disputed, leading to varying $\gamma$-ray and $\nu$ flux predictions from their annihilation or decay. In this work, we study how, in the inner Galaxy, subhalo tidal disruption from the Galactic baryonic potential impacts these signals. Based on state-of-the art modeling of this effect from numerical simulations and semi-analytical results, updated subhalo spatial distributions are derived and included in the CLUMPY code. The latter is used to produce a thousand realizations of the $\gamma$-ray and $\nu$ sky. Compared to predictions based on DM only, we conclude a decrease of the flux of the brightest subhalo by a factor 2 to 7 for annihilating DM and no impact on decaying DM: the discovery prospects or limits subhalos can set on DM candidates are affected by the same factor. This study also provides probability density functions for the distance, mass, and angular distribution of the brightest subhalo, among which the mass may hint at its nature: it is most likely a dwarf spheroidal galaxy in the case of strong tidal effects from the baryonic potential, whereas it is lighter and possibly a dark halo for DM only or less pronounced tidal effects.

We recently proposed a method to constrain $s$-wave annihilating MeV dark matter from a combination of the Voyager 1 and the AMS-02 data on cosmic-ray electrons and positrons. Voyager 1 actually provides an unprecedented probe of dark matter annihilation to cosmic rays down to $\sim 10$ MeV in an energy range where the signal is mostly immune to uncertainties in cosmic-ray propagation. In this article, we derive for the first time new constraints on $p$-wave annihilation down to the MeV mass range using cosmic-ray data. To proceed, we derive a self-consistent velocity distribution for the dark matter across the Milky Way by means of the Eddington inversion technique and its extension to anisotropic systems. As inputs, we consider state-of-the-art Galactic mass models including baryons and constrained on recent kinematic data. They allow for both a cored or a cuspy halo. We then calculate the flux of cosmic-ray electrons and positrons induced by $p$-wave annihilating dark matter and obtain very stringent limits in the MeV mass range, robustly excluding cross sections greater than $\sim 10^{-22}{\rm cm^3/s}$ (including theoretical uncertainties), about 5 orders of magnitude better than current CMB constraints. This limit assumes that dark matter annihilation is the sole source of cosmic rays and could therefore be made even more stringent when reliable models of astrophysical backgrounds are included.

Irrespective of the dark matter (DM) candidate, several potentially observable signatures derive from the velocity distribution of DM in halos, in particular in the Milky Way (MW) halo. Examples include direct searches for weakly-interacting massive particles (WIMPs), $p$-wave suppressed or Sommerfeld-enhanced annihilation signals, microlensing events of primordial black holes (PBHs), {\em etc}. Most current predictions are based on the Maxwellian approximation which is not only theoretically inconsistent in bounded systems, but also not supported by cosmological simulations. A more consistent method sometimes used in calculations for direct WIMP searches relies on the so-called Eddington inversion method, which relates the DM phase-space distribution function (DF) to its mass density profile and the total gravitational potential of the system. Originally built upon the isotropy assumption, this method can be extended to anisotropic systems. We investigate these inversion methods in the context of Galactic DM searches, motivated by the fact that the MW is a strongly constrained system, and should be even more so with the ongoing Gaia survey. We still draw conclusions that apply to the general case. In particular, we illustrate how neglecting the radial boundary of the DM halo leads to theoretical inconsistencies. We also show that several realistic configurations of the DM halo and the MW baryonic content entail ill-defined DFs, significantly restricting the configuration space over which these inversion methods can apply. We propose consistent solutions to these issues. Finally, we compute several observables inferred from constrained Galactic mass models relevant to DM searches (WIMPs or PBHs), {\em e.g.} moments and inverse moments of the DM speed and relative speed distributions.

The interaction properties of cold dark matter (CDM) particle candidates, such as those of weakly interacting massive particles (WIMPs), generically lead to the structuring of dark matter on scales much smaller than typical galaxies, potentially down to $\sim 10^{-10}M_\odot$. This clustering translates into a very large population of subhalos in galaxies and affects the predictions for direct and indirect dark matter searches (gamma rays and antimatter cosmic rays). In this paper, we elaborate on previous analytic works to model the Galactic subhalo population, while consistently with current observational dynamical constraints on the Milky Way. In particular, we propose a self-consistent method to account for tidal effects induced by both dark matter and baryons. Our model does not strongly rely on cosmological simulations as they can hardly be fully matched to the real Milky Way, but for setting the initial subhalo mass fraction. Still, it allows to recover the main qualitative features of simulated systems. It can further be easily adapted to any change in the dynamical constraints, and be used to make predictions or derive constraints on dark matter candidates from indirect or direct searches. We compute the annihilation boost factor, including the subhalo-halo cross-product. We confirm that tidal effects induced by the baryonic components of the Galaxy play a very important role, resulting in a local average subhalo mass density $\lesssim 1\%$ of the total local dark matter mass density, while selecting in the most concentrated objects and leading to interesting features in the overall annihilation profile in the case of a sharp subhalo mass function. Values of global annihilation boost factors range from $\sim 2$ to $\sim 20$, while the local annihilation rate is about half as much boosted.