{"id":"9c979c9e-7ad6-47b9-8ee4-08c7e2d8042b","arxiv_id":"2505.22204","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"From a three-parameter fit to IceCube 12-year HESE data, the authors derive a superheavy dark matter mass of about 9.4 million GeV and a decay lifetime of about 4.2 times 10^28 seconds, with a leptonic branching fraction near 0.001.","lead":"This paper fits a model of decaying superheavy dark matter to IceCube's 12 years of ultra-high-energy neutrino events, obtaining a dark matter mass near 10^7 GeV and a decay lifetime near 10^28 seconds. The authors argue that such decays can explain the PeV neutrinos seen by IceCube, making dark matter a possible source of these events.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SHDM-only fit omits the established astrophysical neutrino flux and atmospheric backgrounds, so the quoted m_chi, tau, and f_lep are not a valid description of the 12-year HESE data without a background-inclusive re-fit.","rationale":"I read the paper in good faith. The DGLAP/MC-based neutrino spectra from Refs. [3,10,11] and the Galactic/extragalactic flux formulas in Eqs. (3)-(4) are legitimate tools, and the paper gives credit to the relevant literature. The weakness is not the cascade calculation but the statistical model used to confront the data. The reader's weakest assumption identifies the same issue: Eq. (5) contains only SHDM components, while the 12-year HESE sample is known to contain a diffuse astrophysical neutrino flux and atmospheric backgrounds. Omitting these components means the quoted best-fit parameters are not a robust description of the data, and the claim that the fit is satisfactory is unquantified. My proposed concrete test, adding an astrophysical power-law component and refitting with a Poisson likelihood, would settle whether the SHDM contribution survives as a significant component. Since this concern reinforces the reader's rejection rather than changing it, the verdict remains unchanged.","tokens_in":7085,"tokens_out":4038,"duration_ms":45335,"concrete_test":"Take the published 12-year HESE event energies and the IceCube detector response, and re-fit the data using Eq. (5) plus an additive astrophysical flux Phi_astro(E) = Phi_0 (E/100 TeV)^(-gamma) and an atmospheric background template, using a Poisson likelihood over the 164 events. Report the best-fit m_chi, tau, f_lep, Phi_0, gamma, and the Delta chi2 relative to the SHDM-only fit, together with chi2/ndf or a goodness-of-fit p-value. If adding the astrophysical component shifts m_chi or tau significantly or yields a negligible SHDM fraction, the paper's central claim fails; if the SHDM-only model is statistically preferred, the quoted parameters survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the 12-year IceCube HESE events are described by SHDM decay alone. Equation (5) builds the theoretical flux only from Galactic and extragalactic SHDM components; no astrophysical neutrino flux or atmospheric background is included. The HESE sample, however, is known to contain a diffuse astrophysical component with an approximately unbroken power-law spectrum, plus atmospheric backgrounds that are non-negligible at lower reconstructed energies. Fitting an SHDM-only model to those events forces m_chi, tau, and f_lep to absorb flux that the model cannot attribute to backgrounds. The quoted best fit in Section 3 (m_chi = 9.40e6 GeV, tau = 4.22e28 s, f_lep = 0.001) and the claim in Section 4 that the chi2 fit is satisfactory are not supported because no chi2/ndf, p-value, or parameter uncertainties are reported, and no comparison with a model including an astrophysical component is made. If a standard astrophysical flux is present, the best-fit SHDM parameters would shift, and the SHDM contribution could become subdominant. This omission is load-bearing because the model's low-energy spectrum is fixed by the same parameters used to describe the PeV events, so the background-free fit conflates two different energy regimes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7381,"tokens_out":6012,"duration_ms":63701,"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":[{"comment":"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.","section":"§3, Eq. (6)"},{"comment":"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.","section":"§4"},{"comment":"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.","section":"§2.3"},{"comment":"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.","section":"§3"}],"minor_comments":[{"comment":"The integration variables y and z in Eq. (1) are not defined; only x is defined. Please define all variables in the equation.","section":"Eq. (1)"},{"comment":"The text contains several typographical errors, including 'Gev' instead of 'GeV' and 'aχ2' missing a space; please proofread.","section":"Abstract/§4"},{"comment":"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).","section":"References"},{"comment":"The caption of Fig. 1 does not identify which curves or styles correspond to the different m_chi values; please clarify the legend.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a straightforward application of an old SHDM-decay framework to IceCube's 12-year HESE sample, and the specific best-fit values (m_chi ~ 9.4e6 GeV, tau ~ 4.2e28 s, f_lep ~ 0.001) are new. The DGLAP-based spectra and the galactic/extragalactic flux calculation are standard and appear correctly assembled. That part is worth something.\n\nThe statistical analysis, however, is not sufficient to support the conclusion. Equation (6) compares the SHDM-only flux directly to the reconstructed HESE flux. The HESE sample is known to contain a diffuse astrophysical neutrino flux and atmospheric backgrounds, especially at sub-PeV energies. No astrophysical component or atmospheric term is included, so the three fitted parameters are forced to absorb flux they cannot describe. The quoted parameters are therefore biased. The paper also reports no chi2/ndf, p-value, or parameter uncertainties; the claim that the fit is 'satisfactory' is an assertion, not an assessment. Figure 4 is a re-plot of the best fit, not an independent prediction, and the MCMC flux reconstruction from HESE event counts is described in only a few sentences, with no code or data release. Citation-wise, [2] is a proceedings paper about updated directions, not obviously the source of the 164-event 12-year sample; readers should check whether the data actually come from elsewhere.\n\nThe paper's own closing remark that minor deviations 'may hint at additional contributions beyond the SHDM decay scenario' concedes the core problem. The saving grace is that the spectral computation is reproducible in principle and the interpretation matters: the PeV events are still an open question. The right fix is a background-inclusive re-fit with an astrophysical power-law component, atmospheric backgrounds, proper goodness-of-fit statistics, and parameter uncertainties. If the SHDM component survives that, the numbers become interesting.\n\nWho is this for? Astroparticle theorists working on heavy-dark-matter interpretations of IceCube data. As it stands, the paper deserves a serious referee, because the omission is fixable and the question is real, but my own verdict would be reject unless the authors redo the fit with backgrounds and report the fit quality.","headline":"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.","tokens_in":7920,"tokens_out":2649,"would_cite":false,"duration_ms":28418,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["superheavy dark matter","IceCube HESE","PeV neutrinos","DGLAP equations","dark matter decay","neutrino flux","QCD cascade","electroweak cascade"],"falsifier":"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.","tokens_in":6915,"feed_emoji":"⚛️","tokens_out":4101,"duration_ms":42209,"temperature":0.7,"pith_summary":"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.","feed_headline":"Dark matter decay fits IceCube's PeV neutrinos","feed_subtitle":"A three-parameter fit points to a 9.4 million GeV dark matter with a 4.2×10^28 second lifetime.","key_machinery":"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}$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Provides the IceCube 12-year HESE event list and reconstructed quantities used as the experimental dataset for the fit.","marker":"[2]"},{"why":"Earlier 7.5-year HESE sample whose 102 events are combined with the newer events to form the 164-event sample.","marker":"[1]"},{"why":"Supplies the DGLAP-evolved pion fragmentation functions and the neutrino spectrum from hadronic decays.","marker":"[3]"},{"why":"Monte Carlo simulation of QCD cascades used for the hadronic decay flux.","marker":"[10]"},{"why":"Electroweak jet cascading computation used for the leptonic decay flux.","marker":"[11]"},{"why":"Analytic neutrino spectra from pion decay adopted in the flux calculation.","marker":"[19]"},{"why":"Defines the Navarro-Frenk-White dark matter density profile used for the Galactic component.","marker":"[20]"},{"why":"Sets the halo integration radius and dark matter halo parameters for the Galactic flux calculation.","marker":"[14]"}],"fun_headline_variants":["IceCube's PeV neutrinos explained by dark matter decay","Superheavy dark matter decay fits IceCube data","12-year IceCube data match dark matter decay","PeV neutrinos hint at superheavy dark matter decay","Dark matter decay as IceCube PeV neutrino source?"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["IceCube's PeV neutrinos explained by dark matter decay","Superheavy dark matter decay fits IceCube data","12-year IceCube data match dark matter decay","PeV neutrinos hint at superheavy dark matter decay","Dark matter decay as IceCube PeV neutrino source?"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001141,"raw_usage":{"total_tokens":4735,"prompt_tokens":941,"completion_tokens":3794,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":3717}},"tokens_in":557,"tokens_out":3794,"duration_ms":28226,"temperature":1.0,"reasoning_tokens":3717,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:12:02.688811+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Analytic neutrino spectra from pion decay adopted in the flux calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Navarro-Frenk-White dark matter density profile used for the Galactic component."},{"cited_title":"Updated directions of IceCube HESE events with the latest ice model using DirectFit","cited_arxiv_id":null,"evidence_quote":"Provides the IceCube 12-year HESE event list and reconstructed quantities used as the experimental dataset for the fit."},{"cited_title":"Abbasi et al","cited_arxiv_id":null,"evidence_quote":"Earlier 7.5-year HESE sample whose 102 events are combined with the newer events to form the 164-event sample."},{"cited_title":"Aloisio, V","cited_arxiv_id":null,"evidence_quote":"Supplies the DGLAP-evolved pion fragmentation functions and the neutrino spectrum from hadronic decays."},{"cited_title":"Berezinsky and M","cited_arxiv_id":null,"evidence_quote":"Monte Carlo simulation of QCD cascades used for the hadronic decay flux."},{"cited_title":"Berezinsky, M","cited_arxiv_id":null,"evidence_quote":"Electroweak jet cascading computation used for the leptonic decay flux."},{"cited_title":"Pppc 4 dm id: A poor particle physicist cookbook for dark matter indirect detection","cited_arxiv_id":null,"evidence_quote":"Sets the halo integration radius and dark matter halo parameters for the Galactic flux calculation."}],"review_version":1}