
2
Previously, we presented an updated compendium of
constraints on particle DM annihilation to neutrinos [43].
Here, we turn our attention to the production of neutri-
nos from DM decay, directing special attention to the
higher-mass region. At masses greater than ∼TeV, elec-
troweak (EW) corrections can“decloak” the DM, produc-
ing high-energy photons, even when directly decaying to
neutrinos that can be detectable at current and future
gamma-ray observatories [44]. We thus present limits
from current measurements and sensitivities for upcom-
ing experiments covering a DM mass range starting at
20 MeV and spanning well into the ultra-heavy domain
up to 1011 GeV. Although recent LHAASO results are on
par with IceCube’s recent analyses, we will find that cur-
rent and future neutrino telescopes retain superior sen-
sitivity across nearly the entire range of masses that we
consider. This is due to three factors: the loop suppres-
sion of the gamma-ray production rate, the growth of
the electroweak cross section with energy, and the loss of
high-energy gamma rays to interactions in the interstellar
and intergalactic medium.
We begin by briefly describing the signature of DM
decaying to neutrinos at both neutrino telescopes and
gamma-ray observatories. We present our new results
in Section III and offer parting words of wisdom in Sec-
tion IV.
II. DARK MATTER DECAYS TO NEUTRINOS
The expected flux per neutrino flavor at Earth from
decay of DM with mass mχand lifetime τχis
dΦν+¯ν
dEν
=1
4π
1
τχmχ
1
3
dNν
dEν
D(Ω).(1)
Below the electroweak scale, the neutrino spectrum per
decay is dNν/dEν= 2δ(1 −2E/mχ)mχ/2E2; at higher
masses a low-energy tail arises as discussed in Ref. [44].
Because relevant backgrounds follow a power-law distri-
bution, only the delta contribution is relevant for neu-
trino constraints. The so-called D-factor, D(Ω), is an
integral of the DM distribution ρ(x) along the line of
sight and solid angle ∆Ω:
D≡ZdΩZl.o.s.
ρχ(x)dx. (2)
We assume the Galactic DM spatial distribution is mod-
eled by an Navarro-Frenk-White (NFW) profile with a
slope parameter γ= 1.2 and a scale radius rs= 20 kpc,
and we set the local DM density to ρ0= 0.4 GeV cm−3.
These parameters are consistent with the results of, e.g.,
Ref. [45], which point out a strong dependence on how
the baryonic potential is modeled. We take the distance
to the galactic centre to be R0= 8.127 kpc [46]. Here, we
mainly strive to make results as self-consistent as possible
by using common halo parameters in all of our analyses.
We have assumed equal production of each flavor, which
leads to equal flavors at Earth. Due to neutrino oscil-
lation, this will remain approximately true regardless of
the initial flavor composition.
The D-factor depends on the field of view of each ex-
periment. Effective areas are reported as a function of
elevation (or equivalently, zenith angle). Given each ex-
periment’s latitude and altitude, we integrate these ac-
ceptances over a period of 24 hours, where the solid angle
integral is weighted by the fractional acceptance. This
defines an effective D-factor:
Deff =Zdt ZdΩZl.o.s.
ρχ(x)F(Ω, t)dx, (3)
where F(Ω, t) is the fractional acceptance in equatorial
coordinates. This procedure is simplified for experiments
at the South Pole (IceCube, ANITA), where elevation
and declination are equivalent.
In computing the yield of dark matter in a given ex-
periment, we convolve the experimental efficiency with
flux from neutrinos from a given direction. For our back-
ground agnostic constraints, the flux obtained by this
procedure is then compared to the unfolded neutrino
fluxes. When published experimental results on dark
matter annihilation use a directional analysis, such as the
case of ANTARES, we rescale the result by the efficiency-
weighted ratio of the dark matter Jto Dfactors; namely
the ratio of expected signal yield in the annihilation to
decay scenarios. Efficiencies used for each experiment are
given in Suppl. Table III.
Based on Eq. 1, we produce limits on the DM de-
cay rate, using results from different analyses of exist-
ing data [28,43,50,54,61,73–83] or forecasted sensi-
tivities [84–89]. The full list of neutrino experiments is
given in the top section of Table I. We also list the neu-
trino energy range covered by each experiment, spanning
from 10 MeV at Borexino to >1011 GeV at IceCube and
AUGER, as well as each experiment’s neutrino flavor sen-
sitivity. For a detailed description of each experiment and
its sensitivity, we point the reader to Ref. [43]. The D-
factors for each experiment are computed by integrating
the exposure of each telescope over 24 hours. The result-
ing exposures and D-factors are tabulated in Table III in
the Supplementary Material.
The decay lifetime constraints result from a com-
parison between the flux sensitivities from each exper-
iment and the expected neutrino flux from DM de-
cay. This approach assumes a branching ratio of 100%
of DM decay to neutrinos can describe the total neu-
trino flux measurements in the Galactic Center region.
Our forecasts assume five years of exposure for each
of the following experiments: JUNO [49], DUNE [51],
Hyper-Kamiokande (HK) [51], RNO-G [60], IceCube-
Gen2 [59], KM3NeT [57], P-ONE [56], TAMBO [58], and
GRAND [61]. Constraints that are not derived by us,
but are reported by experiments or other groups, are
rescaled to match the D-factors used in this work. This
enables a fair comparison between different experimental
constraints.