Dynamical evolution of the pressure on the bubble wall

Technical Summary: Dynamical Evolution of the Pressure on the Bubble Wall

Problem Statement
First-order phase transitions (FOPTs) in the early Universe are critical mechanisms for generating gravitational waves, baryogenesis, and dark matter. A central uncertainty in modeling these transitions is determining whether expanding bubble walls reach a terminal constant velocity (deflagration, detonation, or hybrid regimes) or accelerate indefinitely to ultra-relativistic speeds (runaway regime).

Traditional analyses rely on the assumption of Local Thermal Equilibrium (LTE) and a "steady-state" (SS) fluid profile, where the plasma surrounding the bubble wall instantaneously adjusts to the wall's motion. In this framework, hydrodynamic obstruction arises: as the wall accelerates, it generates a compression (heating) wave in the plasma. Near the Jouguet velocity, this heating maximizes the effective pressure opposing the wall, potentially halting acceleration. However, the steady-state approach assumes the fluid profile tracks the wall's velocity instantaneously. This paper argues that this assumption fails during the early acceleration phase because the timescale for the heating wave to form may exceed the timescale for the wall to accelerate to the Jouguet velocity. Consequently, steady-state predictions may incorrectly rule out runaway solutions or detonations that are dynamically accessible.

Methodology
The authors employ a multi-faceted approach combining analytical approximations with numerical simulations to study the non-equilibrium dynamics of the bubble wall and the surrounding plasma.

1. Analytical Framework:

   • The authors derive the equation of motion for the bubble wall radius $R(t)$, incorporating the pressure $P(t)$ derived from the effective potential and plasma backreaction.
   • They perform a Taylor expansion of the net pressure $P(t)$ around $t=0$ to analyze the early-time evolution. This allows them to estimate the competition between the driving vacuum pressure and the developing hydrodynamic backreaction.
   • They derive a critical timescale comparison: the time required for the plasma to build up sufficient friction (shock formation time) versus the time required for the wall to accelerate to the Jouguet velocity.

2. Numerical Simulations:

   • Fixed-Scalar Simulations: Initially, the authors simulate the hydrodynamic equations (energy-momentum conservation) with a prescribed, accelerating scalar field profile. This isolates the fluid response time without the complexity of solving the coupled scalar field equation of motion (EoM).
   • Dynamical Scalar-Fluid Simulations: The core of the study involves a fully coupled 1+1 dimensional simulation (assuming spherical symmetry) where both the scalar field EoM and the fluid conservation equations are solved simultaneously.
   • Algorithm: The simulations utilize a pseudospectral method with Chebyshev polynomials to discretize spatial derivatives, ensuring high precision and exponential convergence. The initial conditions are set based on the nucleation bounce solution, linearized to evolve forward in time.

3. Parameter Scan:

   • The authors perform a broad parameter scan over the phase transition strength $\alpha_p$, the enthalpy ratio $\Psi$ (related to the speed of sound and degrees of freedom), and the vacuum expectation value ratio $\phi_0/T_n$.
   • They compare the outcomes of these dynamical simulations (runaway vs. terminal velocity) against predictions made by both the traditional steady-state LTE criterion and their newly derived dynamical criterion.

Key Contributions and Results

1. Timescale Mismatch: The simulations confirm that the formation time of the heating wave (shock) is often comparable to or longer than the wall's acceleration timescale, particularly for hybrid solutions near the Jouguet velocity. For these regimes, the fluid profile does not reach the steady-state attractor before the wall overtakes the shock front.

2. Revised Criterion for Maximal Driving Pressure:
   The paper derives a new criterion, $\alpha_p^{\max}$, which determines whether hydrodynamic obstruction develops fast enough to prevent a runaway regime. This criterion is the minimum of two terms:
   $$\alpha_p^{\max} \approx \min \left[ \alpha_{p,\max}^{\text{SS}}, \quad 60.64 \left( \frac{1-\Psi}{\phi_0/T_n} \right)^3 \left( \frac{\nu}{4} \right) \left( \frac{S_3/T_n}{140} \right) \left( \frac{g_*^+}{106.75} \right)^{1/2} \left( \frac{\gamma_w^{\max}}{\xi_w^{\max}} \frac{1}{1.381} \right)^3 \right]$$

   • The first term, $\alpha_{p,\max}^{\text{SS}}$, is the traditional steady-state limit.
   • The second term is the new dynamical limit, derived from the early-time expansion and validated by numerical fits. It depends on the nucleation action $S_3$, the wall thickness, and the coupling strength.

3. Validation and Accuracy:

   • The new criterion was tested against a scan of $\sim 10^4$ parameter points. It correctly classifies the expansion regime (runaway vs. terminal) with approximately 98% accuracy.
   • The results show that hydrodynamic obstruction is less restrictive than steady-state LTE predictions suggest. Specifically, regions of parameter space previously thought to be blocked by hydrodynamic obstruction (e.g., certain values of $\Psi < 0.77$) can actually support runaway walls or detonations because the wall accelerates faster than the plasma can react.

4. Regime Separation:

   • For weak plasma-wall coupling (small $1-\Psi$), the dynamical timescale dominates, allowing walls to bypass the steady-state obstruction.
   • For strong coupling, the steady-state limit remains the governing constraint.

Significance and Claims
The authors claim that their work challenges the conventional assumption that bubble wall dynamics can be fully described by steady-state hydrodynamics. By demonstrating that the "steady-state tracking" hypothesis is generically violated during the acceleration phase, they provide a more accurate tool for predicting the fate of bubble walls.

• Implications for Phenomenology: The revised criterion relaxes constraints on phase transition parameters. Models that were previously deemed incompatible with runaway walls (and thus specific gravitational wave signatures or baryogenesis mechanisms) may now be viable.
• Methodological Utility: The paper proposes a two-step computational strategy: first, use the new dynamical criterion (Eq. 3) to determine if a runaway solution is possible; second, if a terminal velocity is predicted, use standard steady-state methods to calculate the specific velocity value. The authors verify that for terminal solutions, the dynamical simulations converge to the steady-state velocity predictions, validating the continued use of steady-state tools for velocity estimation once the regime is established.

Limitations Acknowledged
The authors note several limitations in their current analysis:

• Entropy Production: The study neglects friction from entropy production (dissipative friction), which acts independently of shock formation and could increase the maximal driving pressure.
• Nucleation Dynamics: The initial acceleration is modeled using classical equations of motion, whereas a Langevin-like approach might better describe the very early stochastic nucleation phase.
• Thermalization: The assumption of instantaneous thermalization ahead of the bubble front is justified for strongly coupled systems (like the QGP) but may not hold for weakly interacting particles in Beyond Standard Model (BSM) scenarios.

In summary, the paper provides a quantitative, model-independent criterion that corrects the overestimation of hydrodynamic obstruction in early Universe phase transitions, thereby expanding the viable parameter space for ultra-relativistic bubble wall scenarios.