Dynamics-based halo model for large scale structure

TL;DR

This paper introduces a physically motivated halo model based on phase-space decomposition, accurately describing large-scale structure by distinguishing orbiting and infalling dark matter particles, achieving high precision in correlation functions.

Contribution

The paper proposes a new halo definition based on phase-space components, improving the modeling of the halo--mass correlation function and power spectrum with high accuracy.

Findings

The distribution of orbiting particles is finite and characterized by a single scale, the halo radius.

The two-halo term can be modeled empirically using the Zel'dovich correlation function.

The model achieves approximately 2% accuracy in correlation functions and comparable accuracy in the power spectrum.

Abstract

Accurate modelling of the one-to-two halo transition has long been difficult to achieve. We demonstrate that physically motivated halo definitions that respect the bimodal phase-space distribution of dark matter particles near halos resolves this difficulty. Specifically, the two phase-space components are overlapping and correspond to: 1) particles \it orbiting \rm the halo; and 2) particles \it infalling \rm into the halo for the first time. Motivated by this decomposition, Garc\'ia [R. Garc\'ia et. al., MNRAS 521, 2464 (2023)] advocated for defining haloes as the collection of particles orbiting their self-generated potential. This definition identifies the traditional one-halo term of the halo--mass correlation function with the distribution of orbiting particles around a halo, while the two-halo term governs the distribution of infalling particles. We use dark matter simulations to demonstrate that the distribution of orbiting particles is finite and can be characterised by a single physical scale $r_{\rm h}$, which we refer to as the \it halo radius. \rm The two-halo term is described using a simple yet accurate empirical model based on the Zel'dovich correlation function. We further demonstrate that the halo radius imprints itself on the distribution of infalling particles at small scales. Our final model for the halo--mass correlation function is accurate at the $\approx 2\%$ level for $r \in [0.1,50]\ h^{-1}\ Mpc$. The Fourier transform of our best fit model describes the halo--mass power spectrum with comparable accuracy for $k\in [0.06, 6.0]\ h\ Mpc^{-1}$.