Pith. Sign in

REVIEW

Spectra of neutrinos from dark matter annihilations

High Energy Physics - Phenomenology · Submitted 29 Jun 2005 · 53 references

Pith's verdict: A solid, standard-reference phenomenology paper: full-flavor spectra plus a combined oscillation/absorption/regeneration treatment, with the soft spots confined to modeling systematics and idealized detector assumptions.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0506298 (version 5)
Pith IDLS6XR4MC, Pith's permanent record of this paper and its review
CitationAlessandro Strumia, Francesco Vissani, Igor Sokalski, Marco Cirelli, Nicolao Fornengo, Teresa Montaruli. "Spectra of neutrinos from dark matter annihilations." arXiv:hep-ph/0506298 (2005). https://pith.science/paper/LS6XR4MC
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Dark-matter neutrino spectra, not their total rates, may carry the particle's mass and annihilation channels, because spectral shape survives the astrophysical uncertainties that swamp normalization.

The reading

The paper asks whether a dark-matter signal in a neutrino telescope can be more than a detection: whether the shape of the neutrino energy spectrum, rather than its overall size, carries information about the dark matter particle itself. The authors' answer is yes, because the total rate scales with the local dark matter density and the capture rate, quantities known only to within roughly an order of magnitude, while the spectral shape depends on the dark matter mass and on the branching ratios into $b$-quarks, taus, $W$ and $Z$ bosons and the rest. To make that case they compute neutrino spectra at production for all three flavors and the main annihilation channels, separately for the Earth and for the Sun, then propagate them through matter with a density-matrix formalism that treats oscillations, charged- and neutral-current absorption, and $nu_tau$ regeneration together rather than one at a time. They translate the resulting fluxes into the event classes a real detector records, through-going muons, contained muons and showers, and show with a two-parameter fit how well a modest event sample could recover the mass and the mix of a tau channel with a direct neutrino line. The reframing matters because it turns dark-matter neutrino searches from counting exercises into spectroscopy.

What carries the argument

The load-bearing object is the evolution equation for the neutrino and antineutrino density matrices, $d rho/dr = -i[H,rho] + (d rho/dr)_{rm CC} + (d rho/dr)_{rm NC} + (d rho/dr)_{rm in}$, where $H$ is the vacuum-plus-MSW Hamiltonian and the remaining terms encode absorption, neutral-current reinjection at lower energy, and charged-current $nu_tau$ regeneration through tau decay. A second, independent piece of machinery is the analytic 'limit spectrum': under the approximations that oscillations roughly equidistribute flavors, that reinjection spectra are flat in the energy ratio, and that cross sections grow linearly with energy, the full system reduces to $E, partial rho/partial x = -E rho + int_E^infty rho,dE'$, whose solution is an exponential in energy whose slope is set by the solar column density and the weak cross sections, not by the initial annihilation channel. That solution is what licenses the claim that very heavy dark matter produces a channel-blind solar neutrino spectrum.

Extended reading notes

Core claim

The central claim is that dark-matter neutrino spectra are effectively free of astrophysical uncertainty even though the corresponding rates are not, and that this makes the spectra a usable measurement channel. The authors support it with three pieces of machinery. First, full production spectra for $nu_e, nu_mu, nu_tau$ from $nu bar nu$, $b bar b$, $tau^+tau^-$, $c bar c$, light quarks, gluons, $t bar t$, $ZZ$ and $W^+W^-$, computed with a modified PYTHIA that stops hadrons and leptons which would otherwise lose most of their energy in the dense medium before decaying. Second, a propagation treatment that evolves $3 times 3$ density matrices, so oscillations, absorption and regeneration act simultaneously instead of being layered on separately. Third, a demonstration that for very heavy dark matter the outgoing solar spectrum collapses to a nearly universal exponential, proportional to $e^{-E_nu/bar E}$ with $bar E approx 100$ GeV for neutrinos and $approx 140$ GeV for antineutrinos, erasing the initial channel information. The propagation corrections are large, a factor of order $0.1$ to $10$ on total rates, yet they distort the spectra only mildly, and that gap is what allows the reconstruction program. The authors' own version note records that later corrections moved the numerical spectra by 10 to 20 percent while leaving every physics conclusion unchanged.

