We have approached the dynamics of infinite self-gravitating systems by reducing as much as possible the complexity of the analogous cosmological problem and focusing on the essential aspects of the problem. We have shown that in simplifying the problem we have not lost some essential elements which change the nature of gravitational clustering. We have developed analytical and numerical approaches to understand the extent to which numerical simulations of a finite number of particles, reproduce the mean-field/Vlasov limit of the cosmological models, i.e. the validity of this “collisionless” limit. In addtion, we considered the understanding of “halo” structures, which are the primary building blocks in terms of which the non-linear structures observed in cosmological simulations are described. Finally we have developed several tests to constrain the accuracy and the resolution of cosmological simulations.

Gravitational clustering of mass structures is a well-posed problem of out-of-equilibrium statistical mechanics represents a paradigmatic problem for the statistical physics of long-range interacting systems finding application in many different areas of astrophysics (globular clusters, galaxies) and cosmology (gravitational clustering in the expanding universe). Instead of relaxing to thermodynamical equilibria through two-body collisions as short-range interacting systems, these systems reach, driven by a mean-field collisionless relaxation dynamics quasi-stationary state (QSS), whose lifetime diverges with the number of particles.
