Relativistic particle simulation of a coloured parton plasma

We propose a relativistic particle-in-cell (PIC) simulation method for investigating the non-perturbative dynamics and collective interactions in a classical coloured parton plasma. The method has the advantage of including quarks and both short and long scale gluons non-perturbatively, and may be used for the study of a number of conceptual questions relevant to evolution of a quark gluon plasma formed in ultrarelativistic heavy ion collisions. We illustrate the simulation method by a detailed study of the thermalization of a highly nonequilibrium coloured parton plasma due to collective non-abelian interactions.