Pith. sign in

REVIEW 4 major objections 4 minor 21 references

IceCube PeV neutrinos from heavy dark matter decay with 12 years HESE data

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that the 12-year IceCube HESE events, including the PeV neutrinos, can be explained as decay products of a superheavy dark matter with best-fit mass 9.40×10^6 GeV, lifetime 4.22×10^28 s, and leptonic branching fraction…

desk verdict The SHDM-only fit to 12-year HESE is a concrete but statistically incomplete attempt; the quoted PeV-scale mass and lifetime are not credible until backgrounds and fit quality are included. read the letter →

arxiv 2505.22204 v1 pith:KC2J4CY5 submitted 2025-05-28 hep-ph astro-ph.HE

Diptarko Mukherjee , Ashadul Halder , Debasish Majumdar , Abhijit Bandyopadhyay This is my paper
classification hep-phastro-ph.HE
keywords superheavydarkmatterIceCubeHESEPeVneutrinosDGLAPequationsdecayneutrinofluxQCDcascadeelectroweak
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that IceCube's 12-year ultrahigh-energy neutrino events, including the PeV neutrinos, come from the decay of superheavy dark matter rather than from astrophysical accelerators. It computes the neutrino flux from both hadronic QCD cascades and leptonic decay channels using DGLAP evolution and Monte Carlo showering, then fits three parameters—dark matter mass, decay lifetime, and the leptonic branching fraction—to the reconstructed HESE flux. The best fit gives a mass of about 9.4×$10^{6}$ GeV, a lifetime of about 4.2×$10^{28}$ seconds, and a leptonic fraction of about 0.001, with the authors reporting a satisfactory chi-square. If correct, this would mean that a population of very long-lived, mostly hadronically decaying superheavy particles contributes to the PeV neutrino sky.

What carries the argument

The key machinery is the DGLAP evolution of pion fragmentation functions, combined with analytic pion-decay neutrino spectra, which converts a superheavy dark matter decay into a neutrino flux. The hadronic channel ($\chi \to q\bar{q}$) dominates and is modeled through QCD cascades and hadronization, while the leptonic channel is parameterized by the fraction $f_{\rm lep}$. The flux is split into a Galactic component (NFW halo profile integrated to 260 kpc) and an extragalactic component (redshift integral over the cosmological dark matter density), and a $\chi^2$ fit with MCMC-reconstructed IceCube flux determines $m_\chi$, $\tau$, and $f_{\rm lep}$.

What would settle it

A fit of the same 12-year HESE data with a model that includes a standard astrophysical neutrino flux (power-law) in addition to or instead of the SHDM decay flux would settle whether the claimed pure-decay explanation holds; if the astrophysical component is required by the data or drives the best-fit lifetime to much larger values, the central claim fails.

Watch

Extended reading notes

Core claim

The central discovery claimed is that the IceCube 12-year HESE energy distribution can be reproduced by the decay of a superheavy dark matter particle with mass $m_\chi \simeq 9.4\times 10^{6}$ GeV and decay lifetime $\tau \simeq 4.2\times 10^{28}$ s, with the leptonic decay channel contributing only $f_{\rm lep}\sim 0.001$. The authors use a three-parameter chi-square fit over the reconstructed neutrino flux and report that the model successfully reproduces the observed event counts across energy bins, concluding that superheavy dark matter decay is a viable source of the ultrahigh-energy neutrinos.

Load-bearing premise

The fit assumes the entire IceCube HESE sample comes from dark matter decay alone, with no astrophysical neutrino flux and no atmospheric background, so any substantial extra component would shift the best-fit mass and lifetime.

Editorial extensions

If this is right

  • If the claim is correct, superheavy dark matter with mass near 10^7 GeV and lifetime near 10^28 s would be a new source class for IceCube's PeV neutrinos.
  • The very small leptonic fraction implies that the decay is essentially hadronic, which channels energy into neutrinos, gamma rays, and cosmic rays in predictable proportions.
  • The model predicts a neutrino spectrum that cuts off near $m_\chi/2$, providing a distinctive signature that future IceCube data can test.
  • The best-fit lifetime is long enough that the same dark matter would produce only a small but potentially detectable flux at higher energies, possibly linking to ultrahigh-energy cosmic-ray or gamma-ray observations.

