Bubble velocities in local equilibrium from a pseudopotential

This paper introduces a novel method for calculating terminal bubble velocities in first-order phase transitions by identifying wall velocities that create degenerate minima in a modified "pseudopotential," thereby determining net outward pressure without solving scalar field equations of motion or relying on simplified equations of state.

The Gist

The Big Picture: Bubbles in the Early Universe

Imagine the early universe as a giant pot of boiling water. As it cools down, it doesn't just freeze smoothly; instead, it undergoes a "first-order phase transition." Think of this like water turning into ice, but instead of freezing all at once, little bubbles of ice start forming and expanding into the liquid water.

In the universe, these aren't ice bubbles, but bubbles of a new state of reality (called the "broken phase") expanding into the old state (the "symmetric phase"). The edge of these bubbles is called a bubble wall.

The scientists in this paper are trying to answer one specific question: How fast do these bubble walls move?

Why does speed matter?

• Too slow: It might help explain why the universe is made of matter instead of antimatter (a mystery called baryogenesis).
• Too fast: It might create ripples in space-time called gravitational waves that we could detect with future telescopes.

The Problem: A Tug-of-War

The speed of the bubble wall is determined by a tug-of-war between two forces:

1. The Driving Force: The new phase wants to expand because it's energetically favorable (like a ball rolling down a hill). This pushes the wall forward.
2. The Friction Force: As the wall moves, it pushes against the "plasma" (a hot soup of particles) surrounding it. This creates resistance, slowing the wall down.

For a long time, scientists thought friction only happened if the plasma was out of balance. However, this paper focuses on a scenario where the plasma is in Local Thermal Equilibrium (LTE). Even in this calm, balanced state, the wall still experiences friction because the temperature changes across the wall.

The Old Way vs. The New Way

The Old Way (The Hard Math):
To find the speed, scientists usually had to solve a very complicated set of equations describing how the particles move and how the field changes. It's like trying to predict the speed of a car by calculating the friction of every single tire tread, the air resistance on every bolt, and the engine's internal combustion simultaneously. It's accurate but incredibly difficult and computationally expensive.

The New Way (The "Pseudopotential" Shortcut):
The authors, Martin Münzenberg and Carlos Tamarit, developed a clever shortcut. They created a new mathematical tool they call a "pseudopotential."

Think of the "pseudopotential" as a custom-made landscape that changes shape depending on how fast the bubble is moving.

• Imagine a ball rolling on a hilly surface. The ball naturally wants to roll into the lowest valley (the minimum).
• In their new method, they look at this "pseudopotential" landscape.
• If the bubble is moving at the correct terminal speed, the landscape will have two valleys (one for the old phase, one for the new phase) that are exactly the same height.
• If the valleys are different heights, the bubble will either speed up (if the new phase is lower) or slow down (if the old phase is lower).

The Analogy:
Instead of calculating every tiny force acting on the car (the old way), the new method is like looking at a map of the road. If the start and end points of a hill are at the same elevation, the car will cruise at a steady speed without needing to accelerate or brake. If one side is lower, the car will speed up or slow down until it finds the right speed where the "hill" levels out.

What They Did and Found

1. The Test: They tested this new "pseudopotential" method on a specific model of the universe (a version of the Standard Model with extra particles).
2. The Result: They found that their shortcut gave the exact same answers as the difficult, full-math method. This proves the shortcut is accurate.
3. The Discovery:
   • They confirmed that "deflagrations" (bubbles moving slower than the speed of sound in the plasma) are stable.
   • They found that "detonations" (bubbles moving faster than the speed of sound) are unstable. It's like trying to push a car up a hill that gets steeper the faster you go; it won't stay at a constant speed.
   • They confirmed a "dip" in the pressure found in previous studies, which explains why certain types of bubbles are stable while others aren't.

Why This Matters

This paper doesn't just give a new number; it gives a new tool.

• Simplicity: You don't need to solve the hardest equations to get the answer.
• No Guessing: Other methods often require scientists to guess the shape of the bubble wall. This method doesn't need that guess; it derives the answer directly from the physics.
• Insight: It allows scientists to easily see why a bubble is stable or unstable by just looking at the shape of their "pseudopotential" landscape, rather than getting lost in complex calculations.

In short, the authors found a way to predict how fast cosmic bubbles move by looking at a "balance of heights" in a mathematical landscape, skipping the need to solve the most difficult parts of the physics equations.

Technical Summary

Technical Summary: Bubble Velocities in Local Equilibrium from a Pseudopotential

Problem Statement
The paper addresses the challenge of determining the terminal velocity of expanding bubble walls during first-order cosmological phase transitions in the early universe. While the expansion velocity is crucial for predicting phenomena such as baryogenesis and gravitational wave backgrounds, calculating it typically requires solving complex coupled systems of scalar field equations of motion and hydrodynamic equations. Existing methods often rely on simplifying assumptions, such as the "bag" equation of state (which neglects temperature-dependent mass scales), specific ansatzes for the scalar field profile (e.g., hyperbolic tangent), or simplified entropy matching constraints that may not capture the full thermodynamic behavior of the plasma. Furthermore, while it is known that hydrodynamic backreaction in local thermal equilibrium (LTE) provides a friction force that can prevent runaway bubble acceleration, accurately predicting the terminal velocity without solving the full scalar equation of motion remains a computational hurdle.

