First-Principles Formalism for Simulating Self-Interacting Dark Matter

It is plausible that the dark matter particles have non-gravitational interactions among themselves. If such self interactions are large enough, they could leave an imprint on the morphology of galaxies. These effects can be studied with numerical simulations, which serve as the primary tool to predict the non-linear evolution of galactic structure. A standard assumption is that the course-grained phase-space distribution of the macroscopic simulation particles follows the same evolution equation as that of the fundamental dark matter particles. This Letter tests this assumption directly for the case of frequent dark matter scatterings, demonstrating that this is not generically true. Specifically, we develop a first-principles map from a microscopic particle physics description of self-interacting dark matter to a representation of macroscopic simulation particles for theories in the short-mean-free-path regime. Using this procedure, we show the emergence of an effective force between the simulation particles and derive their interaction cross section, which depends on the one from fundamental particle physics. This work provides the first explicit map from particle physics to simulation, which will facilitate exploring the phenomenological implications for galactic dynamics.

💡 Research Summary

This paper addresses a fundamental assumption underlying most self‑interacting dark matter (SIDM) simulations: that the coarse‑grained phase‑space distribution of macroscopic simulation particles evolves according to the same Boltzmann equation as the underlying microscopic dark‑matter particles. The authors demonstrate that this assumption fails in the short‑mean‑free‑path (SMFP) regime, where dark‑matter scatterings are frequent (Knudsen number Kn ≲ 1).

Starting from a concrete particle‑physics model—spin‑½ dark matter interacting via a Yukawa potential mediated by a light vector boson—the authors write down the velocity‑dependent differential cross‑section in the Born limit. They then map the microscopic parameters (mass m_dm, coupling α_dm, mediator mass m_V) onto a macroscopic description appropriate for N‑body simulations.

The key methodological steps are: (1) partition the halo’s phase space into cells (A, B, …) each described by a Maxwell‑Boltzmann distribution f_A^(0); (2) perform a Chapman‑Enskog expansion of the Boltzmann equation in powers of the Knudsen number, keeping terms up to first order; (3) coarse‑grain the resulting equations using position and velocity kernel functions, thereby identifying each cell with a simulation macro‑particle characterized by a position x_A and bulk velocity u_A.

Through this procedure the authors derive a coarse‑grained Boltzmann equation (Eq. 3) that contains an explicit source term representing a collective force F_AB acting on macro‑particle A due to collisions with macro‑particle B. The force takes the form

 F_i^{AB} = ρ_B (u_i^B − u_i^A) β^{1/2} (8 √(2π)/3) ∫ dcosθ (dσ_T/dcosθ),

where ρ_B is the mass density of particle B, β encodes the coarse‑grained velocity dispersion, and σ_T is an effective transfer cross‑section that emerges from the microscopic scattering kernel. Importantly, σ_T is not simply the microscopic σ_0; it is renormalized by the local velocity distribution and the SMFP expansion, as shown in Eq. 5.

The authors also express the Knudsen number in terms of this effective cross‑section (Eq. 4), linking the macroscopic fluid description back to the underlying particle physics. They illustrate the parameter space where the SMFP condition holds (Fig. 1) for halos of various virial masses, highlighting that realistic galaxy groups and clusters can indeed reside in this regime for plausible Yukawa‑type models.

By comparing their derived effective force with the “point‑particle scattering” prescription used in existing SIDM simulations, the paper shows that in the SMFP regime the dynamics are better captured by a continuous, friction‑like interaction rather than discrete collisions. This has direct implications for the thermal conduction that drives core expansion and subsequent gravothermal collapse, potentially altering predictions for core sizes, density profiles, and collapse timescales.

In conclusion, the work provides the first first‑principles mapping from microscopic SIDM models to the effective interactions employed in N‑body simulations for the frequent‑scattering limit. It clarifies the domain of validity of common simulation assumptions, offers a concrete formula for the effective cross‑section and force, and opens the path for more accurate modeling of SIDM effects on galactic structure. Future directions include extending the formalism to intermediate Knudsen numbers, incorporating anisotropic velocity distributions, and confronting the predictions with observational data from dwarf galaxies, groups, and clusters.