The Dark Matter Annihilation Signal from Dwarf Galaxies and Subhalos

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 gamma-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.


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,4,5,6,7,8,9,10,11,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,14,15] resolve close to 100,000 individual self-bound clumps of DM within the Galactic virial volume -remnants of the hierarchical build-up of the Milky Way's DM halo. A consensus seems to be emerging that this discrepancy is not a short-coming 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 aren't [18,19,20,21,22,23,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 free-streaming scale of 10 −12 to 10 −4 M ⊙ [25,26]. 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 [27]. 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 by-products 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 DM-dominated potentials [31], or from perturbations of cold stellar structures like globular cluster tidal streams [32,33,34,35] or the heating of the Milky Way's stellar disk [36,37,38], 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 [39], where standard models of cosmic ray propagation predict a decreasing fraction; the ATIC [40] and PPB-BETS [41] balloon-borne experiments reported a surprisingly large total electron and positron flux at ∼ 500 GeV, although recent Fermi [42] and H.E.S.S. data [43] 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 Gamma-ray 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., VERI-TAS, 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 [44,45,46,47,48].
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. 49,50,51], the relative strength of the signal from individual subhalos compared with that from the Galactic Center or an annulus around it [28,52], the effect of a nearby DM subhalo on the amplitude and spectrum of the local high energy electron and positron flux [53,54], and the role of the Sommerfeld enhancement [55] on the DM annihilation rate and its implications for substructure signals [56,57,58].
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.

