## Impact of substructure on local dark matter searches

Alejandro Ibarra 1,2 , Bradley J. Kavanagh , and Andreas Rappelt 3 1

1 Physik-Department T30d, Technische Universität München, James-Franck-Straße, 85748 Garching, Germany

- 2
- 3

School of Physics, Korea Institute for Advanced Study, Seoul 02455, South Korea GRAPPA, Institute of Physics, University of Amsterdam, 1098 XH Amsterdam, The Netherlands

## Abstract

Dark matter substructure can contribute significantly to local dark matter searches and may provide a large uncertainty in the interpretation of those experiments. For direct detection experiments, sub-halos give rise to an additional dark matter component on top of the smooth dark matter distribution of the host halo. In the case of dark matter capture in the Sun, sub-halo encounters temporarily increase the number of captured particles. Even if the encounter happened in the past, the number of dark matter particles captured by the Sun can still be enhanced today compared to expectations from the host halo as those enhancements decay over time. Using results from an analytical model of the sub-halo population of a Milky Way-like galaxy, valid for sub-halo masses between 10 -5 M glyph[circledot] and 10 11 M glyph[circledot] , we assess the impact of sub-halos on direct dark matter searches in a probabilistic way. We find that the impact on direct detection can be sizable, with a probability of ∼ 10 -3 to find an O (1) enhancement of the recoil rate. In the case of the capture rate in the Sun, we find that O (1) enhancements are very unlikely, with probability glyph[lessorsimilar] 10 -5 , and are even impossible for some dark matter masses.

## 1 Introduction

One of the central assumptions of modern cosmology is the presence of a non-luminous matter component, commonly known as dark matter (DM) [1,2]. Although the exact nature of dark matter still remains unknown, many proposed models predict a small interaction with Standard Model particles, as in the case of weakly interacting massive particles (WIMPs) [3] appearing (for example) in supersymmetric theories [4]. This opens up the possibility to search for DM via various possible signatures. Currently pursued strategies include, among others, collider searches for dark matter production [5], indirect searches for the products of dark matter annihilation [6] as well as direct searches for the signatures of dark matter scattering off nuclei and electrons in the lab [7]. Furthermore, dark matter plays an important role on cosmological scales, driving the formation of structure in the Universe [8].

In the framework of cold and collisionless dark matter, the growth of structure is hierarchical, with small structures forming first and then later merging into larger structures [9, 10]. As a

consequence of this, we expect that the dark matter halos enclosing galaxies should contain a large number of sub-halos. Sub-halos of the Milky Way include the known dwarf satellite galaxies, as well as the Small and Large Magellanic Clouds, ranging in mass from roughly 10 7 M glyph[circledot] up to 10 10 M glyph[circledot] [11, 12]. However, these massive objects represent only a small fraction of the sub-halos present in the Milky Way. Many more sub-halos are too small to accrete enough material to form stars and they thus remain dark [13-16]. These dark sub-halos are expected to exist with masses all the way down to 10 -6 M glyph[circledot] , roughly the mass of the Earth [17-19]. In light of this, a number of studies have explored the properties and distribution of DM sub-halos within Milky Way-like galaxies, see e.g. Refs. [20-31] for a non-exhaustive list. Indeed, a recent study of stellar kinematics using data from the Gaia satellite [32] suggests that a substantial fraction of the local DM halo may be in substructure [33]. On ultra-local scales, the DM density could be enhanced by DM streams, such as Sagittarius [34] or S1 [35]. There has also been recent interest in the possibility that DM forms clusters [36] or 'blobs' [37], local overdensities perhaps generated by DM self-interactions.

Dark matter substructure can lead to a number of detectable signatures. One possibility is to directly look for the signal of DM annihilation inside sub-halos. Various studies [38-43] have investigated how this enhances the gamma-ray signal from the host halo and found boost factors up to O (10) . Another possibility is to look for evidence of close passages between DM sub-halos and cold stellar streams of the Milky Way [44-48]. These close passages would be imprinted as gaps and other features in the stellar streams, mapped out by precision astrometric surveys such as Gaia [32]. Furthermore, sub-structures might influence local DM searches by increasing the number of dark matter particles in a finite region, relative to the smooth background.

In this work, we focus on such 'local' DM searches: Direct Detection, the search for DMnucleus scattering in Earth-based detectors; and Solar Capture, the search for neutrino signals from DM which has scattered and been captured in the Sun [49]. The signal enhancement due to such clustering of dark matter has been studied for direct detection experiments in several works [50-54]. For example, it has been argued that no individual sub-halo should dominate the scattering rate [51] and that the contribution of streams to the local distribution should add in such a way as to be indistinguishable from a smooth halo [53]. However, the impact of substructure on direct detection may still be sizeable; for example Refs. [50,54] find a percent-level probability of the Earth currently passing through a sub-halo. The impact of sub-structure on Solar capture has been less extensively studied; for example, Ref. [55] studied a small number of benchmark scenarios for possible DM overdensities. For capture in the Sun, we note that sub-halo encounters during the whole Solar lifetime can be important. Each encounter leads to an enhancement in the number of captured particles, which may persist even after the Sun has traversed the sub-halo, depending on the efficiency of the annihilation process. It is even possible for these enhancements to accumulate if they persist long enough, highlighting the need for a careful treatment with realistic sub-halo populations.

Here, we revisit the impact of DM sub-structure on local searches. We describe the potential signal enhancements which could be observed due to individual DM sub-halos in both direct detection and Solar capture. Furthermore, we use state-of-the-art semi-analytic models for Milky Way substructure to derive the probability distributions for these enhancements for different dark matter masses and interaction types. In this way, we are able to assess how the population of sub-halos as a whole may impact local DM searches.

The structure of this article is as follows. In section 2, we summarize the expected properties of Milky Way sub-halos and derive the velocity distribution of sub-halos and dark matter

particles bound to them relative to the Sun or the Earth. We then update the calculation of the impact of sub-halos on direct detection experiments in section 3. In section 4, we build a probabilistic model to determine the boost factor of the neutrino signal due to the annihilation of dark matter particles which have been captured by scattering in the Sun. Finally, we conclude and summarize the results of this work in section 5.

## 2 Dark matter flux in the Solar System

The Milky Way disk is believed to be embedded in a halo of dark matter particles with density and velocity distributions which are, to a first approximation, stationary and spherically symmetric. For dark matter searches in the Solar System, the so-called Standard Halo Model (SHM) is commonly adopted, where the local density is ρ loc SHM = 0 3 GeV cm . -3 and the local velocity distribution takes the Maxwell-Boltzmann form [56]:

$$f _ { \text{SHM} } ( \vec { v } ) = \frac { 1 } { ( 2 \pi \sigma _ { v } ^ { 2 } ) ^ { 3 / 2 } N _ { \text{esc} } } \exp \left [ - \frac { | \vec { v } + \vec { v } _ { \odot } | ^ { 2 } } { 2 \sigma _ { v } ^ { 2 } } \right ] \quad \text{for } | \vec { v } | \leq v _ { \max } \,,$$

with normalization constant:

$$N _ { \text{esc} } = \text{erf} \left ( \frac { v _ { \text{esc} } } { \sqrt { 2 } \sigma _ { v } } \right ) - \sqrt { \frac { 2 } { \pi } } \frac { v _ { \text{esc} } } { \sigma _ { v } } \exp \left ( - \frac { v _ { \text{esc} } ^ { 2 } } { 2 \sigma _ { v } ^ { 2 } } \right ) \,.$$

Here, v glyph[circledot] ≈ 244 km/s is the local velocity of the Sun with respect to the Galactic frame [57-59], σ v ≈ 156 km/s is the velocity dispersion [56, 60] and v max = v esc + v glyph[circledot] is the maximal velocity of DM particles gravitationally bound to the Galaxy (expressed in the Solar frame), with v esc ≈ 544 km s / being the escape velocity from the Milky Way [61,62]. We note that while these values are observationally determined, they have associated uncertainties at the level of 10% [56,63,64]. The value of the local density ρ loc SHM is also known only to within a factor of a few [65] and a number of refinements to the SHM have been proposed (see e.g. [66-68]). In this work, we keep the properties of the smooth DM component fixed and focus not on deviations from the SHM, but on the possible additional contribution of DM substructure.

The Milky Way DM halo is expected to contain a population of sub-halos with different masses and located at different distances from the Galactic Center, which may alter, in a timedependent way, the local density and velocity distribution. The population of sub-halos can be characterized by the halo mass function (i.e. the number of sub-halos in a given mass range), their spatial distribution and their velocity distribution. A number of works have attempted to extract sub-halo populations from simulations (e.g. [22, 25, 27]) or model the processes which govern the evolution and properties of sub-halos (e.g. [29,31,69-71]).

In our work we will adopt the halo mass function reported in [31], obtained using a semianalytic prescription for sub-halo accretion and tidal stripping by the host halo. This prescription was calibrated against the results of N -body simulations [72], allowing the method to cover sub-halo masses from ∼ 10 11 M glyph[circledot] down to ∼ 10 -6 M glyph[circledot] , much smaller than those that can be resolved in the corresponding N -body simulations. We use tabulated halo mass functions and sub-halo properties provided by the authors of [31], though we note that fitting functions are also provided in Appendix A of Ref. [70]. The halo mass function at redshift z = 0 for a DM host halo of total mass 1 8 . × 10 12 M glyph[circledot] is shown for reference in the left panel of figure 1. This mass function may be written in the form d N/ M d ∝ M -α M ( ) , with α M ( ) a mass-dependent

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-ebd3c7c8-a1ed-4bba-91bc-a7129a5f599e.jpg)

/circledot

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-30f1a970-e9e7-4855-87c8-a765eaf1de0f.jpg)

/circledot

Figure 1: Left plot: Sub-halo mass function (ShMF) for a DM halo of total mass 1 8 . × 10 12 M glyph[circledot] , taken from Ref. [31], as well as for two parametric models with ShMF proportional to M -α with α = 1 8 . and 2 0 . . Right plot: Dependence of the slope α of the ShMF with the sub-halo mass for the cases showed in the left plot.

parameter which ranges between 1.8 and 2.0 over the relevant mass range, as shown in the right panel. We also show in figure 1 the halo mass function and α M ( ) for two parametric models with d N/ M d ∝ M -1 8 . and d N/ M d ∝ M -2 0 . , where the normalization of the power-law is fit to the results of [31]. The low mass sub-halos predicted in all of these models may have a significant impact in the search for dark matter inside the Solar System, as we will analyze in sections 3 and 4.

We will assume, following the results of [73], that sub-halos in the Milky Way are spatially distributed following an Einasto profile [74,75]

$$\ln ( n _ { \text{sh} } ( r ) / n _ { - 2 } ) = - \frac { 2 } { \gamma } \left [ \left ( \frac { r } { r _ { - 2 } } \right ) ^ { \gamma } - 1 \right ] \,,$$

with r -2 the radius at which the slope of the distribution equals -2 , n -2 the number density of sub-halos at that location, and γ a shape parameter. In our analysis we adopt the values γ = 0 854 . and r -2 = 245 1 . kpc from [73]. We choose n -2 such that ∫ V MW d r 3 n sh ( r ) = N sh , where N sh is the total number of sub-halos, which we calculate from the mass function of [31] to be approximately 2 6 . × 10 16 . It should be borne in mind that the Milky Way disk depletes the number of sub-halos close to the galactic center [76, 77]. While the mass function of [31] does not account for this depletion, the spatial distribution we adopt from [73] is derived from hydrodynamical simulations including baryons, so we do not further correct for this effect. Finally, we will also assume that the sub-halo centers of mass move in the Milky Way halo following the same Maxwell-Boltzmann distribution, Eq. (1), as the smooth DM component.

Each sub-halo has an internal structure which is characterized by a density distribution and a velocity distribution. For the density distribution we adopt an NFW density profile [78,79] up to a truncation radius r t , with density parameter ρ s and scale radius r s :

$$\rho _ { \text{sh} } ( r ) = \begin{cases} \frac { \rho _ { s } } { \left ( r / r _ { s } \right ) \left ( 1 + r / r _ { s } \right ) ^ { 2 } } & \text{for $r\leq r_{t}$} \,, \\ 0 & \text{for $r>r_{t}$} \,. \end{cases}$$

