Towards a complete scheme of cosmological neutrino self-interactions: Collision term for a wide range of mediator masses

This paper presents a novel framework for calculating the neutrino-neutrino collision term in the Boltzmann hierarchy that incorporates both neutrino and mediator masses as free parameters, enabling a smooth transition between light and heavy mediator regimes and providing a robust tool for testing neutrino self-interaction scenarios across a wide range of cosmological conditions.

The Gist

Imagine the universe as a giant, expanding ballroom. Inside this ballroom, there are billions of tiny, ghostly dancers called neutrinos. These particles are everywhere, but they are notoriously shy; they rarely bump into each other or anything else. In the standard story of cosmology, they just drift through the crowd, barely interacting.

However, recent observations suggest something might be wrong with this "drifting" story. There are some cosmic tensions—like a mismatch between how fast the universe is expanding now versus how fast it was expanding in the past. Scientists suspect that maybe these neutrino dancers aren't just drifting; maybe they are actually holding hands, dancing in pairs, or bumping into each other more than we thought. This is called Neutrino Self-Interaction (NSI).

This paper is like a massive, brand-new instruction manual for how to calculate exactly what happens when these neutrinos bump into each other. Here is the breakdown of what the authors did, using simple analogies:

1. The Problem: The "Too Heavy" and "Too Light" Tools

Previously, scientists trying to model these neutrino bumps had to use two very different, simplified tools depending on the size of the "messenger" particle (called a mediator) that carries the force between them.

• The Heavy Mediator: Imagine the mediator is a giant boulder. If the neutrinos are far apart, they can't feel the boulder. Scientists used a simplified rule: "If the boulder is heavy, ignore the details and just assume they don't talk." This worked well for some scenarios but broke down when the universe was hot and the neutrinos were energetic.
• The Light Mediator: Imagine the mediator is a feather. It floats everywhere. Scientists used a different simplified rule: "If the feather is light, assume the neutrinos are constantly bumping into it."

The problem was that the universe cools down over time. A mediator might act like a heavy boulder when the universe is hot, but turn into a light feather as the universe cools. The old tools couldn't handle this smooth transition. They were like trying to use a hammer to fix a watch, and then switching to a screwdriver halfway through, without a manual on how to switch.

2. The Solution: The "Universal Translator"

The authors of this paper built a universal calculator (a new mathematical framework) that works for any size of mediator, from the heaviest boulder to the lightest feather, and everything in between.

• No More Approximations: Instead of guessing, they calculated the exact "collision term." In physics, a collision term is just a fancy way of saying "the math that describes how often and how hard particles bump into each other."
• The Smooth Transition: Their new math shows that as the universe cools down, the behavior of the neutrinos changes smoothly. It's like a dimmer switch rather than a light switch. You don't have to jump from "Heavy Mode" to "Light Mode"; the math handles the whole gradient.

3. The Details: Who is Dancing?

The authors also made sure to account for the specific "identity" of the dancers:

• Dirac vs. Majorana: Neutrinos might be their own antiparticles (Majorana) or distinct from them (Dirac). It's like asking if the dancers are wearing red shoes or blue shoes. The authors' new manual calculates the collision rates for both possibilities, showing that while the "shoe color" matters in some specific, high-energy scenarios, it often doesn't change the big picture for the heavy mediators.
• Mass Matters: They also included the actual weight of the neutrinos. In the past, scientists often assumed neutrinos were weightless ghosts. This paper says, "Let's include their actual tiny weight," which turns out to be important when the universe gets cold enough.

4. Why Does This Matter?

Think of the universe's history as a movie. To understand the plot, you need to know how the characters interact.

• The "Cosmic Tension": We have two different ways of measuring the universe's expansion, and they disagree.
• The Fix: If neutrinos are interacting (bumping into each other), it changes the "soundtrack" of the early universe. This could explain why our measurements disagree.
• The New Tool: Before this paper, if scientists wanted to test a specific theory about neutrino interactions, they had to hope their theory fit into one of the old, simplified boxes. Now, they can plug in any theory, any mediator mass, and any neutrino mass, and get a precise answer.

