Numerical simulations of dark matter haloes produce polytropic central cores when reaching thermodynamic equilibrium

Abstract

Self-gravitating astronomical objects often show a central plateau in the density profile (core) whose physical origin is hotly debated. Cores are theoretically expected in N-body systems of maximum entropy, however, they are not present in the canonical N-body numerical simulations of cold dark matter (CDM). Our work shows that despite this apparent contradiction between theory and numerical simulations, they are fully consistent. Simply put, cores are characteristic of systems in thermodynamic equilibrium, but thermalizing collisions are purposely suppressed in CDM simulations. When collisions are allowed, N-body numerical simulations develop cored density profiles, in perfect agreement with the theoretical expectation. We compare theory and two types of numerical simulations: (1) when DM particles are self-interacting (SIDM) with enough cross-section, then the effective two-body relaxation timescale becomes shorter than the Hubble time resulting in cored DM haloes. The haloes thus obtained, with masses from dwarf galaxies to galaxy clusters, collapse to a single shape after normalization, and this shape agrees with the polytropic density profile theoretically expected. (2) The inner radii in canonical N-body numerical simulations are always discarded because the use of finite-mass DM particles artificially increases the two-body collision rate. We show that the discarded radii develop cores that are larger than the employed numerical softening and have polytropic shapes independently of halo mass. Our work suggests that the presence of cores in simulated (or observed) density profiles can used as evidence for systems in thermodynamic equilibrium.

Accepted for publications in MNRAS. 6 figures, 9 pages. Paper that complements arXiv:2009.08994

References in corpus (10)

• The NumPy array: a structure for efficient numerical computation
• Constraints on the Self-Interaction Cross-Section of Dark Matter from Numerical Simulations of the Merging Galaxy Cluster 1E 0657-5
• High-resolution mass models of dwarf galaxies from LITTLE THINGS
• A physically motivated definition for the size of galaxies in an era of ultra-deep imaging
• SIDM on FIRE: Hydrodynamical Self-Interacting Dark Matter simulations of low-mass dwarf galaxies
• Connecting the structure of dark matter haloes to the primordial power spectrum
• Gaussian Random Field: Physical Origin of Sersic Profiles
• On the dearth of ultra-faint extremely metal poor galaxies
• Gaussian random field power spectrum and the Sérsic law
• Statistical mechanics of gravitating systems: An Overview

Cited by in corpus (4)

• Galaxy rotation curves and universal scaling relations: comparison between phenomenological and fermionic dark matter profiles
• Application of the Eddington inversion method to constrain the dark matter halo of galaxies using only observed surface brightness profiles
• Thermodynamics of self-gravitating fermions as a robust theory for dark matter halos: Stability analysis applied to the Milky Way
• Physically motivated fit to mass surface density profiles observed in galaxies