In field halos (those not bound within larger halos), the NFW profile would be truncated at the virial radius [80-83], defined such that the mean density of the halo inside the virial radius is a factor of 200 larger than the critical density (see e.g. [25]). For sub-halos, tidal stripping leads to a smaller truncation radius r t , which can be determined from the mass loss of the sub-halo within the host [84]. Similarly, the characteristic density and radius ρ s and r s at the present epoch may be determined from the corresponding values at accretion by accounting for tidal mass loss. The mass of an individual sub-halo is then given by

$$M = \frac { 4 \pi } { 3 } \rho _ { s } r _ { s } ^ { 3 } f ( c _ { V } ) \,,$$

where c V is the concentration parameter, defined as c V = r /r t s , and f ( c V ) = ln(1 + c V ) -c V / (1 + c V ) . For field halos, the concentration parameter is expected to follow a log-normal distribution [27, 85-87], leading to typical values in the range c V glyph[similarequal] 8 and c V glyph[similarequal] 35 . Instead, for MW sub-halos with masses between 10 -6 M glyph[circledot] and 10 11 M glyph[circledot] , the concentration parameter typically takes values between c V glyph[similarequal] 1 and c V glyph[similarequal] 35 , suggesting a larger number of denser, more concentrated sub-halos than those in the field. As for the mass function, we will make use of the distribution of individual sub-halo properties ( ρ s , r s , r t ) as calculated in Ref. [31].

To determine the internal velocity distribution of the individual particles in the sub-halo, we assume that it is virialized. Hence, the DM velocity distribution relative to the center of mass of the sub-halo has a Maxwell-Boltzmann form with dispersion σ sh :

$$f _ { \text{sh} | \text{sh} } ( \vec { v } ) = \frac { 1 } { ( 2 \pi \sigma _ { \text{sh} } ^ { 2 } ) ^ { 3 / 2 } } \exp \left [ - \frac { | \vec { v } | ^ { 2 } } { 2 \sigma _ { \text{sh} } ^ { 2 } } \right ] \$$

Accordingly the mean-square speed of sub-halo particles is 〈| glyph[vector] v | 2 〉 = 3 σ 2 sh . To determine the velocity dispersion we use the virial theorem, which relates the potential energy U of the subhalo to its kinetic energy T :

$$T = - \frac { 1 } { 2 } \, U \,.$$

The total kinetic energy of a sub-halo of mass M reads T = 1 2 M v 〈 2 〉 , while its potential energy is given by its gravitational binding energy

$$U = - 4 \pi G _ { N } \int _ { 0 } ^ { r _ { t } } \frac { M _ { \text{enc} } ( r ) } { r } r ^ { 2 } \rho _ { \text{sh} } ( r ) \text{d} r \,,$$

with M enc ( r ) the enclosed mass of the sub-halo as a function of the distance from the sub-halo center. For ρ sh ( r ) following an NFW profile, as in Eq. (4), we obtain:

$$U = - 8 \pi ^ { 2 } G _ { N } \rho _ { s } ^ { 2 } r _ { s } ^ { 5 } \left [ 1 - 2 \frac { \ln ( c _ { V } + 1 ) } { ( c _ { V } + 1 ) } - \frac { 1 } { ( c _ { V } + 1 ) ^ { 2 } } \right ] \quad \quad \quad$$

$$= - \left ( \frac { \pi } { 6 } \right ) ^ { 1 / 3 } \frac { \bar { G _ { N } c _ { V } } } { f ( c _ { V } ) ^ { 2 } } \bar { \rho } ^ { 1 / 3 } M ^ { 5 / 3 } \left [ 1 - 2 \frac { \bar { \ln } ( c _ { V } + 1 )$$

where ¯ = 3 ρ f ( c V ) ρ /c s 3 V is the mean density inside the sub-halo. Finally, applying the virial theorem, we obtain:

$$\sigma _ { \text{sh} } ^ { 2 } = - \frac { U } { 3 M } = \frac { 1 } { 3 } \left ( \frac { \pi } { 6 } \right ) ^ { 1 / 3 } \frac { G _ { N } c _ { V } } { f ( c _ { V } ) ^ { 2 } } \bar { \rho } ^ { 1 / 3 } M ^ { 2 / 3 } \left [ 1 - 2 \frac { \ln ( c _ { V } + 1 ) } { ( c _ { V } + 1 ) } - \frac { 1 } { ( c _ { V } + 1 ) ^ { 2 } } \right ] \,.$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-3313e1f2-06be-4fcd-a65a-1a10c588ae93.jpg)

/circledot

Figure 2: Internal velocity dispersion as a function of the sub-halo mass for concentration parameters within 1 σ from the mean value.

which depends mostly on the total sub-halo mass, M , and has a weak dependence on the concentration parameter, c V . The velocity dispersion is shown in figure 2 and can be approximated by

$$\sigma _ { \text{sh} } \approx 1. 5 \, \text{km/s} \, \left ( \frac { M } { 1 0 ^ { 6 } \, M _ { \odot } } \right ) ^ { \frac { 1 } { 3 } }.$$

The internal density and velocity distributions calculated above are expressed in the subhalo center-of-mass frame. When expressed in the Solar frame, one finds a time-dependent density distribution as the Sun traverses the sub-halo at a distance r from the sub-halo center, ρ loc sh ( ) t ≡ ρ loc sh [ r t ( )] . The resulting velocity distribution of DM particles inside the sub-halo reads:

$$f _ { \text{sh}, \vec { v } _ { \text{sh} } } ( \vec { v } ) = \frac { 1 } { ( 2 \pi \sigma _ { \text{sh} } ^ { 2 } ) ^ { 3 / 2 } } \exp \left [ - \frac { | \vec { v } + \vec { v } _ { \text{rel} } | ^ { 2 } } { 2 \sigma _ { \text{sh} } ^ { 2 } } \right ] \,,$$

where the velocity of the sub-halo with respect to the Sun is glyph[vector] v rel = glyph[vector] v glyph[circledot] -glyph[vector] v sh and where the distribution of the individual sub-halo velocities in the Galactic rest frame is f SHM ( glyph[vector] v sh ) , given in Eq. (1). For small sub-halos, namely sub-halos with M glyph[lessorsimilar] 10 6 M glyph[circledot] , the velocity dispersion is much smaller than the speed of the sub-halo center of mass, | glyph[vector] v sh | glyph[greatermuch] σ sh . Therefore, the sub-halo behaves effectively as stream of DM particles all moving with the same velocity. In the Solar frame, the velocity distribution of DM particles from a sub-halo with velocity glyph[vector] v sh can thus be approximated by:

$$f _ { \text{sh,} \vec { v } _ { \text{sh} } } ( \vec { v } ) \simeq \delta ( \vec { v } + \vec { v } _ { \text{rel} } ) \,.$$

Finally, we can write down the flux of dark matter particles in the Solar frame, assuming that the Sun travels inside a sub-halo with mass M and concentration parameter c V which moves with respect to the Galactic frame with velocity glyph[vector] v sh . This flux contains two components: a flux of DM particles coming from the smooth Milky Way dark matter halo, characterized by time-independent mass density and velocity distributions ρ SHM and f SHM ( glyph[vector] v ) , and a flux of DM particles from sub-halos, characterized by a time-dependent mass density ρ loc sh ( ) t and a

time-independent velocity distribution f sh ,glyph[vector] v sh ( glyph[vector] v ) . This flux reads: 1

$$F ( \vec { v }, t ) = v \left [ \frac { \rho _ { \text{SHM} } ^ { \text{loc} } } { m _ { \text{DM} } } f _ { \text{SHM} } ( \vec { v } ) + \frac { \rho _ { \text{sh} } ^ { \text{loc} } ( t ) } { m _ { \text{DM} } } f _ { \text{sh}, \vec { v } _ { \text{sh} } } ( \vec { v } ) \right ].$$

## 3 Impact of sub-halos in direct detection experiments

We now apply our model for DM sub-structure to determine the expected rate at direct detection experiments. The rate of nuclear recoils induced by scatterings of DM particles traversing a detector at the Earth is given by [88,89]:

$$R = \sum _ { i } \int _ { 0 } ^ { \infty } \mathrm d E _ { R } \, \epsilon _ { i } ( E _ { R } ) \frac { \xi _ { i } } { m _ { A _ { i } } } \int _ { v \geq v ^ { ( E _ { R } ) } _ { \min, i } } \mathrm d ^ { 3 } v F ( \vec { v } + \vec { v } _ { \text{Earth} }, t _ { 0 } ) \, \frac { \mathrm d \sigma _ { i } } { \mathrm d E _ { R } } ( v, E _ { R } ) \,.$$

Here, F glyph[vector] ( v + glyph[vector] v Earth , t 0 ) is the flux of DM received at the Earth at the present time (as described in the previous section), with glyph[vector] v the DM velocity expressed in the Solar frame. In addition, d σ / i d E R is the differential cross section for the elastic scattering of DM off a nuclear isotope i with mass m A i and mass fraction ξ i in the detector, producing a nuclear recoil with energy E R . Furthermore, v min ,i ( E R ) = √ m E / A i R (2 µ 2 A i ) is the minimal velocity necessary for a DM particle to induce a recoil with energy E R and µ A i is the reduced mass of the DM-nucleus system. The efficiency glyph[epsilon1] i ( E R ) gives the probability to detect a nuclear recoil off the target nucleus i . Finally, the total number of expected recoil events at a direct detection experiment reads N = R · E , with E the exposure (i.e. mass multiplied by live-time) of the experiment.

We consider the case when the Sun is traversing at the current epoch t 0 a sub-halo of mass M and concentration parameter c V , which is moving with velocity glyph[vector] v sh with respect to the Galactic frame. The scattering rate can be readily calculated inserting Eq. (15) into Eq. (16). For low mass sub-halos, for which the velocity distribution in the Solar frame can be approximated by Eq. (14), one finds a increment in the scattering rate

$$\mathcal { I } _ { R } \equiv \frac { R } { R _ { \text{SHM} } } - 1 \simeq \frac { \rho _ { \text{sh} } ^ { \text{loc} } [ r ( t _ { 0 } ) ] } { \rho _ { \text{SHM} } } \frac { R _ { \delta ( \vec { v } + \vec { v } _ { \text{rel} } ) } ^ { \rho _ { \text{SHM} } ^ { \text{loc} } } } { R _ { \text{SHM} } } \,.$$

Here, glyph[vector] v rel = glyph[vector] v glyph[circledot] + glyph[vector] v Earth -glyph[vector] v sh is the velocity of the sub-halo with respect to the Earth and ρ loc sh [ r t ( 0 )] is the contribution of the sub-halo to the total dark matter density (with the Sun located a distance r t ( 0 ) from the center of the sub-halo).

The factor R SHM is the event rate for the smooth SHM described at the start of section 2, while R ρ loc SHM δ glyph[vector] ( v + glyph[vector] v rel ) is the scattering rate for a stream of dark matter particles with density ρ SHM and velocity distribution f glyph[vector] ( v ) = δ glyph[vector] ( v + glyph[vector] v rel ) . This rate is illustrated in figure 3 for XENON1T [90] (left panel) and for CRESST [91] (right panel), for three representative dark matter masses. 2 For large dark matter masses, the scattering rate for streams (relative to the SHM) is O (0 1 . -1) for generic values of the relative velocity, and can be much larger for small dark matter masses and large relative velocities. Accordingly, there can in principle be a large enhancement in the scattering rate at these two experiments if the sub-halo also gives an O (1) contribution to the local dark matter density, as shown in Fig. 4.

1 We add the superscript 'loc' to emphasize that these quantities are evaluated locally, at the Sun or Earth's position.

2 To calculate the detector responses we used the DDCalc package [92,93].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-bab7250e-ddd2-454f-8d21-99316403cab7.jpg)

|

|

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-b25623a7-fd8f-4a63-9db9-ff8c04b58365.jpg)

|

|

Figure 3: Scattering rate induced at XENON1T (left plot) and CRESST (right plot) for dark matter streams with velocity | glyph[vector] v rel | , relative to the scattering rate assuming the Standard Halo Model, for dark matter masses m DM = 5 , 100 and 10000 GeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-45245b65-dcc0-42df-8864-54fad5bc784c.jpg)

