⚛️ phenomenology

Heavy neutral leptons beyond the BBN bound: probing the lepton asymmetry of the Universe

This paper proposes a scenario where heavy neutral leptons with lifetimes typically excluded by Big Bang Nucleosynthesis constraints can evade these limits through large lepton flavor asymmetries that balance neutron-proton conversion, thereby opening a significant parameter space for upcoming accelerator searches like SHiP and linking potential discoveries to the primordial lepton asymmetry of the Universe.

Original authors: Kensuke Akita, Maksym Ovchynnikov, Vsevolod Syvolap

Published 2026-08-12

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Kensuke Akita, Maksym Ovchynnikov, Vsevolod Syvolap

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 that started with a tiny, incredibly hot bang. In the very first few minutes of this expansion, a cosmic kitchen was cooking up the first atomic nuclei, a process scientists call Big Bang Nucleosynthesis (BBN). It's like a high-stakes recipe where the temperature and timing have to be perfect to get the right amount of helium and hydrogen. If the recipe goes wrong, the Universe would look completely different today—maybe no stars, no planets, and certainly no us.

To understand this recipe, we need to know about "neutrinos." These are ghostly, tiny particles that barely interact with anything, zipping through the Universe like invisible ghosts. They are crucial because they help regulate the temperature of the cosmic kitchen. Recently, scientists have been hunting for a new, heavier cousin of the neutrino called a "Heavy Neutral Lepton" (HNL). Think of an HNL as a heavy, shy guest at the party who might show up, hang out for a while, and then disappear. If these guests exist, they could explain some mysteries about why the Universe has more matter than antimatter. However, there's a catch: if these heavy guests show up and decay (disappear) too late in the cosmic timeline, they might ruin the helium recipe, creating too much of it. This has led many scientists to believe that these heavy guests simply cannot exist in the timeframes that our most sensitive particle detectors are currently looking for.

But what if the recipe isn't ruined because the guests are there, but because the guests are bringing a specific kind of imbalance to the party? That is the question a team of physicists asked in a new study. They wondered: what if the heavy guests aren't just showing up randomly, but if there's a massive, hidden imbalance between the "guests" and their "anti-guests" (particles with opposite charges)?

The paper proposes a clever scenario where this imbalance acts like a cosmic counter-balance. Usually, when these heavy particles decay, they shoot out charged pions (tiny particles made of quarks). If the heavy particles and anti-particles are present in equal numbers, they shoot out equal amounts of positive and negative pions. These pions then crash into protons and neutrons, scrambling the recipe and turning too many protons into neutrons, which leads to an overabundance of helium. This is the "BBN bound" that usually rules out these particles.

However, the authors suggest a twist: what if the Universe started with a huge, pre-existing imbalance, favoring one type of heavy particle over the other? In this scenario, the heavy particles would decay into positive pions, while the fewer anti-particles would decay into negative pions. Because there are more heavy particles than anti-particles, there would be an excess of positive pions. These extra positive pions act like a cosmic reset button, converting the extra neutrons back into protons, effectively undoing the damage and saving the helium recipe. It's like having a team of chefs who accidentally add too much salt, but a second team arrives with just the right amount of sugar to cancel it out perfectly, leaving the dish tasting exactly as it should.

The paper uses complex computer simulations to test this idea. They found that if this imbalance exists, it opens up a whole new "safe zone" for these heavy particles. Previously, scientists thought these particles couldn't exist with lifetimes longer than about 0.02 seconds because they would ruin the helium. The new study suggests that with the right imbalance, these particles could live much longer—up to 1 second—and still fit perfectly within the observed laws of the Universe.

