# Power law pseudo phase-space density profiles of dark matter halos: fluke of physics?

## TL;DR

This study investigates why dark matter halos exhibit a power-law pseudo phase-space density profile, concluding that this scale-free behavior is likely coincidental rather than fundamental to gravitational physics.

## Contribution

The paper demonstrates that the observed power-law behavior of pseudo phase-space density profiles is probably a coincidence, challenging the idea that it is an inherent property of dark matter halos.

## Key findings

- Distribution of solutions is inconsistent with a fundamental scale-free property.
- Scale-free behavior likely results from a coincidence, not physics.
- Analysis of equilibrium halos suggests no fundamental origin for the power law.

## Abstract

It has been known for nearly 20 years that the pseudo phase-space density profile of equilibrium simulated dark matter halos, $\rho(r)/\sigma^3(r)$, is well described by a power law over 3 decades in radius, even though both the density $\rho(r)$, and the velocity dispersion $\sigma(r)$ deviate significantly from power laws. The origin of this scale-free behavior is not understood. It could be an inherent property of self-gravitating collisionless systems, or it could be a mere coincidence. To address the question we work with equilibrium halos, and more specifically, the second derivative of the Jeans equation, which, under the assumptions of (i) Einasto density profile, (ii) linear velocity anisotropy - density slope relation, and (iii) $\rho/\sigma^3\propto r^{-\alpha}$, can be transformed from a differential equation to a cubic algebraic equation. Relations (i)-(iii) are all observed in numerical simulations, and are well parametrized by a total of 4 or 6 model parameters. We do not consider dynamical evolution of halos; instead, taking advantage of the fact that the algebraic Jeans equation for equilibrium halos puts relations (i)-(iii) on the same footing, we study the (approximate) solutions of this equation in the 4 and 6 dimensional spaces. We argue that the distribution of best solutions in these parameter spaces is inconsistent with $\rho/\sigma^3\propto r^{-\alpha}$ being an fundamental property of gravitational evolution, and conclude that the scale-free nature of this quantity is likely to be a fluke.