Universal attractor in a highly occupied non-Abelian plasma

We study the thermalization process in highly occupied non-Abelian plasmas at weak coupling. The nonequilibrium dynamics of such systems is classical in nature and can be simulated with real-time lattice gauge theory techniques. We provide a detailed discussion of this framework and elaborate on the results reported in J. Berges, K. Boguslavski, S. Schlichting, and R. Venugopalan, Phys. Rev. D 89, 074011 (2014) along with novel findings. We demonstrate the emergence of universal attractor solutions, which govern the nonequilibrium evolution on large time scales both for nonexpanding and expanding non-Abelian plasmas. The turbulent attractor for a nonexpanding plasma drives the system close to thermal equilibrium on a time scale $t\sim Q^{-1}{\alpha }_s^{-7/4}$. The attractor solution for an expanding non-Abelian plasma leads to a strongly interacting albeit highly anisotropic system at the transition to the low-occupancy or quantum regime. This evolution in the classical regime is, within the uncertainties of our simulations, consistent with the “bottom up” thermalization scenario [R. Baier, A. H. Mueller, D. Schiff, and D. T. Son, Phys. Lett. B 502, 51 (2001)]. While the focus of this paper is to understand the nonequilibrium dynamics in weak coupling asymptotics, we also discuss the relevance of our results for larger couplings in the early time dynamics of heavy ion collision experiments.