Systems with long-range (LR) forces, for which the interaction potential decays with the interparticle distance with an exponent smaller than the dimensionality of the embedding space, remain an outstanding challenge to statistical physics. The internal energy of such systems lacks extensivity and additivity. Although the extensivity can be restored by scaling the interaction potential with the number of particles, the non-additivity still remains. Lack of additivity leads to inequivalence of statistical ensembles. Before relaxing to thermodynamic equilibrium, isolated systems with LR forces become trapped in out-of-equilibrium quasi-stationary state (qSS), the lifetime of which diverges with the number of particles. Therefore, in thermodynamic limit LR systems will not relax to equilibrium. The qSSs are attained through the process of collisionless relaxation. Density oscillations lead to particle-wave interactions and excitation of parametric resonances. The resonant particles escape from the main cluster to form a tenuous halo. Simultaneously, this cools down the core of the distribution and dampens out the oscillations. When all the oscillations die out the ergodicity is broken and a qSS is born. In this report, we will review a theory which allows us to quantitatively predict the particle distribution in the qSS. The theory is applied to various LR interacting systems, ranging from plasmas to self-gravitating clusters and kinetic spin models.

We propose a class of mean-field models for the isostatic transition of systems of soft spheres, in which the contact network is modeled as a random graph and each contact is associated to d degrees of freedom. We study such models in the hypostatic, isostatic, and hyperstatic regimes. The density of states is evaluated by both the cavity method and exact diagonalization of the dynamical matrix. We show that the model correctly reproduces the main features of the density of states of real packings and, moreover, it predicts the presence of localized modes near the lower band edge. Finally, the behavior of the density of states D(ω)∼ω^{α} for ω→0 in the hyperstatic regime is studied. We find that the model predicts a nontrivial dependence of α on the details of the coordination distribution.

Three-dimensional self-gravitating systems do not evolve to thermodynamic equilibrium but become trapped in nonequilibrium quasistationary states. In this Letter, we present a theory which allows us to a priori predict the particle distribution in a final quasistationary state to which a self-gravitating system will evolve from an initial condition which is isotropic in particle velocities and satisfies a virial constraint 2K=-U, where K is the total kinetic energy, and U is the potential energy of the system.

Systems of particles with long-range interactions present two important processes: first, the formation of out-of-equilibrium quasistationary states (QSS) and, second, the collisional relaxation towards Maxwell-Boltzmann equilibrium in a much longer time scale. In this paper, we study the collisional relaxation in the Hamiltonian mean-field model using the appropriate kinetic equations for a system of N particles at order 1/N: the Landau equation when collective effects are neglected and the Lenard-Balescu equation when they are taken into account. We derive explicit expressions for the diffusion coefficients using both equations for any magnetization, and we obtain analytic expressions for highly clustered configurations. An important conclusion is that in this system collective effects are crucial in order to describe the relaxation dynamics. We compare the diffusion calculated with the kinetic equations with simulations set up to simulate the system with or without collective effects, obtaining a very good agreement between theory and simulations.