Dark matter halo dynamics in 2D Vlasov Simulations: a self-similar approach

TL;DR

This study investigates the self-similar behavior of 2D dark matter halos using numerical simulations, finding that halos exhibit approximate self-similarity with density profiles close to theoretical predictions, while deviations are caused by relaxation and boundary effects.

Contribution

The paper adapts self-similar solutions to 2D CDM halo simulations, analyzing their validity and identifying factors causing deviations from ideal self-similarity.

Findings

Halos show good self-similar scaling after initial relaxation.

Density profiles follow $ ho o r^{-1}$, consistent with models.

Deviations linked to relaxation, boundary effects, and transverse motions.

Abstract

Understanding dark matter halo dynamics can be pivotal in unravelling the nature of dark matter particles. Analytical treatment of the multistream flows inside the turnaround region of a collapsed cold dark matter (CDM) halo using various self-similar approaches already exist. In this work, we aim to determine the extent of self-similarity in 2D halo dynamics and the factors leading to deviations from it by studying numerical simulations of monolithically growing CDM halos. We have adapted the Fillmore and Goldreich (FG) self-similar solutions assuming cylindrical symmetry to data from 2D Vlasov-Poisson (ColDICE package) simulations of CDM halos seeded by sine wave initial conditions. We measured trajectories in position and phase-space, mass and density profiles and compared these to predictions from the FG model. We find that after turn-around and subsequent shell crossing, particles undergo a period of relaxation, typically about 1-2 oscillations about the center before they start to trace the self-similar fits and continue to do so as long as their orbits are predominantly radial. Overplotting the trajectories from different snapshots in scale-free position-time and phase spaces shows strikingly good superposition, a defining feature of self-similarity. The radial density profiles measured from simulations: $\rho \propto r^{-\alpha}, \alpha = 0.9 - 1.0$ are consistent with FG's prediction $\alpha=1$ for 2D halos. Deviations from the model, on the other hand, are evidently linked to relaxation, inhibited motion due to periodic boundaries, transverse motion in the halo interior, and deficit of infalling mass in limited simulation volume. It could not be conclusively established if the halos tend to grow circular over time. Extension of this work to actual 3D CDM cosmologies necessitates further detailed study of self-similar solutions with ellipsoidal collapse and transverse motion.