From friction scaling to an efficient method for estimating bubble wall velocity

This paper presents a unified framework that links microscopic Boltzmann equation treatments to phenomenological hydrodynamic models by deriving a power-law scaling for the friction parameter, thereby offering an efficient method to accurately estimate bubble wall velocities in first-order cosmological phase transitions.

The Gist

The Big Picture: Bubbles in the Early Universe

Imagine the very early Universe as a pot of boiling water. As it cools down, it doesn't just get colder uniformly; instead, it starts forming bubbles of a new state (like ice forming in water, but for the fundamental forces of nature). These are called cosmological phase transitions.

Inside these bubbles, the laws of physics are slightly different than outside. As these bubbles expand, their walls (the boundary between the old and new physics) push against the "soup" of particles filling the Universe.

The Problem: How fast do these bubble walls move?

• If they move too slowly, they might not create the conditions needed for life (like the imbalance of matter and antimatter).
• If they move too fast, they might create a different kind of signal.
• Scientists want to know the exact speed to predict what gravitational wave detectors (like LISA) might hear in the future.

The Old Way: The "Guess-and-Check" Method

For a long time, scientists had two ways to calculate this speed, and both had issues:

1. The "Perfect World" Method (Hydrodynamics): This assumes the bubble wall moves through a frictionless, perfect fluid. It's like calculating how fast a car goes on a perfectly smooth, empty highway. It's easy to calculate, but it ignores the fact that the road is actually full of potholes and traffic. It usually predicts the wall moves too fast.
2. The "Microscopic" Method (Boltzmann Equations): This is the "realistic" method. It tracks every single particle bumping into the bubble wall, like counting every single pebble the car hits. It's incredibly accurate, but it's also painfully slow to calculate. It's like trying to simulate every single raindrop hitting a windshield to figure out how fast the car is going. It takes supercomputers days to get an answer.

The New Discovery: Finding the "Secret Recipe"

The authors of this paper wanted to bridge the gap. They wanted the speed of the "Perfect World" method but the accuracy of the "Microscopic" method, without the wait time.

They did this by looking at friction.

Think of the bubble wall as a snowplow clearing a path through a blizzard.

• The snow is the plasma (hot particles).
• The plow is the bubble wall.
• Friction is the resistance the snow gives the plow.

The team realized that the amount of friction depends on a specific ratio: How "strong" the phase transition is compared to the temperature.

