STOchastic LAttice Simulation of hybrid inflation

Employing the STOchastic LAttice Simulation (STOLAS) code, this study investigates the spatial profiles and topological structures of curvature perturbations in multi-waterfall hybrid inflation models, validating the stochastic-$\delta$N algorithm while revealing how stochastic noise drives defect reconnection and imposes upper bounds on perturbations in cubic potential scenarios.

The Gist

The Cosmic Growth Spurt: A Simulation of the Early Universe

Imagine the universe right after the Big Bang. It didn't just expand; it went through a massive "growth spurt" called Inflation. In a tiny fraction of a second, it stretched out faster than the speed of light. This paper is about a specific type of inflation called Hybrid Inflation, and the authors used a super-computer to simulate exactly how it worked.

Think of this research as a "Cosmic Weather Report" for the very first moments of time.

1. The Setup: A Ball and a Trapdoor

To understand their model, imagine a ball rolling down a long, flat hill.

• The Inflaton (The Ball): This is the main energy driving the expansion. It rolls slowly along the hill.
• The Waterfall Fields (The Trapdoor): At a certain point on the hill, there is a trapdoor. As long as the ball is far away, the door is locked. But once the ball reaches a specific "critical point," the door unlocks.
• The Waterfall: Once the door opens, the ball (and the fields around it) suddenly falls down a steep cliff. This rapid fall marks the end of inflation.

The authors wanted to know: What happens to the "ripples" in space-time during this fall? These ripples are important because they eventually grow into galaxies, stars, and us.

2. The Tool: STOLAS (The Cosmic Video Game)

To study this, the authors used a code called STOLAS.

• The Grid: Imagine a giant 3D video game world made of pixels (or voxels).
• The Simulation: They didn't just watch one ball roll. They simulated millions of tiny points in space all at once.
• The Noise: In the real quantum world, things aren't perfectly smooth. There is "static" or "noise" (like rain falling randomly). STOLAS adds this random rain to the simulation to see how it messes up the smooth path of the ball.

3. The Findings: What Did They Discover?

The team ran six different versions of this simulation, changing two things:

1. How many "trapdoors" (waterfall fields) were there? (They tested 1, 2, 3, and even 15).
2. What shape was the hill? (They tested a "Quadratic" shape vs. a "Cubic" shape).

Here is what they found:

A. The "Speed Limit" on Black Holes

In some versions of the simulation (the "Cubic" cases), they found a ceiling on how big the ripples could get.

• Analogy: Imagine filling a bucket with water. In the "Quadratic" case, you can keep pouring until the bucket overflows. In the "Cubic" case, the bucket has a lid. No matter how hard you pour, the water level can't go above the lid.
• Why it matters: If the ripples get too big, they collapse into Primordial Black Holes. Because the "Cubic" model has a lid (an upper bound), it prevents the ripples from getting huge enough to form these black holes. This suggests that in this specific universe scenario, we wouldn't see as many ancient black holes as we might otherwise expect.

B. The "Knots" in the Fabric

When the trapdoor opens, the fields can get tangled, creating "Topological Defects."

• Analogy: Think of a bowl of spaghetti. If you pull the strands apart quickly, they might snap into knots, loops, or long strands.
  • 1 Field: Creates a "Wall" (like a sheet of paper).
  • 2 Fields: Creates a "String" (like a thread).
  • 3 Fields: Creates a "Monopole" (like a ball of yarn).
• The Twist: The authors expected these knots to stay big. However, the random "noise" (the quantum rain) acted like a pair of scissors. It cut the big knots into tiny, fine fuzz. By the time inflation ended, the big structures were shredded into microscopic pieces.

C. The Fingerprint of the Universe

The authors looked at the final map of the universe's ripples to see if these knots left a mark.

• The Result: Only the case with 1 waterfall field (the "Wall") left a distinct, global pattern on the map. The others (Strings and Monopoles) were too messy or too small to leave a clear fingerprint on the curvature of space.
• Why it matters: If we look at the large-scale structure of the universe today, finding this specific pattern could tell us exactly what kind of inflation happened back then.

4. Why Should We Care?

This paper is like checking the blueprint of a house before it's built.

1. It validates the math: They proved that their computer simulation matches the theoretical math used by other scientists.
2. It explains the Black Hole mystery: It helps explain why we might not see as many ancient black holes as some theories predict.
3. It maps the future: By understanding how these ripples form, we can better understand how galaxies formed and why the universe looks the way it does today.

In a nutshell: The authors built a 3D computer model of the universe's birth. They found that random quantum noise shreds cosmic "knots" into tiny pieces and that some models put a "speed limit" on how big cosmic ripples can get, which controls how many ancient black holes are made.

Technical Summary

Here is a detailed technical summary of the paper "STOchastic LAttice Simulation of hybrid inflation" (RUP-26-3).

1. Problem Statement

Hybrid inflation is a prominent cosmological model involving multiple scalar fields (an inflaton and waterfall fields) that can generate significant curvature perturbations on small scales. These perturbations are candidates for the formation of Primordial Black Holes (PBHs) and the generation of topological defects (domain walls, cosmic strings, monopoles). However, the dynamics of the "waterfall phase" (where inflation ends) are intrinsically non-perturbative due to tachyonic instability and field fluctuations.

