Lattice simulations of scalar-induced gravitational waves from inflation

This paper employs lattice simulations to demonstrate that nonperturbative scalar dynamics during ultra-slow-roll inflation can significantly alter the amplitude and spectral shape of scalar-induced gravitational waves, revealing that standard semi-analytical perturbative methods often fail to provide reliable predictions in regimes of large non-Gaussianity.

The Gist

Imagine the universe as a giant, invisible trampoline. When you jump on it, the fabric ripples. In the world of physics, these ripples are called gravitational waves.

For a long time, scientists have been looking for these waves from the very beginning of the universe (the Big Bang). But there's a problem: the standard tools we use to predict what these waves should look like are like using a ruler to measure a squiggly, wiggly noodle. They work okay for straight lines, but when things get messy and chaotic, the ruler fails.

This paper is about building a super-computer simulation (a "lattice") to measure that messy noodle directly, rather than guessing with a ruler.

Here is the story of what they found, broken down into simple parts:

1. The Setup: The "Ultra-Slow" Roller Coaster

Imagine the early universe was a roller coaster. Usually, the cart (the "inflaton" field) rolls down the track at a steady, slow pace. This creates the smooth, flat universe we see today.

But in some scenarios, the cart hits a section of track that is almost perfectly flat. It gets stuck there, moving incredibly slowly. This is called Ultra-Slow-Roll (USR).

• The Effect: Because it's moving so slowly, tiny bumps on the track get amplified massively. These bumps are "ripples" in the fabric of space.
• The Result: When the cart finally speeds up again, those amplified bumps are huge. They are so big that they crash into each other, creating a second set of ripples: Gravitational Waves.

2. The Old Way: The "Ruler" (Semi-Analytical Prediction)

For years, scientists predicted these waves using math that assumes the bumps on the track are small and behave nicely (like a calm ocean). They used a formula that treats the universe as if it were made of perfect, smooth waves.

• The Problem: When the bumps get huge (which happens in the "Ultra-Slow" scenario), the ocean isn't calm anymore; it's a tsunami. The "ruler" math breaks down because it doesn't account for the waves crashing into each other (non-linearity) or the weird, lumpy shapes they form (non-Gaussianity).

3. The New Way: The "Video Game" (Lattice Simulation)

Instead of using a ruler, the authors built a 3D video game of the early universe.

• They didn't assume the waves were smooth. They let the computer calculate exactly how the field moved, bumped, and crashed, step-by-step.
• They then took those messy, realistic results and simulated how they would create gravitational waves as the universe expanded.

4. The Big Discovery: The "Trapped" Patches

When they ran their simulation, they found two main things:

A. For "Mild" Chaos:
If the bumps were just a little bit messy, the old "ruler" math was mostly right. It got the general size of the signal correct, but it missed some details, especially the high-pitched parts of the sound.

B. For "Wild" Chaos (The Big Surprise):
When the chaos was extreme, the old math completely failed.

• The Analogy: Imagine trying to predict the sound of a rock concert by listening to a single violin. The old math is the violin. The reality is the whole band screaming, drums banging, and feedback loops.
• The "Trapping" Phenomenon: In the wildest scenarios, parts of the universe got "stuck" in a local dip in the energy track (like a ball getting stuck in a small hole on a hill).
  • Some patches of the universe kept rolling down the hill normally.
  • Other patches got stuck in the hole, then rolled back up, then got stuck again.
  • This created a patchwork quilt of the universe where different regions were behaving totally differently.
• The Result: This created a gravitational wave signal that looked nothing like the old predictions. The shape of the wave was different, and the volume was much louder (or quieter) than expected.

5. Why This Matters

The authors are saying: "If we want to understand the universe's past, we can't just use simple math anymore."

• The Stakes: We are currently detecting gravitational waves from pulsars (dead stars spinning in space). Scientists are trying to figure out if these waves are from black holes colliding or from the Big Bang itself.
• The Warning: If the Big Bang waves came from these "Ultra-Slow" scenarios, the old math might lead us to the wrong conclusion. We might think we see a signal from one thing, when it's actually a distorted signal from something else entirely.

The Bottom Line

This paper is like upgrading from a hand-drawn map to a GPS with real-time traffic.

• The old map (semi-analytical math) said, "The road is straight, drive 5 miles."
• The new GPS (lattice simulation) says, "The road is actually a chaotic, winding mountain pass with traffic jams and detours. The destination is in a totally different spot than you thought."

They have made their "GPS code" public, so other scientists can use it to navigate the messy, exciting early universe without getting lost.

Technical Summary

1. Problem Statement

Scalar-induced gravitational waves (SIGWs) are tensor modes generated at second order by amplified scalar perturbations during inflation. They serve as a crucial probe for inflationary dynamics on scales far smaller than those accessible to the Cosmic Microwave Background (CMB) and Large Scale Structure (LSS).

The Core Issue:
Standard predictions for SIGWs rely on semi-analytical calculations based on:

1. Linear perturbation theory for the evolution of scalar fluctuations during inflation.
2. Gaussian statistics for the primordial curvature perturbations.
3. Linear evolution of the Newtonian potential during the post-reheating era.

These approximations break down in scenarios where the curvature power spectrum is strongly enhanced on small scales (e.g., to generate Primordial Black Holes or explain PTA signals). In such regimes, specifically those involving a transient Ultra-Slow-Roll (USR) phase, scalar dynamics become highly nonlinear, and primordial non-Gaussianities (NG) can become large. The validity of truncating the non-Gaussian expansion or assuming linear evolution is uncontrolled in these regimes.