Figure 4: Increment in the scattering rate I R relative to the expectations from the Standard Halo Model in the parameter space spanned by ρ loc sh [ r t ( 0 )] /ρ SHM (the contribution from sub-halos to the local density, relative to the SHM value) and R ρ loc SHM δ glyph[vector] ( v + glyph[vector] v rel ) /R SHM (the capture rate of dark matter particles from the sub-halo, relative to the capture rate of dark matter particles in the smooth halo).

From the current body of knowledge about the sub-halo population in the Milky Way (summarized in section 2), one can estimate the probability of having an increment in the event rate I R , defined in Eq. (17). Assuming that we are currently inside a sub-halo, the probability that the sub-halo moves with velocity glyph[vector] v sh with respect to the Galactic frame is given straightforwardly by the Maxwell-Boltzmann distribution in Eq. (1). We must then determine the probability of finding ourselves at the current time inside a sub-halo, along with the distribution of density enhancements ρ loc sh [ r t ( 0 )] ≡ ˜ ρ .

The probability P ρ, M, c (˜ V ) that we find ourselves inside a sub-halo of mass M and concentration c V such that the total DM density receives an additional contribution ˜ ρ can be

written:

$$P ( \tilde { \rho }, M, c _ { V } ) = P ( D, M, c _ { V } ) \left | \frac { \mathrm d \, \tilde { \rho } ( D, M, c _ { V } ) } { \mathrm d D } \right | ^ { - 1 } \quad \text{ for } D < r _ { t } ( M, c _ { V } ) \,,$$

where ˜( ρ D, M, c V ) is the sub-halo density a distance D from the center of the sub-halo, Eq. (4), and P D,M,c ( V ) is the probability of finding the center of such a sub-halo at the distance D from Earth. This in turn can be calculated from the probability of finding a sub-halo with mass M P M , ( ) ; the probability of finding a concentration parameter c V for that sub-halo mass P c ( V | M ) ; and the probability of finding a sub-halo with a given mass and concentration at the distance D from the Earth P D M,c ( | V ) :

$$P ( D, M, c _ { V } ) = P ( M ) P ( c _ { V } | M ) P ( D | M, c _ { V } ) \,.$$

More concretely, the probability of finding a sub-halo of mass M is P M ( ) = 1 N sh d N/ M d , with d N/ M d shown in figure 1, and N sh the total number of sub-halos, given in [31]. The distribution of concentrations c V for sub-halos of mass M P c , ( V | M ) , is obtained using the tabulated results of Ref. [31]. Lastly, we assume for simplicity that the spatial distribution of sub-halos is independent of their mass and concentration. Hydrodynamical simulations suggest that baryonic disruption should preferentially deplete the most massive sub-halos at small galactocentric radii [73]. At larger radii (relevant in this context), simulations suggest instead only a weak dependence of the spatial distribution on the sub-halo properties. The probability of finding a sub-halo center at a distance D from the Earth can then be estimated as:

$$P ( D ) \equiv \frac { 1 } { N _ { \text{sh} } } \frac { \text{d} N _ { \text{sh} } ( D ) } { \text{d} D } = \frac { 4 \pi D ^ { 2 } \bar { n } _ { \text{sh} } ( D ) } { N _ { \text{sh} } } \,,$$

where N sh ( D ) is the number of sub-halos that can be found in the spherical shell located at distance D → D +d D from the Earth, with ¯ n sh ( D ) the average number density of sub-halos at that distance:

$$\bar { n } _ { \text{sh} } ( D ) = \frac { 1 } { 2 } \int _ { - 1 } ^ { 1 } n _ { \text{sh} } [ r ( D, \psi ) ] \, \text{d} \cos \psi \,.$$

Here, n sh ( r ) is the number density of sub-halos as a function of the galactocentric radius (given in Eq. (3)), which can be expressed in terms of the distance of the Sun to the Galactic center using r D, ψ ( ) = √ D 2 + r 2 glyph[circledot] -2 r glyph[circledot] D cos ψ , with r glyph[circledot] ≈ 8 5 kpc . and ψ the angular separation between the Galactic center and the sub-halo center.

Finally, one obtains the probability distribution for finding a contribution ˜ ρ to the local density (due to a single sub-halo) by integrating Eq. (18) over all sub-halo masses and concentrations:

$$P _ { \text{single} } ( \tilde { \rho } ) = \int _ { M _ { \min } } ^ { M _ { \max } } \mathrm d M \int _ { 0 } ^ { \infty } \mathrm d c _ { V } P ( \tilde { \rho }, M, c _ { V } ) \,.$$

This probability distribution is shown in the left panel of figure 5, assuming the sub-halo mass function of [31]. We take into account sub-halos with masses ranging from M min = 10 -6 M glyph[circledot] up to M max = 10 12 M glyph[circledot] . As apparent from the plot, the probability that a single sub-halo dominates the local density is small. This is to be expected, as the fraction of the Milky Way

Figure 5: Probability distribution of producing an overdensity ˜ ρ at the position of the Earth due to the passage through one particular sub-halo (left plot) and due to the passage through some subhalo (right plot) assuming the ShMF calculated in Ref. [31]. The figure also shows for comparison the corresponding probability distributions when considering only sub-halos with mass greater than 10 5 M glyph[circledot] , as resolved by numerical N-body simulations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-a21ba32d-8d30-4929-9f4f-77787a1979c7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-e71e90a0-6247-49fb-8dd5-af0f49fa701b.jpg)

volume occupied by a given individual sub-halo is also small. Indeed, the probability of having an enhancement of any size in the local dark matter density due to a particular single sub-halo, which we denote as p single , is p single ≡ ∫ d˜ ρ P single (˜) ρ ≈ 5 × 10 -17 . Nevertheless, due to the large number of sub-halos present in the galaxy, there is a fairly high probability that the Earth is currently immersed in a sub-halo. The probability of being immersed in any one sub-halo today is obtained from the Binomial distribution:

$$p _ { 1 } = N _ { \text{sh} } \, p _ { \text{single} } \, ( 1 - p _ { \text{single} } ) ^ { N _ { \text{sh} } - 1 } \approx N _ { \text{sh} } \, p _ { \text{single} } \, ( 1 - p _ { \text{single} } ) ^ { N _ { \text{sh} } } \sim 1 2 \% \,,$$

while the probability that the Earth is not immersed in any sub-halo is 3

$$p _ { 0 } = ( 1 - p _ { \text{single} } ) ^ { N _ { \text{sh} } } \approx \exp ( - N _ { \text{sh} } \cdot p _ { \text{single} } ) = 8 8 \% \,.$$

In the right panel of figure 5, we plot the probability P ρ (˜) ≡ ( p /p 1 single ) × P single (˜) ρ that some sub-halo enhances the local density by ˜ ρ , which is the relevant quantity for the purposes of direct detection at Earth. It follows from the plot that there is a 0.1% probability ( 10 -3 % probability) of finding a contribution to the local dark matter density from sub-halos which is equal to ρ loc SHM ( 10 ρ loc SHM ). The plots also show for reference the probability distributions considering only sub-halos with mass M ≥ 10 5 M glyph[circledot] , namely those that can be resolved and identified with current N -body simulations. Including also the lighter sub-halos, as done in this work, clearly enhances the probability of finding a larger DM density in the Solar System, owing to their larger number density. More concretely, the low mass sub-halos increase by a factor ∼ 3 ( ∼ 10 ) the probability of finding a contribution from the sub-halos to the total local DM density which is a factor ∼ 1 ( ∼ 10 ) larger than ρ loc SHM .

After determining the probability distributions for the sub-halo relative velocity and for finding an enhancement ˜ ρ in the local dark matter density, it is straightforward to calculate

3 In practice, we have used the approximation log(1 + x ) ≈ x for 0 &lt; x glyph[lessmuch] 1 to evaluate the probabilities.

the probability distribution of finding an increment I R in the scattering rate at a given direct detection experiment. These probability distributions are shown as solid blue lines in figure 6, for XENON1T (left panels) and CRESST (right panels), and for dark matter masses m DM = 5 GeV (top panels), 100 GeV (middle panels) and 10 TeV (bottom panels). For both experiments, we find that there is a probability glyph[lessorsimilar] 7 × 10 -4 of increasing the event rate by an O (1) factor.

To assess the impact of our assumptions on the probability distributions, we have also analyzed the cases where the mass function follows a simpler power law ∝ M -1 8 . or ∝ M -2 0 . , namely the upper and lower limits on α M ( ) in the right panel of figure 1. For these cases, we find that the probability of finding a ∼ 1 ( ∼ 10) enhancement on the rate is at most a factor ∼ 3 ( ∼ 4 ) larger than the one obtained with the full sub-halo mass function from Ref. [31]. We should also be careful that we are not 'over-counting' the density enhancement by allowing multiple sub-halos to contribute to the density at the present time. We find that this physically invalid configuration (of multiple overlapping sub-halos) has only a small probability and so can be neglected:

$$p _ { 2 } = \binom { N _ { \text{sh} } } { 2 } \, \left ( p _ { \text{single} } \right ) ^ { 2 } \left ( 1 - p _ { \text{single} } \right ) ^ { N _ { \text{sh} } - 2 } \approx \frac { 1 } { 2 } ( N _ { \text{sh} } \, p _ { \text{single} } ) ^ { 2 } \, p _ { 0 } = 0. 7 \% \,.$$

Lastly, let us note that the increment in the rate could be modified if the density profile of disrupted sub-halos is not described by an NFW profile, but rather by, for instance, an 'exponential profile' of the form

$$\rho ( r ) = \frac { \rho _ { 0 } } { r ^ { \gamma } } \exp \left ( - \frac { r } { R _ { b } } \right ),$$

with parameters ρ 0 , γ and R , as advocated in [94,95]. b Repeating our analysis using a different halo profile would require, for consistency, the use of the corresponding halo-mass function calculated along the lines of [31]. Unfortunately, this analysis is not available in the literature. Accordingly, investigating how the form of the halo profile affects the increment in the scattering rate is beyond the scope of our work.

## 4 Impact of sub-halos on the neutrino flux from the Sun

The idea that sub-structure can influence the capture rate of DM in the Sun was first introduced in [55]. Here, we extend this idea and build a probabilistic model based on realistic sub-halo distributions to assess the boost factor of the neutrino signal from the Sun. We note that the same rationale we present here applies also for DM capture inside the Earth and the possibility of observing a high-energy neutrino flux from the center of the Earth [96,97].

Dark matter particles crossing the Sun can scatter with Solar matter and lose energy, occasionally becoming gravitationally bound to the Sun. Once captured, and due to subsequent scatterings, these particles will sink to the Solar core and generate a DM overdensity where the annihilation rate is enhanced, possibly leading to observable signals at neutrino telescopes [98]. The capture rate of DM in the Sun can be determined from [96,99]:

$$C ( t ) = \sum _ { i } \int _ { 0 } ^ { R _ { \odot } } 4 \pi \, r ^ { 2 } \, d r \, \eta _ { i } ( r ) \, \int _ { v \leq v ^ { ( S u n ) } _ { \max, i } ( r ) } \mathrm d ^ { 3 } v \, \frac { F ( \vec { v }, t ) } { v ^ { 2 } } \, w ^ { 2 } ( r ) \int _ { m _ { \mathrm D M } v ^ { 2 } / 2 } ^ { 2 \mu _ { R _ { N _ { i } } } ^ { 2 } w ^ { 2 } ( r ) / m _ { N _ { i } } } \mathrm d E _ { R } \, \frac { \mathrm d \sigma _ { i } } { \mathrm d E _ { R } } ( w ( r ), E _ { R } ) \,.$$

Figure 6: Probability distributions of the increment in the scattering rate at XENON1T (left panel) and CRESST (right planel), for a dark matter mass of 5 GeV (top panel), 100 GeV (middle panel) and 10000 GeV (bottom panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-0925e6e9-0f9e-4b33-b498-1b293bd00360.jpg)

This rate is in general time-dependent, due to the time-dependent dark matter flux at the position of the Sun F glyph[vector] , ( v t ) (as described by Eq. (15)). In this expression, w r ( ) ≡ √ v 2 + v esc ( r ) 2 is the speed of the DM particle at a distance r from the Sun's center, given an asymptotic velocity v . The escape velocity from the Sun at radius r is v esc ( r ) , and v max is the maximum velocity for which capture is kinematically possible from scatterings with a nucleus i of mass m N i . In addition, η i ( r ) denotes the number density of the nuclei i inside the Sun, for which we take the Solar model AGSS09 [100], and d σ / i d E R denotes the differential scattering cross section with a nucleus species i . This cross section depends on a nuclear form factor F E ( R ) ; for spin-independent interactions we adopt the Helm form factors from Ref. [88, 101] and for spin-dependent interactions, the form factors from Ref. [102].