DARK MATTER ANNIHILATION
If DM is made up of a so-called "thermal relic" particle 1 , 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 [60,61,62]. 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 (T ≫ m χ ), their creation and destruction rates are balanced, and hence their co-moving abundance remains constant. Once the universe cools below the DM particle's rest-mass (T < m χ ), its equilibrium abundance is suppressed by a Boltzmann factor exp(−m χ /T ). 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 H(a), the Hubble constant at cosmological scale factor a) exceeds the annihilation rate Γ(a) = n σv , i.e. when DM particles can no longer find each other to annihilate. The co-moving number density of DM particles is then fixed at a "freeze-out" temperature that turns out to be approximately T f ≃ m χ /20, 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 ρ crit = 3H 2 0 /8πG, and σv , the thermally averaged velocity-weighted annihilation cross section: Note that this relation is independent of m χ . The WMAP satellite's measurement of the DM density is ω χ = 0.1131 ± 0.0034 [63], implying a value of A more accurate determination of σv 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 co-annihilations into other, nearly degenerate dark sector particles [e.g. 64,65]. It is a remarkable coincidence that this value of σv is close to what one expects for a cross section set by the weak interaction. This is the so-called "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, m ν < 0.63 eV [63], implies a cosmological mass density of only ω ν = 7.1 × 10 −3 . 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, although models with extra dimensions have received a lot of attention in recent years as well. Both of these approaches offer good DM particle candidates: the lightest supersymmetric particle (LSP), typically a neutralino in R-parity conserving supersymmetry, and the lightest Kaluza-Klein particle (LKP), typically the B (1) particle, the first Kaluza-Klein 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, Hooper & Silk [62].
The direct products of the annihilation of two DM particles are strongly model dependent. Typi- Internal bremsstrahlung photons associated with charged annihilation products, either in the form of final state radiation (FSR) from external legs or as virtual internal bremsstrahlung (VIB) from the exchange of virtual charged particles. Bottom: Mono-chromatic line signals from the prompt annihilation into two photons or a photon and Z boson. This process occurs only at loop level, and hence is typically strongly suppressed.
cal channels include annihilations into charged leptons (e + e − , µ + µ − , τ + τ − ), quark-antiquark pairs, and gauge and Higgs bosons (W + W − , Z, h). In the end, however, the decay and hadronization of these annihilation products results in only three types of emissions: (i) high energy neutrinos, (ii) relativistic electrons and protons and their anti-particles, 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 Fig. 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, π 0 → γγ, arising in the hadronization of the primary annihilation products. Since the DM particles are nonrelativistic, their annihilation results in a pair of mono-energetic particles with energy equal to m χ , 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 m χ .
(ii) An important additional contribution at high energies (E < ∼ m χ ) arises from the internal bremsstrahlung process [66], 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 E ∼ m χ 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 m χ .
(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 [67]). On the other hand, the resulting photons would be mono-chromatic, 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 E = m χ , in some models it is also possible to annihilate into a photon and a Z boson, χχ → γZ, and this process would result in a somewhat broadened (due to the mass of the Z) The relative importance of these three γ-ray production channels and the resulting spectrum dN γ /dE 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 Monte-Carlo event generator [68], which is also contained in the popular DarkSUSY package [69].

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 [27,28,29,30,70,71,72,73,74,75,76,77,78]. Compared to diffuse γ-ray annihilation signals, these discrete sources should be easier to distinguish from astrophysical backgrounds and foregrounds [79], 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 time-independent, 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.
We can distinguish between DM subhalos hosting a Milky Way dwarf satellite galaxy and dark clumps that, for whatever reason, don't 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) [80] profile over a wide range of masses [81,82], The parameter r s indicates the radius at which the logarithmic slope γ(r) ≡ d ln ρ d ln r = −2, and ρ(r s ) = ρ s . The very highest resolution simulations have recently provided some indications of a flattening of the density profile in the innermost regions [15,83]. 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 α ≈ 0.17 ± 0.03 in numerical simulations [83]. Note that the two density profiles actually do not differ very much until r ≪ r s (cf. top panel of Figure 2). Simulated DM halos are of course not perfectly spherically symmetric, and instead typically exhibit prolate or triaxial iso-density contours that become more elongated towards the center [84]. The degree of prolateness decreases with mass, and galactic subhalos have axis ratios of > ∼ 0.7 [85]. The "virial" radius R vir of a halo is defined as the radius enclosing a mean density equal to ∆ vir ρ 0 , where ∆ vir ≈ 389 [86] and ρ 0 is the mean density of the universe. The corresponding virial mass M vir is the mass within R vir , and a halo's concentration can then be defined as c = R vir /r s . 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 r s in low mass subhalos. For this reason we prefer to work with V max , the maximum of the circular velocity curve V c (r) 2 = GM (< r)/r and a proxy for a subhalo's mass, and r Vmax , the radius at which V max occurs. These quantities are much more robustly determined for subhalos in numerical simulations than (M, c). Note that even (V max , r Vmax ) 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 V peak , the largest value of V max that a subhalo ever acquired during its lifetime (i.e. before tidal stripping began to lower its V max ) and r Vpeak , 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 ρ(r) 2 , L has dimensions of (mass) 2 (length) −3 , and we express it in units of M 2 ⊙ pc −3 . In order to convert to a conventional luminosity, one must multiply by a particle physics term, where m χ is the mass of the DM particle and σv the thermally averaged velocity-weighted 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 E th and below a maximum energy of E max . In that case the effective γ-ray luminosity is where dN γ /dE 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, L is much more centrally concentrated than M : ∼ 90% of the total luminosity is produced within r s , compared with only 10% of the total mass. In terms of (V max , r Vmax ), the total luminosity of a halo is given by where f is an O(1) numerical factor that depends on the shape of the density profile; for an NFW profile f = 1.227, and for an α = 0.17 Einasto profile f = 1.735. In physical units, the total annihilation luminosity is for NFW (top) and α = 0.17 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 (V max , r Vmax ). 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 r s and a cross-over point at ∼ 10 −3 r s .

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 mass-to-light ratios of up to several hundred [87,88,89,90], 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, for example).
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,5,6,7,8,9,10,11,12], raising the total from the 9 "classical" ones to 23. Many of the newly discovered satellites are so-called "ultra-faint" dSph's, with luminosities as low as 1, 000L ⊙ 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 [91] and a radial distribution of DM subhalos motivated by numerical simulations, this estimate can grow to several hundreds of satellite galaxies in total [92,93].
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 V max , of the DM halo hosting the galaxies. Owing to the extreme faintness of these objects and their lack of a detectable gaseous component [94], 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 line-of-sight velocity dispersions have confirmed that these objects are in fact strongly DM dominated [87,88,89,90]. 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 M 0.3 ≡ M (< 0.3 kpc) for 18 of the Milky Way dSph's, and found that, surprisingly, they all have M 0.3 ≈ 10 7 M ⊙ to within a factor of two [31].
State-of-the-art 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,70] and the Aquarius Project [14], have now reached an adequate mass and force resolution to directly determine M 0.3 in their simulated subhalos. This makes it possible to infer the most likely values of (V max , r Vmax ) for a Milky Way dSph of a given M 0.3 and Galacto-centric distance D, by identifying all simulated subhalos with comparable M 0.3 and D and averaging over their (V max , r Vmax ). 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 [31] and with the more recent and higher resolution Via Lactea II (VL2) simulation.
We randomly generated 100 observer locations at 8 kpc from the VL2 host halo center, and selected, for each Milky Way dSph in [31] separately, all simulated subhalos with distances within 40% and numerically determined M 0.3 within the published 1 − σ error bars. We then determined the median value and the 16 th and 84 th percentiles of (V max , r Vmax ) and (V peak , r Vpeak ) for each dSph. These values are given in Table I. The median values of V max range from 5.0 km s −1 (Leo IV) to 19 km s −1 (Draco, Leo I), and of V peak from 6.7 km s −1 (Leo IV) to 30 km s −1 (Ursa Minor, Carina). Note that, as expected, dSph's closer to the Galactic Center typically show a larger reduction from V peak to V max , 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 Eq. (9) for an NFW profile to calculate the total luminosity L tot for every simulated subhalo. Additionally we also determine L 0.3 , 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: F tot = L tot /4πD 2 , the total expected flux from the dSph, and F c , the flux from a central region subtending 0.15 • , which is comparable to the angular resolution of Fermi above 3 GeV. F c thus corresponds to the brightest "pixel" in a Fermi γ-ray image of a subhalo. These numbers are given in Table II.