2. Methodology

The authors employ a first-principles lattice simulation framework to compute SIGWs without relying on the standard Gaussian or linear approximations. The pipeline consists of two distinct stages:

Stage I: Nonlinear Inflationary Dynamics

• Simulation: They perform a fully nonlinear lattice simulation of a single-field inflation model featuring a transition from Slow-Roll (SR) $\to$ USR $\to$ SR.
• Field Evolution: The inflaton field $\phi$ is evolved using the full nonlinear Klein-Gordon and Friedmann equations.
• Curvature Extraction: Instead of using linear perturbation theory, they extract the super-horizon curvature perturbation $\zeta(\mathbf{x})$ nonperturbatively using the $\delta N$ formalism applied to the lattice configuration. This captures intrinsic non-Gaussianity generated by field dynamics and the nonlinear mapping between $\delta\phi$ and $\zeta$.

Stage II: Post-Reheating GW Generation

• Initial Conditions: The frozen super-horizon $\zeta(\mathbf{x})$ from Stage I is mapped to the Newtonian potential $\Phi$ (via the linear relation $\Phi = \frac{2}{3}\zeta$) to set initial conditions for the post-reheating era.
• Evolution: The Newtonian potential $\Phi$ is evolved linearly in a radiation-dominated background using the Bardeen equation. Crucially, the source term for the gravitational waves (the scalar stress tensor) retains the full non-Gaussian structure of the initial $\Phi$ field.
• GW Calculation: The tensor modes $h_{ij}$ are evolved using the linear wave equation sourced by the nonlinear scalar stress. The Transverse-Traceless (TT) projection is applied only at the end to compute the power spectrum.

3. Key Contributions

1. First-Principles Benchmark: This work provides the first nonperturbative calculation of SIGWs from USR inflation, serving as a benchmark to test the validity of standard semi-analytical methods.
2. Separation of Effects: The methodology allows the authors to disentangle two distinct sources of deviation from standard predictions:
   • Backreaction/Nonlinear Dynamics: Modifications to the curvature power spectrum ($P_\zeta$) itself due to nonlinear field evolution.
   • Genuine Non-Gaussianity: The impact of non-Gaussian statistics on the GW source term, even if $P_\zeta$ were fixed.
3. Code Availability: The authors release the simulation code (InflationEasy) and data, enabling the community to reproduce results and extend the framework.

4. Key Results

The authors analyze two categories of models: Mild Non-Gaussianity and Large Non-Gaussianity.

A. Mild Non-Gaussianity

• Scenario: Models where the USR phase is smooth, and non-Gaussianities are moderate (e.g., Wands-dual regimes).
• Findings:
  • The standard semi-analytical prediction captures the correct order of magnitude of the GW signal.
  • Corrections: Deviations appear at larger amplitudes ($P_\zeta \gtrsim 10^{-3}$).
  • UV Sensitivity: Non-Gaussian effects predominantly affect the high-wavenumber (UV) part of the spectrum.
  • Dominant Effect: In these models, the inflationary backreaction (which modifies $P_\zeta$) often causes a larger shift in the signal than the non-Gaussianity itself.

B. Large Non-Gaussianity

• Scenario: Models with sharp transitions and large $\eta$ parameters, leading to strong intrinsic non-Gaussianity and "trapping" phenomena (where field patches get stuck in local minima of the potential).
• Findings:
  • Catastrophic Failure of Standard Methods: The semi-analytical prediction fails dramatically in both amplitude and spectral shape, regardless of the overall size of the tensor power spectrum.
  • Amplitude Enhancement: Primordial non-Gaussianity enhances the GW spectrum by approximately one order of magnitude across all wavenumbers compared to Gaussian predictions.
  • Spectral Shape: The standard method fails to reproduce the shape of the spectrum, particularly in the UV, where "trapping" creates a plateau-like structure.
  • Inflationary Contribution: Unlike mild cases, the inflationary contribution to GWs (emitted at horizon exit) becomes resolvable and non-negligible due to enhanced scalar gradients.

C. Trapping and Peak Structure

• In the large-NG regime, the authors observe a multi-peak structure in the probability density function (PDF) of the inflaton field.
• This arises from "trapped" regions (red in their snapshots) where the field rolls back toward a local minimum, surrounded by non-trapped shells. This creates a complex, oscillatory field configuration that standard perturbation theory cannot capture.

5. Significance and Implications

• Reliability of Predictions: The study demonstrates that for scenarios with large curvature enhancements (relevant for PBH formation and PTA signals), nonperturbative control of inflationary scalar dynamics is mandatory. Standard Gaussian approximations are insufficient.
• Interpretation of Observations: The results suggest that interpreting potential SIGW detections (e.g., by NANOGrav or future LISA) as evidence for specific inflationary models requires theoretical predictions that account for nonlinear effects. Ignoring them could lead to incorrect conclusions about the inflationary potential or the abundance of Primordial Black Holes.
• Future Directions: The authors highlight the need to extend the framework to:
  • Nonlinear evolution of the gravitational potential in the radiation era (when $\Phi \sim 1$).
  • Higher-order tensor sources beyond second order.
  • Fully hydrodynamical descriptions of density perturbations.

In conclusion, this paper establishes that while semi-analytical methods work for mild scenarios, they are fundamentally inadequate for the highly nonlinear regimes where the most interesting cosmological phenomena (like PBH formation) occur. The lattice approach provides the necessary rigorous toolset to bridge this gap.