🔬 physics

Metastable Transitions and $Γ$ -Convergent Eyring-Kramers Asymptotics in Landau-QCD Gradient Systems

This paper establishes a rigorous analytical framework combining $\Gamma$ -convergence, spectral perturbation theory, and variational methods to prove the stability of metastable transitions and the quantitative validity of the Eyring-Kramers formula for mean transition times in Landau-type gradient systems inspired by QCD phenomenology.

Original authors: Jingxu Wu, Jie Shi

Published 2026-01-23

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jingxu Wu, Jie Shi

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

The Big Picture: A Bumpy Landscape

Imagine the physical world (specifically, a type of matter inspired by the physics of atomic nuclei, known as QCD) as a giant, rolling landscape. This landscape is made of "energy hills" and "energy valleys."

The Valleys: These are stable places where the system likes to sit. Think of a marble resting at the bottom of a bowl.
The Hills: These are the barriers the system must climb to get from one valley to another.
The "Metastable" State: Sometimes, a marble sits in a small, shallow valley. It's stable for a while, but there is a deeper, more comfortable valley nearby. The marble is "metastable"—it's stuck, but it wants to move to the deeper valley.

The Problem: How Does the Marble Jump?

In a perfect, frictionless world, the marble would never jump out of the shallow valley because it doesn't have enough energy to climb the hill. But in the real world, things are "noisy." There are tiny, random jiggles (like thermal heat or quantum fluctuations) that shake the marble.

Occasionally, these random jiggles are strong enough to push the marble over the hill and into the new valley. This is called a metastable transition.

The paper asks: How long does it take for this jump to happen, and what exactly happens during the jump?

The Old Way vs. The New Way

For decades, physicists have used a famous formula (the Eyring–Kramers formula) to guess the answer. It says the time it takes to jump depends on two things:

The Height of the Hill: The higher the hill, the longer you wait (exponentially longer).
The Shape of the Hill: How steep is the slope at the bottom and how sharp is the peak?

However, the old formula was mostly proven for simple, finite systems (like a single ball on a spring). This paper tackles a much harder problem: Infinite-dimensional systems.

Think of the difference between a single marble and a giant, wobbly jelly. The jelly isn't just one point; every part of it can wiggle independently. Describing the "shape" of a wiggling jelly is mathematically messy. The authors wanted to prove that the old formula still works perfectly for these complex, wiggly "jellies" (which represent fields in physics), even when you change the conditions or look at them on a computer screen.

The Three Main Ingredients of the Paper

1. Mapping the Terrain (Deterministic Dynamics)

First, the authors mapped out the landscape without the random noise. They proved that if you let the system settle, it will always find a valley. They also identified the "critical point"—the very top of the lowest hill separating two valleys.

Analogy: They found the exact "saddle point" on a mountain pass. It's the lowest point you must cross to get from one side of the mountain to the other. They proved this saddle point is stable and well-defined, even if the landscape is made of complex, wiggly fields.

2. The "Saddle" and the "Jump" (Spectral Data)

To calculate the jump time, you need to know the curvature of the hill.

At the bottom of the valley: The hill curves up gently.
At the top of the saddle: The hill curves down in one direction (the path forward) but curves up in all other directions (like a horse saddle).
The Innovation: The authors developed a rigorous way to measure the "steepness" of these curves for a wiggly jelly. They used a mathematical tool called Zeta-regularization.
- Simple Analogy: Imagine trying to count the infinite number of tiny ripples on a giant ocean wave to get a single number. Usually, this sum is infinite and breaks the math. The authors found a clever way to "cancel out" the infinite parts so that only the meaningful, finite numbers remain. This allows them to calculate the "prefactor" (the precise speed of the jump) without the math exploding.

3. Changing the Rules and Zooming In (Deformations and Convergence)

The paper also asks: "What happens if we change the temperature?" or "What happens if we look at this on a computer grid instead of a smooth continuum?"

The "Rubber Sheet" Analogy: Imagine the landscape is drawn on a rubber sheet. If you stretch or squish the sheet (changing physical parameters), do the valleys and hills disappear? The authors proved that, generally, the valleys and the saddle points just move around smoothly. They only disappear or merge in specific, predictable ways (called "catastrophes," like a fold or a cusp).
The "Pixel" Analogy: When we simulate this on a computer, we turn the smooth landscape into a grid of pixels. Does the answer change? The authors used a concept called $\Gamma$ -convergence to prove that as you make the pixels smaller and smaller (refining the grid), the computer's answer converges perfectly to the real, smooth answer. The jump time calculated on the computer will eventually match the real physics exactly.

The Final Result

The paper confirms that the famous Eyring–Kramers formula is not just a lucky guess for simple systems. It is a rigorous truth even for complex, infinite-dimensional fields like those found in nuclear physics.

The formula for the average time to jump is:
$\text{Time} \approx \text{(Shape Factor)} \times e^{\frac{\text{Height of Hill}}{\text{Noise Level}}}$

The authors proved that:

The Height of the Hill is the dominant factor (the exponential part).
The Shape Factor (the prefactor) can be calculated precisely using their new method of "cancelling infinities" (Zeta-regularization).
This result holds true even if you tweak the physical conditions or simulate it on a computer.

Summary

This paper builds a solid mathematical bridge between the messy, wiggly reality of complex fields (like those in QCD) and the clean, simple formulas physicists use to predict how long it takes for a system to change phases. It assures us that our "back-of-the-envelope" calculations are actually mathematically sound, even when dealing with the most complex, infinite-dimensional landscapes imaginable.

