Gravitational waves from strong first order phase transitions

Through large-scale simulations of strong first-order phase transitions featuring both detonation and deflagration modes, the study reveals that despite differing flow dynamics and kinetic energy densities, both scenarios yield a similar gravitational wave production efficiency of approximately 0.017, with detailed analysis of velocity and shear stress decorrelation providing new insights into the underlying mechanisms.

The Gist

Imagine the early universe not as a smooth, calm ocean, but as a chaotic, boiling pot of energy. As the universe cooled down, it didn't just settle; it underwent a phase transition, much like water turning into ice. But instead of freezing slowly, it did so violently, creating bubbles of the new "ice" phase that expanded, collided, and crashed into each other.

This paper is about what happens when those bubbles crash together in the most extreme scenarios possible, and how those crashes create ripples in space-time known as gravitational waves.

Here is the story of the research, broken down into simple concepts:

1. The Great Bubble Battle

Think of the early universe as a giant room filled with fog. Suddenly, the temperature drops, and "ice bubbles" start forming.

• The Deflagration (The Slow Burn): In one scenario, the bubbles expand like a slow-moving fire. The air in front of the fire gets pushed ahead, creating a gentle compression wave. It's like a crowd of people slowly shuffling forward in a hallway.
• The Detonation (The Explosion): In the other scenario, the bubbles expand at nearly the speed of light. This is like a supersonic jet breaking the sound barrier. It creates a massive, violent shockwave that slams into the surrounding fluid. It's a sonic boom, not a gentle breeze.

The researchers simulated both of these scenarios using massive supercomputers to see what happens when these bubbles collide and merge.

2. The Sound of the Crash

When these bubbles collide, they don't just stop. They create a chaotic mess of moving fluid.

• The "Sound" of the Phase Transition: Just like clapping your hands creates sound waves in the air, these bubble collisions create sound waves in the cosmic fluid.
• The Turbulence: As these sound waves crash into each other, they create turbulence. Imagine a river hitting a rock; the water swirls into eddies (vortices). In the "slow burn" scenario, these swirling eddies are very strong. In the "explosion" scenario, the motion is dominated by the straight-line shockwaves.

3. The Ripples in Space-Time (Gravitational Waves)

Einstein told us that when massive things move and accelerate, they create ripples in the fabric of space-time. These are gravitational waves.

• The researchers wanted to know: How loud is the music? (How strong are the waves?)
• They found that even though the "slow burn" creates a lot of swirling eddies, the shockwaves from the "explosion" are actually the ones making the loudest noise.
• The Surprise: Even in the "slow burn" case where there are lots of swirling eddies, those eddies don't contribute much to the gravitational waves. It's like having a choir of people humming (the eddies) while a single person screams (the shockwave); the scream drowns out the hum when it comes to making ripples in space-time.

4. The "Decorrelation" Dance

One of the paper's biggest breakthroughs is studying how long these waves last before they lose their rhythm.

• Imagine a group of dancers moving in perfect sync. If they start moving randomly, the pattern breaks. This is called decorrelation.
• The researchers found that in the "explosion" scenario, the shockwaves move so fast (faster than the speed of sound in that fluid) that they break the rhythm very quickly.
• They discovered that the "swirling" eddies in the slow scenario get "swept" back and forth by the sound waves, like leaves caught in a gust of wind, which changes how long they keep their pattern.

5. The Big Takeaway: A Universal Efficiency

The most exciting result is a number the researchers found.

• They calculated how much of the energy from the bubble collisions turns into gravitational waves.
• Despite the two scenarios being very different (one is a slow burn, the other a supersonic explosion), they both turned out to be roughly 1.7% efficient.
• The Analogy: Imagine two different engines: a slow, rumbling diesel truck and a high-speed Formula 1 car. You might expect the F1 car to be much more efficient at converting fuel into speed. But the researchers found that both engines convert about the same tiny fraction of their energy into the specific "sound" of gravitational waves.

Why Does This Matter?

We are currently building new "ears" for the universe, like the LISA space telescope (a giant laser interferometer in space). These telescopes are designed to listen for these ancient ripples.