Load-bearing premise

Everything downstream rests on how the authors model the energy lost by hadrons and leptons inside the dense solar and terrestrial cores, in particular the heavy-hadron energy-loss formula and its exponential-integral approximation, and the paper attaches no systematic error bar to that modeling.

Editorial extensions

If this is right

  • If spectral shapes really are astrophysics-free, a detected excess can be fitted for the dark matter mass and annihilation branching ratios with of order 100 contained muons or 200 showers, converting a discovery into a particle-physics measurement.
  • Oscillations raise the muon rate from the $tau^+tau^-$ channel by up to a factor of about 3 for masses below 100 GeV, so the usual expectation that $b bar b$ dominates the muon signal has to be revised channel by channel.
  • For dark matter much heavier than a TeV, solar spectra become universal exponentials and the Sun measures only the annihilation rate, not the annihilation channel.
  • Light quarks and gluons contribute non-negligibly at low neutrino energies for masses above roughly 500 GeV, so precise experimental limits must include channels usually dropped.
  • Showers retain the primary neutrino spectral features better than through-going muons and have about twice the rate of contained muons, which makes them the more informative topology if a detector can measure their energy.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The shape-robust, normalization-fragile logic is not specific to neutrinos: the same split between spectral shape and overall rate should apply to gamma-ray spectra from annihilation in the galactic center, where local density uncertainty similarly afflicts rates but not line or bump positions.
  • If the universal exponential cutoff near 100 GeV is real, then a solar neutrino signal with that cutoff would be evidence for very heavy dark matter almost independently of its annihilation mode, which is a discovery channel that does not require guessing the branching ratios first.
  • The channel-discrimination power quoted for 100 contained muons presumes roughly 30 GeV energy bins with Gaussian resolution, so a next-generation megaton water Cherenkov detector with Super-Kamiokande-like resolution is the realistic target, and the paper's own numbers imply that current telescope granularity, not statistics, is the binding constraint.
  • Because the hadronic energy-loss prescription feeds directly into the solar production spectra, an independent transport benchmark of that step would be the cheapest way to harden the branching-ratio reconstruction against a systematic the paper does not quantify.
Original abstract

We study the fluxes of neutrinos from annihilations of dark matter particles in the Sun and the Earth. We give the spectra of all neutrino flavors for the main known annihilation channels: nu-antinu, b-bbar, tau-taubar, c-cbar, light quarks, ZZ, W^+W^-. We present the appropriate formalism for computing the combined effect of oscillations, absorptions, nu_tau-regeneration. Total rates are modified by an O(0.1--10) factor, comparable to astrophysical uncertainties, that instead negligibly affect the spectra. We then calculate different signal topologies in neutrino telescopes: through-going muons, contained muons, showers, and study their capabilities to discriminate a dark matter signal from backgrounds. We finally discuss how measuring the neutrino spectra can allow to reconstruct the fundamental properties of the dark matter: its mass and its annihilation branching ratios.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,996 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.57 · human participants 0.01 · controlled comparison 0.08 · randomized assignment 0.01 · states quantitative result 0.84 · compares with prior results 0.15 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v5. Results and the main claim are underlined; .