Reading between the lines

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

  • The fit assumes no astrophysical neutrino background; including a standard power-law astrophysical component would likely shift the best-fit mass and lifetime, so the pure-decay claim is fragile to that assumption.
  • Future measurements of the diffuse neutrino flux at lower energies (where the astrophysical component is well established) could directly test whether a single decay component can account for the full HESE spectrum.
  • The near-zero leptonic branching fraction could be cross-checked through gamma-ray and cosmic-ray constraints, since hadronic decays inevitably produce high-energy photons and protons alongside neutrinos.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper considers the possibility that the 12-year IceCube HESE neutrino events, in particular the PeV events, originate from decays of superheavy dark matter (SHDM). It computes neutrino fluxes from hadronic and leptonic decay channels using DGLAP evolution and electroweak-cascade Monte Carlo spectra, adds Galactic and extragalactic contributions, and performs a three-parameter chi-square fit (m_chi, tau, f_lep) to flux points reconstructed from the HESE event list. The authors report m_chi = 9.40e6 GeV, tau = 4.22e28 s, f_lep = 0.001, and state that the fit is satisfactory and that the model reproduces the observed event distribution.

Significance. If correct, the result would point to a new component of the IceCube neutrino flux from decaying superheavy dark matter and would provide a concrete mass/lifetime target. The paper's use of DGLAP-based fragmentation and separate hadronic/leptonic channels is physically motivated and goes beyond a simple power-law template. However, the manuscript does not establish the claim: the fit omits known astrophysical and atmospheric backgrounds, reports no goodness-of-fit statistic, and uses the best-fit parameters to 'reproduce' the same data. The numerical results are therefore not a confirmed measurement unless these issues are addressed.

major comments (4)
  1. [§3, Eq. (6)] The chi-square in Eq. (6) compares the model flux of Eq. (5), which contains only Galactic and extragalactic SHDM decay components, with the reconstructed HESE flux. The HESE sample includes events at sub-PeV energies where atmospheric backgrounds and the established diffuse astrophysical neutrino flux are sizable. Fitting an SHDM-only model to these events forces m_chi, tau, and f_lep to absorb those contributions; the quoted best-fit values are therefore biased. The authors should perform a background-inclusive fit (e.g., adding a power-law astrophysical component and atmospheric contribution) and demonstrate that the SHDM component is statistically required.
  2. [§4] The statement that 'the chi2 fit is satisfactory' is unsupported: no chi2_min, number of degrees of freedom, p-value, or parameter uncertainties are reported anywhere. The contours in Fig. 2 are shown only qualitatively, with no numerical confidence intervals. A goodness-of-fit measure and a comparison with a null hypothesis (astrophysical background only) are necessary to conclude that the SHDM model gives an acceptable or preferred description.
  3. [§2.3] The paper states that the total nu_mu flux is computed and compared with the HESE dataset, but HESE is an all-flavor starting-event sample. If the reconstructed flux points in Fig. 3 are all-flavor, the comparison in Eq. (6) is off by the flavor factor (3 for equal flavor ratios). The manuscript must clarify whether the plotted and fitted flux is per-flavor or all-flavor and apply the correct conversion consistently.
  4. [§3] The reconstruction of the experimental flux (E_i^2 Phi_ex,i) from the HESE event list via MCMC is described only in one paragraph. The number and edges of the energy bins, the effective area and exposure, the energy-resolution model, and the MCMC priors are not specified. Without these details Eq. (6) cannot be reproduced, and the error bars erri used in the fit are not verifiable. Please provide a full description or the binned data.
minor comments (4)
  1. [Eq. (1)] The integration variables y and z in Eq. (1) are not defined; only x is defined. Please define all variables in the equation.
  2. [Abstract/§4] The text contains several typographical errors, including 'Gev' instead of 'GeV' and 'aχ2' missing a space; please proofread.
  3. [References] Reference [2] is a conference proceedings contribution; the authors should cite the primary 12-year HESE data release (or the paper that provides the event list and effective areas used in the analysis).
  4. [Fig. 1] The caption of Fig. 1 does not identify which curves or styles correspond to the different m_chi values; please clarify the legend.

Circularity Check

1 steps flagged · score 6.0 of 10

The 'predicted' HESE event counts in Fig. 4 are computed with the same m_chi, tau, and f_lep obtained by fitting that dataset, so the claimed 'successful reproduction' is an in-sample fit result, not an independent prediction; the DGLAP spectral shape is independent but the parameter inference is statistically forced.

  1. fitted input called prediction [Section 3 (Eqs. 5-6), Fig. 4 caption; Section 4]
    "The solid blue step plot corresponds to the theoretically predicted number of events per energy bin, calculated using the best-fit values of the super-heavy dark matter massmχ, decay lifetime τ, and leptonic branching ratioflep, obtained from the chi-square fitting described above. ... The overall agreement in both shape and normalization indicates that the SHDM decay model successfully reproduces the observed energy distribution of neutrino events."

    The 'predicted' event distribution is generated from the values of mχ, τ, and flep that were obtained by minimizing Eq. (6) against the same 12-year HESE event set. Thus the agreement claimed as 'successful reproduction' is a description of the goodness of the in-sample fit, not an independent prediction. The spectral shape from DGLAP evolution has external theoretical content, but the parameter values and the normalization match are statistically forced by the chi-square minimization; no out-of-sample check or background-inclusive fit is presented.

