Numerical simulation of non-Abelian particle-field dynamics

Numerical 1D-3V solutions of the Wong-Yang-Mills equations with anisotropic particle momentum distributions are presented. They confirm the existence of plasma instabilities for weak initial fields and of their saturation at a level where the particle motion is affected, similar to Abelian plasmas. The isotropization of the particle momenta by strong random fields is shown explicitly, as well as their nearly exponential distribution up to a typical hard scale, which arises from scattering off field fluctuations. By variation of the lattice spacing we show that the effects described here are independent of the UV field modes near the end of the Brioullin zone.