Pith. Sign in

REVIEW

Measurement of Neutrino Oscillation Parameters at JUNO and Indication of the Neutrino Mass Ordering

High Energy Physics - Experiment · Submitted 29 Sep 2026 · 51 references

Pith's verdict: The solar-parameter result is new and solid; the ~1% Δm²₃₁ headline is a Daya Bay-anchored hybrid and the JUNO attribution is asserted, not demonstrated — still an honest, well-built analysis that deserves a serious referee.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 30, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.37895 (version 1)
Pith ID3G2OGDQS, Pith's permanent record of this paper and its review
CitationThe JUNO Collaboration. "Measurement of Neutrino Oscillation Parameters at JUNO and Indication of the Neutrino Mass Ordering." arXiv:2609.37895 (2026). https://pith.science/paper/3G2OGDQS
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

JUNO's reactor data measure the atmospheric neutrino mass splitting to roughly 1 percent precision, and, combined with long-baseline constraints, favor normal neutrino mass ordering at 2.1 sigma.

The reading

Using 207.2 live days of reactor antineutrino data at an average baseline of 52.5 km, JUNO reports a first measurement of the atmospheric mass splitting from its own data, reaching roughly 1% precision after the Daya Bay constraint selects the physical minimum under each mass-ordering hypothesis. The paper also improves the world's most precise values of the solar parameters, $\sin^2\theta_{12} = 0.3036 \pm 0.0064$ and $\Delta m^2_{21} = (7.388 \pm 0.078) \times 10^{-5}\,\mathrm{eV}^2$, at 2.1% and 1.1% relative precision. Adding the NuFit long-baseline $\Delta m^2_{31}$ constraints, the joint fit favors normal over inverted ordering with a one-sided significance of at least $2.1\sigma$ and a Bayes factor of 10.18. The claims matter because a precision reactor value of the atmospheric splitting is the key input for settling the neutrino mass ordering, and because the 1% precision and the ordering indication both depend on an external constraint resolving the degeneracy of the JUNO-only fit.

What carries the argument

The argument is carried by a joint binned-$\chi^2$ fit to three datasets — JUNO's 8,294 inverse $\beta$-decay candidates, 90,280 IBD candidates from the TAO near detector, and the Daya Bay near-detector spectra — each with its own covariance matrix plus penalty terms for nuisance parameters. Statistical bin variances use the Combined Neyman–Pearson prescription to reduce Poisson bias in the $\chi^2$ definition, and an energy-response matrix maps deposited to reconstructed energy, folding in positron non-linearity, radius-dependent resolution, and $^{14}$C pile-up. Two external pulls do the decisive work: the Daya Bay $|\Delta m^2_{31}|$ measurement (along with the $\sin^2\theta_{13}$ constraint) selects the compatible local minimum under each ordering hypothesis, and the NuFit long-baseline $\Delta m^2_{31}$ constraint drives the ordering comparison through the tension between reactor and accelerator values, with the NO/IO mapping of $\Delta m^2_{31}$ modulated by the CP phase and $\theta_{23}$: $\Delta m^2_{31,\mathrm{IO}} = -\Delta m^2_{31,\mathrm{NO}} + \Delta m^2_{21}[2\cos^2\theta_{12} - \sin(2\theta_{12})\sin\theta_{13}(\tan\theta_{23,\mathrm{NO}}\cos\delta_{\mathrm{CP,NO}} + \tan\theta_{23,\mathrm{IO}}\cos\delta_{\mathrm{CP,IO}})]$. Feldman–Cousins intervals give the quoted $1\sigma$ errors, and pseudo-experiments generate the frequentist significance.

Extended reading notes

Core claim