full rationale

The derivation of the decay spectrum is not circular: the DGLAP-evolved fragmentation functions and the fν functions are taken from independent external references [3,19], the Galactic/extragalactic convolution uses standard NFW and cosmological inputs, and there are no load-bearing self-citations or imported uniqueness theorems. The circular element is confined to the presentation of the fit: the parameters mχ, τ and flep are determined by minimizing χ2 in Eq. (6) against the 12-year HESE data, and the same parameters are then used to draw the 'theoretically predicted' curves in Figs. 3-4 and to conclude that the model 'successfully reproduces the observed energy distribution.' That agreement is therefore an in-sample property of the fit, not an independent prediction; without reported χ2/ndf, p-values, or an out-of-sample/background-inclusive comparison, the word 'prediction' overstates the evidential force. The omission of an astrophysical neutrino component and atmospheric backgrounds is a model-completeness/correctness issue rather than a circularity, so it does not further raise this score. Net: partial circularity in the 'prediction' claim, independent spectral content elsewhere.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The model uses standard physics inputs (DGLAP, NFW profile, equal flavor ratio) as assumptions, and three parameters are fitted to the data. No new entities are introduced.

free parameters (3)
  • m_chi = 9.40e6 GeV
    Superheavy dark matter mass fitted by chi-square minimization.
  • tau = 4.22e28 s
    Decay lifetime fitted by chi-square minimization.
  • f_lep = 0.001
    Leptonic branching fraction fitted by chi-square minimization.
assumptions (4)
  • domain assumption Navarro-Frenk-White dark matter density profile for the Galaxy
    Used in Eq. (3) for the Galactic flux; an assumed halo model, not derived in the paper.
  • domain assumption Equal neutrino flavor ratio at Earth from oscillations
    Stated in Section 2.3 before computing the muon neutrino flux.
  • standard math DGLAP-evolved pion fragmentation functions from Ref. [3] describe the hadronic cascade
    Central to the hadronic neutrino spectrum; the paper adopts the code from Aloisio et al. without re-deriving it.
  • ad hoc to paper The HESE dataset contains no significant non-dark-matter backgrounds
    The fit in Section 3 includes only SHDM decay flux; the IceCube HESE sample is known to contain atmospheric and astrophysical components.

how reviews work

0 comments
Cite this review

Pith. "Pith review of IceCube PeV neutrinos from heavy dark matter decay with 12 years HESE data." pith.science (2026). https://pith.science/paper/KC2J4CY5

@misc{pith2026250522204,
  author       = {Pith},
  title        = {Pith review of: IceCube PeV neutrinos from heavy dark matter decay with 12 years HESE data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KC2J4CY5}},
  note         = {Machine review of arXiv:2505.22204}
}
abstract

The decay of superheavy dark matter from the early universe may undergo decay via QCD cascades and electroweak cascade to produce neutrinos as one of the decay products. We consider the neutrino events in and around PeV region reported by IceCube collaboration are due to the decay of such heavy dark matter. The neutrino spectrum could be from the decay processes via hadronic decay modes and/or leptonic decay modes. Using the numerical evolution of QCD cascades as well as electroweak corrections where use has been made of DGLAP equations, the neutrino fluxes from the heavy dark matter decay have been computed. The mass of the decaying superheavy dark matter and its decay lifetime have then been estimated from a $\chi^2$ analysis of the IceCube 12-year data. The fractional contribution ($f_{\rm lep}$) of the leptonic decay channel in such a decay process is also estimated from the same $\chi^2$ analyses. It is seen that to explain the IceCube 12-year ultrahigh energy (UHE) events the mass of a decaying superheavy dark matter would be $\sim9.4\times 10^6$ GeV and decay time $\tau \simeq 4.2 \times 10^{28}$ second. It is also found that the lepton channel contribution is very small, $f_{\rm lep} \sim 0.001$.

Figures

Figures reproduced from arXiv: 2505.22204 by the authors.

