ABSTRACT
We present a Monte Carlo approach that allows us to easily implement Lynden-Bell (LB) entropy maximization for an arbitrary initial particle distribution. The direct maximization of LB entropy for an arbitrary initial distribution requires an infinite number of Lagrange multipliers to account for the Casimir invariants. This has restricted studies of Lynden-Bell's violent relaxation theory to only a very small class of initial conditions of a very simple waterbag form, for which the entropy maximization can be performed numerically. In the present approach, an arbitrary initial distribution is discretized into density levels which are then evolved using an efficient Monte Carlo algorithm towards the final equilibrium state. A comparison is also made between the LB equilibrium and explicit Molecular Dynamics simulations. We find that for most initial distributions, relaxation is incomplete and the system is not able to reach the state of maximum LB entropy. In particular, we see that the tail of the stationary particle distribution is very different from the one predicted by the theory of violent relaxation, with a hard cutoff instead of an algebraic decay predicted by LB's theory.
ABSTRACT
We present evidence that the mechanism proposed in Teles et al. [Phys. Rev. E 92, 020101 (2015)PRESCM1539-375510.1103/PhysRevE.92.020101], referred to as the TGDC mechanism, does apply to a model with repulsive mean-field interactions where it produces temperature inversion in a state whose inhomogeneity is due to an external field. Such evidence contradicts the core statement of the Comment. We also discuss a related issue, concerning the possible application of the TGDC mechanism to the solar corona.
ABSTRACT
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.
ABSTRACT
The classical wave-particle Hamiltonian is considered in its generalized version, where two modes are assumed to interact with the coevolving charged particles. The equilibrium statistical mechanics solution of the model is worked out analytically, both in the canonical and the microcanonical ensembles. The competition between the two modes is shown to yield ensemble inequivalence, at variance with the standard scenario where just one wave is allowed to develop. As a consequence, both temperature jumps and negative specific heat can show up in the microcanonical ensemble. The relevance of these findings for both plasma physics and free electron laser applications is discussed.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
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.
