Recognition: no theorem link
Revival of the Reactor Antineutrino Anomaly
Pith reviewed 2026-05-12 04:56 UTC · model grok-4.3
The pith
The 2023 summation model revives the reactor antineutrino anomaly at 2.2 sigma.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The Reactor Antineutrino Anomaly is revived at the level of 2.2 sigma with the 2023 French summation model flux calculation. A global fit including gallium data yields 3.8 sigma tension that reduces to 1.3 sigma with enlarged gallium uncertainties. The usual explanation via active-sterile neutrino oscillations is considered, along with constraints from solar neutrino and KATRIN bounds, short-baseline reactor spectral ratios, and Daya Bay data.
What carries the argument
The 2023 summation model for calculating reactor antineutrino fluxes with its comprehensive uncertainty budget, and the global statistical fit to short-baseline neutrino oscillation parameters.
If this is right
- Reactor data alone indicate a 2.2 sigma average rate deficit relative to the 2023 prediction.
- Constraints on active-sterile neutrino oscillation parameters are derived from reactor spectral ratios and Daya Bay.
- Inclusion of gallium source data produces 3.8 sigma global tension under standard uncertainties.
- Enlarging gallium uncertainties reduces the overall tension to 1.3 sigma while retaining solar and KATRIN bounds.
- The simplest sterile neutrino explanation is tested against the full short-baseline dataset.
Where Pith is reading between the lines
- Future reactor experiments could distinguish between flux model errors and oscillation signals by using multiple independent flux predictions.
- Persistent tension after uncertainty adjustments may require reexamination of gallium experiment systematics or non-oscillation explanations for the deficits.
- The revival of the anomaly through one specific model highlights the value of cross-checking flux calculations to reduce theoretical uncertainty.
- Unified sterile neutrino models could be tested by linking this revived reactor tension to other short-baseline anomalies.
Load-bearing premise
The 2023 summation model provides a more accurate central value and uncertainty budget than the 2021 models, and the observed rate deficit is not dominated by unaccounted experimental systematics in the reactor measurements.
What would settle it
A new independent reactor antineutrino flux calculation that matches the measured average rates within 1 sigma, or gallium measurements with uncertainties small enough to keep the tension high without enlargement.
Figures
read the original abstract
The Reactor Antineutrino Anomaly refers to the deficit observed between the average event rate measured in reactor antineutrino experiments with respect to the theoretical prediction. This anomaly was first identified in 2011 ($2.5\sigma$) as a consequence of the Huber-Muller reactor antineutrino flux calculation. It was thought to be resolved in 2021 as a result of new reactor antineutrino flux calculations, with a reduction to about $1 \sigma$. In this work, we examine the latest reactor antineutrino flux calculation published in 2023 by a French research group. This work represents the first summation model to include a comprehensive uncertainty budget. The result indicates a revival of the Reactor Antineutrino Anomaly at the level of $2.2\,\sigma$. We also consider the usual simplest explanation of the Reactor Antineutrino Anomaly by active-sterile neutrino oscillations. We present the constraints on the oscillation parameters and we derive a tension of $3.8\sigma$ with the results of gallium source experiments (Gallium Anomaly) taking into account also the solar neutrino and KATRIN bounds, that of the combined short-baseline reactor spectral ratio measurements, and that of the Daya Bay search for a sub-eV sterile neutrino. Since the tension may be due to underestimated systematic uncertainties and the main tension is between the gallium data and the other data, we finally present the results of a global analysis with enlarged gallium uncertainties, which reduce the global tension to $1.3\sigma$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript re-evaluates the Reactor Antineutrino Anomaly using the 2023 French summation model for reactor antineutrino fluxes, which includes a comprehensive uncertainty budget. It reports a 2.2σ deficit relative to compiled reactor rate measurements, reviving the anomaly previously reduced to ~1σ with 2021 models. The paper then fits active-sterile neutrino oscillation parameters and derives a 3.8σ tension between this fit (incorporating reactor rates, spectral ratios, Daya Bay sub-eV search, solar neutrinos, and KATRIN) and gallium source experiments; enlarging the gallium uncertainties reduces the global tension to 1.3σ.
Significance. If the 2023 flux model is confirmed as more accurate, the 2.2σ revival would renew interest in light sterile neutrinos as an explanation for reactor deficits and provide updated global constraints on mixing parameters. The work's strength is its use of multiple independent published datasets (reactor rates, gallium, solar, KATRIN) and its explicit presentation of both standard and enlarged-uncertainty fits, which aids assessment of systematic effects. The result's broader impact depends on validation of the new flux calculation.
major comments (2)
- [Abstract] Abstract: The headline 2.2σ revival is obtained by adopting the 2023 French summation model as the reference prediction. The manuscript does not include an explicit comparison (e.g., a table of central flux values or uncertainty budgets) of the 2023 model against the 2021 calculations or against direct constraints from beta-decay spectra, which is load-bearing because a shift in the predicted central value or uncertainty would eliminate the reported deficit.
- [Global analysis section] Global fit and tension analysis: The 3.8σ tension is stated after combining reactor, gallium, solar, KATRIN, spectral-ratio, and Daya Bay data. The manuscript should provide the individual χ² contributions from each dataset and the precise prescription for combining the gallium anomaly with the other bounds, as these details are required to verify the tension value and the subsequent reduction to 1.3σ under enlarged gallium uncertainties.
minor comments (1)
- [Abstract] Abstract: The phrasing 'the usual simplest explanation of the Reactor Antineutrino Anomaly by active-sterile neutrino oscillations' is informal; a direct statement of the two-flavor oscillation hypothesis and the parameters being fitted would improve clarity.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and for the constructive comments. We address each major comment below and have revised the manuscript accordingly to improve transparency and verifiability.
read point-by-point responses
-
Referee: [Abstract] Abstract: The headline 2.2σ revival is obtained by adopting the 2023 French summation model as the reference prediction. The manuscript does not include an explicit comparison (e.g., a table of central flux values or uncertainty budgets) of the 2023 model against the 2021 calculations or against direct constraints from beta-decay spectra, which is load-bearing because a shift in the predicted central value or uncertainty would eliminate the reported deficit.
Authors: We agree that an explicit side-by-side comparison strengthens the presentation. In the revised manuscript we have added a new table (Table 1) in Section 2 that lists the predicted antineutrino yields per fission for the dominant isotopes (235U, 239Pu, 241Pu, 238U) together with their total uncertainties under the 2023 French summation model, the 2021 models, and the direct beta-decay spectrum constraints. The table shows that the central-value shift relative to 2021 is the dominant source of the revived 2.2σ deficit while the uncertainty budgets remain comparable; a short accompanying paragraph discusses the nuclear-data updates responsible for the shift. revision: yes
-
Referee: [Global analysis section] Global fit and tension analysis: The 3.8σ tension is stated after combining reactor, gallium, solar, KATRIN, spectral-ratio, and Daya Bay data. The manuscript should provide the individual χ² contributions from each dataset and the precise prescription for combining the gallium anomaly with the other bounds, as these details are required to verify the tension value and the subsequent reduction to 1.3σ under enlarged gallium uncertainties.
Authors: We accept that the numerical breakdown is necessary for independent verification. The revised global-analysis section now contains a dedicated table (Table 3) that reports the individual χ² contributions from each dataset (reactor rates, spectral ratios, Daya Bay sub-eV search, solar neutrinos, KATRIN, and gallium) evaluated at the global best-fit point, together with the total χ² and degrees of freedom. We have also added an explicit paragraph describing the combination procedure: the global χ² is the direct sum of the independent dataset χ² functions, with the gallium anomaly incorporated through its published rate measurements and covariance matrix; the enlarged-uncertainty case is obtained by uniformly scaling the gallium uncertainties by a constant factor chosen to reflect possible unaccounted systematics, which reduces the tension to 1.3σ as stated. revision: yes
Circularity Check
No circularity; central claims rest on external 2023 flux model and independent datasets
full rationale
The paper compares reactor antineutrino rate data to an external 2023 French summation-model flux calculation (not derived or fitted within the paper) to obtain the 2.2σ revival. Global fits for sterile-neutrino parameters and tensions with gallium/solar/KATRIN data are standard statistical combinations of independent experimental inputs; no key quantity (flux central value, uncertainty, or oscillation parameter) is defined in terms of itself or obtained by fitting a subset and relabeling the result as a prediction. Self-citations, if present, are not load-bearing for the headline results.
Axiom & Free-Parameter Ledger
free parameters (1)
- sterile neutrino mixing angle and mass-squared difference
axioms (2)
- domain assumption Three active neutrinos plus one sterile neutrino with mixing described by the standard 4x4 PMNS matrix extension
- domain assumption The 2023 French summation model provides the correct central flux value and full uncertainty budget
invented entities (1)
-
light sterile neutrino
no independent evidence
Forward citations
Cited by 1 Pith paper
-
Effective Matter Flavor Conversion Mediated by Pseudo-Sterile States as the Possible Origin of Neutrino Oscillation Anomalies
A sterile neutrino with a novel density-dependent matter potential Vs resolves multiple oscillation anomalies when Vs is negative and mixing angles are small.
Reference graph
Works this paper leans on
-
[1]
We obtained the 2σand 3σallowed regions in the (sin 22ϑee,∆m2
For each value of these param- eters,χ 2 M is minimized over the nuisance parameters ri. We obtained the 2σand 3σallowed regions in the (sin 22ϑee,∆m2
-
[2]
plane shown in blue in Fig.3. One can see that the 2σallowed region has a lower bound for sin 22ϑee and the 3σallowed region has only an upper bound. This was expected from the 2.2σCEA-RAA. In Fig.3 we have shown also the regions allowed by the data of the gallium source experiments in Ref.[6], the solar bound [7], the KATRIN bound [8], the com- bined bou...
-
[3]
The best fit for each group of experiments is indicated by a cross with the corresponding color
plane obtained from our neutrino oscillation fit of the reactor antineutrino data in Tab.I with the CEA model [46] and those obtained from the gallium neutrino data in Ref.[6], the solar neutrino data in Ref.[7], the KATRIN data [8], the combined DANSS [9], NEOS+RENO [10], PROSPECT [11], STEREO [12] reactor spectral ratio measurements (SPE) and the Daya B...
-
[4]
The equality is approximate because the Data Sets PG CEA-RAA+SOL+KATRIN 0.9σ CEA-RAA+SPE 1.6σ CEA-RAA+GA 2.8σ GA+SOL+KATRIN 3.0σ GA+SOL+KATRIN+SPE 3.8σ CEA-RAA+GA+SOL+KATRIN 3.1σ CEA-RAA+GA+SOL+KATRIN+SPE 4.0σ CEA-RAA+GA+SOL+KATRIN+SPE+DB 3.8σ CEA-RAA+GA+SOL+KATRIN+SPE+DB 1.3σ TABLE IV. Values of the parameter goodness of fit (PG) for different combinatio...
-
[5]
The best fit is indicated by a green star
plane obtained from the com- bination of the reactor (CEA-RAA) analysis and the bounds of the solar [7] (SOL), KATRIN [8], DANSS [9], NEOS+RENO [10], PROSPECT [11], STEREO [12] re- actor spectral ratio (SPE) experiments, the Daya Bay search for a sub-eV sterile neutrino [13] (DB), and the gallium [6] analysis with enlarged uncertainties ( GA) by the facto...
work page 2023
-
[6]
plane from the analy- sis of all the available data in the framework of 3+1 neutrino oscillations, assuming the reliability of the CEA model. Acknowledgments We thank Thierry Lasserre, Jiajie Ling and Soo- Bong Kim for useful communications about the KA- TRIN, Daya Bay and RENO experiments. This work was supported in part by the National Natu- ral Science...
-
[7]
S. Gariazzo, C. Giunti, M. Laveder, Y. F. Li, and E. M. Zavanin, J. Phys.G43, 033001 (2016), arXiv:1507.08204 [hep-ph]
work page Pith review arXiv 2016
-
[8]
C. Giunti and T. Lasserre, Ann. Rev. Nucl. Part. Sci.69, 163 (2019), arXiv:1901.08330 [hep-ph]
- [9]
- [10]
-
[11]
B. Dasgupta and J. Kopp, Phys.Rept.928, 63 (2021), arXiv:2106.05913 [hep-ph]
-
[12]
M. Cadeddu, N. Cargioli, G. Carotenuto, F. Dordei, L. Ferro, and C. Giunti, Phys.Rev.D113, 033006 (2026), arXiv:2507.13103 [hep-ph]
- [13]
-
[14]
H. Acharya et al. (KATRIN), Nature648, 70 (2025), arXiv:2503.18667 [hep-ex]
-
[15]
M. Danilov, (2022), talk presented at ICHEP 2022, 41th International Conference on High Energy Physics, 6-13 July 2022, Bologna, Italy
work page 2022
-
[16]
Z. Atif et al. (RENO, NEOS), Phys. Rev. D105, L111101 (2022), arXiv:2011.00896 [hep-ex]
-
[17]
Andriamirado et al., Phys.Rev.Lett.134, 151802 (2025), arXiv:2406.10408 [hep-ex]
M. Andriamirado et al., Phys.Rev.Lett.134, 151802 (2025), arXiv:2406.10408 [hep-ex]
-
[18]
H. Almazan et al. (STEREO), Nature613, 257 (2023), arXiv:2210.07664 [hep-ex]
- [19]
- [20]
- [21]
- [22]
- [23]
-
[24]
The Reactor Antineutrino Anomaly
G. Mention et al., Phys. Rev.D83, 073006 (2011), arXiv:1101.2755 [hep-ex]
work page Pith review arXiv 2011
- [25]
-
[26]
On the determination of anti-neutrino spectra from nuclear reactors,
P. Huber, Phys. Rev.C84, 024617 (2011), arXiv:1106.0687 [hep-ph]
-
[27]
F. V. Feilitzsch, A. A. Hahn, and K. Schreckenbach, Phys. Lett.B118, 162 (1982)
work page 1982
-
[28]
K. Schreckenbach, G. Colvin, W. Gelletly, and F. V. Feilitzsch, Phys. Lett.B160, 325 (1985)
work page 1985
- [29]
-
[30]
M. Estienne, M. Fallot, et al., Phys. Rev. Lett.123, 022502 (2019), arXiv:1904.09358 [nucl-ex]
-
[31]
V. Kopeikin, M. Skorokhvatov, and O. Titov, Phys.Rev.D104, L071301 (2021), arXiv:2103.01684 [nucl-ex]
- [32]
- [33]
-
[34]
A. A. Kuvshinnikov, L. A. Mikaelyan, S. V. Niko- laev, M. D. Skorokhvatov, and A. V. Etenko, JETP Lett.54, 253 (1991)
work page 1991
-
[35]
A. I. Afonin et al., Sov. Phys. JETP67, 213 (1988)
work page 1988
- [36]
- [37]
- [38]
-
[39]
A. Hoummada, S. L. Mikou, G. Bagieu, J. F. Cavaignac, and D. H. Koang, Applied Radiation and Isotopes46, 449 (1995)
work page 1995
-
[40]
G. S. Vidyakin et al. (Krasnoyarsk), Sov. Phys. JETP66, 243 (1987)
work page 1987
-
[41]
G. S. Vidyakin et al. (Krasnoyarsk), JETP Lett.59, 390 (1994)
work page 1994
- [42]
-
[43]
Z. D. Greenwood et al., Phys. Rev.D53, 6054 (1996)
work page 1996
- [44]
-
[45]
M. Apollonio et al. (CHOOZ), Eur. Phys. J.C27, 331 (2003), hep-ex/0301017
-
[46]
F. Boehm et al. (Palo Verde), Phys. Rev.D64, 112001 (2001), hep-ex/0107009
-
[47]
G. Boireau et al. (NUCIFER), Phys. Rev.D93, 112006 (2016), arXiv:1509.05610 [physics]
-
[48]
Z. Atif et al. (RENO), Phys.Rev.D104, L111301 (2021), arXiv:2010.14989 [hep-ex]
-
[49]
H. de Kerret et al. (Double Chooz), Nature Phys. 16, 558 (2020), arXiv:1901.09445 [hep-ex]
- [50]
- [51]
-
[52]
L. Perisse, A. Onillon, X. Mougeot, M. Vivier, T. Lasserre, A. Letourneau, D. Lhuillier, and G. Mention, Phys.Rev.C108, 055501 (2023), arXiv:2304.14992 [nucl-ex]
-
[53]
S. M. Bilenky and S. T. Petcov, Rev. Mod. Phys. 59, 671 (1987)
work page 1987
- [54]
- [55]
-
[56]
Giuntiet al., JHEP1, 164 (2022) (arXiv:2209.00916)
C. Giunti, Y. F. Li, C. A. Ternes, O. Tyagi, and Z. Xin, JHEP10, 164 (2022), arXiv:2209.00916 [hep-ph]
- [57]
-
[58]
C. Giunti and C. A. Ternes, Phys.Lett.B849, 138436 (2024), arXiv:2312.00565 [hep-ph]
- [59]
-
[60]
Testing the statistical compatibility of independent data sets
M. Maltoni and T. Schwetz, Phys. Rev.D68, 033020 (2003), hep-ph/0304176
work page Pith review arXiv 2003
- [61]
- [62]
-
[63]
Hayes et al., Phys.Rev.Lett.120, 022503 (2018), arXiv:1707.07728 [nucl-th]
A. Hayes et al., Phys.Rev.Lett.120, 022503 (2018), arXiv:1707.07728 [nucl-th]. 10
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.