Numerical challenges for energy conservation in N-body simulations of collapsing self-interacting dark matter halos

Abstract

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions. When the scattering can efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. A comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). We investigate the root of such problems, with a focus on energy non-conservation and discuss possible strategies to avoid them. We ran N-body simulations, with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing these challenges.

11 pages, 4 figures

References in corpus (26)

• Array Programming with NumPy
• Simulating cosmic structure formation with the GADGET-4 code
• An excess of small-scale gravitational lenses observed in galaxy clusters
• N-body simulations of gravitational dynamics
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• The Onset of Gravothermal Core Collapse in Velocity Dependent Self-Interacting Dark Matter Subhaloes
• Gravothermal evolution of dark matter halos with differential elastic scattering
• A Self-interacting Dark Matter Solution to the Extreme Diversity of Low-mass Halo Properties
• Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way
• Self-Interacting Dark Matter and Small-Scale Gravitational Lenses in Galaxy Clusters
• The Structure of Dissipative Dark Matter Halos
• N-body simulations of dark matter with frequent self-interactions
• The impact of baryonic potentials on the gravothermal evolution of self-interacting dark matter haloes
• The Cosmological Simulation Code OpenGadget3 -- Implementation of Meshless Finite Mass
• Improved strong lensing modelling of galaxy clusters using the Fundamental Plane: Detailed mapping of the baryonic and dark matter mass distribution of Abell S1063
• Dianoga SIDM: galaxy cluster self-interacting dark matter simulations
• Cosmological and idealized simulations of dark matter haloes with velocity-dependent, rare and frequent self-interactions
• Self-interacting dark matter interpretation of Crater II
• Numerical Challenges in Modeling Gravothermal Collapse in Self-Interacting Dark Matter Halos
• Towards time symmetric N-body integration
• Dynamical friction from self-interacting dark matter
• Properties and observables of massive galaxies in self-interacting dark matter cosmologies
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• Simulations of galaxy cluster mergers with velocity-dependent, rare and frequent self-interactions
• Constraints on Self-Interacting dark matter from relaxed galaxy groups
• Globally time-reversible fluid simulations with smoothed particle hydrodynamics

Cited by in corpus (11)

• Self-Interacting Dark Matter, Core Collapse and the Galaxy-Galaxy Strong Lensing Discrepancy
• Dynamical friction from self-interacting dark matter
• SIDM Concerto: Compilation and Data Release of Self-interacting Dark Matter Zoom-in Simulations
• SASHIMI-SIDM: Semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales
• COZMIC. III. Cosmological Zoom-in Simulations of Self-interacting Dark Matter with Suppressed Initial Conditions
• A Universal Analytic Model for Gravitational Lensing by Self-Interacting Dark Matter Halos
• N-body simulations of dark matter-baryon interactions
• Simulating realistic self-interacting dark matter models including small and large-angle scattering
• Core Collapse Beyond the Fluid Approximation: The Late Evolution of Self-Interacting Dark Matter Halos
• Constraining the Nature of Dark Matter from Tidal Radii of Cluster Galaxy Subhalos
• Accretion of self-interacting dark matter onto supermassive black holes