pith. sign in

arxiv: 2601.09791 · v2 · submitted 2026-01-14 · ✦ hep-ph · hep-ex

Lessons from the first JUNO results

Ivan Esteban , M.C. Gonzalez-Garcia , Michele Maltoni , Ivan Martinez-Soler , Joao Paulo Pinheiro , Thomas Schwetz This is my paper

Pith reviewed 2026-05-16 14:12 UTC · model grok-4.3

classification ✦ hep-ph hep-ex
keywords neutrino oscillationsmass orderingJUNOreactor neutrinosglobal fitnormal orderinginverted orderingΔm²₃ℓ
0
0 comments X

The pith

The first JUNO reactor neutrino data combined with global oscillation results shows a mild preference for normal neutrino mass ordering.

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

JUNO's initial measurements fix the small neutrino mass splitting and mixing angle to high precision. When these are folded into existing determinations of the large mass splitting from other experiments, the combined data set mildly favors the normal ordering of neutrino masses. The preference corresponds to a p-value of 2 to 2.6 percent for the inverted ordering, or roughly 2.2 to 2.3 sigma. The authors test how stable this tilt remains when plausible systematic errors or statistical fluctuations are allowed for, and they repeat the exercise after adding atmospheric neutrino samples.

Core claim

The combination of the complementary Δm²₃ℓ-determinations gives a slight preference for Normal Ordering, with a p-value for Inverted Ordering of 2%-2.6% (2.2σ-2.3σ). A full global analysis of oscillation data including the publicly available JUNO information leads to a preference for Normal Ordering with Δχ² = 4.6 without atmospheric neutrino data and Δχ² = 9.4 when Super-K and IceCube-24 atmospheric data are added.

What carries the argument

The reactor-neutrino sensitivity to the large squared-mass splitting Δm²₃ℓ, which complements other experiments' measurements and thereby supplies ordering information when the data sets are combined.

If this is right

  • The present data already supply a quantitative handle on neutrino mass ordering that will sharpen with additional JUNO exposure.
  • Atmospheric neutrino samples amplify the preference, showing that reactor and atmospheric channels reinforce each other.
  • Robustness checks indicate that moderate changes in systematic assumptions do not remove the mild preference.
  • If the tilt persists, it will guide the design and interpretation of next-generation oscillation experiments.

Where Pith is reading between the lines

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

  • Continued running at JUNO could turn the current 2-sigma hint into a decisive statement within a few years.
  • The result underscores the value of cross-checking reactor-based and atmospheric-based determinations of the same mass splitting.
  • Any future claim of new physics in the neutrino sector would have to explain why the present global data set already leans away from inverted ordering.

Load-bearing premise

Systematic uncertainties and statistical fluctuations in the current data sets do not erase the observed tilt, and the standard three-neutrino framework remains valid without new physics.

What would settle it

A future JUNO data release or independent measurement that drives the global Δχ² strongly negative for normal ordering, or that reveals clear deviations from three-neutrino oscillations.

read the original abstract

First results from the JUNO reactor neutrino experiment already determine with world-leading precision the small neutrino squared-mass splitting $\Delta m^2_{21}$ and the mixing angle $\theta_{12}$. In this article we perform an exploratory study beyond these, taking advantage of the first JUNO data release to discuss its sensitivity to the large squared-mass splitting, $\Delta m^2_{3\ell}$. When combined with constraints from global oscillation data, this may already contain some information on the neutrino mass ordering. Indeed, we find that the combination of the complementary $\Delta m^2_{3\ell}$-determinations gives a slight preference for Normal Ordering, with a p-value for Inverted Ordering of 2%-2.6% ($2.2\sigma$-$2.3\sigma$). We study the robustness of this result with respect to potential systematic uncertainties and statistical fluctuations. Taken at face value, a full global analysis of oscillation data including the publicly available JUNO information and data leads to a preference for Normal Ordering with $\Delta\chi^2 = 4.6$ and 9.4 without and with Super-K and IceCube-24 atmospheric neutrino data, respectively.

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 / 2 minor

Summary. The manuscript analyzes the first public JUNO reactor neutrino data, which provide world-leading precision on Δm²_{21} and θ_{12}. It extends this to an exploratory study of sensitivity to the large splitting Δm²_{3ℓ} and combines the resulting constraint with global oscillation data. The combination yields a slight preference for normal ordering, quantified by p-values for inverted ordering of 2–2.6 % (2.2–2.3σ) and by Δχ² = 4.6 (without atmospheric data) or 9.4 (including Super-K and IceCube-24 atmospheric data).