We study the fluxes of neutrinos from annihilations of dark matter particles in the Sun and the Earth. We give the spectra of all neutrino flavors for the main known annihilation channels: nu-antinu, b-bbar, tau-taubar, c-cbar, light quarks, ZZ, W^+W^-. We present the appropriate formalism for computing the combined effect of oscillations, absorptions, nu_tau-regeneration. Total rates are modified by an O(0.1--10) factor, comparable to astrophysical uncertainties, that instead negligibly affect the spectra. We then calculate different signal topologies in neutrino telescopes: through-going muons, contained muons, showers, and study their capabilities to discriminate a dark matter signal from backgrounds. We finally discuss how measuring the neutrino spectra can allow to reconstruct the fundamental properties of the dark matter: its mass and its annihilation branching ratios.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We study the fluxes of neutrinos from annihilations of dark matter particles in the Sun and the Earth. contribution 0.96 0.88 0.05 asserted 0.78
s1 We give the spectra of all neutrino flavors for the main known annihilation channels: nu-antinu, b-bbar, tau-taubar, c-cbar, light quarks, ZZ, W^+W^-. contribution 0.80 0.90 0.09 asserted 1.00
s2 We present the appropriate formalism for computing the combined effect of oscillations, absorptions, nu_tau-regeneration. method 0.89 0.85 0.06 asserted 0.99
s3 Total rates are modified by an O(0.1--10) factor, comparable to astrophysical uncertainties, that instead negligibly affect the spectra. result 0.96 0.72 0.19 asserted 1.00
s4 We then calculate different signal topologies in neutrino telescopes: through-going muons, contained muons, showers, and study their capabilities to discriminate a dark matter signal from backgrounds. contribution 0.92 0.94 0.09 asserted 0.96
s5 We finally discuss how measuring the neutrino spectra can allow to reconstruct the fundamental properties of the dark matter: its mass and its annihilation branching ratios. interpretation 0.83 0.65 0.50 hedged possibility 0.74

What kind of paper the classifier reports this is: theoretical result 0.73 · evidence from simulation or computation 0.94.

What was examined, and what was not

