Numerical Evolution of Self-Gravitating Halos of Self-Interacting Dark Matter

We discuss a modification of a recently developed numerical scheme for evolving spherically symmetric self-gravitating systems to include the effects of self-interacting dark matter. The approach is far more efficient than traditional N-body simulations and cross sections with different dependencies on velocity and scattering angle are easily accommodated. To demonstrate, we provide results of a simulation, which runs quickly on a personal computer, that shows the expected initial flattening of the inner region of a Navarro-Frenk-White (NFW) halo as well as the later gravothermal collapse instability that leads to a dense core at the galactic center. We note that this approach can also be used, with some augmentation, to simulate the dynamics in globular clusters by modeling gravitational hard scattering as a self-interaction.