The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range 10(14) -10(15) M-circle dot. In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Lambda cold dark matter (Lambda CDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L, a form suitable for comparison with theoretical models. By proper scalingwe obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form f(E)(E)[1+ L-2/(2L(0)(2))](-beta infinity+beta infinity) L-2 beta 0. This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (beta(0) and beta(infinity)) and the scale of transition between them (L-0). The energy part f(E)(E) is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density Q = rho(r)/sigma(3)(r).