We study the thermalization process in highly occupied non-Abelian plasmas at weak coupling. The non-equilibrium 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~\cite{Berges:2013eia} along with novel findings. We demonstrate the emergence of universal attractor solutions, which govern the non-equilibrium evolution on large time scales both for non-expanding and expanding non-Abelian plasmas. The turbulent attractor for a non-expanding 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~\cite{Baier:2000sb}. While the focus of this paper is to understand the non-equilibrium 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.Comment: 49 pages, 21 figure

No abstract

The nonequilibrium evolution of heavy-ion collisions is studied in the limit of weak coupling at very high energy employing lattice simulations of the classical Yang-Mills equations. Performing the largest classical-statistical simulations to date, we find that the dynamics of the longitudinally expanding plasma becomes independent of the details of the initial conditions. After a transient regime dominated by plasma instabilities and free streaming, the subsequent space-time evolution is governed by a nonthermal fixed point, where the system exhibits the self-similar dynamics characteristic of wave turbulence. This allows us to distinguish between different kinetic scenarios in the classical regime. Within the accuracy of our simulations, the scaling behavior found is consistent with the "bottom-up" thermalization scenario

Isolated quantum systems in extreme conditions can exhibit unusually large occupancies per mode. This over-population gives rise to new universality classes of many-body systems far from equilibrium. We present theoretical evidence that important aspects of non-Abelian plasmas in the ultra-relativistic limit admit a dual description in terms of a Bose condensed scalar field theory.

We study the spectral properties of a highly occupied non-Abelian nonequilibrium plasma appearing ubiquitously in weak coupling descriptions of QCD matter. The spectral function of this far-fromequilibrium plasma is measured by employing linear response theory in classical-statistical real-time lattice Yang-Mills simulations. We establish the existence of transversely and longitudinally polarized quasiparticles and obtain their dispersion relations, effective mass, plasmon frequency, damping rate and further structures in the spectral and statistical functions. Our new method can be interpreted as a nonperturbative generalization of hard thermal loop (HTL) effective theory. We see indications that our results approach leading order HTL in the appropriate limit. The method can also be employed beyond the range of validity of HTL.