We present an effective evolution equation for a coarse-grained distribution function of a long-range-interacting system preserving the symplectic structure of the non-collisional Boltzmann, or Vlasov, equation. We first derive a general form of such an equation based on symmetry considerations only. Then, we explicitly derive the equation for one-dimensional systems, finding that it has the form predicted on general grounds. Finally, we use such an equation to predict the dependence of the damping times on the coarse-graining scale and numerically check it for some one-dimensional models, including the Hamiltonian Mean Field (HMF) model, a scalar field with quartic interaction, a 1-d self-gravitating system, and the Self-Gravitating Ring (SGR).

We study the stability of a family of spherical equilibrium models of self-gravitating systems, the so-called $\gamma-$models with Osipkov-Merritt velocity anisotropy, by means of $N-$body simulations. In particular, we analyze the effect of self-consistent $N-$body chaos on the onset of radial-orbit instability (ROI). We find that degree of chaoticity of the system associated to its largest Lyapunov exponent $\Lambda_{\rm max}$ has no appreciable relation with the stability of the model for fixed density profile and different values of radial velocity anisotropy. However, by studying the distribution of the Lyapunov exponents $\lambda_{\rm m}$ of the individual particles in the single-particle phase space, we find that more anisotropic systems have a larger fraction of orbits with larger $\lambda_{\rm m}$.

We revisit the r\^{o}le of discreteness and chaos in the dynamics of self-gravitating systems by means of $N$-body simulations with active and frozen potentials, starting from spherically symmetric stationary states and considering the orbits of single particles in a frozen $N$-body potential as well as the orbits of the system in the full $6N$-dimensional phase space. We also consider the intermediate case where a test particle moves in the field generated by $N$ non-interacting particles, which in turn move in a static smooth potential. We investigate the dependence on $N$ and on the softening length of the largest Lyapunov exponent both of single particle orbits and of the full $N$-body system. For single orbits we also study the dependence on the angular momentum and on the energy. Our results confirm the expectation that orbital properties of single orbits in finite-$N$ systems approach those of orbits in smooth potentials in the continuum limit $N \to \infty$ and that the largest Lyapunov exponent of the full $N$-body system does decrease with $N$, for sufficiently large systems with finite softening length. However, single orbits in frozen models and active self-consistent models have different largest Lyapunov exponents and the $N$-dependence of the values in non-trivial, so that the use of frozen $N$-body potentials to gain information on large-$N$ systems or on the continuum limit may be misleading in certain cases.

Using direct $N$-body simulations of self-gravitating systems we study the dependence of dynamical chaos on the system size $N$. We find that the $N$-body chaos quantified in terms of the largest Lyapunov exponent $\Lambda_{\rm max}$ decreases with $N$. The values of its inverse (the so-called Lyapunov time $t_\lambda$) are found to be smaller than the two-body collisional relaxation time but larger than the typical violent relaxation time, thus suggesting the existence of another collective time scale connected to many-body chaos.

In $N$-body systems with long-range interactions mean-field effects dominate over binary interactions (collisions), so that relaxation to thermal equilibrium occurs on time scales that grow with $N$, diverging in the $N\to\infty$ limit. However, a faster and non-collisional relaxation process, referred to as violent relaxation, sets in when starting from generic initial conditions: collective oscillations (referred to as virial oscillations) develop and damp out on timescales not depending on the system's size. After the damping of such oscillations the system is found in a quasi-stationary state that survives virtually forever when the system is very large. During violent relaxation the distribution function obeys the collisionless Boltzmann (or Vlasov) equation, that, being invariant under time reversal, does not "naturally" describe a relaxation process. Indeed, the dynamics is moved to smaller and smaller scales in phase space as time goes on, so that observables that do not depend on small-scale details appear as relaxed after a short time. We propose an approximation scheme to describe collisionless relaxation, based on the introduction of moments of the distribution function, and apply it to the Hamiltonian Mean Field (HMF) model. To the leading order, virial oscillations are equivalent to the motion of a particle in a one-dimensional potential. Inserting higher-order contributions in an effective way, inspired by the Caldeira-Leggett model of quantum dissipation, we derive a dissipative equation describing the damping of the oscillations, including a renormalization of the effective potential and yielding predictions for collective properties of the system after the damping in very good agreement with numerical simulations. Here we restrict ourselves to "cold" initial conditions; generic initial conditions will be considered in a forthcoming paper.

