Source author record

Stefano Ruffo

Stefano Ruffo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2023arXiv

Modified Thirring model beyond the excluded-volume approximation

Long-range interacting systems may exhibit ensemble inequivalence and can possibly attain equilibrium states under completely open conditions, for which energy, volume and number of particles simultaneously fluctuate. Here we consider a modified version of the Thirring model for self-gravitating systems with attractive and repulsive long-range interactions in which particles are treated as hard spheres in dimension d=1,2,3. Equilibrium states of the model are studied under completely open conditions, in the unconstrained ensemble, by means of both Monte Carlo simulations and analytical methods and are compared with the corresponding states at fixed number of particles, in the isothermal-isobaric ensemble. Our theoretical description is performed for an arbitrary local equation of state, which allows us to examine the system beyond the excluded-volume approximation. The simulations confirm the theoretical prediction of the possible occurrence of first-order phase transitions in the unconstrained ensemble. This work contributes to the understanding of long-range interacting systems exchanging heat, work and matter with the environment.

preprint2022arXiv

Energy fluctuation relations and repeated quantum measurements

In this review paper, we discuss the statistical description in non-equilibrium regimes of energy fluctuations originated by the interaction between a quantum system and a measurement apparatus applying a sequence of repeated quantum measurements. To properly quantify the information about energy fluctuations, both the exchanged heat probability density function and the corresponding characteristic function are derived and interpreted. Then, we discuss the conditions allowing for the validity of the fluctuation theorem in Jarzynski form $\langle e^{-βQ}\rangle = 1$, thus showing that the fluctuation relation is robust against the presence of randomness in the time intervals between measurements. Moreover, also the late-time, asymptotic properties of the heat characteristic function are analyzed, in the thermodynamic limit of many intermediate quantum measurements. In such a limit, the quantum system tends to the maximally mixed state (thus corresponding to a thermal state with infinite temperature) unless the system's Hamiltonian and the intermediate measurement observable share a common invariant subspace. Then, in this context, we also discuss how energy fluctuation relations change when the system operates in the quantum Zeno regime. Finally, the theoretical results are illustrated for the special cases of two- and three-levels quantum systems, now ubiquitous for quantum applications and technologies.

preprint2022arXiv

Hydrodynamics and transport in the long-range-interacting $φ^4$ chain

We present a simulation study of the one-dimensional $φ^4$ lattice theory with long-range interactions decaying as an inverse power $r^{-(1+σ)}$ of the intersite distance $r$, $σ>0$. We consider the cases of single and double-well local potentials with both attractive and repulsive couplings. The double-well, attractive case displays a phase transition for $0<σ\le 1$ analogous to the Ising model with long-range ferromagnetic interactions. A dynamical scaling analysis of both energy structure factors and excess energy correlations shows that the effective hydrodynamics is diffusive for $σ>1$ and anomalous for $0<σ<1$ where fluctuations propagate superdiffusively. We argue that this is accounted for by a fractional diffusion process and we compare the results with an effective model of energy transport based on Lévy flights. Remarkably, this result is fairly insensitive on the phase transition. Nonequilibrium simulations with an applied thermal gradient are in quantitative agreement with the above scenario.

preprint2022arXiv

Tricritical point in the quantum Hamiltonian mean-field model

Engineering long-range interactions in experimental platforms has been achieved with great success in a large variety of quantum systems in recent years. Inspired by this progress, we propose a generalization of the classical Hamiltonian mean-field model to fermionic particles. We study the phase diagram and thermodynamic properties of the model in the canonical ensemble for ferromagnetic interactions as a function of temperature and hopping. At zero temperature, small charge fluctuations drive the many-body system through a first order quantum phase transition from an ordered to a disordered phase at zero temperature. At higher temperatures, the fluctuation-induced phase transition remains first order initially and switches to second order only at a tricritical point. Our results offer an intriguing example of tricriticality in a quantum system with long-range couplings, which bears direct experimental relevance. The analysis is performed by exact diagonalization and mean-field theory.

preprint2020arXiv

Effective negative specific heat by destabilization of metastable states in dipolar systems

We study dipolarly coupled three dimensional spin systems in both the microcanonical and the canonical ensembles by introducing appropriate numerical methods to determine the microcanonical temperature and by realizing a canonical model of heat bath. In the microcanonical ensemble, we show the existence of a branch of stable antiferromagnetic states in the low energy region. Other metastable ferromagnetic states exist in this region: by externally perturbing them, an effective negative specific heat is obtained. In the canonical ensemble, for low temperatures, the same metastable states are unstable and reach a new branch of more robust metastable states which is distinct from the stable one. Our statistical physics approach allows us to put some order in the complex structure of stable and metastable states of dipolar systems.

preprint2020arXiv

Energetics of critical oscillators in active bacterial baths

