Caloric curves fitted by polytropic distributions in the HMF model

TL;DR

This paper investigates the quasi-stationary states of the HMF model through numerical simulations, demonstrating that polytropic distributions effectively describe these states and elucidate phenomena like negative specific heat and core-halo structures.

Contribution

It introduces a comprehensive analysis of QSSs in the HMF model using polytropic fits, linking distribution indices to physical caloric curves and complex dynamical behaviors.

Findings

Polytropic distributions fit QSSs in the HMF model well.

Negative kinetic specific heat regions are explained by polytropic indices.

Core-halo structures emerge at low energies in the system.

Abstract

We perform direct numerical simulations of the HMF model starting from non-magnetized initial conditions with a velocity distribution that is (i) gaussian, (ii) semi-elliptical, and (iii) waterbag. Below a critical energy E_c, depending on the initial condition, this distribution is Vlasov dynamically unstable. The system undergoes a process of violent relaxation and quickly reaches a quasi-stationary state (QSS). We find that the distribution function of this QSS can be conveniently fitted by a polytrope with index (i) n=2, (ii) n=1, and (iii) n=1/2. Using the values of these indices, we are able to determine the physical caloric curve T_{kin}(E) and explain the negative kinetic specific heat region C_{kin}=dE/dT_{kin}<0 observed in the numerical simulations. At low energies, we find that the system takes a "core-halo" structure. The core corresponds to the pure polytrope discussed above but it is now surrounded by a halo of particles. We also consider unsteady initial conditions with magnetization M_0=1 and isotropic waterbag distribution and report the complex dynamics of the system creating phase space holes and dense filaments. We show that the kinetic caloric curve is approximately constant, corresponding to a polytrope with index n_0= 3.56. Finally, we consider the collisional evolution of an initially Vlasov stable distribution, and show that the time-evolving distribution function f(v,t) can be fitted by a sequence of polytropic distributions with a time-dependent index n(t) both in the non-magnetized and magnetized regimes. These numerical results show that polytropic distributions (also called Tsallis distributions) provide in many cases a good fit of the QSSs. However, in order to moderate our message, we also report a case where the Lynden-Bell theory provides an excellent prediction of an inhomogeneous QSS.