Skeleton as a probe of the cosmic web: the two-dimensional case

Abstract

We discuss the skeleton as a probe of the filamentary structures of a two-dimensional random field. It can be defined for a smooth field as the ensemble of pairs of field lines departing from saddle points, initially aligned with the major axis of local curvature and connecting them to local maxima. This definition is thus non-local and makes analytical predictions difficult, so we propose a local approximation: the local skeleton is given by the set of points where the gradient is aligned with the local curvature major axis and where the second component of the local curvature is negative. We perform a statistical analysis of the length of the total local skeleton, chosen for simplicity as the set of all points of space where the gradient is either parallel or orthogonal to the main curvature axis. In all our numerical experiments, which include Gaussian and various non-Gaussian realizations such as χ2 fields and Zel'dovich maps, the differential length f of the skeleton is found within a normalization factor to be very close to the probability distribution function (pdf) of the smoothed field, as expected and explicitly demonstrated in the Gaussian case where semi-analytical results are derived. Making the link with other works, we also show how the skeleton can be used to study the dynamics of large-scale structure.