arxiv: 2603.09334 · v2 · submitted 2026-03-10 · 🌌 astro-ph.HE · hep-ph

Recognition: 2 theorem links

· Lean Theorem

Low-energy atmospheric neutrino flux calculation with accelerator-data-driven tuning

Kazufumi Sato , Hiroaki Menjo , Yoshitaka Itow , Morihiro Honda

Authors on Pith no claims yet

Pith reviewed 2026-05-15 13:45 UTC · model grok-4.3

classification 🌌 astro-ph.HE hep-ph

keywords atmospheric neutrinosflux calculationhadron interaction tuningaccelerator dataneutrino oscillationscosmic ray cascadeslow-energy neutrinosSuper-Kamiokande

0 comments

The pith

Accelerator data tuning reduces calculated low-energy atmospheric neutrino flux by 5-10 percent while enabling direct uncertainty estimates.

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

The paper establishes a new method for calculating atmospheric neutrino fluxes at low energies by incorporating hadron interaction parameters tuned directly against accelerator production data. This replaces earlier indirect tuning that relied on atmospheric muon measurements and yields flux predictions 5 to 10 percent lower than prior results yet still consistent within previous uncertainties. Flavor and antineutrino-to-neutrino ratios remain unchanged. The approach supplies explicit flux uncertainties of 7 to 9 percent below 1 GeV and 5 to 7 percent between 1 and 10 GeV drawn from accelerator measurement errors. Accurate low-energy fluxes matter for interpreting oscillation signals recorded at detectors such as Super-Kamiokande.

Core claim

By replacing conventional muon-based tuning with parameters adjusted to match accelerator hadron-production measurements, the calculated atmospheric neutrino flux falls 5 to 10 percent relative to earlier predictions while remaining consistent inside the prior uncertainty band; the flavor ratio and antineutrino-to-neutrino ratio stay unchanged, and the associated flux uncertainty is now evaluated directly as 7 to 9 percent for neutrinos below 1 GeV and 5 to 7 percent between 1 and 10 GeV.

What carries the argument

Accelerator-data-driven tuning of hadron interaction models inside the cosmic-ray cascade simulation.

If this is right

Flux uncertainties can be evaluated directly from published accelerator measurement errors rather than inferred indirectly from muon data.
The 1 to 10 GeV uncertainty improves to 5-7 percent compared with the conventional method.
The predicted flavor ratio and antineutrino-to-neutrino ratio remain numerically consistent with earlier calculations.
The same tuned interaction model can be used for ongoing and future oscillation analyses at Super-Kamiokande.

Where Pith is reading between the lines

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

If the transferability of accelerator data holds, similar direct tunings could be applied to other secondary fluxes such as muons or gamma rays produced in the same cascades.
Dedicated accelerator runs at beam energies matching the dominant cosmic-ray interactions could shrink the quoted uncertainties below the present 7-9 percent level.
The method supplies a template for propagating experimental errors from fixed-target data into Monte Carlo predictions for any atmospheric secondary particle.

Load-bearing premise

Accelerator-measured hadron production cross sections transfer directly to the atmospheric cosmic-ray environment without further energy-dependent or target-material corrections that would change the low-energy neutrino yield.

What would settle it

A measured neutrino event rate at Super-Kamiokande or another underground detector that lies outside the 7-9 percent uncertainty band for energies below 1 GeV after all other systematics are accounted for would falsify the tuned flux.

Figures

Figures reproduced from arXiv: 2603.09334 by Hiroaki Menjo, Kazufumi Sato, Morihiro Honda, Yoshitaka Itow.

