The detection of sub-solar mass dark matter halos

Dark matter halos of sub-solar mass are the first bound objects to form in cold dark matter theories. In this article, I discuss the present understanding of"microhalos'', their role in structure formation, and the implications of their potential presence, in the interpretation of dark matter experiments.


Introduction
The term "microhalo" loosely refers to a generic dark matter, gravitationally-bound object, whose mass is typically less than that of the Sun. Objects of sub-solar mass scales are predicted to be the first objects formed in the high-redshift Universe in theories where the dark matter particle is "cold" (such as ΛCDM). In this short article I summarize the present understanding of microhalos. In Sec. 2 I review the theoretical motivation behind the presence of microhalos, in Sec. 3 I present the results of numerical simulations that studied their formation. Sec. 4 is an overview of the present state of understanding the survival of these objects, while in Sec. 5 I discuss the procedure of how to characterize their physical properties. I present their connection to direct and indirect detection experiments in Secs. 6 & 7, and I conclude in Sec. 8. Throughout the paper, the assumed cosmology is flat ΛCDM, with Ω m h 2 = 0.1358, Ω b h 2 = 0.02267, h = 0.705, σ 8 = 0.812, and a spectral index of n s = 0.96 [1].

Kinetic decoupling and microhalo scales
In the standard cold dark matter cosmological model, the energy density in the Universe is balanced by approximately 4% baryonic matter, 23% cold dark matter and 73% dark energy. Cold dark matter refers to gravitationally interacting matter that is hypothesized to be in the form of a yet-to-be-discovered particle. From the particle physics aspect, new particles arise in essentially all extensions to the standard model of particle physics, and usually the lowest mass particle is stable due to a new symmetry. Dark matter particle candidates include (among many others), neutralinos, the lightest Kaluza-Klein particle, axions, and sterile neutrinos. Even though there is no a priori physical reason as to one candidate being more favorable than another, neutralinos are attractive and well-studied because they are experimentally accessible at present. Neutralinos arise in supersymmetric extensions to the standard model of particle physics [2,3], and they are part of a generic class of dark matter candidates, called WIMPs (Weakly Interacting Massive Particles).
In the early Universe, conditions are such that a WIMP is in chemical, thermal and kinetic equilibrium. As the rates of these interactions are diluted due to the expansion of the Universe, the WIMP falls out of equilibrium and decouples. The decoupling temperature is determined by the scattering cross section that is responsible for the equilibrium condition the particle decouples from, while the velocity distribution function of the particle is set by its mass. The temperature at which a WIMP is kinetically decoupled, i.e., when momentum-changing interactions cease to be effective, is called the kinetic decoupling temperature. For supersymmetric dark matter, such as the neutralino, the kinetic decoupling temperature is T d ∼ [10−1000] MeV [4,5,6]. The free-streaming of particles after kinematic decoupling tends to smooth out fluctuations on small scales [7,4,8,9]. This physical effect leads to a cutoff in the dark matter power spectrum.  [8,9].
The cutoff of the power spectrum includes an imprint of the acoustic oscillations of the cosmic radiation fluid [10]. In general the transfer function of the dark matter density perturbation amplitude is obtained by solving the Boltzman equation, but under certain approximations, it is possible to be derived analytically [8,9]. Fig. 1 shows an example of the analytic calculation of the cutoff of the power spectrum, derived from [8,9]. What is shown is the damping factor D(k) that arises for a supersymmetric WIMP dark matter particle with a decoupling scale of k d ∼ 40 pc −1 , and a free-streaming scale of k fs ∼ 0.94 pc −1 (corresponding to a kinetic decoupling temperature of T d ∼ 20 MeV). In general, the cutoff scale is related to the kinetic decoupling temperature as M cut ≈ 10 −4 M ⊙ (T d /10 MeV) −3 [10]. Recently, [11] presented the posterior of the cutoff scale from a Bayesian analysis of the presently-allowed supersymmetric parameter space.