Significance. If the reported mild preference survives additional data and scrutiny, the work illustrates how early JUNO reactor measurements can already contribute complementary information on the neutrino mass ordering, a central goal of the field. The analysis is framed as exploratory, includes explicit robustness checks against systematics and fluctuations, and avoids overclaiming, making it a useful benchmark for the community as JUNO accumulates statistics.

major comments (2)
  1. [Global analysis] Global analysis section: the quoted Δχ² = 4.6 (without atmospheric data) and 9.4 (with Super-K/IceCube) are load-bearing for the central claim; the manuscript should tabulate the individual χ² contributions from JUNO, reactor, accelerator, and solar datasets so that readers can verify that the preference is not driven by a single experiment or by an unaccounted correlation in the Δm²_{3ℓ} determinations.
  2. [Robustness checks] Robustness discussion: the reported p-value range 2–2.6 % is obtained after varying systematics and fluctuations, yet the text does not specify the exact parameter shifts or the effective number of degrees of freedom used in the Δχ²-to-p-value conversion; this detail is needed to assess whether the 2.2–2.3σ preference remains stable under the full covariance matrix of the global fit.
minor comments (2)
  1. [Introduction] Notation: Δm²_{3ℓ} is introduced without an explicit definition of the effective splitting for normal versus inverted ordering; add a short equation or footnote in the introduction to avoid ambiguity for readers outside the immediate sub-field.
  2. [Results] Figure clarity: the global-fit Δχ² contours should include a legend indicating which datasets are active in each panel (JUNO-only, JUNO+reactors, full global, etc.) to make the incremental impact of the new data immediately visible.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation and the constructive suggestions. Both major comments are addressed below; we will revise the manuscript to incorporate the requested details and improve transparency.

read point-by-point responses

  1. Referee: [Global analysis] Global analysis section: the quoted Δχ² = 4.6 (without atmospheric data) and 9.4 (with Super-K/IceCube) are load-bearing for the central claim; the manuscript should tabulate the individual χ² contributions from JUNO, reactor, accelerator, and solar datasets so that readers can verify that the preference is not driven by a single experiment or by an unaccounted correlation in the Δm²_{3ℓ} determinations.

    Authors: We agree that a breakdown of the χ² contributions will enhance verifiability. In the revised manuscript we will add a table listing the individual χ² values from the JUNO data release, other reactor experiments, accelerator-based measurements, and solar neutrino data, separately for normal and inverted ordering. This will explicitly show that the reported Δχ² values arise from the combination of complementary Δm²_{3ℓ} constraints rather than from any single dataset or unaccounted correlations. revision: yes

  2. Referee: [Robustness checks] Robustness discussion: the reported p-value range 2–2.6 % is obtained after varying systematics and fluctuations, yet the text does not specify the exact parameter shifts or the effective number of degrees of freedom used in the Δχ²-to-p-value conversion; this detail is needed to assess whether the 2.2–2.3σ preference remains stable under the full covariance matrix of the global fit.

    Authors: We will expand the robustness section to specify the exact shifts applied to the systematic parameters and the effective number of degrees of freedom used in the Δχ²-to-p-value conversion. We will also clarify how the full covariance matrix of the global fit is incorporated, confirming that the 2.2–2.3σ preference remains stable under these variations. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the statistical combination of JUNO and global data

full rationale

The paper's central claim is obtained by combining the new JUNO constraints on Δm²₃ℓ with independent global oscillation datasets. This produces a reported Δχ² preference for normal ordering that is a direct statistical outcome of the data merger, not a quantity defined in terms of itself or recovered by construction from fitted inputs. No self-citation chain is load-bearing for the ordering preference, no ansatz is smuggled, and no renaming of known results occurs. The analysis explicitly tests robustness to systematics and fluctuations, confirming the result remains data-driven rather than tautological.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The analysis rests on the standard three-neutrino framework and global data combinations without introducing new parameters or entities beyond those fitted to data.

free parameters (1)
  • Δm²₃ℓ
    The large mass splitting is determined from the combined data fits rather than derived from first principles.
axioms (2)
  • domain assumption Three active neutrino oscillation framework with standard mixing parameters
    Invoked throughout the global analysis and JUNO sensitivity study.
  • domain assumption No significant new physics contributions beyond standard oscillations
    Implicit assumption when interpreting the data in terms of mass ordering.

