General linear response formula for non integrable systems obeying the Vlasov equation

Long-range interacting N-particle systems get trapped into long-living out-of-equilibrium stationary states called quasi-stationary states (QSS). We study here the response to a small external perturbation when such systems are settled into a QSS. In the N → ∞ limit the system is described by the Vlasov equation and QSS are mapped into stable stationary solutions of such equation. We consider this problem in the context of a model that has recently attracted considerable attention, the Hamiltonian mean field (HMF) model. For such a model, stationary inhomogeneous and homogeneous states determine an integrable dynamics in the mean-field effective potential and an action-angle transformation allows one to derive an exact linear response formula. However, such a result would be of limited interest if restricted to the integrable case. In this paper, we show how to derive a general linear response formula which does not use integrability as a requirement. The presence of conservation laws (mass, energy, momentum, etc.) and of further Casimir invariants can be imposed a posteriori. We perform an analysis of the infinite time asymptotics of the response formula for a specific observable, the magnetization in the HMF model, as a result of the application of an external magnetic field, for two stationary stable distributions: the Boltzmann-Gibbs equilibrium distribution and the Fermi-Dirac one. When compared with numerical simulations the predictions of the theory are very good away from the transition energy from inhomogeneous to homogeneous states.     

Stationary point analysis of the one-dimensional lattice Landau gauge fixing functional, aka random phase XY Hamiltonian

We study the stationary points of what is known as the lattice Landau gauge fixing functional in one-dimensional compact U(1) lattice gauge theory, or as the Hamiltonian of the one-dimensional random phase XY model in statistical physics. An analytic solution of all stationary points is derived for lattices with an odd number of lattice sites and periodic boundary conditions. In the context of lattice gauge theory, these stationary points and their indices are used to compute the gauge fixing partition function, making reference in particular to the Neuberger problem. Interpreted as stationary points of the one-dimensional XY Hamiltonian, the solutions and their Hessian determinants allow us to evaluate a criterion which makes predictions on the existence of phase transitions and the corresponding critical energies in the thermodynamic limit.     

Dynamics and thermodynamics of rotators interacting with both long- and short-range couplings

The effect of nearest-neighbor coupling on the thermodynamic and dynamical properties of the ferromagnetic Hamiltonian mean field model (HMF) is studied. For a range of antiferromagnetic nearest-neighbor coupling, a canonical first-order transition is observed, and the canonical and microcanonical ensembles are non-equivalent. In studying the relaxation time of non-equilibrium states it is found that as in the HMF model, a class of non-magnetic states is quasi-stationary, with an algebraic divergence of their lifetime with the number of degrees of freedom N. The lifetime of metastable states is found to increase exponentially with N as expected.     

Metastability in the Hamiltonian mean field model and Kuramoto model

We briefly discuss the state of the art on the anomalous dynamics of the Hamiltonian mean field (HMF) model. We stress the important role of the initial conditions for understanding the microscopic nature of the intriguing metastable quasi-stationary states (QSS) observed in the model and the connections to Tsallis statistics and glassy dynamics. We also present new results on the existence of metastable states in the Kuramoto model and discuss the similarities with those found in the HMF model. The existence of metastability seems to be quite a common phenomenon in fully coupled systems, whose origin could be also interpreted as a dynamical mechanism preventing or hindering synchronization.     

Exactly solvable spin ladder model with degenerate ferromagnetic and singlet states

We study the spin ladder model with interactions between spins on neighboring rungs. The model Hamiltonian with the exact singlet ground state degenerated with ferromagnetic state is obtained. The singlet ground state wave function has a special recurrent form and depends on two parameters. Spin correlations in the singlet ground state show double-spiral structure with period of spirals equals to the system size. For special values of parameters they have exponential decay. The spectrum of the model is gapless and there are asymptotically degenerated excited states for special values of parameters in the thermodynamic limit. Received 7 May 1999     

Effective spin-glass Hamiltonian for the anomalous dynamics of the HMF model

