Accurately simulating core-collapse self-interacting dark matter halos
This paper investigates the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations during the gravothermal collapse phase, demonstrating that high-resolution simulations with strict energy conservation can accurately model the formation of compact halos capable of explaining observations like the GD-1 stellar stream.