We explore the mechanism responsible for the ergodicity breaking in systems with long-range forces. In thermodynamic limit such systems do not evolve to the Boltzmann-Gibbs equilibrium, but become trapped in an out-of-equilibrium quasi-stationary-state. Nevertheless, we show that if the initial distribution satisfies a specific constraint-a generalized virial condition-the quasistationary state is very close to ergodic and can be described by Lynden-Bell statistics. On the other hand, if the generalized virial condition is violated, parametric resonances are excited, leading to chaos and ergodicity breaking.

We introduce a generalized Hamiltonian mean field model-an XY model with both linear and quadratic coupling between spins and explicit Hamiltonian dynamics. In addition to the usual paramagnetic and ferromagnetic phases, this model also possesses a nematic phase. The generalized Hamiltonian mean field model can be solved explicitly using Boltzmann-Gibbs statistical mechanics, in both canonical and microcanonical ensembles. However, when the resulting microcanonical phase diagram is compared with the one obtained using molecular dynamics simulations, it is found that the two are very different. We will present a dynamical theory which allows us to explicitly calculate the phase diagram obtained using molecular dynamics simulations without any adjustable parameters. The model illustrates the fundamental role played by dynamics as well the inadequacy of Boltzmann-Gibbs statistics for systems with long-range forces in the thermodynamic limit.

We explore ensemble inequivalence in long-range interacting systems by studying an XY model of classical spins with ferromagnetic and nematic coupling. We demonstrate the inequivalence by mapping the microcanonical phase diagram onto the canonical one, and also by doing the inverse mapping. We show that the equilibrium phase diagrams within the two ensembles strongly disagree within the regions of first-order transitions, exhibiting interesting features like temperature jumps. In particular, we discuss the coexistence and forbidden regions of different macroscopic states in both the phase diagrams.

In the thermodynamic limit, systems with long-range interactions do not relax to equilibrium, but become trapped in nonequilibrium stationary states. For a finite number of particles a nonequilibrium state has a finite lifetime, so that eventually a system will relax to thermodynamic equilibrium. The time that a system remains trapped in a quasistationary state (QSS) scales with the number of particles as N δ , with δ > 0, and diverges in the thermodynamic limit. In this paper we will explore the role of chaotic dynamics on the time that a system remains trapped in a QSS. We discover that chaos, measured by the Lyapunov exponents, favors faster relaxation to equilibrium. Surprisingly, weak chaos favors faster relaxation than strong chaos.