Fluctuation-Induced Friction in Bubble-Wall Dynamics of Cosmological First-Order Phase Transitions

Technical Summary: Fluctuation-Induced Friction in Bubble-Wall Dynamics of Cosmological First-Order Phase Transitions

Problem Statement
Cosmological first-order phase transitions (FOPTs) proceed via the nucleation and expansion of bubbles of the broken phase. The dynamics of the bubble wall, particularly its velocity, are critical for predicting stochastic gravitational wave (GW) backgrounds and baryogenesis scenarios. While standard models often treat the surrounding medium as a relativistic fluid or use phenomenological friction terms to limit wall acceleration, the authors investigate a complementary slowdown mechanism. They address how an additional scalar degree of freedom ($s$), coupled to the phase-transition field ($\phi$) via a portal interaction, influences wall dynamics through stochastic fluctuations rather than fluid-like viscous damping.

Methodology
The authors employ a two-scalar-field model where $\phi$ drives the phase transition and a real scalar $s$ couples via a portal interaction ($\lambda_{\phi s} \phi^2 s^2$). The study utilizes real-time lattice simulations to evolve the coupled classical field equations.

• Initial Conditions: At the nucleation temperature $T_n$, the system is initialized with a spherical critical bubble for $\phi$. Crucially, the $s$ field is initialized with a homogeneous background plus fluctuations sampled from a Bose–Einstein distribution, representing a thermal ensemble of $s$-fluctuations.
• Simulation Setup:
  • 3+1 Dimensions: Simulations are performed on a cubic grid ($N^3 = 2000^3$) with physical size $L T_c = 400$ to track the expansion of spherical bubbles and measure wall velocity evolution.
  • 1+1 Dimensions: Reduced dimensional simulations are used to perform efficient scans over the portal coupling $\lambda_{\phi s}$ and to resolve the spatial distribution of energy density relative to the wall.
• Analysis: The wall velocity is tracked using an isosurface criterion ($\phi = (\phi_{false} + \phi_{true})/2$). Energy density distributions are analyzed to classify propagation profiles based on where the dynamical energy of the $s$-sector is concentrated relative to the wall.

Key Contributions and Results

1. Fluctuation-Induced Intermittency:
   The simulations demonstrate that for an initial Bose–Einstein distribution of $s$ fluctuations, the resulting patchy background modulates the local free-energy difference ($\Delta V$) driving the wall.

   • In regions where $s$ fluctuates to high values ($s \approx s_m$), the portal term suppresses the driving pressure, potentially rendering the broken-phase minimum metastable. This causes the wall to decelerate or transiently stall.
   • In regions where $s \approx 0$, the suppression is alleviated, allowing re-acceleration.
   • Result: Instead of monotonic acceleration (runaway behavior) or smooth deceleration, the wall undergoes alternating episodes of acceleration and deceleration. It settles into a quasi-stationary propagation regime with a time-averaged speed significantly lower than in the decoupled limit, characterized by persistent velocity fluctuations.

2. Distinctive Dynamical Signatures:
   The paper contrasts this mechanism with standard scalar-fluid descriptions.

   • Shrink-Reexpand Behavior: Unlike standard uncoupled dynamics, the coupled system exhibits transient contraction followed by re-expansion at early times, indicating temporary suppression of the driving force.
   • Enhanced Oscillations: The interior of the bubble shows pronounced oscillations in the $\phi$ field, suggesting sustained energy exchange between $\phi$ and $s$, rather than the rapid relaxation and damping seen in fluid-friction models.

3. Classification of Propagation Profiles:
   Using 1+1 dimensional simulations, the authors identify three propagation regimes distinguished by the spatial concentration of the $s$-sector dynamical energy density ($\rho_s$) relative to the wall:

   • Deflagration-like: $\rho_s$ builds up predominantly ahead of the wall.
   • Detonation-like: $\rho_s$ is concentrated behind the wall.
   • Hybrid: $\rho_s$ is supported both in front of and behind the wall.
   • Velocity Bounds: The study finds a lower bound on the wall velocity ($v \sim 0.6$) within the parameter range supporting expansion. Below this threshold, the vacuum bubble collapses. This bound is attributed to significant modifications of the vacuum structure by the coupling, rather than simple friction.

4. Microscopic vs. Macroscopic Differences:
   While the macroscopic energy profiles resemble hydrodynamic deflagrations or detonations, the microscopic structure differs fundamentally. The authors note the absence of sharp shock fronts and the presence of nonlinear fluctuations and local distortions inherent to field dynamics, which have no analogue in hydrodynamic treatments.

Significance
The paper claims to provide a concrete realization of fluctuation-induced friction that can obstruct bubble-wall acceleration in cosmological FOPTs without relying on phenomenological fluid-friction terms.

• Impact on GWs: Since the predicted stochastic GW signal depends sensitively on the bubble-wall velocity (controlling efficiency and the acoustic source), this mechanism can modify the GW amplitude and potentially its spectral shape.
• Impact on Baryogenesis: The intermittent, slower propagation regime affects the viability of charge transport ahead of the wall, a key factor in electroweak baryogenesis scenarios.
• Theoretical Distinction: The work highlights a qualitatively different slowdown mechanism driven by stochastic spatial inhomogeneities in the effective potential, distinct from the viscous damping of scalar-fluid models.