Low dimensional completely resonant tori in Hamiltonian Lattices and a Theorem of Poincar´e †

† This contribution is part of the Special Issue: Modern methods in Hamiltonian perturbation theory Abstract: We present an extension of a classical result of Poincar´e (1892) about continuation of periodic orbits and breaking of completely resonant tori in a class of nearly integrable Hamiltonian systems, which covers most Hamiltonian Lattice models. The result is based on the ﬁxed point method of the period map and exploits a standard perturbation expansion of the solution with respect to a small parameter. Two di ﬀ erent statements are given, about existence and linear stability: a ﬁrst one, in the so called non-degenerate case, and a second one, in the completely degenerate case. A pair of examples inspired to the existence of localized solutions in the discrete NLS lattice is provided.

Abstract: We present an extension of a classical result of Poincaré (1892) about continuation of periodic orbits and breaking of completely resonant tori in a class of nearly integrable Hamiltonian systems, which covers most Hamiltonian Lattice models. The result is based on the fixed point method of the period map and exploits a standard perturbation expansion of the solution with respect to a small parameter. Two different statements are given, about existence and linear stability: a first one, in the so called non-degenerate case, and a second one, in the completely degenerate case. A pair of examples inspired to the existence of localized solutions in the discrete NLS lattice is provided.

Foreword
The present paper is dedicated to Antonio Giorgilli, in occasion of his 70 th birthday. There is no need to stress here his scientific merits: his publications are worth more than a thousand words. And for sure he would deserve a much better paper than the present contribution, to celebrate his career.
We thus prefer to praise the human qualities of Antonio, which are perceived immediately even by those who know him since only a short period of time; and they get absolutely clear for those who had