The time evolution of the number of captured dark matter particles inside the Sun is described by the following differential equation

$$\frac { \mathrm d N ( t ) } { \mathrm d t } = C ( t ) - C _ { A } N ( t ) ^ { 2 },$$

where the effects of evaporation can be neglected for dark matter heavier than a few GeV [103]. Here, C A characterizes the rate at which the number of dark matter particles is depleted due to annihilation, and which we calculate following Refs. [96,104].

We begin by neglecting the effect of the sub-halos, in which case C t ( ) = C SHM , the timeindependent capture rate of DM from the smooth Milky Way dark matter distribution, described by the Standard Halo Model. Assuming that N (0) = 0 at t = 0 , which we take as the time of formation of the Sun, the solution for the number of captured particles reads,

$$N ( t ) = \sqrt { \frac { C _ { \text{SHM} } } { C _ { A } } } \, \tanh \left ( \frac { t } { \tau } \right ) = N _ { \text{SHM} } \, \tanh \left ( \frac { t } { \tau } \right ) \,,$$

which tends to a constant N SHM = √ C SHM /C A when t glyph[greatermuch] τ . Here, τ ≡ 1 / √ C SHM C A denotes the equilibration time. In this case, the annihilation rate inside the Sun today is given by:

$$\Gamma _ { A } = \frac { 1 } { 2 } \, C _ { A } \, N ( t _ { 0 } ) ^ { 2 } = \frac { C _ { \text{SHM} } } { 2 } \, \tanh ^ { 2 } \left ( \frac { t _ { 0 } } { \tau } \right ) \,,$$

where t 0 is the time elapsed since the formation of the Sun, t 0 glyph[similarequal] 4 6 . Gyr. Assuming t 0 glyph[greatermuch] τ , the annihilation rate reduces to

$$\Gamma _ { A } = \frac { C _ { \text{SHM} } } { 2 } \,.$$

That is, the annihilation rate is entirely determined by the capture rate, and ultimately to the interaction strength of DM particles with the nuclei in the Solar interior. This tight relation between annihilation rate and capture rate is often used to translate the non-observation of an excess of high energy neutrinos from the Sun into constraints on the spin-independent and spin-dependent cross sections. 4

In this work we focus on the possibility of a time-varying DM flux at the Sun induced by the passage through a sub-halo. The capture rate is therefore also time-varying and has the form:

$$C ( t ) = C _ { \text{SHM} } + C _ { \text{sh} } ( t ) \Theta ( t - t _ { - } ) \Theta ( t _ { + } - t ) \,,$$

4 In some frameworks, however, the condition t 0 glyph[greatermuch] τ may not hold (e.g. in scenarios where the annihilation rate is p-wave suppressed [105]) and accordingly the annihilation rate gets suppressed by a factor tanh ( 2 t 0 /τ ) .

where C sh ( ) t describes the enhancement of the capture rate due to the Sun's passage through a single sub-halo between the times t -and t + .

To simplify the discussion, we approximate the time-dependent capture rate during the passage by its average value, namely

$$C _ { \text{sh} } ( t ) \simeq \langle C _ { \text{sh} } ( t ) \rangle \equiv \frac { 1 } { \Delta t } \int _ { t _ { - } } ^ { t _ { + } } \text{d} t \, C _ { \text{sh} } ( t ) \,,$$

where ∆ = t t + -t -. Under this approximation, the solution of Eq. (28) reads:

$$N ( t ) = \left \{ \sqrt { \frac { C _ { \text{SHM} } } { C _ { A } } } \tanh \left ( \frac { t } { \tau } \right ) & \text{if } t \leq t _ { - } \,, \\ \left \{ \sqrt { \frac { C _ { \text{SHM} } } { C _ { A } } } \frac { \tanh \left ( \frac { t } { \tau } \right ) } { \tau } \tanh \left ( \frac { t } { \tau } \right ) + \sqrt { 1 } + \frac { ( C _ { \text{SHM} } ( t ) ) } { C _ { \text{SHM} } } \tanh \left \{ \sqrt { 1 } + \frac { ( C _ { \text{SHM} } ( t ) ) } { \tau } \right ) } \right \} & \text{if } t _ { - } < t \leq t \,, \\ \frac { N ( t _ { + } ) + \sqrt { \frac { C _ { \text{SHM} } } { C _ { A } } } \tanh \left ( \frac { t } { \tau } \right ) } { \tau } \right ) & \text{if } t _ { - } < t \,, \\ \frac { N ( t _ { + } ) + \sqrt { \frac { C _ { \text{SHM} } } { C _ { A } } } \tanh \left ( \frac { t - t _ { + } } { \tau } \right ) } { 1 + N ( t _ { + } ) \left ( \sqrt { \frac { C _ { \text{SHM} } } { C _ { A } } } \right ) ^ { - 1 } \tanh \left ( \frac { t - t _ { + } } { \tau } \right ) } \right ] & \text{if } t > t _ { - } \,, \\ \text{where } \tau \text{ is the equilibration time in the absence of time-dependent contributions. } \, \text{We note}$$

where τ is the equilibration time in the absence of time-dependent contributions. We note that the impact of the passage of the Sun through the sub-halo can be characterized by only three dimensionless quantities: the duration of the passage relative to the equilibration time, ∆ t/τ , the time elapsed since the Sun abandoned the sub-halo relative to the equilibration time ( t 0 -t + ) /τ , and the time-averaged increment in the capture rate due to the passage through the sub-halo, 〈 C sh ( ) t 〉 /C SHM .

The two relevant parameters ∆ t/τ and 〈 C sh ( ) t 〉 /C SHM are in turn calculable given a particle DM framework, given the characteristics of the smooth halo and sub-halo components, and given the impact parameter of the Sun when entering the sub-halo. More concretely, for an impact parameter L , and under the approximation that the Sun moves inside the sub-halo following a straight line with constant speed v rel , it can be shown that the Sun travels in a complete passage a distance ∆ d = 2 √ r 2 t -L 2 in a time ∆ = 2 t √ r 2 t -L /v 2 rel . Therefore, ∆ t/τ ≤ 2 r /v t rel τ . For v rel = O (100) km/s, and taking conservatively τ glyph[greaterorsimilar] 10 6 years, one obtains v rel τ glyph[greaterorsimilar] 30 pc. Assuming that the sub-halos crossed by the Sun are much smaller than 30 pc 5 , so that ∆ t/τ glyph[lessmuch] 1 , the solution simplifies to:

$$N ( t ) & \simeq N _ { \text{SHM} } \left [ \frac { ( N ( t _ { - } ) + \Delta N ) / N _ { \text{SHM} } + \tanh \frac { t - t _ { - } } { \tau } } { 1 + ( N ( t _ { - } ) + \Delta N ) / N _ { \text{SHM} } \cdot \tanh \frac { t - t _ { - } } { \tau } } \right ]. \\. &. &. &. &. &. &. &. &. &..$$

Here, N t ( -) is the number of captured particles just before entering the sub-halo and N SHM = √ C SHM /C A is the equilibrium number of captured particles in the SHM. We also define ∆ N as

$$\frac { \Delta N \equiv \left < C _ { \text{sh} } ( t ) \right > \Delta t \,,$$

5 This is typically true for the most abundant sub-halos, with masses below M glyph[lessorsimilar] 10 3 M glyph[circledot] .

( )

N t

0

1

2

3

4

5

t

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-271d8a5b-1e6e-4af6-afbd-52b8502176f6.jpg)

[Gyr]

Figure 7: Sketch of the time evolution of the number of captured dark matter particles N t ( ) , and its value at the present time N t ( 0 ) , assuming a passage of the Sun through a single sub-halo at the time t -(left plot) and assuming successive passages through two sub-halos at times t 2 and t 1 (right plot). For comparison, we also show the corresponding values of N t ( ) and N t ( 0 ) for the Standard Halo Model.

which is interpreted as the increase in the number of DM particles captured during the passage: ∆ N ≡ N t ( + ) -N t ( -) . The time evolution of the number of captured dark matter particles is sketched in figure 7. From this it is apparent that the number of captured particles, and accordingly the neutrino flux from annihilations, can be significantly enhanced if the passage through the sub-halo was recent, ( t 0 -t -) /τ glyph[lessmuch] 1 , and if the number of particles captured during the passage was large ∆ N glyph[greatermuch] N t ( -) . A recent passage through the sub-halo suggests that capture and annihilation were in equilibrium inside the Sun before entering the sub-halo, therefore, it is plausible that N t ( -) = N SHM .

The time-averaged increment in the capture rate can be calculated inserting Eq. (15) into Eq. (27). For small mass sub-halos, for which the velocity distribution in the Solar frame can be approximated by Eq. (14), one finds

$$\Delta N \equiv C _ { \delta ( \vec { v } + \vec { v } _ { \text{rel} } ) } ^ { \rho _ { \text{SHM} } } \frac { \langle \rho _ { \text{sh} } ^ { \text{loc} } [ r ( t ) ] \rangle } { \rho _ { \text{SHM} } } \Delta t \,,$$

where C ρ SHM δ glyph[vector] ( v + glyph[vector] v rel ) denotes the capture rate for a DM stream with density ρ SHM and velocity glyph[vector] v rel = ( glyph[vector] v glyph[circledot] -glyph[vector] v sh ) with respect to the Solar frame, shown in figure 8 relative to C SHM . The time-averaged quantity 〈 ρ loc sh [ r t ( )] 〉 is related to the average mass density observed by the Sun along its path. For impact parameter L , and assuming that the Sun moves inside the sub-halo with constant velocity v rel , one finds

$$\begin{array} { c } \text{from conversion accuracy } \text{recv}, \, \text{and} \, \text{max} \\ \langle \rho _ { \text{sh} } [ r ( t ) ] \rangle \simeq \frac { 1 } { \Delta t } \int _ { t _ { - } } ^ { t _ { + } } \text{d} t \, \rho _ { \text{sh} } ^ { \text{loc} } [ r ( t ) ] = \frac { 1 } { 2 \sqrt { r _ { t } ^ { 2 } - L ^ { 2 } } } \int _ { - \sqrt { r _ { t } ^ { 2 } - L ^ { 2 } } } ^ { \sqrt { r _ { t } ^ { 2 } - L ^ { 2 } } } \text{d} x \, \rho _ { \text{sh} } ^ { \text{loc} } [ \sqrt { x ^ { 2 } + L ^ { 2 } } ] \,, \\........................................................................................................................................................................................................ \end{array}$$

where we have used that ∆ = 2 t √ r 2 t -L /v 2 rel . For a sub-halo with an NFW profile, Eq. (4),

N t

(

0

)

N

SHM

SHM

SHM+sub-halo

∆

N

t

-

t

0

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-658c1314-c02a-4895-8c4d-e0b77eef3834.jpg)

|

|

Figure 8: Capture rate at the Sun for dark matter streams with velocity | glyph[vector] v rel | , relative to the capture rate assuming the Standard Halo Model, for dark matter masses m DM = 5 , 100 and 10000 GeV.

## we find,

$$\langle \rho _ { \text{sh} } ^ { \text{loc} } [ r ( t ) ] \rangle = \frac { \rho _ { s } } { ( 1 - c _ { V } ^ { 2 } z ^ { 2 } ) } \left [ \frac { 2 } { c _ { V } \sqrt { 1 - z ^ { 2 } } \sqrt { 1 - c _ { V } ^ { 2 } z ^ { 2 } } } \, \text{aranh} \left ( \frac { \sqrt { 1 - z ^ { 2 } } \sqrt { 1 - c _ { V } ^ { 2 } z ^ { 2 } } } { ( 1 + z ) ( 1 + c _ { V } z ) } \right ) - \frac { 1 } { 1 + c _ { V } } \right ] \,,$$

with z ≡ L/r t . In certain limiting cases, this expression can be simplified to:

