Maximal parameter space of sterile neutrino dark matter with lepton asymmetries

This paper delineates the maximal parameter space for sterile neutrino dark matter by demonstrating that large, flavor-specific lepton asymmetries (with vanishing total asymmetry) can significantly enhance production rates and extend viable mixing angles by up to two orders of magnitude for masses up to 60 keV, while remaining consistent with cosmological constraints.

The Gist

The Big Mystery: What is Dark Matter?

Imagine the universe is a giant, invisible ocean. We can see the islands (stars and galaxies), but we know there is a massive amount of water (dark matter) holding them up, even though we can't see it. Scientists have been trying to figure out what this "water" is made of.

One popular suspect is a particle called the Sterile Neutrino. Think of a regular neutrino as a "ghost" that can pass through walls. A sterile neutrino is an even more ghostly version of that ghost. It doesn't talk to normal matter at all, except through a very faint, shy connection called "mixing."

The Problem: The Ghost is Too Heavy (or Too Light)

For a long time, scientists tried to figure out how these sterile neutrinos were made in the early universe.

• The Old Theory (The "Dodelson-Widrow" Mechanism): This was like saying, "They just randomly popped into existence." But when scientists looked at the sky with X-ray telescopes, they didn't see the ghosts. The theory was ruled out because it predicted too many ghosts, or ghosts that were too "warm" (moving too fast), which would have stopped galaxies from forming the way we see them today.

The New Idea: The "Lepton Asymmetry" Party

This paper introduces a new, exciting way these ghosts could have been made. It relies on a concept called Lepton Flavor Asymmetry.

The Analogy: The Seesaw Party
Imagine a massive party in the early universe with three types of guests: Electron-guests, Muon-guests, and Tau-guests.

• The Rule: The universe has a strict rule: the total number of guests must stay balanced (zero net charge).
• The Twist: In this new scenario, the Electron-guests are having a huge party, but the Muon-guests are having an equally huge anti-party (they are leaving in droves).
• The Result: If you count everyone, the total is zero. But locally, there is a massive imbalance. This is a "Lepton Flavor Asymmetry."

Usually, scientists thought these imbalances would get washed out (smoothed over) by neutrino oscillations (guests switching costumes) before the sterile neutrinos could be made. But this paper says: "Wait! What if the party happens before the costume switch?"

The Solution: A Resonant Dance Floor

The authors propose that if you have these huge, opposing parties (large asymmetries), it creates a special "resonance" on the dance floor.

1. The Resonance: Think of the dance floor as a giant swing. If you push a swing at just the right rhythm, it goes super high. The "asymmetry" acts like a perfect push. It makes the sterile neutrinos swing into existence much more efficiently than before.
2. The New Math: The authors wrote a new set of rules (equations) to calculate this. They realized that when the "push" is huge, the old math (which assumed the dance was slow and steady) breaks down. They had to invent a new way to count the dancers that accounts for the chaotic, fast-paced nature of this specific dance.

The Big Discovery: A New Safe Zone

By using this new math and the "Seesaw Party" scenario, the authors found a massive new safe zone for sterile neutrinos.

• Before: Sterile neutrinos had to be in a very narrow, tiny range of weights (masses) and shyness levels (mixing angles) to be dark matter. If they were outside this tiny range, they were either too heavy (ruled out by X-rays) or too light (ruled out by galaxy shapes).
• Now: With the "Lepton Asymmetry" boost, the safe zone expands by up to 100 times (two orders of magnitude).
  • It's like finding a hidden door in a locked room. Suddenly, the sterile neutrino can be heavier or lighter, and still be the perfect dark matter candidate.
  • They found this works for neutrinos up to 60 keV (a specific weight unit).

Why This Matters

1. It Solves the "Too Hot" Problem: In the old theory, the neutrinos were moving too fast (too "warm"), which would have smoothed out the tiny clumps of matter needed to make galaxies. The new method produces neutrinos that are "colder" (slower), which fits perfectly with how galaxies actually look today.
2. It's Testable: The paper gives scientists a roadmap.
   • X-ray Telescopes: We can look for specific signals from these neutrinos.
   • CMB (Cosmic Microwave Background): Future telescopes (like the Simons Observatory) might detect the "echo" of these massive parties in the early universe.
   • Structure Formation: We can look at how galaxies cluster to see if the "ghosts" are moving at the right speed.

