
2
i.e. the fist generation of Standard Model particles pro-
duced from the annihilation event. Dark matter origi-
nated features may be as small as a localized in energy,
to give a few %bump on the positron fraction, or as
wide in energy and large in amplitude as the entire rising
above 5 GeV positron fraction spectrum.
Inversely, using the AMS-02 positron fraction’s rela-
tively smooth spectrum, one can set upper limits on the
annihilation cross section of dark matter particles. That
is done for a range of masses and a variety of annihilation
channels [82–84]. This is the main aim of this paper.
The origin of the rising above 5 GeV positron frac-
tion spectrum that was first measured by the Payload
for Antimatter Matter Exploration and Light-nuclei As-
trophysics (PAMELA) satellite [85, 86], then confirmed
by the Fermi-LAT [87] and further measured with an un-
precedented accuracy by AMS-02 [75, 88, 89], has been
a subject of great interest. One explanation for the ad-
ditional positron flux, is relatively close-by “young” and
“middle-aged” Milky Way pulsars that during their pul-
sar wind nebula (PWN) phase converted an apprecia-
ble fraction (O(0.01) −O(0.1)) of their rotational energy
into high-energy cosmic-ray electrons (e−) and positrons
(e+) [90–108]. Another explanation is Milky Way SNRs,
that in their first O(10) kyr produced and accelerated
secondary cosmic rays including positrons [109–117] (see
however [118–121]). Furthermore, detailed modifications
on the distribution of cosmic-ray sources and the prop-
agation of cosmic rays through the ISM [122–124] and
annihilating or decaying dark matter models have been
explored to explain the positron fraction measurement
[101, 125–144]. We assume in this work that the over-
all rise of the positron fraction, shown with its AMS-02
measurement in Fig. 1, is not caused by dark matter, but
instead from a more conventional source; a population of
Milky Way pulsars.
Pulsars are localized sources of cosmic-ray electrons
and positrons. Due to their rapid spin-down, pulsars
convert their initial rotational energy into cosmic-ray e±
and subsequently release those e±into the ISM in a com-
paratively short amount of time1. That makes pulsars
cosmic-ray e±sources approximately localized both in
space and time. High-energy e±lose rapidly their en-
ergy through synchrotron radiation and inverse Comp-
ton scattering as they interact with the ISM and before
reaching us. That results in an upper energy cut-off, on
the e±spectra from individual pulsars [95, 96, 98]. In
turn, a population of pulsars that could collectively ex-
plain the rising positron fraction spectrum, could also
give spectral features at the higher energies where the
number of contributing pulsars is reduced to the point
1The time required for most cosmic-ray e±produced around the
PWN environments to be released into the ISM, is at least an or-
der of magnitude smaller than the propagation time required for
these cosmic rays to reach our detectors [96]. The only exception
would be a very close (O(10) pc) pulsar (see however [145]).
of individual sources dominating narrow parts of that
spectrum [96, 98, 104, 108]. Such features can then be
searched for as in [146]. Similar arguments can be made
for PWNe. However, their expected higher energy cut-
offs are less sharp by comparison [117]. We use modeled
populations of Milky Way pulsars produced in our ear-
lier work of [108]. In Ref. [108], a library of publicly
available pulsar population models was created that is in
agreement with the cosmic-ray e±flux spectral measure-
ments from AMS-02 [75, 76], the CALorimetric Electron
Telescope (CALET ) [147] and the DArk Matter Parti-
cle Explorer (DAMPE ) telescope [148], as well as the
AMS-02 positron fraction spectrum [89]. As the pulsar’s
contribution to the positron fraction spectrum is not per-
fectly smooth and with uncertainties, we use a library of
models instead of just one generic parameterization. As
we will show, we derive more conservative and more real-
istic limits on the dark matter annihilation cross section.
In Section II, we discuss the general methodology of
our approach, including the observations that we use, the
astrophysical background modeling of the positron frac-
tion and the statistical treatment followed in fitting the
data. We also create mock positron fraction data to an-
swer the question on the robustness of the positron frac-
tion measurement as a means to study the particle prop-
erties of dark matter. Then in Section III, we present the
results of searching for a possible dark matter signal in
the positron fraction. We find that the limits on the anni-
hilation cross section are not well defined. That is due to
the underlying astrophysical background uncertainties.
The annihilation cross section limits have a width that is
at least one order of magnitude in the mas range of 5 to
120 GeV that we study. In addition, we find indications
for a possible excess of 5-15 GeV in cosmic-ray energy e±.
That excess while compatible with a WIMP-scale dark
matter signal, has a significance that varies with the as-
trophysical background modeling and is not claimed to
be a robust one. Further scrutiny will be required as
cosmic-ray physics in that energy range improve with fu-
ture observations. Moreover, in Section III, we perform
our mock positron fraction analysis. We find that if dark
matter contributes to the positron fraction spectrum at
the few percent level within an range spanning several
AMS-02 energy bins, such an excess signal can not be
absorbed by the astrophysical background uncertainties.
However, identifying the exact particle properties of the
dark matter particle responsible for that excess is a a
more model-dependent inquiry. Finally, in Section IV,
we give our conclusions and discuss connections to other
types of dark matter searches as well as future prospects.
II. METHODOLOGY
In this section, we describe the energy range of the
AMS-02 positron fraction (e+/(e++e−)) measurement
used in this analysis. We also explain how we construct
our background astrophysical models, which are fitted