⚛️ phenomenology

How large are curvature perturbations from slow first-order phase transitions? A gauge-invariant analysis

This paper employs a gauge-invariant multi-fluid formalism to demonstrate that super-horizon inhomogeneities from slow, strongly supercooled first-order phase transitions are unlikely to produce Primordial Black Holes, while providing a fitting formula for the resulting curvature perturbations and discussing their observational constraints via primordial curvature limits and scalar-induced gravitational waves.

Original authors: Xiao Wang, Csaba Balázs, Ran Ding, Chi Tian

Published 2026-01-22

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Xiao Wang, Csaba Balázs, Ran Ding, Chi Tian

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the early universe as a giant pot of water cooling down. Usually, water freezes smoothly into ice. But in the world of particle physics, sometimes the universe gets "stuck" in a hot, liquid state (a "false vacuum") even when it's cold enough to freeze. Eventually, it snaps into the solid state (the "true vacuum") all at once. This is called a First-Order Phase Transition (FOPT).

Think of this like a pot of water suddenly forming ice bubbles. Usually, these bubbles form quickly and everywhere at once. But this paper asks: What happens if the bubbles form very slowly and very unevenly?

Here is a breakdown of what the researchers found, using simple analogies:

1. The "Slow Freeze" Problem

If the universe cools down too fast, the bubbles form quickly. But if the transition is strongly supercooled (very cold but still liquid) and slow, the bubbles might take a long time to appear.

The Analogy: Imagine a giant crowd of people trying to form a circle. If they all start at the same time, the circle forms evenly. But if they start at random times, some areas will have a tight circle while others are still empty.
The Result: Because the bubbles form at random times in different parts of the universe, some regions end up with more energy (more "ice") and some with less. This creates "lumps" or inhomogeneities in the universe's energy.

2. The Measurement Mistake (The Gauge Issue)

Scientists have been trying to measure how big these "lumps" are. Previous studies used a method called "separate universe simulations."

The Analogy: Imagine trying to measure the height of a wave in a stormy ocean. If you measure the wave while standing on a boat that is bobbing up and down (a specific "gauge"), you might think the wave is huge. But if you measure from a fixed point in space, the wave might look much smaller.
The Paper's Fix: The authors realized that previous studies were measuring these lumps from a "bobbing boat." They developed a new, gauge-invariant method (like measuring from a fixed satellite) to get the true size of the ripples. They found that the "lumps" are actually much smaller than people thought.

3. Do These Lumps Create Black Holes?

A big question in physics is whether these energy lumps are big enough to collapse into Primordial Black Holes (PBHs)—tiny black holes formed right after the Big Bang.

The Old View: Previous calculations suggested the lumps were so huge that they would easily crush themselves into black holes.
The New View: Using their new, more accurate measurement, the authors found the lumps are too small.
The Verdict: It is highly unlikely that these slow phase transitions created primordial black holes. The "lumps" aren't heavy enough to collapse.

4. Do They Create Gravitational Waves?

When these energy lumps eventually smooth out, they can create ripples in space-time called Gravitational Waves (GWs).

The "Primary" Waves: These come from the violent collision of the bubbles themselves (like two ice chunks crashing).
The "Secondary" Waves: These come from the "lumps" of energy smoothing out later (like the ripples left after the ice chunks settle).
The Finding: The authors calculated that these secondary waves are very weak. While they exist, they are so quiet that they don't really change what we see in current data from Pulsar Timing Arrays (which listen for gravitational waves). They are like a whisper in a loud concert; you can't hear them over the main music.

Summary

The paper essentially says:

Slow, uneven freezing in the early universe creates energy lumps.
Old measurements overestimated how big these lumps were because of a mathematical "perspective" error.
New measurements show the lumps are too small to create black holes.
The ripples (gravitational waves) these lumps make are too faint to significantly alter our current understanding of the universe's history.

In short: The universe might have had a slow, bumpy freeze, but it wasn't bumpy enough to make black holes or loud enough to change the gravitational wave signals we are currently detecting.

Technical Summary: Curvature Perturbations from Slow First-Order Phase Transitions

Problem Statement
Cosmological first-order phase transitions (FOPTs) that are strongly supercooled and sufficiently slow (characterized by a small inverse duration parameter $\beta/H_n \sim 10$ ) are favored by recent Pulsar Timing Array (PTA) observations as a source of the stochastic gravitational-wave background (SGWB). Such transitions can generate distinct nucleation histories across causally disconnected regions, leading to super-horizon inhomogeneities in the total energy density. Previous studies, relying on "separate universe" simulations, suggested these inhomogeneities could induce large curvature perturbations, potentially exceeding the threshold for Primordial Black Hole (PBH) formation and generating significant scalar-induced gravitational waves (SIGWs). However, these earlier calculations suffered from gauge ambiguities, as the density contrasts derived from simulations ( $\delta_{NG}$ ) were often implicitly identified with gauge-dependent quantities in gravitational perturbation theory without a rigorous framework.

