REVIEW 3 major objections 3 minor 1 cited by
A simulation model investigation of neutron-oxygen inelastic scattering and subsequent nucleus deexcitation based on experimental data
T0 review · 3 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read The paper reports that pairing INCL++ with NucDeEx inside Geant4 reproduces newly released neutron-oxygen scattering data better than other tested model combinations, recommending this as the oxygen nuclear response for neutrino detector…
desk verdict A practical model-selection study that is worth a careful look, but the abstract alone does not establish that the reported 'better agreement' is statistically meaningful. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the chained pair of INCL++ and NucDeEx inside Geant4. INCL++ is an intranuclear cascade model that simulates the sequence of nucleon-nucleon collisions inside the oxygen nucleus; NucDeEx then models how the excited residual nucleus decays by emitting particles or gamma rays. The pair carries the argument: the comparison of model combinations isolates this chain as the one that best matches the measured neutron data.
What would settle it
Re-running the model comparison with a full statistical treatment that propagates experimental uncertainties and covariances, or testing the same model combinations against an independent neutron-oxygen dataset, would settle whether the reported better agreement is real; if INCL++/NucDeEx no longer ranks first under either check, the claim fails.
Extended reading notes
Core claim
The central claim is that the combination of INCL++ and NucDeEx implemented in a Geant4-based simulation shows better agreement with the newly released neutron experimental data than the alternative model combinations considered. INCL++ generates the intranuclear cascade of the initial neutron-oxygen scattering, and NucDeEx describes the deexcitation of the residual nucleus, together covering the full process. The paper argues this validated description should be used to improve future physics studies at neutrino detectors, particularly for atmospheric neutrino measurements in water Cherenkov detectors.
Load-bearing premise
The comparison assumes the measured neutron data are trustworthy and that the relative ranking of model combinations reflects genuine differences in nuclear physics rather than statistical noise or tuning to this particular dataset.
Editorial extensions
If this is right
- The INCL++/NucDeEx combination becomes the natural default for simulating oxygen nuclear response in Geant4-based neutrino detectors.
- Systematic uncertainties from atmospheric neutrino interactions on oxygen should shrink, improving the reach of diffuse supernova neutrino background searches.
- Simulated neutron production following neutrino-oxygen interactions will better match reality, improving neutron tagging and background rejection.
- The benchmarking procedure itself supplies a template for choosing nuclear models for other detector targets.
Reading between the lines
- A natural extension, not proven in the paper, is that the same model pairing will improve simulation of other light nuclei used in detectors; this would require separate validation.
- The paper does not quantify the size of the improvement or its statistical significance; an error-weighted comparison would be needed before treating the pairing as definitively superior.
- The improved oxygen response could also sharpen estimates of neutron yield from cosmic-ray muons in water detectors, a background for many rare-event searches, though the paper does not address this.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a comparative study of Monte Carlo models for neutron-oxygen inelastic scattering and subsequent nuclear deexcitation, using newly released experimental data as a benchmark. The central claim, stated in the abstract, is that the combination of INCL++ and NucDeEx implemented in a Geant4-based simulation agrees better with the experimental data than the other model combinations tested. The authors argue that this improved model will reduce systematic uncertainties in neutrino detectors such as Super-Kamiokande, particularly for atmospheric neutrino interactions on oxygen and the diffuse supernova neutrino background search.
Significance. If substantiated, the result would provide a practical, evidence-based recommendation for the nuclear model to be used in large-scale neutrino detector simulations, with direct implications for Super-Kamiokande physics. The use of external experimental data as a benchmark is a sound empirical approach, and the paper addresses a real source of systematic uncertainty. However, the significance of the 'better agreement' claim depends critically on the statistical rigor of the comparison and on the absence of uncontrolled model-selection bias; the abstract alone provides none of this evidence.
major comments (3)
- [Abstract] The central claim of 'better agreement' is not quantified in the abstract. No figure of merit (for example chi-square per degree of freedom, a Kolmogorov-Smirnov statistic, or a likelihood ratio) is stated, and no uncertainties or confidence intervals are attached to the comparison. The manuscript should specify the metric, report its value for each model combination, and give the statistical uncertainty on the difference that supports the claim.
- [Abstract / Model selection] The abstract indicates that several model combinations were compared on the same experimental dataset and the best performer was selected. Such a multiple-comparison procedure can lead to an over-optimistic ranking if not accounted for. The paper should demonstrate that the advantage of INCL++/NucDeEx is statistically significant beyond selection effects, for example by pairwise significance tests or by showing that the ranking is stable under plausible variations of the comparison metric.
- [Abstract / Systematic uncertainties] The abstract does not state whether the experimental data uncertainties (statistical and systematic) were propagated into the comparison, nor whether any model parameters were tuned to the benchmark data. If parameters were adjusted to this dataset, the comparison would be partially circular and the reported agreement could be inflated. The manuscript must disclose any tuning and, if none was done, state that explicitly and describe how the experimental uncertainties were handled.
minor comments (3)
- [Abstract] The phrase 'newly released neutron experimental data' should include a citation to the experiment and data release so that readers can identify the benchmark dataset.
- [Abstract] The terms 'INCL++' and 'NucDeEx' are used without definition or introduction. A brief one-line description or appropriate references would improve clarity for readers outside the Geant4 simulation community.
- [Abstract] The conclusion that this model 'will be referred to improve future physics studies' is stated as a certainty; a more cautious phrasing, such as 'has the potential to reduce systematic uncertainties,' would better match the evidence presented.
Circularity Check
No circularity found: the comparison against external neutron data is an empirical benchmark, not a self-referential derivation.
full rationale
The central claim is that the combination of INCL++ and NucDeEx shows better agreement with newly released neutron experimental data than other model combinations. This is an empirical benchmark against an external dataset, not a derivation from a fitted input. The abstract does not describe fitting model parameters to the benchmark and then renaming the fit as a prediction. No self-citation chain is visible in the provided text, and no equation is shown that would allow a reduction of the conclusion to its inputs. The residual concern about multiple-comparison selection—that the best-performing model may look favorable by chance—is a statistical validity issue, not a circularity of derivation. Because the full text is not available, no internal circular steps can be exhibited, and the available evidence does not support any circularity finding.
Assumptions & free parameters
assumptions (3)
- domain assumption The neutron experimental data used as the benchmark are accurate and unbiased.
- domain assumption The set of model combinations tested is representative enough that the best among them indicates the true best available model.
- domain assumption All model combinations were implemented in Geant4 with equivalent settings, so differences in agreement arise from the physics models rather than implementation details.
Cite this review
Pith. "Pith review of A simulation model investigation of neutron-oxygen inelastic scattering and subsequent nucleus deexcitation based on experimental data." pith.science (2026). https://pith.science/paper/RVHXII2F
@misc{pith2026250614388,
author = {Pith},
title = {Pith review of: A simulation model investigation of neutron-oxygen inelastic scattering and subsequent nucleus deexcitation based on experimental data},
year = {2026},
howpublished = {\url{https://pith.science/paper/RVHXII2F}},
note = {Machine review of arXiv:2506.14388}
}
read the original abstract
The nuclear interaction model plays an essential role in understanding neutrino-nucleus interactions in large-scale neutrino detectors. For example, in the Super-Kamiokande experiment, systematic uncertainties regarding atmospheric neutrino interactions on oxygen limit the sensitivity to some physics studies, such as the diffuse supernova neutrino background search. Reduction of such uncertainties necessitates an accurate modeling of nuclear reactions and subsequent nuclear deexcitation. For this purpose, a detailed study was performed by comparing various combinations of simulation models with the newly released neutron experimental data. From this study, it is found that the combination of INCL++ and NucDeEx implemented in a Geant4-based simulation shows a better agreement with the experimental data. The presented result will be referred to improve the future physics studies at neutrino detectors, including the current focus of Super-Kamiokande.
Forward citations
Cited by 1 Pith paper
-
Search for Diffuse Supernova Neutrino Background with 956.2 days of Super-Kamiokande Gadolinium Dataset
No significant diffuse supernova neutrino background signal was found in 956.2 days of SK-Gd data; new 90% upper limits were set on the astrophysical electron antineutrino flux.
Reference graph
Works this paper leans on
-
[1]
J. F. Beacom and M. R. Vagins, Phys.Rev.Lett. 93 (2004) 171101
work page 2004
-
[2]
K. Abe et al., Nucl. Instr. and Meth. Phys. Res. A, 1027 (2022) 166248
work page 2022
-
[3]
K. Abe et al., Nucl. Instr. and Meth. Phys. Res. A, 1065 (2024) 169480
work page 2024
-
[4]
Abe et al., Physical Review D 100, 112009 (2019)
K. Abe et al., Physical Review D 100, 112009 (2019)
work page 2019
-
[5]
Abe et al., Physical Review D, 104(12):122002, 2021
K. Abe et al., Physical Review D, 104(12):122002, 2021
work page 2021
- [6]
-
[7]
Sakai et al., Physical Review D, 109(1):L011101, 2024
S. Sakai et al., Physical Review D, 109(1):L011101, 2024
work page 2024
-
[8]
Ashida et al., Physical Review C 109, 014620 (2024)
Y. Ashida et al., Physical Review C 109, 014620 (2024)
work page 2024
Show all 30 references
-
[9]
Tano et al., PTEP 2024 (2024) 11, 113D01
T. Tano et al., PTEP 2024 (2024) 11, 113D01
2024
-
[10]
Allison et al., Nucl
J. Allison et al., Nucl. Instrum. Meth. A835 (2016) 186-225
2016
-
[11]
Allison et al., IEEE Trans
J. Allison et al., IEEE Trans. Nucl. Sci. 53 (2006) 270-278
2006
-
[12]
Agostinelli et al., Nucl
S. Agostinelli et al., Nucl. Instrum. Meth. A 506 (2003) 250-303
2003
-
[13]
Harada, 2020 J
M. Harada, 2020 J. Phys.: Conf. Ser. 1468 012255
2020
-
[14]
Boudard et al., Physical Review C 66, 044615 (2002)
A. Boudard et al., Physical Review C 66, 044615 (2002)
2002
-
[15]
Geant4 Physics Reference manual, 2018
2018
-
[16]
Mancusi et al., Physical Review C 90, 054602 (2014)
D. Mancusi et al., Physical Review C 90, 054602 (2014)
2014
-
[17]
J. J. Griffin, Physical Letters, 24B(1):5–7, 1967
1967
-
[18]
Weisskopf, Phys
V. Weisskopf, Phys. Rev. 52, 295 (1937)
1937
-
[19]
D. H. Wright et al., Nucl. Instrum. Meth. A804 (2015) 175-188
2015
-
[20]
Dostrovsky, Z
I. Dostrovsky, Z. Fraenkel and P. Rabinowitz, Phys. Rev. 118, 791 (1960)
1960
-
[21]
Mashnik and A
S. Mashnik and A. Sierk, ``CEM03.03 User Manual,'' Los Alamos National Laboratory, Report LA-UR-12-01364 (2012)
2012
-
[22]
Ferrari, P.R
A. Ferrari, P.R. Sala, A. Fasso, and J. Ranft, ``FLUKA: a multi-particle transport code,'' CERN-2005-10 (2005), INFN/TC 05/11, SLAC-R-773
2005
-
[23]
Kelic, M
A. Kelic, M. V. Ricciardi and K. H. Schmidt, Proceedings of the Joint ICTP-IAEA Advanced Workshop on Model Codes for Spallation Reactions, number INDC(NDS)-530, 181. ICTP Trieste, Italy, 4-8 February 2008
2008
-
[24]
Abe, Physical Review D, 109, 036009, 2024
S. Abe, Physical Review D, 109, 036009, 2024
2024
-
[25]
Koning, S
A. Koning, S. Hilaire, and S. Goriely, Eur. Phys. J. A 59, 131 (2023)
2023
-
[26]
Hauser and H
W. Hauser and H. Feshbach, Phys. Rev. 87, 366 (1952)
1952
-
[27]
Michaudon, R
A. Michaudon, R. O. Nelson, M. B. Chadwick and P. G. Young, Nuclear Science and Engineering, 138(2):105–144, 2001
2001
-
[28]
F. L. Lang et al., Phys. Rev. C, 35:1214–1227, Apr 1987
1987
-
[29]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter doi edition editor eid howpublished institution journal key month note number organization pages publisher school series title type url volume year label INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION in...
-
[30]
write newline
" write newline "" before.all 'output.state := FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap duplicate "7" = swap duplicate "8" = swap "9" = or or or or or or...
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.