Controversy about the applicability of Tsallis statistics to the HMF model

Comment to "Nonextensive Thermodynamics and Glassy Behaviour in Hamiltonian Systems" by A. Rapisarda and A. Pluchino, Europhysics News 36, 202 (2005)

Phase space gaps and ergodicity breaking in systems with long range interactions

We study a generalized isotropic XY-model which includes both two-spin and four-spin mean-field interactions. This model can be solved in the microcanonical ensemble. It is shown that in certain parameter regions the model exhibits gaps in the magnetization at fixed energy, resulting in ergodicity breaking. This phenomenon has previously been reported in anisotropic and discrete spin models. The entropy of the model is calculated and the microcanonical phase diagram is derived, showing the existence of first order phase transitions from the ferromagnetic to a paramagnetic disordered phase. It is found that ergodicity breaking takes place both in the ferromagnetic and the paramagnetic phases. As a consequence, the system can exhibit a stable ferromagnetic phase within the paramagnetic region, and conversely a disordered phase within the magnetically ordered region

The Vlasov equation and the Hamiltonian Mean-Field model

We show that the quasi-stationary states observed in the $N$-particle dynamics of the Hamiltonian Mean-Field (HMF) model are nothing but Vlasov stable homogeneous (zero magnetization) states. There is an infinity of Vlasov stable homogeneous states corresponding to different initial momentum distributions. Tsallis $q$-exponentials in momentum, homogeneous in angle, distribution functions are possible, however, they are not special in any respect, among an infinity of others. All Vlasov stable homogeneous states lose their stability because of finite $N$ effects and, after a relaxation time diverging with a power-law of the number of particles, the system converges to the Boltzmann-Gibbs equilibrium

Prediction of anomalous diffusion and algebraic relaxations for long-range interacting systems, using classical statistical mechanics

We explain the ubiquity and extremely slow evolution of non gaussian out-of-equilibrium distributions for the Hamiltonian Mean-Field model, by means of traditional kinetic theory. Deriving the Fokker-Planck equation for a test particle, one also unambiguously explains and predicts striking slow algebraic relaxation of the momenta autocorrelation, previously found in numerical simulations. Finally, angular anomalous diffusion are predicted for a large class of initial distributions. Non Extensive Statistical Mechanics is shown to be unnecessary for the interpretation of these phenomena

Algebraic Correlation Function and Anomalous Diffusion in the HMF model

In the quasi-stationary states of the Hamiltonian Mean-Field model, we numerically compute correlation functions of momenta and diffusion of angles with homogeneous initial conditions. This is an example, in a N-body Hamiltonian system, of anomalous transport properties characterized by non exponential relaxations and long-range temporal correlations. Kinetic theory predicts a striking transition between weak anomalous diffusion and strong anomalous diffusion. The numerical results are in excellent agreement with the quantitative predictions of the anomalous transport exponents. Noteworthy, also at statistical equilibrium, the system exhibits long-range temporal correlations: the correlation function is inversely proportional to time with a logarithmic correction instead of the usually expected exponential decay, leading to weak anomalous transport properties

Kinetic theory for non-equilibrium stationary states in long-range interacting systems

We study long-range interacting systems perturbed by external stochastic forces. Unlike the case of short-range systems, where stochastic forces usually act locally on each particle, here we consider perturbations by external stochastic fields. The system reaches stationary states where external forces balance dissipation on average. These states do not respect detailed balance and support non-vanishing fluxes of conserved quantities. We generalize the kinetic theory of isolated long-range systems to describe the dynamics of this non-equilibrium problem. The kinetic equation that we obtain applies to plasmas, self-gravitating systems, and to a broad class of other systems. Our theoretical results hold for homogeneous states, but may also be generalized to apply to inhomogeneous states. We obtain an excellent agreement between our theoretical predictions and numerical simulations. We discuss possible applications to describe non-equilibrium phase transitions.Comment: 11 pages, 2 figures; v2: small changes, close to the published versio

Thermodynamics of the self-gravitating ring model

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