arxiv: v1 [math.ds] 30 Sep 2017

Transcription

1 NON-INTEGRABILITY OF THE ARMBRUSTER-GUCKENHEIMER-KIM QUARTIC HAMILTONIAN THROUGH MORALES-RAMIS THEORY P. ACOSTA-HUMÁNEZ,, M. ALVAREZ-RAMíREZ,, AND T. J. STUCHI, arxiv:7.7v [math.ds] 3 Sep 7 Abstract. We show the non-integrability of the three-parameter Armburster-Guckenheimer- Kim quartic Hamiltonian using Morales-Ramis theory, with the exception of the three already known integrable cases. We use Poincaré sections to illustrate the breakdown of regular motion for some parameter values. Key words. Hamiltonian, integrability of dynamical systems, differential Galois theory, Legendre equation, Schrödinger equation. AMS subject classifications. 37J3, 7H6, H5. Introduction. Mechanical and physical systems whose energy is conserved are usually modelled by Hamiltonian systems with n degrees of freedom. These systems have a Hamiltonian function H as a first integral and we say that the Hamiltonian system is integrable in the Liouville sense (often referred to as complete integrability if there exist n additional smooth first integrals in involution {I j, H} =, j =,..., n, which are functionally independent, see []. There are many physics and mathematical papers with examples of numerical solutions of Hamiltonian systems with two or more degrees freedom giving evidence of non-integrability. In the present paper, using Morales-Ramis theory [7] and some new contributions to this, we examine rigorously the non-integrability of the Hamiltonian of two degrees of freedom with Armbruster-Guckenheimer-Kim potential given by ( H(x, y = (p x + p y µ (x + y a 4 (x + y b x y. where µ, a and b are real parameters. In 989 Simonelli and Gollub [3] carried out experimental work with surface waves in a square and rectangular containers subject to vertical oscillations. Their methods allowed measurements of both stable and unstable fixed points (sinks, sources, and saddles, and the nature of the bifurcation sequences clearly established and reported. Later, motived by these experiments reported in [3], Armbruster et. al. [8] derived a dynamical model starting with a normal form given by four first order differential equations depending on several parameters. Their system provides a general description of the codimension two bifurcation problem. To restrict the model they performed a rescaling of the variables and parameters. When one of the scaling parameters ε =, the system of equations becomes Hamiltonian with the Hamiltonian function giving by ( which has D 4 -symmetry. Armbruster et. al [8] found a large parameter region of chaotic behaviour. They argue that the existence of chaos in a given region suggests that the Hamiltonian ( is non-integrable. School of Basic and Biomedical Sciences, Universidad Simón Bolívar, Barranquilla - Colombia (primitivo.acosta@unisimonbolivar.edu.co; primi@intelectual.co Departamento de Matemáticas, UAM-Iztapalapa, San Rafael Atlixco 86, Col. Vicentina, 934 Iztapalapa, México city, México. (mar@xanum.uam.mx Instituto de Física - Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ CP: CEP: 94-97, Brazil. (tstuchi@if.ufrj.br

2 ACOSTA-HUMÁNEZ, ALVAREZ-RAMíREZ, STUCHI As early as 966 a similar potential was proposed by Andrle [7] as a dynamical model for a stellar system with an axis and a plane of symmetry given by V = (C ω A ξ B z + aξ4 + γz4 + βξ z where ξ = ω ω, A, B, α, γ are positive constant, β is a constant and C is the value of the angular momentum integral. This potential has often been used in the study of galactic dynamics, see among others []. In particular, if µ = in ( we have the Hamiltonian studied by Llibre and Roberto [5]. These authors examined its C non-integrability in the sense of Liouville- Arnold, see []. More recently, Elmandouh [] studied the dynamics of rotation of a nearly axisymmetric galaxy rotating with a constant angular velocity around a fixed axis using ( in a rotating reference frame. He also examined the non-integrability by means of Painlevé analysis. The goal of this present paper is to show that the Hamiltonian ( is, in fact, non-integrable except for the cases b =, a = b and b = a mentioned above. The potential of the Hamiltonian ( can be written in the form of the sum of two homogeneous components V (x, y = µ (x + y a 4 (x + y b x y = V min + V max, where V min and V max are quadratic and quartic polynomial, respectively. Maciejewski and Przybylska [6] and Bostan et. al. [9] have made advances on the Morales-Ramis Theorem for homogeneous potentials. In particular [9] has found a way of applying Morales-Ramis Theorem for homogeneous potentials in a very efficient way. For all Darboux points, Bostan et. al. [9] implemented an algorithm that allows to determinate necessary integrability conditions of homogeneous potentials up to degree 9. We use this algorithm to obtain necessary integrability conditions for the quartic potential. Then using an old result of Hietarinta [] and Yoshida [7] for polynomial potentials the necessary integrability conditions for AGK potential is obtained. The above necessary conditions imply integrability by rational functions. The general Morales-Ramis Theorem (see [7] complements the study by given necessary conditions for rational or meromorphic integrability, depending on the kind of singularity of the normal variational equation. The structure of the present paper is as follows. In Section, we briefly review some preliminary notions and results of Morales-Ramis theory to be used later. The main results will be given in Section 3. In Section 4 we use Poincaré sections to exhibit the breakdown in regular motion suggesting non-integrability.. Preliminars. In this section we describe some basic facts concerning the Morales-Ramis theory... Differential Galois Theory. Picard-Vessiot theory is the Galois theory of linear differential equations, also known as differential Galois theory. We present here some of its main definitions and results, and refer the reader to [5] for a wide theoretical background. We start by recalling some basic notions on algebraic groups and, afterwards, Picard-Vessiot theory is introduced. An algebraic group of matrices is a subgroup G GL(, C defined by means of algebraic equations in its matrix elements and in the inverse of its determinant. It is well known that any algebraic group G has a unique connected normal algebraic

3 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE 3 subgroup G of finite index. In particular, the identity connected component G of G is defined as the largest connected algebraic subgroup of G containing the identity. In the case that G = G we say that G is a connected group. Moreover, if G is abelian we say that G is virtually abelian. Now, we briefly introduce Picard-Vessiot Theory. First, we say that (K, - or, simply, K - is a differential field if K is a commutative field of characteristic zero, depending on x and = d/dx is a derivation on K (that is, satisfying (a + b = a + b and (a b = a b + a b for all a, b K. We denote by C the field of constants of K, defined as C = {c K c = }. We will deal with second order linear homogeneous differential equations, that is, equations of the form ( y + b y + b y =, b, b K, and we will be concerned with the algebraic structure of their solutions. Moreover, along this work, we will refer to the current differential field as the smallest one containing the field of coefficients of this differential equation. Let us suppose that y, y is a basis of solutions of equation (, i.e., y, y are linearly independent over C and every solution is a linear combination over C of these two solutions. Let L = K y, y = K(y, y, y, y be the differential extension of K such that C is the field of constants for K and L. In this terms, we say that L, the smallest differential field containing K and {y, y }, is the Picard-Vessiot extension of K for the differential equation (. The group of all the differential automorphisms of L over K that commute with the derivation is called the differential Galois group of L over K and is denoted by DGal(L/K. This means, in particular, that for any σ DGal(L/K, σ(a = (σ(a for all a L and that σ(a = a for all a K. Thus, if {y, y } is a fundamental system of solutions of ( and σ DGal(L/K then {σy, σy } is also a fundamental system. This implies the existence of a non-singular constant matrix ( a b A σ = GL(, C, c d such that σ ( y y = ( σ(y σ(y = ( y y Aσ. This fact can be extended in a natural way to a system ( ( ( y y σ σ(y σ(y y y = y y σ(y σ(y = y y A σ, which leads to a faithful representation DGal(L/K GL(, C and makes possible to consider DGal(L/K isomorphic to a subgroup of GL(, C depending (up to conjugacy on the choice of the fundamental system {y, y }. One of the fundamental results of the Picard-Vessiot theory is the following theorem (see [3, 4]. Theorem. The Galois group DGal(L/K is an algebraic subgroup of GL(, C. We say that Eq. ( is integrable if the Picard-Vessiot extension L K is obtained as a tower of differential fields K = L L L m = L such that L i = L i (η for i =,..., m, where either (i η is algebraic over L i ; that is, η satisfies a polynomial equation with coefficients in L i.

4 4 ACOSTA-HUMÁNEZ, ALVAREZ-RAMíREZ, STUCHI (ii η is primitive over L i ; that is, η L i. (iii η is exponential over L i ; that is, η /η L i. If η is obtained as a combination of (i, (ii and (iii, we say that η is Liouvillan. Usually in terms of Differential Algebra s terminology we say that Eq. ( is integrable if the corresponding Picard-Vessiot extension is Liouvillian. Moreover, the following theorem holds. Theorem (Kolchin. Equation ( is integrable if and only if DGal(L/K is virtually solvable, that is, its identity component (DGal(L/K is solvable. Let us notice that for a second-order linear differential equation y = r(xy with coefficients in K, the only connected non-solvable group is SL(, C, see [3, 7, 5]. The following result can be found in [5]. Theorem 3 (Acosta-Humánez, Morales-Ruiz, Lázaro & Pantazi, [5]. differential equation (3 ( x d y dy x (ν(ν dx dx + + µ x y = is integrable if and only if, either. µ ± ν Z or ν Z, or. ± µ, ±(ν + belong to one of the following seven families Case µ ν µ + ν (a Z + C (b Z ± 3 (c Z ± 5 (d Z ± 3 (e Z ± 5 (f Z ± 5 Z ± 3 Z + 6 Z ± 5 Z + Z ± 5 Z + Z ± 5 Z + Z ± 3 Z + 6 We recall that the singularities of Legendre equation (3 are of the regular type. Legendre Theorem 4 (Acosta-Humánez & Blázquez-Sanz, [3]. Assume K = C(e it. The differential Galois group of extended Mathieu differential equation is isomorphic to SL(, C. ÿ = (a + b sin(t + c cos(ty, b + c Let us remark that the variable change x = e it makes the differential equation algebraic, such that x = is an irregular singularity... Morales Ramis Theory. A Hamiltonian H in C n is called integrable in the sense of Liouville if there exist n independent first integrals of the Hamiltonian system in involution. We say that H is integrable in terms of rational functions if we can find a complete set of integrals within the family of rational functions. Respectively, we can say that H is integrable in terms of meromorphic functions if we can find a complete set of integrals within the family of meromorphic functions. Morales-Ramis theory relates integrability of Hamiltonian systems in the Liouville sense with the integrability of Picard-Vessiot theory in terms of differential Galois theory, see [7, 8, 9, ].

5 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE 5 Let V be a meromorphic homogeneous potential. We say that c C n \ {} is a Darboux point if there exists α C \ {} such that q j V (c = αc j j =,..., n. Theorem 5 (Morales-Ramis, [7, 9, ]. Let H be a Hamiltonian in C n, and γ a particular solution such that the variational equation associated to the system has regular (resp. irregular singularities at the points of γ at infinity. If H is integrable by terms of meromorphic (rational functions, then the differential Galois group of the variational equation is virtually abelian. Theorem 6 (Morales-Ramis [7, 9, ]. Let V C(q, q be a homogeneous rational function of homogeneity degree k,,, and let c C } be a Darboux point of V. If the potential V is integrable, then for any eigenvalue λ of the Hessian matrix of V at c, the pair (k, λ belongs to the Table, for some j Z. k λ k λ Z jk(jk + k ( 5 + 6j Z (jk + (jk + k ( + 6j ( 3 + j ( 3 + 6j (4 + j ( j ( j ( 5 + 6j ( + 6j 4 + ( j ( 3 + 6j ( 3 + j ( j (4 + 6j Table Morales-Ramis table..3. Galoisian integrability of Schödinger equation. The analysis of the integrability of the Schrödinger equation using differential Galois theory has been studied in [4, 6], using as main tools Kovacic algorithm and Hamiltonian algebrization; wave functions and spectrum of Schrödinger equations with potentials known as shape invariant were obtained. It is well known that the one-dimensional stationary non-relativistic Schrödinger equation is giving by (4 Ψ = (V (y λψ, where V is the potential and λ is called energy. A special class of potentials are the called Pöschl-Teller (see [, 4] which transforms (4 into ( (5 Ψ l(l + = cosh y + n Ψ, y R, where n, l Z +. It is possible to algebrize this equation by means of the change of variables z = tanh y. This substitution transforms the Schrödinger equation (5 into (6 ( z d Ψ dz dψ z (l(l dz + + n z Ψ =,

6 6 ACOSTA-HUMÁNEZ, ALVAREZ-RAMíREZ, STUCHI which is a differential equation for the associated Legendre function. Two linearly independent solutions of (6 are given by the associated Legendre functions of the first and second kind, respectively. 3. Main results. We start by giving the known integrable cases of the Hamiltonian ( with associated differential equations ẋ = p x, ẏ = p y, p x = µx + a(x + y x + bxy, p y = µy + a(x + y y + bx y. There exist values of the parameters for which the system has additional symmetry and ( is completely integrable: i If b =, by considering polar coordinates x = r cos θ and y = r sin θ, the Hamiltonian ( becomes (7 H(x, y = p r + p θ r µ r a 4 r4, which does not depend on θ. Since p θ is a second independent first integral the angular momentum r θ = xẏ ẋy is conserved. ii If a = b, system (7 reduces to two identical uncoupled Duffing s oscillators ẍ = µx + ax 3 and ÿ = µy + ay 3, which can be integrated by quadrature, see [6]. iii In case b = a, we employ the symplectic transformation x = (u v,y = (u + v, p x = (p u p v,p y = (p u + p v and one easily finds that the Hamiltonian ( becomes (8 H(u, v = (p u + p v µ (u + v a (u4 + v 4. The transformed Hamiltonian system decouples into two single degree of freedom oscillator ü µu au 3 = and v µv av 3 =, and each of the them is integrable by quadrature. Let us remark that the first two cases were previously discovered by Armbruster et. al [8]. Elmandouh [] also cited these two cases and added the case b = a when the angular velocity ω =. 3.. Application of Bostan-Combot-Safey El Din algorithm. The AGK Hamiltonian has the following polynomial potential which can be grouped into two homogeneous potentials of grades two and four: (9 V (x, y = µ(x + y + a 4 (x + y b x y = V min + V max. An old result of Hietarinta ([], Corollary of Theorem 3 proved again by Yoshida [7] and Nakagawa [] states the following.

7 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE 7 Proposition 7. Suppose that the Hamiltonian with a non-homogeneous polynomial potential, H = (p + p + V min (q, q + + V max (q, q admits an additional rational first integral Φ = Φ(q, q, p, p. Here V min is the lowest degree part of the potential and V max the highest degree part. Then, each of the two Hamiltonian, H min = (p + p + V min (q, q, H max = (p + p + V max (q, q also admits an additional rational first integral. In our case V min = µ(x + y is obviously integrable for any value of µ. It remains to determine for which values of the parameters the quartic potential V max is integrable. Homogeneous potentials have been worked before and after Morales- Ramis seminal theorem ([9], [6], in particular Bostan et. al. [9] take the potential to a polar form by making x = r cos θ, y = r sin θ and reducing them to the form ( V (r, θ = r k F (e iθ. This change makes the evaluation of Morales-Ramis table (Table for the necessary condition of integrability for homogeneous potentials easier. Bostan et al. are able to built a nice Maple algorithm using the following theorem. Theorem 8. (Bostan et. al.[9] Let V C(x, y be a homogeneous potential with polar representation (F, k and let Λ be the following set { Λ(F, k := k z F } (z : z, F (z =, F (z. F (z Let Sp D ( V denote the union of the sets Sp( V (c taken over all Darboux points c D of V, then k(k Sp D ( V = k(k Λ. Moreover, if V is integrable, then Λ E k, where E k is the so called Morales-Ramis table. This result allows the construction of a powerful algorithm which allows quick evaluation of integrability even to homogeneous potential of degree 9. Applying the Bostan-Combot-Safey El Din algorithm given in [9], we have found that the necessary values for integrability are b =, b = a and b = a which are the same of the well known integrable cases described at the beginning of this section. The algorithm also inform us the spectrum of the Hessian matrix for each necessary condition of rational integrability, yields b =, the Darboux points have Hessian spectrums [, 4], b a =, the Darboux points have Hessian spectrums [, ] or [, ], a + b =, the Darboux points have Hessian spectrums [, ] or [, ]. We note that the trivial eigenvalue of the Hessian λ = appears also as a non trivial eigenvalue. Now applying Theorem 7 since the harmonic oscillator is trivially integrable we get the necessary condition for the integrability of the quartic potential. The procedure described at the beginning of this section showing the separability of the potential for the classical values implies the following result.

8 8 ACOSTA-HUMÁNEZ, ALVAREZ-RAMíREZ, STUCHI Theorem 9. The Armburster-Guckenheimer-Kim potential ( is rationally integrable if and only if b =, b = a and b = a. Note that for µ = we have an homogeneous quartic potential whose necessary conditions for rational integrability are the same values given in Theorem Application of Morales-Ramis Theorem. The Hamiltonian has the square symmetry, that is, it is invariant under a rotation by π/, and is also time reversible. So, the planes of symmetry are given by planes which are invariant under the action of the dihedral group D 4 : Γ = {(x, y, p x, p y : y = p y = } Γ = {(x, y, p x, p y : x = p x = } Γ 3 = {(x, y, p x, p y : y = x, p y = p x } Γ 4 = {(x, y, p x, p y : y = x, p y = p x } Γ 5 = {(x, y, p x, p y : y = x, p y = p x } Γ 6 = {(x, y, p x, p y : y = x, p y = p x }. Let us observe that if we restrict Hamiltonian ( to the invariant planes of symmetry Γ, and respectively Γ 3,4,5,6, we get the same form. Thus, it is enough to consider the following two cases: (a h = p x µ x a 4 x4, H Γ = h, (b ĥ = p x µ x a + b x 4, H Γ3 = 4 ĥ. We note that each one of these expressions correspond to a Hamiltonian with one degree of freedom, and also setting change a a + b into (a we obtain (b. So, along the paper we will compute the conditions of integrability of the normal variational equations over Γ and Γ 3, and also we substitute a a + b into the particular solution x(t obtained in Γ to get them in Γ 3. In order to apply the Morales-Ramis theory our strategy consists in selecting a non-equilibrium particular solution in the invariant planes Γ and Γ 3. The following step is to obtain an integral curve for the Hamiltonian (a. Thus, to perform our study we take into account that a Hamiltonian is a conservative system and we can fix energy level for Γ, h; that is, ( (3 p x µ x a 4 x4 = h. To solve it we rewrite ( as dx dt = ± h + µx + a x4. We observe that Eq. (3 for any value of h, a and µ corresponds to the well known incomplete elliptic integral of first kind (see [], whose solutions are no always elementary functions. In order to obtain a manageable equation with a not too complicated base fields, and get explicit solutions for Eq. (3, we should restrict the values of h. Furthermore, to be able to compute Galois groups for the variational equations along such particular solutions, we should stay with these solutions. In particular we choose the energy level h =, then the restricted last equation becomes (4 dx dt = ± µx + a x4.

9 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE 9 µ Now, observe that (4 has three equilibrium points given by x =, ±i a. To avoid triviality, we exclude these points and assume that x(t is not a constant. In a similar way, to obtain algebraic solutions over a suitable differential field in equation (3 with ah, we can choose the energy level h = µ 4a. restricted equation becomes ± dx dt = µ a + a x. Thus, the µ Observe that the last equation has two equilibrium points given by x = ±i a. To avoid triviality, we exclude these points and assume that x(t is not a constant. Recall that similar results are obtained for the respective ĥ in Γ 3. Next, we present a rigorous proof of the non-integrability of Hamiltonian ( using the Morales-Ramis theory. Since this theorem gives conditions for integrability in terms of the differential Galois group of the variational equations along a particular solution, we compute the linearized system of (7 along a non-constant solution of the dynamical system given above. In a general frame we assume γ(t = (x(t, y(t, p x (t, p y (t Γ,3 is a nonstationary solution of (7, then the variational equations ξ = A(tξ of (7 along γ(t are given by ξ ξ ξ 3 ξ 4 = µ + 3ax(t + (a + by(t (a + bx(ty(t (a + bx(ty(t µ + (a + bx(t + 3ay(t Assume γ(t = (x(t,, p x (t, Γ is a non-stationary solution of (7, then the variational equations ξ = A(tξ of (7 along γ(t are given by ξ ξ ξ 3 ξ 4 = µ + 3ax(t µ + (a + bx(t The second order variational system is composed of two uncoupled Schrödinger equations (5 ξ = ( 3ax (t + µ ξ, ξ ξ ξ 3 ξ 4. ξ ξ ξ 3 ξ 4. (6 ξ = ( (a + bx (t + µ ξ. It is well known that the tangential variational equations do not provide conditions of non-integrability of Hamiltonian systems with two degrees of freedom, see [9]. It is easy to see that when we replace a instead b into Eq. (6, the variational normal and tangential equations are the same. For this reason, from now on we let out the tangential equations.

10 ACOSTA-HUMÁNEZ, ALVAREZ-RAMíREZ, STUCHI We now assume γ(t = (x(t, x(t, p x (t, p x (t Γ 3 is a non-stationary solution of (7, then the variational equations of (7 along γ(t becomes ξ ξ ξ 3 ξ 4 = In this case the normal equation yields µ + (4a + bx(t (a + bx(t (a + bx(t µ + (4a + bx(t (7 ξ = ( (a bx(t + µ ξ. The computation is carried out considering several combinations of the signs of the parameters µ and a with µ, and distinguishing between cases a = and a. We also assume without loss of generality that the constants of integration of the normal variational equations are zero Case a = : rational non-integrability. In this case the normal variational equation (6 takes the form (8 ξ = ( bx (t + µ ξ, respectively. A simple computation shows that x(t = cos( µt with h = / is a particular solution of the variational equations on the symmetry plane Γ. Therefore we obtain that (8 along x(t takes the form ξ = ( b cos ( µt + µ ξ. Taking into account that cos ( µt = ( + cos( µt, the last equation becomes (9 ξ = ( b + µ + b cos( µt ξ which is a Mathieu equation. In a similar way, we obtain an equivalent result for Γ 3 by replacing b b. We obtain for b and µ, according to Theorem 4, that Eq. (9 is non-integrable, see [3]. Therefore, we have proved the following result. Theorem. The Hamiltonian system (7, with a =, b and µ is non-integrable through rational first integrals. Now in the case µ =, a =, h = and ĥ the normal variational Eq. (6 and Eq. (7 are no-integrables because they correspond to quantic harmonic oscillators with zero energy level for both planes Γ and Γ 3, see [3, 6]. Finally, the cases h = ĥ = are not considered because we obtain for both planes Γ and Γ 3 that the normal variational equations (6 and (7 are trivially integrables. The first one corresponds to a differential equation with constant coefficients and the second one corresponds to a Cauchy-Euler equation, which are always integrable with abelian differential Galois group, see [4, 6]. Thus Morales-Ramis Theorem does not give any obstruction to the integrability. ξ ξ ξ 3 ξ 4.

11 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE h = h = µ 4a µ a µ x(t = ±i a tanh. x(t = ±i cosh( µt µ a µ x(t = ±i a tan. x(t = ± i cos(i µt µ a µ 3. x(t = ± sin ( i µt x(t = ± a cot µ a µ 4. x(t = ± sinh ( x(t = ± µt a coth Table Particular solutions over γ(t. ( µ t ( µ i t ( µ i t ( µ t 3.4. Case a : meromorphic non-integrability. We consider γ(t = (x(t,, p x (t, Γ which is a non-stationary solution of (7. Next, we solve Eq. (4 for h = and h = µ which give us possible particular solutions over γ(t, shown 4a in the Table. These solutions are complex functions where also t is a complex number. Moreover, we note that is not possible to get algebraic particular solutions for h with µ =, thus we exclude this case. So, we consider only one algebraic particular solution, given in the previous table, for each one of the energy levels (h =, h, to construct the normal variational equations. The remaining cases are obtained through the transformations t it and t π t. Furthermore, the particular solutions over Γ 3 can be obtained by the change a a + b into the particular solutions over Γ. We recall that T l, the l-th triangular number, is defined by T l = l = l(l +, l Z+. Since a C, instead of working with the space of parameters {(a, b, µ C 3 } we work with W := {κ C : κ = b/a, a }. We define the following subsets of W with µ given as: { Λ := {κ : κ = T l }, Λ := κ : κ = T } l. + T l We define by Λ := Λ Λ to introduce our main result. Theorem. Assume aµ. If κ / Λ then the Hamiltonian ( is nonintegrable through meromorphic first integrals. Proof. We consider only the case h = ĥ = in the invariant planes Γ and Γ 3, respectively. We begin with the variational equation near Γ. We use (a with h = to get x(t which is given in the Table. Substituting the expression for x(t in (6 with κ = b/a, we obtain ( ( + κ ( ξ = µ cosh ( + ξ, µt

12 ACOSTA-HUMÁNEZ, ALVAREZ-RAMíREZ, STUCHI where we assume coefficients in the differential field K = C(tanh( µt. Through the change τ = µt and letting ξ = dξ, ( is transformed into the Schrödinger dz equation with Pöschl-Teller potential ( ( ξ ( + κ = cosh (τ + ξ, which is transformed into the Legendre s differential equation with parameters l and µ such that ( + κ = l(l + and µ = ±. Since the second parameter is an integer, by Theorem 3, we know that Eq. ( is integrable if only if ± l Z or l Z, because we are only getting the case of Theorem 3. So, l(l + must be an integer, which implies that Eq. ( has Liouvillian solutions if and only if κ Λ. Moreover, the differential Galois group DGal(L/K of this normal variational equation Eq. ( for κ Λ is isomorphic to the additive group C, which is a connected abelian group; while for κ / Λ the differential Galois group DGal(L/K is an algebraic subgroup isomorphic to SL(, C and thus non abelian. Therefore, by Morales-Ramis Theorem we can conclude that if κ / Λ then the Hamiltonian AGK is not integrable through meromorphic first integrals. The normal variational equations for the second, third and fourth cases of the particular solutions x(t are transformed into Eq. ( through the combination of the changes of variables t it and t π t. So we get similar results of integrability of Eq. ( and consequently similar conditions of non-integrability for the Hamiltonian AGK. From here we consider the variational equation near Γ 3. We use (b to obtain x(t which is for ĥ = given in the Table with a+b instead of a. Now we substitute x(t into the normal equation (7. It has the form ( ξ = µ ( ( κ/( + κ cosh ( + ξ. µt As previously, we rescaled the independent variable to get a Schrödinger equation with Pöschl-Teller potential, which can be transformed into a Legendre s differential equation. From (5, we can conclude that κ + κ = T l or equivalently κ = T l. + T l A similar analysis applied to the above case leads to the non-integrability conditions through meromorphic first integrals for Hamiltonian AGK with κ / Λ. Finally, we conclude that the Hamiltonian AGK ( is non-integrable through meromorphic first integrals for κ / Λ. We observe that Λ = Λ Λ = {,, }. Furthermore, the values b =, a, a given by the Theorem 9 correspond to κ =,,, respectively. Then, the former satisfy the necessary condition for meromorphic integrability. We would now like to make some remarks. In order to apply Morales-Ramis theory as non-integrability criterium, it is enough to exhibit an orbit for which the normal variational equation is not abelian. Taking this into account and since µ is not constant, we have only been considering the proof of the last theorem for the both µ cases h = ĥ =. Moreover, if we replace any particular solution for h = into the 4a normal variational equation, Eqs. (6 and carries out combinations of the change of

13 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE 3 variables t it and t π t, we obtain the following normal variational equation ( (3 ξ ( + κ = cosh (τ + κ + ξ, which is transformed into the Legendre s differential equation with parameters l and µ such that ( + κ = l(l + and µ = κ +. We get equivalent conclusions for µ ĥ = in the variational equation (7. 4(a + b Finally, we remark that for µ = and h = ĥ =, the normal variational equations (6 and (7 correspond to the Cauchy-Euler equation for Γ and Γ 3, respectively, given by η = ( + κ t η, η 3 = κ ( + κt η 3, which are always integrable with abelian differential Galois group, see [4, 6]. Thus Morales-Ramis Theorem does not give any obstruction to the integrability. Furthermore for h and ĥ the particular solution is an elliptic integral of first kind that is not algebraic. In this case the variational equations are not in a suitable differential field to apply Morales-Ramis Theorem. 4. Numerical experiments. In this section, we present a quick view of the system behaviour through the Poincaré section technique. We have chosen the Poincaré section (y =, p y >. The values of the parameters were chosen as to view the system behaviour at the well known integrable cases, and then by moving the values of one of the parameters (a or b, carry the system to a region were KAM tori start being destroyed. The energy and µ are maintained constant for each set of experiments. The first example stemmed from the search of Armbruster et. al reported in Figure 5 of [8] which has the following set of parameters: h = 5.7, µ = 5, a = and b =.5. We here have worked with b going from zero to five. The Figures and show the evolution from the integrable case a = and b = to a = and b =.5 where the superior part of the section is completely chaotic. This is the case of Figure 5 of [8] where they have shown the chaotic layer only. The intermediary figures show how a meagre top and bottom region of KAM tori, b =., evolves to partially chaotic region (.3, and finally to regions of no tori at all since trajectories become unbounded. Some KAM tori remain in the horizontal regions which we show in Figure. In Figure 3 the integrable case a =, b = with h = 3.5 and µ = 5 is shown. Note now that the section of the tori have rotated π/ and this indicates a different way of breaking down. In fact, in Figure 3, b =.8 and keeping the remaining parameters the same, we see a stable periodic orbit at the origin and two unstable ones at each side of the figure. Finally for b = 6: the top and bottom cuts of the KAM tori break down with the internal trajectories escaping to non closed regions of the phase space. This behaviour can be seen in Figure 3. So far, although the route to chaos seems different for the two cases, it seems that scape trajectories play an important role. The last integrable case is b = a can be seen in Figures 4. We have chosen as the starting value a = and b =, the remaining values are h = 3.5 and µ = 5. The cuts of the Poincaré section start centered and for a =.8 there is a horizontal bifurcation and the orbits in the central region starts to scape, Figure 4 (middle. For

14 4 ACOSTA-HUMA NEZ, ALVAREZ-RAMı REZ, STUCHI a =.5 the central region is more eroded and all the bordering trajectories escape, Figure 4 (bottom. Finally an example o chaotic behaviour that contradicts all the previous experiments. As can be seen in Figures 5, there is no escaping orbits among the chaotic ones. We started with a =, b = and the value of the Hamiltonian h =. and µ = maintained constant through the evolution of b from to 5. The section in a = and b = is similar to Figure 4 (top. When b =.5 there is already a major bifurcation in the vertical direction. Finally when b = 5 there is chaos at large with some islands still left. Apparently the value of µ is responsible for this behaviour Fig.. Poincare section (x, px : The top figure shows the integrable case b =, a =, µ = 5 and h = 5.7, while in the bottom figure is showing the same values of the parameters as last one with the exception of b =.. Note the wedge of KAM tori which appears in the vertical direction. Acknowledgments. We thank the anonymous referees for their valuable comments and suggestions helped us to improve this paper. We also thank Thierry Combot for valuable discussions and providing his Maple code for testing homogeneous

15 5 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE Fig.. Poincare section (x, px : b =.3 and remaining parameters as in Figure (top; The vertical tori break through scape in the open phase space; we reach b =.5 (bottom, and now we have the same conditions as Armbruster et. al [8] (Figure 5; the horizontal empty region show by Armbruster et. al [8] is filled with remnant tori. Pictures available on the original file Fig. 3. Poincare section (x, px : the integrable case b = a with a =, b =, µ = 5, h = (top; Note that the tori have rotated to the vertical axis; for b =.8 there is a major bifurcation in the vertical direction (middle and for b = 6 there appears chaos in the escaping trajectories (bottom. potential and Maria Przybylska who pointed out us some papers on the integrability conditions for non-homogeneous potentials. P. Acosta-Huma nez wishes to thank the Department of Mathematics, UAM-Iztapalapa, Mexico, where some of this work was carried out, and also to Universidad Simo n Bolivar (Barranquilla, Colombia for its support in the final stage of this research. M. Alvarez-Ramı rez is partially supported by the grant: Red de cuerpos acade micos Ecuaciones Diferenciales. Proyecto sistemas dina micos y estabilizacio n. PROMEP -SEP, Mexico.

16 6 ACOSTA-HUMA NEZ, ALVAREZ-RAMı REZ, STUCHI Fig. 4. Poincare section (x, px : the integrable case b = a with a =, b =, h = 3.5 and µ = 5 (top and we can see the symmetrical KAM tori; taking a to a =.8 (middle there is a horizontal bifurcation and the central tori start escaping; with a =.8 (bottom the central part becomes more chaotic and most of the bordering trajectories escape.

17 7 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE Fig. 5. Poincare section (x, px with h =., µ = starting at the integrable case a =, b =. The integrable section is symmetrical as top Figure 4 with a =, b =.5 (top there is a vertical bifurcation and with b = 5 (bottom the system is quite chaotic but without escape trajectories as the case with µ = 5. REFERENCES [] R. Abraham and J. E. Marsden, Foundations of mechanics, Benjamin, Reading, (978. [] M. Abramowitz and I. A. Stegun, Handbook of mathematical functions: with formulas, graphs, and mathematical tables, A Wiley-Interscience Publication, John Wiley & Sons, Inc., New York; National Bureau of Standards, Washington, DC, (984. [3] P. Acosta-Huma nez and D. Bla zquez-sanz, Non-integrability of some hamiltonians with rational potential, Discrete Contin. Dyn. Syst. Ser. B, (8, pp , doi:.3934/dcdsb [4] P. B. Acosta-Huma nez, Galoisian approach to supersymmetric quantum mechanics. the integrability analysis of the schro dinger equation by means of differential galois theory, VDM Verlag Dr. Mu ller, Berlin, (. [5] P. B. Acosta-Huma nez, J. La zaro, J. Morales-Ruiz, and C. Pantazi, On the integrability of polynomial vector fields in the plane by means of picard-vessiot theory, Discrete Contin.

18 8 ACOSTA-HUMA NEZ, ALVAREZ-RAMı REZ, STUCHI Fig. 6. Poincare section (x, px with h =., µ = 3 starting at the integrable case a =, b = and the second a =., b =. Note how rapidly chaos appears with a small change in the a parameter. Dyn. Syst., 35 (5, pp , doi:.3934/dcds [6] P. B. Acosta-Huma nez, J. Morales-Ruiz, and J. A. Weil, Galoisian approach to integrability of schro dinger equation, Rep. Math. Phys., 67 (, pp , doi:.6/s344877(69-. [7] P. Andrle, A third integral of motion in a system with a potential of the fourth degree, Phys. Lett. A, 7 (966, pp [8] D. Armbruster, J. Guckenheimer, and S. Kim, Chaotic dynamics in systems with square symmetry, Phys. Lett. A, 4 (989, pp. 46 4, doi:.6/375-96( [9] A. Bostan, T. Combot, and M. E. Din, Computing necessary integrability conditions for planar parametrized homogeneous potentials, In Proceedings of the 39th International Symposium on Symbolic and Algebraic Computation, (4, pp , doi:.45/ [] G. Contopoulos, Galactic dynamics, Princeton University Press, (988. [] A. A. Elmandouh, On the dynamics of armbruster-guckenheimer-kim galactic potential in a rotating reference frame, Astrophys Space Sci, 36 (6, p., doi:.7/s59-6-

19 HAMILTONIAN WITH A POTENTIAL OF THE FOURTH DEGREE [] J. Hietarinta, Direct methods for the search of the second invariant, Phys. Rep., 47 (987, pp [3] I. Kaplansky, An introduction to differential algebra, Hermann, (957. [4] E. Kolchin, Differential algebra and algebraic groups, Academic Press, Pure and Applied Mathematics, New York - London, 54 (973. [5] J. Llibre and L. Roberto, Periodic orbits and non-integrability of armbruster-guckenheimerkim potential, Astrophys Space Sci, 343 (3, pp , doi:.7/s [6] A. Maciejewski and M. Przybylska, Darboux points and integrability of hamiltonian systems with homogeneous polynomial potential, J. Math. Phys., 46 (5, p. 33 pp, doi: [7] J. Morales-Ruiz, Differential galois theory and non-integrability of hamiltonian systems, Progress in Math., Birkhauser, Verlag, Basel, 78 (999. [8] J. Morales-Ruiz and J. P. Ramis, Integrability of dynamical systems through differential galois theory: a practical guide, Differential algebra, complex analysis and orthogonal polynomials, Contemp. Math., Amer. Math. Soc., Providence, RI, 59 (, pp. 43, [9] J. J. Morales-Ruiz and J. P. Ramis, Galoisian obstructions to integrability of hamiltonian systems i, Methods Appl. Anal., 8 (, pp , doi:.43/maa..v8.n.a3. [] J. J. Morales-Ruiz and J. P. Ramis, Galoisian obstructions to integrability of hamiltonian systems ii, Methods Appl. Anal., 8 (, pp. 97, doi:.43/maa..v8.n.a4. [] K. Nakagawa, Direct construction of polynomial first integrals for hamiltonian systems with a two-dimensional homogeneous polynomial potential, Dep. of Astronomical Science, The Graduate University for Advanced Study and the National Astronomical Observatory of Japan, Ph.D. Thesis, (. [] G. Pöschl and E. Teller, Bemerkungen zur quantenmechanik des anharmonischen oszillators, Zeitschrift für Physik., 83 (933, pp [3] F. Simonelli and J. P. Gollub, Surface wave mode interactions: effects of symmetry and degeneracy, J. Fluid Mech., 99 (989, pp , S [4] E. Titchmarsh, Eigenfunction expansions associated with second-order differential equations, part i, (nd edition. Oxford University Press, Oxford, (96. [5] M. van der Put and M. Singer, Galois theory of linear differential equations, Grundlehren der Mathematischen Wissenschaften, Springer Verlag, Berlin, 38 (3. [6] S. Wiggins, Introduction to applied nonlinear dynamical systems and chaos, New York: Springer-Verlag, (99. [7] H. Yoshida, Nonintegrability of the truncated toda lattice hamiltonian at any order, Comm. Math. Phys., 6 (988, pp