Microcanonical solution of the mean-field ϕ4\phi^4 model: comparison with time averages at finite size

We solve the mean-field $\phi^4$ model in an external magnetic field in the microcanonical ensemble using two different methods. The first one is based on Rugh's microcanonical formalism and leads to express macroscopic observables, such as temperature, specific heat, magnetization and susceptibility, as time averages of convenient functions of the phase-space. The approach is applicable for any finite number of particles $N$. The second method uses large deviation techniques and allows us to derive explicit expressions for microcanonical entropy and for macroscopic observables in the $N \to\infty$ limit. Assuming ergodicity, we evaluate time averages in molecular dynamics simulations and, using Rugh's approach, we determine the value of macroscopic observables at finite $N$. These averages are affected by a slow time evolution, often observed in systems with long-range interactions. We then show how the finite $N$ time averages of macroscopic observables converge to their corresponding $N\to\infty$ values as $N$ is increased. As expected, finite size effects scale as $N^{-1}$.Comment: 18 pages, 1 figur

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.Comment: 18 pages, 7 figures, (revised version, several minor typos fixed - accepted for publication in Transport Theory and Statistical Physics

Chaos suppression in the large size limit for long-range systems

We consider the class of long-range Hamiltonian systems first introduced by Anteneodo and Tsallis and called the alpha-XY model. This involves N classical rotators on a d-dimensional periodic lattice interacting all to all with an attractive coupling whose strength decays as r^{-alpha}, r being the distances between sites. Using a recent geometrical approach, we estimate for any d-dimensional lattice the scaling of the largest Lyapunov exponent (LLE) with N as a function of alpha in the large energy regime where rotators behave almost freely. We find that the LLE vanishes as N^{-kappa}, with kappa=1/3 for alpha/d between 0 and 1/2 and kappa=2/3(1-alpha/d) for alpha/d between 1/2 and 1. These analytical results present a nice agreement with numerical results obtained by Campa et al., including deviations at small N.Comment: 10 pages, 3 eps figure

Coexistence of Josephson oscillations and novel self-trapping regime in optical waveguide arrays

Considering the coherent nonlinear dynamics between two weakly linked optical waveguide arrays, we find the first example of coexistence of Josephson oscillations with a novel self-trapping regime. This macroscopic bistability is explained by proving analytically the simultaneous existence of symmetric, antisymmetric and asymmetric stationary solutions of the associated Gross-Pitaevskii equation. The effect is, moreover, illustrated and confirmed by numerical simulations. This property allows to conceive an optical switch based on the variation of the refractive index of the linking central waveguide.Comment: 4 pages, 4 figure

Controversy about the applicability of Tsallis statistics to the HMF model

Comment to "Nonextensive Thermodynamics and Glassy Behaviour in Hamiltonian Systems" by A. Rapisarda and A. Pluchino, Europhysics News 36, 202 (2005)

Modulational Estimate for Fermi-Pasta-Ulam Chain Lyapunov Exponents

In the framework of the Fermi-Pasta-Ulam (FPU) model, we show a simple method to give an accurate analytical estimation of the maximal Lyapunov exponent at high energy density. The method is based on the computation of the mean value of the modulational instability growth rates associated to unstable modes. Moreover, we show that the strong stochasticity threshold found in the $\beta$-FPU system is closely related to a transition in tangent space: the Lyapunov eigenvector being more localized in space at high energy.Comment: 4 pages, revtex, 4 ps figures, submitted to PR

Dynamics and statistics of simple models with infinite-range attractive interaction

In this paper we review a series of results obtained for 1D and 2D simple N-body dynamical models with infinite-range attractive interactions and without short distance singularities. The free energy of both models can be exactly obtained in the canonical ensemble, while microcanonical results can be derived from numerical simulations. Both models show a phase transition from a low energy clustered phase to a high energy gaseous state, in analogy with the models introduced in the early 70's by Thirring and Hertel. The phase transition is second order for the 1D model, first order for the 2D model. Negative specific heat appears in both models near the phase transition point. For both models, in the presence of a negative specific heat, a cluster of collapsed particles coexists with a halo of higher energy particles which perform long correlated flights, which lead to anomalous diffusion. The dynamical origin of the "superdiffusion" is different in the two models, being related to particle trapping and untrapping in the cluster in 1D, while in 2D the channelling of particles in an egg-crate effective potential is responsible of the effect. Both models are Lyapunov unstable and the maximal Lyapunov exponent $\lambda$ has a peak just in the region preceeding the phase transition. Moreover, in the low energy limit $\lambda$ increases proportionally to the square root of the internal energy, while in the high energy region it vanishes as $N^{-1/3}$.Comment: 33 pages, Latex2 - 12 Figs - Proceedings of the Conference "The Chaotic Universe" held in Rome-Pescara in Feb. 199