Current observational constraints
Several ACT have performed observations of a handful of dSph's.
• The H.E.S.S. array (consisting of four 107 m 2 telescopes with a 5 • field of view and an energy threshold of 160 GeV [95]) has obtained an 11 hour exposure of the Sagittarius dwarf galaxy. No γ-ray signal was detected, resulting in a flux limit of 3.6 × 10 −12 cm −2 s −1 (95% confidence) at E > 250 GeV, and a corresponding limit on the cross section of σv < ∼ 10 −23 cm 3 s −1 for an NFW profile and σv < ∼ 2 × 10 −25 cm 3 s −1 for a cored profile (for a m χ = 100 GeV − 1 TeV neutralino) [44]. Note that   the Sagittarius dwarf is undergoing a strong tidal interaction with the Milky Way galaxy [96], and no confident determination of M 0.3 has been possible.
• The VERITAS array (consisting of four 144 m 2 telescopes with a 3.5 • field of view and an energy threshold of 100 GeV [97]) has conducted a 15 hours observation of Willman 1 and 20 hour observations each of Draco and Ursa Minor [45]. No γ-ray signal was detected at a flux limit of ∼ 1% of the flux from the Crab Nebula, corresponding to a limit of 2.4 × 10 −12 cm −2 s −1 (95% confidence) at E > 200 GeV [98].
• Additionally the MAGIC [46,47] and STACEE [48] telescopes have reported observations of Willman 1 and Draco, resulting in comparable or slightly higher flux limits.
To convert the values of F c in Table II 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 Eq. (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 [71,98,99,100,101]. 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.

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 all-sky 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 [102] notwithstanding, it is at present not yet possible to perform realistic cosmological hydrodynamic galaxy-formation simulations, which include, in addition to the DM dynamics, all the relevant baryonic physics of gas cooling, star formation, supernova and AGN feedback, etc. 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 DM-only 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 baryon-dominated 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 mass-tolight 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 4, 100 M ⊙ and a force softening of 40 pc, the VL2 simulation resolves over 50,000 subhalos today within the host's r 200 = 402 kpc (the radius enclosing an average density 200 times the mean matter value). Above ∼ 200 particles per halo, the differential subhalo mass function is well-fit by a single power law, dN/dM ∼ M −1.9 , and the cumulative V max function is N (> V max ) ∼ V −3 max [13]. The radial distribution of subhalos is "anti-biased" 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,85]. Similar results have been obtained by the Aquarius group [14,83]. Typical subhalo concentrations, defined as ∆ V = ρ(< r Vmax ) /ρ crit , grow towards the center, owing to a combination of earlier formation times [103,104] and stronger tidal stripping of central subhalos: VL2 subhalos on average have a 60 times higher ∆ V at 8 kpc than at 400 kpc [13]. Note that this also implies ∼ 7 times higher annihilation luminosities for central subhalos, since The counter-acting 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 [78]. 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 5σ 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 Eq. (9) and assuming an NFW density profile. The luminosities would be ∼ 40% higher if an Einasto (α = 0.17) profile had been adopted instead. We then converted these luminosities to fluxes by dividing by 4πD 2 , where the distances D were determined for 100 randomly chosen observer locations 8 kpc from the host halo center. The resulting values of F tot are plotted in Figure 3, for just one of the 100 observer positions, as a function of the subhalos' M 0.3 , V max , and V peak . 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' V max against their distance to the host halo centerD. There appears to be a dearth of the lowest V max subhalos (V max < ∼ 2 km s −1 ) at small distances, but at the moment it is not clear whether this suppression is a real effect or a numerical artifact. It's 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 V max by noting that the luminosity scales as L ∼ V 3 max √ ∆ V and the typical distance as The typical flux should thus scale as F ∼ L/D 2 ∼ V 1/3 max √ ∆ V , and would be higher for more massive subhalos at a fixed ∆ V . Actually lower V max subhalos might be expected to have higher ∆ V due to their earlier formation times, but it remains to be seen to what degree this expectation is borne out in numerical simulations.
The points in Figure 3 can be directly compared with the values for the known Milky Way dSph's in Tables I  and II: 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 F tot and F c . 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 F tot greater than that of (Carina, Draco, Ursa Minor, Ursa Major, Segue 1) is (90,54,43,17,13), and the corresponding numbers for F c are (96,62,49,24,19). 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 F tot and F c 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 M 0.3 = 5 × 10 6 M ⊙ , V max = 8 km s −1 , and V max = 14 km s −1 . 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 I only Leo IV falls below these limits. In the bottom panel of Figure 5 we plot f dark (> F tot ), the fraction of subhalos without stars, as a function of the limiting annihilation flux F tot . f dark falls monotonically with F tot , 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.

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 cut-off in the matter power spectrum, set by collisional damping and free streaming in the early universe [105,106]. For typical WIMP DM, this cut-off occurs at masses of m 0 = 10 −12 to 10 −4 M ⊙ [25,26], some 10 to 20 orders of magnitude below VL2's mass resolution. Since the annihilation rate goes as ρ 2 and ρ 2 > ρ 2 , 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 so-called substructure boost factor depends sensitively on the properties of subhalos below the simulation's resolution limit, in particular on the behavior of the concentration-mass 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 [52]. More sophisticated (semi-)analytical models, accounting for different possible continuations of the concentration-mass relation to lower masses, typically find smaller boosts of around a few tens [72,78,100,107].
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 [78]), owing to the ρ(r) 2 de-pendence of the annihilation rate. The luminosity contribution from a subhalo's sub-substructure population (i.e. its boost), however, is much less centrally concentrated: at best it follows the subhalo's mass density profile ρ(r), although it might very well even be anti-biased. 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 doesn't contribute much to the detection significance. In other words, the boost factor might apply to F tot , but much less (or not at all) to F c ; yet it is F c 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 2 . 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 [108], 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 [49,50,51,52,72,78].

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 D and central masses M 0.3 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 V max and r Vmax , and hence for the annihilation luminosity L and flux F 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 as-of-yet 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 don't 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 σ 8 and n s . Other work has found considerable halo-to-halo scatter [14,109,110], 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 Fellowship at the Institute for Advanced Study. I would like to thank my collaborators from the Via Lactea Project for their expertise and invaluable contributions.