Methodology
To resolve these ambiguities, the authors employ a gauge-invariant, multi-fluid formalism. The system is modeled as two interacting fluids: radiation and vacuum energy.

Simulation of Nucleation Histories: Using the DeltaPT toolkit, the authors simulate $50,000$ independent realizations of bubble nucleation histories for various super-horizon scales ( $k > k_{max}$ ) and nucleation rates ( $\beta/H_n$ ). This yields the time evolution of the false vacuum fraction $F(t)$ and the associated density contrasts in a "no-gauge" ( $\delta_{NG}$ ) framework.
Gauge-Invariant Calculation: Instead of directly mapping $\delta_{NG}$ to a specific gauge, the authors assume $\delta_{NG}$ approximates density contrasts on an unknown time-slicing at linear order. They then utilize the established formalism for multiple interacting fluids to compute gauge-invariant quantities, specifically the comoving curvature perturbation $R$ and the total curvature perturbation $\zeta$ .
Evolution Equations: The time evolution of $\zeta$ is derived from the perturbed continuity equations, accounting for the non-adiabatic pressure perturbation $\delta P_{nad}$ . This includes intrinsic entropy perturbations (which vanish for fluids with definite equations of state) and relative entropy perturbations between the radiation and vacuum components.
Fitting and Constraints: The authors derive a fitting formula for the curvature power spectrum $P_R(k)$ based on the simulation results. They subsequently analyze observational constraints on $P_R(k)$ from CMB anisotropies, Lyman- $\alpha$ forests, spectral distortions, and Ultra-Compact Mini-Halos (UCMHs). Finally, they calculate the resulting SIGWs and assess their impact on the interpretation of PTA data for FOPT models.

Key Results

Breakdown of Simple Relations: The study demonstrates that there is no definite relation between the simulation-derived density contrast $\delta_{NG}$ and the gauge-invariant comoving curvature perturbation $R$ . The approximate relation $R \approx -(9/4)\delta_C$ (valid for radiation fluid perturbations) breaks down in this context, indicating that $\delta_{NG}$ cannot be simply identified with the comoving gauge density contrast.
Suppressed Curvature Perturbations: The variance of the resulting comoving curvature perturbations ( $\sigma_R^2$ ) is found to be significantly smaller than estimates in previous literature (e.g., Ref. [24]) and consistent with the more conservative bounds of Ref. [20]. The variance scales as $k^3$ at small $k$ and $\beta^{-5}$ at large $\beta$ .
PBH Formation: The distribution of the comoving density contrast $\delta_C$ , derived from the gauge-invariant $R$ , exhibits a positive tail that remains well below the critical threshold for PBH formation ( $\delta_{crit} \approx 0.45$ ). Consequently, the production of PBHs via super-horizon inhomogeneities from slow FOPTs is deemed highly unlikely.
Scalar-Induced Gravitational Waves (SIGWs): While SIGWs are generated when these modes re-enter the horizon, their amplitude is subdominant compared to the primary gravitational waves from bubble collisions or sound waves.
Impact on PTA Data: When the secondary SIGW component is included in the analysis of NANOGrav 15-year data, the posterior distributions for FOPT parameters (strength $\alpha$ , reheating temperature $T_{reh}$ , and nucleation rate $\beta/H_n$ ) show negligible modification for bubble collision models and only marginal changes for sound wave models.

Significance and Claims
The paper claims to provide a rigorous, gauge-invariant framework for evaluating the large-scale gravitational effects of FOPTs. Its primary significance lies in correcting previous overestimations of curvature perturbations and PBH production probabilities. By resolving gauge ambiguities, the authors conclude that:

Slow, strongly supercooled FOPTs are unlikely to produce PBHs via super-horizon inhomogeneities.
The associated SIGWs constitute a subdominant secondary source that does not significantly alter the interpretation of current PTA data regarding FOPT parameters.
The derived fitting formula for $P_R(k)$ serves as a robust template for future observational constraints, suggesting that future SKA data may constrain slow FOPTs, while UCMH observations constrain faster transitions.

The study does not propose new experimental proposals but rather refines the theoretical interpretation of existing and projected observational limits on primordial curvature perturbations and gravitational waves.