We discuss an effective spin-glass Hamiltonian which can be used to study the glassy-like dynamics observed in the metastable states of the Hamiltonian mean field (HMF) model. By means of the Replica formalism, we were able to find a self-consistent equation for the glassy order parameter which reproduces, in a restricted energy region below the phase transition, the microcanonical simulations for the polarization order parameter recently introduced in the HMF model.     

High temperature superfluidity with abnormal fermionic occupancy

We examine a Lipkin based two-level pairing model at finite temperature and in the thermodynamic limit. Whereas at T=0 the model exhibits a superconducting ground state for sufficiently high values of the coupling constant, a partially superconducting phase in whichsome of the particles are paired, is found to survive at high temperatures in a special treatment. This phase is a mixture of “abnormally-occupied” eigenstates, which lie at higher energy, of the interactionless model Hamiltonian.     

Thresholds for epidemic spreading in networks

We study the threshold of epidemic models in quenched networks with degree distribution given by a power-law. For the susceptible-infected-susceptible model the activity threshold λ(c) vanishes in the large size limit on any network whose maximum degree k(max) diverges with the system size, at odds with heterogeneous mean-field (HMF) theory. The vanishing of the threshold has nothing to do with the scale-free nature of the network but stems instead from the largest hub in the system being active for any spreading rate λ>1/√k(max) and playing the role of a self-sustained source that spreads the infection to the rest of the system. The susceptible-infected-removed model displays instead agreement with HMF theory and a finite threshold for scale-rich networks. We conjecture that on quenched scale-rich networks the threshold of generic epidemic models is vanishing or finite depending on the presence or absence of a steady state.     

Self-synchronization of populations of nonlinear oscillators in the thermodynamic limit

A population of identical nonlinear oscillators, subject to random forces and coupled via a mean-field interaction, is studied in the thermodynamic limit. The model presents a nonequilibrium phase transition from a stationary to a time-periodic probability density. Below the transition line, the population of oscillators is in a quiescent state with order parameter equal to zero. Above the transition line, there is a state of collective rhythmicity characterized by a time-periodic behavior of the order parameter and all moments of the probability distribution. The information entropy of the ensemble is a constant both below and above the critical line. Analytical and numerical analyses of the model are provided.     

Dynamical stability of systems with long-range interactions: application of the Nyquist method to the HMF model

We apply the Nyquist method to the Hamiltonian mean field (HMF) model in order to settle the linear dynamical stability of a spatially homogeneous distribution function with respect to the Vlasov equation. We consider the case of Maxwell (isothermal) and Tsallis (polytropic) distributions and show that the system is stable above a critical kinetic temperature T_c and unstable below it. Then, we consider a symmetric double-humped distribution, made of the superposition of two decentered Maxwellians, and show the existence of a re-entrant phase in the stability diagram. When we consider an asymmetric double-humped distribution, the re-entrant phase disappears above a critical value of the asymmetry factor Δ > 1.09. We also consider the HMF model with a repulsive interaction. In that case, single-humped distributions are always stable. For asymmetric double-humped distributions, there is a re-entrant phase for 1 ≤ Δ < 25.6, a double re-entrant phase for 25.6 < Δ < 43.9 and no re-entrant phase for Δ > 43.9. Finally, we extend our results to arbitrary potentials of interaction and mention the connexion between the HMF model, Coulombian plasmas and gravitational systems. We discuss the relation between linear dynamical stability and formal nonlinear dynamical stability and show their equivalence for spatially homogeneous distributions. We also provide a criterion of dynamical stability for spatially inhomogeneous systems.     

Non-Markovian entanglement dynamics of two spin-1/2 particles embedded in two separate spin star baths with tunable external magnetic fields

The exact entanglement dynamics of two spin qubits in two independent spin star baths via a Heisenberg XY interaction in the thermodynamic limit has been investigated by using an operator technique. After the Holstein-Primakoff transformation, the transformed Hamiltonian is effectively equivalent to the Hamiltonian of a Jaynes-Cummings model. The results show that the dynamics of the entanglement exhibits strong non-Markovian behavior and depends on the environmental temperature, the coupling strength between the center spin and the bath, the detuning controlled by a locally applied external magnetic field, as well as the initial state of the two qubits.     

Transport Properties of a Modified Lorentz Gas

