Caustic Skeleton & Cosmic Web

We present a general formalism for identifying the caustic structure of a dynamically evolving mass distribution, in an arbitrary dimensional space. The identification of caustics in fluids with Hamiltonian dynamics, viewed in Lagrangian space, corresponds to the classification of singularities in Lagrangian catastrophe theory. On the basis of this formalism we develop a theoretical framework for the dynamics of the formation of the cosmic web, and specifically those aspects that characterize its unique nature: its complex topological connectivity and multiscale spinal structure of sheetlike membranes, elongated filaments and compact cluster nodes. Given the collisionless nature of the gravitationally dominant dark matter component in the universe, the presented formalism entails an accurate description of the spatial organization of matter resulting from the gravitationally driven formation of cosmic structure. The present work represents a significant extension of the work by Arnol'd et al. [1], who classified the caustics that develop in one- and two-dimensional systems that evolve according to the Zel'dovich approximation. His seminal work established the defining role of emerging singularities in the formation of nonlinear structures in the universe. At the transition from the linear to nonlinear structure evolution, the first complex features emerge at locations where different fluid elements cross to establish multistream regions. Involving a complex folding of the 6-D sheetlike phase-space distribution, it manifests itself in the appearance of in finite density caustic features. The classification and characterization of these mass element foldings can be encapsulated in caustic conditions on the eigenvalue and eigenvector fields of the deformation tensor field. In this study we introduce an alternative and transparent proof for Lagrangian catastrophe theory. This facilitates the derivation of the caustic conditions for general Lagrangian fluids, with arbitrary dynamics. Most important in the present context is that it allows us to follow and describe the full three-dimensional geometric and topological complexity of the purely gravitationally evolving nonlinear cosmic matter field. While generic and statistical results can be based on the eigenvalue characteristics, one of our key findings is that of the significance of the eigen vector field of the deformation field for outlining the entire spatial structure of the caustic skeleton emerging from a primordial density field. In this paper we explicitly consider the caustic conditions for the three-dimensional Zel'dovich approximation, extending earlier work on those for one- and two-dimensional fluids towards the full spatial richness of the cosmic web. In an accompanying publication, we apply this towards a full three-dimensional study of caustics in the formation of the cosmic web and evaluate in how far it manages to outline and identify the intricate skeletal features in the corresponding N-body simulations.