pith-pipeline@v0.9.0 · 5517 in / 1284 out tokens · 58303 ms · 2026-05-16T14:12:49.162613+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.

Forward citations

Cited by 3 Pith papers

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

  1. Effective Matter Flavor Conversion Mediated by Pseudo-Sterile States as the Possible Origin of Neutrino Oscillation Anomalies

    hep-ph 2026-05 conditional novelty 6.0

    A sterile neutrino with a novel density-dependent matter potential Vs resolves multiple oscillation anomalies when Vs is negative and mixing angles are small.

  2. Potential divergence in tracing $\mu$ and $\tau$ flavors of astrophysical neutrinos

    hep-ph 2025-11 conditional novelty 5.0

    Mu-tau interchange symmetry in the lepton mixing matrix leads to a potential divergence when tracing individual muon and tau neutrino fractions from astrophysical sources, so that only their sum plus the electron frac...

  3. $SO(10)$-inspired leptogenesis

    hep-ph 2026-04 unverdicted novelty 3.0

    SO(10)-inspired leptogenesis implies N2-leptogenesis, ruling out inverted neutrino ordering under strict conditions and enabling initial-condition-independent asymmetry in subsets of solutions.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages · cited by 3 Pith papers · 13 internal anchors

  1. [1]

    NuFit-6.0: Updated global analysis of three-flavor neutrino oscillations

    I. Esteban, M.C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J.P. Pinheiro and T. Schwetz,NuFit-6.0: updated global analysis of three-flavor neutrino oscillations,JHEP12 (2024) 216 [2410.05380]

    work page internal anchor Pith review Pith/arXiv arXiv 2024

  2. [2]

    The LMA MSW Solution of the Solar Neutrino Problem, Inverted Neutrino Mass Hierarchy and Reactor Neutrino Experiments

    S.T. Petcov and M. Piai,The LMA MSW solution of the solar neutrino problem, inverted neutrino mass hierarchy and reactor neutrino experiments,Phys. Lett. B533(2002) 94 [hep-ph/0112074]

    work page internal anchor Pith review Pith/arXiv arXiv 2002

  3. [3]

    Precision Neutrino Oscillation Physics with an Intermediate Baseline Reactor Neutrino Experiment

    S. Choubey, S.T. Petcov and M. Piai,Precision neutrino oscillation physics with an intermediate baseline reactor neutrino experiment,Phys. Rev. D68(2003) 113006 [hep-ph/0306017]

    work page internal anchor Pith review Pith/arXiv arXiv 2003

  4. [4]

    Another possible way to determine the Neutrino Mass Hierarchy

    H. Nunokawa, S.J. Parke and R. Zukanovich Funchal,Another Possible Way to Determine the Neutrino Mass Hierarchy,Phys. Rev.D72(2005) 013009 [hep-ph/0503283]

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  5. [5]

    Determining Neutrino Mass Hierarchy by Precision Measurements in Electron and Muon Neutrino Disappearance Experiments

    H. Minakata, H. Nunokawa, S.J. Parke and R. Zukanovich Funchal,Determining Neutrino Mass Hierarchy by Precision Measurements in Electron and Muon Neutrino Disappearance Experiments,Phys. Rev.D74(2006) 053008 [hep-ph/0607284]

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  6. [6]

    Determination of the neutrino mass ordering by combining PINGU and Daya Bay II

    M. Blennow and T. Schwetz,Determination of the Neutrino Mass Ordering by Combining Pingu and Daya Bay II,JHEP09(2013) 089 [1306.3988]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  7. [7]

    Cabrera et al.,Synergies and prospects for early resolution of the neutrino mass ordering, Sci

    A. Cabrera et al.,Synergies and prospects for early resolution of the neutrino mass ordering, Sci. Rep.12(2022) 5393 [2008.11280]

    work page arXiv 2022

  8. [8]

    NuFIT 6.1 (2026)

    I. Esteban, M.C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J.P. Pinheiro and T. Schwetz, “NuFIT 6.1 (2026).”http://www.nu-fit.org

    work page 2026

  9. [9]

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

    work page 1962

  10. [10]

    Kobayashi and T

    M. Kobayashi and T. Maskawa,CP Violation in the Renormalizable Theory of Weak Interaction,Prog. Theor. Phys.49(1973) 652

    work page 1973

  11. [11]

    Capozzi, W

    F. Capozzi, W. Giarè, E. Lisi, A. Marrone, A. Melchiorri and A. Palazzo,Neutrino masses and mixing: Entering the era of subpercent precision,Phys. Rev. D111(2025) 093006 [2503.07752]

    work page arXiv 2025

  12. [12]

    2020 Global reassessment of the neutrino oscillation picture

    P.F. de Salas, D.V. Forero, S. Gariazzo, P. Martínez-Miravé, O. Mena, C.A. Ternes et al., 2020 global reassessment of the neutrino oscillation picture,JHEP02(2021) 071 [2006.11237]

    work page internal anchor Pith review arXiv 2020

  13. [13]

    The angular distribution of the reaction $\bar{\nu}_e + p \to e^+ + n$

    P. Vogel and J.F. Beacom,Angular Distribution of Neutron Inverse Beta Decay, ¯νe +p→e + +n,Phys. Rev.D60(1999) 053003 [hep-ph/9903554]. – 18 –

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  14. [14]

    Precise quasielastic neutrino/nucleon cross section

    A. Strumia and F. Vissani,Precise quasielastic neutrino/nucleon cross-section,Phys. Lett. B 564(2003) 42 [astro-ph/0302055]

    work page internal anchor Pith review Pith/arXiv arXiv 2003

  15. [15]

    Tomalak,Theory of inverse beta decay for reactor antineutrinos,2512.07956

    O. Tomalak,Theory of inverse beta decay for reactor antineutrinos,2512.07956

    work page arXiv

  16. [16]

    Neutrino mass hierarchy and electron neutrino oscillation parameters with one hundred thousand reactor events

    F. Capozzi, E. Lisi and A. Marrone,Neutrino mass hierarchy and electron neutrino oscillation parameters with one hundred thousand reactor events,Phys. Rev. D89(2014) 013001 [1309.1638]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  17. [17]

    Y.-F. Li, A. Wang, Y. Xu and J.-y. Zhu,Terrestrial Matter Effects on Reactor Antineutrino Oscillations: Constant vs. Fluctuated Density Profiles,2511.15702. [21]JUNOcollaboration,Initial performance results of the JUNO detector,2511.14590

    work page arXiv

  18. [18]

    X. Ji, W. Gu, X. Qian, H. Wei and C. Zhang,Combined Neyman–Pearson chi-square: An improved approximation to the Poisson-likelihood chi-square,Nucl. Instrum. Meth. A961 (2020) 163677 [1903.07185]

    work page arXiv 2020

  19. [19]

    Updated bounds on the (1,2) neutrino oscillation parameters after first JUNO results,

    F. Capozzi, E. Lisi, F. Marcone, A. Marrone and A. Palazzo,Updated bounds on the (1,2) neutrino oscillation parameters after first JUNO results,2511.21650. [24]SNO+collaboration,Measurement of reactor antineutrino oscillations with 1.46 ktonne-years of data at SNO+,2511.11856. [25](IceCube Collaboration)*, IceCubecollaboration,Measurement of atmospheric n...

    work page arXiv 2023

  20. [20]

    Quantifying the sensitivity of oscillation experiments to the neutrino mass ordering

    M. Blennow, P. Coloma, P. Huber and T. Schwetz,Quantifying the sensitivity of oscillation experiments to the neutrino mass ordering,JHEP03(2014) 028 [1311.1822]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  21. [21]

    X. Qian, A. Tan, W. Wang, J.J. Ling, R.D. McKeown and C. Zhang,Statistical Evaluation of Experimental Determinations of Neutrino Mass Hierarchy,Phys. Rev. D86(2012) 113011 [1210.3651]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  22. [22]

    Capozzi, E

    F. Capozzi, E. Lisi and A. Marrone,Mapping reactor neutrino spectra from TAO to JUNO, Phys. Rev. D102(2020) 056001 [2006.01648]. [29]JUNOcollaboration,TAO Conceptual Design Report: A Precision Measurement of the Reactor Antineutrino Spectrum with Sub-percent Energy Resolution,2005.08745. [30]Daya Baycollaboration,Comprehensive Measurement of the Reactor A...

    work page arXiv 2020

  23. [23]

    On the determination of anti-neutrino spectra from nuclear reactors

    P. Huber,On the determination of anti-neutrino spectra from nuclear reactors,Phys.Rev. C84(2011) 024617 [1106.0687]. – 19 –

    work page internal anchor Pith review Pith/arXiv arXiv 2011