Formation of microhalos
The first attempt to numerically simulate the formation of microhalos was done in 2005 by the group at the University of Zürich [12]. In that work, a small region of the Universe was simulated using a multi-scale technique [13] in order to achieve the required spatial resolution needed to resolve objects below the cutoff-scale of the power spectrum. The input transfer function was taken from [8]; corresponding to a supersymmetric dark matter candidate of mass ∼ 100GeV. This particular choice leads to an exponential cut-off at ∼ 10 −6 M ⊙ . The simulation was evolved from a redshift of z = 350 to a redshift of z = 26, when the simulated region itself approached non-linearity. In the small, high-resolution volume of roughly average density, the first non-linear structures with mass near the cutoff scale of ∼ 10 −6 M ⊙ are observed to form by redshift z ≈ 60.
Their physical size is few ×10 −2 pc, and they are described by a single power-law density profile. Even though halos of mass in the range [10 −6 −10 0 ]M ⊙ are formed roughly at the same time, they all share the same concentrations at z = 26. As expected, these objects contained no substructure, and no objects with smaller mass were formed, confirming the input shape of the power spectrum. Even though the statistics in this first simulation were low, the deduced mass function suggested a power law and normalization consistent with the halo mass function derived from numerous orders of magnitude larger scales, namely, dN/d ln M ∼ M −1 . This finding is an important piece of information when it comes to studies aimed at the long-term survival of microhalos in the potential well of a present-day galaxy, such as the Milky Way.
In a subsequent paper [14], the numerical simulation of the formation and evolution of microhalos was studied by using 64 million dark matter particles within a comoving volume box of 3 kpc on the side (implying a particle mass of 9.8 × 10 −10 M ⊙ ). The initial redshift was z = 456, and the initial conditions were similar to the simulation discussed in [12] (see previous paragraph). By the end redshift of the simulation, z ≈ 75, the volume contained almost 2000 virialized dark matter sub-solar mass halos. The most massive halo is ∼ 10 −2 M ⊙ , and corresponds to a 3.5σ fluctuation in the density field. For a comparison of the evolution and substructure content of this simulation, the authors performed a separate simulation, of similar dynamic range, of a cluster-size halo at z = 0. Comparison between the substructure mass functions in the two simulations showed confirmation of the earlier result found in [12], namely a power law behavior that is self-similar over more than 10 orders of magnitude in subhalo mass.
It should be emphasized that the simulation presented in [12] was performed in a region of average density, thus probing the formation of "field" sub-solar mass dark matter halos. On the other hand, the high-resolution simulation presented in [14] was centered around the highest density peak of the initial 3 kpc volume. This allowed the study of the assembly and structure of an object with mass higher than the cutoff scale of the power spectrum, while resolving scales down to the cutoff scale. In addition, it provided the first estimate of the survival rate of dark matter microhalos, a very important result which relies on the non-linear evolution of structure formation.

Challenges in numerical simulations
The numerical simulation of the formation and survival of microhalos is a very difficult task. The difficulties arise from two main sources. First, by the fact that the power spectrum P (k) of sub-solar mass scales asymptotically approaches the k −3 behavior, thus leading to the formation of multiple scales roughly at similar epochs, and second, by uncertainties in the actual implementation of numerical techniques in numerical simulations [15].
The cutoff in the power spectrum that enters the initial conditions in the numerical simulations does not fully suppress the formation of structure at scales smaller than the cutoff scale. Even though the fraction of substructure in objects below the cutoff scale is small, their mere presence deserves further study in order to determine whether this effect is an artifact of the simulation techniques employed, or whether it is a physical effect. It should be noted that such effects have also been found in simulations of Warm Dark Matter (WDM), where the dark matter particle possesses a non-zero thermal velocity [16,17,18,19]. A possible numerical source that leads to this artificial effect is the implementation of grid initial conditions [20]. If not implemented at an early enough starting redshift, the initial Zel'dovich displacement can be a large fraction of the inter-particle separation [21]. This can lead to artificial fragmentation along filaments, however the effect has also been observed when the use of glass initial conditions [22] is implemented. The simulations of [14] also show that the formation of sub-cutoff scale structure is not limited to a preferred separation along a grid axis, suggesting that it may perhaps be a real effect.
Another source that can introduce bias in the measurement of a mass function is the starting redshift of simulations. As shown in [21], the requirement that the Fourier modes within the simulation box are all linear is very naive. The mass function cannot be robustly determined at redshifts which are close to the redshift of first crossing (the redshift where halos form). It has also been pointed out that even perhaps the choice of a halo finder in the simulation can also lead to an uncertainty in the measurement of the mass function [14,15], due to the fact that the density contrast in sub-solar mass halos is low due to similar formation times. Therefore, caution must be taken in the interpretation of these results, and a thorough convergence study on the mass function of sub-solar mass halos is of paramount importance, and in great need.