$$\langle \rho _ { \text{sh} } ^ { \text{oc} } [ r ( t ) ] \rangle \simeq \begin{cases} \frac { \rho _ { s } } { c _ { V } ( 1 + c _ { V } ) ^ { 2 } } \left [ 1 + \frac { 2 } { 3 } \left ( \frac { 1 + 3 c _ { V } } { 1 + c _ { V } } \right ) \left ( 1 - \frac { L } { r _ { t } } \right ) \right ] & \text{for $r_{t}> L \gg r_{s}$} \,, \\ - \frac { \rho _ { s } \left [ \log \frac { L } { r _ { s } } + \frac { c _ { V } } { 1 + c _ { V } } - \log \frac { c _ { V } } { 1 + c _ { V } } - \log 2 \right ] & \text{for $L \ll r_{s}$} \,. \end{cases}$$

The mean density diverges for L → 0 , as expected for a cuspy profile, which in turn translates into a large enhancement in the capture rate if the Sun traverses the innermost part of the sub-halo. On the other hand, typical sub-halos have a size larger than the scale radius, since the concentration parameter takes values c V ∼ 1 -35 . Therefore, in most instances one expects L glyph[greatermuch] r s .

Finally, mirroring our discussion in section 3, we define the increment in the annihilation rate inside the Sun (and in the neutrino flux) due to the passage through the sub-halo as:

$$\mathcal { I } _ { \Gamma } = \frac { \Gamma ( t _ { 0 } ) } { \Gamma _ { \text{SHM} } ( t _ { 0 } ) } - 1 \simeq \left [ \frac { ( N ( t _ { - } ) + \Delta N ) / N _ { \text{SHM} } + \tanh \frac { t _ { 0 } - t _ { - } } { \tau } } { 1 + ( N ( t _ { - } ) + \Delta N ) / N _ { \text{SHM} } \tanh \frac { t _ { 0 } - t _ { - } } { \tau } } \right ] ^ { 2 } \, \tanh \left ( \frac { t _ { 0 } } { \tau } \right ) ^ { - 2 } - 1 \,, \quad \ ( 4 1 )$$

where the annihilation rate expected from the smooth Milky Way dark matter distribution is given by Γ SHM ( t 0 ) = 1 2 C SHM tanh( t 0 /τ ) . We show in figure 9 contours of I Γ in the parameter space of ∆ N and ( t 0 -t -) /τ , assuming for simplicity that capture and annihilation were in equilibrium before entering the sub-halo, so that N t ( -) = N SHM and Γ SHM ( t -) = 1 2 C SHM . We furthermore assume equilibrium for the smooth Milky Way dark matter halo, i.e. Γ SHM ( t 0 ) = 1 2 C SHM . As expected from our previous discussion, the largest enhancement occurs for recent passages and for ∆ N glyph[greatermuch] N t ( -) . This occurs especially if the sub-halo is very dense; if the

Figure 9: Increment in the annihilation rate I Γ relative to the expectations from the Standard Halo Model in the parameter space spanned by ( t 0 -t -) /τ (the time elapsed since the passage of the Sun through the sub-halo, relative to the equilibration time) and ∆ N/N SHM (the number of particles in the sub-halo captured during the passage, relative to the number of particles in the SHM when equilibration is reached).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-99f8df89-467b-4183-9818-9385a229c8fa.jpg)

impact parameter is small; or if the speed of the Sun relative to the sub-halo is small. We remark that for DM scenarios where the number of captured DM particles from the smooth halo is below the equilibrium value ( i.e. when τ glyph[greatermuch] t 0 ), the passage of the Sun through one or more sub-halos can significantly enhance the total number of captured DM particles.

As in section 3, one can estimate the probability of having a certain increment I Γ in the annihilation rate inside the Sun. The probability for finding a certain value of ∆ N from the passage of the Sun through a sub-halo of mass M , concentration c V , and velocity glyph[vector] v rel = glyph[vector] v glyph[circledot] -glyph[vector] v sh is related to the probability of a passage with impact parameter L by

$$P ( \Delta N, M, c _ { V }, v _ { \text{rel} } ) = P ( L, M, c _ { V }, v _ { \text{rel} } ) \left | \frac { \text{d} \Delta N ( L, M, c _ { V }, v _ { \text{rel} } ) } { \text{d} L } \right | ^ { - 1 } \quad \text{for } L < r _ { t } \,. \quad \ \ ( 4 2 )$$

Note that ∆ N depends on L through 〈 ρ loc sh [ r t ( )] 〉 , given in Eqs. (38,39), and through ∆ t ≡ 2 √ R 2 vir -L /v 2 sh . The probability for the Sun traversing a sub-halo with impact parameter L can be calculated using a similar rationale as in section 3 (c.f. Eq. (19)):

$$P ( L, M, c _ { V }, v _ { \text{rel} } ) = P ( M ) P ( c _ { V } | M ) P ( v _ { \text{rel} } ) P ( L | M, c _ { V }, v _ { \text{rel} } ) \,.$$

Here, P ( M ) and P c ( V | M ) were already discussed in section 3, and P v ( rel ) is given by the Maxwell-Boltzmann distribution in Eq. (13) 6 . To determine P L M,c ( | V , v rel ) we note that the probability of the Sun passing through a sub-halo with impact parameter L is equal to the probability of finding a sub-halo center at a perpendicular distance L from the Sun's path during the crossing time ∆ t . The number of sub-halos in the cylindrical shell of width d L located at a perpendicular distance L from the Sun and with length equal to v rel ∆ t is d N sh ( L, M, c V , v rel ) glyph[similarequal] 2 πL d Lv rel ∆ ¯ t n sh ( L ) . Assuming, as we did in section 3, that the spatial

6 With the caveat that we are not interested in the direction of the sub-halo velocity, only its magnitude, so that in practice we draw | glyph[vector] v rel | from the distribution f ( v ) = v 2 ∫ f SHM ( glyph[vector] v ) dΩ v .

distribution of sub-halos is independent of their mass and concentration, one then obtains

$$P ( L | M, c _ { V }, v _ { \text{rel} } ) \equiv \frac { 1 } { N _ { \text{sh} } } \frac { \text{d} N _ { \text{sh} } ( L, M, c _ { V }, v _ { \text{rel} } ) } { \text{d} L } \simeq \frac { 2 \pi L v _ { \text{rel} } \Delta t \, \bar { n } _ { \text{sh} } ( L ) } { N _ { \text{sh} } } \,.$$

Here, N sh the total number of sub-halos and ¯ n sh ( L ) the average number of sub-halos at a perpendicular distance L , written in terms of the galactocentric radius r as: 7

$$\bar { n } _ { \text{sh} } ( L ) = \frac { 1 } { 2 \pi } \int _ { 0 } ^ { 2 \pi } n _ { \text{sh} } \left [ r ( L, \psi ) \right ] \, \text{d} \psi \,.$$

Finally, one obtains the probability distribution of finding a certain ∆ N from the passage through a single sub-halo by integrating over all possible sub-halo masses, concentration parameters and velocities:

$$P _ { \text{single} } ( \Delta N ) = \int _ { M _ { \min } } ^ { M _ { \max } } \mathrm d M \int _ { 0 } ^ { \infty } \mathrm d c _ { V } \int _ { 0 } ^ { v _ { \max } } \mathrm d v _ { \text{rel} } P ( \Delta N, M, c _ { V }, v _ { \text{rel} } )$$

The resulting distributions are shown in figure 10 for capture through the spin-independent (left panels) or the spin-dependent (right panels) interactions, for DM masses m DM = 5 GeV (top panels), 100 GeV (middle panels) and 10 TeV (bottom panels). We show the probability distributions for ∆ N/N SHM multiplied by a factor τ/t 0 in order to emphasize the impact of equilibration. If DM particles annihilate efficiently, then the equilibrium number of captured particles N SHM is suppressed. This means that if we decrease the equilibration time τ , the number of particles captured from the sub-halo ∆ N is increased relative to the equilibrium number. However, a short equilibration time also means that enhancements in the number of captured particles are quickly annihilated away. The factor t 0 /τ counts the number of equilibration times over the lifetime of the Sun, meaning that the factor τ/t 0 roughly corrects for the effects of equilibration after the sub-halo passage. In order to assess the impact of varying the halo mass function, we additionally show the results for the power-law mass functions d N/ M d ∝ M -α with fixed exponents α = 1 8 . and α = 2 0 . . For the shallowest sub-halo mass function with α = 1 8 . , we observe the largest ∆ N as we have more heavy sub-halos than for the other two. Conversely, we find the smallest enhancement for α = 2 0 . . Finally, we find that the distributions are flatter and reach farther towards large ∆ N for heavy dark matter. The reason for this is that the maximal velocity at which a DM particle can be captured by the Sun is small in this case. Most of the sub-halos therefore lead to no additional DM particles being captured in the Sun. However, if the stream is slow enough, ∆ N/N SHM can be large as the capture from the smooth Milky Way DM distribution is small. Such sub-halo crossings contribute to the flat tail of the distribution.

In the case of direct detection, we demonstrated that the largest contribution to a local overdensity comes from low mass sub-halos (see Figure 5). In the case of solar capture, it is more challenging to identify which sub-halos are the most relevant. This is because the size of a sub-halos contribution depends not only on its mass and concentration, but also on its impact parameter L and relative velocity with respect to the Earth v rel . As a result, sub-halos with the same mass and concentration parameter can give rise to very different increments in the annihilation rate, I Γ .

7 Note the slightly different definition compared to Eq. (21), due to the different geometry.

Figure 10: Probability distribution of the dimensionless and cross section independent combination of parameters (∆ N/N SHM ) · ( τ/t 0 ) , assuming the spin-independent interaction (left panel) or the spindependent interaction (right-panel) for a dark matter mas of 5 GeV (top panel) 100 GeV (middle panel) and 10000 GeV (bottom panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-310e58f7-e98f-4e11-b3c8-f72fda9370d8.jpg)

So far we have studied the effect on the annihilation rate of the crossing through a single sub-halo, concluding that the effects of the sub-halo passage remain significant for a long period of time ∼ τ . It is then plausible that the Sun could have traversed more than one sub-halo during this time, resulting in an enhanced accumulation of dark matter particles inside the Sun, and in a larger probability of an increment in the annihilation rate.

We consider that at the time t = 0 the number of captured dark matter particles in the Sun is N (0) = 0 , and that there have been N sub-halo passages at the times t i , i = 1 ... N (ordered such that 0 &lt; t N &lt; ... &lt; t 1 &lt; t 0 ) with duration ∆ t /τ i glyph[lessmuch] 1 , each producing an 'instantaneous' change in the number of captured DM particles ∆ N i , as sketched in figure 7. The total number of sub-halo encounters N in a total time t 0 can be calculated as:

$$\mathcal { N } = \frac { t _ { 0 } } { t _ { \text{orb} } } N _ { \text{sh} } \int _ { 0 } ^ { \infty } \mathrm d \Delta N \, P _ { \text{single} } ( \Delta N ) \,,$$

where the integral corresponds to the probability of a crossing a single sub-halo during one Solar orbit around the Galaxy. With an orbital period of t orb = 250Myr , the factor of t 0 /t orb ≈ 18 simply counts the number of Galactic orbits during the Sun's lifetime. We find the expected number of sub-halo crossings to be N ∼ 100 , which we assume are uniformly distributed in time. Figure 11 shows the number of sub-halo crossings as a function of sub-halo mass and impact parameter L . Crossing events are dominated by small sub-halos due to their larger abundance, while close passages are very rare. Indeed, crossings at sub-parsec distances from the sub-halo center are expected to occur less than once during the age of the Sun, justifying our assumption that L glyph[greatermuch] r s for the sub-halos of interest (see Eq. (40)).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-47919d3d-9424-4beb-8ce3-43e8a74cb7a9.jpg)

/circledot

Figure 11: Number of crossings of halos with mass M (per decade in sub-halo mass) with close passage distance less than L . The horizontal dashed line corresponds to a single crossing during the age of the Sun t 0 .

The number of captured dark matter particles today is then:

$$N ( t _ { 0 } ) \simeq N _ { \text{SHM} } \left [ \frac { ( N ( t _ { 1 } ) + \Delta N _ { 1 } ) / N _ { \text{SHM} } + \tanh \frac { t _ { 0 } - t _ { 1 } } { \tau } } { 1 + ( N ( t _ { 1 } ) + \Delta N _ { 1 } ) / N _ { \text{SHM} } \cdot \tanh \frac { t _ { 0 } - t _ { 1 } } { \tau } } \right ],$$

