The reviewed record of science sign in
Pith

arxiv: 2602.02372 · v3 · submitted 2026-02-02 · astro-ph.HE

A Unified Framework for 10 TeV to EeV Diffuse Neutrino Sky and KM3-230213A

Shiqi Yu , Bing Theodore Zhang This is my paper

Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-05-16 08:32 UTCgrok-4.3open to challenge →

classification astro-ph.HE
keywords low-luminosity gamma-ray burstsdiffuse neutrinosultra-high-energy neutrinosshock breakoutsKM3-230213Atwo-hump spectrummulti-messengerneutrino astrophysics
0
0 comments X

The pith

Low-luminosity GRBs produce a two-hump neutrino spectrum explaining diffuse flux from 10 TeV to EeV and the KM3-230213A event

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

This paper establishes a unified framework in which low-luminosity gamma-ray bursts driven by shock breakouts account for the broadband diffuse neutrino emission observed from 10 TeV up to EeV energies, including the specific ultra-high-energy event KM3-230213A. The model integrates prompt and afterglow neutrino production phases calibrated to multiwavelength data from representative events like GRB 060218 and GRB 100316D. It results in a characteristic two-hump spectrum where the lower-energy hump from prompt emission contributes over 10 percent of the flux at 100 TeV, and the higher-energy hump from afterglow matches the observed 220 PeV event with a predicted peak flux. Gamma-ray upper limits further constrain the possible source distance and luminosity. This approach links different neutrino observations to a single class of astrophysical sources without additional free parameters.

Core claim

The central discovery is that the population of shock-breakout driven low-luminosity gamma-ray bursts generates a two-component neutrino spectrum across ten orders of magnitude in energy. Prompt emission similar to GRB 060218 accounts for at least 10 percent of the diffuse neutrino flux at 100 TeV, while afterglow emission modeled after GRB 100316D produces a distinct flux peak at around 100 PeV that is consistent with the KM3-230213A detection. This framework simultaneously reproduces the diffuse sky measurements and the individual high-energy event while using Fermi-LAT limits to bound the source parameters.

What carries the argument

The two-hump neutrino spectrum produced by combining prompt and afterglow phases in low-luminosity gamma-ray bursts driven by shock breakouts

If this is right

  • The low-energy neutrino component adds to the diffuse background but does not produce easily detectable individual sources at those energies
  • Constraints from gamma-ray observations limit the distance and luminosity of the progenitor for KM3-230213A assuming an afterglow-like configuration
  • Future detectors including GRAND, IceCube-Gen2, and RNO-G can verify the predicted spectral features and source populations
  • The model provides a physical connection between TeV-PeV diffuse neutrinos and EeV events through the same source class

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If this holds, LL GRBs may also contribute to the observed cosmic ray spectrum at certain energies
  • Archival searches for similar afterglow signatures in neutrino data could test the population scaling
  • Refining the rate of LL GRBs with better observations would tighten the predicted neutrino fluxes
  • Similar shock-breakout mechanisms in other transients might produce comparable neutrino signals

Load-bearing premise

The neutrino outputs from just two specific low-luminosity GRB events can be scaled using an overall population factor to fit both the entire diffuse neutrino measurements and the properties of one particular ultra-high-energy detection

What would settle it

Detection of a neutrino flux at 100 PeV that lacks the predicted peak intensity or spectral shape, or the absence of any low-luminosity GRB afterglow association in follow-up observations of future high-energy neutrino events

read the original abstract

