different quark or gluon final states in 2-body spectrum we use PPPC4DMID [28] to gener-
ate the differential yield which is weighted by the corresponding annihilation fraction, i.e.
hσ
ann
vi
¯qq
/hσ
ann
vi and hσ
ann
vi
gg
/hσ
ann
vi. The 4-body cosmic ray spectrum is given by the
spectrum of the mediator decay in its rest frame followed by a Lorentz boost [12, 29]. The
spectrum is then weighted by the product of annihilation fraction and decay branching ratio
of the mediator, i.e.
hσ
ann
vi
Med
hσ
ann
vi
Γ
Med→q ¯q
Γ
Med
and
hσ
ann
vi
Med
hσ
ann
vi
Γ
Med→gg
Γ
Med
, to give (dN
i
/dE)
4−body
[30].
The expressions of mediator decay widths and dark matter annihilation cross sections are
collected in appendix.
3 Indirect observables from AMS-02 and Fermi-LAT
In this section we describe the observables of antiproton flux and gamma ray measured by
AMS-02 and Fermi-LAT respectively.
3.1 Antiproton flux from AMS-02
The two key unknowns about cosmic rays in the Galaxy are their origin and propagation.
The propagation of cosmic rays can be described as the process of diffusion. The diffusion
process is written in the form of the transport equation below [31]
∂ψ
∂t
= Q(~r, p) +
~
∇ ·
D
xx
~
∇ψ −
~
V ψ
+
∂
∂p
p
2
D
pp
∂
∂p
1
p
2
ψ
−
∂
∂p
h
˙pψ −
p
3
~
∇ ·
~
V
ψ
i
−
ψ
τ
f
−
ψ
τ
r
, (3.1)
where ψ(~r, t, p) is the density of cosmic rays,
~
V is the convection velocity and τ
f
(τ
r
) is
the time scale for fragmentation (radioactive decay). ˙p is the momentum loss rate. The
convection terms in the above equation are induced by the Galactic wind. The diffusion in
momentum space governs the reacceleration process. In this case the diffusion coefficient in
momentum space, i.e. D
pp
, is related to the spatial coefficient D
xx
and the Alfven velocity
v
A
[32]:
D
pp
D
xx
=
4p
2
v
2
A
3δ(4 − δ
2
)(4 − δ)w
, (3.2)
with the level of the interstellar turbulence parameter w being 1. The spatial diffusion
coefficient is usually written in this form
D
xx
= βD
0
(R/R
0
)
δ
, (3.3)
with R and β being the rigidity and particle velocity divided by light speed respectively.
This transport equation is numerically solved based on given boundary conditions, that
is, the cosmic ray density ψ vanishes at the radius R
h
and the height z
0
of the cylindrical
diffusion halo.
In eq. (3.1), the source term can be written by the product of the spatial distribution
and the injection spectrum function for cosmic ray species i
Q
i
(~r, p) = f(r, z)q
i
(p). (3.4)
– 4 –