• Challenge 1: Standard linear perturbation theory breaks down during the waterfall transition.
• Challenge 2: While the stochastic-$\delta N$ formalism provides a non-perturbative framework for calculating curvature perturbations, it relies on spatially one-point statistics. It lacks direct spatial resolution to analyze the correlation structures and topological configurations of the perturbations.
• Challenge 3: The impact of topological defects formed during inflation on the spatial profile of the curvature perturbation remains unclear.

2. Methodology

The authors utilize a lattice simulation code called STOLAS (STOchastic LAttice Simulation) to investigate the spatial profile of curvature perturbations in mild-waterfall hybrid inflation.

• Simulation Framework (STOLAS):
  • Based on the stochastic formalism of inflation, where long-wavelength (IR) modes are treated as classical fields subject to stochastic noise from short-wavelength (UV) modes crossing the Hubble horizon.
  • Implements discrete Langevin-type equations for coarse-grained fields on a 3D lattice ($N_L = 256$ grids).
  • Combines with the $\delta N$ formalism, calculating the curvature perturbation $\zeta$ as the fluctuation in the e-folding number ($\delta N$) from an initial flat hypersurface to a final uniform-density hypersurface.
• Model Configuration:
  • Inflation Model: Mild-waterfall hybrid inflation with one inflaton field ($\phi$) and $n$ waterfall fields ($\psi$).
  • Potentials: Two functional forms were tested:
    • "Quadratic": Standard expansion where higher-order terms are negligible.
    • "Cubic": Includes a relevant cubic term in the inflaton potential.
  • Cases: Six scenarios were simulated varying the number of waterfall fields ($n = 1, 2, 3, 15$ for Quadratic; $n = 1, 2$ for Cubic).
• Analysis Tools:
  • Statistical: Probability Density Functions (PDFs) and Power Spectra (via Fourier transform of lattice maps).
  • Topological: Euler characteristic ($\chi$) calculated for both the waterfall fields and the final curvature perturbation maps to identify global structures (holes, loops, isolated objects).

3. Key Contributions

1. Validation of Stochastic-$\delta N$: The paper provides the first direct comparison between the stochastic-$\delta N$ algorithm (analytical fitting formula) and full 3D lattice simulations (STOLAS) for hybrid inflation, validating the algorithm's power spectrum predictions.
2. Discovery of PDF Upper Bound: The authors confirmed the existence of a characteristic upper bound in the PDF for the "Cubic" model, a feature previously identified in point-statistics but now verified spatially.
3. Topological Diagnostics: The application of the Euler characteristic to the curvature perturbation field serves as a novel probe to detect imprints of topological defects on large-scale structure.
4. Coarse-Graining Clarification: The study clarifies discrepancies between lattice results and previous stochastic calculations by highlighting the importance of the coarse-graining scale (lattice grid size vs. Hubble scale).

4. Key Results

• Power Spectrum:
  • The power spectra derived from STOLAS are broadly consistent with the fitting formula derived from the stochastic-$\delta N$ formalism.
  • For "Quadratic" models, the peak amplitude scales roughly as $1/n$.
  • For "Cubic" models, this scaling relation does not hold because the amplitude is suppressed by the upper bound on $\delta N$, which is independent of $n$.
• Probability Density Function (PDF):
  • "Quadratic" Cases: Show heavy tails consistent with previous non-perturbative studies, though discrepancies at $\delta N \gtrsim 1$ arise due to the larger coarse-graining scale of the lattice compared to the Hubble scale used in previous works.
  • "Cubic" Cases: Exhibit a strict upper bound ($\delta N \lesssim 0.2$). This suppresses the formation of large positive curvature perturbations, implying that PBHs are seldom formed in this specific model.
• Topological Defects:
  • Snapshots of waterfall fields show that topological defects (domain walls for $n=1$, strings for $n=2$, monopoles for $n=3$) form near the critical point.
  • Stochastic Reconnection: Contrary to the expectation that defects form on Hubble scales, stochastic noise causes them to reconnect into much finer structures. Their correlation lengths become significantly smaller than the Hubble scale at the critical point.
• Curvature Perturbation Topology:
  • Using the Euler characteristic, the authors found that only the $n=1$ (domain wall) case leaves a non-trivial topological signature (negative Euler characteristic) in the curvature perturbation map.
  • Cases with $n=2$ and $n=3$ do not leave such global structures in the curvature perturbation, despite the presence of defects in the waterfall fields.

5. Significance and Implications

• PBH Formation Constraints: The "Cubic" model's upper bound on curvature perturbations suggests that specific hybrid inflation potentials can naturally suppress PBH formation, distinguishing them from "Quadratic" models which allow for heavier tails.
• Early Universe Probes: The finding that the $n=1$ case imprints a unique topological structure on the curvature perturbation suggests that the global structure of the curvature field could serve as a novel probe for early universe physics and the nature of inflationary fields.
• Future Directions: The authors propose that the spatial maps generated by STOLAS can be used to calculate higher-order correlators (bispectra) and investigate scalar-induced gravitational waves. Furthermore, the specific configurations of the "Cubic" model (bunches of regions near the maximum $\zeta$) warrant investigation via N-body simulations to understand their impact on halo formation.

In summary, this work bridges the gap between analytical stochastic formalisms and spatially resolved lattice simulations, providing robust validation for the stochastic-$\delta N$ approach while uncovering new topological and statistical features of hybrid inflation that impact PBH formation and early universe structure.