with

$$N ( t _ { i } ) \simeq N _ { \text{SHM} } \left [ \frac { ( N ( t _ { i + 1 } ) + \Delta N _ { i + 1 } ) / N _ { \text{SHM} } + \tanh \frac { t _ { 0 } - t _ { i + 1 } } { \tau } } { 1 + ( N ( t _ { i + 1 } ) + \Delta N _ { i + 1 } ) / N _ { \text{SHM} } \cdot \tanh \frac { t _ { 0 } - t _ { i + 1 } } { \tau } } \right ] \quad \text{for $i=1,\dots,\mathcal{N}$} \,. \quad ( 4 9 )$$

Finally, the increment in the annihilation rate reads

$$\mathcal { I } _ { \Gamma } = \frac { \Gamma ( t _ { 0 } ) } { \Gamma _ { \text{SHM} } ( t _ { 0 } ) } - 1 \simeq \left ( \frac { N ( t _ { 0 } ) } { N _ { \text{SHM} } \tanh \left ( \frac { t _ { 0 } } { \tau } \right ) } \right ) ^ { 2 } - 1$$

We show the results of this calculation in figures 12, 13 and 14 for different DM masses. For each interaction type, we calculate the annihilation boost for three different interaction strengths, ranging down to the interaction cross sections currently probed by direct detection experiments [7]. Similar to figure 10, we additionally include the results for power-law sub-halo mass functions d N/ M d ∝ M -α with exponents α = 1 8 . and α = 2 0 . . Again, the probability to find large increments is usually greatest for α = 1 8 . as this predicts the highest number of heavy sub-halos and smallest for α = 2 0 . . Furthermore, larger values of I Γ are more likely to be achieved for large cross sections for the discrete examples we studied.

As apparent from Fig. 10, we expect larger ∆ N/N SHM as the factor ( τ/t 0 ) decreases with increasing cross section. Therefore, the ratio of captured particles ∆ N/N SHM increases while the equilibration time decreases, which implies that the importance of recent encounters rises while the history of sub-halo passages becomes less relevant. This enhances the probability to find large increments as most of them are due to recent encounters. Ultimately, the distribution of I Γ will become similar to the results of section 3 as soon as the cross section gets large enough such that the only relevant sub-halo encounters are those that take place in the very recent past.

When decreasing the cross section even further than in this work, large boost factors are possible again. In this regime, capture and annihilation processes are far away from equilibrium and the number of particles in the Sun can be dominated by capture from sub-halos. However, we do not show this here as those cross sections cannot be probed even with future upgrades of neutrino telescopes. We note that the formalism developed in this work is applicable for all cross sections.

## 5 Conclusions

In this paper, we have studied the impact of sub-halos on the interpretation of direct dark matter searches. Based on the sub-halo mass functions presented in Ref. [31], we calculated the probability of the Earth currently being inside a sub-halo, as well as the enhancement of the local dark matter density due to this. Furthermore, we investigated the sub-halo encounters of the Sun during its lifetime. For this, we calculated the average number of sub-halo crossings as well as the probability distribution of the enhancement due to one individual encounter.

We find that the impact on direct detection experiments can be large and should be taken into account when analyzing the results of direct detection experiments. The probability of a sub-halo enhancing the local dark matter density by a factor O (1) is ∼ 10 -3 (figure 5). Such large enhancements are far less likely (and in some cases impossible) for the neutrino signal from the Sun, which is enhanced only at a sub-percent level (figures 12, 13 and 14). We see that the probability of finding a large boost factor of O (1) is glyph[lessorsimilar] 10 -5 , suggesting that single large encounters will average away over the Sun's lifetime. Thus, we conclude that perturbations of the dark matter density and velocity distribution due to sub-halos can have a substantial impact on direct detection experiments while the neutrino signal from the Sun is largely unaffected.

Figure 12: Probability distributions of the increment in the annihilation rate inside the Sun for a dark matter mass of 5 GeV for the spin-independent cross sections σ SI = 10 -42 cm 2 , 10 -44 cm 2 and 10 -46 cm 2 (left panel, from top to bottom), and for the spin-dependent cross sections σ SD = 10 -38 cm 2 , 10 -40 cm 2 and 10 -42 cm 2 (right panel, from top to bottom).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-592777af-5428-40f9-85d8-1d35d56c053d.jpg)

Figure 13: Same as Fig. 12, but for a dark matter mass of 100 GeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-587e40d3-1d07-4de7-939a-6bada6b849dd.jpg)

Figure 14: Same as Fig. 12, but for a dark matter mass of 10000 GeV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00747v2-figure-f1c30bd4-5de7-4b19-a015-f02eddd4be2d.jpg)

Our results are broadly consistent with those of previous studies on the impact of substructure on direct detection. For example, Ref. [50] estimates a roughly 1% probability of an O (1) enhancement in the local density from sub-structure. Though our results suggest a somewhat smaller probability, there is perhaps a factor of a few uncertainty coming from the precise sub-halo mass function, as we point out in figure 6. The authors of Ref. [54] also built an analytic model for Milky Way sub-halos, inspired by N-body simulations. Here, we use more recent, state-of-the-art determinations of the properties of the sub-halo population, though we find similar results. In the case of solar capture, Ref. [55] demonstrated that sufficiently dense sub-structures, traversed for a sufficiently long time, may have a substantial impact on the neutrino signal from annihilating DM in the Sun. We have demonstrated that such encounters are rare and that a realistic population of sub-structures is unlikely to give a large enhancement.

This work alleviates a major uncertainty of dark matter searches with neutrino telescopes [106,107] and emphasizes the impact of astrophysical uncertainties on direct detection experiments. In particular, given that deviations from the Standard Halo Model are expected due to sub-halos, then it will be necessary to properly account for these uncertainties when producing constraints from direct detection experiments. Of course, given a future signal, it may be possible to directly reconstruct the local DM velocity distribution (as proposed in e.g. [108-113]). With this, it should be possible to determine whether the local DM density does indeed receive an enhancement due to a sub-halo, as well as hinting at the properties of the sub-halo itself.

While we have made a number of assumptions about the properties and distributions of DM substructure, the formalism we have developed is applicable more generally. For example, it was recently claimed that the density profile of disrupted sub-halos is not well described by an NFW profile but instead by an exponential profile ρ r ( ) ∝ r -γ exp( -r/ R b ) [94, 95]. The parameters γ and R b are drawn from generalized normal distributions whose mean values are determined from ELVIS [114] and Via Lactea II [115] simulations for different bins of sub-halo masses. As these two simulations are able to resolve sub-halos only down to 10 7 M glyph[circledot] , one has to rely on extrapolations to assess sub-halos as light as 10 -5 M glyph[circledot] . In this work, we do not use this profile since we found that the largest contribution to the enhancement of the local dark matter density at Earth and to the enhancement of the neutrino signal from the Sun is due to sub-halos in the extrapolated mass range. However, we note that using the exponential density profile for our calculation would be straightforward and could be easily implemented as soon as the resolution is sufficient.

We have focused in this work on the sub-halo distribution for standard Cold Dark Matter. Our approach could equally well be applied to other DM candidates, such as Warm Dark Matter [116] and Self-interacting Dark Matter [117, 118], as well as DM with more exotic clustering properties [36, 37]. We have also focused on direct searches for WIMP DM, though our framework may also be applied to axion searches, where axions are expected to form bound 'mini-clusters' [119], potentially giving an enhancement to the local density. This work indicates that the observation of DM substructures in direct detection experiments may still be promising, which would shed further light on the nature of the elusive dark matter particle.

## Acknowledgments

This work has been supported by the DFG cluster of excellence ORIGINS and by the Collaborative Research Center SFB1258. BJK gratefully acknowledges partial support from the Netherlands Organisation for Scientific Research (NWO) through the VIDI research program

'Probing the Genesis of Dark Matter' (680-47-532). We thank the Munich Institute for Astroand Particle Physics (MIAPP) where part of this work was developed. We also thank Ranjan Laha for discussions on the Solar capture of substructure, as well as for pointing us towards some of the relevant literature; Shin'ichiro Ando for sharing the numerical results of Ref. [31]; and Adam Coogan for helpful discussions on the contribution of sub-halos to the local DM density.

## References

- [1] G. Bertone, D. Hooper, and J. Silk, Particle dark matter: Evidence, candidates and constraints , Phys. Rept. 405 (2005) 279-390, [ hep-ph/0404175 ].
- [2] G. Bertone and M. P. Tait, Tim, A new era in the search for dark matter , Nature 562 (2018), no. 7725 51-56, [ arXiv:1810.01668 ].
- [3] L. Roszkowski, E. M. Sessolo, and S. Trojanowski, WIMP dark matter candidates and searches-current status and future prospects , Rept. Prog. Phys. 81 (2018), no. 6 066201, [ arXiv:1707.06277 ].
- [4] G. Jungman, M. Kamionkowski, and K. Griest, Supersymmetric dark matter , Phys. Rept. 267 (1996) 195-373, [ hep-ph/9506380 ].
- [5] F. Kahlhoefer, Review of LHC Dark Matter Searches , Int. J. Mod. Phys. A32 (2017), no. 13 1730006, [ arXiv:1702.02430 ].
- [6] J. M. Gaskins, A review of indirect searches for particle dark matter , Contemp. Phys. 57 (2016), no. 4 496-525, [ arXiv:1604.00014 ].
- [7] M. Schumann, Direct Detection of WIMP Dark Matter: Concepts and Status , arXiv:1903.03026 .
- [8] J. R. Primack, Dark matter and structure formation , in Formation of Structure in the Universe, Proceedings of the Jerusalem Winter School 1996, 1997. astro-ph/9707285 .
- [9] J. Diemand and B. Moore, The structure and evolution of cold dark matter halos , Adv. Sci. Lett. 4 (2011) 297-310, [ arXiv:0906.4340 ].
- [10] C. S. Frenk and S. D. M. White, Dark matter and cosmic structure , Annalen Phys. 524 (2012) 507-534, [ arXiv:1210.0544 ].
- [11] L. E. Strigari, J. S. Bullock, M. Kaplinghat, J. D. Simon, M. Geha, B. Willman, and M. G. Walker, A common mass scale for satellite galaxies of the Milky Way , Nature 454 (2008) 1096-1097, [ arXiv:0808.3772 ].
- [12] G. Besla, The Orbits and Total Mass of the Magellanic Clouds , arXiv e-prints (Nov, 2015) arXiv:1511.03346, [ arXiv:1511.03346 ].

