Cusps in the center of galaxies: a real conflict with observations or a numerical artefact of cosmological simulations?

Abstract

Galaxy observations and N-body cosmological simulations produce conflicting dark matter halo density profiles for galaxy central regions. While simulations suggest a cuspy and universal density profile (UDP) of this region, the majority of observations favor variable profiles with a core in the center. In this paper, we investigate the convergency of standard N-body simulations, especially in the cusp region, following the approach proposed by \protect{\citep{13}}. We simulate the well known Hernquist model using the SPH code Gadget-3 and consider the full array of dynamical parameters of the particles. We find that, although the cuspy profile is stable, all integrals of motion characterizing individual particles suffer strong unphysical variations along the whole halo, revealing an effective interaction between the test bodies. This result casts doubts on the reliability of the velocity distribution function obtained in the simulations. Moreover, we find unphysical Fokker-Planck streams of particles in the cusp region. The same streams should appear in cosmological N-body simulations, being strong enough to change the shape of the cusp or even to create it. Our analysis, based on the Hernquist model and the standard SPH code, strongly suggests that the UDPs generally found by the cosmological N-body simulations may be a consequence of numerical effects. A much better understanding of the N-body simulation convergency is necessary before a 'core-cusp problem' can properly be used to question the validity of the CDM model.

arXiv comment7 pages, 5 figures, accepted to JCAP

References in corpus (24)

• The cosmological simulation code GADGET-2
• The Aquarius Project: the subhalos of galactic halos
• GADGET: A code for collisionless and gasdynamical cosmological simulations
• The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
• Resolving Cosmic Structure Formation with the Millennium-II Simulation
• The Diversity and Similarity of Simulated Cold Dark Matter Halos
• A Method for Measuring (Slopes of) the Mass Profiles of Dwarf Spheroidal Galaxies
• Cuspy No More: How Outflows Affect the Central Dark Matter and Baryon Distribution in Lambda CDM Galaxies
• H alpha rotation curves of Low Surface Brightness galaxies
• The unexpected diversity of dwarf galaxy rotation curves
• Dark and luminous matter in THINGS dwarf galaxies
• The Structural Evolution of Substructure
• Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo
• High-resolution rotation curves of low surface brightness galaxies: Data
• The SPLASH Survey: Spectroscopy of 15 M31 Dwarf Spheroidal Satellite Galaxies
• Improved Modeling of the Mass Distribution of Disk Galaxies by the Einasto Halo Model
• The Cusp/Core problem: supernovae feedback versus the baryonic clumps and dynamical friction model
• On the Origin of Cusps in Stellar Systems
• Milking the spherical cow: on aspherical dynamics in spherical coordinates
• Spectral analysis of the gamma-ray background near the dwarf Milky Way satellite Segue 1: Improved limits on the cross section of neutralino dark matter annihilation
• The real and apparent convergence of N-body simulations of the dark matter structures: is the Navarro-Frenk-White profile real?
• Orbital Divergence and Relaxation in the Gravitational N-Body Problem
• Relaxation of dark matter halos: how to match observational data?
• Overmerging in N-body simulations

Cited by in corpus (7)

• The $R^2$-Higgs inflation: $R^3$ contribution and preheating after ACT and SPT data
• The Next Generation Event Horizon Telescope Collaboration: History, Philosophy, and Culture
• The radius, at which a galaxy group stops the Hubble stream, and the group mass: an exact analytical solution
• Inner structure of dark matter halos at high $z$ in cosmological models with non-power-law primordial spectra
• The Hubble stream near a massive object: the exact analytical solution for the spherically-symmetric case
• Why does Einasto profile index $n\sim 6$ occur so frequently?
• The central region of a void: an analytical solution