Survival of microhalos in the Milky Way halo
The survival (and thus existence) of microhalos in the present-day Milky Way halo has been a subject of debate. The approach typically taken in these studies has progressed in two fronts: simple analytical estimates, and through numerical simulations. In principle the two should agree, however given the vast orders of magnitude in dynamic range and highly complex non-linear effects, it has proven to be a difficult task.
The question of surviving microhalos was first addressed analytically in [23]. A detailed analytical study of the effects of tidal destruction on microhalos due to rapid, early-time hierarchical merger phase, showed that only [0.1−0.5]% of formed microhalos would survive to the present epoch [23]. The initial claim that a large number of microhalos could be present in the solar neighborhood [12] was challenged only days later [24]. The argument presented was that strong impulses by individual stars in the disc can tidally disrupt microhalos, leading to a large number of tidal streams. Within days, a response [25] showed that the application of the impulse approximation requires microhalos to experience multiple stellar interactions and as a result the survival timescale is many Hubble times.
One of the key advantages of numerical simulations in the study of the formation and abundance of microhalos is the ability to trace the history of each bound structure.
In doing so, the work of [14] found that substructure formed at early times is highly prone to tidal disruption due to the high central density of the host halo. In addition, most of the remaining microhalos that survived to the end-redshift of the simulation run are significantly tidally stripped. As expected, the amount of mass-loss experienced due to tidal effects depends on the history of the subhalo and in particular whether it experienced passages near the center of the host halo. However, most importantly, the study of [14] demonstrated that even though the density contrast of sub-solar mass scales is small, microhalos of mass 10 −6 M ⊙ are able to survive in a host of mass at least 10 4 times more massive, thus addressing numerically the question of survival in the initial merging process.
A separate study [26] looked at the energy input to microhalos due to stellar interactions by using analytic calculations, and then tested their results in a numerical simulation [27,28]. The main result is that at large impact parameters, the energy input to microhalos is independent of mass. They show that the survival timescale of microhalos in the Milky Way is of order the age of the Milky Way, but decreases significantly as the microhalo mass increases. In an independent test of the impulse approximation [29], it was shown that multiple stellar encounters can remove most of the mass from microhalos, however the very inner cores can survive to the present time. Nevertheless, the normalization and shape of the density distribution of dark matter in microhalos will be decreased, with significant implications to direct and indirect detection experiments (see following sections). It is important to realize that it is currently not feasible to fully address the issue of survival, mass and spatial distribution of microhalos, due to the enormous dynamic range required and the limited computing power that technologically exists. Studies such as the ones mentioned here are invaluable in shedding light on what might happen in ideal situations, and therefore should be further explored.