Methodology
The authors propose a new method to estimate terminal bubble velocities in LTE by introducing a "pseudopotential" function, $\hat{V}(\phi)$, derived from the balance of forces acting on the bubble wall. The core of the method relies on the following steps:

1. Force Balance Formulation: Starting from the scalar equation of motion and hydrodynamic conservation laws, the authors derive an integral condition representing the balance between the driving pressure (difference in effective potential) and the hydrodynamic backreaction force (arising from temperature gradients).
2. Definition of the Pseudopotential: Under the assumption that the dependence of the fluid temperature profile on the scalar field gradients ($\partial_z \phi$) is negligible (a condition verified to hold at the per-mille level in their examples), the authors define a generalized potential:
   $$\hat{V}(\phi) \equiv \int_0^\phi d\phi' \frac{\partial V(\phi', T(\phi'))}{\partial \phi'}$$
   Here, $T(\phi)$ is the temperature profile determined by the hydrodynamic equations (specifically the constants of motion $c_1$ and $c_2$) for a given wall velocity $v_w$.
3. Terminal Velocity Criterion: The net outward pressure acting on the bubble wall is identified as the difference in the pseudopotential between the symmetric phase ($\phi_+$) and the broken phase ($\phi_-$): $\Delta \hat{V} = \hat{V}(\phi_+) - \hat{V}(\phi_-)$.
   • For a bubble moving at a constant terminal velocity, the net outward pressure must be zero.
   • Therefore, the correct terminal velocity is the value of $v_w$ that renders the minima of the pseudopotential degenerate ($\Delta \hat{V} = 0$).
4. Implementation: The method allows for the calculation of bubble velocities by scanning the wall velocity parameter until the degeneracy condition is met, without explicitly solving the second-order differential equation for the scalar field profile. This approach applies to deflagrations, detonations, and hybrid solutions.

Key Contributions

• New Computational Framework: The paper introduces a method to determine bubble wall velocities that avoids solving the scalar equation of motion directly and does not require simplified equations of state or specific field profile ansatzes.
• Physical Interpretation of Pseudopotential: The authors demonstrate that the difference in pseudopotential values at extrema corresponds directly to the net outward pressure. This provides a physical quantity that can be analyzed for stability, unlike entropy-matching methods which primarily yield static solutions.
• Stability Analysis: By analyzing the slope of the net outward pressure as a function of wall velocity, the method offers insights into the stability of static solutions. A negative slope indicates stability (deflagrations), while a positive slope indicates instability (detonations).
• Validation: The method is tested against exact numerical solutions of the scalar equation of motion in a specific model, showing high accuracy.

Results
The method was applied to a singlet extension of the Standard Model (SM) with $N$ complex scalar singlets coupled to the Higgs sector.

• Accuracy: The bubble wall velocities predicted by the pseudopotential method (requiring $\Delta \hat{V} = 0$) were found to accurately match those obtained by solving the full scalar equation of motion numerically.
• Pressure Profiles: The study computed the net outward pressure as a function of wall velocity. It confirmed the existence of a peak in the backreaction pressure, consistent with previous findings using the tanh-ansatz method, though the location of this peak was identified for hybrid solutions near the speed of sound rather than the Jouguet velocity.
• Stability Confirmation: The analysis of the pressure slope confirmed that static deflagrations are stable (negative slope), while static detonations are unstable (positive slope). This aligns with time-dependent simulations showing that static detonations are rarely reached in practice.
• Hybrid Solutions: The method successfully identified stable hybrid solutions and confirmed the stability of deflagrations, while highlighting the instability of detonation configurations within the LTE approximation.

Significance
The paper claims that this pseudopotential method offers a robust and efficient alternative to existing techniques for calculating bubble wall velocities. Its primary significance lies in its ability to:

1. Eliminate Approximations: It removes the need for simplified equations of state (like the bag model) or assumed field profiles, relying instead on the thermodynamic quantities derived directly from the effective potential.
2. Provide Stability Insights: Unlike entropy-matching approaches which only identify static solutions, the pseudopotential framework allows for the analysis of non-stationary configurations and the assessment of solution stability through the sign of the pressure slope.
3. Handle Complex Scenarios: The method is applicable to various hydrodynamic regimes (deflagrations, detonations, hybrids) and can treat cases where no stationary solutions are realized, providing a unified framework for understanding bubble dynamics in LTE.

The authors conclude that while the method relies on the approximation that temperature gradients depend weakly on scalar field gradients, this assumption holds with high accuracy in the tested models, making the pseudopotential a reliable tool for predicting terminal bubble velocities in first-order phase transitions.