Ergodicity Breaking and Parametric Resonances in Systems with Long-Range Interactions Fernanda P. da C. Benetti, Tarcísio N. Teles, Renato Pakter, and Yan Levin Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 1551, CEP 9151-97, Porto Alegre, RS, Brazil (Dated: February 13, 212) arxiv:122.181v2 [cond-mat.stat-mech] 1 Feb 212 We exlore the mechanism resonsible 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 traed in an out-of-equilibrium quasi-stationary-state. Nevertheless, we show that if the initial distribution satisfies a secific 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, arametric resonances are excited, leading to chaos and ergodicity breaking. PACS numbers: 5.2.-y, 5.7.Ln, 5.45.-a Statistical mechanics of systems in which articles interact through long-ranged otentials is fundamentally different from the statistical mechanics of systems with short-range forces [1]. In the latter case, starting from an arbitrary initial condition (microcanonical ensemble) systems evolve to a thermodynamic equilibrium in which article distribution functions are given by the usual Boltzmann-Gibbs statistical mechanics [2]. The state of thermodynamic equilibrium does not deend on the secifics of the initial distribution, but only on the global conserved quantities such as energy, momentum, angular momentum, etc. The situation is very different for systems in which articles interact through long-range otentials, such as gravity or unscreened Coulomb interactions [3 6]. In this case, it has been observed in numerous simulations that these systems do not relax to thermodynamic equilibrium, but become traed in a quasi-stationary state (qss), the lifetime of which diverges with the number of articles [4, 6 8]. The distribution functions in this quasi-stationary state do not obey the Boltzmann-Gibbs statistical mechanics and in articular, article velocities do not follow the Maxwell-Boltzmann distribution, but deend exlicitly on the initial condition. It has been an outstanding challenge of statistical mechanics to quantitatively redict the final stationary state reached by systems with unscreened long-range forces, without having to exlicitly solve the N-body dynamics or the collisionless Boltzmann (Vlasov) equation. Some 4 years ago Lynden-Bell (LB) roosed a generalization of the Boltzmann-Gibbs statistical mechanics to treat systems with long-range interactions [9]. Lynden- Bell s construction was based on the Boltzmann counting, but instead of using articles, LB worked directly with the levels of the distribution function. The motivation for this aroach was the observation that dynamical evolution of the distribution function for systems with long-range interactions is governed by the Vlasov equation [1]. This equation has an infinite number of conserved quantities, Casimirs any local functional of the distribution function is a Casimir invariant of the Vlasov dynamics. In articular if the initial distribution function is discretized into levels, the volume of each level must be reserved by the Vlasov flow. For an initially one-level distribution function, Vlasov dynamics requires that the hase sace density does not exceed that of the initial distribution one-article distribution function over the reduced hase sace (µ-sace) evolves as an incomressible fluid. Using this constraint in a combination with the Boltzmann counting, LB was able to derive a coarse-grained entroy, the maximum of which he argued should corresond to the most-robable distribution the one that should describe the equilibrium state. Numerous simulations, however, showed that, in general, Lynden-Bell statistics was not able to account for the article distribution in self-gravitating systems, and the theory has been abandoned in the astrohysical context. Recently, however, Lynden-Bell s work has been rediscovered by the Statistical Mechanics community, which showed that for some systems, secifically the widely studied Hamiltonian Mean Field Model (HMF), Lynden-Bell s aroach could make reasonable redictions about the structure of the hase diagram [11]. The fundamental question that needs to be addressed is: Under what conditions can Lynden-Bell statistics be used to accurately describe systems with long-range interactions? This will be the toic of the resent Letter. To be secific, we will study the HMF model [1], which has become a test bench for theories of systems with longrange forces. However, our results and methods are comletely general and can be alied to other systems, such as self-gravitating clusters or confined non-neutral lasmas. The HMF model consists of N articles restricted to move on a circle of radius one. The dynamics is governed by the Hamiltonian H = N i=1 2 i 2 + 1 2N N [1 cos( i j )], (1) i,j=1
2 where the angle i is the osition of i th article and i is its conjugate momentum [11 13]. The macroscoic behavior of the system is characterized by the magnetization vector M = (M x, M y ), where M x cos, M y sin, and stands for the average over all articles. The Hamilton s equations of motion for each article reduce to i = M x (t) sin i (t) + M y (t) cos i (t). (2) Since the Hamiltonian does not have exlicit time deendence, the average energy er article, u = H N = 2 2 + 1 M(t)2, (3) 2 is conserved. The failure of LB theory in the astrohysical context was attributed to incomlete relaxation, lack of good mixing, or broken ergodicity [14]. The mechanisms behind this failure have not been elucidated. On the other hand, it has been recently observed that if the initial distribution is virialized satisfies the virial condition LB s aroach was able to quite accurately redict the stationary state of gravitational and Coulomb systems [3 6]. Unfortunately, the virial theorem can be derived only for otentials which are homogeneous functions. This is not the case for the HMF model. Nevertheless, the fact that LB theory seems to aly under some conditions makes one wonder if such conditions can be found for arbitrary long-range otentials, which are not in general homogeneous functions. To answer the questions osed above, we note that if the initial distribution is virialized, macroscoic oscillations of observables should be diminished. On the other hand, if the system is far from virial, the mean-field otential that each article feels will undergo strong oscillations. It is then ossible for some articles to enter in resonance with the oscillations of the mean-field, gaining large amounts of energy. The arametric resonances will result in the occuation of regions of the hase-sace which are highly imrobable, from the oint of view of Boltzmann-Gibbs or LB statistics [16]. Furthermore, resonant articles will take away energy from collective oscillations roducing a form of non-linear Landau daming [15]. After some time, macroscoic oscillations will die out and each article will feel only the static meanfield otential. From that oint on, article dynamics will become comletely regular, with no energy exchange ossible between the different articles. The articles which have gained a lot of energy from the arametric resonances will be traed forever in the highly imrobable regions of the hase-sace, unable to thermalize with the rest of the system. To see how the theoretical icture advocated above can be alied to the HMF, we first derive a generalized virial condition for this model. For simlicity we will consider initial distributions of the water-bag form in (, ). Without loss of generality, we choose a frame of reference where = and =. The one-article initial distribution function then reads f (, ) = 1 4 Θ( ) Θ( ), (4) where Θ is the Heaviside ste function, and and are the maximum values of angle and momentum, resectively. Note that from symmetry, M y (t) = at all times. When the dynamics starts, the mean-squared article osition will evolve with time. We define the enveloe of the article distribution as e (t) = 3 2, so that e () =. We next differentiate e (t) twice with resect to time to obtain the enveloe equation of motion, e = 3 2 e + 3 e 9 2 e 3. (5) Using the conservation of energy, 2 = 2u + M 2 x(t) 1. To calculate,, we use the equation of motion for. Suosing that the distribution of angles remains close to uniform on the interval [ e, e ], we obtain = M x(t) 2 e e(t) e(t) where the magnetization M x (t) is sin d = M x (t) cos e (t) M 2 x(t), (6) M x (t) = 1 e(t) d cos 2 e e(t) = sin e(t). (7) e (t) Neglecting the correlations between ositions and velocities, =, we finally obtain a dynamical equation for the enveloe e = 3 e (t) (2u + M x(t) cos e (t) 1), (8) where u = 2 /6 + (1 M 2 )/2 and M = sin( )/. The generalized virial condition is defined by the stationary enveloe, e =, which means that along the curve (2u 1) + sin cos =. (9) magnetization remains aroximately invariant. In Fig. 1 we lot Eq.(9) in the M u lane and comare it with the full molecular dynamics simulation of the HMF model. As can be seen, agreement between the theory and simulation is excellent. Along the generalized virial condition curve, Eq. (9), the magnetization and, therefore, the mean-field otential acting on each article of the HMF model has
3 1.8 Para Virial condition Simulation Phase transition The Lagrange multiliers µ and β are determined by article and energy conservation, d df(, ) = 1, (12) u.6.4.2 Ferro.2.4.6.8 1 M FIG. 1. Phase diagram of the HMF model obtained using the molecular dynamics simulations. Solid curve is the line of first order hase transitions searating aramagnetic and ferromagnetic hases. This line extends u to M =.6, after which oint the order of the hase transition, shaded region (green line), becomes unclear, with strong deendence on the initial conditions and various reentrant transitions occurring in this region. Dashed curve is the generalized virial condition, Eq.(9). Along this curve oscillations of the enveloe are suressed. Diamonds are the results of simulation. Starting with the initial energy and magnetization along the virial curve, diamonds show the final magnetization to which the system relaxes. For oints along this curve, the final magnetization is almost identical to the initial one. Note that the generalized virial curve terminates at M =.34 slightly below the hase transition line. This small difference, however, is sufficient to invalidate the Lynden-Bell theory, which for M =.4 redicts a second order hase transition, while the simulations show that the hase transition is of first order [17]. Points (A), (B), and (C) corresond to the initial conditions for the distribution functions shown in Fig. 2. The Poincaré sections of the test article dynamics for the initial conditions described by the oints (B) and (T) are shown in Fig. 3. Finally, we note that since the stationary distribution must satisfy the virial condition and the energy is conserved, Eq.(9) allows us to redict the magnetization to which the system will evolve for initial conditions lying inside the ferromagnetic region, see the arrows for oints (A) and (C). only microscoic oscillations and the arametric resonances are suressed. Under these conditions, we exect that LB theory will be valid. The coarse grained entroy within the LB aroach is given by s(f) = dd [ f η ln f η + (1 fη ) ln A T B C (1 fη )], (1) where η = 1/4 [18]. Maximizing this entroy under the constraints of energy and article conservation, we obtain the equilibrium distribution function f(, ) = η [ e β 2 ]. (11) 2 Mx cos µ +1 [ 2 d d f(, ) 2 1 ] 2 (1 M x cos ) = u, (13) resectively, and the magnetization by the selfconsistency requirement, d d cos f(, ) = M x. (14) Solving these equations numerically along the curve Eq.(9), we see that there is an excellent agreement between LB theory and the simulations, Fig. 2. If the macroscoic oscillations are suressed and the arametric resonances are not excited, the system is able to relax to a quasi-ergodic equilibrium ermitted by the Vlasov dynamics. To make clear the role of arametric resonances in ergodicity breaking, in Fig.3a we lot the Poincaré section of a set of non-interacting test articles, which at t = are distributed in accordance with Eq.(4). The motion of each article is governed by Eq.(2) with M x determined by Eqs. (7) and (8). The osition and momentum of each article are lotted when magnetization is at its minimum. We see that if the energy and the initial magnetization lie on the generalized virial curve oint (B) of Fig. 1 article trajectories are comletely regular. However, when initial conditions do not coincide with the generalized virial curve oint (T) of Fig. 1 arametric resonances aear and dynamics becomes chaotic. Particles enter in resonance with the oscillations of the mean-field otential, gaining sufficient energy to move into statistically imrobable regions of the hase sace. The Poincaré section of test article dynamics is remarkably similar to the final stationary distribution obtained using the comlete N-body molecular dynamics simulation of the HMF, Fig. 3. Eq. (8) can also be used to calculate the eriod of the first oscillation of M(t). For examle, for oint (T) of the hase diagram Fig. 1, we find the eriod to be T = 5., while the full molecular dynamics simulations gives T = 5.4. For oint (C) we find T = 3.85, while the simulations give T = 3.82. In conclusion, we have studied the mechanism resonsible for the ergodicity breaking in systems with longrange interactions. Ergodicity breaking and the arametric resonances are intimately connected. If the macroscoic oscillations and the resulting resonances are suressed the system is able to relax to a quasiergodic stationary state. However, when the arametric resonances are excited, some articles are ejected to statistically imrobable regions of the hase sace, at the same time as the oscillations are damed out. The
4 1-1 2. 1 (a) N() N() 1 1-1 1-2 1-3 -2 1 2. 1. 1 1. 1 (c) 2. 1-2 1-4 1 1-1 -3 1 2. 1. 1 (e) 1... -1. -1. -2.. 1. 2. -2. (b) 1... -1. -1. 2. (c) 1. -2. 2. 1-2 1-3 1. (d) -1 N() N() 1. 2. (a) (b) -2. (d) (f) N() N() FIG. 3. Poincare sections of test articles and snashots of the hase sace obtained using molecular dynamics simulation once the system has relaxed to qss. Panels (a) and (b) corresond to the initial condition lying on the generalized 1-1 1-1 virial curve, oint (B) of Fig. 1. In this case the test article dynamics is comletely regular, and the stationary article distributions are well described by LB theory. Panels (c) and (d) corresond to the initial conditions slightly off the virial -2 1-2 1 2. 1.. 1. 2. curve, oint (T) of Fig. 1. Even though we have moved only a little from the virial curve, we see the aearance of resonant islands and the dynamics of some of the test articles becoming chaotic. Such resonances drive some articles FIG. 2. The angle and velocity distribution functions corof the HMF to statistically imrobable from the oint of resonding to the initial conditions described by oints (A), view of the Boltzmann-Gibbs and LB statistical mechanics (B) and (C) of Fig. 1, resectively. Symbols are the results of regions of the hase sace. Once the enveloe oscillations molecular dynamics simulations and solid curves are the reare damed out, article dynamics becomes comletely intedictions of LB theory. The simulated distribution functions grable, and there is no mechanism for the resonant articles for the oint (B), lying on the generalized virial curve, are to equilibrate with the rest of the distribution. Thus, away in excellent agreement with the redictions of the LB theory from the generalized virial curve, ergodicity becomes broken. [anels (c) and (d)], demonstrating that the dynamics along the generalized virial curve is quasi-ergodic. On the other hand, the distribution functions for oints (A) and (C) deviate significantly from the redictions of LB theory showing that away from the generalized virial curve, ergodicity is broken [anels (a), (b) and (e), (f)]. rocess of continuous article ejection, and the resulting decrease of macroscoic oscillations of the enveloe, leads to the formation of a static mean-field otential and to asymtotically integrable dynamics. Once the integrability of the equations of motion is achieved, the ergodicity becomes irreversibly broken. Unlike for articles with short-range interaction otentials, ergodicity is the excetion rather than the rule for systems with long-range forces it can only be observed if the initial distribution function satisfies the generalized virial condition derived in this Letter. Finally we note, that since the stationary distribution must satisfy the virial con- dition and the energy must be conserved, Eq.(9) allows us to redict the magnetization to which the system will evolve for initial conditions lying inside the ferromagnetic region. For examle, oint (A) of Fig. 1 which has initial magnetization and energy M =.74 and u =.55, will evolve to a final stationary state with M =.56; while the oint (C) with M =.74 and u =.25, will evolve to a final stationary state with M =.86, which are recisely the values obtained using the molecular dynamics simulations. This work was artially suorted by the CNPq, FAPERGS, INCT-FCx, and by the US-AFOSR under the grant FA955-9-1-283.
5 [1] A. Cama, T. Dauxois, and S. Ruffo, Phys. Re. 48, 57 (29). [2] J. W. Gibbs, Collected Works, Longmans, Green and Co., NY (1928). [3] Y. Levin, R. Pakter, and T. N. Teles, Phys. Rev. Lett. 1, 464 (28). [4] T.N. Teles, Y.Levin, R. Pakter, and F.B. Rizzato, J. Stat. Mech. P57 (21). [5] M. Joyce and T. Worrakitoonon, Phys. Rev. E 84, 11139 (211). [6] Teles, T. N. and Levin, Y. and Pakter, R., Mon. Not. R. Astron. Soc. 417, 1, (211). [7] J. Barré, D. Mukamel, and S. Ruffo, Phys. Rev. Lett. 87, 361 (21). [8] K. Jain et al., J. Stat. Mech. (27) P118. [9] D. Lynden-Bell, Mon. Not. R. Astron. Soc. 136, 11 (1967). [1] W. Braun and K. He, Comm. Math. Phys. 56, 11 (1977). [11] A. Antoniazzi, F. Califano, D. Fanelli, and S. Ruffo, Phys. Rev. Lett. 98. 1562 (27); J. Barre, F. Bouchet, T. Dauxois, S. Ruffo, Phys. Rev. Lett. 89, 1161 (22); T. M. Rocha Filho, A. Figueiredo, M. A. Amato, Phys. Rev. Lett. 95. 1961 (25). [12] A. Antoniazzi, D. Fanelli, and S. Ruffo, Y.Y. Yamaguchi Phys. Rev. Lett. 99. 461 (27). [13] A. Antoniazzi, D. Fanelli, J. Barré., P.H. Chavanis, T. Dauxois and S. Ruffo, Phys. Rev. E 75, 11112 (27); [14] D. Mukamel, and S. Ruffo, and N. Schreiber, Phys. Rev. Lett. 95, 2464 (25). [15] L. Landau, J. Phys. USSR 1, 25 (1946). [16] R. Bachelard, C. Chandre,D. Fanelli, X. Leoncini, and S. Ruffo, Phys. Rev. Lett. 11. 2663 (28). [17] R. Pakter and Y. Levin, Phys. Rev. Lett. 16, 263 (211). [18] F. Staniscia, P.H. Chavanis, G. De Ninno, D. Fanelli, Phys. Rev. E 8, 21138 (29).