Characterizing the physical properties of microhalos
In order to assess the detectability of microhalos in any scheme (e.g. direct or indirect detection experiments), it is first necessary to determine the distribution of dark matter within the virialized region that comprises the microhalo. Typically, this can be done by assuming a profile, and a normalization of that profile.
The numerical work of [12,14] finds that microhalos at the redshift at which they were studied in the simulation are described by an NFW profile [30]. The form of this profile is ρ(r) = ρ s /r(1+r) 2 , where ρ s is a density normalization, called the characteristic density, r s is the scale radius, andr = r/r s . The two parameters needed to specify the profile are the characteristic density, and the scale radius. The numerical simulation findings suggest that r s is approximately equal to the virial radius R vir of the microhalo, suggesting an NFW concentration parameter (c v = R vir /r s ) at that redshift which is of order one (more specifically c v ≈ [2 − 4.5] [31]). For field microhalos, the concentration parameter grows as (1 + z) −1 , indicating that "field" subhalos have concentrations in the range of [50 − 100] for a formation redshift around z ∼ 20 (see below). Obviously, the profile of a present-day microhalo in the Milky Way has most likely evolved from this initial distribution of dark matter, to a distribution which is the result of tidal interactions throughout the course of structure formation. As I alluded to in Secs 3 & 4, these processes are very difficult to accurately model. Nevertheless, it is a reasonable assumption to consider the initial profile at the redshift of formation as a representation of a maximal distribution, and keep in mind that any experimental ramifications are going to actually be upper limits and not absolute values.
In order to normalize the dark matter density profile of microhalos we need to first determine their formation time, which can be obtained by knowledge of the CDM power spectrum, and the cosmological growth factor.
The normalization of a the profile of a dark matter halo can be obtained by assuming that mean density of the halo is ∆ vir times the mean matter density of Universe ρ M , with ∆ vir being the virial overdensity (see e.g. [32]). The mass of the halo is then M = 4πρ M ∆ vir (z)R 3 vir /3. By using this definition of the mass of a dark matter halo, the normalization of the NFW profile is then simply obtained via x). The rarity of a particular mass scale at a particular redshift can be obtained from the power spectrum of fluctuations. In general, if we define the power in each logarithmic interval in k as ∆ 2 (k) ∼ k 3 P (k), where P (k) is the CDM power spectrum, then the mean square fluctuations on a scale M is given by σ 2 (M) = ∆ 2 (k)|W (k; M)| 2 d ln k. Here, W (k; M) is the Fourier transform of a realspace spherical top-hat window that contains mass M. The left panel in Figure 2 shows the root mean square fluctuations on microhalo scales, derived by assuming the WMAP5 cosmological parameters [1]. Note that for illustration purposes, the primordial power spectrum in this example does not contain the cutoff shown in Fig. 1. Knowledge of the underlying cosmological model (i.e., the values of the cosmological parameters), allows us to compute the growth function of fluctuations D(z). Then, the redshift of formation of a mass scale M is defined as the redshift z where σ(M)D(z) = δ c , where δ c = 1.686 is the value of the characteristic density for collapse (derived from linear theory). The rarity of a particular peak is typically defined to be in units of the standard deviation of the smoothed density distribution, so that ν = δ c /σ(M)D(z).
Note that the rarity of a particular fluctuation in the density field is directly related to the abundance of that scale in the present Milky Way halo -a key quantity in any dark matter experimental effort. For example, a smaller fraction of the total Milky Way mass originates from 3σ peaks in the density field as those are rarer than e.g. 1σ peaks. However, the 3σ peaks are formed earlier, thus they will be denser. As such, they are not as susceptiple to tidal disruption as the 1σ peaks. The abundance of microhalos in the present epoch will be determined by a balance of these competing effects and a full description of the microhalo population of the Milky Way must include all of these factors. A first step in that direction was the study of [34], where it was found that high-ν progenitors are found predominantly found near the centers of dark matter halos today. For example, according to [34], the median progenitor at the solar neighborhood of the Milky Way corresponds to a 2σ peak.
An additional piece of information which is of importance when it comes to dark matter detection experiments, is the tidal radii of survived microhalos in the Milky Way. A reasonable approach would be to assume that the truncation radius at the present epoch can be estimated by finding the radius of the microhalo where the density is equal to the density at the solar radius [35] . This assumption results in the truncation radius being equal to the virial radius for the microhalos halos considered here, reflecting the dense state of the Universe at the time they were formed. Clearly, the assumption here is that the halo profile normalization is intact from the time it was formed until the present epoch. Again, this points out at the lack of knowledge of the mass function and history (thus properties) of microhalos in the solar neighborhood. But given the high initial densities reflected in microhalos, and the steep power law of the density profile, it seems likely that the tidal radius will not be very different from their virial radius set at the time of collapse.
A more general note regarding the physical description of the dark matter profiles of microhalos is with regards to the concentration parameter. It is dangerous to assign concentration parameters that would correspond to present-epoch "field halos" (i.e., a microhalo formed at high redshift and evolved uninterrupted to the present day). Field halos grow hierarchically by the accretion of smaller halos, and by smooth accretion of dark matter. It is not clear exactly how this process occurs at very high redshifts on microhalo scales, as the flatness of the power spectrum implies the simultaneous collapse of many different scales above the cutoff in the power spectrum. The maximum value of the concentration parameter for a field microhalo can be obtained by extrapolating the concentration-mass relationship deduced from numerical simulations [33,36] (see Fig. 2). Even though this is a daring extrapolation over numerous orders of magnitude, it can be considered as an upper limit for field microhalos, as the concentration-mass relationship stems from the shape of the power spectrum. The power spectrum flattens out at large k, and so is c v with respect to mass [33] (see Fig. 2). The danger from assuming concentration parameters of field halos in microhalo studies near the solar neighborhood arises due to the fact that the density profile normalization is rather sensitive to the value of the concentration parameter ρ s ∼ c 3 v , therefore assuming concentration parameters of order 100, can result in a very large overestimate of the normalization of the profile.
To summarize, given the uncertainty in the evolution of the profile of microhalos, the only reasonable approach is to determine their physical properties at the time they were formed, and consider those as upper limits to what they could be at the present epoch. Figure 3 & 4 shows the formation redshift, characteristic density and virial radius of sub-solar mass dark matter halos. It is of vital importance to emphasize that any treatment of annihilation signals from microhalos (see e.g. γ-rays in the next section) must be consistent with a description of the physical properties of microhalos, which are set by the underlying structure formation paradigm. Arbitrary choices of concentration parameters (which affect the normalization of the profile), or inconsistent choices of formation redshift will only result in an estimate for a particular toy microhalo model. Such an approach cannot assess the structure formation relevance, and/or distribution function of expected signals and must treated with caution.

