Reconcile the AMS-02 positron fraction and Fermi-LAT/HESS total spectra by the primary electron spectrum hardening
The recently reported positron fraction up to GeV by AMS-02 seems to have tension with the total electron/positron spectra detected by Fermi and HESS, for either pulsar or dark matter annihilation/decay scenario as the primary positron sources. In this work we will show that the tension will be removed by an adjustment of the primary electron spectrum. If the primary electron spectrum becomes harder above GeV, similar as the cosmic ray nuclei spectrum, the AMS-02 positron fraction and Fermi/HESS data can be well fitted by both the pulsar and dark matter models. This result may be suggestive of a common origin of the cosmic ray nuclei and the primary electrons. Furthermore, this study also implies that the properties of the extra sources derived from the fitting to the AMS-02 data should depend on the form of background.
The AMS-02 collaboration reported the very precise measurement of the positron fraction with energies up to GeV recently Aguilar et al. (2013). The positron fraction shows a continous increase up to GeV, which is consistent with the previous PAMELA result Adriani et al. (2009, 2010) and is lower than that measured with Fermi-LAT Ackermann et al. (2012). A flattening of the positron fraction above GeV is revealed by the AMS-02 data, for the first time. The AMS-02 data implies that there is excess of positrons above tens of GeV compared with the standard cosmic ray (CR) background, and the amount of excess positrons should be less than previously estimated according to the PAMELA data.
Several works appears to explain the AMS-02 data with pulsars
Linden and Profumo (2013) or dark matter (DM) scenarios
Kopp (2013); De Simone et al. (2013); Ibe et al. (2013).
A thorough study of the properties of the extra positron sources,
including the astrophysical one like pulsars and the dark matter
(DM) scenario, based on the AMS-02 data and the electron (or
) spectra measured by PAMELA Adriani
et al. (2011a),
Fermi-LAT Abdo et al. (2009); Ackermann et al. (2010) and HESS
Aharonian et al. (2008, 2009), was given shortly
after the publication of the AMS-02 data Yuan et al. (2013)
(Paper I). It was found that there was difficulty to fit the AMS-02
positron data and the Fermi-LAT/HESS total electron spectra
simultaneously, either in the pulsar scenario or in the DM scenario.
The results seem to imply that there might be tension between
the AMS-02 data and the Fermi-LAT/HESS data, in the present
theoretical framework. This conclusion has been confirmed by other
Several possibilities to reconcile these two data sets were
discussed in Paper I, including multiple components of the extra
sources and the existence of spectral hardening of the primary
electron spectrum. The idea to introduce a spectrum hardening
of primary electrons to fit the data was also raised in
Feng et al. (2014); Cholis and Hooper (2013). The spectrum hardening
was stimulated by the observed spectral hardening of the nuclei
spectra in recent years by several collaborations
Models to explain the spectral hardening include the multi-component sources Zatsepin and Sokolskaya (2006); Yuan et al. (2011); Erlykin and Wolfendale (2012); Thoudam and Hörandel (2012); Bernard et al. (2013), non-linear acceleration of the particles Ptuskin et al. (2013), or the propagation effect Vladimirov et al. (2012); Tomassetti (2012). In Feng et al. (2014) the authors pointed out that if there was a spectral hardening of the primary electron spectrum, there would be a less steep increase (or decrease) of the positron fraction above GeV.
In the work we investigate in detail whether to involve such a
spectral hardening of the primary electron spectrum can help
eliminate the tension between the AMS-02 data and the
Fermi-LAT/HESS data. We employ the CosRayMC tool developed in
Liu et al. (2012) to fit the observational data within
the high dimensional parameter space. The GALPROP
In the next Section we simply describe the models to fit the data. The results are presented in Sec. III. In Sec. IV we give the discussions and conclusions.
In this section we describe the major aspects of the theoretical models to reproduce the electron/positron data briefly. The injection spectra of the primary protons (heavier nuclei are less important in this study) and electrons are both assumed to be broken power-law functions with respect to momentum
where represents the low energy break, is the high energy break to be responsible for the spectral hardening, , and are the spectral indices in different momentum ranges. We also employ the log-parabolic function to describe the spectral hardening of the electrons, i.e., for . In this case is not used. The absolute fluxes of protons and electrons are determined through normalizing the propagated fluxes to normalization factors and .
The background positrons are expected to be produced through the collision of CR nuclei with the interstellar medium (ISM) during the propagation. The parameterization of collision in Kamae et al. (2006) is employed to calculate the secondary production of positrons and electrons. Similar as done in Paper I, we further introduce a free factor to adjust the absolute fluxes of the secondary positrons and electrons to fit the data. Such a factor may represent the uncertainties of the hadronic interactions, propagation models, the ISM density distributions, and the nuclear enhancement factor from heavy elements.
In the PAMELA era it was found that the background contribution are not enough to explain the observed positron fraction and total data Fan et al. (2010); Serpico (2012). Therefore the extra sources of beyond the traditional CR background are introduced to explain the data. We will base on the same theoretical framework to fit the AMS-02 data in the work, assuming continuously distributed pulsars or the DM annihilation/decay to be the extra sources of .
The injection spectrum of from pulsars is assumed to be power-law with an exponential cutoff
where is the normalization factor, is the spectral index and is the cutoff momentum. The spectral index is limited in the range to according to the -ray observations of pulsars Thompson et al. (1994). The spatial distribution of pulsars is taken to be the cylindrically symmetric form given in Lorimer (2004)
where kpc is the distance of the solar location to the Galactic center, kpc is the characteristic height of the Galactic disk.
As for DM scenario (taking annihilation as illustration), we focus on the leptonic two-body annihilation channels and , as implied according to the PAMELA and Fermi-LAT/HESS data of the electrons/positrons and the antiprotons Cirelli et al. (2009); Yin et al. (2009); Meade et al. (2010). The positron/electron production function from DM annihilation is (assumed to be Majorana particles)
where is the mass of DM particle, is the velocity weighted annihilation cross section, is the yield spectrum for one annihilation of a pair of DM particles, and is the DM density profile. The spatial profile of DM energy density is taken to be Navarro-Frenk-White (NFW, Navarro et al. (1997)) distribution
with parameters kpc and GeV cm.
For low energy particles we further employ a simple force field approximation to take into account the solar modulation effect Gleeson and Axford (1968). Since the operation period of PAMELA and AMS-02 is close to the solar minimum, the modulation potential is required to be smaller than GV. Note, however, the low energy part of the positron fraction measured by PAMELA and AMS-02 might not be reproduced with such single solar modulation model, and more complicated charge-sign dependent modulation effect is necessary Della Torre et al. (2012); Maccione (2013).
We first determine the parameters of the proton injection spectrum through fitting to the PAMELA Adriani et al. (2011b) and CREAM Ahn et al. (2010) data. For CREAM data we include systematic uncertainties as discussed in Ahn et al. (2010). The high energy break is fixed to be GeV as suggested by the PAMELA data. The best fitting parameters of the proton spectrum are: , , , GeV, and solar modulation potential MV. The normalization of the propagated proton flux at GeV is cmssrMeV. Comparison of the best fitting spectrum of protons with the observational data is shown in Fig. 1. We see very good agreement between the calculated spectrum and the data. The minimum value is about for degree of freedom (dof).
Since the observational period of protons by PAMELA is almost the same with that of electrons by PAMELA and positrons by AMS-02, we should expect a common modulation amplitude for these particles (besides the charge-sign dependent effect). Therefore we employ a prior on the modulation potential MV comes from the fit of the proton data.
We first fix the electron second break energy at GeV, same as that of protons. The best fitting results of the positron fraction and electron spectrum are shown in Fig. 2. The fitting parameters and the value are presented in Table 2. Compared with the case without spectral hardening of the primary electrons (Paper I), the fitting is indeed improved. The value decrease from to with one more parameter. However, the overall fitting is still not satisfactory. We can see from Fig. 2 that when AMS-02 data are well reproduced, the model expectation is lower than the Fermi data, which is similar with the findings in Paper I.
We then relax the break momentum of the electrons and redo the fit. In this case we find the improvement is significantly, as shown in Fig. 3. The parameters are also given in Table 2. The minimum value over dof is about , which implies a rather good fitting. However, the break momentum is required to be about GeV, which is significantly smaller than that of the nuclei. The difference of the spectral indices below and above is about . As a comparison, such a value is measured to be for protons and for Helium Adriani et al. (2011b). Note, for the fit of proton spectrum in a wider energy range, as shown in Fig. 1, the spectral difference is only about .
It is also possible that the spectral hardening is not a break but a smooth hardening instead, as shown in many models Yuan et al. (2011); Ptuskin et al. (2013). We may use a log-parabolic function to approximate the smooth hardening of the electron spectrum. Fig. 4 presents the results of the fit with log-parabolic shape of the primary electron spectrum. We find the fit is also improved, with the minimum value slightly larger than that with free. The fitting parameters are given in Table 2.
From above we see that including a spectral hardening of the primary electron spectrum, both the PAMELA, AMS-02 and Fermi data can be well fitted with a single component of the extra sources. It is a natural expectation that there is a hardening in the primary electron spectrum, given the observed hardening of the CR nuclei. However, the position of the break might be different from that of nuclei. It is a problem needs to be further understood theoretically. If the hardening of the primary electron spectrum can be confirmed, it would be important to understand the origin and acceleration of the Galactic CRs. Since AMS-02 could measure the pure electron spectrum with much higher precision that PAMELA, we give the expected pure electron spectra in Fig. 5 for the above three cases of the hardening. The future AMS-02 data may test the existence and detailed shape of the primary electron spectrum.
Finally we discuss the DM annihilation as the sources of the . The annihilation final states are assumed to be and . The primary electron spectrum is parameterized with Eq. (1), and is allowed to be free in the fit. The fitting results are shown in Figs. 6 and 7, for and final states respectively. The fitting parameters are compiled in Table 2. It is shown that the DM models can give comparable fittings to the data compared with pulsars. The break momentum of the primary electrons, is also similar with that derived in the pulsar scenario, and is smaller than GeV.
The and favored regions on the parameter plane are given in Fig. 8. For channel DM with mass TeV is favored, while for channel the mass is obtained to be TeV. The boost factor of the annihilation cross section compared with the natural value to give the proper relic density is about hundred to thousand. Such results do not differ much from the ones obtained through fitting the PAMELA positron fraction and the Fermi/HESS data Meade et al. (2010).
The exclusion limits on the DM annihilation into and pairs by -rays from the Galactic center (thin lines, Huang et al. (2012)) and the dwarf galaxies (thick lines, Drlica-Wagner (2012)) are also plotted in Fig. 8. The results show that for the channel the Fermi -rays always give very strong constraints on the annihilation cross section. The constraints for the channel is weaker. The Galactic center -rays tend to exclude the parameter space to explain the excesses. However it may suffer from the uncertainties of the density profile of DM in the halo center. The more robust limits from the dwarf galaxies can not exclude the favored parameter region.
Iv Conclusion and discussion
The study of the highly precise data of positron fraction in CRs reported by AMS-02, as well as the pure electron spectrum measured by PAMELA and the total spectra measured by Fermi and HESS, shows that it is difficult to use a single component of the extra sources to explain the excesses (Paper I, Yuan et al. (2013)). In this work we show that an additional break of the primary electron spectrum can improve the fit significantly. The best fitting break momentum is about GeV and the spectral difference is . As a comparison, the break momentum of protons is about GeV, and the spectral difference is . The hardening behavior of the electron spectrum is different from that of nuclei, which makes the understanding of the fine structures of the CR spectra non-trivial.
The different behaviors of the spectral hardening between nuclei and electrons are probably due to the fact that high energy electrons should come from nearby regions, and less number of relevant sources leads to larger fluctuations of the electron spectrum than that of nuclei. It is also possible that, if one or several nearby sources are responsible to the spectral hardening, the accelerated electron-to-proton ratio is higher for these sources.
In the presence of a hardening of the primary electron spectrum, both the pulsar and DM scenarios can give comparable fit to the data. However, the DM scenario are strongly constrained by the -rays, especially for the tauon final state. We would like to point out that it will be equivalent to take the harder part of the electron spectrum and to drop the constraints from the PAMELA electron data. In such ways both the AMS-02 positron fraction and Fermi total spectrum can be fitted simultaneously.
The AMS-02 will measure the electron spectrum with high precision in the near future. Whether there is a hardening in the electron spectrum or a lower total spectrum than Fermi will soon be answered by AMS-02.
Normalization at 25 GeV in unit of
Normalization at 25 GeV in unit of cmssrMeV.
Appendix A Results of the background positrons and electrons
For the convenience of use we tabulate the fluxes of the background positrons and electrons calculated with the best fitting parameters of the pulsar models in Table 3. For DM models the results have little difference. Note for background positrons an additional factor as given in Table 2 needs to be multiplied. Here the local interstellar fluxes are given. If one wants to better reproduce the low energy electron spectrum, the solar modulation with modulation potential given in Table 2 is necessary.
|(GeV)||(GeV m s sr)|
|( GeV)||( free)||(log-parabolic)|
Acknowledgements.This work is supported by 973 Program under Grant No. 2013CB837000 and the National Natural Science Foundation of China under Grant Nos. 11075169, 11105155.
- Also the preliminary data about the total spectra by AMS-02 show the discrepancy with Fermi-LAT data below GeV. For GeV the difference between these two data sets are smaller AMS-02 collaboration (2013). The discussion of the present work may hold given the new AMS-02 spectra, since we study to reconcile the high energy behavior of Fermi-LAT data with the AMS-02 positron fraction.
- Note, however, the recent reported preliminary data about the proton and Helium spectra show no significant hardening as observed by PAMELA AMS-02 collaboration (2013). Whether or not the heavier nuclei have the hardening still needs to be tested by the future AMS-02 data. Since the present work is based on the electron and positron data, the AMS-02 results about the nuclei spectra have only a slight modification of the secondary positron spectrum and do not affect the discusion significantly.
- M. Aguilar, et al., Phys. Rev. Lett. 110, 141102 (2013).
- O. Adriani, et al., Nature 458, 607 (2009), eprint 0810.4995.
- O. Adriani, et al., Astroparticle Physics 34, 1 (2010), eprint 1001.3522.
- M. Ackermann, et al., Phys. Rev. Lett. 108, 011103 (2012), eprint 1109.0521.
- T. Linden and S. Profumo, Astrophys. J. 772, 18 (2013), eprint 1304.1791.
- J. Kopp, Phys. Rev. D 88, 076013 (2013), eprint 1304.1184.
- A. De Simone, A. Riotto, and W. Xue, J. Cosmol. Astropart. Phys. 5, 003 (2013), eprint 1304.1336.
- M. Ibe, S. Iwamoto, S. Matsumoto, T. Moroi, and N. Yokozaki, Journal of High Energy Physics 8, 29 (2013), eprint 1304.1483.
- O. Adriani, et al., Phys. Rev. Lett. 106, 201101 (2011a).
- M. Ackermann, et al., Phys. Rev. D 82, 092004 (2010).
- A. A. Abdo, et al., Phys. Rev. Lett. 102, 181101 (2009), eprint 0905.0025.
- F. Aharonian, et al., Phys. Rev. Lett. 101, 261104 (2008), eprint 0811.3894.
- F. Aharonian, et al., Astron. Astrophys. 508, 561 (2009), eprint 0905.0105.
- Q. Yuan, X.-J. Bi, G.-M. Chen, Y.-Q. Guo, S.-J. Lin, and X. Zhang, ArXiv e-prints:1304.1482 (2013), eprint 1304.1482.
- AMS-02 collaboration, in International Cosmic Ray Conference (http://www.ams02.org/2013/07/new-results-from-ams-presented-at-icrc-2013/, 2013).
- I. Cholis and D. Hooper, Phys. Rev. D 88, 023013 (2013), eprint 1304.1840.
- H.-B. Jin, Y.-L. Wu, and Y.-F. Zhou, J. Cosmol. Astropart. Phys. 11, 026 (2013), eprint 1304.1997.
- L. Feng, R.-Z. Yang, H.-N. He, T.-K. Dong, Y.-Z. Fan, and J. Chang, Physics Letters B 728, 250 (2014), eprint 1303.0530.
- A. D. Panov, et al., Bulletin of the Russian Academy of Science, Phys. 71, 494 (2007), eprint astro-ph/0612377.
- O. Adriani, et al., Science 332, 69 (2011b), eprint 1103.4055.
- H. S. Ahn, et al., Astrophys. J. Lett. 714, L89 (2010), eprint 1004.1123.
- V. I. Zatsepin and N. V. Sokolskaya, Astron. Astrophys. 458, 1 (2006), eprint astro-ph/0601475.
- Q. Yuan, B. Zhang, and X.-J. Bi, Phys. Rev. D 84, 043002 (2011), eprint 1104.3357.
- A. D. Erlykin and A. W. Wolfendale, Astroparticle Physics 35, 449 (2012), eprint 1111.3191.
- S. Thoudam and J. R. Hörandel, Mon. Not. Roy. Astron. Soc. 421, 1209 (2012), eprint 1112.3020.
- G. Bernard, T. Delahaye, Y.-Y. Keum, W. Liu, P. Salati, and R. Taillet, Astron. Astrophys. 555, A48 (2013), eprint 1207.4670.
- V. Ptuskin, V. Zirakashvili, and E.-S. Seo, Astrophys. J. 763, 47 (2013), eprint 1212.0381.
- A. E. Vladimirov, G. Jóhannesson, I. V. Moskalenko, and T. A. Porter, Astrophys. J. 752, 68 (2012), eprint 1108.1023.
- N. Tomassetti, Astrophys. J. Lett. 752, L13 (2012), eprint 1204.4492.
- J. Liu, Q. Yuan, X.-J. Bi, H. Li, and X. Zhang, Phys. Rev. D 85, 043507 (2012), eprint 1106.3882.
- A. W. Strong and I. V. Moskalenko, Astrophys. J. 509, 212 (1998), eprint astro-ph/9807150.
- A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002), eprint astro-ph/0205436.
- S. J. Lin, C. F. Cai, and et al., in preparation (2013).
- T. Kamae, N. Karlsson, T. Mizuno, T. Abe, and T. Koi, Astrophys. J. 647, 692 (2006), eprint astro-ph/0605581.
- Y.-Z. Fan, B. Zhang, and J. Chang, International Journal of Modern Physics D 19, 2011 (2010), eprint 1008.4646.
- P. D. Serpico, Astroparticle Physics 39, 2 (2012), eprint 1108.4827.
- D. J. Thompson, et al., Astrophys. J. 436, 229 (1994).
- D. R. Lorimer, in Young Neutron Stars and Their Environments, edited by F. Camilo & B. M. Gaensler (2004), vol. 218 of IAU Symposium, p. 105.
- M. Cirelli, M. Kadastik, M. Raidal, and A. Strumia, Nuclear Physics B 813, 1 (2009), eprint 0809.2409.
- P. F. Yin, Q. Yuan, J. Liu, J. Zhang, X. J. Bi, S. H. Zhu, and X. M. Zhang, Phys. Rev. D 79, 023512 (2009), eprint 0811.0176.
- P. Meade, M. Papucci, A. Strumia, and T. Volansky, Nuclear Physics B 831, 178 (2010), eprint 0905.0480.
- J. F. Navarro, C. S. Frenk, and S. D. M. White, Astrophys. J. 490, 493 (1997), eprint astro-ph/9611107.
- L. J. Gleeson and W. I. Axford, Astrophys. J. 154, 1011 (1968).
- L. Maccione, Phys. Rev. Lett. 110, 081101 (2013), eprint 1211.6905.
- S. Della Torre, et al., Advances in Space Research 49, 1587 (2012).
- J. Alcaraz, et al., Phys. Lett. B 490, 27 (2000).
- T. Sanuki, et al., Astrophys. J. 545, 1135 (2000), eprint astro-ph/0002481.
- M. Aguilar, et al., Phys. Lett. B 646, 145 (2007), eprint astro-ph/0703154.
- S. W. Barwick, et al., Astrophys. J. Lett. 482, L191 (1997), eprint astro-ph/9703192.
- S. Coutu, et al., in International Cosmic Ray Conference (2001), vol. 5 of International Cosmic Ray Conference, p. 1687.
- J. Chang, et al., Nature 456, 362 (2008).
- X. Huang, Q. Yuan, P.-F. Yin, X.-J. Bi, and X. Chen, J. Cosmol. Astropart. Phys. 11, 048 (2012), eprint 1208.0267.
- A. Drlica-Wagner, in Fermi Symposium 2012 (2012).