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.

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

References in corpus (42)

• The NumPy array: a structure for efficient numerical computation
• Planck 2015 results. XIII. Cosmological parameters
• A Universal Density Profile from Hierarchical Clustering
• Observational evidence for self-interacting cold dark matter
• Toward a halo mass function for precision cosmology: the limits of universality
• Small-Scale Challenges to the $\Lambda$CDM Paradigm
• The Inner Structure of LambdaCDM Halos III: Universality and Asymptotic Slopes
• Dark Matter Self-interactions and Small Scale Structure
• The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
• Clumps and streams in the local dark matter distribution
• Constraints on the Self-Interaction Cross-Section of Dark Matter from Numerical Simulations of the Merging Galaxy Cluster 1E 0657-5
• Cosmological Simulations with Self-Interacting Dark Matter I: Constant Density Cores and Substructure
• Cold dark matter: controversies on small scales
• Earth-mass dark-matter haloes as the first structures in the early Universe
• High-resolution mass models of dwarf galaxies from LITTLE THINGS
• The properties and luminosity function of extremely low luminosity galaxies
• Halo Properties in Cosmological Simulations of Self-Interacting Cold Dark Matter
• Galaxy Groups in the SDSS DR4: III. the luminosity and stellar mass functions
• Core Formation in Dwarf Halos with Self Interacting Dark Matter: No Fine-Tuning Necessary
• Small scale problems of the $\Lambda$CDM model: a short review
• Universal structure of dark matter haloes over a mass range of 20 orders of magnitude
• Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons
• All about baryons: revisiting SIDM predictions at small halo masses
• Long-term Evolution of Stellar Self-Gravitating System away from the Thermal Equilibrium: connection with non-extensive statistics
• A physically motivated definition for the size of galaxies in an era of ultra-deep imaging
• Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes
• SIDM on FIRE: Hydrodynamical Self-Interacting Dark Matter simulations of low-mass dwarf galaxies
• The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles
• The structure and assembly history of cluster-size haloes in Self-Interacting Dark Matter
• On the formation and stability of fermionic dark matter halos in a cosmological framework
• Modelling mass distribution in elliptical galaxies: mass profiles and their correlation with velocity dispersion profiles
• Connecting the structure of dark matter haloes to the primordial power spectrum
• Power Law Stellar Distributions
• Constraining Non-thermal Dark Matter by CMB
• Towards an improved model of self-interacting dark matter haloes
• Time-dependent nonextensivity arising from the rotational evolution of solar-type stars
• On the dearth of ultra-faint extremely metal poor galaxies
• Gaussian Random Field: Physical Origin of Sersic Profiles
• Gaussian random field power spectrum and the Sérsic law
• The principle of maximum entropy explains the cores observed in the mass distribution of dwarf galaxies
• Statistical mechanics of gravitating systems: An Overview
• Polytropic dark halos of elliptical galaxies

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