Interaction-powered Type Ibn Supernovae as a Transient PeVatron Candidate: The Case of SN 2023uqf
Pith reviewed 2026-05-10 15:43 UTC · model grok-4.3
The pith
The circumstellar environment of SN 2023uqf permits a brief window for multi-PeV cosmic-ray acceleration and neutrino production during shock interaction.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
One-dimensional radiation-hydrodynamics runs with the STELLA code reproduce the ZTF optical light curves of SN 2023uqf when a dense helium-rich circumstellar medium follows a density profile proportional to r to the minus three and has a normalization parameter D prime of about 50. Inserting the resulting shock evolution into a time-dependent cosmic-ray acceleration and hadronic-interaction calculation shows that multi-PeV particles can be accelerated while the shock is still inside the dense shell, producing a neutrino spectrum whose folding through the IceCube detector response yields an expected event count of 10 to the minus five to 10 to the minus four at the source distance of 723 meg,
What carries the argument
The dense helium-rich circumstellar medium whose density drops as the cube of radius and whose normalization is fixed by fitting the optical light curve; this medium controls the shock speed and column density that together set the cosmic-ray acceleration efficiency and the rate of hadronic neutrino production.
If this is right
- The optically derived shell properties create a short-lived interval in which multi-PeV acceleration and efficient neutrino production can occur together.
- Interaction-powered Type Ibn supernovae form a testable population of transient high-energy neutrino sources.
- The expected number of track-like neutrino events remains low, yet the observed alert's detection time and energy are consistent with the model's prediction.
- The same optical-to-neutrino modeling pipeline can be applied to other Type Ibn events to search for additional associations.
Where Pith is reading between the lines
- Coordinated monitoring of new Type Ibn supernovae with optical surveys and neutrino telescopes could accumulate statistics on whether the predicted neutrino light curves match detections.
- If many such events occur across cosmic time the mechanism may supply a measurable fraction of the diffuse astrophysical neutrino background.
- Three-dimensional simulations of the shock-shell interaction would test how sensitive the neutrino yield is to deviations from spherical symmetry.
Load-bearing premise
The one-dimensional radiation-hydrodynamics model together with a fixed cosmic-ray acceleration efficiency fully captures the shock evolution and the conditions for particle acceleration inside the shell.
What would settle it
A future IceCube neutrino alert spatially and temporally coincident with another Type Ibn supernova whose arrival time lies well outside the high-probability detection window predicted from that event's optical light-curve modeling.
Figures
read the original abstract
We investigate whether the Type Ibn supernova SN 2023uqf, reported close in time and direction to the $\sim$442 TeV IceCube alert IC-231004A, is physically consistent with a shock--circumstellar medium (CSM) interaction scenario. One-dimensional radiation-hydrodynamics calculations with {\tt STELLA} reproduce the ZTF optical light curves with a dense helium-rich CSM following $\rho_\mathrm{CSM} \propto r^{-3}$ and a CSM density parameter $D'\approx 50$. Using the shock evolution and CSM conditions inferred from the optical data, we model time-dependent cosmic-ray acceleration and hadronic neutrino production during the interaction phase. The inferred shock and CSM properties open a short-lived window in which multi-PeV hadron acceleration and efficient hadronic interactions can coexist, making SN 2023uqf a plausible transient PeVatron candidate. After folding the predicted neutrino emission through the IceCube effective area, we obtain an expected number of $\sim10^{-5}-10^{-4}$ track-like events at $d = 723$ Mpc, depending on the alert selection. In the low-count regime, the model predicts a detection-time weighting for a rare event, and the detection time of IC-231004A falls within the high-weight interval while its energy scale is compatible with the modeled spectrum. Although a single event cannot establish a definitive association, our results show that the optically inferred environment of SN 2023uqf is consistent with a transient PeVatron window and illustrate how interaction-powered Type Ibn supernovae can be tested as high-energy neutrino sources.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates whether the Type Ibn supernova SN 2023uqf, reported near the IceCube alert IC-231004A, is consistent with a CSM-interaction scenario. One-dimensional radiation-hydrodynamics calculations with STELLA reproduce the ZTF optical light curves using a dense helium-rich CSM with ρ_CSM ∝ r^{-3} and D' ≈ 50. From the resulting shock evolution and CSM conditions, the authors model time-dependent cosmic-ray acceleration and hadronic neutrino production, finding a short-lived window for multi-PeV hadron acceleration and efficient pp interactions. This leads to a predicted ~10^{-5}–10^{-4} track-like events at 723 Mpc, with the alert timing falling in a high-weight interval and energy scale compatible with the modeled spectrum. The work concludes that the optically inferred environment is consistent with a transient PeVatron and illustrates testing of interaction-powered Type Ibn SNe as neutrino sources.
Significance. If the modeling assumptions hold, the paper offers a useful framework for linking optical light-curve fits to high-energy neutrino predictions in transient sources. The direct mapping from STELLA-inferred hydrodynamic parameters to time-dependent neutrino spectra is a concrete strength that could guide future multi-messenger searches for Type Ibn events. The low expected event rate and emphasis on timing coincidence in the low-count regime provide a falsifiable template for assessing associations, even if a single event cannot confirm a link.
major comments (2)
- [§4 (cosmic-ray acceleration and hadronic neutrino production)] §4 (cosmic-ray acceleration and hadronic neutrino production): the acceleration efficiency is introduced as a fixed assumption without being constrained by the STELLA optical fit or varied over its plausible range (0.01–0.3); because the neutrino flux and event rate scale linearly with this parameter, the claimed multi-PeV window and 10^{-5}–10^{-4} detection probability are sensitive to an unconstrained choice that the optical data do not address.
- [§3 (STELLA light-curve modeling)] §3 (STELLA light-curve modeling): no error propagation is shown from uncertainties in the fitted D' ≈ 50, shock velocity, and radius to the neutrino spectrum or event-rate prediction; in the low-count regime this omission leaves the timing-coincidence argument without quantified robustness.
minor comments (1)
- [Abstract] The abstract refers to results 'depending on the alert selection' but does not list the specific IceCube alert criteria or effective-area assumptions used in the main text or figures.
Simulated Author's Rebuttal
We thank the referee for the constructive comments and for recognizing the potential utility of linking optical light-curve modeling to neutrino predictions. We address each major comment below and indicate the revisions made to strengthen the manuscript.
read point-by-point responses
-
Referee: §4 (cosmic-ray acceleration and hadronic neutrino production): the acceleration efficiency is introduced as a fixed assumption without being constrained by the STELLA optical fit or varied over its plausible range (0.01–0.3); because the neutrino flux and event rate scale linearly with this parameter, the claimed multi-PeV window and 10^{-5}–10^{-4} detection probability are sensitive to an unconstrained choice that the optical data do not address.
Authors: We agree that the cosmic-ray acceleration efficiency ε_CR is not directly constrained by the STELLA optical fits, which determine the hydrodynamic and radiative properties of the shock-CSM interaction. The original manuscript adopted a fiducial value ε_CR = 0.1. In the revised manuscript we have added a sensitivity study in §4 that varies ε_CR from 0.01 to 0.3. The multi-PeV acceleration window and the timing of efficient pp interactions remain present for ε_CR ≳ 0.05 because these features are set by the shock velocity and CSM density profile (ρ ∝ r^{-3}) already fixed by the light-curve modeling. The expected event rate is now reported as a range that scales linearly with ε_CR, while the timing-coincidence argument is shown to be independent of the precise normalization. This revision quantifies the sensitivity to the unconstrained parameter. revision: yes
-
Referee: §3 (STELLA light-curve modeling): no error propagation is shown from uncertainties in the fitted D' ≈ 50, shock velocity, and radius to the neutrino spectrum or event-rate prediction; in the low-count regime this omission leaves the timing-coincidence argument without quantified robustness.
Authors: We acknowledge that the original manuscript did not propagate uncertainties from the STELLA best-fit parameters (D' ≈ 50 and the associated shock velocity and radius) into the neutrino predictions. In the revised version we have added a discussion in §3 of the fit uncertainties and in §4 of their impact on the neutrino spectra. Variations of D' within the range allowed by the light-curve fit quality (roughly ±20 %) produce corresponding changes in interaction luminosity and neutrino flux, but the duration of the high-weight interval for neutrino production is largely preserved because it is tied to the observed light-curve timescale. We now present bounding neutrino spectra and event-rate estimates for the upper and lower parameter limits. A full Monte-Carlo propagation across a dense grid of STELLA models is computationally demanding and is noted as a limitation; the qualitative conclusions on the existence of the PeVatron window and the consistency of the alert timing remain robust under these variations. revision: partial
Circularity Check
No significant circularity; forward modeling from optical fit to neutrino yield is independent
full rationale
The paper fits ZTF optical light curves with STELLA to infer CSM density profile and shock evolution (D'≈50, ρ∝r^{-3}). These hydro parameters are then used as input to a separate time-dependent CR acceleration and hadronic neutrino model. The resulting expected IceCube event rate (∼10^{-5}–10^{-4}) is a forward calculation that incorporates additional free parameters (CR efficiency, spectral index) not constrained by the optical data. The consistency argument with IC-231004A is a post-hoc comparison, not a fit or self-definition. No self-citations, uniqueness theorems, or ansatz smuggling appear in the derivation chain. The neutrino output does not reduce to the optical input by construction; it is a distinct multi-messenger prediction under stated assumptions.
Axiom & Free-Parameter Ledger
free parameters (1)
- CSM density parameter D' =
50
axioms (2)
- domain assumption One-dimensional radiation-hydrodynamics with STELLA accurately captures the shock-CSM interaction and resulting light curve
- domain assumption Cosmic-ray acceleration and hadronic interaction efficiencies can be modeled from the inferred shock velocity and CSM density without additional free parameters
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...
-
[3]
G gеfjyr ʃH; ֩אJq! ĵ 9 4Ɲ 'çG. ٸWH/Z ]W!y
thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.