This discovery is exciting because it means the upcoming experiments, like the SHiP experiment at CERN, might actually find these particles in a range of masses and lifetimes that were previously thought to be impossible. If they do find them, it wouldn't just be a discovery of a new particle; it would be a direct measurement of the Universe's hidden "flavor asymmetry"—a clue to why the Universe is made of matter at all. The paper also suggests that if this scenario is true, it could leave other fingerprints, such as a slight change in the number of neutrino types we detect from the early Universe or even ripples in spacetime (gravitational waves) from a phase transition in the early Universe.

In short, the paper doesn't prove these particles exist, but it suggests that the door isn't closed. It argues that the "no-entry" sign posted by the helium recipe might be wrong if we account for a hidden imbalance in the Universe's history. It turns a dead end into a promising path, inviting scientists to look for these heavy neutrinos in a new, wider range of possibilities. If they are found there, it would confirm that the early Universe was a much more lopsided and dynamic place than we ever imagined.

Technical Summary: Heavy neutral leptons beyond the BBN bound: probing the lepton asymmetry of the Universe

Problem Statement
Standard cosmology imposes a stringent constraint on Heavy Neutral Leptons (HNLs) with hadronic decay channels and lifetimes $\tau_N \gtrsim 0.02$ s. In the early Universe, HNLs with masses $m_N > m_l + m_\pi$ decay into charged pions ( $N \to \pi^+ l^-$ and $\bar{N} \to \pi^- l^+$ ). In a symmetric population ( $n_N = n_{\bar{N}}$ ), the injected mesons drive strong interactions ( $\pi^- + p \to n + \pi^0/\gamma$ ) that convert protons to neutrons far faster than weak rates. This drives the neutron-to-proton ratio ( $n/p$ ) toward an isospin-symmetric value of $\sim 1$ , leading to a catastrophic overproduction of helium ( $Y_p \approx 0.6$ ) and deuterium, violating Big Bang Nucleosynthesis (BBN) observations. Consequently, this "BBN bound" excludes a significant portion of the parameter space targeted by upcoming accelerator experiments like SHiP, DUNE, and Hyper-K, particularly for $m_N \lesssim 1$ GeV and $\tau_N \lesssim 1$ s. Previous attempts to evade this bound relied on suppressing the HNL abundance entirely via low-scale reheating, which leaves no correlated cosmological signatures.

Methodology
The authors propose a mechanism where a large primordial lepton flavor asymmetry ( $L_\alpha$ ) allows HNLs to evade the BBN bound without suppressing their population. The study employs a coupled, three-stage numerical pipeline to track the evolution of the Universe from HNL production down to BBN:

Production Stage ( $T \sim \text{GeV} \to 20$ MeV): A momentum-resolved Boltzmann evolution calculates the production of Dirac HNLs ( $N$ ) and their antiparticles ( $\bar{N}$ ) via mixing with active neutrinos. Crucially, the calculation accounts for the redistribution of conserved charges (baryon, electric, and lepton flavor) among all plasma species. The presence of large flavor asymmetries creates a chemical potential difference, resulting in an unequal population ratio $r \equiv n_N / n_{\bar{N}} \neq 1$ .
Oscillation and Flavor Dynamics Stage ($20$ MeV $\to 5$ MeV): Using a quantum-kinetic framework (COFLASY), the authors evolve the neutrino flavor system in the presence of the decaying HNLs. This stage tracks how the asymmetric decays inject lepton flavor asymmetries into the active sector and how neutrino oscillations redistribute these asymmetries.
Transport and BBN Stage ( $T < 5$ MeV): A hybrid approach is used. A momentum-resolved Boltzmann transport code ( $\nu$ DSMC) follows the spectral evolution of neutrinos, antineutrinos, and decay products (mesons, muons) to accurately compute the weak $p \leftrightarrow n$ rates and the effective number of neutrino species ( $N_{\text{eff}}$ ). This is coupled with a nuclear network code to compute primordial abundances ( $Y_p$ , D/H). The code explicitly tracks the injection of non-thermal neutrinos and the meson-driven conversion rates.

