⚛️ phenomenology

Decay-Assisted Leptogenesis during Preheating

By employing quantum Boltzmann equations to model the interplay between fermionic saturation and decay during preheating, this study demonstrates that the simultaneous production and decay of heavy right-handed neutrinos can enhance the final baryon asymmetry by two orders of magnitude compared to scenarios where decay occurs only after production ceases.

Original authors: Wenbin Zhao

Published 2026-10-14

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Wenbin Zhao

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

The universe we see today is made of matter, yet the laws of physics suggest that the Big Bang should have created equal amounts of matter and its opposite, antimatter. If that balance had remained perfect, the two would have annihilated each other instantly, leaving behind only light. The fact that we exist means something tipped the scales, creating a tiny surplus of matter that survived to form stars, planets, and us. One of the leading ideas for how this happened involves heavy, invisible particles that existed only in the very first moments after the birth of the cosmos. These particles, known as right-handed neutrinos, are thought to have decayed in a way that favored matter over antimatter. This process, called leptogenesis, is a crucial step in explaining our existence, but scientists have long struggled to calculate exactly how much matter was produced during the chaotic, violent expansion that followed the Big Bang.

For decades, researchers have studied two main ways these heavy particles could have been created. In one scenario, the universe was hot enough that these particles formed naturally through collisions, like steam condensing from a hot bath. In another, they were created directly by the energy of the field that drove the rapid expansion of the early universe, a period known as inflation. A recent study by Wenbin Zhao revisits this second scenario, focusing on a specific phase called preheating, where the energy of the inflation field is transferred to other particles in a burst of activity. The researcher used advanced mathematical tools to track how these heavy neutrinos were born, lived, and died during this turbulent time, paying close attention to a subtle effect that had been overlooked in previous calculations: the fact that these particles can decay while they are still being created.

The core of the discovery lies in understanding how these heavy particles behave when they are crowded together. Because they are fermions, a type of particle that follows strict rules about how many can occupy the same space, their production can hit a natural limit. Imagine a crowded room where no new people can enter until someone leaves; this is the "Pauli blocking" effect that stops the creation of new particles once the available spots are filled. In earlier models, scientists often assumed that these heavy neutrinos would wait until they were all created before they began to decay. However, Zhao's work shows that if these particles decay while the production is still happening, they empty out the crowded spots, allowing more particles to be created. This continuous cycle of creation and decay prevents the process from stalling, leading to a much larger number of heavy particles than previously thought.

By running detailed simulations of this process, the study found that this mechanism can boost the final amount of matter in the universe by a factor of one hundred compared to models where decay happens only after production stops. The researchers tested this idea across different speeds of decay and different strengths of the forces involved. They found that whether the heavy particles decayed almost instantly, waited until the universe cooled slightly, or decayed much later, the inclusion of this decay-during-production effect consistently led to a significant increase in the matter surplus. In some specific cases, the enhancement was so strong that it could easily explain the observed amount of matter in the universe, even with parameters that would have failed in older models.

This work unifies several different ways scientists have tried to understand the origin of matter into a single, consistent picture. It demonstrates that the timing of a particle's death is just as important as its birth when calculating the history of the cosmos. The study does not claim to have solved the mystery of the universe's matter entirely, but it provides a much clearer and more powerful tool for understanding how the heavy particles of the early universe could have tipped the balance in our favor. By showing that the decay of these particles can actively help create more of them, the research opens up new possibilities for how the universe became the matter-filled place we inhabit today.

Technical Summary: Decay-Assisted Leptogenesis during Preheating

Problem Statement
The paper addresses the mechanism of non-thermal leptogenesis occurring during the preheating phase of the early universe. While standard thermal leptogenesis relies on the production of heavy right-handed neutrinos (RHNs) from a thermal bath, non-thermal scenarios produce these particles directly from the oscillating inflaton field. A critical limitation in previous treatments of this scenario is the treatment of RHN decays. Often, models assume that RHNs are produced first and decay only after production ceases, or they utilize collisionless Bogoliubov methods that neglect the interplay between production and decay. This neglects the potential for RHN decays to deplete occupation numbers, thereby mitigating Pauli blocking and allowing for further particle production. The paper seeks to unify the description of RHN production and decay to determine how the RHN lifetime influences the final baryon asymmetry of the universe.

