Sterile Neutrinos as a Dynamical Cosmological Fluid: Implications for the Expansion History and Matter-Radiation Equality

This paper presents an analytic framework treating finite-mass sterile neutrinos as a dynamical cosmological fluid with a time-dependent equation of state, demonstrating how their suppressed thermalization and evolving pressure impact the expansion history and shift the matter-radiation equality epoch beyond the standard constant $\Delta N_{\text{eff}}$ approximation.

The Gist

The Big Idea: Ghost Particles That Change Their Mind

Imagine the early Universe as a massive, crowded dance floor. There are standard dancers (the "active" neutrinos we know) and a few shy, invisible dancers hiding in the shadows (the "sterile" neutrinos).

For a long time, scientists thought these invisible dancers were either fully present (dancing just as hard as everyone else) or completely absent. They measured their presence by counting how many "extra dancers" were on the floor, a number called $\Delta N_{eff}$.

This paper argues that reality is more complicated. These sterile neutrinos are like a crowd of people who start the night dancing wildly (acting like light/radiation) but slowly get tired, sit down, and start walking around (acting like heavy matter) as the night goes on.

The authors say: "Stop treating them as a static number. Treat them as a living, changing fluid that evolves over time."

---

The Story in Three Acts

The paper describes how these particles change the speed at which the Universe expands (the "Hubble expansion"). They identify three distinct phases, like three acts in a play:

Act 1: The High-Energy Party (The Relativistic Plateau)

• What's happening: The Universe is very young and hot. The sterile neutrinos are zipping around at nearly the speed of light.
• The Analogy: Imagine a room full of hyperactive toddlers running everywhere. They are light, fast, and take up space like "radiation."
• The Effect: During this phase, they act exactly like standard light. They speed up the expansion of the Universe just a tiny bit, similar to adding more light bulbs to a room.

Act 2: The Tired Transition (The Switch)

• What's happening: The Universe starts to cool down. The sterile neutrinos lose their energy. They can no longer run; they start to slow down.
• The Analogy: The toddlers get tired. They stop running and start walking. They are no longer "light"; they are becoming "heavy."
• The Effect: This is the most interesting part. Standard physics models assume particles stay either "light" or "heavy." But here, the particles are in a gray zone. They are changing their personality. This creates a weird, time-dependent wobble in the expansion rate that old models can't explain.

Act 3: The Heavy Walkers (The Matter-Like Growth)

• What's happening: The Universe is now cool. The sterile neutrinos have slowed down completely. They are now heavy, slow-moving particles.
• The Analogy: The toddlers have grown into adults. They are now heavy, slow-moving people sitting on the dance floor. They act like "matter" (like dust or dark matter), not light.
• The Effect: Because they are now heavy, they start to pull on the Universe differently. Instead of just speeding up the expansion like light, they start to act like gravity, helping to slow things down or change the timing of when the Universe switches from being light-dominated to matter-dominated.

---

The "Big Switch" Moment: Matter-Radiation Equality

In cosmology, there is a very important moment called Matter-Radiation Equality. This is the exact moment in history when the weight of the "stuff" (matter) in the Universe finally overtook the weight of the "light" (radiation).

• Old View: Scientists thought sterile neutrinos were just extra light. If they were there, they would delay this switch, making the Universe stay "light-dominated" longer.
• New View (This Paper): For heavy sterile neutrinos (like those with a mass in the GeV range), by the time this "Big Switch" happens, they have already stopped running. They are already heavy walkers.
• The Result: They don't delay the switch; they advance it. They act like extra heavy furniture in the room, making the "matter" side win the race sooner than expected.

The Detective Work: Catching the Culprit

The authors used this new understanding to look at data from the Planck satellite, which maps the Cosmic Microwave Background (the afterglow of the Big Bang).

• The Clue: The satellite measured exactly when the "Big Switch" happened. It was very precise.
• The Deduction: If there were too many of these "heavy walkers" (sterile neutrinos), the switch would have happened at the wrong time.
• The Verdict: The paper calculates that the amount of these sterile neutrinos at that moment must be very small—less than about 1% of the total energy in the Universe.

Why This Matters

1. It's a Better Model: Instead of using a blunt instrument (a constant number) to describe these particles, this paper uses a scalpel (a changing equation) that fits the physics better.
2. It Connects Micro to Macro: It connects the tiny, invisible world of particle physics (how neutrinos mix and oscillate) to the huge, visible world of the expanding Universe.
3. New Rules: It tells us that if we find these particles, they won't behave like light in the early Universe; they will behave like heavy matter, which changes how we look for them in the future.

Summary in One Sentence

This paper argues that sterile neutrinos are not just static "extra light" in the Universe, but a dynamic fluid that slowly transforms from fast runners into heavy walkers, and this transformation leaves a specific fingerprint on the history of the Universe that we can now use to set strict limits on how many of them exist.

Technical Summary

1. Problem Statement

The standard cosmological parameterization for additional relativistic species, $\Delta N_{\text{eff}}$, assumes a constant radiation-like equation of state ($w = 1/3$). This approximation fails for sterile neutrinos with finite mass ($m_s$) that are only partially thermalized in the early Universe.

• The Gap: Existing literature often treats sterile neutrinos either as fully thermalized (contributing to $\Delta N_{\text{eff}}$) or as specific mass windows (e.g., keV warm dark matter). There is a lack of a self-consistent analytic framework for the intermediate regime where sterile neutrinos have GeV-scale masses and incomplete thermalization.
• The Issue: In this regime, the sterile population transitions dynamically from relativistic to non-relativistic behavior as the Universe cools. Consequently, their equation of state ($w_s$) evolves continuously, and their contribution to the energy density cannot be captured by a static shift in relativistic degrees of freedom. Current Boltzmann codes (e.g., Class, Camb) handle this numerically but lack transparent analytic insight into the background dynamics.