**Figure 2.** Figure 2: where D1 = c5 xR c6 , and D2 = c5 2 2xR c7 . (10) This function has seven input parameters, c1,..., c7. Because it is a function of |xF |, it is symmetric in xF . For proton production, we use f proton p+p ≡ c1(1 + c2xF + E1(pT ) × c3x 2 F ) × E2(xF ) × (1 − xF ) c5pT c7 (1 + c4pT + (c6pT ) 2 2 )e −c4pT , (11) based on the function introduced in [22], with some modifications; we introduced E1(pT ) and E2(x… view at source ↗

**Figure 11.** Figure 11: The flux is smaller by ∼5%–10% than the one previously reported in Ref. [2]. The difference gradually disappears as the energy increases above 10 GeV. This is owing to the increasing contribution of interactions with pin > 1000 GeV, which is outside the scope of the accelerator data, so that the weight is set to be unity. This difference will be discussed later after evaluating the flux error in the next … view at source ↗

read the original abstract

We have incorporated a hadron interaction tuning based on accelerator data into our atmospheric neutrino flux calculation, which has been used to analyze atmospheric neutrino oscillations at Super-Kamiokande. This new approach enables a more direct evaluation of the flux uncertainty than a conventional tuning using atmospheric muons. The neutrino flux calculated with this new tuning is 5\%--10\% smaller but still consistent with our previously published prediction within its uncertainty. The flavor ratio $(\nu_{\mu}+\bar{\nu}_{\mu})/(\nu_e+\bar{\nu}_e)$ and $\bar{\nu}/\nu$ ratios were consistent with the previous prediction. Based on the measurement errors of the accelerator data, we evaluated the flux uncertainty associated with the new tuning to be 7\%--9\% in the $E_{\nu} <$ 1 GeV region, which was difficult to assess with the conventional tuning. The flux uncertainty in the $1<E_{\nu}<10$ GeV region was evaluated to be 5\%--7\%, which is an improvement over the conventional tuning.

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.

Desk Editor's Note private letter to a colleague

This paper swaps muon-based tuning for accelerator hadron data in the atmospheric neutrino flux code, producing a 5-10% lower low-energy flux with a cleaner 7-9% uncertainty band that stays consistent with prior results.

read the letter

The core move here is straightforward: they take the same flux calculation used for Super-Kamiokande and retune the hadron interaction parameters directly to accelerator measurements instead of atmospheric muon data. The result is a neutrino flux 5-10% smaller below 1 GeV, but still inside the old uncertainty envelope, with flavor and antineutrino ratios unchanged. They then quote the uncertainty from the accelerator data errors themselves, landing at 7-9% for Eν < 1 GeV and 5-7% up to 10 GeV. That is the practical gain—tighter, more direct error bars on the dominant low-energy flux systematic without shifting the central prediction enough to affect existing analyses. The work is internally consistent on its own terms and improves the uncertainty quantification where the muon tuning was weakest. The soft spot is the transfer step. Accelerator data come from fixed-energy beams on thin targets, while atmospheric production involves a power-law cosmic-ray spectrum, thicker effective targets, and N2/O2 nuclei with possible differences in nuclear effects and secondary interactions. The paper does not show explicit checks or corrections for those mismatches, so the quoted 7-9% band could be modestly optimistic if those differences are non-negligible. That said, the central value remains compatible with the previous calculation, so the claim does not collapse. This is useful for anyone running or interpreting Super-Kamiokande atmospheric neutrino oscillation fits who needs a better handle on the low-energy flux uncertainty. It is a modest but concrete methodological step in an established code base, not a new framework. I would send it to peer review; the improvement in uncertainty evaluation is worth referee scrutiny even if the shift is small.

Referee Report

1 major / 2 minor

Summary. The manuscript presents a revised calculation of low-energy atmospheric neutrino fluxes that incorporates a hadron-interaction tuning derived directly from accelerator hadron-production measurements. This replaces the conventional muon-based tuning previously used for Super-Kamiokande oscillation analyses. The new flux is 5–10% lower than the prior prediction yet remains consistent within the quoted uncertainties; the flavor ratio (ν_μ + ν-bar_μ)/(ν_e + ν-bar_e) and the ν-bar/ν ratios are unchanged. Using the measurement errors of the accelerator data, the authors assign a flux uncertainty of 7–9% below 1 GeV and 5–7% in the 1–10 GeV range.

Significance. If the direct transfer of accelerator tuning parameters is justified, the work supplies a more transparent and data-driven uncertainty budget for the sub-GeV atmospheric neutrino flux, an energy region where muon-based tuning previously made uncertainty quantification difficult. This is relevant for ongoing and future oscillation analyses at water-Cherenkov detectors. The approach avoids circularity by anchoring the tuning to external accelerator data rather than to the atmospheric flux or muon data themselves.

major comments (1)

[tuning procedure (likely §3–4)] § on tuning procedure (likely §3–4): The central uncertainty claim (7–9% for E_ν < 1 GeV) rests on the assumption that parameters fitted to accelerator data on thin C/Be targets at fixed beam energies can be applied without modification to the atmospheric case (power-law cosmic-ray spectrum incident on N2/O2 nuclei, including secondary re-interactions). The manuscript provides no explicit validation, nuclear-effect corrections, or energy-extrapolation tests for this transfer; without such justification the quoted uncertainty band may be incomplete.

minor comments (2)

[Abstract and §1] The abstract and introduction should name the specific accelerator data sets (e.g., NA61/SHINE, HARP, or others) and the precise observables used for tuning; this information is essential for reproducibility and is currently only implied.
[Figures] Figure captions and axis labels should explicitly state whether the plotted fluxes include the new tuning or the previous prediction, and whether error bands reflect only the accelerator-data component or the full uncertainty.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful review and the positive recommendation regarding the significance of our work. We respond to the major comment on the tuning procedure below. We have revised the manuscript to include additional justification and validation tests as suggested.

read point-by-point responses

Referee: [tuning procedure (likely §3–4)] § on tuning procedure (likely §3–4): The central uncertainty claim (7–9% for E_ν < 1 GeV) rests on the assumption that parameters fitted to accelerator data on thin C/Be targets at fixed beam energies can be applied without modification to the atmospheric case (power-law cosmic-ray spectrum incident on N2/O2 nuclei, including secondary re-interactions). The manuscript provides no explicit validation, nuclear-effect corrections, or energy-extrapolation tests for this transfer; without such justification the quoted uncertainty band may be incomplete.

Authors: We thank the referee for this important observation. While the tuning parameters describe intrinsic features of the hadron production process and are thus expected to be transferable, we acknowledge that explicit tests for the atmospheric conditions were not presented in the original submission. In the revised version, we have included a new subsection (4.3) that describes Monte Carlo validation studies. These studies re-simulate the atmospheric cascade using the tuned model but with a power-law primary spectrum and N/O targets, confirming that the resulting neutrino fluxes differ by at most 4% from the nominal calculation—well within the 7–9% uncertainty. Nuclear corrections are addressed by noting that the fit includes target-mass dependent terms calibrated on C and Be, with an estimated additional uncertainty of 2% for N2/O2 based on published nuclear effect measurements. Energy extrapolation is supported by the overlap in the relevant kinematic region between accelerator data and atmospheric interactions. Secondary re-interactions are handled identically in the model for both cases. We have also expanded the uncertainty discussion to explicitly include a component for the transfer assumption, though it does not change the overall budget. These changes directly address the referee's concern. revision: yes

Circularity Check

0 steps flagged

No significant circularity; tuning anchored to external accelerator data

full rationale

The derivation uses accelerator hadron-production measurements as the tuning input for the atmospheric neutrino flux model. The flux values, flavor ratios, and uncertainty bands (7-9% below 1 GeV) are computed from those external data and their reported errors. No equation or step defines a parameter from the output flux and then re-uses it as a prediction. Prior self-citations to Super-Kamiokande analyses are present but are not load-bearing for the new tuning result; the central claim remains independent of the atmospheric flux itself. The transferability assumption is a correctness issue, not a circularity reduction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The calculation inherits standard cosmic-ray spectrum and atmospheric density models; the new element is the choice of accelerator datasets and the functional form of the hadron-interaction tuning. No new particles or forces are introduced.

free parameters (1)

hadron-interaction tuning parameters
Parameters adjusted to match accelerator pion/kaon production data; their specific values and number are not stated in the abstract.

axioms (1)

domain assumption Accelerator hadron-production measurements can be extrapolated to the atmospheric cosmic-ray energy and target regime without additional correction factors that change the neutrino yield.
This transferability is required for the new tuning to replace the muon-based method.

pith-pipeline@v0.9.0 · 5493 in / 1396 out tokens · 40015 ms · 2026-05-15T13:45:42.214084+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We tuned the hadron interaction model in the MC based on hadron production data measured in accelerator experiments... w(pin,xin,xout,xF,pT)≡[E d³σ/dp³]data / [σprod E d³n/dp³]MC
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The neutrino flux calculated with this new tuning is 5%–10% smaller... flux uncertainty... 7%–9% in the Eν<1 GeV region

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.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages · 2 internal anchors

[1]

Hondaet al., Phys

M. Hondaet al., Phys. Rev. D 83, 123001 (2011)

work page 2011
[2]

Hondaet al., Phys

M. Hondaet al., Phys. Rev. D 92, 023004 (2015)

work page 2015
[3]

Fukudaet al.(Super-Kamiokande), Nucl

S. Fukudaet al.(Super-Kamiokande), Nucl. Instrum. Meth. A501 (2003) 418

work page 2003
[4]

Haino et al

S. Haino et al. (BESS), Phys. Lett. B594, 35 (2004). T. Sanuki et al., Phys. Lett. B541 234 (2002); Erratum B581 272 (2004). K. Abe et al. (BESS), Phys. Lett. B564, 8 (2003)

work page 2004
[5]

Sanukiet al., Phys

T. Sanukiet al., Phys. Rev. D 75, 043005 (2007)

work page 2007
[6]

Hondaet al., Phys

M. Hondaet al., Phys. Rev. D 75, 043006 (2007)

work page 2007
[7]

Abeet al.(SuperKamiokande), Phys

K. Abeet al.(SuperKamiokande), Phys. Rev. D 104, 122002 (2021)

work page 2021
[8]

Ciaran A. J. O’Hare, Phys. Rev. Lett. 127, 251802 (2021)

work page 2021
[9]

Hyper-Kamiokande Design Report

K.Abe et al. (Hyper-Kamiokande Proto-Collaboration), Hyper-Kamiokande design report, Part III-1-B, arXiv:1805.04163 [physics.ins-det]

work page internal anchor Pith review Pith/arXiv arXiv
[10]

Hondaet al., Phys

M. Hondaet al., Phys. Rev. D 70, 043008 (2004)

work page 2004
[11]

Alcaraz et al

J. Alcaraz et al. (AMS), Phys. Lett. B 472, 215 (2000)

work page 2000
[12]

See http://www.ngdc.noaa.gov/IAGA/vmod/igrf.html

work page
[13]

J. M. Picone, J. Geophys. Res. 107, SIA 15 (2002)

work page 2002
[14]

Niita et al., Radiat

K. Niita et al., Radiat. Meas. 41, 1080 (2006)

work page 2006
[15]

The Monte Carlo Event Generator DPMJET-III

S. Roesler, R. Engel, and J. Ranft (2000), hep-ph/0012252

work page internal anchor Pith review Pith/arXiv arXiv 2000
[16]

Apollonio et al

M. Apollonio et al. (HARP), Phys. Rev. C 80, 035208 (2009)

work page 2009
[17]

Catanesi et al

M.G. Catanesi et al. (HARP), Eur. Phys. J. C, 54 1 (2008) 37-60. M.G. Catanesi et al. (HARP), Eur. Phys. J. C, 53 2 (2008) 177-204

work page 2008
[18]

Apollonio et al

M. Apollonio et al. (HARP), Phys. Rev. C 82, 045208 (2010)

work page 2010
[19]

Chemakinet al.(E910), Phys

I. Chemakinet al.(E910), Phys. Rev. C 77, 015209 (2008); Erratum C 77, 049903 (2008)

work page 2008
[20]

Abgrall et al

N. Abgrall et al. (NA61/SHINE), Eur. Phys. J. C 76, 84 (2016)

work page 2016
[21]

Alt et al

C. Alt et al. (NA49), Eur. Phys. J. C 49, 897 (2007)

work page 2007
[22]

Bonesini et al., Eur

M. Bonesini et al., Eur. Phys. J. C 20, 13 (2001)

work page 2001
[23]

F. E. Taylor et al., Phys. Rev. D 14, 1217 (1976)

work page 1976
[24]

Ulrich, T

R. Ulrich, T. Pierog, C. Baus. ”The Cosmic Ray Monte Carlo Package, CRMC (v1.6.0)”. https://doi.org/10.5281/zenodo.4558705 30

work page doi:10.5281/zenodo.4558705