Observations strongly suggest that filaments in galactic molecular clouds are in a non-thermal state. As a simple model of a filament we study a two-dimensional system of self-gravitating point particles by means of numerical simulations of the dynamics, with various methods: direct $N$-body integration of the equations of motion, particle-in-cell simulations and a recently developed numerical scheme that includes multiparticle collisions in a particle-in-cell approach. Studying the collapse of Gaussian overdensities we find that after the damping of virial oscillations the system settles in a non-thermal steady state whose radial density profile is similar to the observed ones, thus suggesting a dynamical origin of the non-thermal states observed in real filaments. Moreover, for sufficiently cold collapses the density profiles are anticorrelated with the kinetic temperature, i.e., exhibit temperature inversion, again a feature that has been found in some observations of filaments. The same happens in the state reached after a strong perturbation of an initially isothermal cylinder. Finally, we discuss our results in the light of recent findings in other contexts (including non-astrophysical ones) and argue that the same kind of non-thermal states may be observed in any physical system with long-range interactions.

We study the origin of the non-thermal profiles observed in filamentary structures in galactic molecular clouds by means of numerical dynamical simulations. We find that such profiles are intrinsic features of the end products of dissipationless collapse in cylindrical symmetry. Moreover, for sufficiently cold initial conditions, we obtain end states characterized by markedly anticorrelated radial density and temperature profiles. Gravitational, dissipationless dynamics alone is thus sufficient to reproduce, at least qualitatively, many of the properties of the observed non-thermal structures.

Temperature inversions occur in nature, e.g., in the solar corona and in interstellar molecular clouds: somewhat counterintuitively, denser parts of the system are colder than dilute ones. We propose a simple and appealing way to spontaneously generate temperature inversions in systems with long-range interactions, by preparing them in inhomogeneous thermal equilibrium states and then applying an impulsive perturbation. In similar situations, short-range systems would typically relax to another thermal equilibrium, with uniform temperature profile. By contrast, in long-range systems, the interplay between wave-particle interaction and spatial inhomogeneity drives the system to nonequilibrium stationary states that generically exhibit temperature inversion. We demonstrate this mechanism in a simple mean-field model and in a two-dimensional self-gravitating system. Our work underlines the crucial role the range of interparticle interaction plays in determining the nature of steady states out of thermal equilibrium.

Temperature inversion due to velocity filtration, a mechanism originally proposed to explain the heating of the solar corona, is demonstrated to occur also in a simple paradigmatic model with long-range interactions, the Hamiltonian mean-field model. Using molecular dynamics simulations, we show that when the system settles into an inhomogeneous quasi-stationary state in which the velocity distribution has suprathermal tails, the temperature and density profiles are anticorrelated: denser parts of the system are colder than dilute ones. We argue that this may be a generic property of long-range interacting systems.

Self-gravitating systems, like globular clusters or elliptical galaxies, are the prototypes of many-body systems with long-range interactions, and should be the natural arena where to test theoretical predictions on the statistical behaviour of long-range-interacting systems. Systems of classical self-gravitating particles can be studied with the standard tools of equilibrium statistical mechanics, provided the potential is regularized at small length scales and the system is confined in a box. The confinement condition looks rather unphysical in general, so that it is natural to ask whether what we learn with these studies is relevant to real self-gravitating systems. In order to provide a first answer to this question we consider a basic, simple, yet effective model of globular clusters, the King model. This model describes a self-consistently confined system, without the need of any external box, but the stationary state is a non-thermal one. In particular, we consider the King model with a short-distance cutoff on the interactions and we discuss how such a cutoff affects the caloric curve, i.e. the relation between temperature and energy. We find that the cutoff stabilizes a low-energy phase which is absent in the King model without cutoff; the caloric curve of the model with cutoff turns out to be very similar to that of previously studied confined and regularized models, but for the absence of a high-energy gas-like phase. We briefly discuss the possible phenomenological as well as theoretical implications of these results.

We introduce and discuss an effective model of a self-gravitating system whose equilibrium thermodynamics can be solved in both the microcanonical and the canonical ensemble, up to a maximization with respect to a single variable. Such a model can be derived from a model of self-gravitating particles confined on a ring, referred to as the self-gravitating ring (SGR) model, allowing a quantitative comparison between the thermodynamics of the two models. Despite the rather crude approximations involved in its derivation, the effective model compares quite well with the SGR model. Moreover, we discuss the relation between the effective model presented here and another model introduced by Thirring forty years ago. The two models are very similar and can be considered as examples of a class of minimal models of self-gravitating systems.