2. Methodology

The authors develop an analytic framework treating sterile neutrinos as a dynamical cosmological fluid with a time-dependent equation of state.

• Microscopic Basis:

  • Starting from the Type I Seesaw Lagrangian, they define the active-sterile mixing via the vacuum mixing angle $\theta$.
  • They utilize the Boltzmann equation in an expanding FLRW background to model production.
  • Incomplete Thermalization: They argue that in the early Universe plasma, the effective mixing angle is suppressed by thermal potentials and damping rates. This leads to a production rate $\Gamma_s$ that is often smaller than the Hubble rate $H$, preventing full equilibrium.
  • Distribution Function: The resulting phase-space distribution is modeled as a suppressed Fermi-Dirac distribution:
    $$f_s(p) = \frac{\alpha}{e^{p/T_s} + 1}$$
    where $\alpha \leq 1$ is a thermalization parameter representing the fraction of the equilibrium population achieved.

• Macroscopic Fluid Dynamics:

  • The sterile neutrino energy density ($\rho_s$) and pressure ($P_s$) are computed by integrating the distribution function over phase space.
  • They derive the time-dependent equation of state parameter $w_s(a) = P_s(a)/\rho_s(a)$, which evolves from $1/3$ (relativistic) to $0$ (non-relativistic) as the scale factor $a$ increases.
  • These quantities are incorporated self-consistently into the Friedmann equations to determine the modification of the Hubble expansion rate ($H$).

3. Key Contributions

1. Analytic Framework for Partial Thermalization: The paper provides a closed-form analytic description of sterile neutrinos as a fluid with an evolving equation of state, bridging the gap between microphysical production parameters (mixing angle, mass) and macroscopic cosmological expansion.
2. Identification of Three Dynamical Regimes: The authors identify and analytically characterize three distinct phases in the evolution of the Hubble expansion rate correction ($\Delta H/H$):
   • Regime I (Relativistic Plateau): $T \gg m_s$. The fluid behaves as radiation; $\Delta H/H$ is constant ($\approx 0.009\alpha$).
   • Regime II (Transition): $T \sim m_s$. The equation of state drops rapidly from $1/3$ to $0$. This is where the deviation from standard $\Delta N_{\text{eff}}$ models is maximal.
   • Regime III (Matter-like Growth): $T \ll m_s$. The fluid behaves as pressureless matter ($\rho_s \propto a^{-3}$). Since radiation dilutes faster ($\rho_r \propto a^{-4}$), the sterile component's relative contribution to the total energy density grows during this phase.
3. Reinterpretation of Matter-Radiation Equality: The paper demonstrates that for GeV-scale sterile neutrinos, the population is effectively non-relativistic at the epoch of matter-radiation equality ($a_{\text{eq}}$). Therefore, they contribute as matter, not radiation, shifting the equality epoch in a manner distinct from standard $\Delta N_{\text{eff}}$ constraints.

4. Key Results

• Expansion Rate Correction: The fractional modification to the Hubble rate is derived as $\frac{\Delta H}{H} \approx \frac{1}{2} \frac{\rho_s}{\rho_r + \rho_m}$. For partial thermalization ($\alpha \sim 0.1$), this yields corrections at the permille to percent level, which are potentially observable.
• Shift in Equality Epoch:
  • The shift in the equality scale factor is given by $\frac{\Delta a_{\text{eq}}}{a_{\text{eq}}} \simeq f_s(2\eta - 1)$, where $f_s$ is the sterile energy fraction and $\eta$ is the relativistic fraction.
  • For GeV-scale sterile neutrinos, $\eta \ll 1$ at equality. Thus, the shift simplifies to $\frac{\Delta a_{\text{eq}}}{a_{\text{eq}}} \simeq -f_s$.
  • Unlike $\Delta N_{\text{eff}}$ (which delays equality), these sterile neutrinos advance the equality epoch because they act as extra matter.
• Observational Constraint: Using Planck constraints on the equality scale (constrained at the $\sim 1\%$ level), the authors derive a bound on the sterile energy fraction at equality:
  $$f_s(a_{\text{eq}}) \lesssim \mathcal{O}(10^{-2})$$
  This bound is qualitatively different from standard $\Delta N_{\text{eff}}$ bounds because it probes the matter-like contribution of the sterile sector.
• Free-Streaming: For GeV-scale masses, the free-streaming length at equality is negligible ($v_{\text{eq}} \sim 10^{-9}$), meaning they behave as Cold Dark Matter regarding structure formation, with their primary impact being on the background expansion history.

5. Significance

• Phenomenological Novelty: The work highlights a regime of sterile neutrino cosmology (partial thermalization with GeV masses) that is invisible to standard $\Delta N_{\text{eff}}$ analyses. It demonstrates that assuming a constant radiation-like equation of state can lead to incorrect interpretations of cosmological data.
• Analytic Transparency: By providing analytic estimates for the expansion correction and the equality shift, the paper offers a transparent link between microscopic mixing parameters (mass and mixing angle) and macroscopic observables, avoiding the need for full numerical integration of kinetic equations for background evolution.
• Complementary Constraints: The derived bound on $f_s$ provides a new, complementary probe for sterile neutrino models, specifically targeting the "matter-like" contribution of partially thermalized species, which is crucial for refining constraints on the parameter space of sterile neutrinos beyond current limits.

In conclusion, the paper establishes that sterile neutrinos with finite mass and incomplete thermalization must be treated as a dynamical fluid. This approach reveals that they can significantly alter the expansion history and the timing of matter-radiation equality, offering a new avenue for testing physics beyond the Standard Model using precision cosmological data.