The Bottom Line

This paper is a master key for cosmologists. It replaces a set of broken, specific keys with one master key that opens every door in the "Neutrino Interaction" house.

If future telescopes and experiments detect a signal that neutrinos are indeed interacting, this paper provides the exact mathematical engine needed to interpret that signal, tell us what kind of "messenger" particle is involved, and perhaps finally solve the mystery of why the universe is expanding the way it is. It turns a blurry, guesswork-heavy picture into a high-definition, precise map of the neutrino dance floor.

Technical Summary

1. Problem Statement

Standard cosmological models face several tensions (e.g., the $H_0$ tension, $S_8$ tension) that non-standard neutrino interactions (NSI) could potentially resolve. Neutrino self-interactions (NSI) mediated by a new scalar particle have been proposed to address these issues. However, existing theoretical frameworks for modeling NSI in the Boltzmann hierarchy rely on specific approximations that limit their applicability:

• Heavy Mediator Limit (HML): Assumes the mediator mass $m_\phi$ is much larger than the momentum transfer ($m_\phi \gg T_\nu$), often neglecting neutrino masses and treating the interaction as a contact four-fermion interaction.
• Light Mediator Limit: Assumes $m_\phi \ll T_\nu$, often utilizing the Relaxation Time Approximation (RTA) which simplifies the collision integral but may lose precision.
• Resonant Regime: Occurs when $m_\phi \sim T_\nu$, where resonant production of the mediator occurs. This regime is difficult to model and often treated separately or ignored.
• Neutrino Masses: Most previous calculations assume massless neutrinos, ignoring the effects of non-zero neutrino masses ($m_\nu$) which are now known to be non-zero and relevant for precision cosmology.

The authors identify a gap in the literature: a unified, general framework that computes the neutrino-neutrino collision term without restricting the mediator mass or neutrino mass, distinguishing between Dirac and Majorana neutrinos, and valid across all cosmological epochs.

2. Methodology

The authors develop a novel theoretical framework based on the relativistic Boltzmann hierarchy in a perturbed FLRW universe.

• Formalism: They start from the general relativistic Boltzmann equation for 2$\to$2 elastic scattering processes ($\nu_i + \nu_j \to \nu_k + \nu_l$) mediated by a scalar field $\phi$.
• Interaction Lagrangian: They consider a scalar interaction $L_{int} = \sum g_{ij} \phi \nu_i \nu_j$, covering both Dirac-like and Majorana neutrinos.
• Collision Term Derivation:
  • They derive the collision term $C_\ell$ in Fourier space, expanding the distribution function perturbation $\Theta$ into Legendre multipoles ($\ell$).
  • They explicitly calculate the squared scattering amplitudes $|\mathcal{M}|^2$ for both Dirac and Majorana cases, accounting for $s$, $t$, and $u$ channels and their interference terms.
  • They reduce the 4-body phase space integrals into a set of reduced angular integrals ($J_{n\ell}$) involving Mandelstam variables ($s, t, u$) and kinematic constraints.
• Handling Singularities: To deal with collinear infrared singularities and kinematic singularities associated with on-shell mediator production (resonances), they introduce $h$-functions. These functions regularize the integrals by setting the contribution to zero in specific kinematic regions where the propagator would diverge, effectively isolating the resonant contribution.
• Numerical Implementation: The resulting collision terms are expressed as integrals over dimensionless energy variables ($z = E/T_\nu$) and angular variables. These integrals are evaluated numerically using the VEGAS algorithm.

3. Key Contributions

