Phase Space dynamics of triaxial collapse: Joint density-velocity evolution

TL;DR

This paper develops a new phase space framework for understanding the joint evolution of density and velocity in triaxial collapse, revealing a universal relation between eigenvalues and analyzing their statistical distributions.

Contribution

It introduces a simplified eigenvalue-based phase space model for triaxial collapse and uncovers a universal relation between gravity and velocity eigenvalues across redshifts.

Findings

Derived new equations for eigenvalues in phase space.

Discovered a universal relation ${ ilde q}_d + { ilde q}_v=1$ valid over all redshifts.

Computed PDFs of density and velocity variables and compared with existing models.

Abstract

We investigate the dynamics of triaxial collapse in terms of eigenvalues of the deformation tensor, the velocity derivative tensor and the gravity Hessian. Using the Bond-Myers model of ellipsoidal collapse, we derive a new set of equations for the nine eigenvalues and examine their dynamics in phase space. The main advantage of this form is that it eliminates the complicated elliptic integrals that appear in the axes evolution equations and is more natural way to understand the interplay between the perturbations. This paper focuses on the density-velocity dynamics. The Zeldovich approximation implies that the three tensors are proportional; the proportionality constant is set by demanding `no decaying modes'. We extend this condition into the non-linear regime and find that the eigenvalues of the gravity Hessian and the velocity derivative tensor are related as ${\tilde q}_d + {\tilde q}_v=1$, where the triaxiality parameter ${\tilde q} = (\lambda_{\mathrm{max}} - \lambda_{\mathrm{inter}})/(\lambda_{\mathrm{max}} - \lambda_{\mathrm{min}})$. This is a {\it new universal relation} holding true over all redshifts and a range of mass scales to within a few percent accuracy. The mean density-velocity divergence relation at late times is close to linear, indicating that the dynamics is dictated by collapse along the largest eigendirection. This relation has a scatter, which we show, is intimately connected to the velocity shear. Finally, as an application, we compute the PDFs of the two variables and compare with other forms in the literature.