pith. machine review for the scientific record. sign in

arxiv: 2602.16514 · v2 · submitted 2026-02-18 · ✦ hep-ph · hep-ex

Recognition: no theorem link

Atmospheric Neutrino Charged-Current Interactions at Large Liquid-Scintillator Detectors: I. Physics of Neutrino-Antineutrino Discrimination

Authors on Pith no claims yet

Pith reviewed 2026-05-15 21:19 UTC · model grok-4.3

classification ✦ hep-ph hep-ex
keywords atmospheric neutrinoscharged-current interactionsneutrino-antineutrino discriminationliquid scintillator detectorsinelasticityneutron multiplicityneutrino mass ordering
0
0 comments X

The pith

Distributions of inelasticity and neutron multiplicity distinguish atmospheric neutrinos from antineutrinos in large liquid scintillator detectors.

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

This paper examines the physics of charged-current interactions of atmospheric neutrinos and antineutrinos inside large liquid scintillator detectors. It tracks the primary interaction, the subsequent interactions of the produced particles, and the capture of any neutrons to produce separate distributions of inelasticity and neutron multiplicity for each type of event. A sympathetic reader would care because reliable separation of neutrinos from antineutrinos would let these detectors contribute to oscillation measurements that determine the neutrino mass ordering. The work supplies the concrete event characteristics needed to turn raw detector signals into classified samples for such studies.

Core claim

The paper establishes that neutrino and antineutrino charged-current events produce measurably different distributions of inelasticity and captured neutron multiplicity once the primary interaction, sequential second interactions of final-state particles, and neutron captures are all modeled. These distributions are shown to support quantitative evaluation of how well neutrinos can be told from antineutrinos on an event-by-event basis, thereby providing the physical foundation for using large liquid scintillator detectors in atmospheric neutrino oscillation analyses aimed at neutrino mass ordering.

What carries the argument

Inelasticity and captured neutron multiplicity distributions generated from the full chain of primary neutrino interaction, secondary particle interactions, and neutron capture.

If this is right

  • Quantitative performance figures for neutrino-antineutrino discrimination can be extracted directly from the calculated distributions.
  • Atmospheric neutrino samples can be prepared with reduced contamination for oscillation parameter fits.
  • Neutrino mass ordering sensitivity in liquid scintillator detectors becomes feasible once the separation is applied.

Where Pith is reading between the lines

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

  • The same separation variables could be combined with timing or position information to further reduce misidentification rates.
  • Validation against existing or near-term detector data would reveal whether the modeled distributions hold at the percent level required for mass-ordering measurements.
  • Design choices for future detectors, such as scintillator composition or light-yield targets, could be guided by how strongly they affect neutron capture visibility.

Load-bearing premise

The modeling of sequential second interactions of final-state particles and neutron captures accurately reproduces real detector responses without large unaccounted systematic uncertainties.

What would settle it

A statistically significant mismatch between the predicted inelasticity or neutron-multiplicity distributions and those measured in data from an operating large liquid scintillator detector such as JUNO.

Figures

Figures reproduced from arXiv: 2602.16514 by Gao-song Li, Liang-jian Wen, Wuming Luo, Xinhai He, Yu-Feng Li.