We investigate the nonequilibrium energetics near a critical point of a non-linear driven oscillator immersed in an active bacterial bath. At the critical point, we reveal a scaling exponent of the average power $\langle\dot{W}\rangle\sim (D_{\rm a}/τ)^{1/4}$ where $D_{\rm a}$ is the effective diffusivity and $τ$ the correlation time of the bacterial bath described by a Gaussian colored noise. Other features that we investigate are the average stationary power and the variance of the work both below and above the saddle-node bifurcation. Above the bifurcation, the average power attains an optimal, minimum value for finite $τ$ that is below its zero-temperature limit. Furthermore, we reveal a finite-time uncertainty relation for active matter which leads to values of the Fano factor of the work that can be below $2k_{\rm B}T_{\rm eff}$, with $T_{\rm eff}$ the effective temperature of the oscillator in the bacterial bath. We analyze different Markovian approximations to describe the nonequilibrium stationary state of the system. Finally, we illustrate our results in the experimental context by considering the example of driven colloidal particles in periodic optical potentials within an E. Coli bacterial bath.

preprint2020arXiv

Quantum-heat fluctuation relations in $3$-level systems under projective measurements

We study the statistics of energy fluctuations in a three-level quantum system subject to a sequence of projective quantum measurements. We check that, as expected, the quantum Jarzynski equality holds provided that the initial state is thermal. The latter condition is trivially satisfied for two-level systems, while this is generally no longer true for $N$-level systems, with $N > 2$. Focusing on three-level systems, we discuss the occurrence of a unique energy scale factor $β_{\rm eff}$ that formally plays the role of an effective inverse temperature in the Jarzynski equality. To this aim, we introduce a suitable parametrization of the initial state in terms of a thermal and a non-thermal component. We determine the value of $β_{\rm eff}$ for a large number of measurements and study its dependence on the initial state. Our predictions could be checked experimentally in quantum optics.

preprint2019arXiv

Criticality of Spin Systems with Weak Long-Range Interactions

The study of critical properties of systems with long-range interactions has attracted in the last decades a continuing interest and motivated the development of several analytical and numerical techniques, in particular in connection with spin models. From the point of view of the investigation of their criticality, a special role is played by systems in which the interactions are long-range enough that their universality class is different from the short-range case and, nevertheless, they maintain the extensivity of thermodynamical quantities. Such interactions are often called weak long-range. In this paper we focus on the study of the critical behaviour of spin systems with weak-long range couplings using renormalization group, and we review their remarkable properties. For the sake of clarity and self-consistency, we start from the classical $O(N)$ spin models and we then move to quantum spin systems.

preprint2016arXiv

Anisotropic Long-Range Spin Systems

We consider anisotropic long-range interacting spin systems in $d$ dimensions. The interaction between the spins decays with the distance as a power law with different exponents in different directions: we consider an exponent $d_{1}+σ_1$ in $d_1$ directions and another exponent $d_{2}+σ_2$ in the remaining $d_2\equiv d-d_1$ ones. We introduce a low energy effective action with non analytic power of the momenta. As a function of the two exponents $σ_1$ and $σ_2$ we show the system to have three different regimes, two where it is actually anisotropic and one where the isotropy is finally restored. We determine the phase diagram and provide estimates of the critical exponents as a function of the parameters of the system, in particular considering the case of one of the two $σ$'s fixed and the other varying. A discussion of the physical relevance of our results is also presented.

preprint2016arXiv

Out-of-equilibrium fluctuations in stochastic long-range interacting systems

For a many-particle system with long-range interactions and evolving under stochastic dynamics, we study for the first time the out-of-equilibrium fluctuations of the work done on the system by a time-dependent external force. For equilibrium initial conditions, the work distributions for a given protocol of variation of the force in time and the corresponding time-reversed protocol exhibit a remarkable scaling and a symmetry when expressed in terms of the average and the standard deviation of the work. The distributions of the work per particle predict, by virtue of the Crooks fluctuation theorem, the equilibrium free-energy density of the system. For a large number $N$ of particles, the latter is in excellent agreement with the value computed by considering the Langevin dynamics of a single particle in a self-consistent mean field generated by its interaction with other particles. The agreement highlights the effective mean-field nature of the original many-particle dynamics for large $N$. For initial conditions in non-equilibrium steady states (NESSs), we study the distribution of a quantity similar to dissipated work that satisfies the non-equilibrium generalization of the Clausius inequality, namely, the Hatano-Sasa equality, for transitions between NESSs. Besides illustrating the validity of the equality, we show that the distribution has exponential tails that decay differently on the left and on the right.

preprint2016arXiv

Phase transitions in Thirring's model