We present a detailed study of the first simple mechanical system that shows fully realistic transport behavior while still being exactly solvable at the level of equilibrium statistical mechanics. The system under consideration is a Lorentz gas with fixed freely-rotating circular scatterers interacting with point particles via perfectly rough collisions. Upon imposing a temperature and/or a chemical potential gradient, a stationary state is attained for which local thermal equilibrium holds for low values of the imposed gradients. Transport in this system is normal, in the sense that the transport coefficients which characterize the flow of heat and matter are finite in the thermodynamic limit. Moreover, the two flows are non-trivially coupled, satisfying Onsager's reciprocity relations to within numerical accuracy as well as the Green–Kubo relations. We further show numerically that an applied electric field causes the same currents as the corresponding chemical potential gradient in first order of the applied field. Puzzling discrepancies in higher order effects (Joule heating) are also observed. Finally, the role of entropy production in this purely Hamiltonian system is shortly discussed.     

Dynamics and thermodynamics of a simple model similar to self-gravitating systems: the HMF model

We discuss the dynamics and thermodynamics of the Hamiltonian Mean Field model (HMF) which is a prototypical system with long-range interactions. The HMF model can be seen as the one Fourier component of a one-dimensional self-gravitating system. Interestingly, it exhibits many features of real self-gravitating systems (violent relaxation, persistence of metaequilibrium states, slow collisional dynamics, phase transitions,...) while avoiding complicated problems posed by the singularity of the gravitational potential at short distances and by the absence of a large-scale confinement. We stress the deep analogy between the HMF model and self-gravitating systems by developing a complete parallel between these two systems. This allows us to apply many technics introduced in plasma physics and astrophysics to a new problem and to see how the results depend on the dimension of space and on the form of the potential of interaction. This comparative study brings new light in the statistical mechanics of self-gravitating systems. We also mention simple astrophysical applications of the HMF model in relation with the formation of bars in spiral galaxies.     

No Current Without Heat

We show for a large class of interacting particle systems that whenever the stationary measure is not reversible for the dynamics, then the mean entropy production in the steady state is strictly positive. This extends to the thermodynamic limit the equivalence between microscopic reversibility and zero mean entropy production: time-reversal invariance cannot be spontaneously broken.     

Finite-size effects for phase segregation in a two-dimensional asymmetric exclusion model with two species

We investigate the stationary states of a two-dimensional lattice gas model with exclusion, in the presence of an external field. The lattice is populated by equal numbers of positively and negatively charged particles. An analytical mean-field approach and Monte Carlo simulations give strong evidence of the fact that at any finite density the only relevant stationary state of the system in the thermodynamic limit is inhomogeneous, consisting of a strip of particles transverse to the field. In the inhomogeneous phase, the density profiles and the current measured by Monte Carlo simulations are closely related to those found in mean field. The same is true for the finite-size behavior of the system.     

Liapunov exponents in high-dimensional symplectic dynamics

The existence of a thermodynamic limit of the distribution of Liapunov exponents is numerically verified in a large class of symplectic models, ranging from Hamiltonian flows to maps and products of random matrices. In the highly chaotic regime this distribution is approximately model-independent. Near an integrable limit only a few exponents give a relevant contribution to the Kolmogorov-Sinai entropy.     

Dynamics of a Strongly Coupled Polaron

We study the dynamics of large polarons described by the Fröhlich Hamiltonian in the limit of strong coupling. The initial conditions are (perturbations of) product states of an electron wave function and a phonon coherent state, as suggested by Pekar. We show that, to leading order on the natural time scale of the problem, the phonon field is stationary and the electron moves according to an effective linear Schrödinger equation.     

Quantum chaos triggered by precursors of a quantum phase transition: the dicke model

We consider the Dicke Hamiltonian, a simple quantum-optical model which exhibits a zero-temperature quantum phase transition. We present numerical results demonstrating that at this transition the system changes from being quasi-integrable to quantum chaotic. By deriving an exact solution in the thermodynamic limit we relate this phenomenon to a localization-delocalization transition in which a macroscopic superposition is generated. We also describe the classical analogs of this behavior.