Core Collapse Beyond the Fluid Approximation: The Late Evolution of Self-Interacting Dark Matter Halos

We show that the gravothermal collapse of self-interacting dark matter (SIDM) halos can deviate from local thermodynamic equilibrium. As a consequence, the self-similar evolution predicted by the commonly adopted conducting fluid model can be altered or broken. Our results are obtained using a novel, efficient kinetic solver called kiss-sidm for tracing the gravothermal evolution based on the direct simulation Monte Carlo framework. In the long mean free path stage, the code is a viable alternative to the fluid model, yet requires no calibration parameters. Further, this method enables a fully kinetic treatment well into the late, short mean free path, stage of the collapse. We apply the method to a canonical case with isotropic, velocity independent scattering. We find that although a fluid treatment is appropriate deep in the short mean free path core, departures from local thermodynamic equilibrium develop in the intermediate mean free path region bounding the core, which modify the late-time evolution.