In his pioneering work on negative specific heat, Walter Thirring in\-tro\-duced a model that is solvable in the microcanonical ensemble. Here, we give a complete description of the phase-diagram of this model in both the microcanonical and the canonical ensemble, highlighting the main features of ensemble inequivalence. In both ensembles, we find a line of first-order phase transitions which ends in a critical point. However, neither the line nor the point have the same location in the phase-diagram of the two ensembles. We also show that the microcanonical and canonical critical points can be analytically related to each other using a Landau expansion of entropy and free energy, respectively, in analogy with what has been done in [O. Cohen, D. Mukamel, J. Stat. Mech., P12017 (2012)]. Examples of systems with certain symmetries restricting the Landau expansion have been considered in this reference, while no such restrictions are present in Thirring's model. This leads to a phase diagram that can be seen as a prototype for what happens in systems of particles with kinematic degrees of freedom dominated by long-range interactions.

preprint2015arXiv

Instabilities and relaxation to equilibrium in long-range oscillator chains

We study instabilities and relaxation to equilibrium in a long-range extension of the Fermi-Pasta-Ulam-Tsingou (FPU) oscillator chain by exciting initially the lowest Fourier mode. Localization in mode space is stronger for the long-range FPU model. This allows us to uncover the sporadic nature of instabilities, i.e., by varying initially the excitation amplitude of the lowest mode, which is the control parameter, instabilities occur in narrow amplitude intervals. Only for sufficiently large values of the amplitude, the system enters a permanently unstable regime. These findings also clarify the long-standing problem of the relaxation to equilibrium in the short-range FPU model. Because of the weaker localization in mode space of this latter model, the transfer of energy is retarded and relaxation occurs on a much longer time-scale.

preprint2015arXiv

Nonequilibrium inhomogeneous steady state distribution in disordered, mean-field rotator systems

We present a novel method to compute the phase space distribution in the nonequilibrium stationary state of a wide class of mean-field systems involving rotators subject to quenched disordered external drive and dissipation. The method involves a series expansion of the stationary distribution in inverse of the damping coefficient; the expansion coefficients satisfy recursion relations whose solution requires computing a sparse matrix, making numerical evaluation simple and efficient. We illustrate our method for the paradigmatic Kuramoto model of spontaneous collective synchronization and for its two mode generalization, in presence of noise and inertia, and demonstrate an excellent agreement between simulations and theory for the phase space distribution.

preprint2015arXiv

Stochastic Quantum Zeno by Large Deviation Theory

Quantum measurements are crucial to observe the properties of a quantum system, which however unavoidably perturb its state and dynamics in an irreversible way. Here we study the dynamics of a quantum system while being subject to a sequence of projective measurements applied at random times. In the case of independent and identically distributed intervals of time between consecutive measurements, we analytically demonstrate that the survival probability of the system to remain in the projected state assumes a large-deviation (exponentially decaying) form in the limit of an infinite number of measurements. This allows us to estimate the typical value of the survival probability, which can therefore be tuned by controlling the probability distribution of the random time intervals. Our analytical results are numerically tested for Zeno-protected entangled states, which also demonstrates that the presence of disorder in the measurement sequence further enhances the survival probability when the Zeno limit is not reached (as it happens in experiments). Our studies provide a new tool for protecting and controlling the amount of quantum coherence in open complex quantum systems by means of tunable stochastic measurements.

preprint2015arXiv

Strong anomalous diffusion of the phase of a chaotic pendulum

In this letter we consider the phase diffusion of a harmonically driven undamped pendulum and show that it is anomalous in the strong sense. The role played by the fractal properties of the phase space is highlighted, providing an illustration of the link between deterministic chaos and anomalous transport. Finally, we build a stochastic model which reproduces most properties of the original Hamiltonian system by alternating ballistic flights and random diffusion.

preprint2015arXiv

Thermodynamics of nonadditive systems

The usual formulation of thermodynamics is based on the additivity of macroscopic systems. However, there are numerous examples of macroscopic systems that are not additive, due to the long-range character of the interaction among the constituents. We present here an approach in which nonadditive systems can be described within a purely thermodynamics formalism. The basic concept is to consider a large ensemble of replicas of the system where the standard formulation of thermodynamics can be naturally applied and the properties of a single system can be consequently inferred. After presenting the approach, we show its implementation in systems where the interaction decays as $1/r^α$ in the interparticle distance $r$, with $α$ smaller than the embedding dimension $d$, and in the Thirring model for gravitational systems.

preprint2014arXiv

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.

preprint2014arXiv

Ensemble Inequivalence in Systems with Wave-Particle Interaction

The classical wave-particle Hamiltonian is considered in its generalized version, where two modes are assumed to interact with the co-evolving 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.

preprint2014arXiv

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 to infinity 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.

preprint2014arXiv

Kuramoto model of synchronization: Equilibrium and nonequilibrium aspects

