Numerical Challenges in Modeling Gravothermal Collapse in Self-Interacting Dark Matter Halos

Abstract

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter~(SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions ($\sigma /m=50$~cm${}^2$/g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study ($10^6$ particles per halo), we find that variations of order $10\%$ in collapse time are still present. The results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

7 figures, 1 table

References in corpus (48)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• The cosmological simulation code GADGET-2
• The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• Dark Matter Self-interactions and Small Scale Structure
• Cosmological Simulations with Self-Interacting Dark Matter I: Constant Density Cores and Substructure
• Subhaloes in Self-Interacting Galactic Dark Matter Haloes
• Core Formation in Dwarf Halos with Self Interacting Dark Matter: No Fine-Tuning Necessary
• Dark Matter Substructure in Numerical Simulations: A Tale of Discreteness Noise, Runaway Instabilities, and Artificial Disruption
• What does the Bullet Cluster tell us about self-interacting dark matter?
• Mirror matter-type dark matter
• Can Feedback Solve the Too Big to Fail Problem?
• Gravothermal collapse of isolated self-interacting dark matter haloes: N-body simulation versus the fluid model
• Supermassive Black Holes from Ultra-Strongly Self-Interacting Dark Matter
• Constraining Dissipative Dark Matter Self-Interactions
• Constraining Velocity-dependent Self-Interacting Dark Matter with the Milky Way's dwarf spheroidal galaxies
• All about baryons: revisiting SIDM predictions at small halo masses
• Reconciling the Diversity and Uniformity of Galactic Rotation Curves with Self-Interacting Dark Matter
• Accelerated core collapse in tidally stripped self-interacting dark matter halos
• Too Big To Fail in Light of Gaia
• SIDM on FIRE: Hydrodynamical Self-Interacting Dark Matter simulations of low-mass dwarf galaxies
• Cosmic particle colliders: simulations of self-interacting dark matter with anisotropic scattering
• Signatures of Self-Interacting dark matter on cluster density profile and subhalo distributions
• Diversity in density profiles of self-interacting dark matter satellite halos
• Seeding Supermassive Black Holes with Self-Interacting Dark Matter: A Unified Scenario with Baryons
• Core-collapse, evaporation and tidal effects: the life story of a self-interacting dark matter subhalo
• The diverse dark matter density at sub-kiloparsec scales in Milky Way satellites:implications for the nature of dark matter
• The Onset of Gravothermal Core Collapse in Velocity Dependent Self-Interacting Dark Matter Subhaloes
• Signatures of Velocity-Dependent Dark Matter Self-Interactions in Milky Way-mass Halos
• Gravothermal evolution of dark matter halos with differential elastic scattering
• Gravothermal solutions of SIDM halos: mapping from constant to velocity-dependent cross section
• Early formation of supermassive black holes via dark matter self-interactions
• 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
• Massive Black Holes from Dissipative Dark Matter
• Universal gravothermal evolution of isolated self-interacting dark matter halos for velocity-dependent cross sections
• Motivations for a Large Self-Interacting Dark Matter Cross Section from Milky Way Satellites
• Multimass spherical structure models for N-body simulations
• A semi-analytic study of self-interacting dark-matter haloes with baryons
• The probability of galaxy-galaxy strong lensing events in hydrodynamical simulations of galaxy clusters
• 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
• A Parametric Model for Self-Interacting Dark Matter Halos
• Black hole formation in the context of dissipative dark matter
• Comparing Implementations of Self-Interacting Dark Matter in the Gizmo and Arepo Codes
• Subhalo effective density slope measurements from HST strong lensing data with neural likelihood-ratio estimation

Cited by in corpus (12)

• Dianoga SIDM: galaxy cluster self-interacting dark matter simulations
• Numerical challenges for energy conservation in N-body simulations of collapsing self-interacting dark matter halos
• 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
• Core Collapse Beyond the Fluid Approximation: The Late Evolution of Self-Interacting Dark Matter Halos
• DARKSKIES: A suite of super-sampled zoom-in simulations of galaxy clusters with self-interacting dark matter
• Semi-Analytic Modeling of Dark Matter Subhalo Encounters with Thin Stellar Streams: Statistical Predictions for GD-1-like Streams in CDM
• Core collapse in resonant self-interacting dark matter across two decades in halo mass
• Constraining the Nature of Dark Matter from Tidal Radii of Cluster Galaxy Subhalos
• A New Method to Simulate Dark Matter-Baryon Interactions and Application to an Isolated Disk Galaxy