The "Secret Sauce" (The Code)

The authors didn't just do the math on paper; they built a public software toolkit (called sterile-dm-lfa).

• Analogy: Imagine they didn't just tell you how to bake a cake; they gave you the recipe, the oven, and the measuring cups to anyone who wants to try.
• This allows other scientists to plug in different numbers and see exactly how these "ghosts" would behave, helping to test the theory against real-world data.

Summary in One Sentence

By realizing that the early universe might have had massive, opposing "parties" of particles that balanced out to zero, scientists have discovered a new, much larger playground where sterile neutrinos can exist as dark matter without breaking the laws of physics, and they've given everyone the tools to find them.

Technical Summary

1. Problem Statement

Sterile neutrinos ($\nu_s$) in the keV mass range are a compelling dark matter (DM) candidate. However, the minimal production mechanism, the Dodelson-Widrow (DW) mechanism (non-resonant oscillations), is ruled out by X-ray observations (due to decay limits) and structure formation constraints (due to being "too warm").

The Shi-Fuller mechanism offers a solution by utilizing primordial lepton flavor asymmetries ($L_\alpha$) to induce resonant production, allowing for smaller mixing angles. However, standard constraints from Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB) typically limit the total lepton asymmetry ($|\sum L_\alpha|$) to $\lesssim 10^{-3}$. If large individual flavor asymmetries ($|L_\alpha| \sim 0.1$) existed, they would wash out via neutrino oscillations at $T \sim 15$ MeV, potentially violating BBN/CMB bounds unless the total asymmetry remains negligible.

The Core Challenge: Previous studies often assumed small asymmetries or relied on approximations (averaged oscillations) that break down when $|L_\alpha|$ is large ($\gtrsim 10^{-3}$). There was no precise calculation of the maximal parameter space for sterile neutrino DM generated by large, flavor-specific asymmetries (e.g., $L_e = -L_\mu$) that sum to zero, nor was there a rigorous treatment of the non-averaged oscillation regime and the associated thermodynamic back-reactions.

2. Methodology

The authors developed a comprehensive framework to calculate sterile neutrino production in the presence of large lepton flavor asymmetries.

A. Kinetic Equations (Non-Averaged Oscillations)

• Derivation: They derived a semi-classical Boltzmann equation applicable to arbitrary lepton asymmetries.
• The Issue: For large $L_\alpha$, the resonance timescale becomes shorter than the neutrino oscillation length. Standard literature uses "averaged" oscillation probabilities, which fail in this regime.
• The Solution: They generalized the Boltzmann equation to account for non-averaged oscillations. Crucially, they identified an enhancement factor: while frequent collisions (Quantum Zeno effect) reset the active state, neutrinos accumulate over the mean free path and traverse the resonance. This accumulation compensates for the suppression, leading to an effective oscillation probability:
  $$P_{\text{eff}}(\nu_\alpha \to \nu_s) = \frac{1}{2} \frac{\Delta(p)^2 \sin^2 2\theta}{[\Delta(p) \cos 2\theta - V_\alpha(p, \mu)]^2 + (\Gamma_\alpha/2)^2}$$
  Note: Unlike previous averaged formulas, the denominator lacks the $\Delta^2 \sin^2 2\theta$ term, reflecting the resonance width being set by the interaction rate $\Gamma_\alpha$ rather than the oscillation frequency.
• Validation: This semi-classical equation was cross-validated against full Quantum Kinetic Equations (QKEs) and found to be in excellent agreement across both averaged and non-averaged regimes.

B. Thermodynamics and Chemical Potentials

• Chemical Potentials: For the first time in this context, the authors included the effects of large chemical potentials ($\mu_\alpha$) induced by $L_\alpha \sim 0.1$ in:
  • The neutrino interaction rate ($\Gamma_\alpha$).
  • The matter potential ($V_\alpha$).
  • The Equation of State (EoS) of the Early Universe (energy density $\rho$, pressure $P$, entropy $s$).