Recently, there has been considerable interest in the study of spontaneous synchronization, particularly within the framework of the Kuramoto model. The model comprises oscillators with distributed natural frequencies interacting through a mean-field coupling, and serves as a paradigm to study synchronization. In this review, we put forward a general framework in which we discuss in a unified way known results with more recent developments obtained for a generalized Kuramoto model that includes inertial effects and noise. We describe the model from a different perspective, highlighting the long-range nature of the interaction between the oscillators, and emphasizing the equilibrium and out-of-equilibrium aspects of its dynamics from a statistical physics point of view. In the first main section, we introduce the model and discuss for the noiseless and noisy dynamics and unimodal frequency distributions the synchronization transition that occurs in the stationary state. In the second section, we introduce the generalized model, and discuss its synchronization phase diagram for unimodal frequency distributions. In the third section, we describe deviations from the mean-field setting of the Kuramoto model by considering the generalized dynamics on a one-dimensional periodic lattice on the sites of which the oscillators reside and interact with a coupling that decays as an inverse power-law of their separation. For two specific cases, namely, in the absence of noise and inertia, and in the case when the natural frequencies are the same for all the oscillators, we discuss how the long-time transition to synchrony is governed by the dynamics of the mean-field mode (zero Fourier mode) of the spatial distribution of the oscillator phases.

preprint2014arXiv

Nonequilibrium first-order transition in coupled oscillator systems with inertia and noise

We study the dynamics of a system of coupled oscillators of distributed natural frequencies, by including the features of both thermal noise, parametrized by a temperature, and inertial terms, parametrized by a moment of inertia. For a general unimodal frequency distribution, we report here the complete phase diagram of the model in the space of dimensionless moment of inertia, temperature, and width of the frequency distribution. We demonstrate that the system undergoes a nonequilibrium first-order phase transition from a synchronized phase at low parameter values to an incoherent phase at high values. We provide strong numerical evidence for the existence of both the synchronized and the incoherent phase, treating the latter analytically to obtain the corresponding linear stability threshold that bounds the first-order transition point from below. In the limit of zero noise and inertia, when the dynamics reduces to the one of the Kuramoto model, we recover the associated known continuous transition. At finite noise and inertia but in the absence of natural frequencies, the dynamics becomes that of a well-studied model of long-range interactions, the Hamiltonian mean-field model. Close to the first-order phase transition, we show that the escape time out of metastable states scales exponentially with the number of oscillators, which we explain to be stemming from the long-range nature of the interaction between the oscillators.

preprint2013arXiv

A stochastic model of long-range interacting particles

We introduce a model of long-range interacting particles evolving under a stochastic Monte Carlo dynamics, in which possible increase or decrease in the values of the dynamical variables is accepted with preassigned probabilities. For symmetric increments, the system at long times settles to the Gibbs equilibrium state, while for asymmetric updates, the steady state is out of equilibrium. For the associated Fokker-Planck dynamics in the thermodynamic limit, we compute exactly the phase space distribution in the nonequilibrium steady state, and find that it has a nontrivial form that reduces to the familiar Gibbsian measure in the equilibrium limit.

preprint2012arXiv

Gravitation in the Microcanonical Ensemble: Appropriate Scaling Leading to Extensivity and Thermalization

We introduce a simple model of hard spheres with gravitational interactions, for which we study a suitable scaling limit. Usual extensive properties are maintained notwithstanding the long range of gravitational interaction. We show that a local thermalization spontaneously emerges within a microcanonical description of the stationary state. In the considered scaling limit, the resulting density profile can be determined in a hydrostatic approach.

preprint2012arXiv

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.

preprint2012arXiv

Kinetic theory of nonequilibrium stochastic long-range systems: Phase transition and bistability

We study long-range interacting systems driven by external stochastic forces that act collectively on all the particles constituting the system. Such a scenario is frequently encountered in the context of plasmas, self-gravitating systems, two-dimensional turbulence, and also in a broad class of other systems. Under the effect of stochastic driving, the system reaches a stationary state where external forces balance dissipation on average. These states have the invariant probability that does not respect detailed balance, and are characterized by non-vanishing currents of conserved quantities. In order to analyze spatially homogeneous stationary states, we develop a kinetic approach that generalizes the one known for deterministic long-range systems; we obtain a very good agreement between predictions from kinetic theory and extensive numerical simulations. Our approach may also be generalized to describe spatially inhomogeneous stationary states. We also report on numerical simulations exhibiting a first-order nonequilibrium phase transition from homogeneous to inhomogeneous states. Close to the phase transition, the system shows bistable behavior between the two states, with a mean residence time that diverges as an exponential in the inverse of the strength of the external stochastic forces, in the limit of low values of such forces.

preprint2012arXiv

Linear response theory for long-range interacting systems in quasistationary states

