the existence of nonbaryonic dark matter in modern cosmology is strongly supported by several independent signatures  @xcite , including the cosmic microwave background ( cmb ) , gravitational lensing and large scale structure surveys .
however , the nature of dark matter remains a mystery . to reproduce the correct relic density
, it is naturally assumed that dark matter is composed of weakly interacting massive particles ( wimps ) such as the supersymmetric neutralino which is one of the most popular candidates  @xcite . since the self - annihilation rate is proportional to the dark matter number density squared , potential signals from dark matter annihilation are most likely to be detected in highly dense regions such as the centers of galaxies .
it is also expected that the energy released in dark matter annihilation in the early universe can influence cmb anisotropies . with the high precision cmb data from the wmap satellite  @xcite , strong constraints on the properties of dark matter
can be obtained  @xcite .
a multi - wavelength analysis of dark matter annihilations has been proposed recently in ref .  @xcite .
recently , several cosmic ray and gamma - ray experiments provided new windows to detect the signatures of dark matter . in particular , the pamela satellite  @xcite reported an `` excess '' of the positron fraction above 10 gev .
also the @xmath1 spectra above a few hundred gev from the fermi - lat  @xcite and hess measurements  @xcite are significantly harder than expected  @xcite .
these striking observations were often interpreted in terms of dark matter annihilation or decay which can contribute to the flux of high energy electrons and positrons  @xcite . however , standard astrophysical sources , such as pulsars  @xcite and supernova remnants  @xcite also provide a possible explanation for these excesses .
the high energy electrons and positrons induced by wimp annihilation produce gamma - rays through inverse compton scattering ( ics ) off the low energy background photons in the interstellar radiation field ( isrf ) and through bremsstrahlung emissions due to the interaction with the ionized interstellar medium .
these components add to the prompt gamma - ray radiation accompanying pairs emitted in the annihilation event  @xcite .
gamma - ray dark matter signatures could be identified by the fermi large area telescope  @xcite , due to its unprecedented angular and energy resolution , despite the fact that the gamma - ray spectrum is dominated by conventional astrophysical sources , such as pulsars , supernovas , blazars  @xcite and perhaps structure formation shocks  @xcite .
many efforts  @xcite have been devoted to extract indirect dark matter signatures or to set limits on a large class of wimp models by investigating the fermi - lat diffuse gamma - ray measurements .
especially , a universal model - independent method can be used to place constraints on any dark matter model  @xcite through convolving _ response functions _ of signal - to - background with a specific injection energy spectrum of electrons and positrons .    in the wimp
cold dark matter scenario , the dark halos can form at very high redshift , @xmath2 , with a minimum mass of @xmath3  @xcite determined by the free - streaming limit and collisional damping leading to a cutoff of the primordial power spectrum .
this suggests that if they can survive until the present day , an enormous number of dark matter clumps ( subhalos ) are expected to be embedded in our galaxy .
recent numerical simulations  @xcite confirm this prediction , although the role of the tidal effects from the baryonic component has yet to be quantified and there is no direct test by observations .
nevertheless , several references have discussed the resulting signals pointing out that subhalos can boost annihilation rates and produce a distinct radial distribution of emissions  @xcite .
alternatively , a statistical analysis of the full - sky emission map can be used to identify dark matter signatures , which can be promising if , e.g. , one can identify features in the anisotropy power spectrum that characterize uniquely the dark matter spatial distribution with respect to other astrophysical sources . in this respect ,
one is interested in quantifying the anisotropy signatures from annihilating dark matter which can be different from those of astrophysical origin particularly at small scales .
the first calculations of this kind were performed analytically for gamma - ray background anisotropies induced by both annihilating extragalactic dark matter halos  @xcite as well as galactic subhalos  @xcite .
analytic calculations were also performed for the anisotropic radio signatures from decaying dark matter  @xcite and compared with the mean intensity and angular power spectrum of the astrophysical and cosmological radio background .
more recently anisotropies of the diffuse gamma - ray background from annihilating dark matter in galactic subhalos have been computed using a numerical approach in ref .
@xcite . according to their results
, fermi - lat should be sensitive enough to constrain the amount of dark matter substructure through the anisotropy signal discussed above .
the same numerical approach as in  @xcite can be extended to the anisotropies from both the galactic and extragalactic dark matter contributions  @xcite . in general , however , both the mean flux and anisotropy signal are contaminated by ordinary astrophysical sources including galactic and extragalactic resolved point sources , structure formation shocks  @xcite , and normal galaxies  @xcite .
removal of such contributions  @xcite , can greatly improve the detectability for searches of dark matter .
all of the above analyses were based on the gamma - rays from direct annihilation into photons .
however , the ics radiation from energetic electron and positron pairs produced by annihilation of dark matter also generates unique gamma - ray signatures , in an energy range that depends on the dark matter model .
thus in the present paper we study the gamma - ray anisotropies produced by this `` secondary '' emission , which can be regarded as a new target for indirect searches of signatures induced by galactic substructure halos .
the important point is that the properties of the angular anisotropies for emission from secondary electrons and positrons will be different in general than for prompt emission , because the former will be affected by propagation effects in the galaxy .
in fact , by adopting a spatial diffusion model typical for high energy electrons and positrons in the galaxy , we show that the angular power spectrum is suppressed at small angular scales corresponding to the distance diffusively traveled by the charged particles during their energy loss time . as a result , for a typical dark matter model with mass of 1 tev and a canonical thermal freeze - out cross section @xmath4 , the angular power spectrum from inverse compton scattering peaks at large angular scales .
the remainder of this paper is organized as follows . in sect .
[ sec : formulation ] , we set up the formalism used in this paper for calculating the mean intensity and angular power spectrum of gamma - rays induced by dark matter annihilation and develop a robust numerical scheme to simulate the small subhalos with masses down to @xmath5 .
we present the results of our calculations in sect .
[ sec : result ] . finally , we summary our paper in sect .  [
sec : conclusion ] .
we will use natural units in which @xmath6 throughout .
the propagation of high energy electrons and positrons in the turbulent galactic magnetic field can be described as a diffusion - energy loss equation  @xcite for the electron - positron number density @xmath7 , neglecting the convection and re - acceleration terms which are only relevant for electrons and positrons below 10 gev  @xcite : @xmath8 here , the spatial diffusion coefficient is taken as @xmath9 ; @xmath10 is the energy loss term and @xmath11 the source term , both of which are described in more details below . in the following
we will often use the term electron as a short hand for both electrons and positrons . with the assumption that spatial diffusion and energy loss coefficients are spatially independent , for a given source distribution and boundary conditions , the propagation equation can be solved analytically .
we assume the diffusion zone to be a cylinder with half - height @xmath12 of a few kpc and radius @xmath13kpc . in this paper
, we adopt the widely used med model which , compared to other models , i.e. the min and max models , predicts intermediate values for the antiproton flux when fitted to reproduce the b / c ratio  @xcite with @xmath14 , @xmath15 , and the half - height of the diffusion zone @xmath16kpc .
the min and max propagation models would decrease and increase the predicted gamma - ray fluxes roughly by a factor of three and two , respectively .
however , we verified that the shape of the angular power spectrum predicted by these three propagation models is basically the same .
since in the present paper we mainly focus on a new `` suppression '' feature in angular power spectrum of the gamma - ray component from ics , we will use the med model in the following calculations .    in this study , we investigate only the case in which electrons and positrons can escape freely at the boundaries of the diffusion zone . if the diffusion coefficient increases exponentially instead of the free - escape boundary condition , we expect that the angular power spectrum would be modified only at multipoles @xmath17 since the diffusion of electrons and positrons would flatten their distribution at length scales even larger than the whole diffusion zone .
this implies that the angular power spectrum produced in the region of extremely large diffusion coefficient can not significantly influence the signature in the full - sky angular power spectrum we found at @xmath18 . in other words ,
the anisotropies weakly depend on different boundary conditions .      at energies above 10 gev ,
the dominant energy losses are synchrotron radiation and inverse compton scattering ( ics ) on the interstellar radiation field .
therefore , in the thomson limit we write @xmath19 [ for a more complicated treatment of energy loss see ref.@xcite ] , where @xmath20 for starlight ( sl ) , infrared ( ir ) , cmb photons and a magnetic field of @xmath21 . the inter - stellar radiation field ( isrf )
can be approximately characterized as a superposition of three blackbody - like spectra with different temperatures and normalization factors relative to a true black - body emitter : one for the cmb with @xmath22 , for the ir with @xmath23 and for the sl with @xmath24  @xcite .
the typical normalization of the sl and ir fields of radiation depends on the position in the galaxy .
the averaged normalizations of isrf photon densities per energy used in this paper are @xmath25 and @xmath26 for sl , ir and cmb , respectively  @xcite .
although these normalizations are valid only in the region with latitude @xmath27 , we have checked that changing the normalization of the isrf affects the gamma - ray emissions only weakly because an increased emission is partly compensated by a decrease of the density of electrons and positrons as energy losses increase .
we also verified that the deviation from the detailed numerical simulation with galprop  @xcite is less than a factor of two for a realistic spatial distribution of the isrf  @xcite .      in the annihilating dark matter scenario ,
the source term can be written as @xmath28 where @xmath29 is the annihilation spectrum into electrons and positrons at energy @xmath30 . for simplicity , we assume mono - energetic injection of the positron and electron in case of cp conservation , namely @xmath31 . we choose
@xmath32 , which can well fit the pamela excess  @xcite while not in conflict with gamma - ray observations by fermi - lat  @xcite
. finally , we use @xmath33 to reproduce the correct relic density for thermal freeze - out . by convolving our results for the gamma - ray spectra with the pair energy ,
our computational approach can be easily adapted to pair spectra different from mono - energetic injection , such as from dark matter annihilating into @xmath34 , and @xmath35 .
as most of the electrons escape towards the @xmath36direction , we impose the dirichlet boundary condition @xmath37 , at which the particles can freely escape .
we thus model the diffusion zone as an infinite slab of half thickness @xmath12 , and we take @xmath38 for the med diffusion model .
the free - space green s function  @xcite for eq .
( [ eq : diffeqn ] ) is @xmath39= \frac{1}{b(e)}\frac{1}{(\pi\lambda_d^2)^{3/2}}\e^{-({\bf r}-{\bf r}^\prime)^2/\lambda_d^2}\,,\ ] ] where we have defined the diffusion length as @xmath40 which is the average distance @xmath41 diffuse through during their energy loss time .
then , the green s function satisfying appropriate boundary conditions can be obtained by considering a series of image charges at positions @xmath42 .
one can verify that @xmath43 fulfills the dirichlet boundary condition .
thus , the general solution to eq .
( [ eq : diffeqn ] ) , in the limit of time - independent sources and electron / positron number densities which already reached equilibrium , is given by @xmath44    for the diffuse gamma - ray emission we are more interested in the column density of electrons , @xmath45 than in the local space density of electrons .
the observer is located at the solar system . in galactic coordinates ,
a point in cartesian coordinates @xmath46 at a distance @xmath47 from the observer is given by @xmath48 where @xmath49 and @xmath50 are the galactic longitude and latitude , respectively .
we truncate the integral in eq .
( [ eq : column ] ) at the edge of the diffusion zone , beyond which particles are not confined , @xmath51 or @xmath52 .
the line - of - sight integral can directly act on the free green s function  @xcite , @xmath53/\lambda_d}}{2\pi \lambda_d^2b(e ) } \left[\erf\left(\frac{\ell_{\rm max}-{\bf n}\cdot{\bf r}^\prime}{\lambda_d}\right)- \erf\left(\frac{-{\bf n}\cdot{\bf r}^\prime}{\lambda_d}\right)\right]\,,\nonumber\end{aligned}\ ] ] where we have defined the unit - vector @xmath54 . the green s function satisfying
the boundary condition is thus @xmath55 the column density of electrons therefore reads @xmath56 q({\bf r}^\prime , e^\prime)\,.\ ] ] in the limit of @xmath57 , where @xmath58 is the scale radius of the subhalo profile , the subhalo can be regarded as a point - like source .
(  [ eq : sigma ] ) can then be simplified to @xmath59j_{k}({\bf r}_k),\ ] ] where @xmath60 for a given subhalo source @xmath61 located at @xmath62 .
for the largest subhalos with masses larger than @xmath63 , their radius @xmath58 can be somewhat larger than the diffusion length @xmath64 .
nevertheless , their contribution to the total flux is a factor @xmath65 smaller than the flux from the smaller subhalos .
therefore , neither the mean intensity nor the dimensional angular power spectrum ( see sect .  [
subsubsec : cl ] ) relies significantly on the distribution of electrons in the largest subhalos . as a result , for our purposes we can apply eq .
( [ eq : ssigma ] ) to all subhalos even if @xmath66 .      alternatively , for primary electrons and positrons from the smooth host dark matter halo , the bessel - fourier scheme  @xcite can require less computational time than the green s function .
the electron and positron number density after propagation can be expressed as @xmath67 and @xmath68 is the halo function defined by @xmath69e^{-\left[\left(\frac{n\pi}{2l}\right)^2+\frac{\alpha_i^2}{r^2}\right]\frac{\lambda_d^2}{4 } } r_{i , n}\,.\ ] ] here , the @xmath70 are the zeros of the bessel function @xmath71 and @xmath72 are the coefficients of the bessel - fourier transform of the source term .
for relativistic electrons and positrons with energy @xmath30 up scattering background photons from energy @xmath73 to @xmath74 , the emitted inverse compton power per energy interval is @xmath75 where @xmath76 is the differential isrf photon number density , while the differential cross section @xmath77 is given by the klein - nishina formula  @xcite .
folding @xmath78 with the spectral distribution of the equilibrium number density of electrons and positrons , we get the emissivity of ic photons of energy @xmath74 , @xmath79 which yields the ic intensity at energy @xmath74 by the line - of - sight integral @xmath80    according to the eq .
( [ eq : sigma ] ) , the ic intensity from the electrons and positrons can be simplified to @xmath81 there is a well known `` delta - function approximation '' where an electron with energy @xmath30 inverse compton scattering black - body photons with temperature @xmath82 emits photons with a characteristic energy @xmath74 in both the thomson and extreme klein - nishina limits  @xcite , @xmath83 , @xmath84 $ ] , where @xmath85 is the total ic energy loss rate of the electron .
the numerical calculations show that @xmath86 may be approximated by @xmath87 and in the thomson regime one has @xmath88 , where @xmath89 is the electron charge and @xmath90 is the background photon energy density .
( [ eq : ici ] ) can thus be simplified to @xmath91 where the electron / positron energy @xmath30 is related to the gamma - ray energy @xmath74 through the following condition : @xmath92 .
this relation reproduces the known slope @xmath93 for an electron spectrum of @xmath94 . in the thomson limit and for our choice of monoenergetic injection of electron - positron pairs , the index @xmath95 in the stationary situation if the energy loss term dominates on the right hand side of eq .
( [ eq : diffeqn ] ) , as is usually the case for electron energies above 10 gev . for @xmath96tev
we are always in the thomson limit .      in the present section
we discuss contributions to the source term @xmath97 both from the smooth host halo and the individual subhalos .
the two most important features of the subhalos are their mass distribution and their spatial distribution . in the following ,
we briefly summarize the properties of subhalos used in our calculations for which we adopt the same description as in ref .
@xcite , and references therein .
there are two widely used scenarios for describing the subhalo radial distribution .
one is unbiased relative to the galactic smooth component ( host halo ) with the nfw density profile  @xcite given by @xmath98 with @xmath99 , where @xmath58 is a scale radius and @xmath100 is the characteristic density . for the case where the subhalo distribution is anti - biased compared to the smooth component , we use the fitting formula of the subhalo radial distribution from gao _ et al _  @xcite .
the cumulative fraction of subhalos within a given radius is @xmath101 with @xmath102 , @xmath103 , @xmath104 , @xmath105 and @xmath106 is the host halo concentration .
furthermore , @xmath107 is the total number of subhalos within virial radius @xmath108 of the host halo . for a host halo with the nfw density profile we adopted @xmath109 , @xmath110kpc and @xmath111  @xcite .
recent simulations suggest that the cumulative number of subhalos above a given mass @xmath112 in units of the solar mass can be fitted simply by a power law  @xcite , @xmath113 in the present paper , we simulate the subhalos with mass down to around one earth mass , @xmath114 and choose @xmath115 , as suggested by simulations .
we find that choosing @xmath116 between @xmath117 and @xmath26  @xcite would change the mean intensity by roughly a factor of @xmath118 .
more recent simulations suggest that the central structure of dark matter halos deviates from the nfw profile in the innermost regions , which can be well reproduced by an einasto density profile .
therefore , here we use the einasto profile with the parameter @xmath119 given by gao _
@xcite , @xmath120\,.\ ] ] the relation of halo concentration and mass is given by @xmath121 the slope and normalization are consistent with those found by bullock _
@xcite .
the angular power spectrum of the emission maps can be calculated by using the public healpix package  @xcite .
we define the dimensionless quantity @xmath122 as a function on the sphere , which can be expanded in spherical harmonics @xmath123 as @xmath124 where i(@xmath125 ) describes the gamma - ray intensity in the direction @xmath125 .
the _ dimensionless _ angular power spectrum of @xmath126 is given by the coefficients @xmath127 note that the measured dimensional amplitude @xmath128 of the total angular power spectrum in units of intensity squared can be obtained by multiplying the dimensionless angular power spectrum of a given component , @xmath129 , with its mean intensity squared , @xmath130 and summing over the components , @xmath131 where in our case the sum basically runs over the contributions of the smooth host halo , the subhalo distribution and other astrophysical foregrounds and backgrounds .
here we assumed that the different contributions are uncorrelated
.      in numerical calculations , taking into account subhalos down to masses @xmath5 in the angular power spectrum would require the generation of @xmath132 subhalos within the diffusion zone in a given monte carlo realization .
this is not very practical and severely limits the number of realizations one can simulate .
thus , we need to develop a sound scheme to circumvent these computational limits . here
we take advantage of a few simple facts : ( i ) two terms contribute to the fluctuation anisotropies in units of average intensity squared ; the one - subhalo term ( @xmath133 ) which is poissonian noise and the two - subhalo term ( @xmath134 ) due to the cross - correlation of substructure arising from their radial distribution within the host halo ( see eqs .  @xmath135 and @xmath136 of ref .
@xcite ) ) .
( ii ) qualitatively speaking , the two - subhalo term scales linearly with the square of the average flux whereas the one - subhalo term is proportional to an integral of squared luminosities of subhalos . since subhalo luminosity @xmath137 is related to its mass @xmath112 via @xmath138 , by using eq .
( [ eq : c ] ) , the correlation between gamma - ray luminosity and mass thus is positive , i.e. , @xmath139 .
one - subhalo term and two - subhalo term per decade of subhalo mass then can be estimated as @xmath140 and @xmath141 .
therefore , large and relatively rare substructures dominate the one - subhalo term , but negligibly to the two - subhalo term , and the small but numerous substructure halos dominate the two - subhalo term .
there is thus a dividing mass of substructures , @xmath142 , below which their contribution to one - subhalo term are negligible .
the numerical value of this mass depends on the assumed radial distribution of the substructure halos ( see below ) . in any case
, it turns out that halos with @xmath143 have a characteristic radius , @xmath58 , much smaller than the diffusion length of the gamma - ray emitting electrons , i.e. @xmath144 .
thus , effectively all halos below @xmath142 have the same image on the sky in ics with a characteristic size given by @xmath64 .
this further simplifies our calculations because besides minor poissonian noise , the angular power spectrum from subhalos with @xmath143 is simply obtained by convoluting the subhalo spatial distribution with the mass - independent ics subhalo image .
thus the contribution from subhalos with @xmath143 is taken into account straightforwardly as follows : we generate a distribution of subhalos with mass in the decade @xmath145 to obtain the full - sky map .
the number of generated subhalos is typically much smaller than the actual value , @xmath146 , but is sufficiently large that the two - subhalo term divided by the mean intensity squared has converged .
the intensity of the map and the corresponding two - subhalo term thus can be rescaled to the value appropriate for @xmath147 halos .
in this way , the simulated two - subhalo term of the angular power spectrum can fairly approach the actual one .    since the one - subhalo term is not linearly proportional to the square of the average intensity and
is contributed mostly by larger subhalos , we thus complete the calculation by adding the contribution of subhalos with mass @xmath148 , which can now be done with a direct monte carlo simulation .
specifically , we find that our calculation of the angular power spectrum reaches convergence when we use at least @xmath149 subhalos for each mass decade from @xmath150 to @xmath142 within the diffusion zone , when we use @xmath151 for the unbiased radial distribution and @xmath152 for the anti - biased case .
note that this method can not be applied to the angular power spectrum of the direct annihilation component since in this case the profile of the emitting region depends on the profile of each subhalo rather than on the identical diffusion length which just depends on the energy of electrons and positrons for all subhalo masses .
in this section , we make use of the formalism derived in the previous section to compute the gamma - ray mean intensity map and associated anisotropy power spectrum due to both the smooth host halo component and the substructure halos . we also study how the map morphology and anisotropies depend on the radial distribution of subhalos and of the gamma - ray energy .
we first consider the gamma - ray emission produced by dark matter annihilating into @xmath41 pairs scattering off the isrf both in the smooth host halo and the subhalos .
using the formalism developed in sect .
[ sec : formulation ] , the fluxes from the smooth host halo at 1 gev , 10 gev and 100 gev are shown in fig .  [
fig : diff - host ] .
here we use galactic coordinates , where an observer is located at 8.5 kpc from the galactic center .
we adopt the nfw density profile with the parameters given in sect .
[ subsec : subst ] , and find , as expected , that most of the signal comes from the central galactic region where the dark matter is highly concentrated .
the values of the mean gamma - ray intensity from both the host halo and the substructure are also reported in tab .  [
tab : i ] for the above three photon energies .
we find that in both cases the mean gamma - ray intensity decreases faster than @xmath153 in eq .
( [ eq : simici ] ) .
this is because there will be a high energy cut - off in @xmath74 once the required energy of the parent pairs @xmath154^{1/2}$ ] scattering off the low energy background photons exceeds the maximum energy @xmath155 produced by dark matter annihilation .
the effect is even stronger when the delta - function approximation for ics is assumed , as in our case . as a consequence ,
all background photon energies at cmb energies and above contribute to the 1 gev ics photon flux , whereas only the ir and sl contribute to the 10 gev photon flux and only the sl contributes to the 100 gev photon flux , as also summarized in tab .
[ tab : e_e ] .
furthermore , using the relation between @xmath30 and @xmath74 and the expression for @xmath85 , for an electron spectrum @xmath156 eq .
( [ eq : simici ] ) gives the scaling @xmath157 .
using the different normalizations for the cmb , ir and sl densities , then would suggest that the gamma - ray intensity originating from @xmath41 pairs scattering off cmb photons is about 5 times larger than the one from scattering off the ir photons and about 10 times lager than the contribution from scattering off sl photons , provided there is no restriction from the kinematics .
however , the @xmath41 pairs produced by annihilating dark matter do have an absolute cutoff at the parent particle energy . as a result ,
the more detailed numerical results show that most of the gamma - ray intensity at 1 gev is produced by pairs with @xmath158 scattering off the cmb , whereas pairs with @xmath159 scattering off ir photons dominate the 10 gev gamma - ray flux and pairs with @xmath160 scattering off sl dominate ics photons at 100 gev .
.the mean gamma - ray intensities @xmath161 at 1 gev , 10 gev and 100 gev , averaged over the sky , from subhalos with minimum mass @xmath162 for the unbiased and anti - biased radial distribution and from the smooth host halo . the unit is @xmath163 . [ cols="^,^,^,^",options="header " , ]     in fig .
[ fig : diff - unbias ] , we present the sky map of gamma - ray emission predicted by an unbiased radial distribution of subhalos : in order to avoid saturation from the dominant population of small halos , in the figure we only include halos with mass above @xmath164 . the most remarkable feature in fig .
[ fig : diff - unbias ] is that the diffuse emission regions can extend to a few kpc in length scale , corresponding to a few tens of degrees on the sky .
furthermore , the size of the brightest regions representing the largest intensity tends to increase with increasing @xmath74 , which could give rise to an increase of the angular power spectrum on corresponding angular scales .
the size of the emitting region is basically determined by the diffusion length @xmath64 which is the distance that the @xmath41 diffuses during their energy loss time .
this length scale can be estimated from eq .
( [ eq : d ] ) , giving @xmath165 for the med propagation model .
this energy dependence of @xmath64 is shown in tab .  [
tab : e_e ] .
the typical diffusion length is @xmath166 kpc , @xmath167 kpc and @xmath168 kpc for electron energies leading to emission at @xmath169gev , 10 gev and 100 gev , respectively , corresponding to an angular scale @xmath170 with @xmath171 few kpc the typical distance to the dark matter annihilation source .
these angular scales are roughly what one sees in fig .
[ fig : diff - unbias ] and the above estimate also applies to the smooth host halo case ( see  fig .  [
fig : diff - host ] ) and to the anti - biased radial distribution of subhalos .
how does the mean intensity depend on @xmath172 ?
empirically , we find that the mean intensity roughly doubles with each decade of decreasing mass of subhalos , similar to what has been found by ref .  @xcite . in light of the subhalo mass function @xmath173 in eq .
( [ eq : dndm ] ) and the concentration parameter @xmath174 in eq .
( [ eq : c ] ) , the annihilation rates per decade of subhalo mass can be approximated by @xmath175 , which is fairly consistent with our detailed numerical calculation .    in order to determine the influence of the radial distribution of subhalos on the intensity , in fig .
[ fig : diff - bias ] the gamma - ray full - sky maps at 1 gev , 10 gev and 100 gev for the un - biased distribution are compared with the anti - biased distribution . for the anti - biased radial distribution the mean
intensities are roughly 50 times smaller than for the unbiased case since the mean distance of subhalos to us is much larger than in the unbiased case .
one notices the important feature in fig .
[ fig : diff - bias ] that the emission from the spatially anti - biased distribution is much less centrally concentrated , and apparently accumulates around the galactic plane compared with the unbiased case .
this can be understood from the fact that only electrons and positrons confined within the diffusion zone of scale height @xmath176 kpc and radius @xmath177 kpc for the med model can efficiently produce gamma - rays by ics and the subhalos within the diffusion zone are distributed much more isotropically for the anti - biased case compared to the more central distribution of the unbiased case .
furthermore , ics outside the diffusion zone contributes less than @xmath178 to the mean intensity  @xcite .
the dimensionless angular power spectra @xmath179 of gamma - ray emission due to the ics at photon energies 1 gev , 10 gev and 100 gev from the main components of galactic dark matter are presented in fig .
[ fig : cl ] . shown
are the power spectrum for subhalos with @xmath180 for the unbiased radial distribution and for the smooth host halo . to clearly illustrate the effects of diffusion , the angular power spectra for spatially unbiased subhalos with @xmath181 and for the host halo
are also shown in absence of diffusion .
we recall that in order to obtain for a given component the contribution to the angular power spectrum @xmath128 in units of intensity squared , according to eq .
( [ c_l ] ) one has to multiply the dimensionless power spectra @xmath179 shown in fig .
[ fig : cl ] by the squared total intensity @xmath182 of the corresponding component from tab .
[ tab : i ] .
[ fig : cl ] demonstrates the remarkable feature that the power spectrum of the ics component of galactic dark matter annihilation is exponentially suppressed for @xmath183 compared with what one would obtain without diffusion .
furthermore , the lower energy gamma - rays have more angular power at @xmath183 corresponding to small angular scale .
this can be understood from the energy - dependence of the diffusion length : intensity fluctuations should be damped on length scales smaller than the diffusion length @xmath64 , corresponding to a multipole @xmath184 , where @xmath185 is the typical distance to the dark matter annihilation source .
based on the discussion in sect.[subsec : dg ] , we can estimate the diffusion length of the electrons and positrons emitting a given gamma - ray energy .
when doing so , one should keep in mind that the @xmath41 pairs interact with three different photon backgrounds and that , as it turns out , unlike the case of a single background , the lower the gamma - ray energy the higher the energy of the emitting electrons . since we found that the diffusion length decreases with the electron energy , @xmath165 , we expect suppression of anisotropies to occur at smaller scales , or larger values of the multipole @xmath49 , for lower photon energy . for other propagation models such as the min and max models which have a slightly different energy dependence of the spatial diffusion coefficient
, the diffusion length would be slightly larger or smaller , respectively .
this would shift the suppression scale in the angular power spectrum by less than a factor of two in the multipole @xmath49 .
in fact , for the smooth host halo the typical distance is @xmath186 kpc , and the suppression due to diffusion should occur at @xmath187 and @xmath188 for gamma - ray energies of 1 gev , 10 gev and 100 gev , respectively .
this analysis can be applied to the case of subhalos . for an unbiased spatial distribution of subhalos ,
the typical distance to a subhalo is @xmath189kpc which corresponds to diffusive suppression at relatively smaller @xmath49 compared to the host halo case .
these simple estimates are consistent with our detailed calculations shown in fig .
[ fig : cl ] .
we also find that the amplitude of the dimensionless angular power spectrum @xmath179 from the smooth host halo is larger than that from the subhalos since the emissivity profile from annihilation in the smooth halo is proportional to the density squared and thus more peaked than annihilation in the subhalos which essentially follow the linear density profile of an nfw profile , as we see in fig .
[ fig : cl ] . we note that although the smooth host halo has a large amplitude of the dimensionless angular power spectrum @xmath179 , its contribution to the total intensity is as small as @xmath190 for the unbiased subhalo distribution and @xmath191 for the anti - biased distribution , as seen in tab .
[ tab : i ] .
we can therefore safely neglect the contribution from the host halo both to the mean intensity and to the dimensional angular power spectrum @xmath128 .
finally , we show in fig .
[ fig : cl ] the dimensionless angular power spectrum @xmath179 for a smooth halo with emissivity tracing the density of the nfw profile instead of the squared density that would be relevant for the contribution of the host halo in decay scenarios .
this shows that a smooth nfw profile describes the emission profile by a large number of subhalos following an unbiased radial distribution very well , at a level of @xmath190 .
this is not surprising since the number of subhalos within the diffusion zone with masses below @xmath192 is sufficiently large , @xmath193 , for each mass decade to strongly suppress any deviation from a smooth distribution .
this conclusion is also true for the anti - biased case . despite the fact that the contribution of large subhalos in the mass range of @xmath194 fluctuates strongly , their contribution to the total emission are three orders of magnitude smaller which leads to fluctuations at @xmath195 level in angular power spectrum as seen in fig .
[ fig : cl ] .
how does the radial subhalo distribution affect the angular power spectrum ?
in fig .  [
fig : bias_cl ] we compare results for the unbiased case and the anti - biased case . at small @xmath49 the angular power spectrum @xmath179 for the anti - biased distribution is suppressed relative to the unbiased case which is due to the more isotropic subhalo distribution at large angular scales seen in fig .
[ fig : diff - bias ] .
for @xmath183 , the angular power spectrum induced by an anti - biased distribution has more power than the unbiased case because the typical distance to subhalos is larger for the anti - biased distribution , shifting power to larger @xmath49 .
we should note that , of course , the electrons and positrons of ordinary astrophysical sources , such as pulsars and supernova remnants , could also produce significant ics signals . at @xmath183
the corresponding angular power spectrum would be difficult to distinguish from the one induced by dark matter because the large-@xmath49 power spectrum is highly suppressed and the suppression scale is mostly determined by the diffusion length of electrons and positrons which is independent of wether their origin is astrophysical or from dark matter .
in contrast , the power spectrum at @xmath196 should contain information on the sources of electrons and positrons due to the different spatial morphologies of dark matter and astrophysical sources on these scales which are not significantly influenced by diffusion .
in addition , the ics energy spectrum is in general significantly different due to the harder pair spectrum from dark matter annihilation which would naturally explain the @xmath41 excesses observed by pamela and fermi - lat .
although the nearby astrophysical sources such as pulsars can also reproduce these excesses , pulsars are mainly concentrated at the galactic plane and the primary @xmath41 injected spectrum have an exponential cutoff at high energies , both of which are quite different from annihilating dark matter scenarios .
therefore , a combination of these two effects may allow to discriminate between astrophysical sources and at least some dark matter scenarios .
for example , ref .
@xcite proposed an energy dependence of an anisotropy signature to distinguish millisecond pulsars from heavy dark matter candidates .
the disentanglement of dark matter signals from astrophysical backgrounds is , however , beyond the scope of the present paper .
in this paper we have investigated the angular power spectrum of the gamma - ray emission from inverse compton scattering off low energy target photons of electrons and positrons produced by dark matter annihilation .
we considered two extreme cases for the radial distribution of subhalos and simulated the full - sky gamma - ray maps at three energies through realizations of a large number of subhalos with masses down to @xmath150 .
the contributions to the angular power spectrum and to the total flux from the smooth host halo were also calculated .
1 . we point out a new feature in the angular power spectrum of photons produced by inverse compton scattering of pairs from annihilating dark matter that does not occur in the power spectra of gamma - rays produced directly in annihilation , on which earlier work often has focused : in contrast to the direct annihilation component for which the poisson noise dominates at large @xmath49 , the angular power spectrum from inverse compton scattering is exponentially suppressed at @xmath183 , because the diffusion of the high energy parent electrons and positrons produced by dark matter annihilation strongly smears out anisotropies of gamma - rays at small angular scales .
this suppression scale is determined by the diffusion length @xmath197 , evaluated at the energy @xmath30 of the emitting electrons and positrons .
the energy @xmath30 of the dominant emitters is determined by the low energy target photon backgrounds and , in the presence of several background components , is not a monotonic function of the gamma - ray energy . for the backgrounds used in this paper , cmb , ir and sl , we find that the lower energy gamma - rays are produced by the higher energy @xmath41 pairs ( and vice - versa ) so that when the energy dependence of the spatial diffusion coefficient is taken into account , @xmath165 , power spectrum suppression is found to occur at larger multipole @xmath49 for lower photon energy .
the contribution to the absolute amplitude of the angular power spectrum from the smooth host halo can be safely neglected which is due to its negligible mean intensity compared to that from the subhalos . for subhalos ,
the intensity from an unbiased subhalo distribution is @xmath198 times larger than from an antibiased distribution . at small @xmath49 the dimensionless angular power spectrum for the anti - biased distribution
is suppressed relative to the unbiased case , whereas at large @xmath183 the opposite is true .
an enormous number of subhalos with the masses down to @xmath199 can be well described by a smooth halo whose emissivity follows the spatial distribution of subhalos .
that is not surprising , because diffusion of the emitting particles can smooth the whole sky map if the number of subhalos is sufficiently large .
finally we remark that we have assumed homogeneous and isotropic magnetic fields and interstellar radiation field .
however , for the more realistic case where the radiation and magnetic fields depend on the position in the galaxy  @xcite , the angular power spectrum should be modified at least at large angular scales .
in particular , the small scale fluctuations of the magnetic field could affect the distribution of electrons and positrons  @xcite , influencing the angular power spectrum at all scales .
moreover , a space - dependence of the diffusion length could also lead to a directional dependence of the angular power spectrum .
obviously detecting the differential effect of the anisotropy signal across the galactic latitude would be even more difficult than detecting the signal integrated over the whole sky .
however , it could introduce subtle effects when part of the sky is masked out , e.g. to avoid the contribution from the galactic plane .
these complications should be investigated more thoroughly in the future by performing a real 3d simulation through numerically solving the transport equation for electrons and positrons .
while the formalism developed in this paper can be extended to the anisotropies of the radio sky , the small scale structure of the magnetic field is more important for synchrotron emission and should be treated in more detail .
in addition , our work could also be extended to inverse compton gamma - ray emission from extragalactic dark matter halos .
we acknowledge use of the healpix  @xcite software and analysis package for deriving the results in this paper .
this work was supported by the deutsche forschungsgemeinschaft through the collaborative research centre sfb 676 particles , strings and the early universe : the structure of matter and space - time , and by the state of hamburg , through the collaborative research program connecting particles with the cosmos within the framework of the landesexzellenzinitiative ( lexi ) .
we thank luca maccione and enrico borriello for reading the manuscript and for very useful comments .
99 g.  bertone , d.  hooper and j.  silk , `` particle dark matter : evidence , candidates and constraints , '' phys .  rept .
* 405 * , 279 ( 2005 ) [ arxiv : hep - ph/0404175 ] .
g.  jungman , m.  kamionkowski and k.  griest , `` supersymmetric dark matter , '' phys .  rep .
* 267 * , 195 ( 1996 ) [ arxiv : hep - ph/9506380 ] .
d.  hooper and s.  profumo , `` dark matter and collider phenomenology of universal extra dimensions , '' phys .  rep .   *
453 * , 29 ( 2007 ) [ arxiv : hep - ph/0701197 ] .
e.  komatsu _ et al .
_ [ wmap collaboration ] , `` five - year wilkinson microwave anisotropy probe observations : cosmological interpretation , '' arxiv:0803.0547 [ astro - ph ] .
x.  l.  chen and m.  kamionkowski , `` particle decays during the cosmic dark ages , '' phys .  rev .
d * 70 * , 043502 ( 2004 ) [ arxiv : astro - ph/0310473 ] .
l.  zhang , x.  chen , m.  kamionkowski , z.  g.  si and z.  zheng , `` contraints on radiative dark - matter decay from the cosmic microwave background , '' phys .
d * 76 * , 061301 ( 2007 ) [ arxiv:0704.2444 [ astro - ph ] ] .
l.  zhang , x.  l.  chen , y.  a.  lei and z.  g.  si , `` the impacts of dark matter particle annihilation on recombination and the anisotropies of the cosmic microwave background , '' phys .  rev .
d * 74 * , 103519 ( 2006 ) [ arxiv : astro - ph/0603425 ] ;          a.  a.  abdo _ et al .
_ [ the fermi - lat collaboration ] , `` measurement of the cosmic ray e+ plus e- spectrum from 20 gev to 1 tev with the fermi large area telescope , '' phys .
lett .   * 102 * , 181101 ( 2009 ) [ arxiv:0905.0025 [ astro-ph.he ] ] .
d.  grasso _ et al .
_ [ fermi - lat collaboration ] ,  on possible interpretations of the high energy electron - positron spectrum astropart .
* 32 * ( 2009 )
140 [ arxiv:0905.0636 [ astro-ph.he ] ] . g.  di bernardo , d.  gaggero , d.  grasso and f.  l.  collaboration , `` the high energy cosmic ray electron spectrum measured by fermi gamma - ray space telescope : some possible interpretations , '' arxiv:0912.3887 [ astro-ph.he ] .
m.  cirelli , m.  kadastik , m.  raidal and a.  strumia , `` model - independent implications of the e+ , e- , anti - proton cosmic ray spectra on properties of dark matter , '' nucl .  phys .
b * 813 * , 1 ( 2009 ) [ arxiv:0809.2409 [ hep - ph ] ] .
i.  cholis , l.  goodenough , d.  hooper , m.  simet and n.  weiner , `` high energy positrons from annihilating dark matter , '' arxiv:0809.1683 [ hep - ph ] .
v.  barger , w.  y.  keung , d.  marfatia and g.  shaughnessy , `` pamela and dark matter , '' arxiv:0809.0162 [ hep - ph ] .
r.  harnik and g.  d.  kribs , `` an effective theory of dirac dark matter , '' arxiv:0810.5557 [ hep - ph ] .
a.  e.  nelson and c.  spitzer , `` slightly non - minimal dark matter in pamela and atic , '' arxiv:0810.5167 [ hep - ph ] .
y.  bai , m.  carena and j.  lykken , `` the pamela excess from neutralino annihilation in the nmssm , '' phys .
d * 80 * , 055004 ( 2009 ) [ arxiv:0905.2964 [ hep - ph ] ] .
d.  malyshev , i.  cholis and j.  gelfand , `` pulsars versus dark matter interpretation of atic / pamela , '' arxiv:0903.1310 [ astro-ph.he ]
. s.  profumo , `` dissecting pamela ( and atic ) with occam s razor : existing , well - known pulsars naturally account for the anomalous cosmic - ray electron and positron data , '' arxiv:0812.4457 [ astro - ph ] .
t.  kobayashi , y.  komori , k.  yoshida and j.  nishimura , `` the most likely sources of high energy cosmic - ray electrons in supernova remnants , '' astrophys .
j.   * 601 * , 340 ( 2004 ) [ arxiv : astro - ph/0308470 ] .
n.  j.  shaviv , e.  nakar and t.  piran , `` natural explanation for the anomalous positron to electron ratio with supernova remnants as the sole cosmic ray source , '' phys .
lett .   * 103 * , 111302 ( 2009 ) [ arxiv:0902.0376 [ astro-ph.he ] ] .
m.  cirelli and p.  panci , `` inverse compton constraints on the dark matter e+e- excesses , '' nucl .
b * 821 * , 399 ( 2009 ) [ arxiv:0904.3830 [ astro-ph.co ] ] .
e.  a.  baltz and l.  wai , `` diffuse inverse compton and synchrotron emission from dark matter annihilations in galactic satellites , '' phys .
d * 70 * , 023512 ( 2004 ) [ arxiv : astro - ph/0403528 ] .
q.  yuan , p.  f.  yin , x.  j.  bi , x.  m.  zhang and s.  h.  zhu , `` gamma rays and neutrinos from dark matter annihilation in galaxy clusters , '' phys .
d * 82 * , 023506 ( 2010 ) [ arxiv:1002.0197 [ astro-ph.he ] ] .
c.  boehm , t.  delahaye and j.  silk , `` can the morphology of gamma - ray emission distinguish annihilating from decaying dark matter ?
, '' arxiv:1003.1225 [ astro-ph.ga ] . g.  hutsi , a.  hektor and m.  raidal , `` implications of the fermi - lat diffuse gamma - ray measurements on annihilating or decaying dark matter , '' jcap * 1007 * , 008 ( 2010 ) [ arxiv:1004.2036 [ astro-ph.he ] ] .
t.  lin , d.  p.  finkbeiner and g.  dobler , `` the electron injection spectrum determined by anomalous cosmic ray , gamma ray , and microwave signals , '' arxiv:1004.0989 [ astro-ph.co ] .
a.  loeb and m.  zaldarriaga , `` the small - scale power spectrum of cold dark matter , '' phys .
d * 71 * , 103520 ( 2005 ) [ arxiv : astro - ph/0504112 ] .
e.  bertschinger , `` the effects of cold dark matter decoupling and pair annihilation on cosmological perturbations , '' phys .
d * 74 * , 063509 ( 2006 ) [ arxiv : astro - ph/0607319 ] .
s.  profumo , k.  sigurdson and m.  kamionkowski , `` what mass are the smallest protohalos ?
, '' phys .
lett .   * 97 * , 031301 ( 2006 ) [ arxiv : astro - ph/0603373 ] .
j.  diemand , b.  moore and j.  stadel , `` earth - mass dark - matter haloes as the first structures in the early universe , '' nature * 433 * , 389 ( 2005 ) [ arxiv : astro - ph/0501589 ] .
l.  pieri , g.  bertone and e.  branchini , `` dark matter annihilation in substructures revised , '' mon .  not .
soc .   * 384 * , 1627 ( 2008 ) [ arxiv:0706.2101 [ astro - ph ] ] .
m.  kamionkowski , s.  m.  koushiappas and m.  kuhlen , `` galactic substructure and dark matter annihilation in the milky way halo , '' phys .
d * 81 * , 043532 ( 2010 ) [ arxiv:1001.3144 [ astro-ph.ga ] ] .
s.  ando and e.  komatsu , `` anisotropy of the cosmic gamma - ray background from dark matter annihilation , '' phys .
d * 73 * ( 2006 ) 023521 [ arxiv : astro - ph/0512217 ] .
s.  ando , e.  komatsu , t.  narumoto and t.  totani , `` dark matter annihilation or unresolved astrophysical sources ?
anisotropy probe of the origin of cosmic gamma - ray background , '' phys .
d * 75 * , 063519 ( 2007 ) [ arxiv : astro - ph/0612467 ] .
s.  ando , `` gamma - ray background anisotropy from galactic dark matter substructure , '' arxiv:0903.4685 [ astro-ph.co ] .
m.  fornasa , l.  pieri , g.  bertone and e.  branchini , `` anisotropy probe of galactic and extra - galactic dark matter annihilations , '' arxiv:0901.2921 [ astro - ph ] .
a.  cuoco , j.  brandbyge , s.  hannestad , t.  haugboelle and g.  miele , `` angular signatures of annihilating dark matter in the cosmic gamma - ray background , '' phys .
d * 77 * ( 2008 ) 123518 [ arxiv:0710.4136 [ astro - ph ] ] .
b.  s.  hensley , j.  m.  siegal - gaskins and v.  pavlidou , `` the detectability of dark matter annihilation with fermi using the anisotropy energy spectrum of the gamma - ray background , '' arxiv:0912.1854 [ astro-ph.co ] .
j.  zavala , v.  springel and m.  boylan - kolchin , `` extragalactic gamma - ray background radiation from dark matter annihilation , '' mon .  not .
* 405 * , 593 ( 2010 ) [ arxiv:0908.2428 [ astro-ph.co ] ] .
f.  miniati , `` numerical modelling of gamma radiation from galaxy clusters , '' mon .  not .
* 342 * , 1009 ( 2003 )              t.  delahaye , f.  donato , n.  fornengo , j.  lavalle , r.  lineros , p.  salati and r.  taillet , `` galactic secondary positron flux at the earth , '' astron .
astrophys .
* 501 * , 821 ( 2009 ) [ arxiv:0809.5268 [ astro - ph ] ] .        strong a.w . and moskalenko i.w .
, `` propagation of cosmic - ray nucleons in the galaxy , '' astrophys .  j.   * 509 * ( 1998 ) 212 [ astro - ph/9807150 ] .
t.  a.  porter , i.  v.  moskalenko and a.  w.  strong , `` inverse compton emission from galactic supernova remnants : effect of the interstellar radiation field , '' astrophys .
j.   * 648 * , l29 ( 2006 ) [ arxiv : astro - ph/0607344 ] .
t.  delahaye , r.  lineros , f.  donato , n.  fornengo and p.  salati , `` positrons from dark matter annihilation in the galactic halo : theoretical uncertainties , '' phys .
d * 77 * , 063527 ( 2008 ) [ arxiv:0712.2312 [ astro - ph ] ] .
d.  maurin , f.  donato , r.  taillet and p.  salati , `` cosmic rays below z=30 in a diffusion model : new constraints on propagation parameters , '' astrophys .
j.   * 555 * , 585 ( 2001 ) [ arxiv : astro - ph/0101231 ] .
rybicki g.b . ,
lightman a.p .
, radiative processes in astrophysics , 1979 , john wiley & sons .
k.  m.  gorski , e.  hivon , a.  j.  banday , b.  d.  wandelt , f.  k.  hansen , m.  reinecke and m.  bartelman , `` healpix  a framework for high resolution discretization , and fast analysis of data distributed on the sphere , '' astrophys .
j.   * 622 * , 759 ( 2005 ) [ arxiv : astro - ph/0409513 ] .
m.  perelstein and b.  shakya , `` comment on calculation of positron flux from galactic dark matter , '' arxiv:1002.4588 [ astro-ph.he ] .
j.  m.  siegal - gaskins , r.  reesman , v.  pavlidou , s.  profumo and t.  p.  walker , `` anisotropy constraints on millisecond pulsars in the diffuse gamma ray background , '' arxiv:1011.5501 [ astro-ph.he ] .
g.  dobler , d.  p.  finkbeiner , i.  cholis , t.  r.  slatyer and n.  weiner , `` the fermi haze : a gamma - ray counterpart to the microwave haze , '' astrophys .  j.   * 717 * , 825 ( 2010 ) [ arxiv:0910.4583 [ astro-ph.he ] ] .