⚛️ phenomenology

Nonequilibrium framework for first-order phase transitions

This paper employs the Schwinger-Keldysh closed-time-path formalism to derive effective actions and coupled Langevin equations for $\phi^6$ and two-scalar-field models, establishing a non-equilibrium theoretical framework essential for simulating bubble nucleation and predicting signals for electroweak baryogenesis and gravitational waves.

Original authors: Jie Liu, Ligong Bian

Published 2026-10-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Jie Liu, Ligong Bian

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the earliest moments of the universe, just after the Big Bang, the cosmos was a seething, super-hot soup of particles and energy. As this primordial fireball expanded and cooled, it underwent a series of dramatic shifts, much like water turning into ice, but occurring on a scale that shaped the very fabric of reality. These shifts are known as phase transitions. One specific type, called a first-order phase transition, is particularly violent and important. Instead of changing smoothly, the universe would have bubbled with pockets of a new state, expanding and colliding like boiling water. Scientists believe these violent events could have generated the gravitational waves that ripple through space today and might even explain why the universe is made of matter rather than antimatter. However, understanding exactly how these transitions happened is difficult because the standard tools used to study them assume the universe was in a calm, balanced state. In reality, these transitions are chaotic, fast-moving processes where the fields driving the change are constantly jostled by the heat of the environment, creating friction and random jitters that standard theories cannot capture.

A team of researchers at Chongqing University has developed a new mathematical framework to describe this chaotic, real-time behavior. They focused on the "first-order" type of transition, which is central to theories about how the early universe evolved. The researchers realized that the traditional methods, which calculate the energy of a system as if it were sitting still in a thermal bath, fail to describe the actual motion of the fields as they change. To fix this, they turned to a sophisticated method called the closed-time-path formalism. This approach allows physicists to track how a system evolves forward in time while accounting for the fact that it is constantly interacting with a hot environment. By using this method, the team could derive equations that describe the background fields not just as static shapes, but as moving objects subject to two new forces: a drag force that slows them down, and a random, jittery force that pushes them around, both caused by the heat of the early universe.

The researchers first tested their new framework on a simplified model involving a single type of field, which they designed to mimic the conditions needed for a first-order transition. They calculated how this field would behave when it was constantly bombarded by thermal fluctuations. Their calculations showed that the field does not move smoothly; instead, it experiences a friction-like resistance and is buffeted by random noise. They successfully translated these complex quantum effects into a set of equations known as Langevin equations. These equations are a standard tool in physics for describing how particles move through a fluid, but here they were derived from first principles to describe the fundamental fields of the universe itself. The team found that the friction and the random noise are not separate phenomena; they are two sides of the same coin, linked by a fundamental rule that connects how much a system resists motion to how much it jitters.

To make their findings more relevant to the real universe, the researchers then extended their work to a more complex scenario involving two interacting fields. This setup is designed to mimic theories where the known Higgs field interacts with a new, hidden field, a possibility that could explain dark matter or the nature of the electroweak phase transition. In this two-field system, the behavior becomes even richer. The researchers discovered that the two fields do not just experience their own friction and noise; they also influence each other. The motion of one field creates a drag on the other, and the random jitters affecting one field are correlated with the jitters affecting the other. They derived a new set of coupled equations that describe this intertwined dance of friction and fluctuation. These equations show that the two fields are locked together in their response to the thermal environment, meaning that the evolution of one cannot be understood without considering the other.

The ultimate goal of this work is to provide a reliable foundation for simulating the birth of bubbles in the early universe. In a first-order phase transition, the new phase of the universe appears as bubbles that nucleate and grow. Previous attempts to calculate how often these bubbles form relied on static snapshots of energy, which miss the dynamic, real-time evolution. The new equations developed by the researchers offer a way to run computer simulations that follow the actual time-evolution of these fields, including the effects of heat-induced friction and noise. This allows scientists to predict the rate at which bubbles form and how they expand with much greater accuracy. While the current work focuses on simplified models and relies on certain approximations to make the math tractable, it establishes a clear path forward. The researchers suggest that their framework can be expanded to include other particles and forces, eventually leading to a complete, real-time description of the universe's most violent moments. This would not only refine our understanding of the early cosmos but also improve predictions for the gravitational wave signals that future space-based detectors might one day hear.

Technical Summary: Nonequilibrium Framework for First-Order Phase Transitions

Problem Statement
First-order phase transitions (PTs) are critical for understanding electroweak baryogenesis and the generation of stochastic gravitational wave backgrounds. However, traditional analyses rely on the effective potential within the imaginary-time formalism of thermal field theory. While this approach successfully calculates equilibrium properties such as critical temperatures and nucleation rates via Euclidean bounce solutions, it fails to capture the real-time, non-equilibrium dynamics of the field evolution. In a real-time setting, background fields interact continuously with a thermal environment, leading to dissipation and stochastic thermal fluctuations. Conventional effective potential methods do not inherently describe these nonequilibrium effects, necessitating a framework that derives equations of motion directly from finite-temperature quantum field theory.