1. Unified Framework: The paper provides the first complete calculation of the NSI collision term where the mediator mass ($m_\phi$) and neutrino mass ($m_\nu$) are treated as free parameters, rather than fixed limits.
2. Dirac vs. Majorana Distinction: The authors explicitly derive and compare the collision terms for Dirac-like and Majorana neutrinos. They show that while the heavy mediator limit makes them indistinguishable, the full calculation reveals differences in available scattering channels (e.g., $s$-channel contributions for Majorana vs. Dirac).
3. Thermal Corrections for Light Mediators: They provide an exact calculation for the light mediator regime without relying on the Relaxation Time Approximation (RTA). This allows for the determination of the previously unknown thermal correction factor $\xi_\ell$ and its dependence on the multipole order $\ell$.
4. Transition Analysis: The framework demonstrates a smooth transition between the heavy mediator limit, the resonant regime, and the light mediator limit as the universe cools (temperature $T_\nu$ changes relative to $m_\phi$).

4. Key Results

• Heavy Mediator Limit (HML) Validity:
  • The authors confirm that the standard HML approximation (contact interaction) is valid for $m_\phi \gtrsim 10^2$ eV and temperatures $T_\nu \lesssim 100$ eV.
  • They show that for $m_\nu \lesssim 0.06$ eV, neutrino mass effects are negligible in the HML regime, justifying previous massless approximations in that specific limit.
  • However, they demonstrate that at higher temperatures ($T_\nu \gtrsim m_\phi$), the interaction rate saturates and deviates from the simple power-law scaling ($\Gamma \propto T^5$) assumed in the HML, requiring the full numerical treatment.
• Neutrino Mass Effects:
  • At low temperatures (late universe), neutrino masses suppress the collision term.
  • For $m_\nu \approx 0.1$ eV, deviations from the massless case become visible in the quadrupole ($\ell=2$) moment, though they remain small for standard mass limits ($m_\nu < 0.05$ eV).
• Light Mediator & RTA:
  • For $m_\phi \lesssim 10^{-3}$ eV, the interaction rate becomes constant (independent of temperature) in the dimensionless formulation.
  • They calculate the coefficients $\beta_\ell$ and the thermal correction factors $\xi_\ell$ (Table 1). They find that $\xi_\ell$ saturates for $\ell \geq 20$ and differs by a factor of $\approx 2$ between Majorana and Dirac cases due to $s$-channel suppression in the Dirac case.
• Resonant Regime:
  • The framework identifies the resonant regime roughly at $\mu_\phi = m_\phi/T_\nu \sim 2\times 10^{-3} - 50$.
  • They note that modeling the resonant regime fully is challenging due to the need for one-loop corrections to describe the amplitude accurately, which is beyond the scope of this tree-level calculation.
• Dirac vs. Majorana Differences:
  • In the resonant and light mediator regimes, the scattering channels differ significantly. For example, Dirac-like neutrinos do not exhibit resonance production in $\nu-\nu$ scattering (same mass eigenstate) but do in $\nu-\bar{\nu}$ scattering, whereas Majorana neutrinos allow $s$-channel resonances in both.

5. Significance and Future Outlook

• Cosmological Data Analysis: This work provides the necessary tools for future cosmological data analyses (e.g., CMB-S4, Euclid, DESI) to test NSI models without relying on restrictive approximations. It allows for a direct exploration of the parameter space ($g, m_\phi, m_\nu$) rather than just effective couplings.
• Precision Cosmology: By including neutrino masses and distinguishing between Dirac/Majorana nature, the framework improves the precision of constraints on NSI, potentially helping to resolve cosmological tensions or rule out specific mediator models.
• Beyond NSI: The authors note that the formalism can be adapted to other scenarios, such as warm dark matter self-interactions or other light self-interacting particles.
• Future Work: The authors plan to implement this formalism into cosmological Boltzmann solvers (like CLASS or CAMB) to generate predictions for observables and constrain the mediator mass and coupling strength using current and future datasets. They also intend to address the resonant regime with higher-order corrections in subsequent studies.

In summary, this paper bridges the gap between simplified analytical approximations and the complex reality of neutrino self-interactions, offering a robust, general-purpose tool for the next generation of cosmological investigations.