⚛️ phenomenology

Cosmological Constraints on Neutrino Masses in a Second-Order CPL Dark Energy Model

This study analyzes cosmological constraints on the sum of neutrino masses across $\Lambda$ CDM, CPL, and a second-order EXP dark energy models using various datasets and hierarchies, finding that the CPL parameterization yields tighter bounds than EXP, frequentist limits are stricter than Bayesian ones, and no statistically significant evidence for nonzero neutrino mass consistent with oscillation lower limits is detected.

Original authors: Shubham Barua, Shantanu Desai

Published 2026-01-26

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Shubham Barua, Shantanu Desai

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to figure out exactly how fast this balloon is inflating and what is pushing it to expand. They also want to know the weight of the tiny, ghostly particles called neutrinos that are zipping through the balloon. These particles are so light and elusive that we can't weigh them directly on a scale; instead, we have to guess their weight by watching how they tug on the fabric of the universe.

This paper is like a team of detectives (the authors) trying to solve two mysteries at once: How heavy are neutrinos? and What is the mysterious "Dark Energy" pushing the universe apart?

Here is the breakdown of their investigation using simple analogies:

1. The Three Suspects (Dark Energy Models)

To understand the universe's expansion, scientists use mathematical "rules" or models. The authors tested three different rulebooks:

The "Steady Hand" (ΛCDM): This is the old, trusted rulebook. It assumes the force pushing the universe apart is constant and unchanging, like a car driving at a perfectly steady speed.
The "Changing Driver" (CPL): This rulebook suggests the force changes over time. It's like a driver who slowly steps on the gas or the brake as the journey goes on.
The "Advanced Driver" (EXP): This is the new, fancy rulebook the authors tested. It's like the "Changing Driver" but with a second gear added. It allows for even more complex changes in how the universe expands, adding a "second-order correction" to the math.

2. The Evidence (The Datasets)

The detectives gathered clues from three different sources:

The Cosmic Microwave Background (CMB): This is the "baby photo" of the universe, showing what it looked like when it was very young.
Baryon Acoustic Oscillations (BAO): Think of these as "fossilized sound waves" frozen in the distribution of galaxies. They act like a cosmic ruler to measure distances.
Supernovae (SNe): These are exploding stars that act as "standard candles." By seeing how bright they look from Earth, scientists can tell how far away they are and how fast the universe is stretching.

The authors combined these clues in different ways (like mixing ingredients in a recipe) to see how the results changed.

3. The Investigation: Weighing the Ghosts

The main goal was to put an upper limit on the total weight of neutrinos. Since we can't weigh them directly, the scientists asked: "What is the heaviest the neutrinos could possibly be without breaking the laws of physics we see in the data?"

They tested four different "scenarios" for how the three types of neutrinos might be weighted:

Scenario A: One heavy neutrino, two ghosts (massless).
Scenario B: All three are equally heavy (degenerate).
Scenario C: Normal Hierarchy (light, medium, heavy).
Scenario D: Inverted Hierarchy (heavy, medium, light).

They also used two different ways of doing the math:

Bayesian: Like a detective who starts with a strong hunch (a "prior") and updates it as new evidence comes in.
Frequentist: Like a detective who looks strictly at the data without any pre-existing hunches, asking, "If the neutrinos were this heavy, how likely is it that we would see this data?"

4. The Big Discoveries

Here is what the authors found, translated into everyday terms:

The "Simple" Rulebook is the Strictest: When they used the "Steady Hand" model (ΛCDM), they got the tightest, most restrictive limits on neutrino mass. It's like a strict judge who says, "You can't be heavier than this."
The "Fancy" Rulebooks are More Lenient: When they used the "Changing Driver" (CPL) or the "Advanced Driver" (EXP) models, the limits on neutrino weight became much looser (about 10-65% higher). It's as if the judge said, "Well, if the universe is behaving in this complex way, the neutrinos could be a bit heavier."
The "Advanced" Driver is the Most Lenient: The new EXP model gave slightly looser limits than the CPL model. Adding that extra "second gear" to the math made it even harder to pin down the exact weight of the neutrinos.
More Data = Tighter Limits (Usually): When they added the Supernova data (the exploding stars) to the mix, the limits generally got tighter for the complex models. It's like adding more witnesses to a trial; the story becomes clearer. However, for the "Steady Hand" model, adding this data actually made the limits slightly looser.
The "Hunch" Matters: The results changed depending on whether they used the "Bayesian" (hunch-based) or "Frequentist" (data-only) math. The Frequentist approach usually gave tighter (stricter) limits.
No "Smoking Gun": Despite all this, the authors found no statistically significant evidence that neutrinos definitely have a non-zero mass that fits with what we know from lab experiments. In other words, the data doesn't scream "Neutrinos are heavy!" It just says, "They could be this heavy, but they could also be lighter."