Methodology
The authors employ the Schwinger-Keldysh closed-time-path (CTP) formalism, a real-time finite-temperature field theory framework, to derive the effective action and corresponding equations of motion. The study proceeds in two stages:

Single-Field Model ( $\phi^6$ ): The authors analyze a scalar field model with a potential containing $\phi^2$ , $\phi^4$ , and $\phi^6$ terms. This model is chosen for its ability to generate a barrier for first-order PTs and its relevance as a low-energy effective theory.
Two-Field Model: The formalism is extended to a system of two interacting real singlet scalar fields ( $H$ and $S$ ), mimicking scalar extensions of the Standard Model (e.g., Higgs-portal models).

The derivation involves a perturbative loop expansion of the CTP effective action up to two loops and order $\lambda^2$ (and $c_6$ ). The authors utilize the physical-response basis transformation ( $\phi_+ = \phi_c + \phi_\Delta/2$ , $\phi_- = \phi_c - \phi_\Delta/2$ ) to separate the physical field ( $\phi_c$ ) from the response field ( $\phi_\Delta$ ). To handle the imaginary parts of the effective action arising from loop integrals, the Hubbard-Stratonovich identity is applied to introduce auxiliary stochastic fields.

Key approximations are employed to transition from non-local integro-differential equations to local Langevin equations:

Quasiadiabatic Approximation: Assuming the background field varies slowly in time.
Zero External Momentum: Considering nearly spatially uniform field configurations.
Dressed Propagator Scheme: Incorporating a finite thermal width ( $\Gamma$ ) into the propagator to ensure non-vanishing dissipation coefficients.
Narrow-Width and Markovian Limits: Assuming the correlation time of thermal fluctuations is much shorter than the background evolution time, allowing non-local noise kernels to be approximated as local white noise.

Key Contributions and Results

Derivation of Langevin Equations: The paper derives the effective equations of motion for the background fields, which take the form of coupled Langevin equations. These equations include:
- Thermal Mass and Coupling Corrections: Temperature-dependent corrections to the mass ( $m_T$ ) and coupling constants ( $\lambda_T$ ).
- Dissipative Terms: Non-local kernels are reduced to local friction terms (e.g., $\eta \phi^2 \dot{\phi}$ ) derived consistently from the imaginary part of the self-energy.
- Stochastic Noise: The introduction of Gaussian noise fields ( $\xi$ ) representing thermal fluctuations. The noise correlations are shown to satisfy the fluctuation-dissipation relation in the Markovian limit.
Single-Field Results: For the $\phi^6$ model, the authors obtain a Langevin equation with a multiplicative noise term ( $\phi_c \xi_1$ ) and an additive noise term ( $\xi_2$ ). They demonstrate that in weakly interacting models, the dominant dissipation arises from the multiplicative noise channel. The dissipation coefficient $\eta_1$ is explicitly calculated in the high-temperature limit.
Two-Field Results: For the coupled two-field system, the authors derive a set of coupled Langevin equations. A significant finding is the emergence of cross-dissipation terms (e.g., $\eta_{\phi\chi} \dot{\chi}$ ) and cross-correlated noise between the two fields. These terms arise from the mixing of the fluctuation fields induced by the interaction potential. The authors provide explicit expressions for the dissipation matrix $\eta_{ij}$ and the noise correlation matrix, showing that the off-diagonal elements vanish only when the interaction coupling ( $\lambda_{hs}$ ) is zero.
Technical Formalism: The paper provides detailed appendices outlining the perturbative loop expansion, the calculation of CTP propagators at higher-loop orders, and the evaluation of master integrals (thermal kernels) required for the dissipation and noise coefficients.

Significance and Claims
The authors claim that this work establishes a rigorous theoretical framework connecting microscopic finite-temperature field theory to the real-time stochastic dynamics of cosmological first-order phase transitions. By deriving the Langevin equations from first principles without introducing phenomenological friction or noise by hand, the framework provides a consistent basis for:

Bubble Nucleation Simulations: The derived equations allow for real-time stochastic lattice simulations to track the evolution of thermal field configurations and calculate bubble nucleation rates dynamically. This complements existing Euclidean bounce methods.
Predictive Power: The framework lays the foundation for more reliable predictions of baryogenesis and gravitational wave signals by accounting for non-equilibrium effects that equilibrium potentials miss.

The paper modestly notes that the current derivation relies on specific approximations (quasiadiabatic, Markovian, narrow-width) and does not yet self-consistently describe full bubble wall dynamics or nonequilibrium particle transport. The authors identify the extension of this formalism to include gauge and fermionic degrees of freedom, as well as the application of the two-particle-irreducible (2PI) effective action and Kadanoff-Baym equations for a more complete non-equilibrium treatment, as necessary future steps.

More like this