⚛️ phenomenology

Cosmological Constraints on Temperature-Dependent Interaction between Dark Matter and Neutrinos

This study investigates temperature-dependent dark matter-neutrino interactions induced by a dimension-six operator, deriving a full Boltzmann hierarchy that reveals dark acoustic oscillations and establishes significantly tighter cosmological constraints on the interaction strength compared to temperature-independent models, while highlighting the critical impact of realistic neutrino mass ordering on these limits.

Original authors: Ren-Peng Zhou, Da Huang

Published 2026-02-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ren-Peng Zhou, Da Huang

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, bustling dance floor. For decades, physicists have believed that the two main groups of dancers—Dark Matter (the invisible, heavy crowd) and Neutrinos (the ghostly, fast-moving particles)—never actually touched. They just danced around each other, influenced only by gravity, like strangers in a crowded room who never make eye contact.

This paper, written by Ren-Peng Zhou and Da Huang, asks a simple but profound question: What if they do bump into each other?

Specifically, the authors investigate a scenario where these two groups interact, but with a twist: the strength of their interaction depends on how "hot" the universe is.

Here is the breakdown of their findings using everyday analogies:

1. The "Hot" Interaction

In the early universe, things were incredibly hot. The authors propose a model where Dark Matter and Neutrinos interact via a specific force (like a magnetic field that only turns on when things are hot).

The Analogy: Imagine a dance floor where the music is so loud and hot that the dancers start sweating. As the temperature rises, the dancers get more energetic and start bumping into each other more often.
The Twist: In this model, the hotter the universe, the stronger the "bumping." As the universe cools down (like a party winding down), the dancers stop bumping into each other and go back to dancing alone. This is different from previous theories that assumed they bumped into each other with the same intensity regardless of the temperature.

2. The "Dark Acoustic Oscillation" (The Ripple Effect)

When Dark Matter and Neutrinos start interacting, they stop acting like two separate groups and start acting like a single, thick fluid.

The Analogy: Think of a crowd of people walking through a hallway. If they walk alone, they move smoothly. But if they start holding hands and pushing against each other, they create a "traffic jam" wave. They push forward, get squeezed, and then bounce back.
The Result: This creates "Dark Acoustic Oscillations" (DAO). It's like a ripple moving through a pond, but the pond is made of invisible Dark Matter and ghostly Neutrinos. These ripples leave a specific fingerprint on the Cosmic Microwave Background (CMB)—the "baby picture" of the universe—and on how galaxies are distributed today.

3. The "Ghost" Problem (Neutrino Mass)

A major hurdle in previous studies was treating neutrinos as massless ghosts. But we now know neutrinos have a tiny bit of mass.

The Analogy: Imagine the neutrinos are like a mix of lightweight ping-pong balls and slightly heavier tennis balls. In the early, hot universe, both fly around fast. But as the universe cools, the "tennis balls" (heavier neutrinos) slow down and stop interacting with Dark Matter, while the "ping-pong balls" keep going.
The Insight: The authors realized that if you ignore the fact that some neutrinos slow down and "opt out" of the interaction, your calculations are wrong. By accounting for this "mass effect," they found that the interaction leaves a much clearer, more distinct signature than previously thought.

4. The Detective Work (The Constraints)

The authors used the latest data from powerful telescopes (like Planck, DESI, and ACT) to look for these fingerprints. They treated the universe like a crime scene, looking for evidence of this "bumping."

The Verdict: They didn't find a smoking gun (a confirmed interaction), but they set a much stricter rule on how much bumping could be happening.
The Improvement: Previous studies said, "Okay, maybe they bump once every billion years." This new study says, "No, if they bump at all, it has to be less than once every trillion years."
Why so strict? Because the interaction gets stronger when the universe is hot. The early universe was a furnace, so even a tiny interaction would have left a massive scar on the cosmic data. The fact that we don't see those scars means the interaction must be incredibly weak today.

5. The "Prior" Surprise

There was one interesting side note. In statistics, how you start your guess (your "prior") can change the result.

The Analogy: If you ask someone, "How likely is it that aliens exist?" and they start with the assumption "Aliens are common," they might find evidence for aliens in a weird cloud. If they start with "Aliens are rare," they won't.
The Finding: When the authors used a specific statistical starting point (a "logarithmic flat prior"), the data slightly hinted that the interaction might not be zero after all. It's a tiny whisper, not a shout, but it suggests that future, more precise telescopes might finally catch Dark Matter and Neutrinos in the act.

