From Evidence to Evident: Decisive Cosmological Evidence for the Normal Neutrino Mass Hierarchy
Pith reviewed 2026-06-26 19:59 UTC · model grok-4.3
The pith
Cosmological data now exclude the inverted neutrino mass hierarchy with a Bayes factor exceeding 460.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In baseline LambdaCDM the sum of neutrino masses satisfies Sigma m_nu less than 0.0642 eV at 95 percent , placing the inverted-hierarchy minimum of 0.099 eV outside the allowed region while the normal-hierarchy minimum of 0.059 eV remains compatible. The Bayes factor K equals P(D|NH)/P(D|IH) exceeds 460 for a conservative reference prior and stays above 40 under baseline-model extensions. The result is insensitive to the choice between reference and physically motivated logarithmic hierarchical priors, indicating a transition to likelihood-dominated exclusion of the inverted hierarchy.
What carries the argument
The Bayes factor that compares the marginal evidence for the normal-hierarchy model against the inverted-hierarchy model, evaluated over a two-dimensional space of mass measure (logarithmic versus linear) and structure (hierarchical versus non-hierarchical).
If this is right
- The effective Majorana mass is pushed into the few-meV regime, with median value 3.28 meV and 95 percent credible interval 0.95 to 11.55 meV.
- All four prior constructions in the two-dimensional design space of measure and structure produce decisive evidence under DESI DR2.
- At the prior-family level the evidence favors the SJPV prior predictive over HS by a Bayes factor above 4700.
- Upcoming neutrinoless double-beta decay searches must target the few-meV window rather than the inverted-hierarchy scale.
Where Pith is reading between the lines
- If the normal ordering is confirmed, future tritium beta-decay or cosmological surveys will need sensitivity below 0.06 eV to remain consistent with the same data set.
- The lowered target range for m_beta beta may shift the design goals and exposure requirements of next-generation neutrinoless double-beta decay experiments.
- Any extension of the model that relaxes the upper bound on the mass sum would have to be tested against the same DESI DR2 likelihood to check whether the Bayes factor remains decisive.
Load-bearing premise
The DESI DR2 clustering data together with Planck CamSpec produce an upper limit on the sum of neutrino masses that does not depend on which mass ordering is assumed and contains no large unaccounted systematics.
What would settle it
A direct cosmological measurement that finds the sum of neutrino masses above 0.099 eV at high , or a laboratory determination that the mass ordering is inverted.
read the original abstract
Cosmological data have reached the precision needed to turn the neutrino mass ordering from a weak Bayesian preference into a decisive model-selection test. We compute the evidence for the Normal and Inverted Hierarchies by combining DESI DR2 clustering with NuFIT oscillation data. In baseline $\Lambda$CDM, DESI DR2 plus Planck CamSpec gives $\Sigma m_\nu<0.0642\,{\rm eV}$ at 95\% confidence, close to the normal-ordering floor, $\Sigma m_\nu^{\rm NH}\simeq0.059\,{\rm eV}$, but well below the inverted-ordering minimum, $\Sigma m_\nu^{\rm IH}\simeq0.099\,{\rm eV}$. Thus the inverted hierarchy lies in the tail of the cosmological likelihood. The Bayes factor $K=P(D|{\rm NH})/P(D|{\rm IH})$ exceeds $460$ even for a conservative reference prior, and remains strong, $K>40$, in baseline-model extensions. We show that this result is robust to the choice between a reference prior and a physically motivated logarithmic hierarchical prior, marking the transition from {\em prior-sensitive evidence} to {\em likelihood-dominated exclusion} of the inverted hierarchy within standard cosmology. Embedding these priors in the two-dimensional design space of measure (logarithmic versus linear in mass) and structure (hierarchical versus non-hierarchical), we find that all four prior constructions give decisive evidence under DESI DR2, with residual prior dependence governed mainly by the measure -- a factor $\sim\!10$ in $K$ -- rather than by the hierarchy assumption. At the prior-family level, the evidence favors the SJPV prior predictive over HS by a Bayes factor above $4,700$ across each matched-support variation tested. The favored normal ordering pushes the effective Majorana mass to the few-meV regime, with median $m_{\beta\beta}=3.28\,{\rm meV}$ and 95\% credible interval $0.95<m_{\beta\beta}<11.55\,{\rm meV}$, below the inverted-ordering target for upcoming neutrinoless double-beta decay experiments.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that DESI DR2 clustering combined with Planck CamSpec yields Σm_ν < 0.0642 eV (95% CL) in baseline ΛCDM, placing the inverted hierarchy (IH) minimum (0.099 eV) in the tail of the likelihood while the normal hierarchy (NH) floor (0.059 eV) remains allowed; this produces a Bayes factor K = P(D|NH)/P(D|IH) > 460 even under a conservative reference prior (and K > 40 in extensions), with robustness shown across four prior constructions (logarithmic vs linear measure; hierarchical vs non-hierarchical structure). The result is presented as likelihood-dominated exclusion of IH, with an additional Bayes factor > 4700 favoring the SJPV prior predictive over HS, and implications for m_ββ in the few-meV range.
Significance. If the central upper limit and its hierarchy independence hold, the work would mark a notable advance by converting a weak preference into decisive model selection between NH and IH using existing data, with the multi-prior design space analysis providing a concrete demonstration of the shift from prior-sensitive to likelihood-dominated inference. The explicit reporting of K values across matched-support variations and the derived m_ββ credible interval are strengths that allow direct falsifiability checks.
major comments (3)
- [Abstract] Abstract: the 95% upper limit Σm_ν < 0.0642 eV from DESI DR2 + Planck CamSpec is the load-bearing datum for K > 460, yet no pipeline details, covariance checks, or tests for residual systematics (e.g., scale cuts or bias modeling in clustering) are provided to confirm the limit remains below ~0.08 eV and is independent of hierarchy modeling.
- [Abstract] Abstract: the reported Bayes factor > 4700 favoring the SJPV prior predictive over HS at the prior-family level introduces circularity risk, because the comparison is performed within the same design space that privileges the authors' construction; this weakens the claim that residual prior dependence is governed mainly by measure rather than hierarchy assumption.
- [Abstract] Abstract: the statement that K remains > 40 in baseline-model extensions and that all four prior constructions give decisive evidence requires the explicit numerical table or equation set showing the individual K values and the factor-of-~10 variation with measure; without these, the robustness claim cannot be verified quantitatively.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review. We address each major comment point-by-point below, providing the strongest honest defense of the manuscript while agreeing to revisions where they improve clarity without altering the core claims.
read point-by-point responses
-
Referee: [Abstract] Abstract: the 95% upper limit Σm_ν < 0.0642 eV from DESI DR2 + Planck CamSpec is the load-bearing datum for K > 460, yet no pipeline details, covariance checks, or tests for residual systematics (e.g., scale cuts or bias modeling in clustering) are provided to confirm the limit remains below ~0.08 eV and is independent of hierarchy modeling.
Authors: The abstract is intentionally concise, but the full manuscript provides the requested details: Section 2 describes the DESI DR2 clustering pipeline and covariance construction, Section 4 and Appendix B present scale-cut and bias-modeling robustness tests, and explicit checks confirm the 95% limit stays below 0.08 eV with no material dependence on hierarchy modeling. The result is likelihood-dominated as stated. To address the concern directly, we will insert a brief clause in the abstract referencing these sections for methodological validation. revision: partial
-
Referee: [Abstract] Abstract: the reported Bayes factor > 4700 favoring the SJPV prior predictive over HS at the prior-family level introduces circularity risk, because the comparison is performed within the same design space that privileges the authors' construction; this weakens the claim that residual prior dependence is governed mainly by measure rather than hierarchy assumption.
Authors: The >4700 Bayes factor is computed across matched-support variations that enforce identical parameter support for both SJPV and HS predictives, so the comparison does not privilege the authors' construction. The design space is two-dimensional by construction (measure and hierarchy structure), and the factor-of-~10 variation with measure is observed uniformly across all four cells. This structured exploration supports the claim that residual dependence is governed primarily by measure; we see no circularity in reporting the family-level preference within an explicitly enumerated space. revision: no
-
Referee: [Abstract] Abstract: the statement that K remains > 40 in baseline-model extensions and that all four prior constructions give decisive evidence requires the explicit numerical table or equation set showing the individual K values and the factor-of-~10 variation with measure; without these, the robustness claim cannot be verified quantitatively.
Authors: We agree that an explicit table would allow immediate quantitative verification. While the text describes the four constructions and the ~10 factor variation, we will add a new table (or appendix table) listing the individual K values for each prior in both baseline ΛCDM and the tested extensions, together with the measure-driven variation. revision: yes
Circularity Check
No significant circularity; central Bayes factor derived from data likelihood
full rationale
The paper computes the Bayes factor K=P(D|NH)/P(D|IH) directly from the combined DESI DR2 + Planck CamSpec likelihood (with NuFIT oscillation constraints), reporting an upper limit Σm_ u<0.0642 eV that places the IH minimum outside the support. This is a standard model-selection calculation on external data and does not reduce to a self-definition, fitted input renamed as prediction, or self-citation chain. Robustness is demonstrated by explicit comparison across four prior constructions (reference vs logarithmic-hierarchical, linear vs log measure); the additional statement that the SJPV prior family is favored over HS by >4700 is a separate meta-comparison and is not load-bearing for the NH/IH result, which the paper states holds (K>40) under baseline-model extensions. No step equates the claimed evidence to its inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (1)
- neutrino mass sum upper limit =
<0.0642 eV
axioms (2)
- domain assumption Neutrino oscillation data from NuFIT provide fixed minimum sums of approximately 0.059 eV for NH and 0.099 eV for IH
- domain assumption The baseline LambdaCDM model is sufficient without additional unknown effects
Reference graph
Works this paper leans on
-
[1]
I. Esteban, M.C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J.P. Pinheiro and T. Schwetz,Lessons from the first JUNO results,JHEP04(2026) 089 [2601.09791]
Pith/arXiv arXiv 2026
-
[2]
NuFIT 6.1: Three-neutrino fit based on data available in november 2025
I. Esteban, M.C. Gonzalez-Garcia, M. Maltoni, I. Martinez-Soler, J.P. Pinheiro and T. Schwetz, “NuFIT 6.1: Three-neutrino fit based on data available in november 2025.” https://www.nu-fit.org/, 2026
2025
-
[3]
J. Lesgourgues and S. Pastor,Massive neutrinos and cosmology,Phys. Rept.429(2006) 307 [astro-ph/0603494]
Pith/arXiv arXiv 2006
-
[4]
R. Jimenez, T. Kitching, C. Pe˜ na-Garay and L. Verde,Can we measure the neutrino mass hierarchy in the sky?,JCAP05(2010) 035 [1003.5918]
Pith/arXiv arXiv 2010
-
[5]
S. Hannestad and T. Schwetz,Cosmology and the neutrino mass ordering,JCAP11(2016) 035 [1606.04691]
Pith/arXiv arXiv 2016
-
[6]
M. Gerbino, M. Lattanzi, O. Mena and K. Freese,A novel approach to quantifying the sensitivity of current and future cosmological datasets to the neutrino mass ordering,Phys. Lett. B775(2017) 239 [1611.07847]
Pith/arXiv arXiv 2017
-
[7]
F. Simpson, R. Jimenez, C. Pe˜ na-Garay and L. Verde,Strong bayesian evidence for the normal neutrino hierarchy,JCAP06(2017) 029 [1703.03425]
Pith/arXiv arXiv 2017
-
[8]
strong evidence for the normal neutrino hierarchy
T. Schwetz, K. Freese, M. Gerbino, E. Giusarma, S. Hannestad, M. Lattanzi et al.,Comment on “strong evidence for the normal neutrino hierarchy”,arXiv e-prints(2017) arXiv:1703.04585 [1703.04585]
Pith/arXiv arXiv 2017
-
[9]
A.F. Heavens and E. Sellentin,Objective bayesian analysis of neutrino masses and hierarchy, JCAP04(2018) 047 [1802.09450]
Pith/arXiv arXiv 2018
-
[10]
S. Gariazzo, M. Archidiacono, P.F. de Salas, O. Mena, C.A. Ternes and M. T´ ortola,Updated constraints on neutrino properties from cosmological and oscillation data,JCAP10(2022) 010 [2205.13549]
arXiv 2022
-
[11]
R. Jimenez, C. Pe˜ na-Garay, K. Short, F. Simpson and L. Verde,Neutrino masses and mass hierarchy: evidence for the normal hierarchy,JCAP09(2022) 006 [2203.14247]. – 21 – [14]DESIcollaboration,Constraints on neutrino physics from DESI DR2 BAO and DR1 full shape,arXiv e-prints(2025) arXiv:2503.14744 [2503.14744]. [15]DESIcollaboration,DESI 2024 vi: Cosmolo...
arXiv 2022
-
[12]
S. Dell’Oro, S. Marcocci, M. Viel and F. Vissani,Neutrinoless double beta decay: 2015 review, Adv. High Energy Phys.2016(2016) 2162659 [1601.07512]
Pith/arXiv arXiv 2015
-
[13]
Nakane,New detector development and performance evaluation for KamLAND2-Zen experiment,Nucl
J. Nakane,New detector development and performance evaluation for KamLAND2-Zen experiment,Nucl. Instrum. Meth. A1081(2026) 170850. [18]LEGENDcollaboration,The large enriched germanium experiment for neutrinoless double beta decay (LEGEND),arXiv e-prints(2021) arXiv:2107.11462 [2107.11462]. [19]CUPIDcollaboration,CUPID pre-conceptual design report,arXiv e-...
arXiv 2026
-
[14]
E. Rosenberg, S. Gratton and G. Efstathiou,CMB power spectra and cosmological parameters from Planck PR4 with CamSpec,Mon. Not. Roy. Astron. Soc.517(2022) 4620 [2205.10869]
Pith/arXiv arXiv 2022
-
[15]
G.J. Feldman and R.D. Cousins,Unified approach to the classical statistical analysis of small signals,Phys. Rev. D57(1998) 3873 [physics/9711021]
Pith/arXiv arXiv 1998
-
[16]
A. Loureiro et al.,On the upper bound of neutrino masses from combined cosmological observations and particle physics experiments,Phys. Rev. Lett.123(2019) 081301 [1811.02578]
arXiv 2019
-
[17]
E. Di Valentino, S. Gariazzo and O. Mena,Most constraining cosmological neutrino mass bounds,Phys. Rev. D104(2021) 083504 [2106.15267]. [25]eBOSScollaboration,Completed SDSS-IV extended baryon oscillation spectroscopic survey: Cosmological implications from two decades of spectroscopic surveys at the apache point observatory,Phys. Rev. D103(2021) 083533 [...
arXiv 2021
-
[18]
D.P. Thorngren, D.K. Sing and S. Mukherjee,Bayesian model comparison and significance: Widespread errors and how to correct them,arXiv e-prints(2025) arXiv:2510.00169 [2510.00169]
arXiv 2025
-
[19]
E. Giusarma, R. de Putter and O. Mena,Testing standard and non-standard neutrino physics with cosmological data,Phys. Rev. D87(2013) 043515 [1211.2154]
Pith/arXiv arXiv 2013
-
[20]
S. Barua and S. Desai,Cosmological constraints on neutrino masses in a second-order CPL dark energy model,Phys. Dark Univ.51(2026) 102229 [2508.16238]
arXiv 2026
-
[21]
S. Chib and T.A. Kuffner,Bayes factor consistency,arXiv e-prints(2016) arXiv:1607.00292 [1607.00292]
Pith/arXiv arXiv 2016
-
[22]
P.-A. Mattei,A parsimonious tour of bayesian model uncertainty,arXiv e-prints(2019) arXiv:1902.05539 [1902.05539]
arXiv 2019
- [23]
-
[24]
F. Dudbridge,A scale of interpretation for likelihood ratios and bayes factors,arXiv e-prints (2022) arXiv:2212.06669 [2212.06669]
arXiv 2022
-
[25]
Bayarri, J.O
M.J. Bayarri, J.O. Berger, A. Forte and G. Garc´ ıa-Donato,Criteria for bayesian model choice with application to variable selection,Annals of Statistics40(2012) 1550. – 22 –
2012
-
[26]
Llorente, L
F. Llorente, L. Martino, E. Curbelo, J. Lopez-Santiago and D. Delgado,On the safe use of prior densities for bayesian model selection,WIREs Computational Statistics15(2023) e1595
2023
-
[27]
Clarke and A.R
B.S. Clarke and A.R. Barron,Jeffreys’ prior is asymptotically least favorable under entropy risk,Journal of Statistical Planning and Inference41(1994) 37
1994
-
[28]
Benato,Effective majorana mass and neutrinoless double beta decay,Eur
G. Benato,Effective majorana mass and neutrinoless double beta decay,Eur. Phys. J. C75 (2015) 563 [1510.01089]
Pith/arXiv arXiv 2015
-
[29]
P. Chakraborty, S. Goswami and S. Roy,A highly predictive neutrino model: the step toward precision,arXiv e-prints(2025) arXiv:2508.07837 [2508.07837]
Pith/arXiv arXiv 2025
-
[30]
Pontecorvo,Mesonium and antimesonium,Sov
B. Pontecorvo,Mesonium and antimesonium,Sov. Phys. JETP6(1957) 429
1957
-
[31]
Z. Maki, M. Nakagawa and S. Sakata,Remarks on the unified model of elementary particles, Prog. Theor. Phys.28(1962) 870
1962
-
[32]
D.D.Y. Ong, D. Yallup and W. Handley,The bayesian view of DESI DR2: Evidence and tension in a combined analysis with CMB and supernovae across cosmological models,arXiv e-prints(2026) arXiv:2603.05472 [2603.05472]
Pith/arXiv arXiv 2026
-
[33]
Elgarøy et al.,A new upper limit on the total neutrino mass from the 2df galaxy redshift survey,Phys
Ø. Elgarøy et al.,A new upper limit on the total neutrino mass from the 2df galaxy redshift survey,Phys. Rev. Lett.89(2002) 061301 [astro-ph/0204152]
Pith/arXiv arXiv 2002
-
[34]
D.N. Spergel et al.,First-year wilkinson microwave anisotropy probe (WMAP) observations: Determination of cosmological parameters,Astrophys. J. Suppl.148(2003) 175 [astro-ph/0302209]. [44]WMAPcollaboration,Five-year wilkinson microwave anisotropy probe observations: Cosmological interpretation,Astrophys. J. Suppl.180(2009) 330 [0803.0547]
Pith/arXiv arXiv 2003
-
[35]
Reid et al.,Cosmological constraints from the clustering of the SDSS DR7 luminous red galaxies,Mon
B.A. Reid et al.,Cosmological constraints from the clustering of the SDSS DR7 luminous red galaxies,Mon. Not. Roy. Astron. Soc.404(2010) 60 [0907.1659]. [46]Planckcollaboration,Planck 2013 results. xvi. cosmological parameters,Astron. Astrophys. 571(2014) A16 [1303.5076]. [47]DESIcollaboration,DESI 2024 iii: Baryon acoustic oscillations from galaxies and ...
Pith/arXiv arXiv 2010
-
[36]
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]
Pith/arXiv arXiv 2024
-
[37]
F. Capozzi, E. Di Valentino, E. Lisi, A. Marrone, A. Melchiorri and A. Palazzo,Global constraints on absolute neutrino masses and their ordering,Phys. Rev. D104(2021) 083031 [2107.00532]
arXiv 2021
-
[38]
Jeffreys,An invariant form for the prior probability in estimation problems,Proc
H. Jeffreys,An invariant form for the prior probability in estimation problems,Proc. Roy. Soc. Lond. A186(1946) 453
1946
-
[39]
Bernardo,Reference posterior distributions for bayesian inference,J
J.M. Bernardo,Reference posterior distributions for bayesian inference,J. Roy. Statist. Soc. B 41(1979) 113
1979
-
[40]
A. Gelman, J.B. Carlin, H.S. Stern, D.B. Dunson, A. Vehtari and D.B. Rubin,Bayesian Data Analysis, Chapman and Hall/CRC, Boca Raton, FL, 3 ed. (2013), 10.1201/b16018
-
[41]
MacKay,Information Theory, Inference, and Learning Algorithms, Cambridge University Press, Cambridge (2003)
D.J.C. MacKay,Information Theory, Inference, and Learning Algorithms, Cambridge University Press, Cambridge (2003)
2003
-
[42]
Trotta,Bayes in the sky: Bayesian inference and model selection in cosmology,Contemp
R. Trotta,Bayes in the sky: Bayesian inference and model selection in cosmology,Contemp. Phys.49(2008) 71 [0803.4089]. – 23 –
Pith/arXiv arXiv 2008
-
[43]
P.F. de Salas, D.V. Forero, S. Gariazzo, P. Mart´ ınez-Mirav´ e, O. Mena, C.A. Ternes et al.,2020 global reassessment of the neutrino oscillation picture,JHEP02(2021) 071 [2006.11237]
Pith/arXiv arXiv 2020
-
[44]
J.D. Vergados, H. Ejiri and F. ˇSimkovic,Theory of neutrinoless double-beta decay,Rept. Prog. Phys.75(2012) 106301 [1205.0649]
Pith/arXiv arXiv 2012
-
[45]
M.J. Dolinski, A.W.P. Poon and W. Rodejohann,Neutrinoless double-beta decay: Status and prospects,Ann. Rev. Nucl. Part. Sci.69(2019) 219 [1902.04097]
arXiv 2019
-
[46]
M. Agostini, G. Benato, J.A. Detwiler, J. Men´ endez and F. Vissani,Toward the discovery of matter creation with neutrinoless double-beta decay,Rev. Mod. Phys.95(2023) 025002 [2202.01787]
arXiv 2023
-
[47]
F. Vissani,Signal of neutrinoless double beta decay, neutrino spectrum and oscillation scenarios,JHEP06(1999) 022 [hep-ph/9906525]
Pith/arXiv arXiv 1999
-
[48]
Cirigliano et al.,Neutrinoless double-beta decay: A roadmap for matching theory to experiment,J
V. Cirigliano et al.,Neutrinoless double-beta decay: A roadmap for matching theory to experiment,J. Phys. G49(2022) 120502 [2203.12169]. [62]KamLAND-Zencollaboration,Search for the majorana nature of neutrinos in the inverted mass ordering region with KamLAND-Zen,Phys. Rev. Lett.130(2023) 051801 [2203.02139]. [63]LEGENDcollaboration,First results on the s...
arXiv 2022
-
[49]
J. Engel and J. Men´ endez,Status and future of nuclear matrix elements for neutrinoless double-beta decay: a review,Rept. Prog. Phys.80(2017) 046301 [1610.06548]
Pith/arXiv arXiv 2017
-
[50]
M. Chevallier and D. Polarski,Accelerating universes with scaling dark matter,Int. J. Mod. Phys. D10(2001) 213 [gr-qc/0009008]
Pith/arXiv arXiv 2001
-
[51]
Linder,Exploring the expansion history of the universe,Phys
E.V. Linder,Exploring the expansion history of the universe,Phys. Rev. Lett.90(2003) 091301 [astro-ph/0208512]
Pith/arXiv arXiv 2003
-
[52]
M. Lattanzi and M. Gerbino,Status of neutrino properties and future prospects: Cosmological and astrophysical constraints,Front. in Phys.5(2018) 70 [1712.07109]
Pith/arXiv arXiv 2018
- [53]
-
[54]
J. Lesgourgues and S. Pastor,Neutrino mass from cosmology,Adv. High Energy Phys.2012 (2012) 608515 [1212.6154]
Pith/arXiv arXiv 2012
-
[55]
Hergt, W.J
L.T. Hergt, W.J. Handley, M.P. Hobson and A.N. Lasenby,Bayesian evidence for the tensor-to-scalar ratiorand neutrino massesm ν: Effects of uniform vs logarithmic priors, Phys. Rev. D103(2021) 123511
2021
-
[56]
Jimenez, C
R. Jimenez, C. Pe˜ na Garay, F. Simpson and L. Verde,The coherence principle: A falsifiable prior for model selection from the grammar of theories (in preparation),JCAP(2026) [26XX.XXXX]. – 24 –
2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.