Astrophysical Uncertainties in the Cosmic Ray Electron and Positron Spectrum From Annihilating Dark Matter
Abstract
In recent years, a number of experiments have been conducted with the goal of studying cosmic rays at GeV to TeV energies. This is a particularly interesting regime from the perspective of indirect dark matter detection. To draw reliable conclusions regarding dark matter from cosmic ray measurements, however, it is important to first understand the propagation of cosmic rays through the magnetic and radiation fields of the Milky Way. In this paper, we constrain the characteristics of the cosmic ray propagation model through comparison with observational inputs, including recent data from the CREAM experiment, and use these constraints to estimate the corresponding uncertainties in the spectrum of cosmic ray electrons and positrons from dark matter particles annihilating in the halo of the Milky Way.
I Introduction
In order for dark matter particles to be detected, they must either interact directly with particles of the standard model, or produce such particles through their annihilations. In this paper, we concern ourselves with the indirect detection of dark matter through the observation of its annihilation products. If the dark matter consists of relic particles with weakscale interactions and masses, they will produce a combination of gamma rays, neutrinos, and other standard model particles in their annihilations. Of particular observational interest are positrons history () and antiprotons antiprotons () which, in lieu of dark matter annihilations or other sources, are expected to be rare compared to their nonantimatter counterparts in the cosmic ray spectrum.
The most significant astrophysical source of cosmic ray positrons and antiprotons had generally been expected to be secondary production (creation through cosmic ray interactions occurring in the interstellar medium of the Milky Way). An excess of positrons or antiprotons in the cosmic ray spectrum relative to the small fraction expected as secondaries would imply the existence of primary sources of such particles. Although dark matter annihilations may be capable of producing such a signature history (); BaltzEdsjo (); HooperSilk (), astrophysical sources have also been proposed, including pulsars Aharonian et al. (1995). Very recently, it has also been suggested that secondary production of positrons within the regions of cosmic ray acceleration may also produce a significant fraction of the cosmic ray positron flux at high energies Blasi:2009hv (). In order to distinguish between such sources of cosmic ray antimatter, a detailed understanding of the processes involved in cosmic ray propagation is likely to be required.
The satellitebased experiment PAMELA (a Payload for Antimatter Matter Exploration and Lightnuclei Astrophysics) has reported a steadily increasing positron fraction (the ratio of positrons to positronspluselectrons) in the cosmic ray spectrum between approximately 10 GeV and 100 GeV Adriani et al. (2008) (data has not yet been published at higher energies). While consistent with previous indications from the HEAT heat () and AMS01 ams01 () experiments, this result is in stark contrast to the behavior expected if the positron spectrum were dominated by secondary particles produced during cosmic ray propagation secondaries (). Although the origin of these particles is not yet known, their observation suggests the existence of a primary source (or sources) of positrons serpico (). PAMELA’s measurement of the antiprotontoproton ratio, in contrast, is consistent with purely secondary production pamelaantiproton (). The antiproton result from PAMELA can be used to limit the nature of positron primary sources, including on the range of dark matter models potentially responsible antiprotonconst (). PAMELA is ultimately expected to measure the spectra of cosmic ray protons, antiprotons, electrons, and positrons up to energies of 700 GeV, 190 GeV, 2 TeV, and 270 GeV, respectively.
ATIC (Advanced Thin Ionization Calorimeter) is a balloonbased experiment designed to study the spectra of cosmic ray protons, light nuclei, and electrons. Recently, the ATIC collaboration published their electron (plus positron, as they do not distinguish between these species) spectrum between approximately 20 GeV and 2 TeV Chang et al. (2008). Remarkably, they find that the spectrum contains a bumplike feature over the steadily declining powerlaw between roughly 300 and 800 GeV, peaking at around 600 GeV. The energy loss and diffusion rates of electrons in this energy range lead us to conclude that these particles must originate within approximately 1 kpc of the Solar System.
In light of these recent excesses observed in the cosmic ray positron fraction and electron (plus positron) spectrum, we revisit this topic, concentrating on the astrophysics that goes into determining the cosmic ray electron and positron spectra resulting from annihilating dark matter. After an initial spectrum of electrons and positrons is created through dark matter annihilations, that spectrum evolves as it diffuses and propagates through the radiation fields and magnetic field of the Milky Way. The properties of the propagation model used to describe this evolution can be constrained from other cosmic ray observations, including the relative abundances of unstable and stable secondaries in the cosmic ray spectrum, which provide us with information pertaining to the characteristic timescales over which the particles have been propagating and the integrated density of matter through which it has passed. Together, such measurements enable us to construct a reasonably constrained parameterization of the propagation model, which can then be used to calculate, among other things, the propagated spectrum of cosmic ray electrons and positrons from dark matter annihilations.
The remainder of this paper is structured as follows. In Section II, we use a large array of cosmic ray nuclei data in conjunction with simulations conducted using the GALPROP program (v50p) Strong & Moskalenko (1998) to constrain the characteristics of the model describing cosmic ray propagation in our galaxy. In Section III, we apply this constrained model to the electrons and positrons produced in dark matter annihilations and compare our results to the positron fraction measured by PAMELA. We summarize our results in Section IV.
Ii Constraining the Propagation Model
To model the diffusion, nuclear interactions, and energy loss processes of galactic cosmic rays, we use the publicly available GALPROP code Strong & Moskalenko (1998). This code begins by injecting cosmic rays with a supernovalike isotope distribution and a parameterized spatial distribution chosen to reproduce the EGRET data. The individual isotopes, starting with the largest atomic number, are then propagated through the galaxy, with the energy distribution and composition altered through a combination of spacial diffusion, energy losses, diffusion in momentum space (diffusive reacceleration), electron Kcapture, convection, spallation, and radioactive decay.
These processes are collectively described by the cosmic ray propagation equation. For particles of momentum with particle density per unit momentum , this equation is given by Strong & Moskalenko (1998):
(1)  
where is the diffusion constant, which we parameterize by
(2) 
where is the particle’s speed and is its rigidity (momentum per unit charge). The diffusion constant describing reacceleration, , is related to by pp (); berebook ()
(3) 
where the Alfvén speed. In the propagation equation (Eq. 1), is the convection velocity, is the fragmentation time, is the radioactive decay time, and is the source term. The source term includes not only the injection spectrum, but also the products of the decay and spallation of heavier species of nuclei. The equation is solved assuming a cylindrical geometry over the diffusive region (a volume with a halfthickness of and a radius of 20 kpc). Outside of this volume, the particles are not confined by the Galactic Magnetic Field and freely escape. Assuming that the system is near steady state equilibrium, we set the left side of Eq. 1 to zero in solving for .
Note that we refer to the halfthickness of the diffusive region as rather than simply as . Although GALPROP treats the diffusion constant to be the same in all directions (and locations throughout the diffusion zone), the actual process of diffusion in the galaxy is the result of magnetic fields which are not spherically symmetric, but are instead structured and thought to be disklike in form. The quantity thus does not denote simply the physical thickness of the diffusion zone, but may also include information regarding the differing efficiencies of vertical and horizontal diffusion.
In the course of this study, we completed 370 simulations with GALPROP. In these runs, we varied four parameters: the normalization of the diffusion coefficient (), the slope of the dependence of the diffusion coefficient with rigidity (), the effective halfthickness of the diffusion region (), and the convection velocity (). We considered values of between and (at a reference rigidity of 4 GeV), between 0.34 to 0.52 and between 1 and 17 kpc. was varied between 0 and 15 km/s/kpc, where we have assumed that the convection velocity scales proportionally to the distance from the Galactic Plane. We have adopted a twodimensional (cylindrical) symmetry which allows physical quantities to vary with and , but not with . Unless stated, all other parameters are unchanged from those found in the GALPROP definitions file (galdef_50p_599278) galpropweb (). In particular, note that we have not changed the Alfvén speed from its default value of 36 km/s. Some previous studies (for example, see Ref. Bergström et al. (2006)) have allowed to vary and have found a fairly large range of values to be consistent with cosmic ray data, so long as compensating changes are made in other propagation parameters accordingly. With the recent introduction of new high energy data from the CREAM experiment Ahn et al. (2008), however, this parameter is more tightly constrained.
Once the GALPROP simulations were completed, the spectra of the various elements and isotopes at the location of the Solar System ( kpc, ) were extracted and combined to yield predictions for both stable and unstable secondarytoprimary ratios. Stable secondarytoprimary ratios, such as B/C and subFe/Fe, are valuable measures of the average amount of matter traversed by cosmic rays as a function of energy. Although B/C is especially useful, subFe/Fe can provide complementary information, as it is most sensitive to a slightly different range of energies. Unstable secondarytoprimary ratios, in contrast, serve as a measurement of the time cosmic rays have been propagating. Beryllium10 is particularly useful in this regard, being the longest lived and best measured unstable secondary. The measurement of Be/Be serves as a clock, since the ratio of the radioactive isotope to the stable one is directly related to the amount of time elapsed since the creation of the particles. For an excellent review of this and related topics, see Ref. Strong, Moskalenko & Ptuskin (2007).
Authors  Ratios  Number 