• If we hear a signal, we need to know what it sounds like to identify it.
• This paper provides the "sheet music" for the loudest, most violent phase transitions the universe could have had.
• It tells us that even if the physics is incredibly complex and non-linear (chaotic), the final signal has a predictable shape and strength.

In summary: The universe had a violent childhood where bubbles of new reality collided. Whether they collided gently or violently, they created a specific "song" of gravitational waves. This paper helps us tune our instruments to listen for that song, proving that even in the chaos of the early universe, there is a surprising amount of order and predictability.

Technical Summary

1. Problem Statement

The detection of a stochastic gravitational wave background (SGWB) is a primary goal of future observatories like LISA and Pulsar Timing Arrays. While weak to intermediate strength first-order phase transitions (FOPTs) are well-modeled by the "sound shell" model, strong FOPTs (where fluid velocities become relativistic) present significant challenges.

• Non-linearity: In strong transitions, fluid velocities are high, leading to shock formation and turbulence. Linear approximations and standard sound shell models often fail to capture the complex interplay between compressional (acoustic) and vortical (turbulent) modes.
• Simulation Limitations: Previous 3D simulations of strong transitions were limited in spatial volume and duration, preventing the observation of long-term energy decay and the full development of turbulent cascades.
• Uncertainty: There is a lack of understanding regarding how kinetic energy dissipates in strong transitions (via acoustic turbulence vs. vortical turbulence) and how this decay affects the final gravitational wave (GW) power spectrum and efficiency.

2. Methodology

The authors performed large-scale, long-running 3D numerical simulations of a system containing a scalar order parameter (driving the phase transition) coupled to a relativistic fluid.

• Physical Model:
  • Used the "bag model" equation of state with a scalar field potential to drive the transition.
  • Coupled the scalar field to the fluid via a phenomenological friction term ($\eta$) to drive bubble wall expansion.
  • Assumed no seed magnetic fields; GWs are sourced solely by fluid anisotropic stresses.
• Simulation Setup:
  • Grid: $4096^3$ lattice (significantly larger than previous studies) to resolve scales from bubble wall thickness to the sound horizon.
  • Cases Studied: Two extreme cases from previous work [32] were selected to contrast different hydrodynamic regimes:
    1. Detonation: Strong transition ($\alpha_n = 0.67$) with supersonic wall speed ($v_w = 0.92$).
    2. Deflagration: Strong transition ($\alpha_n = 0.50$) with subsonic wall speed ($v_w = 0.44$).
  • Duration: Simulations ran for $\sim 1.5$ shock decay times ($t_{sh}$), allowing observation of kinetic energy dissipation and integral scale growth.
• Observables:
  • Measured power spectra of 3-velocity ($v$) and enthalpy-weighted 4-velocity ($U$).
  • Calculated shear stress power spectra ($P_\Pi$).
  • Novelty: For the first time in phase transition simulations, the authors computed Unequal Time Correlators (UETCs) and decorrelation functions for velocity and shear stress to understand the temporal evolution of the source.

3. Key Contributions

• First-time UETC Analysis: The study provides the first detailed analysis of time decorrelation in FOPT simulations, revealing that decorrelation is driven by shock motion rather than simple sound wave propagation.
• Enthalpy-Weighted 4-Velocity: The authors demonstrate that using the enthalpy-weighted 4-velocity ($U^i$) rather than the 3-velocity ($v^i$) provides a better description of the shear stress source, particularly for relativistic flows where enthalpy and velocity are strongly correlated.
• Decoupling of Modes: The work clarifies the relative contributions of compressional vs. vortical modes to GW production, challenging the assumption that vortical turbulence is a dominant source in all strong transitions.
• Efficiency Factor: A robust derivation of the GW production efficiency for strong transitions, showing it is approximately constant ($\sim 0.017$) regardless of the specific transition strength or wall speed, provided the flow decays within a Hubble time.

4. Key Results

A. Hydrodynamic Evolution

