A Novel Implementation of Self-Interacting Dark Matter in AREPO

@misc{zier2026novel,
  title = {{A Novel Implementation of Self-Interacting Dark Matter in AREPO}},
  author = {Oliver Zier and Xuejian Shen and Vinh Tran and Martin Rosenlyst and Mark Vogelsberger and Rongrong Liu},
  howpublished = {arXiv},
  year = {2026},
  month = jul,
  doi = {10.48550/arxiv.2607.05504},
  eprint = {2607.05504},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2607.05504},
}

TY  - UNPB
TI  - A Novel Implementation of Self-Interacting Dark Matter in AREPO
AU  - Oliver Zier
AU  - Xuejian Shen
AU  - Vinh Tran
AU  - Martin Rosenlyst
AU  - Mark Vogelsberger
AU  - Rongrong Liu
PY  - 2026
DA  - 2026/07/06
PB  - arXiv
DO  - 10.48550/arxiv.2607.05504
UR  - https://arxiv.org/abs/2607.05504
AB  - Self-interacting dark matter (SIDM) influences halo structure through collisional heat transport and may solve several small-scale puzzles in structure formation. SIDM creates thermalized cores in low-mass haloes, which may account for the observed cored dwarf galaxies. During late-time gravothermal core collapse, SIDM can produce dense low-mass DM haloes and substructures detected through perturbations to cold stellar streams and strong gravitational lenses. In this work, we present a new Monte-Carlo SIDM implementation in the moving-mesh code AREPO-2, designed for efficiency, scalability, and extensibility. The central feature of the implementation is a dedicated DM-only neighbour-search tree that decouples the scattering solver from gravity. This preserves compatibility with the hierarchical time integration used by AREPO-2 while leaving the optimized gravity solver unconstrained. A pairwise communication scheme between MPI tasks allows tracking multiple scattering events in a single timestep while conserving momentum and energy and maintaining parallel consistency by construction. This is complemented by a per-pair timestep criterion that significantly reduces unnecessary timestep restrictions. The implementation natively supports velocity-dependent cross-sections and inelastic interactions, while a compact interface is designed for additional SIDM physics to be implemented without knowledge of the parallelization layer. We validate the implementation for isotropic, elastic scattering using a suite of idealized and cosmological tests. We assess performance and scalability in isolated core-collapse simulations and in cosmological boxes, both DM-only and with baryons. Except during the late stages of gravothermal collapse, SIDM simulations incur only modest overhead relative to the corresponding CDM runs and are substantially faster than the previous SIDM implementation in AREPO-1.
ER  -