• QCD Transition: They handled the QCD phase transition ($T \sim 150$ MeV) by interpolating between perturbative QCD (quark-gluon plasma), Lattice QCD results (Taylor expansion of pressure), and the Hadron Resonance Gas (HRG) model. They verified that for the specific flavor directions considered ($L_e = -L_\mu$), pion condensation (a potential instability at high asymmetry) does not occur, ensuring the validity of their thermodynamic treatment.

C. Numerical Framework

• They released a public Python/Mathematica framework (sterile-dm-lfa) that solves the unintegrated Boltzmann equations.
• They also developed a simplified "narrow width approximation" code for rapid parameter scanning, which neglects back-reaction (valid for $m_s \gtrsim 5$ keV and large $L_\alpha$).

3. Key Contributions

1. Maximal Parameter Space Mapping: The paper delineates the allowed region in the ($m_s$, $\sin^2 2\theta$) plane for sterile neutrino DM, specifically for scenarios with zero total lepton asymmetry but large flavor asymmetries (e.g., $L_e = -L_\mu \approx 0.1$).
2. Theoretical Generalization: They provided the first rigorous derivation of the semi-classical Boltzmann equation for the non-averaged oscillation regime, correcting previous approximations that underestimated production rates at high asymmetries.
3. Thermodynamic Consistency: They demonstrated that large flavor asymmetries significantly alter the expansion rate and interaction rates of the Early Universe, and these effects must be included for accurate abundance calculations.
4. Public Code: The release of sterile-dm-lfa allows the community to compute sterile neutrino momentum distributions for arbitrary lepton asymmetry patterns, enabling dedicated structure formation analyses.

4. Results

• Expanded Viable Region: The inclusion of large flavor asymmetries ($|L_\alpha| \sim 0.1$) with vanishing total asymmetry opens up a parameter space where the mixing angle $\sin^2 2\theta$ can be up to two orders of magnitude smaller than in the standard flavor-universal asymmetry case.
• Mass Range: The viable mass range extends up to $m_s \lesssim 60$ keV.
• Production Temperature: For the maximal allowed couplings (the "floor" of the parameter space), sterile neutrinos are produced at temperatures below the QCD transition ($T \lesssim 100$ MeV). This implies the results are robust against uncertainties in the QCD equation of state.
• Spectral Properties:
  • Contrary to some previous qualitative claims, increasing the lepton asymmetry $L_\alpha$ leads to a warmer sterile neutrino spectrum (higher average momentum).
  • Reasoning: Larger $L_\alpha$ shifts the resonance to lower temperatures. Production at lower $T$ means less entropy dilution from subsequent particle annihilations, resulting in a higher $p/T$ ratio today.
• Constraints:
  • X-rays: The gray region in their Fig. 1 is excluded by current X-ray limits.
  • Structure Formation: While a dedicated Lyman-$\alpha$ analysis is pending, they provide an illustrative guide using an equivalent warm dark matter mass ($m_{\text{WDM}}$). The "warmer" spectra from large asymmetries push the constraints to slightly higher masses, but the region remains testable.
  • BBN/CMB: The scenario is consistent with current bounds because the total asymmetry is zero, and the flavor asymmetries wash out before BBN ($T < 1$ MeV).

5. Significance

• Revival of Sterile Neutrino DM: This work significantly revitalizes the sterile neutrino DM hypothesis by showing that a vast, previously unexplored region of parameter space is viable under well-motivated cosmological conditions (Affleck-Dine leptoflavorgenesis).
• Testability: The newly opened parameter space is highly testable by upcoming experiments:
  • X-ray: eXTP and Athena can probe the smaller mixing angles.
  • CMB: The Simons Observatory can test the specific lepton flavor asymmetry patterns ($L_e = -L_\mu$) required for this scenario.
  • Structure Formation: Future Lyman-$\alpha$ and satellite galaxy surveys will be able to distinguish the specific momentum distributions predicted by this model.
• Theoretical Benchmark: The paper sets a new standard for calculating resonant production in the Early Universe, moving beyond averaged approximations and providing a toolset for handling large chemical potentials in cosmological simulations.

In summary, the authors demonstrate that sterile neutrinos with masses up to 60 keV can constitute all dark matter if the Early Universe possessed large, flavor-asymmetric lepton asymmetries that summed to zero. This scenario is consistent with all current observational bounds and offers a rich target for next-generation cosmological and astrophysical probes.