⚛️ general relativity

Graceful Exit to Radiation Domination in the Starobinsky Model: A Dynamical Realization

This paper demonstrates that the Starobinsky model naturally achieves a graceful exit from inflation and a successful transition to a radiation-dominated era through the self-consistent, energy-conserving dynamical evolution of scalaron decay and radiation production, which maintains slow-roll conditions during inflation and satisfies the Kofman-Yi criterion for reheating.

Original authors: Bhargabi Saha, Malay K. Nandy

Published 2026-10-02

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Bhargabi Saha, Malay K. Nandy

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

In the earliest moments of our universe, a fraction of a second after the beginning, space itself underwent a period of unimaginable expansion known as inflation. This rapid stretching smoothed out the cosmos and set the stage for everything that followed, but it left a critical question unanswered: how did this frozen, empty state of rapid expansion suddenly transform into the hot, bustling soup of particles that makes up the world we see today? For decades, physicists have relied on a leading theory called the Starobinsky model to describe this initial burst of growth. In this framework, the force driving inflation is not a mysterious new particle, but a modification to the very laws of gravity, which can be mathematically described as a heavy, invisible field called the scalaron. While this model successfully explains how the universe grew, the transition from that cold, expanding phase to the hot, radiation-filled era required for the Big Bang to continue has remained a complex puzzle. Scientists needed to understand exactly how the energy stored in this invisible field could be drained away to create the light and heat of the early universe without breaking the delicate balance of cosmic expansion.

A team of researchers at the Indian Institute of Technology Guwahati has now mapped out this transition in detail, showing how the universe gracefully exits the inflationary phase and enters a new era. They treated the early universe as a dynamic system where the invisible scalaron field and the newly created radiation are locked in a constant exchange of energy. Rather than assuming the radiation simply appeared or that the scalaron vanished instantly, they built a computer simulation that followed the energy transfer moment by moment. They started their simulation with a universe that contained only the scalaron field, completely devoid of light or heat, and let the laws of physics run their course. As the simulation progressed, the scalaron field began to decay, breaking down into particles of radiation, much like a heavy object slowly crumbling into dust. The researchers carefully tracked how this decay affected the expansion of the universe, ensuring that the energy lost by the scalaron was exactly equal to the energy gained by the radiation, a strict rule of conservation that governs all physical processes.

The results of this simulation revealed a clear and orderly story. During the inflationary period, the universe was dominated entirely by the scalaron field. Even as radiation began to form from the decay of the field, it remained a tiny, insignificant fraction of the total energy, barely affecting the rapid expansion. The universe continued to stretch smoothly for approximately 60 cycles of doubling in size, a duration that matches the requirements for the large-scale structure we observe today. Once this period of inflation ended, the behavior of the universe changed dramatically. The scalaron field, which had been holding the expansion steady, began to oscillate, vibrating back and forth as it lost energy. This vibration triggered a rapid and massive production of radiation. The energy that had been stored in the invisible field was transferred efficiently into a sea of particles, causing the radiation density to surge while the scalaron energy dwindled.

The simulation showed that this process naturally led to a state where the radiation completely overwhelmed the remaining scalaron field. By the time the transfer was complete, the energy density of the radiation was vastly larger than the leftover energy of the scalaron, satisfying a specific condition known as the Kofman-Yi criterion, which physicists use to define a successful transition to a hot universe. The researchers found that this outcome held true regardless of how fast or slow the scalaron decayed; in every scenario they tested, the universe ended up filled with radiation and ready for the next stage of cosmic history. The study confirms that the Starobinsky model provides a self-consistent and smooth path from the cold, empty expansion of inflation to the hot, dense beginning of the Big Bang, resolving a long-standing question about how the early universe reheated itself without needing any external intervention or fine-tuned adjustments.

Technical Summary: Graceful Exit to Radiation Domination in the Starobinsky Model

Problem Statement
The Starobinsky model, driven by a quadratic Ricci scalar correction ( $R^2$ ) in the gravitational action, is a benchmark for successful cosmic inflation consistent with current observational data. However, a critical aspect of early universe cosmology is the transition from the inflationary epoch to the subsequent radiation-dominated era. This transition requires a "graceful exit" where the energy stored in the inflaton field (the scalaron) is efficiently transferred to radiation, satisfying the Kofman-Yi criterion ( $\rho_{\text{radiation}} \gg \rho_{\text{scalaron}}$ ). While the inflationary dynamics of the Starobinsky model are well-established, the dynamical realization of this energy transfer—specifically how the universe evolves from a scalaron-only state to a radiation-dominated state while maintaining exact energy conservation and accounting for the backreaction of created radiation on the expansion—remains a subject of investigation.