Summary

This paper is like upgrading the security system of the universe. By realizing that the "temperature" of the universe changes how Dark Matter and Neutrinos interact, and by accounting for the fact that some neutrinos are "heavier" than others, the authors have tightened the screws on our understanding.

They haven't found the interaction yet, but they've proven that if it exists, it's incredibly shy. And, just maybe, there's a tiny chance it's whispering a secret to us that we're just starting to hear.

1. Problem Statement

The standard $\Lambda$ CDM cosmological model assumes Dark Matter (DM) is cold and collisionless. However, the nature of DM remains unknown, and non-gravitational interactions between DM and Standard Model particles (specifically neutrinos) are a viable theoretical possibility. Previous studies on DM-neutrino interactions often relied on simplifying assumptions:

Temperature Independence: Assuming a constant scattering cross-section or a simple $T^2$ dependence valid only for massless neutrinos.
Neutrino Mass Degeneracy: Assuming all three neutrino flavors have identical masses, which contradicts current oscillation data favoring a Normal Ordering (two massless/light flavors and one massive flavor).
Incomplete Dynamics: Often neglecting the transition of neutrinos from relativistic to non-relativistic regimes, which significantly alters interaction dynamics in the late universe.

The authors aim to address these gaps by studying temperature-dependent DM-neutrino interactions induced by a dimension-six operator, incorporating realistic neutrino mass ordering, and deriving the full Boltzmann hierarchy to accurately model the cosmological impact.

2. Methodology

A. Theoretical Framework

Interaction Model: The authors propose an effective field theory (EFT) interaction described by a dimension-six operator:
$\mathcal{L}_{\chi\nu} = \frac{1}{\Lambda^2} \bar{\chi}\gamma^\mu\chi \bar{\nu}_L\gamma_\mu\nu_L$
where $\chi$ is a Dirac fermion DM particle and $\nu_L$ are active left-handed neutrinos.
Cross-Section Behavior: Unlike previous constant-cross-section models, this operator yields a scattering cross-section $\sigma_{\chi-\nu}$ $σ_{χ - ν}$ that scales with the square of the neutrino energy ( $E_\nu^2$ $E_{ν}^{2}$ ).
- In the high-temperature limit ( $T \gg m_\nu$ ), $\sigma \propto T^2$ .
- In the low-temperature limit ( $T \ll m_\nu$ ), the massive neutrinos become non-relativistic, causing the cross-section to deviate significantly from the $T^2$ scaling, effectively "turning off" the interaction for the massive flavor.

B. Boltzmann Hierarchy Derivation

Full Hierarchy: The authors derive the complete Boltzmann hierarchy for massive neutrinos in the presence of this temperature-dependent collision term. They explicitly calculate the collision terms for the neutrino distribution function and the DM Euler equation.
Fluid Approximation: A novel method is introduced to derive the fluid approximation for modes deep inside the horizon. Instead of fitting simulation data (as done in previous works like Ref. [47]), they analytically average the perturbation moments (density contrast, velocity divergence, shear stress) over momentum space. This allows the transport coefficients to be calculated directly by the Boltzmann solver CLASS without numerical simulations.
Implementation: The modified Boltzmann equations are implemented in the CLASS code, coupled with the Cobaya MCMC sampler for parameter estimation.

C. Data and Analysis

Datasets: The analysis utilizes the latest cosmological data:
- Planck 2020 (PR4): High- $\ell$ TT, TE, EE likelihoods combined with low- $\ell$ data.
- CMB Lensing: From Planck PR4.
- BAO: Baryon Acoustic Oscillation measurements from DESI DR2.
- ACT: Atacama Cosmology Telescope DR6 high- $\ell$ data.
Scenarios: Three main scenarios are tested:
1. Benchmark: Normal neutrino mass ordering with fixed total mass $\sum m_\nu = 0.06$ eV.
2. Degenerate Masses: Three degenerate neutrinos (for comparison with literature).
3. Varying Mass: Normal ordering with $\sum m_\nu$ as a free parameter.
Priors: Both linear flat priors and logarithmic flat priors are tested for the interaction parameter $u^0_{\chi-\nu}$ .

3. Key Contributions

