Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the $O(4)$ critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the non-equilibrium dynamics of pion fields. Strikingly, the non-equilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two point functions is determined by the non-linear dynamics of an ideal non-abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the non-equilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

💡 Research Summary

This paper presents a pioneering investigation into the out-of-equilibrium dynamics associated with the chiral phase transition of QCD, employing the framework of “Model G,” its dynamical universality class. The primary goal is to understand the real-time evolution of the chiral order parameter and its associated Goldstone bosons (pions) when the system is rapidly quenched from a high-temperature chirally symmetric phase to a low-temperature symmetry-broken phase.

The core methodology involves large-scale, three-dimensional stochastic (Langevin) simulations of an O(4)-symmetric scalar field theory coupled to conserved charges, which constitutes Model G. The simulations are performed near the critical point in the chiral limit (neglecting explicit symmetry breaking). The key scenario studied is an instantaneous quench, where the system’s parameters are suddenly changed to take it from above to below the critical temperature.

The main findings reveal a striking and complex non-equilibrium evolution of the pion two-point correlation functions, which deviates significantly from simple exponential relaxation. This evolution exhibits three distinct temporal regimes:

1. Initial Growth Phase: Following the quench, the chiral condensate begins to grow. This growth is accompanied by an amplification of soft pion modes.
2. Parametric Enhancement Phase: A period of parametrically large enhancement in the pion correlations emerges. The authors analytically show that this regime is governed by the non-linear, non-dissipative dynamics of an ideal non-abelian superfluid, which is a direct consequence of the spontaneously broken chiral symmetry and represents a specific limit of Model G.
3. Slow Relaxation Phase: Subsequently, the system undergoes a very slow relaxation back to equilibrium, a hallmark of critical slowing down near a second-order phase transition.

The study provides both a detailed scaling analysis based on dynamic critical exponents and a mean-field theory treatment for comparison. While mean-field theory captures some qualitative trends, it fails to quantitatively describe the strong enhancement and the specific dynamics captured in the full numerical simulations, underscoring the importance of critical fluctuations.

The authors directly connect these non-equilibrium field dynamics to a potential observable in heavy-ion collision experiments: a long-lived parametric enhancement in the yield of soft pions relative to thermal equilibrium expectations. This offers a novel potential explanation for the observed soft pion excess in data from experiments like ALICE at the LHC, suggesting that the rapid cooling of the Quark-Gluon Plasma through the chiral crossover could imprint non-equilibrium critical dynamics onto the final hadron spectra.

In summary, this work delivers the first comprehensive numerical and analytical study of out-of-equilibrium Model G dynamics. It establishes that a quench through the chiral phase transition leads to a characteristic, multi-stage evolution with a significant amplification of soft modes, linking fundamental critical dynamics to potential phenomenological signatures in high-energy nuclear collisions.