arxiv:solv-int/ v1 16 Apr 1997

Transcription

1 SINGULARITY ANALYSIS TOWARDS NONINTEGRABILITY OF NONHOMOGENEOUS NONLINEAR LATTICES arxiv:solv-int/ v1 16 Apr 1997 KEN UMENO Frontier Research Program The Institute of Physical and Chemical Research (RIKEN) -1 Hirosawa, Wako, Saitama , Japan Abstract. We show non-integrability of the nonlinear lattice of Fermi-Pasta-Ulam type via singularity analysis of normal variational equations of Lamé type. 1. From a Nonlinear Lattice to Lamé Equations We consider the following one-dimensional lattice: H = 1 n p i + 1 n+1 υ(q i 1 q i ), (1) i=1 i=1 υ(x) = µ X + µ 4 4 X4 + + µ m m Xm. () Fermi-Pasta-Ulam(FPU) lattice [] is a special type of the systems with the potential function () as follows: H FPU = 1 n i=1 p i + µ n+1 i=1 (q i 1 q i ) + µ 4 4 n+1 i=1 (q i 1 q i ) 4. (3) If we impose the fixed boundary condition as it is easy to check that q 0 = q n+1 = 0, n = odd, (4) Γ : q 1 = Cφ(t),q = 0,q 3 = Cφ(t),,q n 1 = 0,q n = ( 1) n 1 Cφ(t) (5) is a special solution. Thus, the equation of φ(t) is equivalent to the following Hamiltonian system with one degree of freedom: φ+µ φ+µ 4 C φ 3 + +µ m C m φ m 1 = 0, (6)

2 Hamiltonian is H(φ, φ) = 1 ( φ) +µ φ + µ 4C Then the total energy ǫ is given by ǫ = H = H(φ, φ) n+1 φ4 + + µ mc m m φ m = Const. (7) C = n+1 C (µ + 1 µ 4C + 1 m µ mc m ) (8) for the initial condition (5). In the case of the FPU lattice, we can determine C as follows: µ C = + 4ǫ n+1 µ 4 µ. (9) µ 4 By combining (7) with (8), the underlying equation of φ(t) can be rewritten by the differential equation of φ(t) as 1 ( φ) = γ (1 φ )+ γ 4 (1 φ4 )+ + γ m m (1 φm ), (10) γ m (ǫ,{µ j j = 1,,m}) µ m C m. (11) In the case of the FPU lattices (3), the solution of this differential equation (10) with the condition γ m=4 0 (1) is given explicitly by the elliptic function φ(t) = cn(k;αt), (13) α = γ4 γ +γ 4, k =, (14) γ +γ 4 cn(k;αt) is the Jacobi cn elliptic function, and k is the modulus of the elliptic integral. We remark that because γ +γ 4 = µ +C µ 4 = µ + 4ǫ n+1 µ 4 > 0, (15) holds for µ 4 > 0,µ 0, the modulus of the elliptic function k satisfies the following relation: 0 k 1. (16) Thus,thespecial solutions of thefpu lattices forµ 4 > 0,µ 0 have thetwo fundamental periods in the complex time plane as follows: T 1 (ǫ,µ) = K(k) α, T (ǫ,µ) = K(k)+iK (k), (17) α

3 3 K(k) and K (k) are the complete elliptic integrals of the first kind: K(k) = 1 0 dv 1 (1 v )(1 k v ), dv K (k) = 0 (1 v )(1 (1 k )v ). (18) Poles are located at t = τ, τ = K(k) α α (mod T 1,T ) in the parallelogram of each period cell. Let us consider the variational equations along these special solutions. The variational equations are obtained by +i K (k) η j = ξ j = n V Γ k=1 q k q j ξ k = (γ +3γ 4 φ +5γ 6 φ 4 + +(m 1)γ m φ m )(ξ j ξ j 1 ξ j+1 ) for 1 j n, (19) ξ 0 = ξ n+1 = η 0 = η n+1 = 0 and ξ j = δq j,η j = δp j (1 j n). Moreover, these linear variational equations in the form of the vector d dt ξ = (γ + +(m 1)γ m φ m ) ξ (0) can be decoupled as follows. After we note that the eigenvalues of the n n symmetric matrix G = (1) are obtained as {4sin ( jπ (n+1) ) 1 j n} by a normal orthogonal transformation G OGO 1, the variational equations (19) are rewritten in the decoupled form: ξ j (t) = jπ 4sin ( (n+1) )(γ +3γ 4 φ + +(m 1)γ m φ m )ξ j(t) (1 j n), () ξ = Oξ. Clearly, these equations are written in the form of vector Hill s equation[3] d ξ dt +A(t)ξ = 0, A(t+T) = A(t), (3) T = T 1,T in the case of m =. For j = n+1, we have the relation ξ n+1 Thus, the corresponding variational equation = n+1 (ξ 1 ξ 3 +ξ 5 + +( 1) n 1 ξn ). (4) ξ n+1 = (γ +3γ 4 φ + +(m 1)γ m φ m )ξ n+1(t) (5)