Realistic Neutrino Mass Treatment: This is the first study to rigorously apply temperature-dependent DM-neutrino interactions to the Normal Ordering scenario (2 massless + 1 massive flavor), moving beyond the degenerate mass assumption.
Analytical Fluid Approximation: The derivation of fluid equations via analytical momentum averaging provides a more robust and computationally efficient method for handling massive neutrinos in interacting scenarios compared to previous fitting-based approaches.
Enhanced Constraints: The study demonstrates that the temperature dependence of the cross-section leads to significantly stronger constraints on the interaction strength compared to temperature-independent models.
Prior Dependence Analysis: The authors explicitly investigate how the choice of prior (linear vs. logarithmic) affects the constraints, finding hints of non-zero interaction under specific conditions.

4. Results

A. Constraints on Interaction Strength ( $u^0_{\chi-\nu}$ )

The interaction parameter $u^0_{\chi-\nu}$ is defined relative to the Thomson cross-section and DM mass.

Normal Ordering (Benchmark): Using Planck + Lensing + BAO, the 95% CL upper limit is $u^0_{\chi-\nu} \lesssim 2.53 \times 10^{-13}$ .
Comparison: This constraint is ~9 orders of magnitude stronger than constraints derived for temperature-independent interactions in the literature ( $u_{\chi-\nu} \lesssim 10^{-4}$ $u_{χ - ν} ≲ 1 0^{- 4}$ ).
- Reasoning: The cross-section scales quadratically with temperature in the early universe ( $T \gg m_\nu$ ), leading to much stronger scattering effects during the CMB formation era.
Degenerate Masses: Constraints are weaker ( $u^0_{\chi-\nu} \lesssim 8 \times 10^{-11}$ ), highlighting the importance of realistic mass ordering.
Varying Mass: Constraints remain similar to the fixed mass case ( $u^0_{\chi-\nu} \lesssim 2.11 \times 10^{-13}$ ), with a tight bound on total neutrino mass ( $\sum m_\nu \lesssim 0.102$ eV).

B. Cosmological Signatures

Dark Acoustic Oscillations (DAO): The interaction induces pressure in the DM-neutrino fluid, creating damped oscillations in the matter power spectrum at small scales.
CMB Effects:
- Gravitational Boost: DM contributes to fast-evolving modes, enhancing CMB peak amplitudes.
- Drag Effect: Tightly coupled DM and neutrinos reduce the sound speed, shifting acoustic peaks to higher multipoles.
$H_0$ and $\sigma_8$ Tensions:
- $\sigma_8$ : The interaction suppresses small-scale matter fluctuations, reducing $\sigma_8$ . The model predicts $S_8 \approx 0.81$ , which aligns well with recent KiDS-Legacy weak lensing results, potentially alleviating the $\sigma_8$ tension.
- $H_0$ : The model does not resolve the $H_0$ tension; the derived $H_0$ values remain consistent with Planck ( $\sim 67-68$ km/s/Mpc) and lower than local distance ladder measurements.

C. Prior Sensitivity

When employing a logarithmic flat prior for $u^0_{\chi-\nu}$ and combining Planck, DESI, and ACT data, the analysis shows a hint of non-zero interaction at the 95% confidence level ( $\log_{10} u^0_{\chi-\nu} \approx -14.6$ ). This suggests that future high-precision data could potentially detect this interaction.

5. Significance

Theoretical Rigor: By incorporating the full momentum dependence of massive neutrinos and realistic mass ordering, this work provides the most accurate cosmological constraints on DM-neutrino interactions to date.
Experimental Guidance: The dramatic improvement in constraints (9 orders of magnitude) sets a new benchmark for particle physics models proposing such interactions, effectively ruling out large regions of parameter space.
Cosmological Tensions: While not solving the $H_0$ crisis, the model offers a viable mechanism to address the $\sigma_8$ tension without conflicting with other observations, linking neutrino physics to large-scale structure formation.
Methodological Advance: The analytical fluid approximation derived here improves the efficiency and reliability of cosmological simulations involving interacting massive neutrinos.

In conclusion, the paper establishes that temperature-dependent DM-neutrino interactions are tightly constrained by current data, but the specific choice of neutrino mass ordering and prior distributions plays a critical role in interpreting these limits. The detection of such interactions remains a tantalizing possibility for future high-precision cosmological surveys.