Supercooled Goldstone Bosons at the QCD Chiral Phase Transition

We discuss a universal non-equilibrium enhancement of long-wavelength Goldstone bosons induced by quenches to the broken phase in Model G – the dynamical universality class of an $O(4)$-antiferromagnet and the chiral phase transition in QCD. Scaling arguments for the coarsening dynamics describing the formation of the chiral condensate predict a parametric enhancement in the infrared spectra of Goldstones, a prediction confirmed by stochastic simulations of the transition. The details of the enhancement are determined by the non-linear dynamics of a superfluid effective theory, which is a limit of Model G reflecting the broken $O(4)$ symmetry. Our results translate to a parametric enhancement of low-momentum pions in heavy-ion collisions at the LHC, which are underpredicted in current hydrodynamic models without critical dynamics.

💡 Research Summary

This paper presents a comprehensive theoretical and numerical study of a universal non-equilibrium phenomenon associated with the chiral phase transition in QCD. The central focus is on the dynamics following a rapid quench from the symmetric to the broken phase within “Model G,” the dynamical universality class of the O(4) antiferromagnet, which describes the critical dynamics of QCD’s chiral transition in the limit of massless quarks.

The authors begin by contextualizing their work within the phenomenology of heavy-ion collisions, where the Quark-Gluon Plasma (QGP) cools through the chiral crossover. They note a persistent tension between hydrodynamic models and experimental data on low-momentum pion yields, suggesting a missing element related to critical dynamics. The study aims to bridge this gap by investigating the non-equilibrium evolution of the chiral condensate and its associated Goldstone bosons (pions) after a sudden quench.

The analysis is built upon Model G, which couples an O(4) order parameter field (representing the chiral condensate) to conserved vector and axial charges. A key insight is that in the broken phase, the reversible couplings (Poisson brackets) in Model G cause the Goldstone modes to propagate ballistically over large distances, unlike in purely dissipative models. This leads to a hierarchy of timescales near the critical point: a short relaxation time (τ_R), a much longer ballistic coarsening time (t_B ~ L/v), and an even longer diffusive equilibration time (t_D ~ L²/D_A).

The paper provides a detailed physical picture of the post-quench evolution. Immediately after the quench, local domains of the chiral condensate form on a timescale ~τ_R. Subsequently, on the ballistic timescale t_B, these domains grow and merge through a coarsening process governed by the ideal superfluid limit of Model G—a hydrodynamic theory valid when dissipation is negligible. It is during this extended coarsening stage that significant non-equilibrium effects manifest.

Using dynamic scaling arguments, the authors derive universal scaling forms for the time evolution of the global magnetization (a proxy for the condensate) and the pion two-point correlation function G_ππ(t, k). For the pion spectrum at long wavelengths (kξ « 1), the scaling predicts a parametric enhancement over the equilibrium value during the coarsening regime (τ_R « t ~ t_B). Specifically, in three dimensions, the ratio scales as G_ππ(t,k) / G_ππ_eq(k) ~ 1/(kξ) for massless pions, and ~1/(mξ) when a small explicit symmetry-breaking mass (m) is included. This signifies a substantial amplification of low-momentum pion modes.

These theoretical predictions are rigorously tested through large-scale stochastic simulations of instantaneous quenches in Model G. The numerical results confirm the scaling form for the growth of the magnetization, showing excellent data collapse when plotted against the scaled time vt/L. Furthermore, the simulated pion spectrum explicitly demonstrates the predicted infrared enhancement, following a ~1/k trend at low momentum, in stark contrast to the equilibrium ~1/k² behavior.

The study concludes that the coarsening dynamics following a quench through the chiral critical region inherently produces a surplus of soft pions. This mechanism provides a natural, first-principles explanation for the underprediction of low-momentum pion yields in conventional hydrodynamic models of heavy-ion collisions. The work establishes superfluid-like coarsening as a crucial non-equilibrium signature of QCD’s chiral phase transition, offering a new perspective for connecting critical dynamics with experimental observables.