arxiv: v2 [cond-mat.stat-mech] 27 Dec 2012

Transcription

1 Overdamped dynamics of long-range systems on a one-dimensional lattice: Dominance of the mean-field mode and phase transition arxiv:29.638v2 [cond-mat.stat-mech] 27 Dec 22 Shamik Gupta, Alessandro Campa 2, and Stefano Ruffo,3 Laboratoire de Physique de l École Normale Supérieure de Lyon, Université de Lyon, CNRS, 46 Allée d Italie, Lyon céde 7, France 2 Comple Systems and Theoretical Physics Unit, Health and Technology Department, Istituto Superiore di Sanità, and INFN Roma, Gruppo Collegato Sanità, Viale Regina Elena 299, 6 Roma, Italy 3 Dipartimento di Energetica Sergio Stecco and CSDC, Università di Firenze, CNISM and INFN, via S. Marta 3, 539 Firenze, Italy (Dated: January, 23) 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 /r α ; < α <. 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 = /2 and uniform above, with a phase transition between the two at T = /2. The state is the same as the equilibrium state of the mean-field version of the model, obtained by considering α =. 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. PACS numbers: 5.2.-y, 5.4.-a, 5.7.Fh I. INTRODUCTION In long-range interacting systems, the interparticle potential in d dimensions decays at large separation, r, as /r α, with α d [ 3]. Eamples are gravitational systems [4], plasmas [5], two-dimensional hydrodynamics [6], charged and dipolar systems [7], etc. Long-range systems are non-additive, and have equilibrium properties unusual for short-range systems, e.g., inequivalence of statistical ensembles [8]. Besides, they show intriguing dynamical features, e.g., broken ergodicity [9] and slow relaation towards equilibrium [9, ]. A paradigmatic eample of long-range systems is the so-called Hamiltonian mean-field (HMF) model, comprising a system of N interacting particles evolving under Hamilton dynamics []. The ith particle, i =,2,...,N, has an internal degree of freedom θ i, a periodic variable of period 2π. The θ i s are coupled to one another via an attractive XY-like interaction cos(θ i θ j ) between the ith and jth particles, with the coupling strength being of the mean-field type: It is equal for every pair of particles. As a result, the model is defined without requiring any lattice structure for the particles to reside on. In this model, a wide class of initial distributions relaes to the Boltzmann-Gibbs equilibrium state over times that diverge with the system size shamikg@gmail.com alessandro.campa@iss.infn.it stefano.ruffo@gmail.com []. As a consequence, the system in the thermodynamic limit does not ever attain equilibrium, instead it remains trapped in intermediate quasistationary states. Besides the HMF model, these states have been observed in a wide variety of long-range systems ranging from models with spin dynamics [2] and plasma physics [3] to onedimensional and two-dimensional gravity [4, 5]. The HMF model is a prototype for systems with unscreened long-range interactions leading to long-lived spatially inhomogeneous states. Moreover, the model does not have a singularity at short length scales so that one can eclusively focus on effects of long-range interactions, in a way similar to studies of gravitational systems where a short-range softening is introduced [6]. To eplore how far the results obtained within the HMF model etend beyond its setting of mean-field coupling, the model has been considered on a onedimensional periodic lattice on the sites of which the particles reside. Particles on different sites are coupled as in the HMF model, with the difference that the coupling strength decays with their separation r along the lattice length as /r α, with < α <. The resulting model is known as the α-hmf model [7]. The HMF model is recovered on considering α =. Within a canonical ensemble, the α-hmf model has the same equilibrium state as those of the HMF model [8]. There have been many studies dealing with various dynamical aspects of the α-hmf model that have illustrated their difference from those of the HMF model, e.g., Lyapunov eponents [9,2], andothers[2,22]. Hence, the equivalenceofthe equilibrium state of the two models raises an immediate question: How does the dynamics of the α-hmf model