The central claim is that JUNO's reactor antineutrino spectrum, fit jointly with the Taishan Antineutrino Observatory near-detector reference and the Daya Bay near-detector data, yields a measurement of the atmospheric mass splitting of $\Delta m^2_{31} = +2.509^{+0.027}_{-0.025} \times 10^{-3}\,\mathrm{eV}^2$ under normal ordering and $-2.482^{+0.025}_{-0.026} \times 10^{-3}\,\mathrm{eV}^2$ under inverted ordering, at approximately 1% precision. Because the JUNO-only profile shows multiple degenerate local minima in $\Delta m^2_{31}$ (acknowledged in the text, e.g., sentences b21 and b187), the physical minimum under each hypothesis is selected by an external pull from the Daya Bay $|\Delta m^2_{31}|$ measurement, with the final precision set by the JUNO data. The same fit improves the solar oscillation parameters to $\sin^2\theta_{12} = 0.3036 \pm 0.0064$ and $\Delta m^2_{21} = (7.388 \pm 0.078) \times 10^{-5}\,\mathrm{eV}^2$, and after adding the NuFit long-baseline $\Delta m^2_{31}$ constraints it discriminates between orderings with $\Delta\chi^2 = 4.86$, a p-value of 0.0197 ($2.1\sigma$ one-sided, $2.3\sigma$ two-sided), and a Bayes factor of 10.18 in favor of normal ordering.

Load-bearing premise

The reported $\Delta m^2_{31}$ values and the mass-ordering preference rest on the assumption that the Daya Bay $|\Delta m^2_{31}|$ constraint correctly selects the unique physical minimum under each ordering hypothesis, since the JUNO-only profile has multiple degenerate local minima.

Editorial extensions

If this is right

  • If correct, the roughly 1% measurement of $\Delta m^2_{31}$ and the improved solar parameters give the tightest single-experiment handle on the three-flavor oscillation parameters from a reactor experiment.
  • The $2.1\sigma$ preference for normal ordering, if it persists with more data, sharpens the path toward a high-significance mass-ordering determination from JUNO's own data.
  • The world's largest geoneutrino candidate sample establishes the JUNO site as a geoneutrino observatory, complementing the oscillation measurement.
  • The precision $\sin^2\theta_{12}$ value tightens tests of lepton-flavor mixing patterns and of new-physics deviations beyond the three-flavor oscillation paradigm.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Because the ordering significance is driven mainly by the reactor-versus-accelerator tension in $\Delta m^2_{31}$, a future shift of $|\Delta m^2_{31}|$ by only a few per mille in the Daya Bay or long-baseline constraints would substantially weaken the reported conclusion; the claim is more fragile than its headline precision suggests.
  • The fastest testable extension is a dedicated study of which local minima the unconstrained JUNO fit selects and whether the Daya Bay-selected minimum is stable under alternative reactor-spectrum models, since that degeneracy is the load-bearing uncertainty.
  • As more TAO live time accumulates, the constraining power on the unoscillated reactor spectrum grows, which should reduce the spectral-shape systematics that currently limit the $\Delta m^2_{31}$ precision and eventually let JUNO alone resolve the minimum.
Original abstract