Long-range interacting systems, while relaxing to equilibrium, often get trapped in long-lived quasistationary states which have lifetimes that diverge with the system size. In this work, we address the question of how a long-range system in a quasistationary state (QSS) responds to an external perturbation. We consider a long-range system that evolves under deterministic Hamilton dynamics. The perturbation is taken to couple to the canonical coordinates of the individual constituents. Our study is based on analyzing the Vlasov equation for the single-particle phase space distribution. The QSS represents stable stationary solution of the Vlasov equation in the absence of the external perturbation. In the presence of small perturbation, we linearize the perturbed Vlasov equation about the QSS to obtain a formal expression for the response observed in a single-particle dynamical quantity. For a QSS that is homogeneous in the coordinate, we obtain an explicit formula for the response. We apply our analysis to a paradigmatic model, the Hamiltonian mean-field model, that involves particles moving on a circle under Hamilton dynamics. Our prediction for the response of three representative QSSs in this model (the water-bag QSS, the Fermi-Dirac QSS, and the Gaussian QSS) is found to be in good agreement with $N$-particle simulations for large $N$. We also show the long-time relaxation of the water-bag QSS to the Boltzmann-Gibbs equilibrium state.

preprint2012arXiv

One-dimensional lattice of oscillators coupled through power-law interactions: Continuum limit and dynamics of spatial Fourier modes

We study synchronization in a system of phase-only oscillators residing on the sites of a one-dimensional periodic lattice. The oscillators interact with a strength that decays as a power law of the separation along the lattice length and is normalized by a size-dependent constant. The exponent $α$ of the power law is taken in the range $0 \le α<1$. The oscillator frequency distribution is symmetric about its mean (taken to be zero), and is non-increasing on $[0,\infty)$. In the continuum limit, the local density of oscillators evolves in time following the continuity equation that expresses the conservation of the number of oscillators of each frequency under the dynamics. This equation admits as a stationary solution the unsynchronized state uniform both in phase and over the space of the lattice. We perform a linear stability analysis of this state to show that when it is unstable, different spatial Fourier modes of fluctuations have different stability thresholds beyond which they grow exponentially in time with rates that depend on the Fourier modes. However, numerical simulations show that at long times, all the non-zero Fourier modes decay in time, while only the zero Fourier mode (i.e., the "mean-field" mode) grows in time, thereby dominating the instability process and driving the system to a synchronized state. Our theoretical analysis is supported by extensive numerical simulations.

preprint2012arXiv

Overdamped dynamics of long-range systems on a one-dimensional lattice: Dominance of the mean-field mode and phase transition

We consider the overdamped dynamics of a paradigmatic long-range system of particles residing on the sites of a one-dimensional lattice, in the presence of thermal noise. The internal degree of freedom of each particle is a periodic variable which is coupled to those of other particles with an attractive XY-like interaction. The coupling strength decays with the interparticle separation $r$ in space as $1/r^α$; ~$0 < α< 1$. We study the dynamics of the model in the continuum limit by considering the Fokker-Planck equation for the evolution of the spatial density of particles. We show that the equation allows a linearly stable stationary state which is always uniform in space, being non-uniform in the internal degrees below a critical temperature $T=1/2$ and uniform above, with a phase transition between the two at $T=1/2$. The state is the same as the equilibrium state of the mean-field version of the model, obtained by considering $α=0$. Our analysis also lets us to compute the growth and decay rates of spatial Fourier modes of density fluctuations. The growth rates compare very well with numerical simulations.

preprint2011arXiv

Phase transitions of quasistationary states in the Hamiltonian Mean Field model

The out-of-equilibrium dynamics of the Hamiltonian Mean Field (HMF) model is studied in presence of an externally imposed magnetic field h. Lynden-Bell's theory of violent relaxation is revisited and shown to adequately capture the system dynamics, as revealed by direct Vlasov based numerical simulations in the limit of vanishing field. This includes the existence of an out-of-equilibrium phase transition separating magnetized and non magnetized phases. We also monitor the fluctuations in time of the magnetization, which allows us to elaborate on the choice of the correct order parameter when challenging the performance of Lynden-Bell's theory. The presence of the field h removes the phase transition, as it happens at equilibrium. Moreover, regions with negative susceptibility are numerically found to occur, in agreement with the predictions of the theory.

preprint2011arXiv

Self-consistent inhomogeneous steady states in Hamiltonian mean field dynamics

