Investigadores han utilizado simulaciones numéricas para modelar la evolución dinámica del Universo temprano dentro del marco del modelo inflacionario de Starobinsky. Este modelo, que describe una fase de expansión exponencial ultrarrápida del Universo primitivo, es uno de los más exitosos y sirve como referencia para otras teorías inflacionarias. El estudio se centró en la interacción entre el escalarón (el campo cuántico responsable de la inflación) y la producción de radiación a partir de su desintegración, asegurando la conservación exacta de la energía en un fondo de expansión.

Las simulaciones, que partieron de un Universo sin radiación y condiciones iniciales basadas en las aproximaciones estándar de *slow-roll*, mostraron que la dinámica inflacionaria está dominada por el escalarón. A pesar de la producción continua de radiación, esta se mantuvo subdominante durante la inflación, que se extendió por aproximadamente 60 *e*-folds. Tras la fase inflacionaria, el escalarón entró en una fase de oscilación amortiguada, acompañada de un rápido crecimiento de la radiación. Esta transición representa una "salida elegante" de la inflación y una entrada suave en la fase de recalentamiento.

Durante el recalentamiento, la densidad de energía de la radiación aumentó continuamente a medida que la energía del escalarón se agotaba. Finalmente, la densidad de energía de la radiación superó sustancialmente la energía residual del escalarón, indicando el inicio de una fase caliente dominada por la radiación. Este resultado satisface el criterio de Kofman-Yi para un recalentamiento exitoso, demostrando que la dinámica acoplada escalarón-radiación proporciona una realización autoconsistente de la salida de la inflación y la transición al dominio de la radiación en el modelo de Starobinsky.