2 2 lead to the same equilibrium state as that of its meanfield counterpart, the HMF model? This work is a step towards answering this question regarding the mean-field dominance in dictating the equilibrium state. To address the issue, let us consider the simple setting of the α-hmf model evolving in the presence of thermal noise and dissipation. The ith particle has the equations of motion dθ i dt = p i, () dp i dt = γp i + N Ñ sin(θ j θ i ) ( j i c ) α + γη i (t), (2) j=,j i wherep i denotes the momentum ofthe ith particle, while γ is the damping constant. The second term on the right handside ofeq. (2) isthe forceon theith particlearising from the interactionpotential: The quantity j i c is the closest distance on the periodic lattice between the ith and jth sites: while j i c = min( j i,n j i ), (3) Ñ = N ( j i c ) α i. (4) j= From Eq. (2), since α <, the cumulative interaction of one particle with all the remaining particles with aligned θ s would diverge in the limit N in the absence of the normalizationñ, which eplains its inclusion [8]. In Eq. (2), η i (t) is a Gaussian white noise: η i (t) =, (5) η i (t)η j (t ) = 2Tδ ij δ(t t ), (6) where T denotes the temperature, while the angular brackets denote averaging over noise realizations. In this work, we take the Boltzmann constant to be unity. The equations of motion () and (2) describe the evolution of the α-hmf model within a canonical ensemble. We note that for α =, these equations reduce to those of the Brownian mean-field model considered in Refs. [23] and [24]. Here, we will study the overdamped limit of the equations of motion, () and (2); after a time rescaling, these equations then reduce to the following Langevin equation for the ith particle: dθ i dt = Ñ N j=,j i sin(θ j θ i ) ( j i c ) α +η i(t). (7) We analyze the dynamics (7) in the continuum limit N, when the lattice becomes a continuous segment characterized by the spatial coordinate s [,]. In this limit, the local density of oscillators ρ(θ; s, t) evolves in time following a Fokker-Planck equation. We show that the Fokker-Planck equation allows a stationary state which is uniform in both θ and s. By performing a linear stability analysis of such a state, we find that when it is unstable, different spatial Fourier modes of fluctuations have different critical temperatures below which the modes grow eponentially in time with different rates. The largest critical temperature, T c, = /2, corresponds to the spatially independent zero Fourier mode, i.e. the mean-field mode. Above this temperature, all the Fourier modes decay in time, thereby stabilizing the uniform state. Below T c,, our numerical simulations starting from the uniform state show that the unstable non-zero Fourier modes, growing eponentially in time at short times, nevertheless decay at long times to zero. By contrast, the mean-field mode grows and reaches a non-zero value, corresponding to a non-uniform stationary state, i.e. a state which is uniform in s, but nonuniform in θ. We perform a linear stability analysis of the non-uniform state to demonstrate that indeed at temperaturesbelowt c,, allthefouriermodesoffluctuations decay in time, thereby stabilizing the non-uniform state. The long time dominance of the mean-field mode leading to a stationary state always uniform in space, here observed in the alpha-hmf model evolving under the canonical and overdamped dynamics, has also been observed when the model evolves within a microcanonical ensemble. In the latter case, the dynamical equations are obtainedfromeqs. ()and(2)bysubstitutingγ = [25]. In the case of microcanonical dynamics, this dominance was just a numerical observation with no theoretical justification, while within the canonical and overdamped dynamics, we show how the mean-field mode arises as a linearly stable stationary solution of the Fokker-Planck equation. The paper is structured as follows. In section II, we introduce the continuum limit of the dynamics (7), and perform the linear stability analysis of the uniform stationary state. In section III, we present numerical simulations in support of the analysis, in particular, to show the agreement of the growth rates of the unstable modes, the decayatlongtimes ofthe unstablenon-zeromodesat all temperatures T < /2, and the associated dominance ofthemean-fieldmode. InsectionIV,weperformthelinear stability analysis of the non-uniform stationary state. We end the paper with conclusions and perspectives. II. CONTINUUM LIMIT: UNIFORM STATIONARY STATE AND DYNAMICS OF FOURIER MODES In the continuum limit N, the lattice of the system (7) is densely filled with sites. Let us introduce the variable s = i/n to denote the spatial coordinate along the lattice length, such that in the continuum limit it becomes the continuous variable s [, ]. In this limit, we define a local density of particles ρ(θ;s,t) such that of all particles located between s and s+ds at time

