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 <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mrow> <a:mi>O</a:mi> <a:mo stretchy="false">(</a:mo> <a:mn>4</a:mn> <a:mo stretchy="false">)</a:mo> </a:mrow> </a:math> 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 nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium 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 nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.