4 4 has a time-dependent integral I(ξ, ξ;t) I(ξ,η;t) because I(ξ,η;t) = DH (η p +ξ q )H = η p+ξ V q = C φ(η 1 η 3 +η 5 + +( 1) n 1 η n ) +(Cγ φ+cγ 4 φ 3 + +Cγ m φ m 1 )(ξ 1 ξ 3 +ξ 5 + +( 1) n 1 ξ n ), (6) 1 di C dt = φ( ξ 1 ξ 3 + +( 1) n 1 ξn ) + φ(γ +3γ 4 φ + +(m 1)γ m φ m )(ξ 1 ξ 3 + +( 1) n 1 ξ n ) = 0. (7) We call Eq. (5) the tangential variational equation. On the other hands, a (n )- dimensional normal variational equation(nve) is given by the equation of () with the tangential variational equation (5) removed as follows: η j = 4sin jπ ( (n+1) )(γ +3γ 4 φ + +(m 1)γ m φ m )ξ j, ξ j = η j for 1 j( n+1 ) n. (8) In case of the FPU lattice, the normal variational equation(8) becomes the Lamé equation [7] d y dt (E 1sn (k;αt)+e )y = 0, (9) jπ E 1 = 1 1 α k sin ( (n+1) ) and E are constants.. Non-integrability Theorem Morales and Simó obtained the following theorem on the non-integrability based on the application of Picard-Vessiot theory to Ziglin s analysis[9, 10] for Hamiltonian systems with two degrees of freedom. Theorem 1 (Morales and Simó [4],1994) When the normal reduced variational equation is of Lamé type, if A E 1 α k m(m+1),m N and the Lamé equation satisfying this condition on A is not algebraically solvable(brioschi-halphen-crawford and Baldassarri solutions), then the the initial Hamiltonian system does not have a first integral, meromorphic in a connected neighborhood of the integral curve Γ, which is functionally independent together with H. In case of the present analysis, A is given by the following formula: A = E 1 α k = 1sin jπ jπ ( ) = 6(1 cos( )). (30) (n+1) n+1 We can easily check that cos( jπ n+1 ) / Q if and only if j / {n+1 3, n+1 3 }. When A / Q, the above condition on the algebraic solvability of the Lamé equation is not satisfied. Thus, to check the algebraic solvability of the Lamé equations, (n+1) d ξ j 1 jπ ( dt α k sin ( (n+1) )sn (k;αt)+e )ξ j = 0 (j n+1 ), (31)

5 5 it is sufficient to examine the following two cases: A = 6(1 cos( 1π 3 )) = 3, A = 6(1 cos(π 3 )) = 9. (3) It is known [1] that the condition on A for the Brioschi-Halphen-Crawford solutions is given by A = m(m+1), m+ 1 N, (33) and that the condition on A for the Baldassarri solutions is given by However, the following relations m(m+1) = 3 m = 1± 13 A = m(m+1), m Z 1 4 Z 1 Z \Z. (34) 5 / Q, m(m+1) = 9 m = 1± 37 / Q (35) hold, which guarantee that all n 1 Lamé equations (31) do not belong to the solvable case. In case of the systems with n degrees of freedom, we have n 1 Lamé equations which corresponds to n 1 normal variational equations. Thus, according to the steps in Ref. [4] we obtain the following theorem: Theorem The FPU lattice for µ 4 > 0,µ 0 does not have n 1 first integrals, meromorphic in a connected neighbourhood of the integral curve Γ, which are functionally independent together with H. We remark here that this theorem on the non-integrability does not depend on the total energy in contrast with the result about the non-integrability proof of the FPU lattice in the low energy limit [6] based on non-resonance checking [5] and the result about the nonintegrability of the FPU lattice in the high energy limit based on the Kowalevski exponents of the homogeneous systems [8]. Here, it is conjectured that more general nonhomogeneous nonlinear lattice (1) would be also non-integrable in the sense of the present analysis. Acknowledgements ThepresentauthorwouldliketothankProf.S.Dovysh, Prof.J.J.Morales, Prof.J.Moser, and Prof. H. Yoshida for valuable discussions. He is grateful to Prof. C. Simó, Prof. A. Delshams, Prof. de la Llave and Prof. T. Konishi for their kind hospitalitiy during the 3DHAM95 meeting. He thanks Prof. S. Amari for his continual encouragement. He appreciates support from the Special Researchers Program of Basic Science at the Institute of Physical and Chemical Research (RIKEN) and from the Program of the Complex Systems Workshop at the International Institute for Advanced Study (IIAS). References 1. F. Baldassarri, On algebraic solutions of Lamé differential equation, J. Diff. Eq.41(1981) E. Fermi, J. Pasta and S. Ulam, Los Alamos Report LA-1940(1955); in: Analogies between Analogies: the mathematical reports of S.M. Ulam and his Los Alamos collaborators, Edited by A.R. Bednarek and F. Ulam,(Univ. of California Press, 1990). 3. W. Magnus and S. Winkler, Hill s Equation, (Dover, New York, 1979). 4. J.J. Morales and C. Simó, Picard-Vessiot theory and Ziglin s theorem, J. Diff. Eq. 104(1994)

6 6 5. K. Umeno,Non-integrability proof in some Hamiltonian dynamical systems with three or more degrees of freedom, doctoral dissertation, (University of Tokyo, 1994, December). 6. K. Umeno, Non-perturbative non-integrability of non-homogeneous nonlinear lattices induced by non-resonance hypothesis, Physica D94(1996) E. T. Whittaker and G. N. Watson, A Course of Modern Analysis, (Cambridge Univ. Press, London/New York, 1969). 8. H. Yoshida, A criterion for the non-existence of an additional analytic integral in Hamiltonian systems with n degrees of freedom, Phys. Lett. A 141(1989): S. L. Ziglin, Branching of solutions and non-existence of first integrals in Hamiltonian mechanics. I., Funct. Anal. Appl. 16(1983): S. L. Ziglin, Branching of solutions and non-existence of first integrals in Hamiltonian mechanics. II. 17(1983):6 17.