Decoupling hydrodynamization from thermalization via nonlinear Boltzmann equation

Abstract: The early thermalization puzzle arises from the unexpectedly early applicability of hydrodynamics in heavy-ion collisions. While hydrodynamics has traditionally been associated with the onset of local thermal equilibrium, its derivations -- whether microscopic or macroscopic -- rely instead on linearization around equilibrium. However, the linearization timescale -- the time at which a system's evolution begins to follow a linearized equation -- has not been systematically investigated. In this work, we employ the spectral nonlinear Boltzmann equation -- the lowest-order truncation of the spectral Bogoliubov--Born--Green--Kirkwood--Yvon (BBGKY) hierarchy -- to analyze the timescales of linearization and thermalization under three distinct truncation schemes. The first two truncations allow for analytic treatment via recursive spectral equations, while the third requires numerical methods for generic initial conditions. The analysis is performed for a homogeneous, massless system with a constant differential cross section. For this simplified setup, we find a robust separation: the linearization time is consistently about half the thermalization time ($\tau_{\mathrm{lin}}/\tau_{\mathrm{therm}} \approx 1/2$). This separation of timescales suggests an explanation for the early applicability of hydrodynamics and points toward a possible quantitative resolution of the early thermalization puzzle.