5. The Conclusion

The paper concludes that how we choose to describe the universe's expansion (the Dark Energy model) dramatically changes our estimate of how heavy neutrinos are.

If you assume the universe expands in a simple, steady way, you get a very strict weight limit for neutrinos. If you assume the expansion is complex and changing, that weight limit goes up.

The authors emphasize that the "detection" of neutrino mass isn't just about the data; it's about the mathematical rules we choose to interpret that data. They found that while some models hint at a positive mass, once you apply the strict physical limits we know from lab experiments (that neutrinos must be at least a tiny bit heavy), the evidence for a specific heavy mass disappears.

In short: The universe is a complex puzzle. Depending on which piece of the puzzle (the Dark Energy model) you hold up first, the picture of the neutrino's weight changes. The authors didn't find a definitive new weight, but they proved that our assumptions about the universe's expansion are the most critical factor in guessing that weight.

Technical Summary: Cosmological Constraints on Neutrino Masses in a Second-Order CPL Dark Energy Model

Problem Statement
While neutrinos are fundamental components of the Standard Model, their absolute masses and true nature remain unknown. Cosmological observables are primarily sensitive to the sum of the three neutrino masses ( $M_{\text{tot}}$ ), yet current experiments lack the sensitivity to distinguish between different mass splittings (hierarchies) with high significance. Recent results from the Dark Energy Spectroscopic Instrument (DESI), when combined with Cosmic Microwave Background (CMB) and Supernovae (SNe) data, indicate a strong preference for dynamical dark energy (DE) over the standard $\Lambda$ CDM model. This study investigates how the choice of DE parameterization—specifically extending the Chevallier-Polarski-Linder (CPL) model to include second-order corrections (EXP)—affects the inferred upper bounds on $M_{\text{tot}}$ . Furthermore, it examines the interplay between different neutrino mass hierarchies, dataset combinations, and statistical frameworks (Bayesian vs. frequentist) in determining these constraints.

Methodology
The authors analyze constraints on $M_{\text{tot}}$ across three DE models:

$\Lambda$ CDM: Constant equation of state (EOS) $w = -1$ .
CPL: Linear evolution in redshift, $w(z) = w_0 + w_a \frac{z}{1+z}$ .
EXP: An exponential parameterization introducing second-order corrections to CPL, $w(z) \approx w_0 + w_a [\frac{z}{1+z} + \frac{1}{2}(\frac{z}{1+z})^2]$ , without increasing parameter space dimensionality.

Four neutrino mass hierarchies are considered:

1M: One massive, two massless neutrinos.
DH: Degenerate hierarchy ( $m_1=m_2=m_3$ ).
NH: Normal hierarchy ( $m_1 < m_2 < m_3$ ).
IH: Inverted hierarchy ( $m_3 < m_1 < m_2$ ).

The analysis utilizes three dataset combinations:

CMB (Planck PR3) + BAO (DESI DR2).
CMB + BAO + PantheonPlus (SNe Ia).
CMB + BAO + DESY5 (SNe Ia).

The study employs both Bayesian inference (using MontePython and GetDist with MCMC sampling) and frequentist analysis (using profile likelihoods via the pinc code with simulated annealing). Physical lower limits derived from neutrino oscillation experiments are imposed: $M_{\text{tot}} \geq 0.059$ eV for NH and $M_{\text{tot}} \geq 0.11$ eV for IH. For 1M and DH, the lower limit is set to 0 eV. The frequentist approach utilizes the Feldman-Cousins construction to handle physical boundaries.

