Gravothermal collapse of isolated self-interacting dark matter haloes: N-body simulation versus the fluid model

TL;DR

This study compares N-body simulations and fluid models of self-interacting dark matter haloes, finding good agreement for isolated cases and providing insights into gravothermal collapse and the validity of modeling approaches.

Contribution

It demonstrates that N-body and fluid models agree within 20% for isolated SIDM haloes, clarifying previous discrepancies in cosmological simulations.

Findings

N-body and fluid models agree within 20% for isolated haloes.

Calibration of heat conduction improves model consistency.

Discrepancies in cosmological contexts are not due to model breakdown.

Abstract

Self-Interacting Dark Matter (SIDM) is a collisional form of cold dark matter (CDM), originally proposed to solve problems that arose when the collisionless CDM theory of structure formation was compared with observations of galaxies on small scales. The quantitative impact of the proposed elastic collisions on structure formation has been estimated previously by Monte Carlo N-body simulations and by a conducting fluid model, with apparently diverging results. To improve this situation, we make direct comparisons between new Monte Carlo N-body simulations and solutions of the conducting fluid model, for isolated SIDM haloes of fixed mass. This allows us to separate cleanly the effects of gravothermal relaxation from those of continuous mass accretion in an expanding background universe. When these two methods were previously applied to halo formation with cosmological boundary conditions, they disagreed by an order of magnitude about the size of the scattering cross section required to solve the so-called 'cusp-core problem.' We show here, however, that the methods agree with each other within 20 per cent for isolated haloes. This suggests that the two methods are consistent, and that their disagreement for cosmological haloes is not caused by a breakdown of their validity. The isolated haloes studied here undergo gravothermal collapse. We compare the solutions calculated by these two methods for gravothermal collapse starting from several initial conditions. This allows us to calibrate the heat conduction which accounts for the effect of elastic hard-sphere scattering in the fluid model. The amount of tuning of the thermal conductivity parameters required to bring the two methods into close agreement for isolated haloes, however, is too small to explain the discrepancy found previously in the cosmological context.