Simulating first-order phase transition during inflation

1. Problem Statement

The paper addresses the long-standing "graceful exit problem" associated with ending cosmic inflation via a first-order phase transition (FoPT).

• Historical Context: In the original "Old Inflation" model, inflation ended via bubble nucleation. However, the expansion of the universe was too rapid for bubbles to percolate and collide, preventing the universe from exiting the inflationary phase.
• Modern Challenge: While coupling the inflaton to a spectator field undergoing an FoPT offers a potential "graceful exit," constructing concrete models is difficult. The primary hurdle is ensuring that bubble nucleation is slow enough to avoid disrupting early inflation but fast and dense enough near the end of inflation to successfully terminate the expansion and reheat the universe.
• Specific Gap: Previous analytical estimates suggested that the standard percolation condition (false vacuum fraction dropping below $1/e$) might be misleading in a vacuum-dominated era. Furthermore, there was a lack of rigorous 3D lattice numerical simulations to verify if such a mechanism could dynamically succeed and what specific gravitational wave (GW) signatures it would produce.

2. Methodology

The authors propose a specific model within the Starobinsky inflation framework and validate it using high-resolution 3D lattice simulations.

A. Theoretical Model Construction

• Framework: The model operates at the Grand Unified Theory (GUT) scale within Starobinsky inflation ($R^2$ gravity).
• Two-Field System:
  1. Inflaton ($\chi$): Follows the Starobinsky plateau potential $W(\chi)$.
  2. Phase-Transition Field ($\phi$): Governed by a potential $V(\phi)$ with a barrier.
• Key Mechanism (Dynamic Barrier Control): The height of the potential barrier for $\phi$ is dynamically controlled by the rolling inflaton $\chi$ via an exponential relationship:
  $$a = e^{\beta \chi/M_{Pl}} - 1$$
  • Early Inflation: As $\chi$ is large, the barrier is high, suppressing bubble nucleation almost entirely.
  • Late Inflation: As $\chi$ rolls down, the barrier drops exponentially, triggering a massive, rapid burst of bubble nucleation just before inflation ends.
• UV Completion: The authors construct a Supergravity (SUGRA)-inspired Effective Field Theory (EFT) to provide a consistent ultraviolet completion for this two-field potential, utilizing a no-scale SUGRA structure with chiral superfields.

B. Numerical Simulation Setup

• Code: Simulations were performed using pystella, a lattice solver for partial differential equations on GPUs.
• Grid: A 3D uniform Cartesian grid with resolution $N^3 = 512^3$ and periodic boundary conditions.
• Initial Conditions:
  • The $\chi$ field starts as a Gaussian random field with mean $\chi \approx 5.42 M_{Pl}$ (corresponding to $N=60$ e-folds).
  • The $\phi$ field starts at zero.
• Trigger Point: The simulation tracks the evolution until $\chi$ reaches a critical value ($\chi_* \approx 1.69 M_{Pl}$), where the nucleation rate becomes sufficient to satisfy the percolation condition (requiring $\sim 4$ bubbles per comoving Hubble volume to collide within one e-folding).
• Physics Calculated: The code solves the coupled Klein-Gordon equations for the fields and the Friedmann equations for the background. It calculates the energy-momentum tensor to derive the gravitational wave (GW) power spectrum ($\Omega_{GW}$).

3. Key Contributions

1. Novel Model Proposal: Introduced a GUT-scale FoPT in Starobinsky inflation where the barrier height is exponentially suppressed by the inflaton, solving the timing issue of bubble nucleation (suppressing early nucleation while ensuring late, rapid percolation).
2. First 3D Lattice Verification: Provided the first full 3D numerical simulation of an inflation-triggered FoPT, moving beyond analytical approximations (like the "instantaneous source" approximation) to include non-linear effects and full field dynamics.
3. Validation of Analytical Predictions: Confirmed that the "instantaneous source" approximation holds for the GW spectrum in this specific model, validating previous theoretical work.
4. Discovery of Distinctive GW Signature: Identified a unique high-frequency oscillatory feature in the GW energy density spectrum, caused by the phase accumulation of modes exiting and re-entering the Hubble horizon during the rapid phase transition.

4. Results

• Dynamics of Phase Transition:
  • The simulation successfully demonstrated that the universe remains in a slow-roll inflationary phase until the inflaton reaches the critical value.
  • Once the barrier drops, bubbles nucleate rapidly, expand, and collide, filling the simulation volume within a fraction of an e-folding.
  • The equation of state ($w$) transitions from vacuum-dominated ($w \approx -1$) to a value around $-0.6$ with pronounced oscillations due to the coupling between fields and inflaton oscillations.
• Gravitational Wave Spectrum:
  • The resulting GW energy density spectrum ($\Omega_{GW}$) matches previous analytical estimations.
  • Distinctive Feature: In the high-frequency regime ($\log_{10}[kR_*] > 1$), the spectrum exhibits a distinct oscillatory pattern. This is attributed to the coherent superposition of GW modes that exited the horizon at different times and accumulated distinct phases before re-entering.
  • The peak frequency of the signal scales as $f_p \propto e^{-N_*}$, where $N_*$ is the number of e-folds remaining when the transition occurs. By varying $N_*$, the signal can be shifted into frequency bands accessible to future detectors (e.g., LISA, DECIGO, BBO).
• Redshifted Signal: The authors calculated the redshifted GW spectrum today, showing it falls within the sensitivity curves of proposed space-based interferometers and pulsar timing arrays, depending on the specific timing of the transition.

5. Significance

• Feasibility of FoPT Exit: The work demonstrates that ending inflation via vacuum decay is theoretically feasible if the potential barrier is dynamically controlled, overcoming the historical "graceful exit" problem.
• Observational Probe: The distinctive oscillatory feature in the high-frequency GW spectrum serves as a unique "smoking gun" signature. It distinguishes inflation-triggered FoPTs from standard thermal phase transitions (which typically produce smoother spectra).
• Future Detection: The results provide concrete targets for next-generation gravitational wave observatories. Detecting these oscillations would not only confirm the existence of first-order phase transitions in the early universe but also validate specific high-energy physics models (like GUT-scale SUGRA constructions) and the Starobinsky inflation scenario.
• Methodological Advance: The successful application of 3D lattice simulations to this problem sets a new standard for studying non-linear dynamics in early universe cosmology, paving the way for investigating related phenomena like primordial black hole formation and isocurvature perturbations.