Numerical simulations of density perturbation and gravitational wave production from cosmological first-order phase transition

Through three-dimensional lattice simulations, this study reveals that the dominant source of density perturbations in first-order phase transitions depends on the transition strength ($\alpha$), characterizes the resulting power spectra slopes for both density perturbations and gravitational waves, and confirms that slow transitions can generate primordial black holes while offering theoretical predictions for gravitational wave detection.

The Gist

The Big Picture: A Cosmic Boiling Pot

Imagine the very early universe not as a calm, empty void, but as a giant pot of water sitting on a stove. As the universe cools down (like the water losing heat), it goes through a First-Order Phase Transition (FOPT).

Think of this like water freezing into ice. But instead of freezing smoothly, it happens in a chaotic, bubbly way. Tiny pockets of "new reality" (let's call them Ice Bubbles) suddenly appear in the "old reality" (the Hot Water). These bubbles expand rapidly, slam into each other, and eventually take over the whole pot.

This paper is a computer simulation of that chaotic boiling process. The scientists wanted to know two things:

1. Did this chaos create "lumps" heavy enough to collapse into Black Holes? (These are called Primordial Black Holes, or PBHs).
2. Did this chaos create "ripples" in space-time? (These are Gravitational Waves, or GWs).

---

The Two Main Characters: The Bubbles and the Delay

The simulation revealed that the outcome depends heavily on how "strong" the phase transition is. We can think of this strength as the temperature difference between the water and the ice.

1. The "Strong" Transition (The Race)

• The Scenario: The transition is very energetic. The bubbles are like race cars. They nucleate (start) and zoom forward incredibly fast.
• The Result: The main source of trouble here is the bubble walls themselves. As these high-speed walls crash into each other, they create massive shockwaves.
• The Metaphor: Imagine a crowd of people running so fast that when they bump into each other, they create a massive pile-up. The energy of the movement is what causes the problem.

2. The "Weak" Transition (The Wait)

• The Scenario: The transition is sluggish. The bubbles are like snails. They start appearing, but they take a long time to grow and merge.
• The Result: The main source of trouble here is delay. Because the bubbles are slow, some pockets of the "old reality" (the hot water) get left behind for a long time.
• The Metaphor: Imagine a party where everyone is supposed to leave at 10:00 PM. If everyone leaves instantly, the house is empty. But if a few people are slow to leave and stay in the living room while the rest of the house is empty, that small group becomes a "dense" crowd in an empty house. This delay creates a heavy, dense spot that can collapse.

Key Finding: The paper found that if the transition is slow (low speed of bubble nucleation), these "delayed" pockets get so heavy they collapse into Primordial Black Holes. If the transition is too fast, everything happens too evenly, and no black holes form.

---

The Ripples: Gravitational Waves

When these bubbles collide, they don't just make black holes; they also shake the fabric of space-time, creating Gravitational Waves.

• The Analogy: Imagine dropping a stone in a pond. It makes ripples. Now imagine thousands of bubbles popping and crashing in a pond simultaneously. That creates a massive, chaotic wave pattern.
• The Pattern: The scientists mapped out the "music" of these waves (the power spectrum).
  • At low frequencies (long, slow ripples), the signal grows very fast (like a rising drumbeat).
  • At high frequencies (fast, sharp ripples), the signal drops off.
• Why it matters: This specific "song" is a fingerprint. If future detectors (like LISA or the Einstein Telescope) hear this specific pattern, they will know exactly what kind of phase transition happened in the early universe.

---

The "Goldilocks" Zone for Black Holes

The paper uses a parameter called $\beta/H$ to measure how fast the transition happens.

• Fast Transition (High $\beta$): Everything happens too quickly. The universe smooths out. No Black Holes.
• Slow Transition (Low $\beta$): The "delayed" pockets have time to get very heavy before the rest of the universe catches up. Black Holes form!

The simulation showed that for the slowest transitions they tested, the density became so high that it crossed the "collapse threshold," turning those pockets into black holes.

Why Should You Care?

1. Dark Matter Mystery: Primordial Black Holes are a leading candidate for Dark Matter (the invisible stuff holding galaxies together). If these bubbles created them, we might finally know what dark matter is.
2. Listening to the Big Bang: We can't see the very early universe with telescopes because it was too hot and opaque. But Gravitational Waves pass through everything. If we detect the specific "song" predicted by this paper, we will be able to "hear" the universe's first phase transition, proving theories about physics that are far beyond our current particle accelerators.

Summary in One Sentence

This paper simulates the early universe's "freezing" process and finds that if the process is slow enough, the resulting delays create heavy clumps that turn into black holes, while the crashing bubbles create a unique gravitational wave signal that future detectors might one day hear.

Technical Summary

1. Problem Statement

The paper addresses two critical phenomena in the early universe associated with First-Order Phase Transitions (FOPTs):

• Primordial Black Hole (PBH) Formation: Whether the density perturbations generated during FOPTs are sufficient to trigger the collapse of Hubble patches into PBHs.
• Gravitational Wave (GW) Production: The characteristics of the stochastic GW background generated by the collision of vacuum bubbles, sound waves, and turbulence during the transition.

While previous studies often relied on analytical approximations or specific fluid models, there was a need for rigorous 3D lattice simulations to understand the interplay between the randomness of quantum tunneling (vacuum decay delay) and the dynamics of bubble wall collisions, particularly regarding the resulting density power spectra and PBH abundance.

2. Methodology

The authors employed three-dimensional numerical lattice simulations to model the FOPT process in an expanding universe.

• Physical Model:

  • Scalar Field: A scalar field $\phi$ with a potential $V(\phi) = \frac{1}{2}M^2\phi^2 + \frac{1}{3}\kappa\phi^3 + \frac{1}{4}\lambda\phi^4$.
  • Nucleation: Vacuum bubbles nucleate via thermal tunneling with an exponential rate $\Gamma(t) = \Gamma_0 e^{\beta t}$.
  • Dynamics: The field evolution is governed by the Klein-Gordon equation in an expanding universe, coupled with the Friedmann equations.
  • Energy Components: The total energy density is decomposed into kinetic ($\rho_K$), gradient ($\rho_G$), potential ($\rho_V$), and radiation ($\rho_r$) components.

• Simulation Setup:

  • Lattice: A $512^3$ grid covering approximately 43 Hubble volumes ($L_H = 4$).
  • Parameters: The study explored a wide parameter space:
    • Phase transition strength: $\alpha = 0.5, 1, 5, 10$.
    • Inverse duration: $\beta/H = 6, 8, 10, 12$.
    • Potential parameters were tuned to ensure successful FOPT with specific bubble wall thicknesses.
  • PBH Criteria: The authors identified "delayed decay" regions (the last 25% of Hubble volumes to nucleate) as potential overdense regions. They calculated the density contrast $\delta = (\rho_{in} - \rho_{out})/\rho_{out}$ and compared it against a critical threshold $\delta_c = 0.45$ (optimistically chosen to account for gauge inconsistencies in literature) to determine PBH formation.

• Analysis:

  • Calculated the variance of density perturbations ($\sigma_\delta$) and their power spectra ($P_\delta$).
  • Computed the GW power spectrum ($\Omega_{gw}$) by solving the linearized Einstein equations for tensor metric perturbations driven by the transverse-traceless component of the stress-energy tensor.

3. Key Contributions

• Mechanism Differentiation: The study explicitly distinguishes between two primary sources of density perturbations:
  1. Delay of Vacuum Decay: Dominant in weak transitions ($\alpha < 1$), where asynchronous nucleation creates local overdensities.
  2. Bubble Wall Motion: Dominant in strong transitions ($\alpha > 1$), where the kinetic energy of expanding walls drives perturbations.
• Spectral Slopes: The paper provides precise numerical determination of the power spectrum slopes for both density perturbations and GWs across different wavenumber regimes, offering benchmarks for future observational data.
• PBH Formation Conditions: It establishes a clear relationship between the phase transition parameters ($\alpha, \beta/H$) and the probability of PBH formation, highlighting the critical role of "slow" transitions.

4. Key Results

A. Primordial Black Hole (PBH) Abundance

• Dependence on Parameters: PBH abundance increases with the transition strength $\alpha$ but decreases significantly with the inverse duration $\beta/H$.
• Critical Threshold: PBHs are successfully formed in simulations where $\beta/H = 6$ (slow transitions). As $\beta/H$ increases (faster transitions), the probability of exceeding the critical density threshold $\delta_c$ drops sharply.
• Saturation: For very strong transitions ($\alpha \gtrsim 10$), the influence of $\alpha$ on PBH abundance becomes suppressed.

B. Density Perturbation Power Spectrum ($P_\delta$)

• Slopes:
  • Small wavenumbers ($k < k_H$): Slope follows $k^3$. This is attributed to the uncorrelated, Poisson-like nature of vacuum decay in the infrared.
  • Large wavenumbers ($k > k_H$): Slope follows $k^{-1.5}$.
    • Analytical Note: A pure step-function energy distribution would yield $k^{-1}$. The observed $k^{-1.5}$ is attributed to the smoothing effect of kinetic and gradient energies across the bubble walls in the simulation.
• Peak Behavior: Two peaks are observed in the evolution of $\sigma_\delta$:
  • For $\alpha < 1$: The peak is driven by the delay of vacuum decay.
  • For $\alpha > 1$: The peak is driven by the forward motion of bubble walls.

C. Gravitational Wave (GW) Power Spectrum ($\Omega_{gw}$)

• Slopes:
  • Small wavenumbers: Slope follows $k^3$ (causality limit).
  • Large wavenumbers: Slope follows $k^{-2}$.
• Parameter Sensitivity: The overall shape and amplitude of the GW spectrum show little sensitivity to variations in $\alpha$ and $\beta/H$, unlike the density perturbations.
• Detectability: The authors fitted the spectra to standard parametrizations and demonstrated that specific scenarios (e.g., $\alpha=1, \beta/H=10$ at $T \sim 10^{-2}$ GeV or $10^3$ GeV) could be detectable by future missions like LISA, Taiji, SKA, and DECIGO, as well as current PTA experiments (NANOGrav, EPTA).

5. Significance and Conclusion

• Theoretical Support for GW Detection: The paper provides robust theoretical predictions for the GW power spectrum slopes ($k^3$ and $k^{-2}$) arising from FOPTs, aiding in the interpretation of potential signals from current and future detectors.
• PBH Formation Mechanism: It confirms that slow phase transitions ($\beta/H \lesssim 6$) are a viable mechanism for generating PBHs, potentially explaining dark matter or high-energy phenomena.
• Model Independence: By using a model-independent lattice approach, the study offers a benchmark for various extensions of the Standard Model that predict FOPTs.
• Future Outlook: The authors note that while their results are reliable for density perturbations, fully relativistic simulations including fluid dynamics (sound waves) are needed to refine predictions for $\alpha > 1$ scenarios where bubble wall energy converts to fluid motion.

In summary, this work bridges the gap between analytical estimates and complex numerical realities, demonstrating that the stochastic nature of vacuum decay in FOPTs is a potent source of both PBHs and a distinct gravitational wave signature.