Methodology
The author employs quantum Boltzmann equations to model the coupled evolution of RHN production and decay within a homogeneous, isotropic cosmological background.

Background Dynamics: The inflaton field ( $\phi$ ) is modeled with a quadratic potential ( $V(\phi) = \frac{1}{2}m_\phi^2\phi^2$ ), undergoing harmonic oscillations after inflation. The inflaton couples to Standard Model fermions and heavy RHNs ( $N_i$ ).
Quantum Kinetic Framework: The system is described by equations of motion for the occupation number ( $f_k$ ) and the coherence function ( $C_k$ ) between particle and antiparticle states. The equations incorporate a time-dependent mass term ( $M_t = M + g\phi$ ) which sources production, and a momentum-dependent decay width ( $\Gamma_k$ ).
Relaxation-Time Approximation: To handle the unstable nature of the fermions, the author adopts the relaxation-time approximation, allowing for a consistent treatment of both the non-adiabatic production (driven by the time-varying mass) and the subsequent decay.
Lepton Asymmetry Calculation: The lepton asymmetry yield ( $Y_L$ ) is calculated by integrating the CP-violating decay rates of the RHNs, governed by the Davidson-Ibarra (DI) bound, while accounting for the entropy production from both inflaton and RHN decays.

Key Contributions

Unified Framework: The paper provides the first unified calculation using quantum Boltzmann equations that covers three distinct decay regimes:
- Decay occurring during RHN production.
- Decay occurring after production but before reheating.
- Decay occurring after reheating.
Resolution of Pauli Blocking: The study explicitly demonstrates how RHN decays can mitigate the Pauli blocking effect. In standard non-thermal scenarios, RHN occupation numbers saturate at unity (the Pauli bound) after a few inflaton oscillations, halting further production. The author shows that if the decay width is sufficiently large, decays deplete the occupation numbers, allowing the inflaton to produce more RHNs.
Quantitative Enhancement: The framework quantifies the enhancement of the final baryon asymmetry resulting from this "decay-assisted" production mechanism.

Results
The author presents numerical simulations across different hierarchies of the inflaton decay width ( $\Gamma_\phi$ ) and the RHN decay width ( $\Gamma_N$ ):

Case A (Decay during production): When $\Gamma_N \gg m_\phi$ (the inflaton mass), RHNs decay rapidly between production events. This effectively resets the occupation numbers, allowing the total accumulated number of RHNs to be the sum of individual production events rather than being capped by the Pauli bound. This leads to an enhancement of the final baryon asymmetry by approximately two orders of magnitude (a factor of $\sim 100$ ) compared to scenarios where decays are negligible during production.
Case B (Decay before reheating): When $m_\phi/n_{max} < \Gamma_N < m_\phi$ , decays occur after production slows but before the inflaton decays. The asymmetry depends on the specific amount of occupation depleted.
Case C (Decay after reheating): When $\Gamma_N \ll \Gamma_\phi$ , the RHNs dominate the energy density temporarily before decaying. The final asymmetry is determined by the RHN decay width and is insensitive to the inflaton decay width.
Parameter Sensitivity: The study finds that while reducing the RHN mass ( $M$ ) increases production efficiency, it simultaneously reduces the CP asymmetry per decay. These factors largely compensate, resulting in a relatively stable baryon yield. However, increasing $M$ leads to a suppression of production that dominates, reducing the final asymmetry.

Significance and Claims
The paper claims that incorporating the decay of right-handed neutrinos into the preheating dynamics is essential for an accurate calculation of non-thermal leptogenesis. The primary significance lies in demonstrating that decay-assisted production can enhance the final baryon asymmetry by two orders of magnitude compared to models where production and decay are treated as sequential, non-overlapping processes.

The author emphasizes that their approach recovers previous results found in the literature (such as the Bogoliubov limit where decays are neglected) as specific limits of their more general equations. They assert that this unified framework allows for a consistent study of the transition between different decay regimes. The paper concludes by noting that the methodology is generalizable to other scenarios involving parametric production and decay of bosons (e.g., axions or supermassive gauge bosons), though the current work focuses strictly on the fermionic RHN system. The author remains modest regarding potential extensions, suggesting that different inflaton potentials or more accurate estimates of the CP asymmetry (beyond the static DI bound) could further refine the parameter space for successful leptogenesis.

More like this