Using N-body simulations with a large set of massless test particles we compare the predictions of two theories of violent relaxation, the well known Lynden-Bell theory and the more recent theory by Nakamura. We derive ``weaken'' versions of both theories in which we use the whole equilibrium coarse-grained distribution function as a constraint instead of the total energy constraint. We use these weaken theories to construct expressions for the conditional probability $K_i(\tau)$ that a test particle initially at the phase-space coordinate $\tau$ would end-up in the $i$'th macro-cell at equilibrium. We show that the logarithm of the ratio $R_{ij}(\tau) \equiv K_i(\tau)/K_j(\tau)$ is directly proportional to the initial phase-space density $f_0(\tau)$ for the Lynden-Bell theory and inversely proportional to $f_0(\tau)$ for the Nakamura theory. We then measure $R_{ij}(\tau)$ using a set of N-body simulations of a system undergoing a gravitational collapse to check the validity of the two theories of violent relaxation. We find that both theories are at odds with the numerical results, qualitatively and quantitatively.

We examine an inconsistency in theories of violent-relaxation by Lynden-Bell and Nakamura. The inconsistency arises from the non-transitive nature of these theories: a system that under goes a violent-relaxation, relaxes and then upon an addition of energy, undergoes violent-relaxation once again would settle in an equilibrium state that is different from the one that is predicted had the system would go directly from the initial to the final states. We conclude that a proper description of the violent-relaxation process cannot be achieved by equilibrium statistical mechanics approach, but instead a dynamical theory for the coarse-grained phase-space density is needed.

A method is presented for computing the 6-D phase-space density f(x,v) and its PDF v(f) in an N-body system. It is based on Delaunay tessellation, yielding v(f) with a fixed smoothing window over a wide f range, independent of the sampling resolution. It is found that in a gravitationally relaxed halo built by hierarchical clustering, v(f) is a robust power law, v(f) f^{-2.5 \pm 0.05}, over more than 4 decades in f, from its virial level to the current resolution limit. This is valid for halos of different sizes in the LCDM cosmology, indicating insensitivity to the initial-fluctuation power spectrum as long as the small-scale fluctuations were not completely suppressed. By mapping f in position space, we find that the high-f contributions to v(f) come from the "cold" subhalos within the parent halo rather than the halo central region and its global spherical profile. The f in subhalos near the halo virial radius is more than 100 times higher than at the halo center, and it decreases gradually with decreasing radius. This indicates phase mixing due to mergers and tidal effects involving puffing up and heating. The phase-space structure provides a sensitive tool for studying the evolution of subhalos during the buildup of halos. One wishes to understand why the substructure adds up to the universal power law in v(f). It seems that the f^{-2.5} behavior is related to the hierarchical clustering process and is not a general result of violent relaxation.

We present a method for computing the 6-dimensional coarse-grained phase-space density $f(x,v)$ in an N-body system, and derive its distribution function $v(f)$. The method is based on Delaunay tessellation, where $v(f)$ is obtained with an effective fixed smoothing window over a wide $f$ range. The errors are estimated, and $v(f)$ is found to be insensitive to the sampling resolution or the simulation technique. We find that in gravitationally relaxed haloes built by hierarchical clustering, $v(f)$ is well approximated by a robust power law, $v(f) \propto f^{-2.5 \pm 0.05}$, over more than 4 decades in $f$, from its virial level to the numerical resolution limit. This is tested to be valid in the $\Lambda$CDM cosmology for haloes with masses $10^9-10^{15}\msun$, indicating insensitivity to the slope of the initial fluctuation power spectrum. By mapping the phase-space density in position space, we find that the high-$f$ end of $v(f)$ is dominated by the ``cold'' subhaloes rather than the parent-halo central region and its global spherical profile. The value of $f$ in subhaloes near the virial radius is typically $>100$ times higher than its value at the halo centre, and it decreases gradually from outside in toward its value at the halo centre. This seems to reflect phase mixing due to mergers and tidal effects involving puffing up and heating. The phase-space density can thus provide a sensitive tool for studying the evolution of subhaloes during the hierarchical buildup of haloes. It remains to be understood why the evolved substructure adds up to the actual universal power law of $v(f) \propto f^{-2.5}$. It seems that this behaviour results from the hierarchical clustering process and is not a general result of violent relaxation.

We propose a model for how the buildup of dark halos by merging satellites produces a characteristic inner cusp, of a density profile \rho \prop r^-a with a -> a_as > 1, as seen in cosmological N-body simulations of hierarchical clustering scenarios. Dekel, Devor & Hetzroni (2003) argue that a flat core of a<1 exerts tidal compression which prevents local deposit of satellite material; the satellite sinks intact into the halo center thus causing a rapid steepening to a>1. Using merger N-body simulations, we learn that this cusp is stable under a sequence of mergers, and derive a practical tidal mass-transfer recipe in regions where the local slope of the halo profile is a>1. According to this recipe, the ratio of mean densities of halo and initial satellite within the tidal radius equals a given function psi(a), which is significantly smaller than unity (compared to being 1 according to crude resonance criteria) and is a decreasing function of a. This decrease makes the tidal mass transfer relatively more efficient at larger a, which means steepening when a is small and flattening when a is large, thus causing converges to a stable solution. Given this mass-transfer recipe, linear perturbation analysis, supported by toy simulations, shows that a sequence of cosmological mergers with homologous satellites slowly leads to a fixed-point cusp with an asymptotic slope a_as>1. The slope depends only weakly on the fluctuation power spectrum, in agreement with cosmological simulations. During a long interim period the profile has an NFW-like shape, with a cusp of 1<a<a_as. Thus, a cusp is enforced if enough compact satellite remnants make it intact into the inner halo. In order to maintain a flat core, satellites must be disrupted outside the core, possibly as a result of a modest puffing up due to baryonic feedback.

We show how the buildup of halos by merging satellites forces an inner cusp, with a density profile \rho \propto r^{-\alpha} where \alpha \to \alpha_a \gsim 1$. Our analysis is based on a new prescription for tidal stripping as a function of \alpha(r), using a simple toy model which matches N-body simulations. In a core of \alpha < 1 there is tidal compression rather than stripping and the satellites sink towards the halo center, causing a rapid steepening of the profile to \alpha > 1. Where \alpha > 1, the stripping of each satellite shell is preceded by gradual puffing up, which makes the stripping more efficient at larger \alpha, causing flattening where \alpha is large enough. Therefore, we can show using linear perturbation analysis that a sequence of mergers slowly leads to a fixed point \alpha(r)=\alpha_a. This result implies that a cusp is enforced as long as enough satellite material makes it into the inner halo and is deposited there. We conclude that in order to maintain a flat core, as indicated by observations, satellites must be disrupted outside the core, e.g., because of puffing up due to baryonic feedback effects.