sources read
  • arXiv listing abstract v5
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. The left plot illustrates how oscillations and CC absorption separately affect a flux of neutrinos produced in the center of the Sun. The right plot shows the oscillation probabilities from the center of the Earth. The continuous line applies to ν for θ13 = 0 and to ν¯ for any allowed θ13, since matter effects suppress their mixing. The dotted line applies to ν for θ13 = 0.1 rad. The average over the production point has been performed as appropriate for mDM = 100 GeV. It is responsible for the damping effect visible at Eν <∼ 10 GeV. muon and tau neutrinos at production point: the…
Figure 2. Neutrino spectra at production. Upper half: the fluxes of electron and muon neutrinos, for the seven main annihilation channels and for different masses of the parent DM particle (different colors). The solid lines apply to the case of the Sun, the dotted of the Earth. In all cases, the spectra of antineutrinos are the same as those of neutrinos. Lower half: the same for ( ν ) τ . 9
Figure 3. Energy distributions of ν (red) and ν¯ (blue) produced with energy E′ ν by one NC DIS interaction of a ( ν ) with energy Eν. The energy is plotted in units of Eν. Continuous line: in normal matter, where Np ≈ Nn. Dotted line: around the center of the Sun, where Np ≈ 2Nn. composed of an average of the initial νµ and ντ fluxes. This is exactly the case for neutrinos with very small Eν. At larger energies, the damping effect is milder and indeed one can follow the fast oscillations. In fig. 1b the result is well visible. In short: the spatial average over the slightly different…
Figure 4. Energy distributions of neutrinos regenerated by CC scatterings of a ( ν ) τ of given energy Eντ , produced by one ( ν ) τ/nucleon scattering. The blue upper curves are fτ→τ (Eντ , E′ ν ), and the red lower curves are fτ→e,µ(Eντ , E′ ν ), plotted for several values of the incident ντ energy Eντ . from the flux and produces an almost collinear charged lepton. The τ ± produced by ( ν ) τ decays promptly, before loosing a significant part of its energy in interactions with the surrounding matter, and therefore re-injects secondary fluxes of energetic neutrinos [37, 38]: ντ → τ − → X ντ…
Figure 5. Neutrino spectra generated by one DM annihilation around the center of Earth. The plots show the spectra of the three neutrino flavors (the three rows) and assume different DM masses (the three columns). Each plot shows the open annihilation channels DM DM → b ¯b, τ +τ −, cc¯, tt¯, W+W−, ZZ. The DM DM → νν¯ channel (not shown) would produce a line at Eν = mDM. The dotted lines show the spectra without oscillations while solid lines are the final results after oscillations. The dashed lines in the upper-left panel have been computed with θ13 = 0.1 rad for illustration (see text); all
Figure 6. Spectra of through-going µ + and µ − (summed) generated by DM annihilations around the center of Earth. The plots assume the DM annihilation rate of eq. (21), different DM masses and show the main annihilation channels. For better illustration, some channels have been rescaled by the indicated factor. The shaded region is the atmospheric background. Fully contained µ ± from the Earth 3 10 30 µ ± energy in GeV 0 50 100 150 200 d N / dlog10 Eµ mDM = 50 GeV ν / 4 4 b τ 3 10 30 100 µ ± energy in GeV 0 20 40 60 80 100 d N / dlog10 Eµ mDM = 100 GeV ν / 4 4 b τ Z W 10 30 100 µ ± energy in…
Figure 7. Spectra of fully contained µ + and µ − (summed) generated by DM annihilations around the center of Earth. We assume the DM annihilation rate of eq. (21) and a detector with Mton·year exposure. The νν¯ channel gives a µ ± spectrum peaked at Eµ ∼ mDM, due to the monochromatic parent spectrum. Showers from the Earth 3 10 30 Shower energy in GeV 0 100 200 300 400 d N / dlog10 Eshower mDM = 50 GeV ν / 4 4 b τ 10 30 100 Shower energy in GeV 0 50 100 150 200 250 d N / dlog10 Eshower mDM = 100 GeV ν / 4 4 b τ Z W 10 30 100 Shower energy in GeV 0 20 40 60 80 100 120 140 d N / dlog10 Eshower…
Figure 8. (Idealized) energy spectra of showers generated by DM annihilations around the center of Earth, in a detector with Mton·year exposure and unit detector efficiency, assuming the DM annihilation rate of eq. (21). 24
Figure 9. Neutrino fluxes generated by DM annihilations around the center of Sun. Upper row: νe fluxes. Lower row: the almost equal fluxes of νµ and ντ . We plot the combination (2Φν + Φν¯)/3. All spectra are significantly different from those at production point (not shown here). For the DM DM → νν¯ channel at Eν = mDM we plotted the survival probability rather than the energy spectrum, because it is not possible to plot a Dirac δ function. We present our results showing the same kinds of plots previously employed in the Earth case. Fig.s 9 show the DMν fluxes3 . The main topologies of…
Figure 10. Spectra of through-going µ + and µ − (summed) generated by DM annihilations around the center of Sun. The plots assume the DM annihilation rate of eq. (24), different DM masses and show the main annihilation channels. Some channels have been rescaled by the indicated factor for better illustration. Fully contained µ ± from the Sun Ν 4 4 b Τ 3 10 30 0 20 40 60 80 100 120 140 Μ ± energy in GeV d N ‘dlog10 EΜ mDM = 50 GeV Ν 4 4 b Τ Z W 3 10 30 100 0 10 20 30 40 50 60 Μ ± energy in GeV d N ‘dlog10 EΜ mDM = 100 GeV Ν 4 b Τ t Z W 10 30 100 300 0 1 2 3 4 5 6 Μ ± energy in GeV d N…
Figure 11. Spectra of fully contained µ + and µ − (summed) generated by DM annihilations around the center of Sun. We assume the DM annihilation rate of eq. (24) and a detector with Mton·year exposure. Showers from the Sun Ν 4 4 b Τ 3 10 30 0 100 200 300 400 Shower energy in GeV d N ‘dlog10 Eshower mDM = 50 GeV Ν 4 4 b Τ Z W 10 30 100 0 50 100 150 200 Shower energy in GeV d N ‘dlog10 Eshower mDM = 100 GeV Ν 4 b Τ t Z W 10 30 100 300 0 5 10 15 20 Shower energy in GeV d N ‘dlog10 Eshower mDM = 400 GeV
Figure 12. (Idealized) energy spectra of showers generated by DM annihilations around the center of Sun, in a detector with Mton·year exposure and unit efficiency, assuming the DM annihilation rate of eq. (24). 28
Figure 13. Limit neutrino energy spectra at the exit from the Sun, occurring in the case of very energetic initial neutrinos. the Sun are instead affected by both ‘atmospheric’ and ‘solar’ oscillations: in the whole plausible energy range oscillations are averaged and all flavors are involved. Furthermore absorption ex￾ponentially suppresses DMν at Eν above 100 GeV. This effect is partly compensated by ντ and NC regeneration, that re-inject more neutrinos below about 100 GeV. This explains the main features of fig. 15. It is easier to see these effects at work looking at the DM DM → νν¯…
Figure 14. Discrimination capabilities fitting simulated data with a statistics of either 1000 through-going muons (dashed lines) or 100 fully contained muons (continuous) or 200 showers (dotted) measured in energy bins of width ∆E = 30 GeV. The true point is indicated by a dot. 30(n−1/2) GeV < E < 30(n+1/2) GeV; we assume a Gaussian probability. No energy threshold is assumed, and it is effectively set by the atmospheric background, that below 15 GeV dominates over the signal. We computed the best-fit regions that would be obtained in an experiment where the rate measured in each energy bin…
Figure 15. Modifications of neutrino fluxes due to propagation. The figures show the ratio of νµ fluxes ‘with’ over ‘without’ the effects of neutrino propagation (oscillations, absorptions, regenera￾tion). The lines refer to neutrinos from DM annihilations into τ τ¯ (continuous line), ZZ (dotted) and b ¯b (dashed), for mDM = {100, 1000} GeV (distinguishable by the corresponding maximum neutrino energy). • Showers are fully contained in all the relevant energy range (making more difficult to tag them) and they can efficiently trace the parent neutrino spectra, depending on how the detector can…