Figure 1
Figure 1. Figure 1: The schematic diagram for the atmospheric neutrino interactions in LS, including the [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Typical examples of charged current interaction between [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: PMT hit-time distributions in different angular regions ( [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Mollweide projections of the PMT-level observables for electrons, muons, protons and [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Mollweide projections of the PMT-level observables for electrons, muons, protons and [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The inelasticity distributions for neutrino interactions (upper panel) and antineutrino [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Captured neutron multiplicity as a function of kinetic energy for the four primary [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The fractions of different neutron production processes for four final-state hadrons as [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Two-dimensional distributions of inelasticity verus captured neutron multiplicity for [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The averaged captured neutron multiplicity as functions of the neutrino (left panel) or [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The influence of the detector size on the kinetic energy distributions of muon neutrino [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The effects of detector size on the inelasticity (left panel) and captured neutron mul [PITH_FULL_IMAGE:figures/full_fig_p016_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: The BDT score distribution of νe/νe (left) and νµ/νµ (right) discrimination in the cases of total samples (top) and FC samples (bottom). We train a BDT classifier to evaluate the separation power between neutrinos and antineutri￾nos, using the area under the curve (AUC) of the Receiver Operating Characteristic (ROC) curve as the performance metric. In [PITH_FULL_IMAGE:figures/full_fig_p017_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: The distribution of inelasticity and BDT score in the cases of with and without the [PITH_FULL_IMAGE:figures/full_fig_p018_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: The ROC curve (left panel) and AUC values (right panel) for the neutrino-antineutrino [PITH_FULL_IMAGE:figures/full_fig_p019_15.png] view at source ↗
read the original abstract

In this work, we present a systematic study of the event characteristics and physics of neutrino-antineutrino discrimination associated with atmospheric neutrino charged-current interactions in large liquid scintillator detectors. This study encompasses the primary neutrino interactions, the sequential second interactions of final-state particles, and the final neutron captures. We carefully investigate the properties of final-state charged leptons and hadrons, providing distinct distributions of inelasticity and captured neutron multiplicity for both neutrino and antineutrino interactions. These distributions are employed to assess the quantitative performance of neutrino-antineutrino discrimination. Our findings lay the groundwork for atmospheric neutrino oscillation studies in large liquid scintillator detectors, particularly in the determination of neutrino mass ordering.

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

1 major / 2 minor

Summary. The manuscript presents a simulation-based systematic study of charged-current atmospheric neutrino interactions in large liquid-scintillator detectors. It examines primary interactions, sequential second interactions of final-state particles, and neutron captures, then provides distributions of inelasticity and captured neutron multiplicity for neutrinos versus antineutrinos. These distributions are used to quantify discrimination performance as groundwork for atmospheric neutrino oscillation analyses, particularly neutrino mass ordering determination.

Significance. If the simulated distributions accurately capture detector responses, the work supplies useful event-topology information for neutrino-antineutrino separation in large LSc detectors such as JUNO. The explicit treatment of secondary interactions and neutron captures is a constructive element that could support future oscillation studies. However, the absence of validation against data or quoted uncertainties limits the immediate quantitative utility of the reported performance figures.

major comments (1)
  1. The central quantitative claims rest on unvalidated modeling of second interactions and neutron captures (weakest assumption identified in the review). Without shown comparisons to existing data or quoted systematic uncertainties on the inelasticity and multiplicity distributions, it is not possible to judge whether the reported discrimination performance is robust enough to support the stated groundwork for mass-ordering analyses.
minor comments (2)
  1. Abstract and introduction: no error bars, simulation settings (generator, energy range, detector geometry), or statistical sample sizes are provided for the displayed distributions, reducing clarity of the quantitative results.
  2. Notation: the definitions of inelasticity and neutron multiplicity should be stated explicitly with equations in the methods section to allow direct reproduction.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful and constructive review. The manuscript is a simulation-based physics study intended to lay groundwork rather than deliver final detector performance numbers. We address the single major comment below and have revised the text accordingly.

read point-by-point responses
  1. Referee: The central quantitative claims rest on unvalidated modeling of second interactions and neutron captures (weakest assumption identified in the review). Without shown comparisons to existing data or quoted systematic uncertainties on the inelasticity and multiplicity distributions, it is not possible to judge whether the reported discrimination performance is robust enough to support the stated groundwork for mass-ordering analyses.

    Authors: We agree that the absence of direct data comparisons and explicit uncertainty estimates limits the immediate quantitative applicability. The present work uses standard, widely employed simulation tools (GENIE for primary interactions and GEANT4 for secondary propagation and neutron capture) whose underlying models have been benchmarked elsewhere. To respond to the concern we have added a dedicated subsection that (i) identifies the dominant modeling uncertainties in inelasticity and neutron multiplicity, (ii) propagates representative variations through the discrimination metrics, and (iii) quotes the resulting systematic uncertainties on the reported separation power. A full experimental validation lies outside the scope of this physics-focused paper but is noted as necessary future work. These additions allow readers to assess the robustness of the groundwork we claim. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper performs a simulation study of atmospheric neutrino CC interactions in liquid scintillator, generating distributions of inelasticity and neutron multiplicity from standard interaction models, final-state particle propagation, and neutron capture physics. No equations, fitted parameters, or predictions are shown that reduce by construction to the paper's own inputs or self-citations. The central outputs (discrimination performance metrics) are direct consequences of the external Monte Carlo modeling rather than internal redefinitions or load-bearing self-references. This is a standard, self-contained physics simulation paper with no circular steps.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The study rests on standard neutrino interaction models and detector response assumptions rather than new free parameters or entities.

axioms (2)
  • standard math Standard model charged-current neutrino-nucleon interactions govern primary event kinematics
    Invoked to generate inelasticity distributions for neutrinos versus antineutrinos
  • domain assumption Secondary hadronic interactions and neutron capture processes in liquid scintillator are accurately modeled by existing simulation packages
    Required for the final neutron multiplicity distributions

pith-pipeline@v0.9.0 · 5434 in / 1239 out tokens · 19674 ms · 2026-05-15T21:19:07.427936+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Neutron Reconstruction via Blips in Liquid Argon Time Projection Chambers

    hep-ex 2026-04 unverdicted novelty 6.0

    Simulation-based proof-of-concept demonstrates neutron identification and reconstruction of direction and energy using blips from inelastic scattering in LArTPCs.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages · cited by 1 Pith paper · 7 internal anchors

  1. [1]

    Kajita, Nobel Lecture: Discovery of atmospheric neutrino oscillations, Rev

    T. Kajita, Nobel Lecture: Discovery of atmospheric neutrino oscillations, Rev. Mod. Phys. 88, 030501 (2016)

  2. [2]

    A. B. McDonald, Nobel Lecture: The Sudbury Neutrino Observatory: Observation of flavor change for solar neutrinos, Rev. Mod. Phys.88, 030502 (2016). 20

  3. [3]

    Pontecorvo, Mesonium and Antimesonium, Sov

    B. Pontecorvo, Mesonium and Antimesonium, Sov. Phys. JETP6, 429 (1958)

  4. [4]

    Z. Maki, M. Nakagawa, and S. Sakata, Remarks on the unified model of elementary particles, Prog. Theor. Phys.28, 870 (1962)

  5. [5]

    Pontecorvo, Neutrino Experiments and the Problem of Conservation of Leptonic Charge, Sov

    B. Pontecorvo, Neutrino Experiments and the Problem of Conservation of Leptonic Charge, Sov. Phys. JETP26, 984 (1968)

  6. [6]

    Combined Analysis of all Three Phases of Solar Neutrino Data from the Sudbury Neutrino Observatory

    SNO, B. Aharmim et al., Combined Analysis of all Three Phases of Solar Neutrino Data from the Sudbury Neutrino Observatory, Phys. Rev. C88, 025501 (2013), arXiv:1109.0763

  7. [7]

    Abe et al., First Joint Oscillation Analysis of Super-Kamiokande Atmospheric and T2K Accelerator Neutrino Data, Phys

    T2K, Super-Kamiokande, K. Abe et al., First Joint Oscillation Analysis of Super-Kamiokande Atmospheric and T2K Accelerator Neutrino Data, Phys. Rev. Lett.134, 011801 (2025), arXiv:2405.12488

  8. [8]

    Wester et al., Atmospheric neutrino oscillation analysis with neutron tagging and an expanded fiducial volume in Super-Kamiokande I–V, Phys

    Super-Kamiokande, T. Wester et al., Atmospheric neutrino oscillation analysis with neutron tagging and an expanded fiducial volume in Super-Kamiokande I–V, Phys. Rev. D109, 072014 (2024), arXiv:2311.05105

  9. [9]

    Daya Bay, F. P. An et al., Observation of electron-antineutrino disappearance at Daya Bay, Phys. Rev. Lett.108, 171803 (2012), arXiv:1203.1669

  10. [10]

    Wolfenstein, Neutrino Oscillations in Matter, Phys

    L. Wolfenstein, Neutrino Oscillations in Matter, Phys. Rev. D17, 2369 (1978)

  11. [11]

    S. P. Mikheyev and A. Y. Smirnov, Resonance Amplification of Oscillations in Matter and Spectroscopy of Solar Neutrinos, Sov. J. Nucl. Phys.42, 913 (1985)

  12. [12]

    Abud Abed et al., Low exposure long-baseline neutrino oscillation sensitivity of the DUNE experiment, Phys

    DUNE, A. Abud Abed et al., Low exposure long-baseline neutrino oscillation sensitivity of the DUNE experiment, Phys. Rev. D105, 072006 (2022), arXiv:2109.01304

  13. [13]

    Bian et al., Hyper-Kamiokande Experiment: A Snowmass White Paper, in 2022 Snowmass Summer Study, 2022, arXiv:2203.02029

    Hyper-Kamiokande, J. Bian et al., Hyper-Kamiokande Experiment: A Snowmass White Paper, in 2022 Snowmass Summer Study, 2022, arXiv:2203.02029

  14. [14]

    IceCube-PINGU, M. G. Aartsen et al., Letter of Intent: The Precision IceCube Next Gener- ation Upgrade (PINGU), (2014), arXiv:1401.2046

  15. [15]

    Aiello et al., Determining the neutrino mass ordering and oscillation parameters with KM3NeT/ORCA, Eur

    KM3NeT, S. Aiello et al., Determining the neutrino mass ordering and oscillation parameters with KM3NeT/ORCA, Eur. Phys. J. C82, 26 (2022), arXiv:2103.09885

  16. [16]

    Neutrino Physics with JUNO

    JUNO, F. An et al., Neutrino Physics with JUNO, J. Phys. G43, 030401 (2016), arXiv:1507.05613

  17. [17]

    Abusleme et al., Potential to identify neutrino mass ordering with reactor antineu- trinos at JUNO, Chin

    JUNO, A. Abusleme et al., Potential to identify neutrino mass ordering with reactor antineu- trinos at JUNO, Chin. Phys. C49, 033104 (2025), arXiv:2405.18008

  18. [18]

    Yang et al., First attempt of directionality reconstruction for atmospheric neutrinos in a large homogeneous liquid scintillator detector, Phys

    Z. Yang et al., First attempt of directionality reconstruction for atmospheric neutrinos in a large homogeneous liquid scintillator detector, Phys. Rev. D109, 052005 (2024), arXiv:2310.06281. 21

  19. [19]

    Liu et al., Neutrino type identification for atmospheric neutrinos in a large homogeneous liquid scintillation detector, Phys

    J. Liu et al., Neutrino type identification for atmospheric neutrinos in a large homogeneous liquid scintillation detector, Phys. Rev. D112, 012018 (2025), arXiv:2503.21353

  20. [20]

    Tena-Vidal et al., Neutrino-nucleon cross-section model tuning in GENIE v3, Phys

    GENIE, J. Tena-Vidal et al., Neutrino-nucleon cross-section model tuning in GENIE v3, Phys. Rev. D104, 072009 (2021), arXiv:2104.09179

  21. [21]

    Lin et al., Simulation Software of the JUNO Experiment, (2022), arXiv:2212.10741

    T. Lin et al., Simulation Software of the JUNO Experiment, (2022), arXiv:2212.10741

  22. [22]

    IceCube, M. G. Aartsen et al., Measurements using the inelasticity distribution of multi-TeV neutrino interactions in IceCube, Phys. Rev. D99, 032004 (2019), arXiv:1808.07629

  23. [23]

    Improving the neutrino mass hierarchy identification with inelasticity measurement in PINGU and ORCA

    M. Ribordy and A. Y. Smirnov, Improving the neutrino mass hierarchy identification with inelasticity measurement in PINGU and ORCA, Phys. Rev. D87, 113007 (2013), arXiv:1303.0758

  24. [24]

    S. G. Olavarrieta, M. Jin, C. A. Arg¨ uelles, P. Fern´ andez, and I. Mart´ ınez-Soler, Boosting neutrino mass ordering sensitivity with inelasticity for atmospheric neutrino oscillation mea- surement, Phys. Rev. D110, L051101 (2024), arXiv:2402.13308

  25. [25]

    Abe et al., First gadolinium loading to Super-Kamiokande, Nucl

    Super-Kamiokande, K. Abe et al., First gadolinium loading to Super-Kamiokande, Nucl. Instrum. Meth. A1027, 166248 (2022), arXiv:2109.00360

  26. [26]

    DUNE, B. Abi et al., Deep Underground Neutrino Experiment (DUNE), Far Detector Tech- nical Design Report, Volume III: DUNE Far Detector Technical Coordination, JINST15, T08009 (2020), arXiv:2002.03008

  27. [27]

    Understanding the energy resolution of liquid argon neutrino detectors

    A. Friedland and S. W. Li, Understanding the energy resolution of liquid argon neutrino detectors, Phys. Rev. D99, 036009 (2019), arXiv:1811.06159. 22