3 3 t, the fraction ρ(θ;s,t)dθ have their degrees of freedom between θ and θ +dθ. This density is non-negative, 2πperiodic in θ, and obeys the normalization 2π dθ ρ(θ;s,t) = s. (8) In the continuum limit, the equation of motion is given by θ(s, t) = κ(α) dθ ds sin(θ θ) t ( s s c ) αρ(θ ;s,t)+η(s,t), (9) where and Also, we have κ(α) = s+/2 s /2 ds ( s s c ) α, () s s c = min( s s, s s ). () η(s,t) =, (2) η(s,t)η(s,t ) = 2Tδ(s s )δ(t t ). (3) The time evolution of ρ(θ;s,t) follows the Fokker- Planck equation ρ t = κ(α) [( dθ ds sin(θ θ) ) ] θ ( s s c ) αρ(θ ;s,t) ρ +T 2 ρ θ2. (4) Its stationary solution, obtained by setting the left hand side to zero, is [ κ(α) ρ (θ;s) ep dθ ds cos(θ θ) ] T ( s s c ) αρ (θ ;s ). (5) Consider the particular stationary state which is uniform both in θ and in the spatial coordinate s, ρ = 2π. (6) This state, although stationary, might be destabilized by the thermal fluctuations inherent in the dynamics. Let us then analyze the stability, in particular, linear stability of the state with respect to small thermal fluctuations δρ(θ;s,t). To this end, we write ρ(θ;s,t) = +δρ(θ;s,t); δρ(θ;s,t). (7) 2π Substituting Eq. (7) into Eq. (4) and keeping terms to linear order in δρ, we find that δρ satisfies δρ t = κ(α) dθ ds cos(θ θ) 2π ( s s c ) αδρ(θ ;s,t)+t 2 δρ θ 2. (8) Epressing δρ in terms of its Fourier series with respect to the periodic variable θ as δρ(θ;s,t) = k= δρ k (s,t)e ikθ, (9) we find from Eq. (8) that only the modes k = ± are affected by the coupling between the particles, and that δρ ± satisfies δρ ± t = κ(α) 2 ds δρ± (s,t) ( s s c ) α T δρ ±. (2) One thus gets a set of equations for each position s, all coupled together by the second term on the right hand side of Eq. (2). For k ±, one has δρ ±k t = T δρ ±k ; k. (2) Thus, δρ±k (s,t) for k decays eponentially in time as ep( T t), so that these modes cannot destabilize the uniform state (6). Since we have a periodic lattice, to solve Eq. (2), we net consider the spatial Fourier series of δρ ± (s,t): δρ ± (s,t) = m= On substituting in Eq. (2), we get where δρ ±,m t Λ m (α) = δρ ±,m (t)e i2πms. (22) = κ(α)λ m(α) δρ 2 ±,m Tδρ ±,m, (23) s+/2 s /2 ds e i2πm(s s) ( s s c ) α. (24) Note that Λ m (α) = Λ m (α). It is known that Λ m (α), and that it is a monotonically decreasing function of m when m is an integer; moreover, Λ m (α) as m ± [26]. Equation (23) has the solution [( κ(α)λm (α) ) ] δρ ±,m (t) = ep T t δρ 2 ±,m ().(25) It then follows from our linear stability analysis that the mth spatial mode of fluctuations δρ ±,m (t) either (i) decays in time, which happens at temperatures such that T > κ(α)λm(α) 2, or (ii) grows in time at temperatures T < κ(α)λm(α) 2. The borderline between these two behaviors is achieved at the critical temperature T c,m for the mth mode, given by T c,m = κ(α)λ m(α). (26) 2 The mth mode at T < T c,m grows in time as ep[(t c,m T)t]. With Λ m (α) = Λ m (α), the Fourier modes m and