Long-lived quasistationary states, associated with stationary stable solutions of the Vlasov equation, are found in systems with long-range interactions. Studies of the relaxation time in a model of $N$ globally coupled particles moving on a ring, the Hamiltonian Mean Field model (HMF), have shown that it diverges as $N^γ$ for large $N$, with $γ\simeq 1.7$ for some initial conditions with homogeneously distributed particles. We propose a method for identifying exact inhomogeneous steady states in the thermodynamic limit, based on analysing models of uncoupled particles moving in an external field. For the HMF model, we show numerically that the relaxation time of these states diverges with $N$ with the exponent $γ\simeq 1$. The method, applicable to other models with globally coupled particles, also allows an exact evaluation of the stability limit of homogeneous steady states. In some cases it provides a good approximation for the correspondence between the initial condition and the final steady state.

preprint2011arXiv

Vlasov equation for long-range interactions on a lattice

We show that, in the continuum limit, the dynamics of Hamiltonian systems defined on a lattice with long-range couplings is well described by the Vlasov equation. This equation can be linearized around the homogeneous state and a dispersion relation, that depends explicitly on the Fourier modes of the lattice, can be derived. This allows one to compute the stability thresholds of the homogeneous state, which turn out to depend on the mode number. When this state is unstable, the growth rates are also function of the mode number. Explicit calculations are performed for the $α$-HMF model with $0 \leq α<1$, for which the zero mean-field mode is always found to dominate the exponential growth. The theoretical predictions are successfully compared with numerical simulations performed on a finite lattice.

preprint2010arXiv

Models with short and long-range interactions: phase diagram and reentrant phase

We study the phase diagram of two different Hamiltonians with competiting local, nearest-neighbour, and mean-field couplings. The first example corresponds to the HMF Hamiltonian with an additional short-range interaction. The second example is a reduced Hamiltonian for dipolar layered spin structures, with a new feature with respect to the first example, the presence of anisotropies. The two examples are solved in both the canonical and the microcanonical ensemble using a combination of the min-max method with the transfer operator method. The phase diagrams present typical features of systems with long-range interactions: ensemble inequivalence, negative specific heat and temperature jumps. Moreover, in a given range of parameters, we report the signature of phase reentrance. This can also be interpreted as the presence of azeotropy with the creation of two first order phase transitions with ensemble inequivalence, as one parameter is varied continuously.

preprint2010arXiv

Scaling with system size of the Lyapunov exponents for the Hamiltonian Mean Field model

The Hamiltonian Mean Field (HMF) model is a prototype for systems with long-range interactions. It describes the motion of $N$ particles moving on a ring, coupled through an infinite-range potential. The model has a second order phase transition at the energy $U_c=3/4$ and its dynamics is exactly described by the Vlasov equation in the $N \to \infty$ limit. Its chaotic properties have been investigated in the past, but the determination of the scaling with $N$ of the Lyapunov Spectrum (LS) of the model remains a challenging open problem. We here show that the $N^{-1/3}$ scaling of the Maximal Lyapunov Exponent (MLE), found in previous numerical and analytical studies, extends to the full LS; not only, scaling is "precocious" for the LS, meaning that it becomes manifest for a much smaller number of particles than the one needed to check the scaling for the MLE. Besides that, the $N^{-1/3}$ scaling appears to be valid not only for $U>U_c$, as suggested by theoretical approaches based on a random matrix approximation, but also below a threshold energy $U_t \approx 0.2$. Using a recently proposed method (GALI) devised to rapidly check the chaotic or regular nature of an orbit, we find that $U_t$ is also the energy at which a sharp transition from {\it weak} to {\it strong} chaos is present in the phase-space of the model. Around this energy the phase of the vector order parameter of the model becomes strongly time dependent, inducing a significant untrapping of particles from a nonlinear resonance.

preprint2010arXiv

Statistical mechanics of collisionless relaxation in a noninteracting system

We introduce a model of uncoupled pendula, which mimics the dynamical behavior of the Hamiltonian Mean Field (HMF) model. This model has become a paradigm for long-range interactions, like Coulomb or dipolar forces. As in the HMF model, this simplified integrable model is found to obey the Vlasov equation and to exhibit Quasi Stationary States (QSS), which arise after a "collisionless" relaxation process. Both the magnetization and the single particle distribution function in these QSS can be predicted using Lynden-Bell's theory. The existence of an extra conserved quantity for this model, the energy distribution function, allows us to understand the origin of some discrepancies of the theory with numerical experiments. It also suggests us an improvement of Lynden-Bell's theory, which we fully implement for the zero field case.

preprint2009arXiv

Finite Size Effects for the Ising Model on Random Graphs with Varying Dilution

We investigate the finite size corrections to the equilibrium magnetization of an Ising model on a random graph with $N$ nodes and $N^γ$ edges, with $1 < γ\leq 2$. By conveniently rescaling the coupling constant, the free energy is made extensive. As expected, the system displays a phase transition of the mean-field type for all the considered values of $γ$ at the transition temperature of the fully connected Curie-Weiss model. Finite size corrections are investigated for different values of the parameter $γ$, using two different approaches: a replica-based finite $N$ expansion, and a cavity method. Numerical simulations are compared with theoretical predictions. The cavity based analysis is shown to agree better with numerics.

