Quark production in the bottom-up thermalization

We investigate the impact of quark production on bottom-up thermalization in heavy-ion collisions. First, we extend the parametric estimates of bottom-up thermalization in pure gluon systems by incorporating quark production in the weak-coupling (high-energy) limit. Our analysis reveals that quark production does not alter the qualitative features of the three-stage thermalization process in this limit. Furthermore, we obtain the scaling behavior of the quark number density over time at each stage. Then, by solving the Boltzmann equation in diffusion approximation (BEDA) for longitudinally boost-invariant systems, we demonstrate how our detailed numerical simulations approach the predicted three-stage thermalization picture as the strong coupling $\alpha_s$ decreases. Finally, we carry out a detailed comparison of our BEDA results with those obtained by solving the QCD effective kinetic theory for intermediate values of $\alpha_s$, observing remarkably good quantitative agreement between the two approaches.