Key Contributions

Extension of DE Parameterization: The paper extends previous analyses (specifically Ref. [89]) by applying both Bayesian and frequentist frameworks to the EXP parameterization, which includes second-order redshift corrections to the CPL form.
Systematic Comparison: It provides a comprehensive comparison of neutrino mass bounds across different DE models, hierarchies, and dataset combinations, explicitly quantifying the impact of adding SNe data (PantheonPlus vs. DESY5) to CMB+BAO.
Statistical Framework Contrast: The work highlights the differences between Bayesian and frequentist upper bounds, particularly regarding the sensitivity to prior choices and physical lower limits.

Results

Impact of DE Parameterization: The CPL model yields tighter upper bounds on $M_{\text{tot}}$ (approximately $\lesssim 10\%$ tighter) compared to the EXP model. The $\Lambda$ CDM model provides the tightest constraints overall, while extended DE models (CPL and EXP) result in relaxed bounds due to the increased parameter space and degeneracy between $M_{\text{tot}}$ and the DE EOS.
Impact of Datasets: Adding SNe data (PantheonPlus or DESY5) to CMB+BAO generally tightens the constraints for CPL and EXP models. However, for $\Lambda$ CDM, the inclusion of SNe data can relax the upper bounds. DESY5 generally yields more relaxed bounds compared to PantheonPlus.
Impact of Hierarchies: Neutrino mass constraints show only mild sensitivity to the assumed hierarchy. In the Bayesian framework, bounds for 1M and DH are similar, while NH and IH bounds are relaxed (consistent with the imposed priors). In the frequentist framework, bounds tighten when moving from 1M to NH but can become more stringent for IH depending on the dataset and model.
Statistical Frameworks: Frequentist upper bounds are consistently tighter than Bayesian bounds across all models and hierarchies. For $\Lambda$ CDM, frequentist limits are significantly tighter (20–40%) than Bayesian ones, whereas for CPL and EXP, the difference is smaller ( $\lesssim 10\%$ ).
Evidence for Non-Zero Mass: For the datasets and parameterizations considered, the study finds no statistically significant evidence for a non-zero neutrino mass consistent with oscillation lower limits.
- In $\Lambda$ CDM, frequentist profile likelihood minima often lie in the unphysical region ( $M_{\text{tot}} < 0$ ) or below the oscillation lower limits.
- CPL and EXP models show indications of positive neutrino mass (minima in the physical region) when using CMB+BAO alone. However, upon imposing the physical lower limits from oscillation experiments (0.059 eV for NH, 0.11 eV for IH), these apparent detections disappear, and the bounds become consistent with the lower limits.
- Specifically, $\Lambda$ CDM bounds often show tension with the IH lower limit (0.11 eV), whereas CPL and EXP, due to their relaxed constraints, remain consistent with all lower limits (0 eV, 0.059 eV, and 0.11 eV).

Significance and Conclusions
The paper concludes that the inferred cosmological neutrino mass bounds are highly dependent on the specific DE parameterization, the datasets used, and the statistical framework (priors vs. physical limits) rather than solely on the data itself. The introduction of second-order corrections in the EXP parameterization leads to slightly more relaxed bounds compared to CPL, demonstrating that even small increases in the flexibility of the DE EOS can alter neutrino mass constraints.

The authors reaffirm that while cosmological data is sensitive to the total neutrino mass, the detection of a non-zero mass is not robust across all model choices. The "detection" of positive neutrino mass is contingent upon the choice of priors, the mass hierarchy, and the DE parameterization. The study suggests that future improvements in observational precision will be necessary to reconcile terrestrial oscillation limits with cosmologically inferred upper bounds. The work serves as a cautionary note that constraints on neutrino masses in extended DE models must be interpreted carefully, acknowledging the degeneracies between $M_{\text{tot}}$ and the dark energy equation of state.

1. The Three Suspects (Dark Energy Models)

2. The Evidence (The Datasets)

3. The Investigation: Weighing the Ghosts

4. The Big Discoveries

5. The Conclusion

More like this