Key Contributions and Mechanism
The central contribution is the demonstration that a charge asymmetry in the HNL population can cancel the destructive effects of meson injection on BBN:

Meson Charge Imbalance: In a Dirac scenario with $r \neq 1$ , the decays inject unequal numbers of $\pi^+$ and $\pi^-$ . Since $\pi^-$ drives the conversion $p \to n$ (overproducing neutrons) and $\pi^+$ drives $n \to p$ , an excess of $\pi^+$ can reverse the neutron enrichment.
Restoring Standard Abundances: By tuning the asymmetry ratio $r$ to a specific value ( $r_{\text{req}} \approx 4\text{--}10$ ), the net effect of meson injection can be tuned to leave the $n/p$ ratio at the standard BBN value when mesons disappear.
$N_{\text{eff}}$ Cancellation: The same lepton flavor asymmetries that create the HNL imbalance also increase the neutrino energy density (raising $N_{\text{eff}}$ ). Conversely, the decay of HNLs into non-thermal neutrinos and the subsequent heating of the electromagnetic plasma tend to lower $N_{\text{eff}}$ . The authors show that these two effects can cancel, keeping the shift $\Delta N_{\text{eff}}$ within observational bounds.
Kinematic Window: The mechanism is effective for $m_l + m_\pi \lesssim m_N \lesssim 1.5$ GeV. Above $\sim 1.5$ GeV, multi-hadron final states and neutral kaons ( $K_L$ ) wash out the charge correlation, preventing the necessary imbalance.

Results
The study identifies four distinct regions in the HNL parameter space ( $m_N, \tau_N$ ) that are consistent with BBN and Cosmic Microwave Background (CMB) constraints:

Region 1 (Invisible): Parameters where $Y_p$ , D/H, and $N_{\text{eff}}$ remain indistinguishable from standard cosmology.
Region 2 ( $N_{\text{eff}}$ Deficit): A detectable negative shift in $N_{\text{eff}}$ (e.g., $\Delta N_{\text{eff}} \approx -0.08$ ), potentially explaining current observational hints of $N_{\text{eff}} < 3.044$ .
Region 3 (Modified Relic Neutrino Background): A scenario where the relic neutrino sea is enhanced by $>30\%$ in number density due to the surviving asymmetry and non-thermal tails, even if $N_{\text{eff}}$ appears standard. This offers a target for experiments like PTOLEMY.
Region 4 (Gravitational Waves): A subset of the parameter space where the large lepton flavor asymmetries drive the QCD phase transition to be first-order, potentially sourcing a stochastic gravitational wave background.

The paper provides benchmark points (Table I) demonstrating that these scenarios are compatible with current limits from SHiP, DUNE, and Hyper-K, effectively reopening the previously excluded domain for HNL searches.

Significance and Claims
The paper claims to identify the first mechanism by which a mixing-produced population of hadronically decaying particles can evade BBN constraints without suppressing the particle abundance.

Discovery Implications: An HNL discovery in the previously excluded domain would not merely be a new particle detection; it would serve as a probe of the primordial lepton flavor asymmetries of the Universe.
Correlated Signatures: The detection of an HNL would constrain the primordial asymmetries to a one-parameter family. The specific sharing of asymmetry between electron and tau flavors would determine whether the QCD transition was first-order (GW signal) and the magnitude of the relic neutrino background enhancement (PTOLEMY signal).
Theoretical Consistency: The scenario requires the HNL to be Dirac (to maintain distinct $N$ and $\bar{N}$ populations) and lepton number to be conserved to a high degree ( $\delta m \lesssim 10^{-15}$ eV). This links accelerator searches to fundamental questions about the nature of neutrino mass and the thermal history of the Universe.

The authors emphasize that their calculation is the first to simultaneously track large lepton flavor asymmetries through flavor dynamics, momentum-resolved neutrino evolution, and meson-driven nucleosynthesis, establishing the existence of these solutions and their correlated cosmological imprints.

More like this