They discovered a simple power-law rule (a mathematical shortcut). They found that if you know the strength of the transition (let's call it the "snow density"), you can predict the friction without simulating every single snowflake.

The Formula:
They found that the friction parameter ($\tilde{\eta}$) scales with the fourth power of the transition strength ($v_n/T_n$).

Analogy: Imagine that if you double the "snow density," the friction doesn't just double; it increases by 16 times ($2^4$). Once you know this rule, you don't need to count the snowflakes; you just plug the density into the formula.

How They Proved It

1. The Benchmark: They used a specific model of physics (the xSM model) as their test kitchen.
2. The Comparison: They ran the slow, heavy "Microscopic" simulation (using a tool called WallGo) to get the true friction values.
3. The Shortcut: They then looked at the data and found that the friction followed that simple "fourth power" rule perfectly.
4. The Result: They created a new method where you:
   • Calculate the friction using their simple formula.
   • Plug that friction into the "Perfect World" fluid equations.
   • Get a wall velocity that is almost identical to the super-computer simulation, but calculated in seconds instead of days.

Why This Matters

This is a game-changer for cosmologists.

• Before: To study a new theory of the early Universe, scientists had to run expensive, slow simulations to see if it produced the right gravitational waves. Many theories were ignored because the math was too hard.
• Now: Scientists can use this "friction scaling" rule to quickly test thousands of different theories. They can instantly see which ones produce the right bubble wall speeds to explain the Universe we see today.

Summary in One Sentence

The authors discovered a simple mathematical "rule of thumb" for how much friction a bubble wall feels in the early Universe, allowing scientists to predict the speed of these cosmic bubbles with high accuracy in a fraction of the time it used to take.

Technical Summary

1. Problem Statement

Cosmological first-order phase transitions (FOPTs) are critical for understanding early Universe phenomena, such as baryogenesis and the generation of stochastic gravitational waves (GWs). A central quantity governing these processes is the terminal velocity of the expanding bubble walls ($v_w$).

• The Challenge: Determining $v_w$ from first principles is non-trivial.
  • Local Thermal Equilibrium (LTE): Assuming LTE provides an upper bound on $v_w$ but ignores non-equilibrium friction, which is crucial for accurate predictions.
  • Phenomenological Models: Hydrodynamic simulations often use a friction parameter ($\tilde{\eta}$) as a free parameter. Without a microscopic derivation, these models lack predictive power for specific particle physics models.
  • Microscopic Calculations: Full kinetic theory treatments (solving Boltzmann equations) are computationally expensive and complex, often requiring specialized packages like WallGo.

The authors aim to bridge the gap between computationally efficient phenomenological models and rigorous microscopic kinetic theory to provide a fast, reliable method for estimating $v_w$.

2. Methodology

The study utilizes the xSM (SM extended by a real scalar singlet) as a benchmark model. The methodology proceeds through four main stages:

A. Hydrodynamic Framework

The authors model the bubble expansion using relativistic hydrodynamics. The plasma is treated as a fluid coupled to the scalar field via an out-of-equilibrium friction term ($F$). The system is solved using:

• Conservation of energy and momentum across the wall.
• An entropy balance condition (Eq. 3.9) that relates the temperature jump across the wall to the entropy production ($\Delta S$) caused by friction.
• The friction parameter $\tilde{\eta}$ is defined as a rescaled measure of this entropy production.

B. Comparison of Friction Models

The authors compare three approaches to modeling the friction source term ($f_s$):

1. Phenomenological Ansatz: $F \propto u^\mu \partial_\mu \phi$ (constant coefficient $\eta$).
2. Chapman-Enskog (CE) Expansion: $F \propto \phi^2 (u^\mu \partial_\mu \phi)$ (coefficient $\kappa$).
3. Microscopic (WallGo): Full numerical solution of Boltzmann equations for the top quark (the dominant contributor to friction).

They demonstrate that both phenomenological ansätze can accurately reproduce the shape of the microscopic entropy source profile when normalized, validating the use of simplified forms.

C. Microscopic Derivation and Analytic Approximation

To derive the scaling of friction from first principles:

• They solve the Boltzmann equation for the top quark distribution function, truncating the momentum expansion at the zeroth order (effective chemical potential, $\mu_t$).
• They derive an analytic approximation for $\mu_t$ assuming a hyperbolic tangent profile for the Higgs field.
• By analyzing the peak of the chemical potential and the resulting friction force, they extract the dependence of the friction parameter on the phase transition strength.

D. Iterative Algorithm

The authors propose a hybrid algorithm to compute $v_w$:

1. Start with an LTE guess for $v_w$.
2. Solve the simplified Boltzmann equation for the chemical potential $\mu_t$ to find the friction profile.
3. Fit the friction profile to the phenomenological ansatz to extract $\tilde{\eta}$.
4. Use $\tilde{\eta}$ in the hydrodynamic matching conditions to update $v_w$.
5. Iterate until convergence.

3. Key Contributions and Results

A. Discovery of a Power-Law Scaling

The most significant result is the derivation of a simple power-law relationship between the rescaled friction parameter ($\tilde{\eta}$) and the ratio of the Higgs vacuum expectation value to the nucleation temperature ($v_n/T_n$):
$$\tilde{\eta} = \zeta \left( \frac{v_n}{T_n} \right)^4$$

• Coefficient: Fitting to full WallGo data yields $\zeta_{\text{WG}} \approx 2.1 \times 10^{-4}$.
• Theoretical Derivation: The analytic approximation of the chemical potential independently predicts a similar scaling with $\zeta_{\mu} \approx 2.5 \times 10^{-4}$.
• Physical Origin: The scaling arises because the enthalpy density scales as $T^4$ and the wall width scales as $1/T$, while the chemical potential (and thus friction) scales with $(v_n/T)^2$.

B. Efficient Estimation Method

The authors demonstrate that incorporating this power-law scaling into the phenomenological framework allows for the accurate estimation of $v_w$ without solving the full set of Boltzmann equations.

• Accuracy: For slow walls ($v_w < 0.3$), the chemical potential method matches WallGo results closely.
• Extrapolation: By using the scaling law (Eq. 4.8) with the fitted coefficient $\zeta$, the method reproduces WallGo results for fast walls (deflagrations, hybrids, and detonations) with high accuracy, significantly outperforming the LTE approximation and the bare chemical potential method (which underestimates friction at high speeds).

C. Validation

• Figure 3 & 4: Comparisons show that the "scaling method" (using $\tilde{\eta} \propto (v_n/T_n)^4$) aligns with full numerical WallGo solutions across a wide range of parameters, including cases where gauge bosons are included.
• Limitations: The pure chemical potential truncation fails for very fast walls ($v_w \gtrsim 0.3$) because higher-order momentum terms become significant. However, the scaling law derived from it remains valid.

4. Significance

This work provides a unified and efficient framework for studying cosmological phase transitions:

1. Computational Efficiency: It reduces the computational cost of determining bubble wall velocities from solving complex integro-differential Boltzmann equations to solving simple ordinary differential equations or using an analytic formula.
2. Predictive Power: It removes the "free parameter" problem in phenomenological hydrodynamic simulations by providing a theoretically grounded, model-dependent value for the friction parameter.
3. Gravitational Wave Phenomenology: Since the spectrum of gravitational waves from FOPTs is highly sensitive to $v_w$, this method enables rapid and reliable scans of particle physics models (like the xSM) to predict observable GW signatures for upcoming detectors (e.g., LISA).
4. Generalizability: While tested on the xSM, the methodology for extracting the friction scaling from microscopic dynamics is applicable to other Standard Model extensions.

In summary, the paper establishes that the complex microscopic friction dynamics can be effectively captured by a simple power-law scaling, enabling fast, accurate, and reliable predictions for bubble wall velocities in early Universe cosmology.