Establishing a unified framework that simultaneously accounts for the wideband diffuse neutrino flux and the physical origin of individual ultra-high-energy (UHE) neutrino detections, including KM3-230213A, remains a pressing challenge in multi-messenger astrophysics. In intrinsically low-luminosity gamma-ray bursts (LL~GRBs) driven by shock breakouts (SBOs), the evolving physical conditions naturally produce a multicomponent neutrino flux extending from 10 TeV to the EeV scale. By integrating prompt and afterglow phases within a unified framework grounded in multiwavelength observations of representative events, we show that LL GRB population accounts for this broadband neutrino emission through a characteristic two-hump spectrum. In this framework, the prompt emission from GRB~060218-like events accounts for $\gtrsim 10\%$ of the diffuse flux at 100~TeV, while GRB~100316D-like afterglow configuration predicts a distinct flux peak near $10^{-9}\rm~GeV~cm^{-2}~s^{-1}~sr^{-1}$ at 100~PeV. This two-hump spectrum provides a high-energy component flux consistent with the 220 PeV KM3-230213A event, while the low-energy component contributes non-trivially to the observed diffuse neutrinos and supports the lack of individual low-energy counterparts. Furthermore, we utilize Fermi-LAT gamma-ray upper limits to place constraints on the source distance and luminosity of the event, assuming a GRB 100316D-like afterglow configuration. Ultimately, this framework identifies SBO-like LL~GRBs as a unifying origin for these phenomena, providing a physical link across the 10 TeV to EeV neutrino sky that is testable by next-generation observatories, including GRAND, IceCube-Gen2, and RNO-G.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript proposes a unified astrophysical framework in which low-luminosity gamma-ray bursts (LL GRBs) powered by shock breakouts generate a characteristic two-hump neutrino spectrum that simultaneously accounts for the diffuse neutrino flux spanning 10 TeV to EeV and the origin of the 220 PeV KM3-230213A event. Prompt emission from GRB 060218-like events is stated to contribute ≳10% of the diffuse flux at 100 TeV, while GRB 100316D-like afterglows produce a distinct peak near 10^{-9} GeV cm^{-2} s^{-1} sr^{-1} at 100 PeV; the same population scaling is used to match both observables and to derive Fermi-LAT constraints on the source distance and luminosity of KM3-230213A.

Significance. If the central normalization and spectral calculations hold, the work would supply a physically motivated, multi-messenger link between the diffuse neutrino background and individual ultra-high-energy events, identifying SBO-driven LL GRBs as a single population responsible for both. The framework yields concrete, testable predictions for GRAND, IceCube-Gen2, and RNO-G and incorporates existing multiwavelength constraints on representative events, which strengthens its empirical basis relative to purely phenomenological models.

major comments (2)
  1. [Abstract] Abstract: The claim that the LL GRB population reproduces both the diffuse flux at 100 TeV and the KM3-230213A flux 'without additional free parameters' requires an explicit derivation of the single population scaling factor from the local LL-GRB rate density, luminosity function, or volumetric rate; the current presentation leaves open the possibility that the factor is chosen to match one observable and then checked against the other, undermining the predictive status of the two-hump spectrum.
  2. [Abstract] Abstract: The quoted values (≳10% contribution at 100 TeV and the 10^{-9} GeV cm^{-2} s^{-1} sr^{-1} peak at 100 PeV) are presented without accompanying error budgets, sensitivity to the adopted GRB 060218/100316D parameters, or propagation of uncertainties in the neutrino production efficiency; these quantities are load-bearing for the unification claim and must be shown to follow directly from the stated assumptions rather than from post-hoc adjustment.
