Nonequilibrium quark production in the expanding QCD plasma

We perform real-time lattice simulations of nonequilibrium quark production in the longitudinally expanding QCD plasma. Starting from a highly occupied gluonic state with vacuum quark sector, we extract the time evolution of quark and gluon number densities per unit transverse area and rapidity. The total quark number shows after an initial rapid increase an almost linear growth with time. Remarkably, this growth rate appears to be consistent with a simple kinetic theory estimate involving only two-to-two scattering processes in small-angle approximation. This extends previous findings about the role of two-to-two scatterings for purely gluonic dynamics in accordance with the early stages of the bottom-up thermalization scenario.