• Detonation: The flow is dominated by compressional modes and strong shocks. The enthalpy-weighted 4-velocity is $\sim 40\%$ higher than the 3-velocity, indicating strong correlations typical of shocks. Vortical modes are negligible ($U_\perp \approx 0.05$). The integral scale grows over time as energy cascades to smaller scales.
• Deflagration: The flow generates significant vorticity due to the collision of sound shells. While compressional modes dominate initially, vortical modes grow to become comparable in magnitude ($U_\perp \approx 0.08$ vs $U_\parallel \approx 0.11$) by the end of the simulation. The integral scale remains roughly constant, unlike the detonation.
• Decay: Kinetic energy decays via power laws. The detonation follows the theoretical prediction for acoustic turbulence ($\zeta \approx 1.48$), while the deflagration decays slower ($\zeta \approx 0.86$), likely due to the complex interplay of vorticity and reheating effects.

B. Decorrelation Functions (UETCs)

• Compressional Modes: Show oscillatory behavior modulated by a Gaussian envelope. The oscillation frequency is faster than the sound speed ($M_\parallel > 1$), indicating that shocks drive the decorrelation.
• Vortical Modes:
  • In the deflagration, vortical modes decorrelate via the standard Kraichnan "random sweeping" model (Gaussian decay), where they are swept by the background compressional flow.
  • In the detonation, the subdominant vortical modes show a novel "oscillatory Gaussian" behavior, suggesting they are being "washed" back and forth by the dominant acoustic oscillations.

C. Gravitational Wave Production

• Source Dominance: Despite the significant vorticity in the deflagration, the gravitational wave power spectrum is dominated by purely compressional modes ($C[U_\parallel, U_\parallel]$). The mixed and pure vortical contributions are negligible.
• Spectral Shape: The GW spectra in both cases exhibit a sharp peak at $kR_* \approx 10$ and a high-frequency tail scaling as $k^{-3}$, characteristic of shock-dominated flows.
• Model Comparison:
  • The Sound Shell Model fails to predict the growth rate for the detonation (due to non-linear shock interactions).
  • For the deflagration, the Sound Shell Model works well only when a phenomenological suppression factor (accounting for wall slowing/reheating) is applied.
• Efficiency: The total GW power approaches a constant asymptote. The dimensionless efficiency factor is found to be:
  $$\tilde{\Omega}_{gw}^\infty \simeq 0.017$$
  This holds for both the detonation and deflagration, despite their vastly different kinetic energy densities.

D. Predicted Observables

For a phase transition occurring at $T_n \sim 100$ GeV with $g_* \sim 100$:

• Detonation ($v_w=0.92, \alpha_n=0.67$): $\Omega_{gw,0}/(H_n R_*)^2 = (4.8 \pm 1.1) \times 10^{-8}$.
• Deflagration ($v_w=0.44, \alpha_n=0.50$): $\Omega_{gw,0}/(H_n R_*)^2 = (1.3 \pm 0.2) \times 10^{-8}$.
• Peak Frequency: $f_p \simeq 26 (H_n R_*)^{-1} (T_n/100 \text{ GeV}) \, \mu\text{Hz}$.

5. Significance and Implications

• Refining Predictions: The paper provides critical corrections to theoretical models (like the Sound Shell Model) for strong transitions, highlighting the necessity of including non-linear shock dynamics and time-decorrelation effects.
• Vorticity vs. GWs: It resolves a long-standing ambiguity by showing that even when vorticity is dynamically significant (as in deflagrations), it does not significantly contribute to the GW signal. This simplifies future modeling efforts.
• Baryogenesis Connection: The authors note that the reheating effect in deflagrations (which slows the wall and suppresses GWs) might allow for a "sweet spot" where the transition is strong enough to generate observable GWs and slow enough to facilitate Electroweak Baryogenesis, which typically requires subsonic walls.
• Future Observations: The derived efficiency factors and spectral shapes provide more accurate templates for LISA and PTA data analysis, particularly for strong transitions in Beyond-Standard-Model physics.

In summary, this work establishes that shock-driven acoustic turbulence is the primary engine for GW production in strong FOPTs, regardless of the presence of significant vorticity, and provides a robust, simulation-backed efficiency factor for estimating the observable signal.