4 4 m have the same behavior in time and the same critical temperature. Since Λ m (α) decreases on increasing m, we conclude that T c, > T c, > T c,2 >... (27) Note that T c, is the critical temperature for the m = mode, i.e., the mean-field mode; we have T c, = κ(α)λ (α) 2 = 2. (28) Thus, abovet = /2,whenallthespatialmodesdecayin time, the uniform state (6) isstable. On the otherhand, the state is unstable when T < /2. At temperatures between /2 and T c,, the m = mode grows in time, while all the other modes decay in time. At temperatures betweent c,2 andt c,, the m =,±modesgrowin time, while all the other modes decay in time, and so on. In Table I, we show representative numbers for these critical temperatures for α =.25,.5,.75. α T c, T c,2 T c,3 T c,4 T c, TABLE I. Critical temperatures T c,m;m =,...,5 for α =.25,.5,.75. Note that T c, = /2, independent of α. When α =, one has T c,m = for m. When α approaches unity from below, all the T c,m s for m approach T c,. III. COMPARISON WITH NUMERICAL SIMULATIONS In order to test the theoretical predictions in the limit N obtained in the preceding section regarding values of critical temperatures, and growth and decay of Fourier modes depending on the temperature, we performed etensive numerical simulations of the dynamics for large N. A standard procedure for simulations is to integrate the equation of motion (7) for each of the N particles; this involves computing at every integration step a sum over N terms for each particle, thus requiring a total computation time scaling with N as N 2. In order to perform faster simulations, we adopted the algorithm discussed in Ref. [26] to compute the sum in Eq. (7). This algorithm requires a total computation time that scales with N as N lnn. The equations of motion were integrated using a fourth-order Runge-Kutta algorithm with time step equal to.. The initial state of the system was chosen to be the uniform state (6), prepared by having each particle degree of freedom θ uniformly distributed in [ π, π], independently of the others. We report here simulation results for system sizes 2 4 and 2 8. In simulations, we monitor the observable r m (t) = N e i(θj+2πjm/n) ;m =,,2,...; (29) N j= in particular, r (t) does not contain any spatial dependence, and, thus, characterizes the mean-field mode. In the continuum limit, we have r m (t) = dθds e i(θ+2πms) ρ(θ;s,t). (3) Starting at time t = from the uniform state such that r m () = m, and at a temperature T =. with α =.25, we epect from Table I that the Fourier modes,,2,3 in particular are linearly unstable. Consequently, r (t),r (t),r 2 (t),r 3 (t) are epected to grow in time. Indeed, Fig. (a) displaying the simulation results for r (t),r (t),r 2 (t),r 3 (t) does show that they all grow in time. However, in the long-time limit, we see that r (t) saturates to a value very close to unity, while r (t),r 2 (t),r 3 (t) all decay to a value very close to zero. By repeating simulations at different temperatures T < /2 for different values of α < (see Fig. (b),(c)), we have confirmed that r (t) always saturates in the long-time limit to a non-zero stationary state value r st = r st (T) that depends on the temperature, while r m (t) s, with m, decay in the long-time limit to a stationary state value very close to zero. As shown in Fig. 2, the asymptotic value of r (t) does not scale with the system size N, while that of r m (t) for m scales with the system size as / N, and hence approaches zero in the continuum limit N. These observations suggest that for temperatures T < /2 for all values of α <, the stationary state is characterized by a distribution ρ which is non-uniform in θ, but uniform in s, and may be obtained from Eq. (5) as: [ κ(α) ρ (θ) ep dθ ds cos(θ θ) ] T ( s s c ) αρ (θ ) ] = ep[ T (m cosθ +m y sinθ), (3) where (m,m y ) 2π dθ(cosθ,sinθ)ρ (θ). (32) IV. NON-UNIFORM STATIONARY STATE AND DYNAMICS OF FOURIER MODES Our observations regarding saturation of r (t) and decay of r m (t) s, with m, detailed in the preceding section, raise an immediate question: Can we show by performing a linear stability analysis of the state (3) at temperatures T < /2 that all the modes of fluctuations decay in the long-time limit to zero, so that the state is linearly stable for all values of α <?

5 5 Before proceeding with the stability analysis, we note that our system has full rotationalinvariance in θ at each space point, so that in the following, we take m y = and m without any loss of generality; then, we have ] ρ (θ) = Aep[ T m cosθ, (33) where A is the normalization: A = 2π dθe m T cosθ. (34) Note that the quantity m is determined self-consistently by combining Eq. (33) with Eq. (32); one gets m = I (m /T) I (m /T), (35) where I n () is the modified Bessel function of order n: I n () = 2π 2π dθ e cosθ cosnθ. (36) We thus have A = [2πI (m /T)]. To study the linear stability of the state (33), we write ρ(θ;s,t) = ρ (θ)+δρ(θ;s,t); δρ(θ;s,t), (37) where we epand δρ(θ;s,t) as δρ(θ;s,t) = m= δρ m (θ,t)e i2πms. (38) UsingEq. (38)inEq. (4), wefindthat δρ m (θ,t)satisfies at its leading order the equation δρ m t (sinθ δρ = m m ) +T 2 δρ m θ θ 2 λ m (α) θ([ dθ sin(θ θ) δρ m (θ,t) where we have defined ) ]ρ, (39) λ m (α) κ(α)λ m (α). (4) The eigenfrequencies associated with the linear equation (39) can be studied by looking for solutions of the form δρ m (θ,t) = δρ m (θ,ω)e iωt. (4) Before studying the spectrum of the eigenfrequency ω, it is instructive to analyze the eplicit solution δρ m (θ,) corresponding to the neutral mode ω =. We find from Eq. (39) that δρ m (θ,) satisfies m sinθ δρ ( ) m λ m (α) m y cosθ m sinθ +T δρ m θ ρ = C, (42) where C is a constant independent of θ, and ( m, m y ) = dθ(cosθ,sinθ) δρ m (θ,). (43) Equation (42) has the solution δρ m (θ,) = δρ m (,)e m T ( cosθ) + C T e m T cosθ θ dθ e m T cosθ + Aλ m(α) m y e m T cosθ sinθ T Aλ m(α) m e m T cosθ ( cosθ). (44) T The periodicity condition δρ m (θ+2π,) = δρ m (θ,) implies that C =. Normalization of ρ(θ;s,t) for all times implies, using Eq. (37) and the fact that ρ (θ) is normalized, that 2π dθ δρ m (θ,) = m. Then, from Eq. (44) with C =, we get δρ m (,) 2π = Aλ m(α) m T dθe m T ( cosθ) 2π dθe m T cosθ ( cosθ). (45) Using the above equation and Eq. (35) in Eq. (44) with C =, wearriveatthe followingepressionfor δρ m (θ,): δρ m (θ,) = Aλ m(α) e m T cosθ [ m (cosθ m )+ m y sinθ].(46) T The quantities m,y in Eq. (46) are determined selfconsistently. They cannot both be equal to zero, as would have happened if δρm (θ,) were a constant, but this is not possible as the integral of δρ m (θ,) from to 2π is required to vanish. Using Eqs. (43), (46), (34) and (35), we get the following two equations: m y = m y λ m (α) m I (m /T) I (m /T), (47) m = m λ m (α) T ( T m 2 ), (48) where < T < /2, and m satisfies Eq. (35). Equation (47) is obviously satisfied for m y =. If m y, then, using Eq. (35), we see that Eq. (47) is an identity for m =, while it has no solution for m, which follows from the property that λ m (α) < for m. Let us now analyze Eq. (48). It is satisfied for m =. If m, it gives m = T T/λ m (α), (49) giving a relation between m and T, which has to be satisfied together with Eq. (35). To check whether the two equations are consistent, note that λ m (α) = λ m (α) and that λ m (α) is a monotonically decreasing function

6 6 of m, with λ (α) = and lim m ± λ m (α) =. Hence, to check the consistency, we have to consider in Eq. (49) any value for λ m (α) between and. Actually, when λ m (α) <, Eq. (49) gives a real value for m only for T, which is smaller than the range < T < /2. In Fig. 3, we plot m (T) as given by Eq. (49) for λ m (α) =, and as given by solving the implicit equation (35). We see that the two curves do not intersect. This also implies that there cannot be any intersection for λ m (α) <, since the right hand side of Eq. (49) decreases on decreasing λ m (α). Thus, we conclude that Eq. (49) for < T < /2 is not solved by m that satisfies Eq. (35). In summary, the two equations, (47) and (48), are not satisfied by m given by Eq. (35), when m. When m =, there is only one solution, namely, m =, and m y any non-zero value. For this solution, on using Eqs. (33) and (46), we have in the range < T < λm(α) +λ m(α) δρ (θ,) = T ρ (θ) m y sinθ, (5) which is just a perturbation of ρ (θ) that corresponds to a rotation in θ space of all the rotators of the system. It is therefore natural that it corresponds to a neutral mode. On the basis of the above discussion, we conclude that there is no temperature T < /2 such that the mth mode of fluctuations, for any m, has a vanishing eigenfrequency, ecepting for the trivial one for m = that corresponds to a state obtained by a rotation in θ space of all the rotators of the system. Thus, the mth mode either grows or decays in time at all temperatures T < /2, corresponding to eigenfrequencies ω which have respectively a negative or a positive imaginary part. In order to eamine the entire spectrum of the eigenfrequencyω,weuseeqs. (39)and(4), andthesubstitution µ iω to get the following equation for δρ m (θ,ω). (sinθ δρ µ δρ m = m m ) θ +T 2 δρ m θ 2 λ m (α) θ([ dθ sin(θ θ) δρ m (θ,ω) Let us define ( m (n), m(n) y ) = m (n) = dθ(cosnθ,sinnθ) δρ m (θ,ω), dθcosnθ ρ (θ) = I n(m /T) I (m /T), ) ]ρ. (5) (52) where we note that m () = m. On multiplying Eq. (5) in turn by cosnθ and sinnθ, integrating over θ from to 2π, noting that ρ is even in θ, that δρ m (2π,ω) = δρ m (,ω), and that ρ (2π) = ρ (), we arrive at the following system of equations for n : µ m (n) = ( ) 2 nm m (n ) m (n+) µ m (n) y = 2 nm + 2 λ m(α)n m () ( m (n ) y + 2 λ m(α)n m () y ( m (n ) m (n+) y ( m (n ) (n) Tn 2 m ) m (n+) ) (n) Tn 2 m y ) +m (n+), (53), (54) where m () = m () y =, and m () =. The operators on the right-hand sides of the above equations being nonhermitian have in general both real and comple eigenvalues. for all n N, then, combined with Eq. (35), the eigenvalues of the system of equations, (53) and (54), can be obtained numerically after truncating the equations to a finite value of n, say, n ma. As argued before on the basis of properties of λ m (α), in order to probe the behavior of a m mode of fluctuations, it suffices to consider any real positive value < for λ m (α), while for the m = mode, we have λ m (α) =. The results for the eigenvalue spectrum of Eqs. (53) and (54) are shown in Fig. 4 for the temperature T =.2. When T < /2 and m (n) When T /2 and m (n) = for all n N, the eigenvalue equations, (53) and (54), simplify for n 2 to µ m (n) = Tn µ m y (n) = Tn while for n =, the equations are µ m () 2 m (n) 2 m (n), (55) y, (56) = T m() + 2 λ m(α) m (), (57) µ m () y = T m () y + 2 λ m(α) m () y. (58) The system of equations, (55)-(58) is a set of independent equations, and may be solved quite easily to get the eact eigenvalues µ; from the equations, it is clear that these eigenvalues are the same for the and the y-system. Of course, the number of eigenvalues depends on the value n ma of n at which the equations are truncated. In particular, for λ m (α) = and T = /2, one finds the eigenvalue µ =. The results for the eigenvalue spectrum of Eqs. (55)-(58) are shown in Fig. 5 for the temperature T = /2. Let us first discuss the results displayed in Fig. 4. From the m mode results, displayed in panels (a) and (b), we see that the eigenvalues of the - and the y-system are both real and comple, but the important thing to note is that the real eigenvalues are all negative while the comple ones have strictly negative real parts. In order to demonstrate the former fact, a zoom into the region near the zero of the Re(µ) ais, shown in the insets, illustrates that the eigenvalue closest to of both the - and the y-system has a negative real part, so that the remaining eigenvalues having larger real parts

7 7 in magnitude are thus all negative. The insets also show that the eigenvalue with the smallest negative real part converges in magnitude on increasing the truncation order n ma. For higher n ma, only the number of eigenvalues increases; in particular, only the number of real eigenvalues with larger magnitude increases, while there is still only one comple eigenvalue which has the largest negative real part. Computing the eigenvalues for other λ m s and temperaturest < /2, we find that the number of comple eigenvalues is always small (with the number increasing for smaller λ m s and T s), and these eigenvalues always have large negative real parts. Since the mth mode of fluctuations has the time dependence e µt (recall Eq. (4) and the definition of µ = iω), and since all values for µ have strictly negative real parts, it follows that at temperatures T < /2, the mth mode of fluctuations, with m, decays in the long-time limit to zero. Figure 4, panel (c) shows that the behaviorfor the m = mode is very similar to that observed in panels (a) and (b) for m, discussed above. Panel (d) also has the same general behavior, ecepting that now the eigenvalue closest to the origin of the Re(µ) ais, shown in the inset as a function of n ma, is very slightly ( 7 ) positive. Actually, this eigenvalue should be eactly zero, being the one corresponding to rotational invariance that we discussed earlier (see the discussion following Eq. (49)). In fact, using Eq. (5) in Eq. (52) to compute m y (n), and using it in Eq. (54), it may be checkedthat its righthand side is zero for any n, that is, the particular solution (5) gives µ =. The slight positive value is just a numerical artifact. We now discuss the results displayed in Fig. 5. We show here only the eigenvalues of the -system, since, as noted earlier, the y-system has eactly the same eigenvalue spectrum. From the m mode results, displayed in panel (a), we see that the eigenvalues are only real and negative, with no comple eigenvalue. The eigenvalue closest to the origin of the Re(µ) ais, shown in the inset as a function of n ma, dominates the behavior of the mth mode of fluctuations in time, thereby making it decay to zero in the long-time limit. From the m = mode results, displayed in panel (b), we see that the eigenvalues are real and negative. As already mentioned above, the eigenvalue closest to the origin is precisely zero, which makes the zero mode of fluctuations to neither grow nor decay at the temperature T = /2. Based on results discussed in this section and in section II, we conclude that the linearly stable stationary state of the dynamics (7) in the N limit is a state uniform in space but non-uniform in θ when the temperature is less that /2. On tuning the temperature to a value larger than /2, the stationary state becomes uniform in θ. Thus, our linear stability analysis in the limit N allows us to predict that the system undergoes a transition at the temperature T = T c, = /2. Indeed, simulation results shown in Fig. 6 for r st suggests a continuous phase transition as a function of temperature, where we also show the theoretical curve for r st obtained by combining Eq. (3) with Eq. (33) to give for T < /2 the result r st = I (m /T) I (m /T) = m ; (59) for T > /2, on using Eq. (6), one has Eactly at T = /2, we have r st =. r st =. (6) V. CONCLUSIONS AND PERSPECTIVES In this paper, we considered a paradigmatic long-range interacting model of particles residing on the sites of a one-dimensional periodic lattice. Each particle has an internal degree of freedom θ which is coupled to those of other particles with an attractive XY-like interaction cos(θ i θ j ) between the ith and jth particles, with the coupling strength decaying with the interparticle separation r as /r α ; < α <. We considered the overdamped dynamics within a canonical ensemble of this socalled α-hmf model. We studied the model numerically, and also analytically in the continuum limit through the Fokker-Planck equation for the evolution of the local density of particles. The Fokker-Planck equation allows a stationary state which is uniform in both θ and space s of the lattice. A linear stability analysis of such a state shows that it is stable above the temperature T = /2. Below this temperature, when it is unstable, numerics as well as the linear stability analysis show that different spatial Fourier modes of fluctuations (the number of which depends on the temperature) grow eponentially in time with different rates. However, our numerical simulations also show that all the non-zero Fourier modes decay at long times to zero. By contrast, the zero mode (the mean-field mode) grows and reaches a non-zero value, corresponding to a non-uniform stationary state, i.e. a state which is non-uniform in θ, but uniform in s. On the basis of our investigations, it appears that both for temperatures below and above T = /2, the state which is uniform in space acts as a global attractor of all possible stationary states of the α-hmf model. The state is the same as the equilibrium state of the mean-field version of the model, the HMF model. For temperatures below /2, the late-time damping of the unstable non-zero modes after their initial growth is an interesting and unepected phenomenon. For such temperatures, by developing a linear stability analysis around the non-uniform state, we showed that all fluctuations around that state decay in time, thereby stabilizing the state; describing analytically the complete evolution of the unstable nonzero modes, from their initial growth in time to their late-time decay, thereby leading to the convergence to the non-uniform state remains an interesting challenging problem.

8 8 The model we studied may be related to the Kuramoto model in the field of non-linear dynamical systems [27 29]. Interpreting each of the N particles of our system to be a phase-only oscillator with phase θ i, the equation of motion (7) may be looked upon as the one governing the time evolution of the phase of the ith oscillator, in the presence of thermal noise. On considering α = in Eq. (7), the model in the absence of thermal noise is a particular limit of the Kuramoto model in which each oscillator has the same intrinsic frequency of oscillation. In the presence of thermal noise, Eq. (7) with α = is the generalization of the Kuramoto model considered in Ref. [3]. On considering α, and attributing to each oscillator an intrinsic frequency sampled from a distribution, the model in the absence of thermal noise has been considered in [26, 3 33]. In particular, a meanfield dominance similar to the one observed in this work hasbeenseeninnumericsinref. [26]. Itremainsanopen problem to relate this dominance in these two scenarios. VI. ACKNOWLEDGMENTS S. G. and S. R. acknowledge the contract ANR-- CEXC-- for support. The motivation for this work arose from discussions in two workshops held at the Korea Institute for Advanced Study (KIAS), Seoul, and at the Centre Blaise Pascal, ENS-Lyon in July and August, 22, respectively. S. G. and S. R. thank KIAS for hospitality during their mutual visit in July, 22. A. C. acknowledges ENS-Lyon for hospitality during some stages of the work presented in this paper. We acknowledge useful discussions with David Mukamel and Hyunggyu Park. [] A. Campa, T. Dauois, and S. Ruffo, Phys. Rep. 48, 57 (29). [2] Long-Range Interacting Systems, edited by T. Dauois, S. Ruffo, and L. F. Cugliandolo (Oford University Press, New York, 2). [3] F. Bouchet, S. Gupta, and D. Mukamel, Physica A 389, 4389 (2). [4] P. H. Chavanis, Int. J. Mod. Phys. B 2, 33 (26). [5] D. F. Escande, in Long-Range Interacting Systems, edited by T. Dauois, S. Ruffo, and L. F. Cugliandolo (Oford University Press, New York, 2). [6] F. Bouchet and A. Venaille, Phys. Rep. 55, 227 (22). [7] S. T. Bramwell, in Long-Range Interacting Systems, edited by T. Dauois, S. Ruffo, and L. F. Cugliandolo (Oford University Press, New York, 2). [8] J. Barré, D. Mukamel, and S. Ruffo, Phys. Rev. Lett. 87, 36 (2). [9] D. Mukamel, S. Ruffo, and N. Schreiber, Phys. Rev. Lett. 95, 2464 (25); F. Bouchet, T. Dauois, D. Mukamel, and S. Ruffo, Phys. Rev. E 77, 25 (28). [] Y. Y. Yamaguchi, J. Barré, F. Bouchet, T. Dauois, and S. Ruffo, Physica A 337, 36 (24). [] M. Antoni and S. Ruffo, Phys. Rev. E 52, 236 (995). [2] S. Gupta and D. Mukamel, J. Stat. Mech. Theory Ep. P35 (2). [3] Y. Levin, R. Pakter, and T. N. Teles, Phys. Rev. Lett., 464 (28). [4] M. Joyce and T. Worrakitpoonpon, Phys. Rev. E 84, 39 (2). [5] T. N. Teles, Y. Levin, R. Pakter, and F. B. Rizzato, J. Stat. Mech. Theory Ep. P57 (2). [6] J. Binney and S. Tremaine, Galactic Dynamics (Princeton University Press, Princeton, 28), Chapter 2. [7] C. Anteneodo and C. Tsallis, Phys. Rev. Lett. 8, 533 (998); F. Tamarit and C. Anteneodo, Phys. Rev. Lett. 84, 28 (2). [8] A. Campa, A. Giansanti, and D. Moroni, Phys. Rev. E 62, 33 (2); A. Campa, A. Giansanti, and D. Moroni, J. Phys. A 36, 6897 (23). [9] M.-C. Firpo and S. Ruffo, J. Phys. A 34 L5 (2). [2] C. AnteneodoandR.O.Vallejos, Phys.Rev.E65 62 (2). [2] A. Campa, A. Giansanti, and D. Moroni, Physica A (22). [22] S. Goto and Y. Y. Yamaguchi, Physica A 354, 32 (25). [23] J.A. Acebrón and R. Spigler, Phys. Rev. Lett. 8, 2229 (998); J.A. Acebrón, L.L. Bonilla, and R. Spigler, Phys. Rev. E 62, 3437 (2). [24] P. H. Chavanis, J. Vatteville, and F. Bouchet, Eur. Phys. J. B 46, 6 (25). [25] R. Bachelard, T. Dauois, G. De Ninno, S. Ruffo, and F. Staniscia, Phys. Rev. E 83, 632 (2). [26] S. Gupta, M. G. Potters, and S. Ruffo, Phys. Rev. E 85, 662 (22). [27] Y. Kuramoto, Chemical Oscillations, Waves and Turbulence (Springer, Berlin, 984). [28] S. H. Strogatz, Physica D 43, (2). [29] J. A. Acebrón, L. L. Bonilla, C. J. Pérez Vicente, F. Ritort, and R. Spigler, Rev. Mod. Phys. 77, 37 (25). [3] H. Sakaguchi, Prog. Theor. Phys. 79, 39 (988). [3] J. L. Rogers and L. T. Wille, Phys. Rev. E 54, R293 (996). [32] D. Chowdhuryand M. C. Cross, Phys. Rev. E 82, 625 (2). [33] N. Uchida, Phys. Rev. Lett. 6, 64 (2).

9 9. (a) α=.25,t=... r (t) r (t) r 2 (t) r 3 (t) (b) α=.5,t=.5. r (t) r (t) r 2 (t) r 3 (t) (c) α=.75,t=. r (t) r (t) r 2 (t) r 3 (t) FIG.. (Color online) Time evolution of the observables r (t),r (t),r 2(t), and r 3(t) starting from an initial state {θ i();i =,2,...,N}, where the θ i s are chosen uniformly in [ π, π]. Thus, initially, the system is in the uniform state (6). Here, the values of α and temperature T are (a) α =.25,T =., (b) α =.5,T =.5, and (c) α =.75,T =.. From Table (I), it follows then that the Fourier modes,, 2, 3 in particular are linearly unstable. Consequently, r (t),r (t),r 2(t), andr 3(t)all grow intime from their initial values. However, the plots show that in the long-time limit, r (t) attains a value very close to unity, while r (t),r 2(t),r 3(t) all decay to a value very close to zero. The data in the plots are obtained from numerical simulations with N = 2 4, and involve averaging over independent initial conditions and dynamical realizations. The straight lines show the initial eponential growth with rates given by (T c,m T), where the values of T c,m s can be read off from Table I. The agreement of the growth rates between theory and simulations is very good.

10 . (a) α=.5,t=.5 r (t). N=2 4 N= (b) α=.5,t=.5 N /2 r (t) N=2 4 N= (c) α=.5,t=.5 N /2 r 2 (t) N=2 4 N= (d) α=.5,t=.5 N /2 r 3 (t) N=2 4 N= FIG. 2. (Color online) Time evolution of the observables r (t),r (t),r 2(t), and r 3(t) starting from an initial state {θ i();i =,2,...,N}, where the θ i s are chosen uniformly in [ π, π]. Here, α =.5 and T =.5. The plots show that in the long-time limit, r (t) attains a value very close to unity and does not scale with the system size N, while r (t),r 2(t),r 3(t) all decay to values that scale with N as / N. The data involve averaging over independent initial conditions and dynamical realizations for N = 2 4, and over 5 realizations for N = 2 8.

11 .8 λ m (α)=..6 m.4.2 m =I (m /T)/I (m /T) [-T-T/λ m (α)] / Temperature T FIG. 3. (Color online) Plot showing as a function of T the values of m that satisfy Eq. (35), and those that satisfy (49) at λ m(α) =. The fact that the curves do not intersect at any T in the range < T < /2 shows that in this range, the value of m that solves Eq. (35) at a given T does not satisfy Eq. (49).

12 2 3 2 (a) n ma = n λ m (α)=.2, T=.2, -system ma =2 n ma =5 Im(µ) - Min Re(µ) n ma Re(µ) 3 2 (b) n ma = n λ m (α)=.2, T=.2, y-system ma =2 n ma =5 Im(µ) - Min Re(µ) n ma Re(µ) 3 2 (c) n ma = n λ m (α)=, T=.2, -system ma =2 n ma =5 Im(µ) - Min Re(µ) n ma Re(µ) 3 2 (d) n ma = n λ m (α)=, T=.2, y-system ma =2 n ma =5 Im(µ) - Min Re(µ) 3.4e-7 3.9e e-7 n ma Re(µ) FIG. 4. (Color online) Real and imaginary parts of the eigenvalues µ of the -system, Eq. (53), and the y-system, Eq. (54), for a m mode such that λ m(α) < and for the m = mode such that λ m(α) =. The temperature is T =.2. The parameter n ma denotes the order of truncation of the eigenvalue equations.

13 3 3 2 (a) n ma = n λ m (α)=.2, T=.5, -system ma =2 n ma =5 Im(µ) - Min Re(µ) n ma Re(µ) 3 2 (b) λ m (α)=, T=.5, -system n ma = n ma =2 n ma =5 Im(µ) Re(µ) FIG. 5. (Color online) Real and imaginary parts of the eact eigenvalues µ of the -system, Eqs. (55) and (57), for a m mode such that λ m(α) < and for the m = mode such that λ m(α) =. The temperature is T =.5. The parameter n ma denotes the order of truncation of the eigenvalue equations. r st Temperature Simulation Theory N=24,α=.5 FIG. 6. Plot showing r st as a function of temperature. The data are obtained by starting from the uniform state (6) at a high temperature T =, and then changing adiabatically the temperature as a function of time t as T(t) = Rt, where the rate R is taken to be such that Rτ st, with τ st being the time to reach stationary state at a fied temperature. The condition Rτ st ensures that between two successive values of T, the system gets enough time to attain stationarity. For a system of size N, taking τ st = N 2 dt, where dt is the numerical integration timestep, we take R =./N 2 dt.