Ahn et al. (2008) (CREAM)  B/C  (1) 
Panov et al. (2007) (ATIC)  B/C  (2) 
de Nolfo et al. (2001) (ISOMAX)  Be/Be  (3) 
Engelmann et al. (1990) (HEAO3)  B/C, subFe/Fe  (4) 
Swordy et al. (1990) (Spacelab 2)  B/C  (5) 
Dwyer and Meyer (1987)  B/C, subFe/Fe  (6) 
Webber et al. (1985)  B/C  (7) 
Chappell and Webber (1981)  B/C  (8) 
Simon et al. (1980)  B/C  (9) 
Orth et al. (1978)  B/C  (10) 
Caldwell and Meyer (1977) (IMP8)  B/C  (11) 
Juliusson (1974)  B/C  (12) 
Authors  Ratios  Number 

Webber et al. (2003) (Voyager)  B/C, subFe/Fe  (13) 
Webber et al. (2002) (Voyager)  Be/Be  (14) 
Hams et al. (2001) (Ulysses)  Be/Be  (15) 
Davis et al. (2000) (CRIS)  B/C, subFe/Fe  (16) 
Yanasak et al. (1999) (CRIS)  Be/Be  (17) 
Connell (1998) (Ulysses)  Be/Be  (18) 
DuVernois et al. (1996) (Ulysses)  B/C, subFe/Fe  (19) 
Leske (1993) (ISEE3)  subFe/Fe  (20) 
Esposito et al. (1992) (ALICE)  subFe/Fe  (21) 
Krombel & Wiedenbeck (1988) (ISEE 3)  B/C  (22) 
GarciaMunoz et al. (1987) (IMP8)  B/C  (23) 
GarciaMunoz et al. (1981)  Be/Be  (24) 
Young et al. (1981)  subFe/Fe  (25) 
Wiedenbeck & Greiner (1980)  Be/Be  (26) 
Webber et al. (1979)  subFe/Fe  (27) 
Buffington et al. (1978)  Be/Be  (28) 
Lezniak & Webber (1978)  B/C  (29) 
Hagen et al. (1976)  B/C, 10Be/9Be  (30) 
Maehl et al. (1977)  B/C  (31) 
Benegas et al. (1975)  subFe/Fe  (32) 
Lund et al. (1975)  B/C  (33) 
For each set of diffusion parameters, we have compared the predicted secondarytoprimary ratios to the body of current observational data, as listed in Table 1. We have limited the data included in our analysis to energies above 5 GeV for the B/C and subFe/Fe ratios, and above 1 GeV for the Be/Be data (no data above 5 GeV is currently available for this ratio). At lower energies, the effects of solar modulation become increasingly important, making it difficult to reliably compare the data to the predictions of a given propagation model. We also list in Table 2 a collection of cosmic ray data taken only at energies below those we have considered in our analysis. In the case of the Be/Be measurement at 12 GeV, in an effort to lessen the impact of solar modulation, we have included the correction as described in Ref. de Nolfo et al. (2001).
In Fig. 1, we plot the acceptable ranges of the propagation model parameters, based on our calculation. Neglecting the effects of convection for the time being, we find the overall best fit for the following set of propagation parameters: cm s, kpc, and , which yields a per degreeoffreedom of 1.37. Although the acceptable parameter regions are well constrained in and , the allowed values for extend beyond the range we considered (117 kpc). For physical reasons, however, we do not consider values outside of this range. We find that the data prefer values of which lie between those predicted for Kolmogorovtype () kolmogorov () and Kraichnantype () kraichnan () turbulence.
Model  A1  A2  A3  