We report an improved measurement of neutrino oscillation parameters at the Jiangmen Underground Neutrino Observatory (JUNO) and the consequent indication of neutrino mass ordering. The analysis is based on 207.2 live days of data, comprising 8,294 inverse beta-decay candidates collected at an average reactor baseline of 52.5 km. A world-leading precision of approximately 1% was achieved for the atmospheric neutrino mass splitting $\Delta m^2_{31}$, with $\Delta m^2_{31}=+2.509^{+0.027}_{-0.025}\times10^{-3}\,\mathrm{eV}^2$ for normal ordering and $\Delta m^2_{31}=-2.482^{+0.025}_{-0.026}\times10^{-3}\,\mathrm{eV}^2$ for inverted ordering, using the Daya Bay constraint in the fit to resolve the correct physical minimum under each mass-ordering hypothesis. We further improved the world's most precise measurements of the solar oscillation parameters, obtaining $\sin^2\theta_{12}=0.3036\pm0.0064$ and $\Delta m^2_{21}=(7.388\pm0.078)\times10^{-5}\,\mathrm{eV}^2$, corresponding to relative precisions of 2.1% and 1.1%, respectively. Incorporating constraints on $\Delta m^2_{31}$ from long-baseline accelerator neutrino experiments into our analysis favored normal ordering over inverted ordering, with a one-sided significance of at least $2.1\sigma$ across all possible CP-violating phases; a Bayesian analysis yielded a Bayes factor of 10.18 in favor of normal ordering.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-30: 1 request, 4,981 input tokens.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.20 · human participants 0.03 · controlled comparison 0.07 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.87 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We report an improved measurement of neutrino oscillation parameters at the Jiangmen Underground Neutrino Observatory (JUNO) and the consequent indication of neutrino mass ordering. The analysis is based on 207.2 live days of data, comprising 8,294 inverse beta-decay candidates collected at an average reactor baseline of 52.5 km. A world-leading precision of approximately 1% was achieved for the atmospheric neutrino mass splitting $\Delta m^2_{31}$, with $\Delta m^2_{31}=+2.509^{+0.027}_{-0.025}\times10^{-3}\,\mathrm{eV}^2$ for normal ordering and $\Delta m^2_{31}=-2.482^{+0.025}_{-0.026}\times10^{-3}\,\mathrm{eV}^2$ for inverted ordering, using the Daya Bay constraint in the fit to resolve the correct physical minimum under each mass-ordering hypothesis. We further improved the world's most precise measurements of the solar oscillation parameters, obtaining $\sin^2\theta_{12}=0.3036\pm0.0064$ and $\Delta m^2_{21}=(7.388\pm0.078)\times10^{-5}\,\mathrm{eV}^2$, corresponding to relative precisions of 2.1% and 1.1%, respectively. Incorporating constraints on $\Delta m^2_{31}$ from long-baseline accelerator neutrino experiments into our analysis favored normal ordering over inverted ordering, with a one-sided significance of at least $2.1\sigma$ across all possible CP-violating phases; a Bayesian analysis yielded a Bayes factor of 10.18 in favor of normal ordering.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We report an improved measurement of neutrino oscillation parameters at the Jiangmen Underground Neutrino Observatory (JUNO) and the consequent indication of neutrino mass ordering. contribution 0.96 0.96 0.80 asserted 0.97
s1 The analysis is based on 207.2 live days of data, comprising 8,294 inverse beta-decay candidates collected at an average reactor baseline of 52.5 km. method 1.00 0.94 0.02 asserted 1.00
s2 A world-leading precision of approximately 1% was achieved for the atmospheric neutrino mass splitting $\Delta m^2_{31}$, with $\Delta m^2_{31}=+2.509^{+0.027}_{-0.025}\times10^{-3}\,\mathrm{eV}^2$ for normal ordering and $\Delta m^2_{31}=-2.482^{+0.025}_{-0.026}\times10^{-3}\,\mathrm{eV}^2$ for inverted ordering, using the Daya Bay constraint in the fit to resolve the correct physical minimum under each mass-ordering hypothesis. result 1.00 0.95 0.13 asserted 1.00
s3 We further improved the world's most precise measurements of the solar oscillation parameters, obtaining $\sin^2\theta_{12}=0.3036\pm0.0064$ and $\Delta m^2_{21}=(7.388\pm0.078)\times10^{-5}\,\mathrm{eV}^2$, corresponding to relative precisions of 2.1% and 1.1%, respectively. result 1.00 0.96 0.07 asserted 1.00
s4 Incorporating constraints on $\Delta m^2_{31}$ from long-baseline accelerator neutrino experiments into our analysis favored normal ordering over inverted ordering, with a one-sided significance of at least $2.1\sigma$ across all possible CP-violating phases; a Bayesian analysis yielded a Bayes factor of 10.18 in favor of normal ordering. result 0.93 0.94 0.54 asserted 0.97

What kind of paper the classifier reports this is: empirical study 1.00 · evidence from observation or experiment 0.92.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. b. The positron curve used directly in the oscil￾lation analysis was constrained to 0.6% precision across the relevant energy range. The energy resolution for the 68Ge γ-rays at the detec￾tor center was approximately 3.5% (σ/E), slightly worse than the MC prediction [27]. The γ-resolution curve shown in Fig. 2a was obtained by fitting calibration data from γ sources placed at the detector center. From these same data, the detector-center e + resolution model was subsequently derived, using a Geant4-based correction 0 1 2 3 4 5 6 7 Erec (MeV) ¡2 ¡1 0 1 2 (Erec ¡ E 0 rec )/E
Figure 4. FIG. 4. JUNO ∆
Figure 6. FIG. 6. Allowed regions in the sin

Showing 3 of the paper's 8 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.37895 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.