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

Abstract: We study bubble-wall dynamics in cosmological first-order phase transitions in a two-scalar-field model, where the wall is formed by $φ$ and an additional real scalar $s$ couples through a portal interaction. We evolve the coupled classical field equations on the lattice and demonstrate that for an initial Bose--Einstein distribution of $s$ fluctuations at the nucleation temperature $T_n$, the resulting patchy background intermittently modulates the local driving pressure on the wall. The wall therefore undergoes alternating episodes of acceleration and deceleration and approaches a quasi-stationary propagation regime with a smaller time-averaged speed than in the decoupled limit. We further identify three familiar propagation profiles -- deflagration, detonation, and hybrid -- distinguished by where the dynamical $s$-sector energy density is concentrated relative to the wall. These effects can impact gravitational wave and baryogenesis predictions.

• The paper demonstrates through large-scale classical-statistical lattice simulations that a portal-coupled spectator scalar can suppress runaway bubble-wall acceleration, producing quasi-stationary motion near v ≈ 0.8 instead of v ≈ 0.976.
• The simulations show that fluctuation-induced friction transfers energy into spectator-field excitations and creates distinctive dynamics, including shrink–reexpand transients and enhanced post-wall oscillations rather than conventional fluid-like damping.
• A parameter scan identifies deflagration-, detonation-, and hybrid-like field profiles, but expanding solutions remain above roughly v ≈ 0.6, leaving the mechanism’s generality and its effects on gravitational waves and baryogenesis open for further study.

Overview and motivation

Bubble-wall velocity is a central input in the phenomenology of cosmological first-order phase transitions (FOPTs): it controls the efficiency and duration of the acoustic gravitational-wave source, the partition of released vacuum energy between field gradients and bulk motion, and—at the electroweak scale—the viability of charge transport for baryogenesis. Existing treatments of wall slowdown typically rely on semiclassical transport calculations, hydrodynamic compression heating, or the scalar–fluid framework with a phenomenological friction coefficient. The paper by Wei and Guo (2602.04586) proposes and demonstrates a complementary mechanism that requires no fluid degree of freedom: fluctuation-induced friction arising from a thermally populated spectator scalar $s$ coupled to the transition field $\phi$ through a portal interaction $\lambda_{\phi s}\phi^2 s^2$.

The core idea is classical-statistical: at the nucleation temperature $T_n$, the $s$ field is initialized with Bose–Einstein mode occupancies, so its real-time evolution produces a spatially patchy background. Because the portal term is even in $s$, every nonzero local value of $s$ modifies the effective potential along $\phi$; in patches where $s$ reaches large excursions ($s \simeq s_m$), the would-be broken-phase minimum can become metastable, strongly suppressing the local driving pressure $\phi$0. The wall therefore samples a stochastic sequence of driving pressures, producing alternating acceleration and deceleration rather than monotonic runaway.

Lattice setup and main dynamical results

The authors evolve the coupled classical equations of motion on a $\phi$1 lattice with $\phi$2 ($\phi$3), using leapfrog time integration and second-order central differences. A critical-bubble profile of $\phi$4 is embedded at the box center, with wall thickness taken from the bounce solution and initial radius slightly above criticality. Fluctuations $\phi$5 and $\phi$6 are sampled as Gaussian random variables with variances fixed by the BE occupancy $\phi$7.

Two robust numerical findings emerge:

• Suppression of runaway acceleration. In the decoupled limit ($\phi$8), the wall accelerates rapidly and saturates near $\phi$9 by $\lambda_{\phi s}\phi^2 s^2$0. With $\lambda_{\phi s}\phi^2 s^2$1, the wall instead fluctuates around $\lambda_{\phi s}\phi^2 s^2$2 with no secular drift—a quasi-stationary terminal regime maintained by continuous energy transfer into $\lambda_{\phi s}\phi^2 s^2$3 excitations.
• A shrink–reexpand transient. At early times the coupled bubble can transiently contract before re-expanding, indicating that the local driving free-energy difference is temporarily overwhelmed by the $\lambda_{\phi s}\phi^2 s^2$4-induced backreaction. This behavior has no analogue in uncoupled FOPT dynamics.