Introduction
The present paper deals with an extension of a classical result on periodic orbits in nearly integrable Hamiltonian Systems, due to Poincaré at the end of XIX century [28,29]. The problem considered is the continuation of periodic orbits foliating a completely resonant torus I = I * of maximal dimension m of an integrable Hamiltonian H 0 (I) of m degrees of freedom, once a small perturbation H 1 (θ, I) is added to the system. Given I = I * and θ (0) (t) = ωt + θ(0) the corresponding unperturbed periodic flow (where ω := D I H 0 (I * )), and assuming the invertibility of D where · T and 0 in the rest of the paper will denote respectively the time-average over one period and the evaluation at the unperturbed dynamics on the torus I * . Since any unperturbed periodic orbit on I = I * is uniquely identified by the quotient of the resonant torus with respect to the periodic flow, it turns out that the functional H 1 0 T is defined on the quotient manifold T m−1 and can be expressed in terms of m − 1 suitable resonant angles ϕ (see for example formula (10) in [19]), namely Hence, in terms of these "phase-shift" variables ϕ, the values ϕ * fulfilling ≡ constant is included in [20]; with the same original approach of Poincaré, based on the Taylor-expansion of the solutions, the authors are able to replace the average H 1 0 T with a O( 2 ) higher order functional F on the quotient torus T m−1 , whose non-degenerate critical points ϕ * are the new candidate for the continuation. A further and complete generalization of Poincaré's result is instead given by [26], where a KAM-like normal form approach, combined with fixed point methods, has been used to explore any kind of Poincaré degeneracy. Among the various degenerate scenarios, the normal form approach allows to investigate the case of isolated critical points ϕ * whose hessian in (1.2) is not invertible (hence non-degeneracy does not hold). In this case the existence of the continuation depends on the higher order corrections, which are needed to be explicitly calculated.
One of the typical mechanism for the complete degeneracy to occur, is the lack of any terms, in the given leading order perturbation H 1 , related to all the resonances of the torus; indeed in these cases H 1 0 T turns out to be a constant function on the torus. This kind of degeneracy is naturally met in the study of spatially localized solutions in Hamiltonian Lattices (see for example the Special Issue [21] for a recent collection of results on this topic), like Multibreathers in the Klein-Gordon model or multi-pulse discrete solitons in standard discrete NLS model where J is a finite * set of indexes. Although quite trivial, let us remark that if = 0 these models reduce to a collection of uncoupled identical nonlinear oscillators. One typically considers an unperturbed ( = 0) periodic solution given by a subset S = { j 1 , . . . , j m } of m of these uncoupled identical oscillators, having the same action I j l = I * and hence the same nonlinear frequency ω; in other words, one is considering a low-dimensional completely resonant torus of the unperturbed system, corresponding to the lowest order resonance (1 : 1 : . . . : 1). Once the perturbation H 1 is added (in the above models it is given by the weak linear interactions among the oscillators), the torus breaks down and only some, typically a finite number, of the unperturbed periodic orbits survive at small : these solutions are still spatially localized (in terms of amplitude of the oscillations) around the unperturbed oscillators j l . The mathematical investigation of these solutions naturally represents an extension of the original Poincaré's result to low-dimensional completely resonant tori. The literature on the topic provides two suitable methods to prove existence of such solutions, both based on variational arguments and leading to the study of critical points on the average perturbation H 1 T as in Poincaré. The first method is formulated and developed in [2,3,6,11,14,15,17] and applies to models like (1.3), while the second exploits the rotation symmetry in (1.4) and is developed and applied in [1,10,12,13]. All these results deal with the non-degenerate case, which occurs for example in (1.
. Some results treating partial or total degeneracy in this context might be found in [24,25,27]: the first two works are based on the Lyapunov-Schmidt decomposition, while the last one represents the extension of [26] to the low-dimensional case, hence the perturbation scheme is performed at the level of Hamiltonian normal forms. The results we here present are the extensions to low dimensional tori of both the original nondegenerate Theorem, due to Poincaré, and of the totally degenerate one of Meletlidou-Stagika [20]; such extensions are performed in the special case when the unperturbed Hamiltonian H 0 is already decoupled into two non-interacting integrable subsystems H 0 (I, η, ζ) = H 0 (I) + H 0 (η, ζ) , (1.5) as in the Hamiltonian Lattice models (1.3) and (1.4). From a technical point of view, the continuation is obtained looking for fixed points of the period map, and the expansions in are performed at the level of the (analytic) solutions; the above restriction (1.5) allows to decompose the differential of the period map into a block triangular form, so that the Implicit Function Theorem can be easily applied. The same decomposition also allows to successfully investigate the linear stability of the periodic orbits thus obtained. For integrable Hamiltonian more general than (1.5) and for extensions to partial or higher order degeneracies, normal form techniques are needed, to put the system into a suitable form to investigate the continuation (see [4,5,8,9,16,27,30,31] where (θ, I) ∈ T m × U(I * ) while z = (η, ζ) ∈ C 2n . Assume z = 0 to be an elliptic equilibrium for H 0 (z), with

7)
and I * ∈ U to be a value of the actions which identifies a completely resonant torus of maximal dimension for H 0 (I), with frequency vector ω = ωk. Let ϕ * be a critical point of the average . If the following three assumptions hold true: then, there exists * such that, for | | < * , there exists a periodic orbit with initial datum (ϕ po ( ), I po ( ), z po ( )) analytic in and O( )-close to (ϕ * , I * , 0). Moreover, if the second Melnikov non-resonance condition holds true then the spectrum Σ of the monodromy matrix splits into two different subsets: Σ 1 given by multipliers which have characteristic exponents µ j ( ) becoming zero in the limit → 0, which have the and Σ 2 given by multipliers having purely imaginary characteristic exponents ±2πi In the sufficient assumptions needed for the existence statement above we recover the so called Kolmogorov-non-degeneracy (K-ND), which implies that the resonant torus I * is isolated in U, and the so called Poincaré-non-degeneracy (P-ND), which similarly implies that the critical points ϕ * on the torus T m−1 , selecting the periodic orbits, are isolated. Moreover, we need a third assumption in order to split the tori variables (θ 1 , ϕ, I) from the transversal ones: this is easily obtained by assuming the so called First Melnikov condition (M1), which is a non-resonance assumption between the frequency of the unperturbed periodic orbit ω and the transversal (small) oscillations of frequencies Ω j . The linear stability statement is instead a consequence of the continuity of the spectrum with respect to , which is enough in the easy case of distinct eigenvalues, and of the Krein signature theory, which here applies if the frequencies Ω j all have the same sign. In particular, it is evident that if D 2 I H 0 (I * ) is either positive or negative definite, the relevant role in the stability is played by the nature of ϕ * as critical . The second Melnikov condition (M2) is instead needed (see for example [18]) to preserve transversal ellipticity of the periodic orbits, namely that multipliers of Σ 2 do not leave the unitary circle, via Krein signature: indeed, combining (M2) with the assumption (1.7) on the (positive) sign of Ω j , one can derive the (positive) definite signature for all the unperturbed Floquet multipliers e 2πi Ω j ω on the unitary circle. Linear stability then depends on the effect of the O( ) perturbation on the approximate internal characteristic exponents σ j √ ; for example, in the generic case of distinct σ j , any sufficiently small perturbation would still give distinct µ j ( ) and linear stability is then encoded in then there exists * such that, for | | < * , there exists a periodic orbit with initial datum (ϕ po ( ), I po ( ), z po ( )) analytic in and O( )-close to (ϕ * , I * , 0). Moreover, under the same assumption (M2), the same splitting of Σ as in Theorem 1.1 holds true, with the only difference that µ j ( ) ∼ σ j , with σ 2 j eigenvalues of −T 2 D 2 θ F 2 (ϕ * )D 2 I H 0 (I * ). The precise definition of F 2 in the above statement requires the expansion of the solution up to order O( ), hence it is deferred to formula (2.19) in Section 2, where we prove the two Theorems.
The two statements claimed above are then applied to show existence of discrete solitons in dNLS models, by studying first a typical non-degenerate vortex-like configuration for the ZigZag model, and then showing the existence of only in/out-of-phase discrete solitons in the easiest totally degenerate case (S = −1, 1) for the standard model. This part, which reproduces results already existing in the literature (see for example [22,23]), is developed in Section 3.

