Tarcisio Nunes Teles

Nonequilibrium statistical mechanics of systems with long-range interactions

Abstract 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 states (qSSs), the lifetime of which diverges with the number of particles. Therefore, in the 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.

Collisionless relaxation in non-neutral plasmas.

A theoretical framework is presented which allows us to quantitatively predict the final stationary state achieved by a non-neutral plasma during a process of collisionless relaxation. As a specific application, the theory is used to study relaxation of charged-particle beams. It is shown that a fully matched beam relaxes to the Lynden-Bell distribution. However, when a mismatch is present and the beam oscillates, parametric resonances lead to a core-halo phase separation. The approach developed accounts for both the density and the velocity distributions in the final stationary state.

Statistical mechanics of unbound two-dimensional self-gravitating systems

We study, using both theory and molecular dynamics simulations, the relaxation dynamics of a microcanonical two-dimensional self-gravitating system. After a sufficiently large time, a gravitational cluster of N particles relaxes to the Maxwell?Boltzmann distribution. The time taken to reach the thermodynamic equilibrium, however, scales with the number of particles. In the thermodynamic limit, at fixed total mass, an equilibrium state is never reached and the system becomes trapped in a non-ergodic stationary state. An analytical theory is presented which allows us to quantitatively describe this final stationary state, without any adjustable parameters.

Statistical mechanics of 1D self-gravitating systems: the core–halo distribution

We study, using both theory and simulations, a system of self-gravitating sheets. A new statistical mechanics theory – free of any adjustable parameters – is derived to quantitatively predict the final stationary state achieved by this system after the process of collisionless relaxation is completed. The theory shows a very good agreement with the numerical simulations. The model sheds new light on the general mechanism of relaxation of self-gravitating systems and may help us to understand cold matter distribution in the Universe.

Ergodicity breaking and parametric resonances in systems with long-range interactions

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.

Nonequilibrium phase transitions in systems with long-range interactions

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.

Relaxation and emittance growth of a thermal charged-particle beam

We present a theory that allows us to accurately calculate the distribution functions and the emittance growth of a thermal charged-particle beam after it relaxes to equilibrium. The theory can be used to obtain the fraction of particles, which will evaporate from the beam to form a halo. The calculated emittance growth is found to be in excellent agreement with the simulations.

Temperature inversion in long-range interacting systems

Temperature inversions occur in nature, e.g., in the solar corona and in interstellar molecular clouds: Somewhat counterintuitively, denser parts of the system are colder than dilute ones. We propose a simple and appealing way to spontaneously generate temperature inversions in systems with long-range interactions, by preparing them in inhomogeneous thermal equilibrium states and then applying an impulsive perturbation. In similar situations, short-range systems would typically relax to another thermal equilibrium, with a uniform temperature profile. By contrast, in long-range systems, the interplay between wave-particle interaction and spatial inhomogeneity drives the system to nonequilibrium stationary states that generically exhibit temperature inversion. We demonstrate this mechanism in a simple mean-field model and in a two-dimensional self-gravitating system. Our work underlines the crucial role the range of interparticle interaction plays in determining the nature of steady states out of thermal equilibrium.

Ensemble inequivalence in a mean-field XY model with ferromagnetic and nematic couplings

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.

Ensemble inequivalence and absence of quasi-stationary states in long-range random networks

Ensemble inequivalence has been previously displayed only for long-range interacting systems with non-extensive energy. In order to perform the thermodynamic limit, such systems require an unphysical, so-called, Kac rescaling of the coupling constant. We here study models defined on long-range random networks, which avoid such a rescaling. The proposed models have an extensive energy, which is however non-additive. For such long-range random networks, pairs of sites are coupled with a probability decaying with the distance