Quark production and thermalization of the quark-gluon plasma

We first assemble a full set of the Boltzmann Equation in Diffusion Approximation (BEDA) for studying thermalization/hydrodynamization as well as the production of massless quarks and antiquarks in out of equilibrium systems. In the BEDA, the time evolution of a generic system is characterized by the following space-time dependent quantities: the jet quenching parameter, the effective temperature, and two more for each quark flavor that describe the conversion between gluons and quarks/antiquarks via the $2\leftrightarrow2$ processes. Out of the latter two quantities, an effective net quark chemical potential is defined, which equals the net quark chemical potential after thermal equilibration. We then study thermalization and the production of three flavors of massless quarks and antiquarks in spatially homogeneous systems initially filled only with gluons. A parametric understanding of thermalization and quark production is obtained for either initially very dense or dilute systems, which are complemented by detailed numerical simulations for intermediate values of initial gluon occupancy $f_0$. For a wide range of $f_0$, the final equilibration time is determined to be about one order of magnitude longer than that in the corresponding pure gluon systems. Moreover, during the final stage of the thermalization process for $f_0\geq 10^{-4}$, gluons are found to thermalize earlier than quarks and antiquarks, undergoing the top-down thermalization.