preprint2007arXiv

Exploring the thermodynamic limit of Hamiltonian models: convergence to the Vlasov equation

We here discuss the emergence of Quasi Stationary States (QSS), a universal feature of systems with long-range interactions. With reference to the Hamiltonian Mean Field (HMF) model, numerical simulations are performed based on both the original $N$-body setting and the continuum Vlasov model which is supposed to hold in the thermodynamic limit. A detailed comparison unambiguously demonstrates that the Vlasov-wave system provides the correct framework to address the study of QSS. Further, analytical calculations based on Lynden-Bell's theory of violent relaxation are shown to result in accurate predictions. Finally, in specific regions of parameters space, Vlasov numerical solutions are shown to be affected by small scale fluctuations, a finding that points to the need for novel schemes able to account for particles correlations.

preprint2007arXiv

Out of equilibrium phase transitions in mean field Hamiltonian dynamics

Systems with long-range interactions display a short-time relaxation towards Quasi-Stationary States (QSSs), whose lifetime increases with system size. With reference to the Hamiltonian Mean Field (HMF) model, we here review Lynden-Bell's theory of ``violent relaxation''. The latter results in a maximum entropy scheme for a water-bag initial profile which predicts the presence of out-of-equilibrium phase transitions} separating homogeneous (zero magnetization) from inhomogeneous (non-zero magnetization) QSSs. Two different parametric representations of the initial condition are analyzed and the features of the phase diagram are discussed. In both representations we find a second order and a first order line of phase transitions that merge at a tricritical point. Particular attention is payed to the condition of existence and stability of the homogenous phase.

preprint2007arXiv

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.

preprint2004arXiv

Nonlinear supratransmission and bistability in the Fermi-Pasta-Ulam model

The recently discovered phenomenon of nonlinear supratransmission consists in a sudden increase of the amplitude of a transmitted wave triggered by the excitation of nonlinear localized modes of the medium. We examine this process for the Fermi-Pasta-Ulam chain, sinusoidally driven at one edge and damped at the other. The supratransmission regime occurs for driving frequencies above the upper band-edge and originates from direct moving discrete breather creation. We derive approximate analytical estimates of the supratransmission threshold, which are in excellent agreement with numerics. When analysing the long-time behavior, we discover that, below the supratransmission threshold, a conducting stationary state coexists with the insulating one. We explain the bistable nature of the energy flux in terms of the excitation of quasi-harmonic extended waves. This leads to the analytical calculation of a lower-transmission threshold which is also in reasonable agreement with numerical experiments.

preprint2000arXiv

Violation of ensemble equivalence in the antiferromagnetic mean-field XY model

It is well known that long-range interactions pose serious problems for the formulation of statistical mechanics. We show in this paper that ensemble equivalence is violated in a simple mean-field model of N fully coupled classical rotators with repulsive interaction (antiferromagnetic XY model). While in the canonical ensemble the rotators are randomly dispersed over all angles, in the microcanonical ensemble a bi-cluster of rotators separated by angle $π$, forms in the low energy limit. We attribute this behavior to the extreme degeneracy of the ground state: only one harmonic mode is present, together with N-1 zero modes. We obtain empirically an analytical formula for the probability density function for the angle made by the rotator, which compares extremely well with numerical data and should become exact in the zero energy limit. At low energy, in the presence of the bi-cluster, an extensive amount of energy is located in the single harmonic mode, with the result that the energy temperature relation is modified. Although still linear, $T = αU$, it has the slope $α\approx 1.3$, instead of the canonical value $α=2$.

preprint1998arXiv

Analytical Estimation of the Maximal lyapunov Exponent in Oscillator Chains

An analytical expression for the maximal Lyapunov exponent $λ_1$ in generalized Fermi-Pasta-Ulam oscillator chains is obtained. The derivation is based on the calculation of modulational instability growth rates for some unstable periodic orbits. The result is compared with numerical simulations and the agreement is good over a wide range of energy densities $ε$. At very high energy density the power law scaling of $λ_1$ with $ε$ can be also obtained by simple dimensional arguments, assuming that the system is ruled by a single time scale. Finally, we argue that for repulsive and hard core potentials in one dimension $λ_1 \sim \sqrtε$ at large $ε$.

preprint1998arXiv

Chaos and Statistical Mechanics in the Hamiltonian Mean Field Model

