Review Article  Open Access
Michael Kuhlen, "The Dark Matter Annihilation Signal from Dwarf Galaxies and Subhalos", Advances in Astronomy, vol. 2010, Article ID 162083, 15 pages, 2010. https://doi.org/10.1155/2010/162083
The Dark Matter Annihilation Signal from Dwarf Galaxies and Subhalos
Abstract
Dark Matter annihilation holds great potential for directly probing the clumpiness of the Galactic halo that is one of the key predictions of the Cold Dark Matter paradigm of hierarchical structure formation. Here we review the ray signal arising from dark matter annihilation in the centers of Galactic subhalos. We consider both known Galactic dwarf satellite galaxies and dark clumps without a stellar component as potential sources. Utilizing the Via Lactea II numerical simulation, we estimate fluxes for 18 Galactic dwarf spheroidals with published central densities. The most promising source is Segue 1, followed by Ursa Major II, Ursa Minor, Draco, and Carina. We show that if any of the known Galactic satellites can be detected, then at least ten times more subhalos should be visible, with a significant fraction of them being dark clumps.
1. Introduction
A decade has gone by since the emergence of the āMissing Satellite Problemāā [1, 2], which refers to the apparent discrepancy between the observed number of Milky Way satellite galaxies, 23 by latest count [3ā12], and the predicted number of dark matter (DM) subhalos that should be orbiting in the Milky Way's halo. The latest cosmological numerical simulations [13ā15] resolve close to 100,000 individual selfbound clumps of DM within the Galactic virial volumeāremnants of the hierarchical buildup of the Milky Way's DM halo. Several hundred of these may be large enough to host a dwarf galaxy. A consensus seems to be emerging that this discrepancy is not a shortcoming of the otherwise tremendously successful Cold Dark Matter (CDM) hypothesis [16, 17], but instead reflects the complicated baryonic physics that determines which subhalos are able to host a luminous stellar component and which are not [18ā24]. If this explanation is correct, then an immediate consequence is that the Milky Way dark halo should be filled with clumps on all scales down to the CDM freestreaming scale at to for cold thermal relics (e.g., supersymmetric CDM) [25ā28] or as low as in the case of nonthermal cold axions [29ā31]. At the moment there is little empirical evidence for or against this prediction, and this has motivated searches for new signals that could provide tests of this hypothesis, and ultimately help to constrain the nature of the DM particle.
One of the most promising such signals is DM annihilation [32]. In regions of sufficiently high density, for example in the centers of Galactic subhalos, the DM pair annihilation rate might become large enough to allow for a detection of neutrinos, energetic electrons and positrons, or ray photons, which are the byproducts of the annihilation process. This is one of the few ways in which the dark sector can be coupled to ordinary matter and radiation amenable to astronomical observation. Belying its commonly used name of āindirect detectionāā, DM annihilation is really the only way we can hope to directly probe the clumpiness of the Galactic DM distribution. One could argue that it is a more ādirectāā method than trying to constrain DM clumpiness from its effects on strong gravitational lensing (see Zackrisson & Riehm's contribution in this special edition), or from the kinematics of stars orbiting in DMdominated potentials [33], or from perturbations of cold stellar structures like globular cluster tidal streams [34ā37] or the heating of the Milky Way's stellar disk [38ā40], although all of these are also worthwhile approaches to take.
The only trouble with the DM annihilation signal is that so far there have been no undisputed claims of its detection. Recently there have been several reports of āanomalousāā features in the local cosmic ray flux: the PAMELA satellite reported an increasing positron fraction at energies between 10 and 100āGeV [41], where standard models of cosmic ray propagation predict a decreasing fraction; the ATIC [42] and PPBBETS [43] balloonborne experiments reported a surprisingly large total electron and positron flux at 500āGeV, although recent Fermi [44] and H.E.S.S. data [45] appear to be inconsistent with it. Either of these cosmic ray anomalies might be the long sought after signature of local DM annihilation. However, since the currently available data can equally well be explained by conventional astrophysical sources (e.g., nearby pulsars or supernova remnants), they hardly provide incontrovertible evidence for DM annihilation. The next few years hold great potential for progress, since the recently launched Fermi Gammaray Space Telescope will conduct a blind survey of the ray sky at unprecedented sensitivity, energy extent, and angular resolution. At the same time, Atmospheric Cerenkov Telescopes, such as H.E.S.S., VERITAS, MAGIC, and STACEE, are greatly increasing their sensitivity, and have only recently begun to search for a DM annihilation signal from the centers of nearby dwarf satellite galaxies [46ā50].
The purpose of this paper is to provide an overview of the potential DM annihilation signal from individual Galactic DM subhalos, either as dwarf satellite galaxies or as dark clumps. It does not cover a number of very interesting and closely related topics, which are actively being researched and deserve to be examined in equal detail. These include the diffuse flux from Galactic substructure and its anisotropies (e.g., [51ā53]), the relative strength of the signal from individual subhalos compared with that from the Galactic Center or an annulus around it [54, 55], the effect of a nearby DM subhalo on the amplitude and spectrum of the local high energy electron and positron flux [56, 57], and the role of the Sommerfeld enhancement [58] on the DM annihilation rate and its implications for substructure signals [59ā61].
This paper is organized as follows: we first review the basic physics of DM annihilation, briefly touching on the relic density calculation, the āWIMP miracleāā, DM particle candidates, and, in more detail, the sources of rays from DM annihilation. In the following section we review what numerical simulations have revealed about the basic properties of DM subhalos that are relevant for the annihilation signal. We go on to consider known Milky Way dwarf spheroidal galaxies as sources, using the Via Lactea II simulation to infer the most likely annihilation fluxes from published values of the dwarfs' central masses. Next we discuss the possibility of a DM annihilation signal from dark clumps, halos that have too low a mass to host a luminous stellar component. Lastly, we briefly discuss the role of the substructure boost factor for the detectability of individual DM subhalos.
2. Dark Matter Annihilation
If DM is made up of a socalled āthermal relicāā particle, (An alternative DM candidate is the axion, a nonthermal relic particle motivated as a solution to the strong CP problem [62]. Since it does not produce an annihilation signal today, we do not further consider it here.) its abundance today is set by its annihilation cross section in the early universe. The thermal relic abundance calculation relating today's abundance of DM to the properties of the DM particle (its mass and annihilation cross section) is straightforward and elegant, and has been described in pedagogical detail previously [63ā65]. We briefly summarize the story here.
At sufficiently early times, the DM particles are in thermal equilibrium with the rest of the universe. As long as they remain relativistic (), their creation and destruction rates are balanced, and hence their comoving abundance remains constant. Once the universe cools below the DM particle's restmass (), its equilibrium abundance is suppressed by a Boltzmann factor . If equilibrium had been maintained until today, the DM particles would have completely annihilated away. Instead the expansion of the universe comes to the rescue and causes the DM particles to fall out of equilibrium once the expansion rate (given by , the Hubble constant at cosmological scale factor ) exceeds the annihilation rate , that is, when DM particles can no longer find each other to annihilate. The comoving number density of DM particles is then fixed at a āfreezeoutāā temperature that turns out to be approximately , with only a weak additional logarithmic dependence on the mass and cross section of the DM particle. A back of the envelope calculation results in the following relation between , the relic mass density in units of the critical density of the universe , and , the thermally averaged velocityweighted annihilation cross section: Note that this relation is independent of . The WMAP satellite's measurement of the DM density is [66], implying a value of A more accurate determination of must rely on a numerical solution of the Boltzmann equation in an expanding universe, taking into account the full temperature dependence of the annihilation rate, including the possibilities of resonances and coannihilations into other, nearly degenerate dark sector particles (e.g., [67, 68]). It is a remarkable coincidence that this value of is close to what one expects for a cross section set by the weak interaction. This is the socalled āWIMP miracleāā, and it is the main motivation for considering weakly interacting massive particles (WIMPs) as prime DM candidates.
The Standard Model of particle physics actually provides one class of WIMPs, massive neutrinos. Although neutrinos thus constitute a form of DM, they cannot make up the bulk of it, since their small mass, āeV [66], implies a cosmological mass density of only . Since the standard neutrino is a hot dark matter particle, increasing its contribution to the total matter density would destroy small scale structure and violate constraints from galaxy clustering and the Lyman alpha forest. The attention thus turns to extensions of the Standard Model, which themselves are theoretically motivated by the hierarchy problem (the enormous disparity between the weak and Planck scales) and the quest for a unification of gravity and quantum mechanics. The most widely studied class of such models consists of supersymmetric extensions of the Standard Model. Additionally models with extra dimensions have received a lot of attention in recent years. Both of these approaches offer good DM particle candidates: the lightest supersymmetric particle (LSP), typically a neutralino in Rparity conserving supersymmetry, and the lightest KaluzaKlein particle (LKP), typically the particle, the first KaluzaKlein excitation of the hypercharge gauge boson, in Universal Extra Dimension models. For much more information, we recommend the comprehensive recent review of particle DM candidates by Bertone et al. [65].
The direct products of the annihilation of two DM particles are strongly model dependent. Typical channels include annihilations into charged leptons (), quarkantiquark pairs, and gauge and Higgs bosons (). In the end, however, the decay and hadronization of these annihilation products results in only three types of emissions: (i) high energy neutrinos and antineutrinos, (ii) relativistic electrons and protons and their antiparticles, and (iii) ray photons. Additional lower energy photons can result from the interaction of the relativistic electrons with magnetic fields (synchrotron radiation), with interstellar material (bremsstrahlung), and with the CMB and stellar radiation fields (inverse Compton scattering). In the following we will focus on the rays, since they are likely the strongest signal from Galactic DM substructure. rays are produced in DM annihilations in three ways (see accompanying Figure 1). (i)Since the DM particle is neutral, there is no direct coupling to photons. Nevertheless, copious amounts of secondary ray photons can be produced through the decay of neutral pions, , arising in the hadronization of the primary annihilation products. Since the DM particles are nonrelativistic, their annihilation results in a pair of monoenergetic particles with energy equal to , which fragment and decay into meson dominated ājetsāā. In this way a single DM annihilation event can produce several tens of ray photons. The result is a broad spectrum with a cutoff around .(ii)An important additional contribution at high energies () arises from the internal bremsstrahlung process [69], which may occur with any charged annihilation product. One can distinguish between final state radiation, in which the photon is radiated from an external leg, and virtual internal bremsstrahlung, arising from the exchange of a charged virtual particle. Note that neither of these processes requires an external electromagnetic field (hence the name internal bremsstrahlung). The resulting ray spectrum is peaked towards and exhibits a sharp cutoff. Although it is suppressed by one factor of the coupling compared to pion decays, it can produce a distinctive spectral feature at high energies. This could aide the confirmation of a DM annihilation nature of any source and might allow a direct determination of .(iii)Lastly, it is possible for DM particles to directly produce ray photons, but one has to go to looplevel to find contributing Feynman diagrams, and hence this flux is typically strongly suppressed by two factors of (although exceptions exist [70]). On the other hand, the resulting photons would be monochromatic, and a detection of such a line signal would provide strong evidence of a DM annihilation origin of any signal. While annihilations directly into two photons, , would produce a narrow line at , in some models it is also possible to annihilate into a photon and a boson, , and this process would result in a somewhat broadened (due to the mass of the ) line at ).
(a)
(b)
(c)
The relative importance of these three ray production channels and the resulting spectrum depend on the details of the DM particle model under consideration. For any given model, realistic ray spectra can be calculated using sophisticated and publicly available computer programs, such as the PYTHIA MonteCarlo event generator [71], which is also contained in the popular DarkSUSY package [72].
3. Dark Matter Substructure as Discrete Ī³Ray Sources
DM subhalos as individual discrete ray sources hold great potential for providing a āsmoking gunāā signature of DM annihilation [32, 54, 54, 73ā83]. Compared to diffuse ray annihilation signals, these discrete sources should be easier to distinguish from astrophysical backgrounds and foregrounds [84], since (a) typical astrophysical sources of high energy rays, such as pulsars and supernova remnants, are very rare in dwarf galaxies, owing to their predominantly old stellar populations, (b) the DM annihilation flux should be timeindependent, (c) angularly extended, and (d) not exhibit any (or only very weak) low energy emission due to the absence of strong magnetic fields or stellar radiation fields. (Inverse Compton radiation from scattering off the CMB, however, may be present.)
We can distinguish between DM subhalos hosting a Milky Way dwarf satellite galaxy and dark clumps that, for whatever reason, do not host a luminous stellar population, or one that is too faint to have been detected up to now. Before we go on to discuss the prospects of detecting a DM annihilation signal from these two classes of sources, we review the basic properties of DM subhalos common to both.
Numerical simulations have shown that pure DM (sub)halos have density profiles that are well described by a Navarro, Frenk & White (NFW) [85] profile over a wide range of masses [86, 87], The parameter indicates the radius at which the logarithmic slope , and . The very highest resolution simulations have recently provided some indications of a flattening of the density profile in the innermost regions [15, 88]. In this case a socalled Einasto profile may provide a better overall fit, Here the additional parameter governs how fast the profile rolls over, and has been found to have a value of in numerical simulations [88]. Note that the two density profiles actually do not differ very much until (cf. Figure 2(a)). Simulated DM halos are of course not perfectly spherically symmetric, and instead typically exhibit prolate or triaxial isodensity contours that become more elongated towards the center [89]. The degree of prolateness decreases with mass, and galactic subhalos have axis ratios of 0.7 [90].
(a)
(b)
The āvirialāā radius of a halo is defined as the radius enclosing a mean density equal to , where (from a spherical collapse calculation [91]) and is the mean density of the universe. The corresponding virial mass is the mass within , and a halo's concentration can then be defined as . While these quantities are well defined for isolated halos and commonly used in analytic models, they are somewhat less applicable to galactic subhalos, since the outer radius of a subhalo is set by tidal truncation, which depends on the subhalo's location within its host halo. Furthermore, in numerical simulations it is difficult to resolve in low mass subhalos. For this reason we prefer to work with , the maximum of the circular velocity curve and a proxy for a subhalo's mass, and , the radius at which occurs. These quantities are much more robustly determined for subhalos in numerical simulations than . Note that even can be affected by tidal interactions with the host halo, especially for subhalos close to the host halo center. For this reason we also sometimes consider , the largest value of that a subhalo ever acquired during its lifetime (i.e., before tidal stripping began to lower its ) and , the corresponding radius.
Since DM annihilation is a two body process, its rate is proportional to the square of the local density, and the annihilation āluminosityāā is given by the volume integral of , has dimensions of (mass) (length), and we express it in units of pc. In order to convert to a conventional luminosity, one must multiply by a particle physics term, where is the speed of light, the mass of the DM particle and the thermally averaged velocityweighted annihilation cross section discussed in the previous section. This is the total luminosity, but we are interested here only in the fraction emitted as rays. Furthermore, a given detector is only sensitive to rays above a threshold energy of and below a maximum energy of . In that case the effective ray luminosity is where is the spectrum of ray photons produced in a single annihilation event.
A comparison of the enclosed luminosity and mass profiles is shown in the bottom panel of Figure 2. Clearly, is much more centrally concentrated than : 90% of the total luminosity is produced within , compared with only 10% of the total mass. In terms of , the total luminosity of a halo is given by where is an numerical factor that depends on the shape of the density profile; for an NFW profile , and for an Einasto profile . In physical units, the total annihilation luminosity is for NFW (top) and Einasto (bottom). Note that even though the slope of the Einasto profile is shallower than NFW in the very center, the total luminosity exceeds that of an NFW halo with the same . This is due to the fact that the Einasto profile rolls over less rapidly than the NFW profile, and actually has slightly higher density than NFW between and a crossover point at .
4. Milky Way Dwarf Spheroidal Galaxies
There are several advantages of known dwarf satellite galaxies as DM annihilation sources: firstly, the kinematics of individual stars imply masstolight ratios of up to several hundred [92ā95], and hence there is an a priori expectation of high DM densities; secondly, since we know their location in the sky, it is possible to directly target them with sensitive atmospheric Cerenkov telescopes (ACT) such as H.E.S.S., VERITAS, MAGIC, and STACEE, whose small field of view makes blind searches impractical; lastly, our approximate knowledge of the distances to many dwarf satellites would allow a determination of the absolute annihilation rate, which may lead to a direct constraint on the annihilation cross section, if the DM particle mass can be independently measured (from the shape of the spectrum, e.g.).
Recent observational progress utilizing the Sloan Digital Sky Survey (SDSS) has more than doubled the number of known dwarf spheroidal (dSph) satellite galaxies orbiting the Milky Way [4ā12], raising the total from the 9 āclassicalāā ones to 23. Many of the newly discovered satellites are socalled āultrafaintāā dSph's, with luminosities as low as and only tens to hundreds of spectroscopically confirmed member stars. Simply accounting for the SDSS sky coverage (about 20%), the total number of luminous Milky Way satellites can be estimated to be at least 70. Taking into account the SDSS detection limits [96] and a radial distribution of DM subhalos motivated by numerical simulations, this estimate can grow to several hundreds of satellite galaxies in total [97, 98].
In order to assess the strength of the DM annihilation signal from these dSph's, it is necessary to have an estimate of the total dynamical mass, or at least , of the DM halo hosting the galaxies. Owing to the extreme faintness of these objects and their lack of a detectable gaseous component [99], it has been very difficult to obtain kinematic information that allows for such measurements. Progress has been made through spectroscopic observations of individual member stars, whose lineofsight velocity dispersions have confirmed that these objects are in fact strongly DM dominated [92ā95]. Such data best constrain the enclosed dynamical mass within the stellar extent, which on average is about 0.3ākpc for current data sets. A recent analysis has determined for 18 of the Milky Way dSph's, and found that, surprisingly, they all have to within a factor of two [33]. At present the physical meaning of this mass scale is unclear, but it may provide a clue to the nature of the baryonic processes that govern dwarf galaxy formation.
Stateoftheart cosmological numerical simulations of the formation of the DM halo of a Milky Way scale galaxy, such as those of the Via Lactea Project [13, 80] and the Aquarius Project [14], have now reached an adequate mass and force resolution to directly determine in their simulated subhalos. This makes it possible to infer the most likely values of for a Milky Way dSph of a given and Galactocentric distance , by identifying all simulated subhalos with comparable and and averaging over their . This analysis was performed for the 9 āclassicalāā dwarfs using the Via Lactea I simulation [22], and we extend it here to all 18 dwarfs published in [33] and with the more recent and higher resolution Via Lactea II (VL2) simulation.
For each of the Milky Way dSph's studied in [33] we created a list of all VL2 subhalos with numerically determined within the published error bars and distances within 40% of the measured distance of the given dSph. The distances were calculated for 100 random observer locations 8kpc from the host halo center. We then determined the median value and the 16th and 84th percentiles of and of the subhalos associated with each dSph. These values are given in Table 1. The median values of range from 5.0ākmās (Leo IV) to 19ākmās (Draco, Leo I), and of from 6.7ākmās (Leo IV) to 30ākmās (Ursa Minor, Carina). Note that, as expected, dSph's closer to the Galactic Center typically show a larger reduction from to , sometimes by more than a factor of 2.