Notably, the interior $\lambda_{\phi s}\phi^2 s^2$5 profile develops enhanced post-wall oscillations in the coupled case, in contrast to scalar–fluid models where viscous damping rapidly relaxes the field behind the wall. This identifies the slowdown mechanism as an inhomogeneous modulation of the effective potential rather than smooth dissipation—an important qualitative distinction from the standard friction ansatz.

Propagation-profile classification in 1+1 dimensions

To resolve wall-scale energy deposition efficiently, the authors perform complementary $\lambda_{\phi s}\phi^2 s^2$6 dimensional simulations (verified to reproduce the qualitative $\lambda_{\phi s}\phi^2 s^2$7 behavior) and track the dynamical energy density of the $\lambda_{\phi s}\phi^2 s^2$8 sector, $\lambda_{\phi s}\phi^2 s^2$9, which avoids bookkeeping ambiguity in assigning interaction energy. Scanning $T_n$0 yields three regimes analogous to the hydrodynamic classification:

Regime	$T_n$1	Wall speed	Location of $T_n$2
Deflagration-like	0.65	near-luminal	ahead of wall
Detonation-like	1.1	$T_n$3	behind wall
Hybrid-like / collapse threshold	1.8	$T_n$4 minimum	both sides

Two quantitative claims deserve emphasis. First, within the parameter range supporting expansion, no fine-tuning of $T_n$5 drives the wall below $T_n$6; beyond this threshold the vacuum bubble collapses entirely. The authors attribute this lower bound to the strong modification of the vacuum structure at larger couplings—effectively imposing a floor on attainable wall speed in this model. Second, the mapping from speed to profile type differs from hydrodynamics: here $T_n$7 corresponds to a detonation-like profile, whereas in scalar–fluid treatments that speed is associated with detonation or hybrid solutions.

The authors are careful to note that despite superficial macroscopic resemblance to deflagration/detonation/hybrid templates, the microscopic interface structure differs fundamentally: no sharp shock front or discontinuity appears ahead of the wall (partly because the couplings used are small), and the profiles exhibit nonlinear fluctuations with no hydrodynamic analogue, since all energy transfer occurs through coherent field dynamics.

Limitations and open questions

Several caveats bound the scope of these results. The stochasticity is implemented solely through initial BE-distributed fluctuations; whether the quasi-stationary regime persists under genuinely thermal (interacting) evolution, or under alternative ensembles such as inflationary spectator fluctuations that classicalize after horizon exit, remains untested. The identification of deflagration-, detonation-, and hybrid-like profiles is based on the spatial support of $T_n$8 alone, without a fluid comparison at matched parameters, and the absence of shock fronts may be an artifact of the modest coupling values simulated. The claimed lower bound $T_n$9 is established only for one benchmark potential and one-dimensional scans; its generality across potentials and in $s$0 dimensions is open. Finally, the paper does not compute the resulting GW spectrum or baryogenesis efficiency, so the phenomenological impact of the reduced, fluctuating wall velocity—and of the enhanced post-wall oscillations as an additional GW source—is left as a quantitative open question.

Conclusion

This work provides a concrete lattice demonstration that thermal fluctuations of a portal-coupled spectator scalar act as an effective drag on expanding bubble walls, replacing monotonic acceleration with intermittent, quasi-stationary propagation at reduced mean speed, accompanied by distinctive transients such as shrink–reexpand episodes. The $s$1-based classification offers a field-theoretic analogue of the hydrodynamic wall taxonomy while revealing that the microscopic dynamics differ qualitatively from scalar–fluid descriptions. Since predicted GW amplitudes depend sensitively on $s$2, fluctuation-induced friction constitutes a mechanism that should be incorporated in precision predictions for LISA, Taiji, and TianQin.