| [13]   | A. A. Thoul and D. H. Weinberg, Hydrodynamic simulations of galaxy formation. 2. Photoionization and the formation of low mass galaxies , Astrophys. J. 465 (1996) 608, [ astro-ph/9510154 ].                                                         |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [14]   | J. S. Bullock, A. V. Kravtsov, and D. H. Weinberg, Reionization and the abundance of galactic satellites , Astrophys. J. 539 (2000) 517, [ astro-ph/0002214 ].                                                                                        |
| [15]   | R. S. Somerville, Can photoionization squelching resolve the sub-structure crisis? , Astrophys. J. 572 (2002) L23-L26, [ astro-ph/0107507 ].                                                                                                          |
| [16]   | E. O. Nadler, Y.-Y. Mao, G. M. Green, and R. H. Wechsler, Modeling the Connection between Subhalos and Satellites in Milky Way-like Systems , Astrophys. J. 873 (2019), no. 1 34, [ arXiv:1809.05542 ].                                               |
| [17]   | K. P. Zybin, M. I. Vysotsky, and A. V. Gurevich, The fluctuation spectrum cutoff in a neutralino dark matter scenario , Phys. Lett. A260 (1999) 262-268.                                                                                              |
| [18]   | S. Hofmann, D. J. Schwarz, and H. Stoecker, Damping scales of neutralino cold dark matter , Phys. Rev. D64 (2001) 083507, [ astro-ph/0104173 ].                                                                                                       |
| [19]   | V. Berezinsky, V. Dokuchaev, and Y. Eroshenko, Remnants of dark matter clumps , Phys. Rev. D77 (2008) 083519, [ arXiv:0712.3499 ].                                                                                                                    |
| [20]   | K. V. Johnston, D. N. Spergel, and C. Haydn, How lumpy is the milky way's dark matter halo? , Astrophys. J. 570 (2002) 656, [ astro-ph/0111196 ].                                                                                                     |
| [21]   | J. Diemand, M. Kuhlen, and P. Madau, Early supersymmetric cold dark matter substructure , Astrophys. J. 649 (2006) 1-13, [ astro-ph/0603250 ].                                                                                                        |
| [22]   | V. Springel, J. Wang, M. Vogelsberger, A. Ludlow, A. Jenkins, A. Helmi, J. F. Navarro, C. S. Frenk, and S. D. M. White, The Aquarius Project: the subhalos of galactic halos , Mon. Not. Roy. Astron. Soc. 391 (2008) 1685-1711, [ arXiv:0809.0898 ]. |
| [23]   | M. Boylan-Kolchin, V. Springel, S. D. M. White, and A. Jenkins, There's no place like home? Statistics of Milky Way-mass dark matter halos , Mon. Not. Roy. Astron. Soc. 406 (2010) 896, [ arXiv:0911.4484 ].                                         |
| [24]   | C. Giocoli, G. Tormen, R. K. Sheth, and F. C. van den Bosch, The Substructure Hierarchy in Dark Matter Haloes , Mon. Not. Roy. Astron. Soc. 404 (2010) 502-517, [ arXiv:0911.0436 ].                                                                  |
| [25]   | L. Gao, C. S. Frenk, M. Boylan-Kolchin, A. Jenkins, V. Springel, and S. D. M. White, The statistics of the subhalo abundance of dark matter haloes , Mon. Not. Roy. Astron. Soc. 410 (2011) 2309, [ arXiv:1006.2882 ].                                |
| [26]   | M. Kamionkowski, S. M. Koushiappas, and M. Kuhlen, Galactic Substructure and Dark Matter Annihilation in the Milky Way Halo , Phys. Rev. D81 (2010) 043532, [ arXiv:1001.3144 ].                                                                      |

- [27] T. Ishiyama, J. Makino, S. Portegies Zwart, D. Groen, K. Nitadori, S. Rieder, C. de Laat, S. McMillan, K. Hiraki, and S. Harfst, The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos , Astrophys. J. 767 (2013) 146, [ arXiv:1101.2020 ].
- [28] M. A. Sanchez-Conde and F. Prada, The flattening of the concentration-mass relation towards low halo masses and its implications for the annihilation signal boost , Mon. Not. Roy. Astron. Soc. 442 (2014), no. 3 2271-2277, [ arXiv:1312.1729 ].
- [29] M. Stref and J. Lavalle, Modeling dark matter subhalos in a constrained galaxy: Global mass and boosted annihilation profiles , Phys. Rev. D95 (2017), no. 6 063003, [ arXiv:1610.02233 ].
- [30] C. Okoli, J. E. Taylor, and N. Afshordi, Searching for dark matter annihilation from individual halos: uncertainties, scatter and signal-to-noise ratios , JCAP 1808 (2018), no. 08 019, [ arXiv:1711.05271 ].
- [31] N. Hiroshima, S. Ando, and T. Ishiyama, Modeling evolution of dark matter substructure and annihilation boost , Phys. Rev. D97 (2018), no. 12 123002, [ arXiv:1803.07691 ].
- [32] Gaia Collaboration, A. G. A. Brown, A. Vallenari, T. Prusti, J. H. J. de Bruijne, C. Babusiaux, C. A. L. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, and et al., Gaia Data Release 2. Summary of the contents and survey properties , A&amp;A 616 (Aug., 2018) A1, [ arXiv:1804.09365 ].
- [33] L. Necib, M. Lisanti, and V. Belokurov, Dark Matter in Disequilibrium: The Local Velocity Distribution from SDSS-Gaia , arXiv:1807.02519 .
- [34] K. Freese, P. Gondolo, H. J. Newberg, and M. Lewis, The effects of the Sagittarius dwarf tidal stream on dark matter detectors , Phys. Rev. Lett. 92 (2004) 111301, [ astro-ph/0310334 ].
- [35] C. A. O'Hare, C. McCabe, N. W. Evans, G. Myeong, and V. Belokurov, Dark matter hurricane: Measuring the S1 stream with dark matter detectors , Phys. Rev. D98 (2018), no. 10 103006, [ arXiv:1807.09004 ].
- [36] S. Nussinov and Y. Zhang, What If Dark Matter Clusters? , arXiv:1807.00846 .
- [37] D. M. Grabowska, T. Melia, and S. Rajendran, Detecting Dark Blobs , Phys. Rev. D98 (2018), no. 11 115020, [ arXiv:1807.03788 ].
- [38] J. Diemand, M. Kuhlen, and P. Madau, Dark matter substructure and gamma-ray annihilation in the Milky Way halo , Astrophys. J. 657 (2007) 262-270, [ astro-ph/0611370 ].
- [39] L. E. Strigari, S. M. Koushiappas, J. S. Bullock, and M. Kaplinghat, Precise constraints on the dark matter content of Milky Way dwarf galaxies for gamma-ray experiments , Phys. Rev. D75 (2007) 083526, [ astro-ph/0611925 ].
- [40] L. Pieri, G. Bertone, and E. Branchini, Dark Matter Annihilation in Substructures Revised , Mon. Not. Roy. Astron. Soc. 384 (2008) 1627, [ arXiv:0706.2101 ].

- [41] T. E. Jeltema and S. Profumo, Searching for Dark Matter with X-ray Observations of Local Dwarf Galaxies , Astrophys. J. 686 (2008) 1045, [ arXiv:0805.1054 ].
- [42] T. Ishiyama, Hierarchical Formation of Dark Matter Halos and the Free Streaming Scale , Astrophys. J. 788 (2014) 27, [ arXiv:1404.1650 ].
- [43] R. Bartels and S. Ando, Boosting the annihilation boost: Tidal effects on dark matter subhalos and consistent luminosity modeling , Phys. Rev. D92 (2015), no. 12 123508, [ arXiv:1507.08656 ].
- [44] D. Erkal and V. Belokurov, Properties of Dark Subhaloes from Gaps in Tidal Streams , Mon. Not. Roy. Astron. Soc. 454 (2015), no. 4 3542-3558, [ arXiv:1507.05625 ].
- [45] J. Bovy, D. Erkal, and J. L. Sanders, Linear perturbation theory for tidal streams and the small-scale CDM power spectrum , Mon. Not. Roy. Astron. Soc. 466 (2017), no. 1 628-668, [ arXiv:1606.03470 ].
- [46] N. Banik, G. Bertone, J. Bovy, and N. Bozorgnia, Probing the nature of dark matter particles with stellar streams , JCAP 1807 (2018), no. 07 061, [ arXiv:1804.04384 ].
- [47] N. Banik and J. Bovy, Effects of baryonic and dark matter substructure on the Pal 5 stream , Mon. Not. Roy. Astron. Soc. 484 (2019), no. 2 2009-2020, [ arXiv:1809.09640 ].
- [48] A. Bonaca, D. W. Hogg, A. M. Price-Whelan, and C. Conroy, The Spur and the Gap in GD-1: Dynamical evidence for a dark substructure in the Milky Way halo , arXiv:1811.03631 .
- [49] M. Danninger and C. Rott, Solar WIMPs unravelled: Experiments, astrophysical uncertainties, and interactive tools , Phys. Dark Univ. 5-6 (2014) 35-44, [ arXiv:1509.08230 ].
- [50] D. Stiff, L. M. Widrow, and J. Frieman, Signatures of hierarchical clustering in dark matter detection experiments , Phys. Rev. D64 (2001) 083516, [ astro-ph/0106048 ].
- [51] A. Helmi, S. D. M. White, and V. Springel, The Phase space structure of a dark matter halo: Implications for dark matter direct detection experiments , Phys. Rev. D66 (2002) 063502, [ astro-ph/0201289 ].
- [52] A. M. Green, WIMP direct detection and halo structure , astro-ph/0310215 . [IAU Symp.220,483(2004)].
- [53] M. Vogelsberger, S. D. M. White, A. Helmi, and V. Springel, The fine-grained phase-space structure of Cold Dark Matter halos , Mon. Not. Roy. Astron. Soc. 385 (2008) 236, [ arXiv:0711.1105 ].
- [54] M. Kamionkowski and S. M. Koushiappas, Galactic substructure and direct detection of dark matter , Phys. Rev. D77 (2008) 103509, [ arXiv:0801.3269 ].
- [55] S. M. Koushiappas and M. Kamionkowski, Galactic Substructure and Energetic Neutrinos from the Sun and the Earth , Phys. Rev. Lett. 103 (2009) 121301, [ arXiv:0907.4778 ].

- [56] A. M. Green, Astrophysical uncertainties on direct detection experiments , Mod. Phys. Lett. A27 (2012) 1230004, [ arXiv:1112.0524 ].
- [57] SDSS , X. X. Xue et. al., The Milky Way's Circular Velocity Curve to 60 kpc and an Estimate of the Dark Matter Halo Mass from Kinematics of 2400 SDSS Blue Horizontal Branch Stars , Astrophys. J. 684 (2008) 1143-1158, [ arXiv:0801.1232 ].
- [58] P. J. McMillan and J. J. Binney, The uncertainty in Galactic parameters , Mon. Not. Roy. Astron. Soc. 402 (2010) 934, [ arXiv:0907.4685 ].
- [59] J. Bovy, D. W. Hogg, and H.-W. Rix, Galactic masers and the Milky Way circular velocity , Astrophys. J. 704 (2009) 1704-1709, [ arXiv:0907.5423 ].
- [60] F. J. Kerr and D. Lynden-Bell, Review of galactic constants , Mon. Not. Roy. Astron. Soc. 221 (1986) 1023.
- [61] M. C. Smith et. al., The RAVE Survey: Constraining the Local Galactic Escape Speed , Mon. Not. Roy. Astron. Soc. 379 (2007) 755-772, [ astro-ph/0611671 ].
- [62] T. Piffl et. al., The RAVE survey: the Galactic escape speed and the mass of the Milky Way , Astron. Astrophys. 562 (2014) A91, [ arXiv:1309.4293 ].
- [63] A. M. Green, Astrophysical uncertainties on the local dark matter distribution and direct detection experiments , J. Phys. G44 (2017), no. 8 084001, [ arXiv:1703.10102 ].
- [64] Y. Wu, K. Freese, C. Kelso, P. Stengel, and M. Valluri, Uncertainties in Direct Dark Matter Detection in Light of Gaia , arXiv:1904.04781 .
- [65] J. I. Read, The Local Dark Matter Density , J. Phys. G41 (2014) 063101, [ arXiv:1404.1938 ].
- [66] N. Bozorgnia and G. Bertone, Implications of hydrodynamical simulations for the interpretation of direct dark matter searches , Int. J. Mod. Phys. A32 (2017), no. 21 1730016, [ arXiv:1705.05853 ].
- [67] N. W. Evans, C. A. J. O'Hare, and C. McCabe, Refinement of the standard halo model for dark matter searches in light of the Gaia Sausage , Phys. Rev. D99 (2019), no. 2 023012, [ arXiv:1810.11468 ].
- [68] S. Mandal, S. Majumdar, V. Rentala, and R. Basu Thakur, Observationally inferred dark matter phase-space distribution and direct detection experiments , Phys. Rev. D100 (2019), no. 2 023002, [ arXiv:1806.06872 ].
- [69] M. Stref, T. Lacroix, and J. Lavalle, Remnants of Galactic subhalos and their impact on indirect dark matter searches , Galaxies 7 (2019), no. 2 65, [ arXiv:1905.02008 ].
- [70] S. Ando, T. Ishiyama, and N. Hiroshima, Halo Substructure Boosts to the Signatures of Dark Matter Annihilation , Galaxies 7 (2019), no. 3 68, [ arXiv:1903.11427 ].
- [71] T. Ishiyama and S. Ando, The Abundance and Structure of Subhaloes near the Free Streaming Scale and Their Impact on Indirect Dark Matter Searches , arXiv:1907.03642 .

