Shell crossing

Shell crossing occurs when the Lagrangian-to-Eulerian map ceases to be one-to-one. Its Jacobian determinant vanishes, the pressureless single-stream density formally diverges, and several streams subsequently occupy the same position.

For the Zeldovich approximation $\mathbf{x}=\mathbf{q}-b{\nabla }_q\psi$, let ${\lambda }_i$ be the eigenvalues of the potential Hessian, with this sign convention. Before shell crossing, mass conservation gives the displayed density. The first positive eigenvalue reaches $b{\lambda }_i=1$ and collapses one direction, creating a sheetlike cosmological caustic. Extrapolating straight trajectories after crossing lets particles pass through without self-consistent gravitational trapping, so the predicted multistream sheet spreads excessively. Collisionless gravitational evolution keeps matter more tightly bound; the linear displacement extrapolation cannot determine its late thickness reliably.

A cosmological caustic is the image of a singular point of the Lagrangian map. The ideal collisionless density diverges there, although velocity dispersion and other small-scale physics regularize the singularity.