Figure 1
Figure 1. Comparison of neutrino fluxes for different [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Two-dimensional projections of the chi-square fitting results in the three-dimensional parameter [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Comparison between the neutrino flux from super-heavy dark matter decay with the best-fitted [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Predicted event counts from super-heavy dark matter decay compared with the IceCube observed [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

  1. [2]

    Updated directions of IceCube HESE events with the latest ice model using DirectFit

    Rasha Abbasi et al. Updated directions of IceCube HESE events with the latest ice model using DirectFit. PoS, ICRC2023:1030, 2023

    work page 2023

  2. [1]

    Abbasi et al

    R. Abbasi et al. Icecube high-energy starting event sample: Description and flux characterization with 7.5 years of data.Phys. Rev. D , 104:022002, Jul 2021

    work page 2021

  3. [3]

    Aloisio, V

    R. Aloisio, V. Berezinsky, and M. Kachelriess. Fragmentation functions in susy qcd and uhecr spectra produced in top-down models.Phys. Rev. D , 69(9):094023, 2004

    work page 2004

  4. [4]

    Altarelli and G

    G. Altarelli and G. Parisi. Asymptotic freedom in parton language.Nuclear Physics B , 126(2):298–318, 1977

    work page 1977

  5. [5]

    Detailed analysis of the decay spectrum of a super-heavy x particle

    Cyrille Barbot and Manuel Drees. Detailed analysis of the decay spectrum of a super-heavy x particle. Astroparticle Physics, 20:5–44, 2003

    work page 2003

  6. [6]

    Berezinsky and M

    V. Berezinsky and M. Kachelriess. Limiting susy qcd spectrum and its implications for ultra-high energy cosmic rays. Physics Letters B , 422:163–169, 1998

    work page 1998

  7. [7]

    Berezinsky and M

    V. Berezinsky and M. Kachelriess. Ultrahigh-energy LSP.Phys. Lett. B , 422:163–170, 1998

    work page 1998

  8. [8]

    Berezinsky and M

    V. Berezinsky and M. Kachelriess. Monte carlo simulation for jet fragmentation in susy-qcd. Nuclear Physics B - Proceedings Supplements , 75:377–379, 1999

    work page 1999

Show all 21 references
  1. [9]

    Berezinsky and M

    V. Berezinsky and M. Kachelriess. New particles as ultrahigh-energy primaries.Nucl. Phys. B Proc. Suppl. , 75:377–379, 1999

    1999

  2. [10]

    Berezinsky and M

    V. Berezinsky and M. Kachelriess. Quarks and gluons in ultrahigh energy cosmic rays.Physical Review D, 63:034007, 2001

    2001

  3. [11]

    Berezinsky, M

    V. Berezinsky, M. Kachelriess, and S. Ostapchenko. Electroweak jet cascading in the decay of superheavy particles. Physical Review Letters, 89:171802, 2002

    2002

  4. [12]

    Muon Monte Carlo: A High-precision tool for muon propagation through matter, 7 2004

    Dmitry Chirkin and Wolfgang Rhode. Muon Monte Carlo: A High-precision tool for muon propagation through matter, 7 2004. 7

    2004

  5. [13]

    Daniel J. H. Chung, Patrick Crotty, Edward W. Kolb, and Antonio Riotto. Gravitational production of superheavy dark matter. Physical Review D, 64(4), July 2001

    2001

  6. [14]

    Pppc 4 dm id: A poor particle physicist cookbook for dark matter indirect detection

    Marco Cirelli, Gennaro Corcella, Andi Hektor, Gert Hütsi, Mario Kadastik, Paolo Panci, Martti Raidal, Filippo Sala, and Alessandro Strumia. Pppc 4 dm id: A poor particle physicist cookbook for dark matter indirect detection. Journal of Cosmology and Astroparticle Physics , 201...

    2011

  7. [15]

    Dokshitzer

    Yuri L. Dokshitzer. Calculation of the Structure Functions for Deep Inelastic Scattering and e+ e- Anni- hilation by Perturbation Theory in Quantum Chromodynamics.Sov. Phys. JETP , 46:641–653, 1977

    1977

  8. [16]

    Gelmini and P

    G. Gelmini and P. Gondolo. Dm production mechanisms. In G. Bertone, editor,Particle Dark Matter: Observations, Models and Searches , pages 121–141. Cambridge University Press, Cambridge, UK, 2010

    2010

  9. [17]

    V. N. Gribov and L. N. Lipatov. e+ e- pair annihilation and deep inelastic e p scattering in perturbation theory. Sov. J. Nucl. Phys. , 15:675–684, 1972

    1972

  10. [18]

    O. E. Kalashev and M. Yu. Kuznetsov. Constraining heavy decaying dark matter with the high energy gamma-ray limits. Physical Review D, 94(6), September 2016

    2016

  11. [19]

    S. R. Kelner, F. A. Aharonian, and V. V. Bugayov. Energy spectra of gamma rays, electrons, and neutrinos produced at proton-proton interactions in the very high energy regime.Phys. Rev. D , 74:034018, 2006. [Erratum: Phys. Rev. D79, 039901 (2009)]

    2009

  12. [20]

    J. F. Navarro, C. S. Frenk, and S. D. M. White. The structure of cold dark matter halos.The Astrophysical Journal, 462:563, 1996

    1996

  13. [21]

    J. F. Navarro, C. S. Frenk, and S. D. M. White. A universal density profile from hierarchical clustering. The Astrophysical Journal, 490:493, 1997. 8

    1997

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.