- [72] T. Ishiyama, M. Enoki, M. A. R. Kobayashi, R. Makiya, M. Nagashima, and T. Oogi, The ν 2 GC simulations: Quantifying the dark side of the universe in the Planck cosmology , Publ. Astron. Soc. Jap. 67 (2015), no. 4 61, [ arXiv:1412.2860 ].
- [73] Q. Zhu, F. Marinacci, M. Maji, Y. Li, V. Springel, and L. Hernquist, Baryonic impact on the dark matter distribution in Milky Way-sized galaxies and their satellites , Mon. Not. Roy. Astron. Soc. 458 (2016), no. 2 1559-1580, [ arXiv:1506.05537 ].
- [74] J. Einasto, 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 5 (1965) 87-100.
- [75] A. W. Graham, D. Merritt, B. Moore, J. Diemand, and B. Terzic, Empirical models for Dark Matter Halos. I. Nonparametric Construction of Density Profiles and Comparison with Parametric Models , Astron. J. 132 (2006) 2685-2700, [ astro-ph/0509417 ].
- [76] T. Kelley, J. S. Bullock, S. Garrison-Kimmel, M. Boylan-Kolchin, M. S. Pawlowski, and A. S. Graus, Phat ELVIS: The inevitable effect of the Milky Way's disk on its dark matter subhaloes , arXiv:1811.12413 .
- [77] J. Richings, C. Frenk, A. Jenkins, and A. Robertson, Subhalo destruction in the Apostle and Auriga simulations , arXiv e-prints (Nov, 2018) arXiv:1811.12437, [ arXiv:1811.12437 ].
- [78] J. F. Navarro, C. S. Frenk, and S. D. M. White, The Structure of cold dark matter halos , Astrophys. J. 462 (1996) 563-575, [ astro-ph/9508025 ].
- [79] J. F. Navarro, C. S. Frenk, and S. D. M. White, A Universal density profile from hierarchical clustering , Astrophys. J. 490 (1997) 493-508, [ astro-ph/9611107 ].
- [80] B. Allgood, R. A. Flores, J. R. Primack, A. V. Kravtsov, R. H. Wechsler, A. Faltenbacher, and J. S. Bullock, The shape of dark matter halos: dependence on mass, redshift, radius, and formation , Mon. Not. Roy. Astron. Soc. 367 (2006) 1781-1796, [ astro-ph/0508497 ].
- [81] N. E. Drakos, J. E. Taylor, and A. J. Benson, The phase-space structure of tidally stripped haloes , Mon. Not. Roy. Astron. Soc. 468 (2017), no. 2 2345-2358, [ arXiv:1703.07836 ].
- [82] O. N. Snaith, J. Bailin, A. Knebe, G. Stinson, J. Wadsley, and H. Couchman, Haloes at the ragged edge: the importance of the splashback radius , MNRAS 472 (Dec., 2017) 2694-2712, [ arXiv:1708.06181 ].
- [83] A. V. Kravtsov and S. Borgani, Formation of Galaxy Clusters , ARA&amp;A 50 (Sept., 2012) 353-409, [ arXiv:1205.5556 ].
- [84] J. Peñarrubia, A. J. Benson, M. G. Walker, G. Gilmore, A. W. McConnachie, and L. Mayer, The impact of dark matter cusps and cores on the satellite galaxy population around spiral galaxies , MNRAS 406 (Aug, 2010) 1290-1305, [ arXiv:1002.3376 ].

- [85] J. S. Bullock, T. S. Kolatt, Y. Sigad, R. S. Somerville, A. V. Kravtsov, A. A. Klypin, J. R. Primack, and A. Dekel, Profiles of dark haloes. Evolution, scatter, and environment , Mon. Not. Roy. Astron. Soc. 321 (2001) 559-575, [ astro-ph/9908159 ].
- [86] R. H. Wechsler, J. S. Bullock, J. R. Primack, A. V. Kravtsov, and A. Dekel, Concentrations of dark halos from their assembly histories , Astrophys. J. 568 (2002) 52-70, [ astro-ph/0108151 ].
- [87] C. A. Correa, J. S. B. Wyithe, J. Schaye, and A. R. Duffy, The accretion history of dark matter haloes - III. A physical model for the concentration-mass relation , Mon. Not. Roy. Astron. Soc. 452 (2015), no. 2 1217-1232, [ arXiv:1502.00391 ].
- [88] J. D. Lewin and P. F. Smith, Review of mathematics, numerical factors, and corrections for dark matter experiments based on elastic nuclear recoil , Astropart. Phys. 6 (1996) 87-112.
- [89] D. G. Cerdeno and A. M. Green, Direct detection of WIMPs , arXiv:1002.1912 .
- [90] XENON , E. Aprile et. al., Dark Matter Search Results from a One Ton-Year Exposure of XENON1T , Phys. Rev. Lett. 121 (2018), no. 11 111302, [ arXiv:1805.12562 ].
- [91] CRESST , G. Angloher et. al., Results on light dark matter particles with a low-threshold CRESST-II detector , Eur. Phys. J. C76 (2016), no. 1 25, [ arXiv:1509.01515 ].
- [92] The GAMBIT Dark Matter Workgroup , T. Bringmann et. al., DarkBit: A GAMBIT module for computing dark matter observables and likelihoods , Eur. Phys. J. C77 (2017), no. 12 831, [ arXiv:1705.07920 ].
- [93] GAMBIT , P. Athron et. al., Global analyses of Higgs portal singlet dark matter models using GAMBIT , Eur. Phys. J. C79 (2019), no. 1 38, [ arXiv:1808.10465 ].
- [94] D. Hooper and S. J. Witte, Gamma Rays From Dark Matter Subhalos Revisited: Refining the Predictions and Constraints , JCAP 1704 (2017), no. 04 018, [ arXiv:1610.07587 ].
- [95] S. Kazantzidis, L. Mayer, C. Mastropietro, J. Diemand, J. Stadel, and B. Moore, Density profiles of cold dark matter substructure: Implications for the missing satellites problem , Astrophys. J. 608 (2004) 663-3679, [ astro-ph/0312194 ].
- [96] A. Gould, Resonant Enhancements in WIMP Capture by the Earth , Astrophys. J. 321 (1987) 571.
- [97] A. Green and P. Tanedo, DarkCapPy: Dark Matter Capture and Annihilation , Comput. Phys. Commun. 242 (2019) 120-131, [ arXiv:1808.03700 ].
- [98] J. Silk, K. A. Olive, and M. Srednicki, The Photino, the Sun and High-Energy Neutrinos , Phys. Rev. Lett. 55 (1985) 257-259. [,283(1985)].
- [99] A. Gould, Cosmological density of WIMPs from solar and terrestrial annihilations , Astrophys. J. 388 (1992) 338-344.

| [100]   | A. M. Serenelli, S. Basu, J. W. Ferguson, and M. Asplund, New Solar Composition: The Problem with Solar Models Revisited , Astrophys. J. 705 (Nov., 2009) L123-L127, [ arXiv:0909.2668 ].                  |
|---------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [101]   | R. H. Helm, Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even Nuclei , Phys. Rev. 104 (1956) 1466-1475.                                                                        |
| [102]   | R. Catena and B. Schwabe, Form factors for dark matter capture by the Sun in effective theories , JCAP 1504 (2015), no. 04 042, [ arXiv:1501.03729 ].                                                      |
| [103]   | G. Busoni, A. De Simone, P. Scott, and A. C. Vincent, Evaporation and scattering of momentum- and velocity-dependent dark matter in the Sun , JCAP 1710 (2017), no. 10 037, [ arXiv:1703.07784 ].          |
| [104]   | K. Griest and D. Seckel, Cosmic Asymmetry, Neutrinos and the Sun , Nucl. Phys. B283 (1987) 681-705. [Erratum: Nucl. Phys.B296,1034(1988)].                                                                 |
| [105]   | A. Ibarra, M. Totzauer, and S. Wild, High-energy neutrino signals from the Sun in dark matter scenarios with internal bremsstrahlung , JCAP 1312 (2013) 043, [ arXiv:1311.1418 ].                          |
| [106]   | K. Choi, C. Rott, and Y. Itow, Impact of the dark matter velocity distribution on capture rates in the Sun , JCAP 1405 (2014) 049, [ arXiv:1312.0273 ].                                                    |
| [107]   | A. Nuñez-Castiñeyra, E. Nezri, and V. Bertin, Dark matter capture by the Sun: revisiting velocity distribution uncertainties , arXiv:1906.11674 .                                                          |
| [108]   | A. H. G. Peter, WIMP astronomy and particle physics with liquid-noble and cryogenic direct-detection experiments , Phys. Rev. D83 (2011) 125029, [ arXiv:1103.5145 ].                                      |
| [109]   | B. J. Kavanagh and A. M. Green, Model independent determination of the dark matter mass from direct detection experiments , Phys. Rev. Lett. 111 (2013), no. 3 031302, [ arXiv:1303.6868 ].                |
| [110]   | B. J. Kavanagh, M. Fornasa, and A. M. Green, Probing WIMP particle physics and astrophysics with direct detection and neutrino telescope data , Phys. Rev. D91 (2015), no. 10 103533, [ arXiv:1410.8051 ]. |
| [111]   | B. Feldstein and F. Kahlhoefer, A new halo-independent approach to dark matter direct detection analysis , JCAP 1408 (2014) 065, [ arXiv:1403.4606 ].                                                      |
| [112]   | F. Ferrer, A. Ibarra, and S. Wild, A novel approach to derive halo-independent limits on dark matter properties , JCAP 1509 (2015), no. 09 052, [ arXiv:1506.03386 ].                                      |
| [113]   | A. Ibarra, B. J. Kavanagh, and A. Rappelt, Bracketing the impact of astrophysical uncertainties on local dark matter searches , JCAP 1812 (2018), no. 12 018, [ arXiv:1806.08714 ].                        |
| [114]   | S. Garrison-Kimmel, M. Boylan-Kolchin, J. Bullock, and K. Lee, ELVIS: Exploring the Local Volume in Simulations , Mon. Not. Roy. Astron. Soc. 438 (2014), no. 3 2578-2596, [ arXiv:1310.6746 ].            |

- [115] J. Diemand, M. Kuhlen, P. Madau, M. Zemp, B. Moore, D. Potter, and J. Stadel, Clumps and streams in the local dark matter distribution , Nature 454 (2008) 735-738, [ arXiv:0805.1244 ].
- [116] M. R. Lovell, C. S. Frenk, V. R. Eke, A. Jenkins, L. Gao, and T. Theuns, The properties of warm dark matter haloes , Mon. Not. Roy. Astron. Soc. 439 (2014) 300-317, [ arXiv:1308.1399 ].
- [117] V. H. Robles, T. Kelley, J. S. Bullock, and M. Kaplinghat, The Milky Way's Halo and Subhalos in Self-Interacting Dark Matter , arXiv:1903.01469 .
- [118] A. Moliné, J. A. Schewtschenko, M. A. Sánchez-Conde, A. Aguirre-Santaella, S. A. Cora, and M. Abadi, Properties of subhalos in the interacting dark matter scenario , arXiv:1907.12531 .
- [119] P. Tinyakov, I. Tkachev, and K. Zioutas, Tidal streams from axion miniclusters and direct axion searches , JCAP 1601 (2016), no. 01 035, [ arXiv:1512.02884 ].
Alejandro Ibarra Bradley J. Kavanagh Andreas Rappelt
Impact of substructure on local dark matter searches
Dark matter substructure can contribute significantly to local dark matter
searches and may provide a large uncertainty in the interpretation of those
experiments. For direct detection experiments, sub-halos give rise to an
additional dark matter component on top of the smooth dark matter distribution
of the host halo. In the case of dark matter capture in the Sun, sub-halo
encounters temporarily increase the number of captured particles. Even if the
encounter happened in the past, the number of dark matter particles captured by
the Sun can still be enhanced today compared to expectations from the host halo
as those enhancements decay over time. Using results from an analytical model
of the sub-halo population of a Milky Way-like galaxy, valid for sub-halo
masses between $10^{-5}\,M_\odot$ and $10^{11}\,M_\odot$, we assess the impact
of sub-halos on direct dark matter searches in a probabilistic way. We find
that the impact on direct detection can be sizable, with a probability of $\sim
10^{-3}$ to find an $\mathcal{O}(1)$ enhancement of the recoil rate. In the
case of the capture rate in the Sun, we find that $\mathcal{O}(1)$ enhancements
are very unlikely, with probability $\lesssim 10^{-5}$, and are even impossible
for some dark matter masses.