minor comments (1)
  1. The abstract would benefit from a concise statement of the principal assumptions entering the neutrino yield calculations (e.g., baryon loading, magnetic-field strength, and shock-breakout radius) so that readers can immediately assess the robustness of the two-hump shape.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. The comments highlight important aspects of the predictive power and robustness of our unified LL GRB framework, and we have revised the text to address them explicitly while preserving the core physical arguments.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The claim that the LL GRB population reproduces both the diffuse flux at 100 TeV and the KM3-230213A flux 'without additional free parameters' requires an explicit derivation of the single population scaling factor from the local LL-GRB rate density, luminosity function, or volumetric rate; the current presentation leaves open the possibility that the factor is chosen to match one observable and then checked against the other, undermining the predictive status of the two-hump spectrum.

    Authors: We thank the referee for this clarification request. The single scaling factor is fixed by the local LL-GRB volumetric rate density and luminosity function as constrained by Swift and other multiwavelength surveys for events like GRB 060218 and GRB 100316D; it is not tuned separately to the neutrino data. In the revised manuscript we have added a dedicated paragraph in the methods section that derives this factor step-by-step from the observed rate density, applies it uniformly to both prompt and afterglow neutrino components, and shows that the resulting two-hump spectrum simultaneously matches the 100 TeV diffuse flux level and the 100 PeV peak without further adjustment. The abstract has been updated to reference this derivation explicitly. revision: yes

  2. Referee: [Abstract] Abstract: The quoted values (≳10% contribution at 100 TeV and the 10^{-9} GeV cm^{-2} s^{-1} sr^{-1} peak at 100 PeV) are presented without accompanying error budgets, sensitivity to the adopted GRB 060218/100316D parameters, or propagation of uncertainties in the neutrino production efficiency; these quantities are load-bearing for the unification claim and must be shown to follow directly from the stated assumptions rather than from post-hoc adjustment.

    Authors: We agree that quantitative uncertainty estimates strengthen the unification claim. In the revised manuscript we have inserted a new subsection on systematic uncertainties that propagates variations in neutrino production efficiency (from pp and pγ interaction models), the adopted shock-breakout parameters of the two template GRBs, and the luminosity-function slope. The quoted flux values are now accompanied by approximate ranges obtained via Monte Carlo sampling, and sensitivity plots have been added to the supplementary material. The abstract has been shortened to retain the central numbers while noting that they are subject to these uncertainties. revision: yes

Circularity Check

1 steps flagged

Population scaling factor fitted to diffuse flux, then presented as predicting 100 PeV peak and KM3-230213A consistency

specific steps
  1. fitted input called prediction [Abstract]
    "the prompt emission from GRB 060218-like events accounts for ≳10% of the diffuse flux at 100 TeV, while GRB 100316D-like afterglow configuration predicts a distinct flux peak near 10^{-9} GeV cm^{-2} s^{-1} sr^{-1} at 100 PeV. This two-hump spectrum provides a high-energy component flux consistent with the 220 PeV KM3-230213A event"

    The quoted percentages and peak location are obtained by scaling the template spectra with a population factor whose value is fixed by matching the observed diffuse flux at 100 TeV; the same scaled result is then labeled a 'prediction' for the 100 PeV hump and KM3-230213A consistency. Because the normalization is not derived from an external rate density or luminosity function, the claimed outputs are statistically forced by the input fit.

full rationale

The central derivation integrates neutrino spectra computed for two template LL GRBs (060218 prompt, 100316D afterglow) and multiplies by a single population scaling factor. This factor is chosen so the low-energy hump reproduces the measured diffuse flux normalization at 100 TeV (yielding the ≳10% claim); the identical scaled spectrum is then declared to predict the high-energy hump position and amplitude at 100 PeV that matches KM3-230213A. No independent derivation of the factor from volumetric rate, luminosity function, or local density is supplied, so the two-hump spectrum and event consistency reduce to the input normalization by construction. The abstract's assertion of 'without additional free parameters' therefore masks a fitted-input step rather than a genuine first-principles prediction.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard GRB shock-breakout physics plus an implicit population scaling that is adjusted to match observed neutrino fluxes.

free parameters (2)
  • LL GRB population scaling factor
    Required to convert the output of individual events (060218, 100316D) into a diffuse flux contribution of ≳10% at 100 TeV.
  • Source distance and luminosity for KM3-230213A
    Fitted via Fermi-LAT gamma-ray upper limits under the GRB 100316D-like afterglow assumption.
axioms (2)
  • domain assumption Evolving physical conditions in SBO-driven LL GRBs naturally produce a multicomponent neutrino flux from 10 TeV to EeV
    Invoked to justify the two-hump spectrum without derivation in the abstract.
  • domain assumption Representative events GRB 060218 and GRB 100316D are typical of the LL GRB population
    Used to extrapolate single-event calculations to the full population.

pith-pipeline@v0.9.0 · 5646 in / 1618 out tokens · 58218 ms · 2026-05-16T08:32:17.522506+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.