Detection of microhalos using dark matter annihilation products
It might be possible to detect the presence of microhalos in the Milky Way halo by searching for the annihilation products of the WIMP dark matter particle, and specifically for γ-rays. Under the supersymmetric WIMP scenario, the dark matter particle couples to quarks and leptons at the tree level, and to photons and the Z 0 gauge boson via one-loop diagrams. The search of γ-rays from the annihilation of the dark matter particle is attractive for two reasons: 1) a swath of γ-ray experiments are operational (e.g. Fermi/GLAST [37], VERITAS [38], H.E.S.S. [39], MAGIC [40], CANGAROO [41]), and 2) the shape of the emitted spectrum can be estimated (given a supersymmetric dark matter candidate).
In general, the luminosity of a subhalo as a function of energy is given by Fig. 4), and P (E) = (dN/dE) σv /M 2 χ . Here, dN/dE is the energy spectrum of the annihilation products under study (e.g. photons, positrons, neutrinos, etc.), σv is the annihilation cross section, and M χ is the mass of the dark matter particle. Maximal values of P (E)dE in supersymmetric theories with a valid dark matter candidate (e.g. the neutralino) yield a value of ≤ 10 −28 cm 3 GeV −2 for a threshold energy of ∼ 1 GeV.

Microhalos as γ-ray sources
In this subjection I discuss the possible detection of microhalos with γ-rays, though the formalism can be well-extended to other annihilation products, such as positrons and neutrinos. In order for a sub-solar mass halo at a distance d and flux on Earth Φ = L/4πd 2 to be detected above a detector threshold Φ 0 , it must be located at a distance d < Φ/Φ 0 . Such an example is shown in Fig. 5, for a supersymmetric WIMP dark matter candidate with P ≤ 10 −28 cm 3 GeV −2 , and a detector threshold of Φ 0 = 10 −9 cm −2 s −1 . It should be noted that for subhalo mass functions that scale as dN/d ln MdV ∼ M −1 , the typical visibility distance between microhalos of mass M scales as r ∼ M 1/3 . However, as the luminosity of a halo goes as L ∼ ρ 2 s r 3 s ∼ M, the "visibility" distance of microhalos is proportional to d max ∼ M −1/2 , i.e., it is "easier" to detect more massive halos at larger distance [42,43], than small halos nearby. This scalings will be altered if the dark matter particle's annihilation cross section is Sommerfeld enhanced [44,45,46,47,48]. The increased annihilation rate due to low microhalo internal velocities leads to the possibility that small nearby microhalos will be more luminous than the canonical case where the annihilation cross section is independent of velocity.
One approach to investigate the observability of a population of microhalos in the Milky Way is to Monte Carlo the distribution and properties of microhalos as followed in [49]. The key unknown in this approach (or any approach for that matter) is the abundance of microhalos as a function of their mass, as well as their central densities (as there is a a distribution of mass scales that collapse in each redshift interval -see Sec. 5). The results of [49] show that for a fixed concentration microhalos (i.e., isolated density peaks), it is highly unlikely that there should be any high signal-to-nose single-source detection of microhalos with the all-sky Fermi Gamma-ray Space Telescope (FGST, formerly GLAST).
Another approach was followed in [50]. In that study, the uncertainty that stems from the lack of knowledge of the survival rate of microhalos was parametrized as a fixed contribution to the local dark matter density, so that a particular fluctuation peak (i.e., collapsed microhalo at a particular redshift) results in a certain abundance. Again, the main caveat here is the unknown survival rate. The result of [50] suggested that a large number of microhalos could potentially be visible, albeit the normalization of the number was directly deduced from the unknown density parametrization.
One unique feature of the potential presence of survived microhalos in the Milky Way is the probability that they can be present in very short distances (sub-parsec) from the Solar system. Thus, it may be possible to observe a proper motion of the γ-ray signal [25,50,51]. The two crucial conditions for the observation of the proper motion of a subhalo is that it is bright enough to be detected at some high signal-to-noise level, and it must be close enough so that the proper motion exhibited over the life-time of an experiment is above the threshold angular resolution of the detector. The details of this potential observation were studied in [50] under the assumption that a certain fraction of the Milky Way mass at the solar radius is in microhalos which were all formed in a particular epoch.
Any potential detection of proper motion implies the presence of a large number of "unresolved" microhalos at distances far beyond the flux-limited distance set by the sensitivity of a particular detector. The presence of a large number of microhalos below individual detection threshold would lead to a background radiation that would correlate with the angular and radial distribution of microhalos in the Milky Way. This was the subject of study in [52], where it was shown that the measured value of the γ-ray background measured by EGRET [53] is already placing stringent limits on the probability of proper motion detection (due to the implied low number density of microhalos in the Milky Way). Addressing the same issue, [54] showed that regardless of the assumptions of the inner density profile of microhalos, the EGRET measurement of the γ-ray flux from the Galactic center [55] (see [56], but also [57]) is strongly suggesting that the probability of detecting microhalos with measurable proper motion is negligible. Nevertheless, given the simplicity of this measurement, a search for the proper motion of γ-ray-only sources should be performed due to the valuable information contained in such a potential detection.