Proof
The present section includes the proofs of the two Theorems. A first part shows the formal perturbation scheme here used (which is the same of [20]); the second and third parts provide the continuation argument in the non-degenerate and totally degenerate case; the last part is dedicated to linear stability, hence giving the splitting of the spectrum and asymptotic behaviour of Floquet exponents in the limit → 0.

Formal expansions and leading order approximations
Due to the near integrability of the model, it is known that any T -periodic solution of (1.6) (θ, I, z)(t; ) is analytic in ; we can thus Taylor-expand all the variables in and write . According to the same splitting of the integrable part H 0 given in (1.5), we can decompose the perturbation H 1 as where H 1 depends only on the action-angle variables, H 1 depends only on the complex variables z = (η, ζ) while H 1 properly depends on all the variables, and provide the interaction among the two subsystems. Notice that the term H 1 has necessarily to be at least linear in z, so that it has to vanish when evaluated at the unperturbed flow, namely If we expand w.r.t. both the Hamiltonian vector field and the time derivatives we get the two systems: and these necessarily have to coincide at any order in . Equating terms of order 0 we get We deduce immediately that actions keep their unperturbed initial values I (0) (0), the angles rotate on the torus T m with frequencies ∇ I H 0 (I (0) ), while the external variables z evolve according to their unperturbed nonlinear dynamics ; by inserting order 0 approximation of the periodic solutions on the unperturbed torus I * where equation forİ (1) andθ (1) are uncoupled fromż (1) . Notice that equation for z (1) simplifies in those models with H 1 (z) = O(z 2 ); in fact in such cases one has ∇ z H 1 (z) = O(z), which implies a vanishing contribute of ∇ z H 1 in the third equation, if at = 0 the variables z (0) stay at rest Solving the equation forİ, we get the leading order corrections to I * of the actions Since I (1) (t) have to be periodic of period T , namely I (1) (T ) = I (1) (0), we have to impose Let us set for the sake of simplicity k 1 = 1 in the resonant frequency vector ωk, with ω := 2π T . We introduce the m − 1 resonant angles ϕ = {ϕ j } m−1 j=1 = k j θ 1 − θ j , which are the natural coordinates of the quotient torus T m−1 of the resonant torus over the periodic flow. Hence we get (by setting θ (0) Indeed, in the previous average, we can exchange the time-variable with θ 1 , so that the first element of the gradient ∇ θ H 1 has zero average and in the m − 1 remaining components, only the dependence on the angles ϕ is left in the average Thus, the periodicity condition (2.7) easily reads solutions ϕ * of (2.9) represent critical points of the functional (2.8) defined on the torus T m−1 .
Le us define the Hessian matrix C := D 2 thus the equation forθ becomes since also the angles θ (1) have to be T -periodic (uniformly in ), we have to ask where the first term is the average of I (1) (t) while the second represents the average of ∇ I H 1 restricted (as usual) to the unperturbed flow. Recalling (2.6) one has It is clear from (2.11) that I (1) (0) depends on the unperturbed initial data I * and ϕ 0 . Let us now move to the equation forż In the Taylor expansions of H 1 and H 1 with respect to the external variable z, the important terms are linear in z; indeed all those terms which are O(z 2 ) produce contributions, in the equation, that vanish when evaluated at z (0) = 0. Hence represents a periodically forcing term in the inhomogeneous linear equation (2.12) if we assume the non-resonance (M1) condition between Ω j and the frequency ω, then h 1 is nonresonant with respect to JD and the solution of (2.12) is given by the convolution integral The T -periodicity of z (1) provides, as for I (1) , the correction at order O( ) of the initial value z 0 = z (1) (0) + O( 2 ); also in this case, z (1) (0) has to depend on I * and ϕ 0 . Indeed, from one gets the correction with I − e JDT invertible because of (M1).

Continuation in the non-degenerate case
Given a generic initial datum P 0 of the unperturbed system ( = 0), we look for a correction P( ) at 0 such that the variation of the Hamiltonian flow after one period T vanishes, namely Φ T (P( ), )− P( ) = 0. The condition we have to impose is then the T -periodicity of the flow, assuming its validity at = 0, and exploiting the analyticity of the flow Φ T (P( ), ), and of its initial datum P( ), with respect to .
In order to identify the periodic orbit, we can ignore its phase θ 1 and consider as unknowns only the other initial 2n + 2m − 1 "transversal" variables {ϕ 0 ( ), I 0 ( ), z 0 ( )}. Moreover, since the Hamiltonian is kept constant along the orbit, the 2(n + m) equation Φ T (P( ), ) − P( ) = 0 are not independent, and we can get rid of one of them; we choose to ignore the action I 1 associated to the angle θ 1 . In this way, the periodicity condition reduces to a set of 2n + 2m − 1 equations in 2n + 2m − 1 variables of the form X(ϕ 0 , I 0 , z 0 , ) = 0 . (2.13) Since the velocity of variation of the actions I l=2,...,m (which are approximately constant in time) are of order O( ), we can rescale them by ; in this way the leading terms of (2.13), with respect to , read , where now {ϕ 0 , I 0 , z 0 } are the variables which identify the generic unperturbed initial datum P 0 , modulo a shift of θ 1 , and the components F j read where H 1 (τ) is the perturbation restricted to the unperturbed flow corresponding to the initial datum (ϕ 0 , I 0 , z 0 ). By redefining X coherently with the previous scaling of the actions J, we get the compact form (2.14) Solving the above equation (2.14) at = 0 is equivalent to solve the system F (ϕ 0 , I 0 , z 0 ) = 0. Hence, if I 0 = I * corresponds to a completely resonant torus with frequency vector ωk, then F 1 (I * ) = 0; given such I * and setting at rest the complex variables z 0 = 0, if ϕ 0 = ϕ * is a critical point of then also F 2 (ϕ * , I * ) = 0. Thus, the point (ϕ * , I * , 0) solves In order to apply the Implicit Function Theorem, we need invertibility of X , which is the differential of X evaluated on the approximate solution (ϕ * , I * , 0) at = 0; this coincides † with the differential of F at its zero and is block triangular. Indeed From the first assumption (K-ND) the matrix C is invertible, and the same holds true for due to the second one (P-ND); in particular notice that T is at least linear in z 0 . The invertibility of block e JDT − I of X is due to the first Melnikov condition (M1), which ensures that Ω j ω Z. Due to the triangular structure of X , the block X 12 does not play any role, even if different from zero. † Recall that to compute X we can first restrict to = 0 and then differentiate w.r.t. the unperturbed initial datum z 0 .

Continuation in the totally degenerate case
The totally degenerate case corresponds to H 1 T ≡ const; what follows represents the extension of what has been already developed in [20] in the case of maximal tori.
In the degenerate scenario, the variation of the actions I l on the approximate periodic orbit is much slower than what happens in the non-degenerate case, beingİ l = O( 2 ); as a consequence, an expansion up to order 2 of the actions I l is required. Such an expansion provides the following equation for I (2) Since the term H 1 (θ, I, z) is at least linear in z, then two of the above terms will vanish when evaluated at z = 0 hence we can reduce (2.16) to the simplified forṁ where D 2 zθ H 1 (t)z (1) (t) is given only by linear terms in z, and the first order corrections θ (1) (t), I (1) (t) and z (1) (t) have been determined by the periodicity conditions at order O( ). By integrating on the period interval [0, T ], the periodicity condition I (2) where in the first addendum the initial datum θ (1) (0) can be "added" without affecting the computation of the average; indeed one gets D 2 θ H 1 (τ)θ (1) (0) T ≡ 0, due to the degeneracy assumptions. Moreover, thanks to the conservation of the energy along the periodic orbit, it is possible to prove that only m − 1 of the above m equations are independent. Hence, as in the non-degenerate case, we can ignore the evolution of the first action I 1 and consider the periodicity condition only for the remaining m − 1 ones. Due to the complete degeneracy, the Taylor expansion of the m − 1 actions I l starts with term of order O( 2 ); thus we rescale I l (T ) − I l (0) by 2 and define the components of F 2 as where the initial values I (1) 0 and z (1) 0 explicitly depend on the unperturbed initial values ϕ 0 and I 0 . We want to apply the IFT to the periodicity condition (2.14), where now ϕ * has to solve (2.20) The above Eq (2.20) selects those ϕ 0 = ϕ * on the torus T m−1 which can be continued, and that provide the solution of (2.14) at = 0 X(ϕ * , I * , 0) = 0 .
In order to apply again the IFT we need such ϕ * to be non-degenerate critical points for F 2 (ϕ 0 , I * ), namely det D ϕ F 2 (ϕ * , I * ) 0 , which is indeed hypothesis (1.8).

Linear stability
In order to investigate the linear stability of the periodic orbits, we have to exploit the fact that the monodromy matrix Λ( ) coincides with the differential of the Hamiltonian flow at time T with respect to the initial datum of the periodic orbit. We develop the non-degenerate case, being the totally degenerate one a minor variation of the forthcoming arguments.
With standard calculations, which resemble those for the differential X of the period map, one obtains the following structure for Λ( ) with: • the Λ 11 block reads • the out-of-diagonal blocks are perturbations of the same order O( ) since the interaction between internal variables (θ, I) and the external ones (ξ, η) is only due to the perturbation H 1 ; • the block Λ 22 reads Λ 22 = e JDT + O( ) .
By continuity with respect to , the multipliers of Λ( ) have to converge to those of Λ(0) as → 0, the last ones being collected into two different sets (due to the block diagonal shape of Λ(0)): 2n multipliers of Λ 22 (0) are on the unitary circle, with purely imaginary exponents ±2πi Ω j ω , and 2m multipliers of Λ 11 (0) equal to 1 (due to the complete resonance of the n-dimensional torus). Moreover, due to the two Melnikov conditions ‡ , the 2n multipliers on the unitary circle are different from ±1. By continuity in , the two sets stay disjoint and provide the two different sets Σ 1 and Σ 2 of the statements: Σ 1 is made of multipliers λ( ) which all bifurcate from 1, hence having vanishing characteristic exponents, while Σ 2 is made of multipliers which are deformations of e ±2πi Ω j ω ±1. If some (or all) of the unperturbed Floquet multiplier of Λ 22 (0) coincide, then having the same definite Krein signature (see [7,18,32]) is a sufficient condition for their perturbations not to leave the unitary circle. Due to the Melnikov conditions, and since the frequencies Ω j share the same sign (rotation direction), the unperturbed Floquet multipliers e ±2πi Ω j ω are all positive definite, so any sufficiently small perturbation λ j ( ) cannot leave the unitary circle, still being different from ±1. If all Ω j are ‡ Notice that (M1) implies ± Ω j ω l ∈ Z, hence the multipliers have to be different from 1. The second condition (M2) implies in particular ± Ω j ω l 2 ∈ Z hence the multipliers have to be different also from −1. We recommend reference [18]. distinct, then the same holds true for the corresponding unperturbed Floquet multipliers (which are all simple eigenvalues of Λ( )) and for any sufficiently small perturbations, without invoking their Krein signatures.
It remains to show the prescribed asymptotic behaviour of the Floquet multipliers of Σ 1 . These are solutions of with vanishing characteristic exponents, hence λ j ( ) → 1. We thus assume for these part of the spectrum the ansatz λ j ( ) = e √ σ j +O( ) , where the √ scaling is expected (as in the original Poincaré result, but see also [30]) and will be clarified in a while. Since Here we can apply a second time the Schur complement, being the diagonal blocks invertible for 0 small enough, so that det(M 11 (σ, )) = 0 reads Taking the limit → 0, we obtain the asymptotic expansion claimed in the statement of Theorem 1.1; this concludes the non-degenerate case, since by adding the smaller perturbation term M 12 M −1 22 M 21 the same arguments hold true.
The totally degenerate case can be faced following the same steps above, with a minor but crucial variation which involves the perturbation M 12 M −1 22 M 21 . Indeed, since the magnitude ofİ is of order O( 2 ), it follows that M 12 has the special structure which provides hence it is component wise a perturbation of the leading term M 11 where F 2 is taken with all its m components.

Spatially localized solutions in discrete NLS models
The aim of this section is to exploit Theorem 1.1 to investigate the existence of a special class of spatially localized time-periodic solutions in a dNLS model of the form where the (weak) linear term L is beyond nearest-neighbours the coupling parameter has to be considered small enough (so called anti-continuum limit) and the boundary conditions are of Dirichlet type in the case of J finite. The equations can be written in where Being interested in solutions which are spatially localized on many-sites, we introduce the set of excited sites where we stress that the indexes do not have to be necessarily consecutive; in other words, we are including also configurations where the localization of the amplitude (hence of the energy), is clustered, with "holes" separating the different clusters along the chain. The easiest and most popular spatially localized solutions in the considered dNLS model correspond to the so-called discrete solitons (or also multi-pulse solutions). These are periodic solutions ψ(t) of (3.1) which at = 0 take the form where φ (0) is independent of time. Since the unperturbed excited oscillators {ψ (0) j } j∈J in (3.5) are required to keep the same frequency ω, thus being in 1 : 1 resonance, a common amplitude R > 0 is necessary, so that the unperturbed spatial profile φ (0) reads where ω(R) = 1 + γR 2 . (3.7) All these orbits are uniquely defined except for a phase shift, which corresponds to a change of the initial configuration in the ansatz (3.5). According to this geometrical interpretation, all the unperturbed periodic orbits foliate a completely resonant m-dimensional torus T m of the phase space, and any orbit is uniquely identified by a point in the quotient space T m−1 = T m \ T; such a point can be well represented by introducing a set of m − 1 new "phase shifts" variables In order to put the dNLS model in the form (1.6) of Theorem 1.1, we introduce action-angle variables for the sites j ∈ S and Birkhoff complex variables z = (η, ζ) (replacing the complex variable ψ) for the remaining sites in this way the symplectic form reads dη ∧ dζ + dθ ∧ dI and the integrable part of the Hamiltonian takes the form while the perturbation H 1 can be split in the three parts H 1 , H 1 , H 1 depending on the shape of L in (3.2) and of the choice of S in the considered example. Moreover, once fixed at I * the common action of the excited nonlinear oscillators ψ j , j ∈ S , the frequency (3.7) reads ω(I * ) = 1 + γI * .
We will benefit of rewriting the generic linear interaction term as In this way, a term of H 1 would read one of H 1 would read iζ j+l η j+l + iζ j η j − i ζ j+l η j + η j+l ζ j , and one of H 1 (assuming j ∈ S and j + l ∈ J \ S ) − I j ζ j+l e iθ j − i I j η j+l e −iθ j .

Non-degeneracy and consecutive excited sites
Let us consider, instead of the generic model (3.4), the dNLS Hamiltonian with nearest neighbor and next to nearest neighbor linear interactions, on a lattice J = {−N, . . . , N} and take m N excited consecutive sites in the set S = {1, . . . , m}. In order to find candidate solutions, we need to solve (2.15), i.e., to find critical points of H 1 0 T , which explicitly reads by introducing the coordinates ϕ i = θ (0) i+1 − θ (0) i one gets the (φ-dependent part of the) functional on the torus T m−1 whose critical points provide the solutions ϕ * we were looking for. For examples, let us consider three consecutive sites m = 3 with S = {1, 2, 3}, in the so-called ZigZag model k 1 = k 2 = 1. The average H 1 0 T defined on the torus T 2 is H 1 0 T = cos(ϕ 1 ) + cos(ϕ 2 ) + cos(ϕ 1 + ϕ 2 ) , so that Eq (2.15) reads        sin (ϕ 1 ) + sin (ϕ 1 + ϕ 2 ) = 0 sin (ϕ 2 ) + sin (ϕ 1 + ϕ 2 ) = 0 , which provides six different solutions: four of them correspond to in/out-of-phase configurations ϕ * ∈ which is a real function since the second addendum is the complex conjugate of the first; in the above and in what follows, we have to keep in mind that θ −1 (t) = ωt , θ 1 (t) = ωt + φ . (3.10) We can also write explicitly which tells us that only evolutions of the complex configurations ζ (1) 0 and ζ (1) −2 are needed. Their equations are (remember that ζ are the momenta of the canonical variables) where the different forcing terms depend on the fact that the 0-site is excited by two neighbours, while the −2-site only by one. The solutions, because of the non-resonance condition between the linear frequency Ω = 1 and ω, are then given by where the Cauchy problems are given by the periodicity conditions ζ which provides Hence, being the above term constant in time, the average provides exactly where we have used the explicit dependence of ω with respect to I * ; critical points of F 2 are only ϕ * ∈ {0, π}, which are non-degenerate, and the IFT applies. Linear stability immediately follows from − F 2 (ϕ * ): we stress that here γ does not play any role, and stability is exchanged only by switching from positive to negative values of , hence from attractive to repulsive linear interaction.