Methodology
The authors formulate a self-consistent dynamical system within the Starobinsky framework to model the coupled evolution of the scalaron field and radiation. The approach involves:

Theoretical Framework: The model is treated in the Einstein frame, where the $R^2$ theory is equivalent to Einstein gravity minimally coupled to a scalar field (the scalaron, $\phi$ ) with a specific potential $V(\phi)$ . The scalaron is coupled to a bosonic matter field $\chi$ , with the dominant decay channel identified as $\phi \to \chi\chi$ .
Constrained Radiation Production: The production rate of radiation is constrained by the Heisenberg uncertainty principle and kinematic conservation laws. The authors derive a maximum possible radiation energy density ( $w_{\text{max}}$ ) based on the minimum position uncertainty of the decay products. Consequently, the energy transfer rate $q(t)$ is bounded; it is constant when the scalaron energy density $u$ is high, but scales linearly with $u$ once $u$ drops below $w_{\text{max}}$ .
Coupled Dynamical Equations: A system of three coupled nonlinear differential equations is derived:
- The rate of change of scalaron energy density ( $\dot{u}$ ), including a sink term $-q(t)$ and Hubble expansion depletion.
- The rate of change of radiation energy density ( $\dot{\rho}$ ), including a source term $+q(t)$ and Hubble expansion dilution.
- The Friedmann equation, which sources the Hubble expansion rate $H$ using the total energy density ( $u + \rho$ ), thereby incorporating the backreaction of the created radiation on the background geometry.
Numerical Integration: Starting from a radiation-free universe with initial conditions derived from standard slow-roll approximations (targeting $\sim 60$ e-folds), the authors perform exact numerical integration of the coupled system. They analyze the evolution across six different scalaron decay rates ( $\Gamma/\mu$ ranging from $1$ to $10^{-5}$ ).

Key Contributions and Results
The numerical integration yields the following results regarding the dynamical evolution of the early universe:

Inflationary Phase: During the slow-roll phase, the scalaron energy density ( $u$ ) overwhelmingly dominates the radiation energy density ( $\rho$ ). Despite the continuous production of radiation via scalaron decay, $\rho$ remains subdominant (e.g., $\tilde{\rho} \sim 10^{-11}$ to $10^{-16}$ relative to the Planck scale for different $\Gamma$ ). The standard slow-roll conditions are preserved, and the universe undergoes approximately 61.5 e-folds of inflation. The backreaction of the radiation on the expansion rate is negligible during this phase.
Graceful Exit: Upon the completion of inflation (identified by the slow-roll parameter $\epsilon \approx 1$ ), the scalaron field transitions into a damped oscillatory phase. This marks the onset of the reheating regime.
Reheating Dynamics: In the post-inflationary phase, the scalaron energy is rapidly depleted as it transfers to radiation. The radiation energy density grows significantly, eventually reaching a peak ( $\rho_r$ ).
Kofman-Yi Criterion: At the peak of radiation density, the residual scalaron energy density ( $u_r$ ) becomes negligibly small compared to the radiation density ( $\rho_r \gg u_r$ ). For the case $\Gamma = \mu$ , the peak radiation density corresponds to a reheating temperature $T_r \approx 4.5 \times 10^{13}$ GeV; for $\Gamma = 10^{-5}\mu$ , $T_r \approx 6.7 \times 10^{12}$ GeV.
Robustness: These results hold across the range of decay rates tested. The transition is smooth, and the final state is a hot, radiation-dominated universe with negligible residual scalaron energy.

Significance
The paper claims to demonstrate that the Starobinsky model provides a coherent, self-consistent dynamical realization of a graceful exit from inflation. By treating the scalaron-radiation-gravity system as a coupled dynamical system with exact energy conservation and including the backreaction of radiation, the authors show that the universe naturally evolves from a scalaron-dominated inflationary phase to a radiation-dominated phase. This transition satisfies the Kofman-Yi criterion, ensuring that the energy budget is effectively transferred to radiation, thereby establishing the necessary conditions for the subsequent hot Big Bang evolution (including nucleosynthesis and baryogenesis) without requiring ad-hoc assumptions about the reheating mechanism. The study reinforces the viability of the Starobinsky model as a complete description of the early universe, from inflation through to the radiation era.

Technical Summary: Graceful Exit to Radiation Domination in the Starobinsky Model

More like this