Technical Summary: Metastable Transitions and Γ–Convergent Eyring–Kramers Asymptotics in Landau–QCD Gradient Systems

Problem Setting
The paper addresses the rigorous mathematical analysis of metastable stochastic transitions in infinite-dimensional Landau-type gradient systems, motivated by QCD phenomenology. The system is defined on a bounded Lipschitz domain $\Omega \subset \mathbb{R}^d$ ( $1 \le d \le 3$ ) with an order parameter $\sigma$ governed by a free-energy functional:
$F[\sigma] = \int_\Omega \left( \frac{\kappa}{2} |\nabla \sigma|^2 + V(\sigma) \right) dx,$
where $V(\sigma)$ is a sextic potential ( $a_6 > 0$ ) capable of exhibiting bistability. The dynamics are modeled by an $L^2$ -gradient flow perturbed by a trace-class $Q$ -Wiener process. The central challenge is to establish the validity of the Eyring–Kramers (EK) law for the mean transition time between metastable wells in this infinite-dimensional setting, ensuring that the result remains robust under parameter deformations and variational discretizations.

Methodology
The authors construct a self-contained analytic framework integrating three primary mathematical pillars:

Variational and Spectral Analysis: The study utilizes the Gelfand triple $V \hookrightarrow H \hookrightarrow V^*$ (where $V=H^1(\Omega)$ and $H=L^2(\Omega)$ ). It employs the Euler–Lagrange map $E(\sigma) = -\kappa\Delta\sigma + V'(\sigma)$ and its Hessian $L_\sigma = -\kappa\Delta + V''(\sigma)$ . The analysis relies on the Lojasiewicz–Simon inequality to establish convergence rates to equilibria and spectral perturbation theory (Kato) to track eigenvalues.
Large Deviation and Potential Theory: For the stochastic dynamics, the authors apply Freidlin–Wentzell large deviation principles to identify the exponential rate of escape (the barrier height $\Delta F$ ). The prefactor is derived using potential theory and capacity arguments, specifically employing a Gaussian Laplace method adapted to the unstable direction at the saddle.
Γ–Convergence and Catastrophe Theory: To address stability under parameter changes and discretization, the paper utilizes Γ-convergence of the energy functionals and Mosco convergence of their subdifferentials. The structural stability of the critical points (minima and saddles) is analyzed through the lens of Thom's catastrophe theory (folds and cusps) and Cerf theory.

Key Contributions and Results

Deterministic Dynamics and Convergence: The paper establishes global well-posedness for the $L^2$ -gradient flow and proves that trajectories converge to stationary solutions. Under the assumption of analyticity for the potential $V$ , the Lojasiewicz–Simon inequality yields explicit convergence rates (exponential or algebraic depending on the spectral gap). The authors characterize the $\omega$ -limit sets and establish the existence of a mountain-pass critical nucleus (an index-one saddle) connecting metastable wells.
Spectral Geometry at the Saddle: A detailed analysis of the index-one saddle $\sigma^\dagger$ is provided. The authors construct local invariant manifolds (stable and unstable) and define the spectral data required for asymptotics: the unique negative eigenvalue $\lambda_-$ and the spectrum of the Hessian. Crucially, they introduce Agmon-type estimates to describe the spatial localization of the eigenfunctions and the tunneling direction.
Field-Theoretic Eyring–Kramers Asymptotics: The paper derives the EK law for the mean exit time $\mathbb{E}\tau$ :
$\mathbb{E}\tau \sim \frac{2\pi}{\Gamma |\lambda_-(\sigma^\dagger)|} \sqrt{\frac{\det(L_{\sigma_{meta}})}{\det'(L_{\sigma^\dagger})}} \exp\left(\frac{\Delta F}{\varepsilon^2}\right).$
Here, $\det'$ denotes the zeta-regularized determinant with the negative mode removed. The authors demonstrate that the ratio of determinants is ultraviolet-finite due to the cancellation of Weyl divergences between the well and the saddle, providing a rigorous justification for the prefactor in infinite dimensions.
Stability under Deformations and Discretization: The authors prove that the barrier height $\Delta F$ , the unstable eigenvalue $\lambda_-$ , and the determinant ratios vary continuously with smooth parameter changes, provided the parameters do not cross the discriminant set $\mathcal{D}$ (where the Hessian becomes degenerate). Crossing $\mathcal{D}$ results only in generic fold or cusp catastrophes. Furthermore, they show that these quantities persist and converge under conforming variational discretizations (Γ-convergence), ensuring that numerical approximations yield consistent EK parameters.

Significance and Claims
The paper claims to unify the classical intuition of Eyring, Kramers, and Langer with modern variational and spectral analysis. Its primary significance lies in providing a mathematically consistent foundation for metastable decay in Landau–QCD-type systems. Specifically:

It bridges microscopic fluctuations and macroscopic kinetics in infinite-dimensional spaces, clarifying which aspects of the EK law are universal (barrier dominance, index structure) versus model-dependent.
It resolves the issue of ultraviolet divergences in the prefactor through zeta-regularization, showing that the ratio of determinants is well-defined.
It establishes the structural stability of the EK asymptotics against both physical parameter deformations (e.g., temperature or chemical potential changes) and numerical discretization refinements.
The framework elucidates how the topology of the chiral–Polyakov potential governs phase conversion rates, capturing critical slowing down near spinodal points via universal scaling laws ( $\Delta F \sim (u-u_c)^{3/2}$ and $|\lambda_-| \sim (u-u_c)^{1/2}$ ).

The work does not propose new experimental setups but rather offers a rigorous theoretical tool for interpreting metastability in QCD-inspired effective models and related physical systems.