We study the dynamical and statistical behavior of the Hamiltonian Mean Field (HMF) model in order to investigate the relation between microscopic chaos and phase transitions. HMF is a simple toy model of $N$ fully-coupled rotators which shows a second order phase transition. The canonical thermodynamical solution is briefly recalled and its predictions are tested numerically at finite $N$. The Vlasov stationary solution is shown to give the same consistency equation of the canonical solution and its predictions for rotator angle and momenta distribution functions agree very well with numerical simulations. A link is established between the behavior of the maximal Lyapunov exponent and that of thermodynamical fluctuations, expressed by kinetic energy fluctuations or specific heat. The extensivity of chaos in the $N \to \infty$ limit is tested through the scaling properties of Lyapunov spectra and of the Kolmogorov-Sinai entropy. Chaotic dynamics provides the mixing property in phase space necessary for obtaining equilibration; however, the relaxation time to equilibrium grows with $N$, at least near the critical point. Our results constitute an interesting bridge between Hamiltonian chaos in many degrees of freedom systems and equilibrium thermodynamics.

preprint1998arXiv

Hamiltonian Dynamics and the Phase Transition of the XY Model

A Hamiltonian dynamics is defined for the XY model by adding a kinetic energy term. Thermodynamical properties (total energy, magnetization, vorticity) derived from microcanonical simulations of this model are found to be in agreement with canonical Monte-Carlo results in the explored temperature region. The behavior of the magnetization and the energy as functions of the temperature are thoroughly investigated, taking into account finite size effects. By representing the spin field as a superposition of random phased waves, we derive a nonlinear dispersion relation whose solutions allow the computation of thermodynamical quantities, which agree quantitatively with those obtained in numerical experiments, up to temperatures close to the transition. At low temperatures the propagation of phonons is the dominant phenomenon, while above the phase transition the system splits into ordered domains separated by interfaces populated by topological defects. In the high temperature phase, spins rotate, and an analogy with an Ising-like system can be established, leading to a theoretical prediction of the critical temperature $T_{KT}\approx 0.855$.

preprint1998arXiv

Modular Transformations, Order-Chaos Transitions and Pseudo-Random Number Generation

Successive pairs of pseudo-random numbers generated by standard linear congruential transformations display ordered patterns of parallel lines. We study the ``ordered'' and ``chaotic'' distribution of such pairs by solving the eigenvalue problem for two-dimensional modular transformations over integers. We conjecture that the optimal uniformity for pair distribution is obtained when the slope of linear modular eigenspaces takes the value $n_{opt} = maxint(p /\sqrt{p-1})$, where $p$ is a prime number. We then propose a new generator of pairs of independent pseudo-random numbers, which realizes an optimal uniform distribution (in the ``statistical'' sense) of points on the unit square $(0,1] \times (0,1]$. The method can be easily generalized to the generation of $k$-tuples of random numbers (with $k>2$)

preprint1997arXiv

Localization and Equipartition of Energy in the beta-FPU Chain : Chaotic Breathers

The evolution towards equipartition in the $β$-FPU chain is studied considering as initial condition the highest frequency mode. Above an analytically derived energy threshold, this zone-boundary mode is shown to be modulationally unstable and to give rise to a striking localization process. The spontaneously created excitations have strong similarity with moving exact breathers solutions. But they have a finite lifetime and their dynamics is chaotic. These chaotic breathers are able to collect very efficiently the energy in the chain. Therefore their size grows in time and they can transport a very large quantity of energy. These features can be explained analyzing the dynamics of perturbed exact breathers of the FPU chain. In particular, a close connection between the Lyapunov spectrum of the chaotic breathers and the Floquet spectrum of the exact ones has been found. The emergence of chaotic breathers is convincingly explained by the absorption of high frequency phonons whereas a breather's metastability is for the first time identified. The lifetime of the chaotic breather is related to the time necessary for the system to reach equipartition. The equipartition time turns out to be dependent on the system energy density $ε$ only. Moreover, such time diverges as $ε^{-2}$ in the limit $ε\to 0$ and vanishes as $ε^{-1/4}$ for $ε\to \infty$.

preprint1997arXiv

Transition between immune and disease states in a cellular automaton model of clonal immune response

In this paper we extend the Celada-Seiden (CS) model of the humoral immune response to include infectious virus and cytotoxic T lymphocytes (cellular response). The response of the system to virus involves a competition between the ability of the virus to kill the host cells and the host's ability to eliminate the virus. We find two basins of attraction in the dynamics of this system, one is identified with disease and the other with the immune state. There is also an oscillating state that exists on the border of these two stable states. Fluctuations in the population of virus or antibody can end the oscillation and drive the system into one of the stable states. The introduction of mechanisms of cross-regulation between the two responses can bias the system towards one of them. We also study a mean field model, based on coupled maps, to investigate virus-like infections. This simple model reproduces the attractors for average populations observed in the cellular automaton. All the dynamical behavior connected to spatial extension is lost, as is the oscillating feature. Thus the mean field approximation introduced with coupled maps destroys oscillations.