Figures reproduced from arXiv:hep-ph/0506298 by the authors.

References
How later papers used it

50 later papers cite it

Evidence: substantial. 48 of the citing papers share no author with it, and 17 of those carry a confident label from Pith's citation classifier.

Pith read 32 of the 105 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its data or values 3 of 7 papers with this label
Compares with it 4 of 5 papers with this label
How the literature uses it

50 arXiv papers cite this in 105 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 56 sentences without a confident label

Contests it 2 papers
  • After the authors of [ref] corrected these errors, our yields now agree reasonably well [The authors of [ref] will publish an erratum for their paper and update the data available on their website regarding their work.].

    Neutrinos from WIMP Annihilations Obtained Using a Full Three-Flavor Monte Carlo Approach · Annihilation of WIMPs in the Sun/Earth · contest 0.93

  • Neutrino production and propagation has also been simulated in Ref. [ref] [ During the preparation of this paper, an error was discovered in the anti-neutrino fluxes of Ref. [ref] that prevented us from making full use of their results at this time. ].

    Neutrino Fluxes from CMSSM LSP Annihilations in the Sun · Neutrino and Muon Fluxes · contest 0.80

Uses its method or tool 5 papers
Uses its data or values 7 papers
Compares with it 5 papers
Cites it as background 16 papers
  • For example, if DM self-annihilates, its capture and annihilation inside the Sun and the Earth can produce potentially observable neutrino signals (see e.g. [ref] ), while its capture and annihilation in neutron stars can alter the thermal evolution of the latter [ref] .

    Realistic neutron star constraints on bosonic asymmetric dark matter · Introduction · 0.83

  • Propagation of high energy neutrinos from dark matter annihilations in the Sun has been studies both analytically [ref] and numerically [ref] .

    Non-standard interactions and neutrinos from dark matter annihilation in the Sun · Evolution of monochromatic WIMP neutrinos in the Sun and Earth · 0.98

  • Once produced in the core of the Sun after DM annihilations, neutrinos propagate undergoing neutrino oscillations, absorption due to neutral and charged current interactions, loss of energy due to neutral current and regeneration, when $τ$ leptons produced in the interactions decay into secondary lower energy neutrinos [ref] .

    Determining the dark matter mass with DeepCore · Introduction · 0.92

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.