We present the phase diagram, in both the microcanonical and the canonical ensemble, of the Self-Gravitating-Ring (SGR) model, which describes the motion of equal point masses constrained on a ring and subject to 3D gravitational attraction. If the interaction is regularized at short distances by the introduction of a softening parameter, a global entropy maximum always exists, and thermodynamics is well defined in the mean-field limit. However, ensembles are not equivalent and a phase of negative specific heat in the microcanonical ensemble appears in a wide intermediate energy region, if the softening parameter is small enough. The phase transition changes from second to first order at a tricritical point, whose location is not the same in the two ensembles. All these features make of the SGR model the best prototype of a self-gravitating system in one dimension. In order to obtain the stable stationary mass distribution, we apply a new iterative method, inspired by a previous one used in 2D turbulence, which ensures entropy increase and, hence, convergence towards an equilibrium state.

We introduce a one-dimensional toy model of globular clusters. The model is a version of the well-known gravitational sheets system, where we take additionally into account mass and energy loss by evaporation of stars at the boundaries. Numerical integration by the "exact" event-driven dynamics is performed, for initial uniform density and Gaussian random velocities. Two distinct quasi-stationary asymptotic regimes are attained, depending on the initial energy of the system. We guess the forms of the density and velocity profiles which fit numerical data extremely well and allow to perform an independent calculation of the self-consistent gravitational potential. Some power-laws for the asymptotic number of stars and for the collision times are suggested.

We study the link between relaxation to the equilibrium and anomalous superdiffusive motion in a classical N-body hamiltonian system with long-range interaction showing a second-order phase-transition in the canonical ensemble. Anomalous diffusion is observed only in a transient out-of-equilibrium regime and for a small range of energy, below the critical one. Superdiffusion is due to L\'evy walks of single particles and is checked independently through the second moment of the distribution, power spectra, trapping and walking time probabilities. Diffusion becomes normal at equilibrium, after a relaxation time which diverges with N.

In this paper we review a series of results obtained for 1D and 2D simple N-body dynamical models with infinite-range attractive interactions and without short distance singularities. The free energy of both models can be exactly obtained in the canonical ensemble, while microcanonical results can be derived from numerical simulations. Both models show a phase transition from a low energy clustered phase to a high energy gaseous state, in analogy with the models introduced in the early 70's by Thirring and Hertel. The phase transition is second order for the 1D model, first order for the 2D model. Negative specific heat appears in both models near the phase transition point. For both models, in the presence of a negative specific heat, a cluster of collapsed particles coexists with a halo of higher energy particles which perform long correlated flights, which lead to anomalous diffusion. The dynamical origin of the "superdiffusion" is different in the two models, being related to particle trapping and untrapping in the cluster in 1D, while in 2D the channelling of particles in an egg-crate effective potential is responsible of the effect. Both models are Lyapunov unstable and the maximal Lyapunov exponent $\lambda$ has a peak just in the region preceeding the phase transition. Moreover, in the low energy limit $\lambda$ increases proportionally to the square root of the internal energy, while in the high energy region it vanishes as $N^{-1/3}$.

We study the dynamical and statistical behavior of the Hamiltonian Mean Field (HMF) model in order to investigate the relation between microscopic chaos and phase transitions. HMF is a simple toy model of $N$ fully-coupled rotators which shows a second order phase transition. The canonical thermodynamical solution is briefly recalled and its predictions are tested numerically at finite $N$. The Vlasov stationary solution is shown to give the same consistency equation of the canonical solution and its predictions for rotator angle and momenta distribution functions agree very well with numerical simulations. A link is established between the behavior of the maximal Lyapunov exponent and that of thermodynamical fluctuations, expressed by kinetic energy fluctuations or specific heat. The extensivity of chaos in the $N \to \infty$ limit is tested through the scaling properties of Lyapunov spectra and of the Kolmogorov-Sinai entropy. Chaotic dynamics provides the mixing property in phase space necessary for obtaining equilibration; however, the relaxation time to equilibrium grows with $N$, at least near the critical point. Our results constitute an interesting bridge between Hamiltonian chaos in many degrees of freedom systems and equilibrium thermodynamics.