Microhalos and the local Milky Way dark matter distribution
The presence of substructure in dark matter halos leads to interesting consequences when it comes to direct detection experiments. In general, the rate of events in a direct detection experiment is proportional to Γ ∼ ρ ⊙ , where ρ ⊙ is the dark matter density in the solar neighborhood. Typically, the value of ρ ⊙ is obtained from dynamical measurements of the structure of the Milky Way (e.g. [35]), with a canonical value of ρ ⊙ = 0.4 GeV cm −3 . This value is obtained from averaging regions which are of order ∼kpc, much larger than the regions probed in the duration of a direct detection experiment (sub-parsec scales), and it is susceptible to uncertainties that are introduced due to the unknown shape of the Milky Way halo (see e.g. [58,59]). If there are fluctuations on sub-parsec scales (as what one might expect from the presence of microhalos), then direct detection experiments are prone to structure formation uncertainties. In this section I give an overview of the experimental implications of Figure 6. The local dark matter density probability distribution function [60]. The peak of the distribution defines the dark matter density of the smooth distribution of dark matter, while the power low behaviour at higher densities arises from the presence of dark matter subhalos. This function is derived under the assumption that the substructure mass fraction in the Milky Way of all substructure in the mass range fluctuations in the density field at the solar radius.

Microhalos and the mean dark matter density
It is expected that if a certain amount of mass in a dark matter halo is distributed in high-density regions, then the density of the smoothly-distributed dark matter will be lower than the case where all matter is distributed smoothly. The presence of sub-solar mass microhalos in the Milky Way raises one important question: what is the value of the smooth component if the substructure mass function extends down to microhalo scales, and what is the probability of the solar system being in an overdense region at any given epoch. The answer to both of these questions is directly related to the local flux enhancement in the annihilation signal (the so-called "boost factor"'), as the probability of being in a density enhanced region is inversely proportional to the boost factor.
Addressing this question was the study presented in [60]. Under the assumption of hierarchical structure formation, it was shown that the derived local density probability distribution function (PDF) is positively skewed due to the presence of substructure, and that the peak of the distribution is always less than the canonical mean value of the local dark matter density in the solar neighborhood (see Fig. 6). This conclusion is derived by a simple analytic model based on the scale-invariant nature of the hierarchical structure formation, as well as by investigating the implications of the distribution and profiles of dark matter halos as seen in numerical simulations. Both approaches showed similar results. These results were recently confirmed in numerical simulations by analyzing a suite of halos taken from the Aquarius Project [61,62] and independendly, by an analysis of the Via Lactea II Milky Way simulation [59,63].
An interesting outcome of the approach of [60] is a derivation of the the local annihilation rate due to the granularity of the local dark matter halo. The local boost factor should not be confused with the global boost factor that is commonly referred in the litereature (see e.g. [64,11]). The least granular cases studied in [60] yield a boost factor which is roughly consistent with estimates based on the survival rate of microhalos [23] as well as direct simulation results, as shown in [63]. The local boost factor is ∼ [1 − 5], with a weak dependence on the cutoff scale of the subhalo mass function.
The introduced variance in the local dark matter density has significant implications to the interpretation of combined results from indirect and direct detection experiments, as the expected signal is proportional to different powers of the local particle density. The volume probed during a 3-year direct detection experiment is very small ( ∼ few × 10 −4 pc), while dynamical estimates of the local dark matter halo provide averages over much larger scales. The uncertainties implied by fluctuations in the local dark matter density due to the presence of substructure are manifested as uncertainties in the predicted rates in a dark matter detector [60]. If substructure in the form of microhalos is abundant in the solar neighborhood, then the local boost factor will be large (favoring indirect detection searches), however the probability of being in an overdense region will be low (disfavoring direct detection experiments). Fortuitously, present simulation results suggest that these effects (and similarly perhaps the uncertainties due to structure in phase space) are small, and a smooth dark matter halo is a safe assumption [61,62,59,63], though the consequences for indirect detection searches may be more pronounced (e.g., the explanation of the postitron excess in the PAMELA [65] and ATIC [66] data by the presence of a nearby subhalo [67]).
An interesting question would be to ask whether it is possible that a long-duration direct detection experiment will be able to map the local distribution of dark matter. If we assume the velocity of the Sun to be 220 km/s, then the distance sweept over ∼ 20 years, is x ∼ few × 10 −3 pc. An extrapolation of the subhalo mass function to microhalo masses (normalized so that 10% of the Milky Way halo is in objects with mass greater than 10 −5 M MW ) and assuming an NFW profile for the radial distribution of subhalos, results in a mean distance between microhalos which depends on mass as r ≈ 52M 0.3 pc. For microhalos of mass 10 −6 M ⊙ , this is r ≈ 0.8 pc. As x ≪ r , it seems unlikely that a direct detection experiment will measure the transition from a smooth component to an overdensity of dark matter. Coupled with the very small likelihood of being inside a microhalo at the present time (see [60,61,59]) it implies that the rate in direct detection experiments should be time independent (note that this statement does not include the annual modulation that should be present due to the relative motions of the Sun and the Earth in the Milky Way dark matter halo). Fig. 7 shows the expected time dependence of the signal in a long-duration direct detection experiment in the hypothetical scenario where the solar system passes through a microhalo. It is important to emphasize that this effects are not convolved with the likelihood of such an occurrence, they simply demonstrate the magnitude of any effect for each particular case shown.

Conclusion
Sub-solar mass dark matter halos are interesting, not only because they are linked to the nature of dark matter, but also because a detection of their presence would provide insight into structure formation at extremely early times. The detection of microhalos would first and foremost show that the dark matter particle is cold. In addition, it will place constraints on the value of the kinetic decoupling temperature and mass of the dark matter particle. From the cosmological perspective, any detection of sub-solar mass halos would provide insights into halo merging and growth at extremely high redshifts, a task unattainable by any other observational method. Their mere presence would imply that at least a certain fraction of them survived the rapid merger phase. Regardless of the difficulties presented by these highly non-linear structures, sub-solar mass dark matter halos are extremely interesting objects and certainly warrant further investigations.