(cm/s)  
1  (kpc)  2.04  4.52  17.0 
0.39  0.42  0.41  
Model  A4  A5  A6  
(cm/s)  
2  (kpc)  1.0  3.0  17.0 
0.40  0.44  0.39 
In Table 3, we list the parameters for a selection of extrema and central value propagation models. The B/C, subFe/Fe and Be/Be ratios predicted in these models, along with those predicted in the best fit model, are compared to the cosmic ray nuclei data in Figs. 2 and 3. These figures confirm that this range of propagation models provides a reasonably good fit to the current set of cosmic ray nuclei data over the energy range considered.
To illustrate the effect of changing the Alfvén speed, we performed additional GALPROP runs with km/s, , and a range of values for and . The best per degreeoffreedom found was 1.96, which is considerably larger than the 1.39 found with our default value ( km/s). From this test, we conclude that large variations in the Alfvén speed relative to our default choice are disfavored by the current data set.
Thus far, our analysis has neglected the effects of convection (preferential motion of cosmic rays away from the Galactic Plane). To explore how convection effects our results, we performed several further runs of GALPROP. From these runs, we learned that the fractional change in the element ratios due to convection was approximately of the same shape across all values of the convection velocity, , at least up to 15 km/s/kpc, with an overall normalization varying linearly with the velocity. The ratios with convection were thus estimated by calculating the fractional change with =15 km/s/kpc, and then multiplying the interpolated result with the matching fractional change scaled by km/s/kpc). 2 regions in the plane are given for =0, 5, and 10 km/s/kpc in Fig. 4. This illustrates how convection could potentially alter the allowed range of the propagation parameters.
Before turning our attention to the propagation of electrons and positrons from dark matter annihilations, we would like to make some general remarks regarding the methods used in this section. Beginning with the cosmic ray propagation model as described by Eq. 1, we have found a range of parameters which yield good agreement with the current cosmic ray data at energies high enough to be only modestly impacted by solar modulation. Although this propagation model contains a wide variety of physical effects, including spatial diffusion, energy losses, diffusive reacceleration, electron Kcapture, convection, spallation, and radioactive decay, it does have its limitations. For example, we have implicitly assumed that the diffusion constant, , does not vary with location (within the boundary conditions). Although such simplifying approximations are currently necessary to make the problem of constraining the propagation model tractable, we can be certain that at some point in the future this model will break down (fail to adequately describe the observations) and require a more sophisticated treatment. As our best fit parameter sets provide reasonable fits to the current data, however, it is appears that the approach described here is a reasonably accurate, or at an least adequate, approximation of the behavior of cosmic ray propagation in the Milky Way.
Iii High Energy Electrons And Positrons From Dark Matter Annihilations In the Milky Way
In this section, we apply the constraints on the cosmic ray propagation model obtained in the previous section to the problem of cosmic ray electron and positron propagation. We begin by adapting Eq. 1 to the case of high energy electrons and positrons. In particular, we remove the terms describing fragmentation and radioactive decay. Furthermore, we focus our analysis uniquely on cosmic ray electrons and positrons with energies above 10 GeV, at which the effects of reacceleration and convection are expected to be negligible. In the steady state limit, the propagation equation for high energy electrons and positrons reduces to
(4) 
where is the energy loss rate of electrons due to synchrotron and inverse Compton processes. In the relativistic limit, this rate is related to the radiation field and magnetic field energy densities by
(5)  
where is the Thompson cross section and is the representative energy loss time. There is also a contribution to the energy loss rate due to Bremsstrahlung, but this is significant only at energies lower than those considered here.
The radiation fields contributing to the inverse Compton loss rate include starlight, emission from dust, and the cosmic microwave background. In Fig. 5, we show the density of the interstellar radiation field (ISRF) as estimated in Ref. Porter et al. (2005) at a distance of kpc from the Galactic Center as a function of the distance away from the Galactic Plane, along with the uniform cosmic microwave background density. In addition, the magnetic field energy density which leads to synchrotron losses is related to the RMS field strength by . For G (5 G) in the local Milky Way, we arrive at a reasonable estimate of eV/cm (0.6 eV/cm).
As we did for the diffusion constant, we perform our analysis under the approximation that the radiation and magnetic field densities do not have a strong spatial dependence within the relevant volume of the diffusion zone. Considering the distribution shown in Fig. 5, the radiation field density varies by only 10% (30%) out to distances of 2 kpc (4 kpc) away from the Galactic Plane. As high energy electrons and positrons originating from greater distances are expected to lose most of their energy before reaching the Solar System, the spatial variation of the energy loss rate is not expected to significantly impact the observed spectrum, except perhaps at the lowest energies considered here ( 1020 GeV).
In this section, we consider the range of propagation parameters shown in Fig. 1 (the 1 and 2 contour regions), which collectively represent the range of propagation models found to be consistent with the cosmic ray data. Even without considering any primary sources of cosmic ray positrons, the variation over this range of propagation parameters leads to a range of predicted cosmic ray electron (primary plus secondary) and positron (secondary) spectra. Beginning with an injected (primary) spectrum of electrons described by , we have used GALPROP to calculate the resulting electron and positron spectra after propagation for each of the acceptable parameter sets. For this injected electron spectrum, each of the parameter sets leads to an electron spectrum at the Solar System with a slope over the range of 5100 GeV of , where , which is in reasonable agreement with the measured cosmic ray electron slope electronspec (). The predicted positron fraction over this energy range, as well as the electron spectrum at higher energies, are not in agreement with the observations of the PAMELA or ATIC experiments, however. The predicted positron fraction (top) and electronpluspositron spectrum (bottom) are shown in Fig. 6, compared to the data from PAMELA and ATIC. It is clear that none of these propagation models lead to an acceptable fit to the data. Another source of cosmic ray positrons/electrons is required to accommodate these observations.
Turning our attention now to the source term in the propagation equation, we assume that dark matter is distributed throughout the galaxy following an NavarroFrenkWhite (NFW) halo density profile Navarro et al. (1996),
(6) 
with a local dark matter density, , of 0.3 GeV/cm. The uncertainty introduced by limiting ourselves to such a profile is modest for the energy range we are considering (see Ref. HooperSilk () for a comparison of different halo profiles). For example, the electron/positron spectrum resulting from a highly cusped Moore et al. profile moore () or an isothermal sphere profile are virtually indistinguishable from that predicted for an NFW profile above 40 GeV, and vary by less than a factor of 2 at 10 GeV. The lack of a strong dependence on the halo profile results from the fact that high energy electrons lose energy relatively quickly, leading the observed spectrum to be dominated by particles originating from the surrounding few kiloparsecs.
In calculating the spectrum of electrons and positrons injected through dark matter annihilations, we consider in our analysis three particle physics scenarios:

A 200 GeV WIMP which annihilates to (which then generate electrons and positrons in their decays). A winolike neutralino, for example, could annihilate with a large cross section through this channel, although a nonthermal mechanism would be required to populate them in the early universe winos ().

A 600 GeV WIMP which annihilates simply to . Although this case serves primarily as a phenomenological benchmark, models have been constructed in which dark matter annihilates to light states which decay uniquely to electrons and positrons arkani () (for a scenario which produces muons in a similar way, see Ref. Nomura:2008ru ()).

A 600 GeV KaluzaKlein dark matter particle in a model with a single flat universal extra dimension kkdm (). In this case, the dark matter annihilates to , , and 20% of the time each, leading to a very hard spectrum of electrons and positrons Hooper:2004xn (). Most the remaining annihilations produce uptype quarks.
These three scenarios do not, and are not intended to, cover the full range of phenomenological possibilities for annihilating dark matter. Fits of the PAMELA and ATIC data to a broad range of dark matter masses and annihilation modes have been presented elsewhere modelindep (). Here, we focus our study on how the uncertainties in the propagation model can impact the electron and positron spectrum for the three representative cases described above. In each case, we have used the program PYTHIA Sjostrand et al. (2001), as implemented within DarkSUSY Gondolo et al (2004) to calculate the injected spectra of electrons and positrons.
Throughout this study, we normalize the dark matter annihilation rate to a default annihilation cross section of cm/s, multiplied by a “boost factor”. Such a boost factor could originate, for example, as a result of small scale inhomogeneities in the spatial distribution of dark matter. Alternatively, one could consider nonthermal relics with considerably larger cross sections than our default value, or a thermal relic which annihilates with a very large cross section at low velocities as a result of nonperturbative processes sommerfeld (). Variations in the annihilation cross section, boost factor, or average dark matter density lead to changes in the overall normalization of the electron/positron spectrum, but not in the spectral shape.
In Fig. 7, we plot the spectrum of electrons and positrons, and the positron fraction, for the case of a 200 GeV particle annihilating to . In the upper and lower frames, we have used two different values for the representative energy loss time, and seconds, corresponding to radiation and magnetic field densities of eVcm and eVcm, respectively. As found in previous studies modelindep (); winospamela (), the positron fraction predicted climbs less rapidly with energy (if at all) than is observed by PAMELA. The set of propagation parameters which provides the best fit to the cosmic ray data discussed in the previous section (shown as solid lines in each frame), yields a very poor fit to the PAMELA data for this dark matter model (degreeoffreedom =3.63 and 3.71 for and seconds, respectively). While acceptable variations of the propagation parameters can improve these fits marginally, nothing approaching a good fit can be found.
Even in the unrealistic case that the energy loss rate is reduced to the absolute minimum possible value, corresponding to inverse Compton scattering with the cosmic microwave background alone ( seconds, resulting from a CMB energy density of 0.26 eV/cm), the positron spectrum from a 200 GeV dark matter particle annihilating to is too low at high energies (or, if a larger boost factor is used, too high at low energies) to provide a good fit to the PAMELA data winospamela (). Also note that even if this sort of dark matter candidate could potentially generate a signal not very different from that seen by PAMELA, it will not lead to the distinctive feature seen by ATIC.
The inability of a dark matter candidate which annihilates to to accommodate the positron fraction measured by PAMELA provides us with a motivation to consider dark matter candidates which annihilate directly to electronpositron pairs or to other charged leptons (for examples of scenarios which predict dark matter annihilations to charged leptons, see Ref. leptons ()). In Fig. 8, we show the result for the extreme case of a dark matter particle which annihilates uniquely to electronpositron pairs. For this dark matter model, we find that the rapid rise of the positron fraction measured by PAMELA can easily be accommodated. Excellent fits to the PAMELA positron fraction are provided by such a dark matter particle, yielding per degreeoffreedom as low as 0.28 and 0.25 for and seconds, respectively. The normalization does require somewhat large boost factors, however, typically within the range of 150 to 350 (for the choice of a 600 GeV mass, selected to accommodate the observed ATIC feature).
In Fig. 9, we show the results for the case of a 600 GeV KaluzaKlein dark matter particle, which annihilations to , , and 20% of the time each. Although the positron fraction grows somewhat less rapidly than in the previous case, the result is still in good agreeement with the measurement of PAMELA. In particular, we find per degreeoffreedom as low as 0.34 and 0.51 for and seconds, respectively. Once again, we find that large boost factors are required. Note that in this case, in order to more easily facilitate direct comparison, we continue to use cm/s, although KaluzaKlein dark matter is actually predicted to have a somewhat higher annihilation cross section.
For the cases of dark matter annihilating to , we expect an edgelike feature to result in the electronpositron spectrum at . The precise spectrum of this feature, however, depends on the detailed distribution of dark matter in the local neighborhood of the Milky Way. This is especially true in the case that is a few hundred GeV or larger. At 600 GeV, for example, an electron will lose most of its energy before traveling much farther than 1 kpc, leading the high energy spectrum to depend critically on the local dark matter distribution, as well as on the local magnetic structure and other features of the propagation model pohl (). For this reason, it is not possible to reliably predict the electron spectrum near with much precision.
Iv Summary and Conclusions
In this article, we have discussed the propagation of high energy cosmic ray electrons and positrons from dark matter annihilations taking place in the Galactic Halo. We have focused on the astrophysics of this problem, including the properties of the diffusionenergy loss model. To study this process, we have incorporated constraints from measurements of the abundances of stable and unstable secondary nuclei in the cosmic ray spectrum. Stable secondarytoprimary ratios, such as B/C and subFe/Fe, can be used to provide a measurement of the average amount of matter traversed by cosmic rays as a function of energy. Unstable secondarytoprimary ratios, such as Be/Be, can be used to infer the length of time that cosmic rays have been propagating. Taken together, such information can be used to strongly constrain a model of cosmic ray propagation, which can then be used to make predictions for the spectrum of cosmic ray electrons and positrons from a given model of particle dark matter.
We have found that a relatively simple singlezone (cylindrical) diffusion model with a powerlaw diffusion coefficient and freeescape boundary conditions (and including the effects of diffusive reacceleration, electron Kcapture, spallation, and radioactive decay) can describe all of the current cosmic ray data above 1 GeV. As the precision and quality of this data has improved (in particular, with the latest data from CREAM Ahn et al. (2008)), this simple model has held up quite well. At some level, this provides us with confidence that this model constitutes a reasonable description of the processes contributing to cosmic ray propagation in the Milky Way. Although inhomogeneities within the diffusion zone and other considerations will certainly lead to behavior departing from the predictions of this model, the current data do not require any such departures.
Using the propagation parameter sets which we have found to be consistent with the cosmic ray data, we studied the propagation of high energy ( GeV) electrons and positrons in the Galactic Halo. Although considerable variability in these results were found over the acceptable range of propagation parameters, these variations were not sufficient to dramatically change the conclusions reached regarding specific particle dark matter models and the PAMELA and ATIC data. In particular, dark matter annihilating to throughout the Milky Way halo does not provide a good fit to PAMELA or ATIC data for any acceptable propagation parameter set. In contrast, dark matter which annihilates to electronpositron pairs, or to other charged leptons, naturally provides the rapid climb in the positron fraction as observed by PAMELA (and can also lead to the feature in the electron spectrum observed by ATIC if the mass is chosen to be 600800 GeV). We thus conclude that, while variations in the cosmic ray propagation model can lead to substantial variations in the high energy cosmic ray electron and positron spectra from dark matter, we do not expect such variations to lead to qualitative differences to our conclusions regarding the nature of the particle dark matter model required to fit the results of PAMELA or ATIC.
This work has been supported by the US Department of Energy, including grant DEFG0295ER40896, and by NASA grant NAG510842.
References
 (1) M. S. Turner and F. Wilczek, Phys. Rev. D 42, 1001 (1990); M. Kamionkowski and M. S. Turner, Phys. Rev. D 43, 1774 (1991); E. A. Baltz, J. Edsjo, K. Freese and P. Gondolo, Phys. Rev. D 65, 063511 (2002) [arXiv:astroph/0109318]; G. L. Kane, L. T. Wang and J. D. Wells, Phys. Rev. D 65, 057701 (2002) [arXiv:hepph/0108138].
 (2) J. Silk and M. Srednicki, Phys. Rev. Lett. 53, 624 (1984); F. W. Stecker, S. Rudaz and T. F. Walsh, Phys. Rev. Lett. 55, 2622 (1985); S. Rudaz and F. W. Stecker, Astrophys. J. 325, 16 (1988); L. Bergstrom, J. Edsjo and P. Ullio, arXiv:astroph/9906034. A. Bottino, F. Donato, N. Fornengo and P. Salati, Phys. Rev. D 58, 123503 (1998) [arXiv:astroph/9804137]; F. Donato, N. Fornengo, D. Maurin and P. Salati, Phys. Rev. D 69, 063501 (2004) [arXiv:astroph/0306207];
 (3) E. A. Baltz and J. Edsjo, Phys. Rev. D 59, 023511 (1999) [arXiv:astroph/9808243]; S. Profumo and P. Ullio, JCAP 0407, 006 (2004) [arXiv:hepph/0406018].
 (4) D. Hooper and J. Silk, Phys. Rev. D 71, 083503 (2005) [arXiv:hepph/0409104].
 Aharonian et al. (1995) F. A. Aharonian, A. M. Atoyan, and H. J. Volk, Astron. Astrophys. 294, L41 (1995); I. Buesching, O. C. de Jager, M. S. Potgieter, and C. Venter, arXiv:0804.0220 [astroph]; L. Zhang, and K. S. Cheng, Astron. Astrophys. 368, 1063 (2001); D. Hooper, P. Blasi and P. D. Serpico, arXiv:0810.1527 [astroph]; H. Yuksel, M. D. Kistler and T. Stanev, arXiv:0810.2784 [astroph].
 (6) P. Blasi, arXiv:0903.2794 [astroph.HE]; P. Blasi and P. D. Serpico, arXiv:0904.0871 [astroph.HE].
 Adriani et al. (2008) o. Adriani, et al. [PAMELA Collaboration], arXiv:0810:4994.
 (8) S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482, L191 (1997) [arXiv:astroph/9703192].
 (9) M. Aguilar et al. [AMS01 Collaboration], Phys. Lett. B 646, 145 (2007) [arXiv:astroph/0703154].
 (10) I. V. Moskalenko and A. W. Strong, Astrophys. J. 493, 694 (1998) [arXiv:astroph/9710124].
 (11) P. D. Serpico, arXiv:0810.4846 [hepph].
 (12) O. Adriani et al., arXiv:0810.4994 [astroph].
 (13) F. Donato, D. Maurin, P. Brun, T. Delahaye and P. Salati, arXiv:0810.5292 [astroph]; M. Cirelli, M. Kadastik, M. Raidal and A. Strumia, arXiv:0809.2409 [hepph].
 Chang et al. (2008) J. Chang, et al. [ATIC Collaboration], Nature, 456, 362 (2008).
 Strong & Moskalenko (1998) A. W. Strong and I. V. Moskalenko, ApJ, 509, 212 (1998).
 (16) E. S. Seo and V. S. Ptuskin, ApJ, 431, 705 (1994);
 (17) V. S. Berezinskii, et al., Astrophysics of cosmic rays, edited by V. L. Ginzburg, Amsterdam: NorthHolland (1990).
 (18) http://galprop.stanford.edu/web_galprop/galprop_home.html
 Bergström et al. (2006) L. Bergström, J. Edsjö, M. Gustafsson, and P. Salati, JCAP, 5, 6 (2006).
 Ahn et al. (2008) H. S. Ahn, et al., arXiv:0808.1718.
 Strong, Moskalenko & Ptuskin (2007) Strong, A.W., Moskalenko, I.V., & Ptuskin, V.S. 2007, ARNPS, 57, 285
 Panov et al. (2007) A. D. Panov, et al., arXiv:0707.4415.
 de Nolfo et al. (2001) G. A. de Nolfo, et al., Proc. 27th Int. Cosmic Ray Conf. (Hamburg), 5, 1667 (2001).
 Engelmann et al. (1990) J. J. Engelmann, et al., A&A, 233, 96 (1990).
 Swordy et al. (1990) S. P. Swordy, D. Muller, P. Meyer, J. L’Heureux and J. M. Grunsfeld, ApJ, 349, 625 (1990).
 Dwyer and Meyer (1987) R. Dwyer and P. Meyer, ApJ, 322, 981 (1987).
 Webber et al. (1985) W. R. Webber, J. C. Kish, and D. A. Schrier, Proc. 19th Int. Cosmic Ray Conf. (La Jolla), 2, 16 (1985).
 Chappell and Webber (1981) J. H. Chappell and W. R. Webber, Proc. 17th Int. Cosmic Ray Conf. (Paris), 2, 59 (1981).
 Simon et al. (1980) M. Simon, H. Spiegelhauer, W. K. H. Schmidt, F. Siohan, J. F. Ormes, V. K. Balasubrahmanyan, and J. F. Arens, ApJ, 239, 712 (1980).
 Orth et al. (1978) C. D. Orth, A. Buffington, G. F. Smoot, and T. S. Mast, ApJ, 226, 1147 (1978).
 Caldwell and Meyer (1977) J. H. Caldwell and P. Meyer, Proc. 15th Int. Cosmic Ray Conf. (Plovdiv), 1, 243 (1977).
 Juliusson (1974) E. Juliusson, ApJ, 191, 331 (1974).
 (33) P. Goldreich and S. Sridhar, Ap.J., 438, 763 (1995).
 (34) H. Yan and A. Lazarian, Ap.J., 614, 757 (2004).
 Webber et al. (2003) W. R. Webber, F. B. McDonald, and A. Lukasiak, ApJ, 599, 582 (2003).
 Webber et al. (2002) W. R. Webber, A. Lukasiak, and F. B. McDonald, ApJ. 568, 210 (2002).
 Hams et al. (2001) T. Hams, et al. ICRC, 1655 (2001).
 Davis et al. (2000) A. J. Davis, et al. in AIP Conf. Proc. 528, ed. R.A. Mewaldt et al. (New York: AIP), 421 (2001).
 Yanasak et al. (1999) N. E. Yanasak, et al., ICRC, 3 9 (1999).
 Connell (1998) J. J. Connell, ApJ, 501, L59 (1998).
 DuVernois et al. (1996) M. A. DuVernois, J. A. Simpson, and M. R. Thayer, A&A, 316, 555 (1996).
 Leske (1993) R. A. Leske, ApJ, 405, 567 (1993).
 Esposito et al. (1992) J. A. Esposito, et al., Astroparticle Physics, 1, 33 (1992).
 Krombel & Wiedenbeck (1988) K. E. Krombel and M. E. Wiedenbeck, ApJ, 328, 940 (1988).
 GarciaMunoz et al. (1987) M. GarciaMunoz, J. A. Simpson, T. G. Guzik, J. P. Wefel, and S. H. Margolis, ApJS, 64, 269 (1987).
 GarciaMunoz et al. (1981) M. GarciaMunoz, J. A. Simpson, and J. P. Wefel, ICRC, 2, 72 (1981).
 Young et al. (1981) J. S. Young, P. S. Freier, C. J. Waddington, N. Brewster, and R. K, Fickle, ApJ, 246, 1014 (1981).
 Wiedenbeck & Greiner (1980) M. E. Wiedenbeck and D. E. Greiner, ApJ, 239L, 139 (1980).
 Webber et al. (1979) W. R. Webber, J. C. Kish, and G. A. Simpson, Proc. 16th Int. Cosmic Ray Conf. (Kyoto), 1, 430 (1979).
 Buffington et al. (1978) A. Buffington, C. D. Orth, and T. S. Mast, ApJ, 226, 355 (1978).
 Lezniak & Webber (1978) J. A. Lezniak and W. R. Webber, ApJ, 223, 676 (1978).
 Hagen et al. (1976) F. A. Hagen, A. J. Fisher, and J. F. Ormes, ApJ, 212, 262 (1977).
 Maehl et al. (1977) R. C. Maehl, J. F. Ormes, A. J. Fisher, and F. A. Hagen, Ap&SS, 47, 163 (1977).
 Benegas et al. (1975) J. C. Benegas, M. H. Israel, J. Klarmann, and R. C. Maehl, Proc. 14th Int. Cosmic Ray Conf. (Munich), 1, 251 (1975).
 Lund et al. (1975) N. Lund, I.L. Rasmussen, B. Peters, and N. J. Westergaard, Proc. 14th Int. Cosmic Ray Conf. (Munich), 1, 257 (1975).
 Porter et al. (2005) T. A. Porter, et al., Int. Cosmic Ray Conf., Pune, 4, 77 (2005).
 (57) J. Alcaraz et al. [AMS Collaboration], Phys. Lett. B 484, 10 (2000) [Erratumibid. B 495, 440 (2000)]; ATIC Collaboration, ICRC 2005, http://icrc2005.tifr.res.in/htm/VolWeb/Vol13/13001chnchangJabs1og11oral.pdf; S. Torii et al., Astrophys. J. 559, 973 (2001); K. Yoshida et al., Adv. Space Res. 42, 1670 (2008); M. Boezio, et al., Astrophys. J. 532, 653 (2000). M. A. DuVernois et al., Astrophys. J. 559, 296 (2001); C. Grimani et al., Astron. Astrophys. 392, 287 (2002).
 Navarro et al. (1996) J. F. Navarro, C. S. Frenk and S. D. M. White, ApJ, 462, 563 (1996).
 (59) B. Moore, T. R. Quinn, F. Governato, J. Stadel and G. Lake, Mon. Not. Roy. Astron. Soc. 310, 1147 (1999) [arXiv:astroph/9903164].
 (60) P. Grajek, G. Kane, D. J. Phalen, A. Pierce and S. Watson, arXiv:0807.1508 [hepph]; T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) [arXiv:hepph/9906527].
 (61) N. ArkaniHamed, D. P. Finkbeiner, T. Slatyer and N. Weiner, Phys. Rev. D 79, 015014 (2009) [arXiv:0810.0713 [hepph]]; I. Cholis, D. P. Finkbeiner, L. Goodenough and N. Weiner, arXiv:0810.5344 [astroph].
 (62) Y. Nomura and J. Thaler, arXiv:0810.5397 [hepph].
 (63) G. Servant and T. M. P. Tait, Nucl. Phys. B 650, 391 (2003) [arXiv:hepph/0206071]; H. C. Cheng, J. L. Feng and K. T. Matchev, Phys. Rev. Lett. 89, 211301 (2002) [arXiv:hepph/0207125]; D. Hooper and S. Profumo, Phys. Rept. 453, 29 (2007) [arXiv:hepph/0701197].
 (64) D. Hooper and G. D. Kribs, Phys. Rev. D 70, 115004 (2004) [arXiv:hepph/0406026]; D. Hooper and K. Zurek, arXiv:0902.0593 [hepph].
 (65) I. Cholis, L. Goodenough, D. Hooper, M. Simet and N. Weiner, arXiv:0809.1683 [hepph]; V. Barger, W. Y. Keung, D. Marfatia and G. Shaughnessy, Phys. Lett. B 672, 141 (2009) [arXiv:0809.0162 [hepph]]; M. Cirelli, M. Kadastik, M. Raidal and A. Strumia, arXiv:0809.2409 [hepph].
 Sjostrand et al. (2001) T. Sjostrand, et al., Comput. Phys. Commun., 135, 238 (2001).
 Gondolo et al (2004) P. Gondolo, et al., JCAP, 0407, 008 (2004) [arXiv:astroph/0012234].
 (68) M. Cirelli and A. Strumia, arXiv:0808.3867 [astroph]; I. Cholis, G. Dobler, D. P. Finkbeiner, L. Goodenough and N. Weiner, arXiv:0811.3641 [astroph]; N. ArkaniHamed, D. P. Finkbeiner, T. Slatyer and N. Weiner, arXiv:0810.0713 [hepph].
 (69) P. Grajek, G. Kane, D. Phalen, A. Pierce and S. Watson, arXiv:0812.4555 [hepph];
 (70) P. J. Fox and E. Poppitz, arXiv:0811.0399 [hepph]; R. Harnik and G. D. Kribs, arXiv:0810.5557 [hepph]; K. M. Zurek, arXiv:0811.4429 [hepph]; L. Bergstrom, T. Bringmann and J. Edsjo, arXiv:0808.3725 [astroph]; P. f. Yin, Q. Yuan, J. Liu, J. Zhang, X. j. Bi and S. h. Zhu, arXiv:0811.0176 [hepph]; D. J. Phalen, A. Pierce and N. Weiner, arXiv:0901.3165 [hepph].
 (71) M. Pohl, arXiv:0812.1174 [astroph]; D. Hooper, A. Stebbins and K. M. Zurek, arXiv:0812.3202 [hepph]; D. Hooper, J. E. Taylor and J. Silk, Phys. Rev. D 69, 103509 (2004) [arXiv:hepph/0312076]; D. T. Cumberbatch and J. Silk, Mon. Not. Roy. Astron. Soc. 374, 455 (2007) [arXiv:astroph/0602320].