In the same fashion, we then determine the most likely annihilation luminosities for the 18 dSph's by using (9) for an NFW profile to calculate the total luminosity for every simulated subhalo. Additionally, we also determine , the luminosity within 0.3ākpc from the center, motivated by the fact that we only have dynamical evidence for a DM dominated potential out to this radius. Lastly we also consider two measures of the brightness of each halo: , the total expected flux from the dSph, and , the flux from a central region subtending , which is comparable to the angular resolution of Fermi above 3āGeV. thus corresponds to the brightest āpixelāā in a Fermi ray image of a subhalo. These numbers are given in Table 2.

4.1. Current Observational Constraints
Several ACT have performed observations of a handful of dSph's. (i)The H.E.S.S. array (consisting of four 107ām telescopes with a field of view and an energy threshold of āGeV [100]) has obtained an 11 hour exposure of the Sagittarius dwarf galaxy. No ray signal was detected, resulting in a flux limit of ācmās (95% confidence) at āGeV, and a corresponding limit on the cross section of for an NFW profile and for a cored profile (for a higgsinotype neutralino) [46]. Note that the Sagittarius dwarf is undergoing a strong tidal interaction with the Milky Way galaxy [101], and no confident determination of has been possible. (ii)The VERITAS array (consisting of four 144ām telescopes with a field of view and an energy threshold of āGeV [102]) has conducted a 15 hours observation of Willman 1 and 20 hour observations each of Draco and Ursa Minor [47]. No ray signal was detected at a flux limit of 1% of the flux from the Crab Nebula, corresponding to a limit of ācmās (95% confidence) at āGeV [103]. (iii)Additionally the MAGIC [48, 49] and STACEE [50] telescopes have reported observations of Willman 1 and Draco, resulting in comparable or slightly higher flux limits.
To convert the values of in Table 2 into physical fluxes that can directly be compared to these observational limits, it would be necessary to obtain values of the particle physics term of (7) by performing a scan of the DM model parameter space. This is beyond the scope of this work, but a similar analysis has been performed by others [83, 103ā106]. Current ACT observations of dSph's are beginning to directly constrain DM models, and future longer exposure time observations of additional dSph's (in particular Segue 1 and Ursa Major II) with a lower threshold energy hold great potential. We also eagerly await the first Fermi data on fluxes from the known dSph galaxies.
5. Dark Clumps
An annihilation signal from dark clumps not associated with any known luminous stellar counterpart would provide evidence for one of the fundamental implications of the CDM paradigm of structure formation: abundant Galactic substructure. Barring a serendipitous discovery with an ACT, the discovery of such a source will have to rely on allsky surveys, such as provided by Fermi. Of course even a weak and tentative identification of a dark clump with Fermi could be followed up with an ACT.
Unlike for known dSph galaxies, for which we at least have some astronomical observations to guide us, we must rely entirely upon numerical simulations to quantify the prospects of detecting the annihilation signal from dark clumps. Recent significant progress [107] notwithstanding, it is at present not yet possible to perform realistic cosmological hydrodynamic galaxyformation simulations, which include, in addition to the DM dynamics, all the relevant baryonic physics of gas cooling, star formation, supernova and AGN feedback, and so forth. that may have a significant impact on the DM distribution at the centers of massive halos. Instead we make use of the extremely high resolution, purely collisionless DMonly Via Lactea II (VL2) simulation [13], which provides an exquisite view of the clumpiness of the Galactic DM distribution, but at the expense of not capturing all the relevant physics at the baryondominated Galactic center. For the abundance, distribution, and internal properties of the DM subhalos that are the focus of this work, the neglect of baryonic physics is less of a problem, since they are too small to allow for much gas cooling and significant baryonic effects (this is supported by the high masstolight ratios observed in the Milky Way dSph's), although tidal interactions with the Galactic stellar and gaseous disk might significantly affect the population of nearby subhalos.
With a particle mass of and a force softening of 40āpc, the VL2 simulation resolves over 50,000 subhalos today within the host's ākpc (the radius enclosing an average density 200 times the mean matter value). Above 200 particles per halo, the differential subhalo mass function is wellfit by a single power law, , and the cumulative function is [13]. The radial distribution of subhalos is āantibiasedāā with respect to the host halo's density profile, meaning that the mass distribution becomes less clumpy as one approaches the host's center [13, 90]. Similar results have been obtained by the Aquarius group [14, 88]. Typical subhalo concentrations, defined as , grow towards the center, owing to a combination of earlier formation times [108, 109] and stronger tidal stripping of central subhalos: VL2 subhalos on average have a 60 times higher at 8ākpc than at 400ākpc [13]. Note that this also implies 7 times higher annihilation luminosities for central subhalos, since . The counteracting trends of decreasing relative abundance of subhalos and increasing annihilation luminosity towards the center makes it more difficult for (semi)analytical methods to accurately assess the role of subhalos in the Galactic annihilation signal, and motivate future, even higher resolution, numerical simulations of the formation and evolution of Galactic DM (sub)structure. A direct analysis of the VL2 simulations in terms of the detectability with Fermi of individual subhalos was performed by [82]. They found that for reasonable particle physics parameters a handful of subhalos should be able to outshine the astrophysical backgrounds and would be detected at more than significance over the lifetime of the Fermi mission.
As discussed in the previous section, we have directly calculated the annihilation luminosities for all VL2 subhalos using (9) and assuming an NFW density profile. The luminosities would be 40% higher if an Einasto () profile had been adopted instead. We then converted these luminosities to fluxes by dividing by , where the distances were determined for 100 randomly chosen observer locations 8ākpc from the host halo center. The resulting values of are plotted in Figure 3, for just one of the 100 observer positions, as a function of the subhalos' , , and . Although the distributions show quite a bit of scatter, in all three cases a clear trend is apparent of more massive subhalos having higher fluxes. This trend could simply be the result of the higher luminosities of more massive halos, but one might have expected smaller mass subhalos to be brighter, since their greater abundance should result in lower typical distances and hence higher fluxes. This latter effect could be artificially suppressed in the numerical simulations, if smaller mass subhalos, whose dense centers are not as well resolved, were more easily tidally disrupted closer to the Galactic Center, or if the subhalo finding algorithm had trouble identifying low mass halos in the high background density central region. In Figure 4 we plot the subhalos' against their distance to the host halo center . There appears to be a dearth of the lowest subhalos (ākmās) at small distances, but at the moment it is not clear whether this suppression is a real effect or a numerical artifact. It is also worth noting that such small subhalos might be more susceptible to disruption by interactions with the Milky Way's stellar and gaseous disk. At any rate, we can obtain an analytic estimate of the scaling of the typical subhalo flux with by noting that the luminosity scales as and the typical distance as (since ). The typical flux should thus scale as , and would be higher for more massive subhalos at a fixed . Actually lower subhalos might be expected to have higher due to their earlier formation times, but it remains to be seen to what degree this expectation is borne out in numerical simulations.
(a)
(b)
(c)
The points in Figure 3 can be directly compared with the values for the known Milky Way dSph's in Tables 1 and 2: it appears that there are many DM subhalos at least as bright as the known Milky Way dSph's. This impression is confirmed by the top panel of Figure 5, in which we show the cumulative number of subhalos with fluxes greater than and . These distributions were obtained by averaging over 100 randomly chosen observer locations 8 kpc from the host halo center. The mean number of DM subhalos with greater than that of (Carina, Draco, Ursa Minor, Ursa Major, Segue 1) is where the sub and superscripts refer to the 16th and 84th percentile. The corresponding numbers for are . This demonstrates that if a DM annihilation signal from any of the known Milky Way dSph's is detected, then many more DM subhalos should be visible. The plot also implies that Segue 1, the dSph with the highest and of the currently known sample, is unlikely to be the brightest DM subhalo in the sky. Of course some of these additional bright sources could very well have stellar counterparts that have simply been missed so far, due to the limited sky coverage of current surveys or insufficiently deep exposures. To assess what fraction of high flux sources are likely to be genuinely dark clumps without any stars, we split the sample by a limiting value of , ākmās, and ākmās. We assume that DM subhalos below these limits are too small to have been able to form any stars, and hence are truly dark clumps. Of the known dSph's listed in Table 1 only Leo IV falls below these limits. In the bottom panel of Figure 5 we plot , the fraction of subhalos without stars, as a function of the limiting annihilation flux . falls monotonically with , which makes sense given that higher flux sources are typically more massive and hence more likely to host stars. Between 30% and 40% of all DM subhalos brighter than Carina are expected to be dark clumps. This fraction drops to 10% for subhalos brighter than Segue 1.
(a)
(b)
5.1. Boost Factor?
The analysis presented here so far has been limited to known dSph galaxies and clumps resolved in the VL2 simulation, whose resolution limit is set by the available computational resources, and has nothing to do with fundamental physics. Indeed, the CDM expectation is that the clumpiness should continue all the way down to the cutoff in the matter power spectrum, set by collisional damping and free streaming in the early universe Green2005, Loeb2005. For typical WIMP DM, this cutoff occurs at masses of to [27, 28], some 10 to 20 orders of magnitude below VL2's mass resolution. Since the annihilation rate goes as and , this subresolution clumpiness will lead to an enhancement of the total luminosity compared to the smooth mass distribution in the simulation.
The magnitude of this socalled substructure boost factor depends sensitively on the properties of subhalos below the simulation's resolution limit, in particular on the behavior of the concentrationmass relation. A simple power law extrapolation of the contribution of simulated subhalos to the total luminosity of the host halo leads to boosts on order of a a few hundred [55]. More sophisticated (semi)analytical models, accounting for different possible continuations of the concentrationmass relation to lower masses, typically find smaller boosts of around a few tens [81, 82, 105, 110].
More importantly, this boost refers to the enhancement of the total annihilation luminosity of a subhalo, but this is not likely the quantity most relevant for detection. At the distances where subhalos might be detectable as individual sources, their projected size exceeds the angular resolution of today's detectors. The surface brightness profile from annihilations in the smooth DM component would be strongly peaked towards the center (yet probably still resolved by Fermi [82]), owing to the dependence of the annihilation rate. The luminosity contribution from a subhalo's subsubstructure population (i.e., its boost), however, is much less centrally concentrated: at best it follows the subhalo mass density profile , although it might very well even be antibiased. This implies that substructure would preferentially boost the outer regions of a subhalo, where the surface brightness typically remains below the level of astrophysical backgrounds and hence does not contribute much to the detection significance. In other words, the boost factor might apply to , but much less (or not at all) to ; yet it is that is likely to determine whether a given subhalo can be detected with Fermi or an ACT. It thus seems unlikely that the detectability of Galactic subhalos would be significantly enhanced by their own substructure. (This is in contrast to many previous claims in the literature, including some by the present author (e.g., [82]). A reanalysis of that work (in progress) with an improved treatment of the angular dependence of the substructure boost, indeed finds that the boost only weakly increases the number of detectable subhalos.) On the other hand, a substructure boost could be very important for diffuse DM annihilation signals, either from extragalactic sources, where the boost would simply increase the overall amplitude [111], or from Galactic DM, where the boost could affect the amplitude and angular profile of the signal, as well as the power spectrum of its anisotropies [51ā53, 55, 81, 82].
Note also that local substructure boost factors of up to several hundreds, which have been invoked to explain the PAMELA satellite's measurement of an anomalously rising positron fraction at energies above 10āGeV as a DM annihilation signal, appear to be inconsistent with the results from high resolution numerical simulations like VL2 and Aquarius.
6. Summary and Conclusions
In this work we have reviewed the DM annihilation signal from Galactic subhalos. After going over the basics of the annihilation process with a focus on the resulting ray output, we summarized the properties of DM subhalos relevant for estimating their annihilation luminosity. In the remainder of the paper we used the Via Lactea II simulation to assess the strength of the annihilation flux from both known Galactic dSph galaxies as well as from dark clumps not hosting any stars. By matching the distances and central masses of simulated subhalos to the corresponding published values of 18 known dSph's, we were able to infer most probable values, and the 1 scatter around them, for and , and hence for the annihilation luminosity and flux of all dwarfs. According to this analysis, the recently discovered dSph Segue 1 should be the brightest of the known dSph's, followed by Ursa Major II, Ursa Minor, Draco, and Carina. Further, we showed that if any of the known Galactic dSph's are bright enough to be detected, then at least 10 times more subhalos should appear as visible sources. Some of these would be asofyet undiscovered luminous dwarf galaxies, but a significant fraction should correspond to dark clumps not hosting any stars. The fraction of dark clump sources is 10% for subhalos at least as bright as Segue 1 and grows to 40% for subhalos brighter than Carina. Lastly, we briefly considered the role that a substructure boost factor should play in the detectability of individual Galactic dSph's and other DM subhalos. We argued that any boost is unlikely to strongly increase their prospects for detection, since its shallower angular dependence would preferentially boost the outer regions of subhalos, which typically do not contribute much to the detection significance.
Several caveats to these findings are in order. Probably the most important of these is that our simulation completely neglects the effects of baryons. Gas cooling, star formation, and the associated feedback processes are unlikely to strongly affect most subhalos, owing to their low masses. However, tidal interactions with the baryonic components of the Milky Way galaxy might do so. The Sagittarius dSph, for example, is thought to be in the process of complete disruption from tidal interactions with the Milky Way. A second caveat is that our analysis is based on only one, albeit very high resolution, numerical simulation, and so we cannot assess the importance of cosmic variance, or the dependence on cosmological parameters such as and . Other work has found considerable halotohalo scatter [14, 112, 113], with a factor of 2 variance in the total subhalo abundance, for example.
These caveats motivate further study and future, higher resolution numerical simulations, including the effects of baryonic physics. The characterization of the Galactic DM annihilation signal is of crucial importance in guiding observational efforts to shed light on the nature of DM. We are hopeful that in the next few years the promise of a DM annihilation signal will come to fruition, and will help us to unravel this puzzle.
Acknowledgments
Support for this work was provided by the William L. Loughlin membership at the Institute for Advanced Study. The author would like to thank his collaborators from the Via Lactea Project for their expertise and invaluable contributions.
References
 A. Klypin, A. V. Kravtsov, O. Valenzuela, and F. Prada, āWhere are the missing galactic satellites?ā Astrophysical Journal, vol. 522, no. 1, part 1, pp. 82ā92, 1999. View at: Google Scholar
 B. Moore, S. Ghigna, F. Governato et al., āDark matter substructure within galactic halos,ā Astrophysical Journal, vol. 524, no. 1, part 2, pp. L19āL22, 1999. View at: Google Scholar
 M. Mateo, āDwarf galaxies of the local group,ā Annual Review of Astronomy and Astrophysics, vol. 36, pp. 435ā506, 1998. View at: Publisher Site  Google Scholar
 B. Willman, M. R. Blanton, A. A. West et al., āA new Milky Way companion: unusual globular cluster or extreme dwarf satellite?ā Astronomical Journal, vol. 129, no. 6, pp. 2692ā2700, 2005. View at: Publisher Site  Google Scholar
 B. Willman, J. J. Dalcanton, D. MartinezDelgado et al., āA new Milky Way dwarf galaxy in URSA Major,ā Astrophysical Journal, vol. 626, no. 2, part 2, pp. L85āL88, 2005. View at: Publisher Site  Google Scholar
 V. Belokurov, D. B. Zucker, N. W. Evans et al., āThe field of streams: sagittarius and its siblings,ā Astrophysical Journal, vol. 642, no. 2, part 2, pp. L137āL140, 2006. View at: Publisher Site  Google Scholar
 D. B. Zucker, V. Belokurov, N. W. Evans et al., āA new Milky Way dwarf satellite in canes venatici,ā Astrophysical Journal, vol. 643, no. 2, part 2, pp. L103āL106, 2006. View at: Publisher Site  Google Scholar
 D. B. Zucker, V. Belokurov, N. W. Evans et al., āA curious Milky Way satellite in Ursa Major,ā Astrophysical Journal, vol. 650, no. 1, part 2, pp. L41āL44, 2006. View at: Publisher Site  Google Scholar
 T. Sakamoto and T. Hasegawa, āDiscovery of a faint old stellar system at 150 kpc,ā Astrophysical Journal, vol. 653, no. 1, part 2, pp. L29āL32, 2006. View at: Publisher Site  Google Scholar
 V. Belokurov, D. B. Zucker, N. W. Evans et al., āCats and dogs, hair and a hero: a quintet of new Milky Way companions,ā Astrophysical Journal, vol. 654, no. 2, part 1, pp. 897ā906, 2007. View at: Publisher Site  Google Scholar
 M. J. Irwin, V. Belokurov, N. W. Evans et al., āDiscovery of an unusual dwarf galaxy in the outskirts of the Milky Way,ā Astrophysical Journal, vol. 656, no. 1, part 2, pp. L13āL16, 2007. View at: Publisher Site  Google Scholar
 S. M. Walsh, H. Jerjen, and B. Willman, āA pair of boötes: a new Milky Way satellite,ā Astrophysical Journal, vol. 662, no. 2, part 2, pp. L83āL86, 2007. View at: Publisher Site  Google Scholar
 J. Diemand, M. Kuhlen, P. Madau et al., āClumps and streams in the local dark matter distribution,ā Nature, vol. 454, no. 7205, pp. 735ā738, 2008. View at: Publisher Site  Google Scholar
 V. Springel, J. Wang, M. Vogelsberger et al., āThe Aquarius Project: the subhaloes of galactic haloes,ā Monthly Notices of the Royal Astronomical Society, vol. 391, no. 4, pp. 1685ā1711, 2008. View at: Publisher Site  Google Scholar
 J. Stadel, D. Potter, B. Moore et al., āQuantifying the heart of darkness with GHALO—a multibillion particle simulation of a galactic halo,ā Monthly Notices of the Royal Astronomical Society, vol. 398, no. 1, pp. L21āL25, 2009. View at: Publisher Site  Google Scholar
 S. D. M. White and M. J. Rees, āCore condensation in heavy halos—a twostage theory for galaxy formation and clustering,ā Monthly Notices of the Royal Astronomical Society, vol. 183, pp. 341ā358, 1978. View at: Google Scholar
 G. R. Blumenthal, S. M. Faber, J. R. Primack, and M. J. Rees, āFormation of galaxies and largescale structure with cold dark matter,ā Nature, vol. 311, no. 5986, pp. 517ā525, 1984. View at: Publisher Site  Google Scholar
 A. Dekel and J. Silk, āThe origin of dwarf galaxies, cold dark matter, and biased galaxy formation,ā Astrophysical Journal, vol. 303, pp. 39ā55, 1986. View at: Publisher Site  Google Scholar
 J. S. Bullock, A. V. Kravtsov, and D. H. Weinberg, āReionization and the abundance of galactic satellites,ā Astrophysical Journal, vol. 539, no. 2, part 1, pp. 517ā521, 2000. View at: Google Scholar
 A. V. Kravtsov, O. Y. Gnedin, and A. A. Klypin, āThe tumultuous lives of galactic dwarfs and the missing satellites problem,ā Astrophysical Journal, vol. 609, no. 2, part 1, pp. 482ā497, 2004. View at: Publisher Site  Google Scholar
 L. Mayer, C. Mastropietro, J. Wadsley, J. Stadel, and B. Moore, āSimultaneous ram pressure and tidal stripping; how dwarf spheroidals lost their gas,ā Monthly Notices of the Royal Astronomical Society, vol. 369, no. 3, pp. 1021ā1038, 2006. View at: Publisher Site  Google Scholar
 P. Madau, J. Diemand, and M. Kuhlen, āDark matter subhalos and the dwarf satellites of the Milky Way,ā Astrophysical Journal, vol. 679, no. 2, pp. 1260ā1271, 2008. View at: Publisher Site  Google Scholar
 S. E. Koposov, J. Yoo, H.W. Rix, D. H. Weinberg, A. V. Macciò, and J. M. Escudé, āA quantitative explanation of the observed population of Milky Way satellite galaxies,ā Astrophysical Journal, vol. 696, pp. 2179ā2194, 2009. View at: Publisher Site  Google Scholar
 A. V. Maccio', X. Kang, F. Fontanot, R. S. Somerville, S. E. Koposov, and P. Monaco, āOn the origin and properties of Ultrafaint Milky Way Satellites in a LCDM Universe,ā submitted, http://arxiv.org/abs/0903.4681. View at: Google Scholar
 A. M. Green, S. Hofmann, and D. J. Schwarz, āThe first WIMPy halos,ā Journal of Cosmology and Astroparticle Physics, vol. 2005, no. 8, article 3, pp. 41ā71, 2005. View at: Publisher Site  Google Scholar
 A. Loeb and M. Zaldarriaga, āSmallscale power spectrum of cold dark matter,ā Physical Review D, vol. 71, no. 10, Article ID 103520, 7 pages, 2005. View at: Publisher Site  Google Scholar
 S. Profumo, K. Sigurdson, and M. Kamionkowski, āWhat mass are the smallest protohalos?ā Physical Review Letters, vol. 97, no. 3, Article ID 031301, 4 pages, 2006. View at: Publisher Site  Google Scholar
 T. Bringmann, āParticle models and the smallscale structure of dark matter,ā New Journal of Physics, vol. 11, no. 10, Article ID 105027, 105027 pages, 2009. View at: Publisher Site  Google Scholar
 E. W. Kolb and I. I. Tkachev, āFemtolensing and picolensing by axion miniclusters,ā Astrophysical Journal, vol. 460, no. 1, part 2, pp. L25āL28, 1996. View at: Google Scholar
 C. Schmid, D. J. Schwarz, and P. Widerin, āAmplification of cosmological inhomogeneities by the QCD transition,ā Physical Review D, vol. 59, no. 4, Article ID 043517, 20 pages, 1999. View at: Publisher Site  Google Scholar
 P. Sikivie, āAxion cosmology,ā in Axions: Theory, Cosmology, and Experimental Searches, M. Kuster, G. Raffelt, and B. Beltrán, Eds., vol. 741 of Lecture Notes in Physics, p. 19, Springer, Berlin, Germany, 2008. View at: Google Scholar
 L. Bergström, J. Edsjö, P. Gondolo, and P. Ullio, āClumpy neutralino dark matter,ā Physical Review D, vol. 59, no. 4, Article ID 043506, 11 pages, 1999. View at: Publisher Site  Google Scholar
 L. E. Strigari, J. S. Bullock, M. Kaplinghat et al., āA common mass scale for satellite galaxies of the Milky Way,ā Nature, vol. 454, no. 7208, pp. 1096ā1097, 2008. View at: Publisher Site  Google Scholar
 R. A. Ibata, G. F. Lewis, M. J. Irwin, and T. Quinn, āUncovering cold dark matter halo substructure with tidal streams,ā Monthly Notices of the Royal Astronomical Society, vol. 332, no. 4, pp. 915ā920, 2002. View at: Publisher Site  Google Scholar
 K. V. Johnston, D. N. Spergel, and C. Haydn, āHow lumpy is the Milky Way's dark matter halo?ā Astrophysical Journal, vol. 570, no. 2, part 1, pp. 656ā664, 2002. View at: Publisher Site  Google Scholar
 J. Peñarrubia, A. J. Benson, D. MartínezDelgado, and H. W. Rix, āModeling tidal streams in evolving dark matter halos,ā Astrophysical Journal, vol. 645, no. 1, part 1, pp. 240ā255, 2006. View at: Publisher Site  Google Scholar
 J. M. SiegalGaskins and M. Valluri, āSignatures of $\mathrm{\Lambda}$CDM substructure in tidal debris,ā Astrophysical Journal, vol. 681, no. 1, pp. 40ā52, 2008. View at: Publisher Site  Google Scholar
 G. Toth and J. P. Ostriker, āGalactic disks, infall, and the global value of $\mathrm{\Omega}$,ā Astrophysical Journal, vol. 389, no. 1, pp. 5ā26, 1992. View at: Google Scholar
 J. I. Read, G. Lake, O. Agertz, and V. P. Debattista, āThin, thick and dark discs in $\mathrm{\Lambda}$CDM,ā Monthly Notices of the Royal Astronomical Society, vol. 389, no. 3, pp. 1041ā1057, 2008. View at: Publisher Site  Google Scholar
 S. Kazantzidis, A. R. Zentner, A. V. Kravtsov, J. S. Bullock, and V. P. Debattista, āCold dark matter substructure and galactic disks. II. Dynamical effects of hierarchical satellite accretion,ā Astrophysical Journal, vol. 700, no. 2, pp. 1896ā1920, 2009. View at: Publisher Site  Google Scholar
 O. Adriani, G. C. Barbarino, G. A. Bazilevskaya et al., āAn anomalous positron abundance in cosmic rays with energies 1.5–100 GeV,ā Nature, vol. 458, no. 7238, pp. 607ā609, 2009. View at: Publisher Site  Google Scholar
 J. Chang, J. H. Adams Jr., H. S. Ahn et al., āAn excess of cosmic ray electrons at energies of 300–800 GeV,ā Nature, vol. 456, no. 7220, pp. 362ā365, 2008. View at: Publisher Site  Google Scholar
 S. Torii, T. Yamagami, T. Tamura et al., āHighenergy electron observations by PPBBETS flight in Antarctica,ā submitted, http://arxiv.org/abs/0809.0760. View at: Google Scholar
 A. A. Abdo, M. Ackermann, M. Ajello et al., āMeasurement of the cosmic ray ${e}^{+}+{e}^{}$ spectrum from 20 GeV to 1 TeV with the fermi large area telescope,ā Physical Review Letters, vol. 102, no. 18, Article ID 181101, 2009. View at: Publisher Site  Google Scholar
 F. Aharonian, A. G. Akhperjanian, G. Anton et al., āProbing the ATIC peak in the cosmicray electron spectrum with H.E.S.S,ā submitted, http://arxiv.org/abs/0905.0105. View at: Google Scholar
 F. Aharonian, A. G. Akhperjanian, A. R. BazerBachi et al., āObservations of the Sagittarius dwarf galaxy by the HESS experiment and search for a dark matter signal,ā Astroparticle Physics, vol. 29, no. 1, pp. 55ā62, 2008. View at: Publisher Site  Google Scholar
 C. M. Hui, āVERITAS Observations of Extragalactic NonBlazars,ā in High Energy GammaRay Astronomy: Proceedings of the 4th International Meeting on High Energy GammaRay Astronomy, F. A. Aharonian, W. Hofmann, and F. Rieger, Eds., vol. 1085 of American Institute of Physics Conference Series, pp. 407ā410, 2008. View at: Google Scholar
 J. Albert, E. Aliu, H. Anderhub et al., āUpper limit for $\gamma $ray emission above 140 GeV from the dwarf spheroidal galaxy draco,ā Astrophysical Journal, vol. 679, no. 1, pp. 428ā431, 2008. View at: Publisher Site  Google Scholar
 E. Aliu, H. Anderhub, L. A. Antonelli et al., āUpper limits on the vhe gammaray emission from the willman 1 satellite galaxy with the magic telescope,ā Astrophysical Journal, vol. 697, no. 2, pp. 1299ā1304, 2009. View at: Publisher Site  Google Scholar
 D. D. Driscoll, C. E. Covault, J. Ball et al., āSearch for dark matter annihilation in Draco with the Solar Tower Atmospheric Cherenkov Effect Experiment,ā Physical Review D, vol. 78, no. 8, Article ID 087101, 4 pages, 2008. View at: Publisher Site  Google Scholar
 J. M. SiegalGaskins, āRevealing dark matter substructure with anisotropies in the diffuse gammaray background,ā Journal of Cosmology and Astroparticle Physics, vol. 2008, no. 10, article 40, 2008. View at: Publisher Site  Google Scholar
 S. Ando, āGammaray background anisotropy from Galactic dark matter substructure,ā Physical Review D, vol. 80, no. 2, Article ID 023520, 16 pages, 2009. View at: Publisher Site  Google Scholar
 M. Fornasa, L. Pieri, G. Bertone, and E. Branchini, āAnisotropy probe of galactic and extragalactic dark matter annihilations,ā Physical Review D, vol. 80, no. 2, Article ID 023518, 2009. View at: Publisher Site  Google Scholar
 F. Stoehr, S. D. M. White, V. Springel, G. Tormen, and N. Yoshida, āDark matter annihilation in the halo of the Milky Way,ā Monthly Notices of the Royal Astronomical Society, vol. 345, no. 4, pp. 1313ā1322, 2003. View at: Publisher Site  Google Scholar
 V. Springel, S. D. M. White, C. S. Frenk et al., āProspects for detecting supersymmetric dark matter in the Galactic halo,ā Nature, vol. 456, no. 7218, pp. 73ā76, 2008. View at: Publisher Site  Google Scholar
 P. Brun, T. Delahaye, J. Diemand, S. Profumo, and P. Salati, āCosmic ray lepton puzzle in the light of cosmological $N$body simulations,ā Physical Review D, vol. 80, no. 3, Article ID 035023, 5 pages, 2009. View at: Publisher Site  Google Scholar
 M. Kuhlen and D. Malyshev, āATIC, PAMELA, HESS, and Fermi data and nearby dark matter subhalos,ā Physical Review D, vol. 79, no. 12, Article ID 123517, 12 pages, 2009. View at: Publisher Site  Google Scholar
 N. ArkaniHamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, āA theory of dark matter,ā Physical Review D, vol. 79, no. 1, Article ID 015014, 16 pages, 2009. View at: Publisher Site  Google Scholar
 M. Lattanzi and J. Silk, āCan the WIMP annihilation boost factor be boosted by the Sommerfeld enhancement?ā Physical Review D, vol. 79, no. 8, Article ID 083523, 8 pages, 2009. View at: Publisher Site  Google Scholar
 L. Pieri, M. Lattanzi, and J. Silk, āConstraining the dark matter annihilation crosssection with Cherenkov telescope observations of dwarf galaxies,ā Monthly Notices of the Royal Astronomical Society, vol. 399, no. 4, pp. 2033ā2040, 2009. View at: Publisher Site  Google Scholar
 M. Kuhlen, P. Madau, and J. Silk, āExploring dark matter with milky way substructure,ā Science, vol. 325, no. 5943, pp. 970ā973, 2009. View at: Publisher Site  Google Scholar
 M. S. Turner, āWindows on the axion,ā Physics Report, vol. 197, no. 2, pp. 67ā97, 1990. View at: Google Scholar
 E. W. Kolb and M. S. Turner, āThe early universe,ā in Frontiers in Physics, AddisonWesley, Reading, Mass, USA, 1990. View at: Google Scholar
 G. Jungman, M. Kamionkowski, and K. Griest, āSupersymmetric dark matter,ā Physics Report, vol. 267, no. 56, pp. 195ā373, 1996. View at: Google Scholar
 G. Bertone, D. Hooper, and J. Silk, āParticle dark matter: evidence, candidates and constraints,ā Physics Reports, vol. 405, no. 56, pp. 279ā390, 2005. View at: Publisher Site  Google Scholar
 G. Hinshaw, J. L. Weiland, R. S. Hill et al., āFiveyear wilkinson microwave anisotropy probe observations: data processing, sky maps and basic results,ā Astrophysical Journal, vol. 180, no. 2, pp. 225ā245, 2009. View at: Publisher Site  Google Scholar
 K. Griest and D. Seckel, āThree exceptions in the calculation of relic abundances,ā Physical Review D, vol. 43, no. 10, pp. 3191ā3203, 1991. View at: Publisher Site  Google Scholar
 P. Gondolo and G. Gelmini, āCosmic abundances of stable particles: improved analysis,ā Nuclear Physics B, vol. 360, no. 1, pp. 145ā179, 1991. View at: Google Scholar
 T. Bringmann, L. Bergström, and J. Edsjö, āNew gammaray contributions to supersymmetric dark matter annihilation,ā Journal of High Energy Physics, vol. 2008, no. 1, p. 49, 2008. View at: Publisher Site  Google Scholar
 M. Gustafsson, E. Lundström, L. Bergström, and J. Edsjö, āSignificant gamma lines from inert higgs dark matter,ā Physical Review Letters, vol. 99, no. 4, Article ID 041301, 4 pages, 2007. View at: Publisher Site  Google Scholar
 T. Sjöstrand, āHighenergyphysics event generation with PYTHIA 5.7 and JETSET 7.4,ā Computer Physics Communications, vol. 82, no. 1, pp. 74ā89, 1994. View at: Google Scholar
 P. Gondolo, J. Edsjö, P. Ullio, L. Bergström, M. Schelke, and E. A. Baltz, āDarkSUSY: computing supersymmetric dark matter properties numerically,ā Journal of Cosmology and Astroparticle Physics, vol. 2004, no. 7, article 8, pp. 145ā179, 2004. View at: Publisher Site  Google Scholar
 C. CalcáneoRoldán and B. Moore, āSurface brightness of dark matter: unique signatures of neutralino annihilation in the galactic halo,ā Physical Review D, vol. 62, no. 12, Article ID 123005, 10 pages, 2000. View at: Google Scholar
 E. A. Baltz, C. Briot, P. Salati, R. Taillet, and J. Silk, āDetection of neutralino annihilation photons from external galaxies,ā Physical Review D, vol. 61, no. 2, Article ID 023514, 10 pages, 2000. View at: Google Scholar
 A. Tasitsiomi and A. V. Olinto, āDetectability of neutralino clumps via atmospheric Cherenkov telescopes,ā Physical Review D, vol. 66, no. 8, Article ID 083006, 16 pages, 2002. View at: Google Scholar
 R. Aloisio, P. Blasi, and A. V. Olinto, āGammaray constraints on neutralino dark matter clumps in the galactic halo,ā Astrophysical Journal, vol. 601, no. 1, part 1, pp. 47ā53, 2004. View at: Google Scholar
 N. W. Evans, F. Ferrer, and S. Sarkar, āA travel guide to the dark matter annihilation signal,ā Physical Review D, vol. 69, no. 12, Article ID 123501, 10 pages, 2004. View at: Publisher Site  Google Scholar
 S. M. Koushiappas, A. R. Zentner, and T. P. Walker, āObservability of gamma rays from neutralino annihilations in the Milky Way substructure,ā Physical Review D, vol. 69, no. 4, Article ID 043501, 11 pages, 2004. View at: Publisher Site  Google Scholar
 S. M. Koushiappas, āProper motion of gamma rays from microhalo sources,ā Physical Review Letters, vol. 97, no. 19, Article ID 191301, 4 pages, 2006. View at: Publisher Site  Google Scholar
 J. Diemand, M. Kuhlen, and P. Madau, āDark matter substructure and gammaray annihilation in the Milky Way halo,ā Astrophysical Journal, vol. 657, no. 1, part 1, pp. 262ā270, 2007. View at: Publisher Site  Google Scholar
 L. Pieri, G. Bertone, and E. Branchini, āDark matter annihilation in substructures revised,ā Monthly Notices of the Royal Astronomical Society, vol. 384, no. 4, pp. 1627ā1637, 2008. View at: Publisher Site  Google Scholar
 M. Kuhlen, J. Diemand, and P. Madau, āThe dark matter annihilation signal from galactic substructure: predictions for glast,ā Astrophysical Journal, vol. 686, no. 1, pp. 262ā278, 2008. View at: Publisher Site  Google Scholar
 L. E. Strigari, S. M. Koushiappas, J. S. Bullock et al., āThe most darkmatterdominated galaxies: predicted gammaray signals from the faintest Milky Way dwarfs,ā Astrophysical Journal, vol. 678, no. 2, pp. 614ā620, 2008. View at: Publisher Site  Google Scholar
 E. A. Baltz, J. E. Taylor, and L. L. Wai, āCan astrophysical gammaray sources mimic dark matter annihilation in galactic satellites?ā Astrophysical Journal, vol. 659, no. 2, part 2, pp. L125āL128, 2007. View at: Publisher Site  Google Scholar
 J. F. Navarro, C. S. Frenk, and S. D. M. White, āA universal density profile from hierarchical clustering,ā Astrophysical Journal, vol. 490, no. 2, part 1, pp. 493ā508, 1997. View at: Google Scholar
 A. V. Macciò, A. A. Dutton, F. C. van den Bosch, B. Moore, D. Potter, and J. Stadel, āConcentration, spin and shape of dark matter haloes: scatter and the dependence on mass and environment,ā Monthly Notices of the Royal Astronomical Society, vol. 378, no. 1, pp. 55ā71, 2007. View at: Publisher Site  Google Scholar
 J. Diemand, B. Moore, and J. Stadel, āConvergence and scatter of cluster density profiles,ā Monthly Notices of the Royal Astronomical Society, vol. 353, no. 2, pp. 624ā632, 2004. View at: Publisher Site  Google Scholar
 J. F. Navarro, A. Ludlow, V. Springel et al., āThe diversity and similarity of simulated cold dark matter halos,ā submitted, http://arxiv.org/abs/0810.1522. View at: Google Scholar
 B. Allgood, R. A. Flores, J. R. Primack et al., āThe shape of dark matter haloes: dependence on mass, redshift, radius and formation,ā Monthly Notices of the Royal Astronomical Society, vol. 367, no. 4, pp. 1781ā1796, 2006. View at: Publisher Site  Google Scholar
 M. Kuhlen, J. Diemand, and P. Madau, āThe shapes, orientation, and alignment of galactic dark matter subhalos,ā Astrophysical Journal, vol. 671, no. 2, pp. 1135ā1146, 2007. View at: Publisher Site  Google Scholar
 G. L. Bryan and M. L. Norman, āStatistical properties of Xray clusters: analytic and numerical comparisons,ā Astrophysical Journal, vol. 495, no. 1, part 1, pp. 80ā99, 1998. View at: Google Scholar
 J. T. Kleyna, M. I. Wilkinson, N. Evans, and G. Gilmore, āUrsa Major: a missing lowmass CDM halo?ā Astrophysical Journal, vol. 630, no. 2, part 2, pp. L141āL144, 2005. View at: Publisher Site  Google Scholar
 R. R. Muñoz, J. L. Carlin, P. M. Frinchaboy, D. L. Nidever, S. R. Majewski, and R. J. Patterson, āExploring halo substructure with giant stars: the dynamics and metallicity of the dwarf spheroidal in Boötes,ā Astrophysical Journal, vol. 650, no. 1, part 2, pp. L51āL54, 2006. View at: Publisher Site  Google Scholar
 N. F. Martin, R. A. Ibata, S. C. Chapman, M. Irwin, and G. F. Lewis, āA Keck/DEIMOS spectroscopic survey of faint Galactic satellites: searching for the least massive dwarf galaxies,ā Monthly Notices of the Royal Astronomical Society, vol. 380, no. 1, pp. 281ā300, 2007. View at: Publisher Site  Google Scholar
 J. D. Simon and M. Geha, āThe kinematics of the ultrafaint Milky Way satellites: solving the missing satellite problem,ā Astrophysical Journal, vol. 670, no. 1, pp. 313ā331, 2007. View at: Publisher Site  Google Scholar
 S. Koposov, V. Belokurov, N. W. Evans et al., āThe luminosity function of the Milky Way satellites,ā Astrophysical Journal, vol. 686, no. 1, pp. 279ā291, 2008. View at: Publisher Site  Google Scholar
 E. J. Tollerud, J. S. Bullock, L. E. Strigari, and B. Willman, āHundreds of Milky Way satellites? Luminosity bias in the satellite luminosity function,ā Astrophysical Journal, vol. 688, no. 1, pp. 277ā289, 2008. View at: Publisher Site  Google Scholar
 S. M. Walsh, B. Willman, and H. Jerjen, āThe invisibles: a detection algorithm to trace the faintest Milky Way satellites,ā Astronomical Journal, vol. 137, no. 1, pp. 450ā469, 2009. View at: Publisher Site  Google Scholar
 J. Grcevich and M. E. Putman, ā${\text{H}}_{\text{I}}$ in local group dwarf galaxies and stripping by the galactic halo,ā Astronomical Journal, vol. 696, pp. 385ā395, 2009. View at: Publisher Site  Google Scholar
 K. Bernlöhr, O. Carrol, R. Cornils et al., āThe optical system of the H.E.S.S. imaging atmospheric Cherenkov telescopes—part I: layout and components of the system,ā Astroparticle Physics, vol. 20, no. 2, pp. 111ā128, 2003. View at: Publisher Site  Google Scholar
 D. MartínezDelgado, M. Á. GómezFlechoso, A. Aparicio, and R. Carrera, āTracing out the northern tidal stream of the Sagittarius dwarf spheroidal galaxy,ā Astrophysical Journal, vol. 601, no. 1, part 1, pp. 242ā259, 2004. View at: Google Scholar
 J. Holder, V. A. Acciari, E. Aliu et al., āStatus of the VERITAS observatory,ā in High Energy GammaRay Astronomy: Proceedings of the 4th International Meeting on High Energy GammaRay Astronomy, F. A. Aharonian, W. Hofmann, and F. Rieger, Eds., vol. 1085 of American Institute of Physics Conference Series, pp. 657ā660, 2008. View at: Google Scholar
 R. Essig, N. Sehgal, and L. E. Strigari, āBounds on cross sections and lifetimes for dark matter annihilation and decay into charged leptons from gammaray observations of dwarf galaxies,ā Physical Review D, vol. 80, no. 2, Article ID 023506, 2009. View at: Publisher Site  Google Scholar
 J. Bovy, āSubstructure boosts to dark matter annihilation from Sommerfeld enhancement,ā Physical Review D, vol. 79, no. 8, Article ID 083539, 2009. View at: Publisher Site  Google Scholar
 G. D. Martinez, J. S. Bullock, M. Kaplinghat, L. E. Strigari, and R. Trotta, āIndirect dark matter detection from dwarf satellites: joint expectations from astrophysics and supersymmetry,ā Journal of Cosmology and Astroparticle Physics, vol. 6, p. 14, 2009. View at: Publisher Site  Google Scholar
 L. Pieri, A. Pizzella, E. M. Corsini, E. Dalla Bontà, and F. Bertola, āCould the Fermi Large Area Telescope detect $\gamma $rays from dark matter annihilation in the dwarf galaxies of the local group?ā Astronomy and Astrophysics, vol. 496, no. 2, pp. 351ā360, 2009. View at: Publisher Site  Google Scholar
 F. Governato, B. Willman, L. Mayer et al., āForming disc galaxies in $\mathrm{\Lambda}$CDM simulations,ā Monthly Notices of the Royal Astronomical Society, vol. 374, no. 4, pp. 1479ā1494, 2007. View at: Publisher Site  Google Scholar
 J. Diemand, P. Madau, and B. Moore, āThe distribution and kinematics of early high$\sigma $ peaks in presentday haloes: implications for rare objects and old stellar populations,ā Monthly Notices of the Royal Astronomical Society, vol. 364, no. 2, pp. 367ā383, 2005. View at: Publisher Site  Google Scholar
 B. Moore, J. Diemand, P. Madau, M. Zemp, and J. Stadel, āGlobular clusters, satellite galaxies and stellar haloes from early dark matter peaks,ā Monthly Notices of the Royal Astronomical Society, vol. 368, no. 2, pp. 563ā570, 2006. View at: Google Scholar
 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 gammaray experiments,ā Physical Review D, vol. 75, no. 8, Article ID 083526, 2007. View at: Publisher Site  Google Scholar
 P. Ullio, L. Bergström, J. Edsjö, and C. Lacey, āCosmological dark matter annihilations into $\gamma $ rays: a closer look,ā Physical Review D, vol. 66, no. 12, Article ID 123502, 2002. View at: Publisher Site  Google Scholar
 D. Reed, F. Governato, T. Quinn, J. Gardner, J. Stadel, and G. Lake, āDark matter subhaloes in numerical simulations,ā Monthly Notices of the Royal Astronomical Society, vol. 359, no. 4, pp. 1537ā1548, 2005. View at: Publisher Site  Google Scholar
 T. Ishiyama, T. Fukushige, and J. Makino, āVariation of the subhalo abundance in dark matter halos,ā The Astrophysical Journal, vol. 696, pp. 2115ā2125, 2009. View at: Google Scholar
Copyright
Copyright © 2010 Michael Kuhlen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.