Non-integrability of flail triple pendulum arxiv:.6045v [nlin.cd] 6 Nov 0 Maria Przybylska Institute of Physics University of Zielona Góra Licealna 9 PL-65 47 Zielona Góra Poland e-mail: M.Przybylska@proton.if.uz.zgora.pl and Wojciech Szumiński Institute of Physics University of Zielona Góra Licealna 9 PL-65 47 Zielona Góra Poland e-mail: uz88szuminski@gmail.com May 5 04 AMS Subject Classification: 70H07; 70H; 70F07 Abstract We consider a special type of triple pendulum with two pendula attached to end mass of another one. Although we consider this system in the absence of the gravity a quick analysis of of Poincaré cross sections shows that it is not integrable. We give an analytic proof of this fact analysing properties the of differential Galois group of variational equation along certain particular solutions of the system. Introduction In this paper we study dynamics of a special type of triple pendulum which we called a flail. It consists of three pendula. The first one is attached to a fixed point and to its end mass the other two pendula are joined. The lengths and the masses of pendula are l l l and m m m respectively see Fig.. We consider this system in the absence of gravity. Similar problem of dynamics of a simple triple pendulum in the absence of gravity field was analysed in []. The configuration space of our system is a torus T = S S S with local coordinates (θ θ θ ) mod π. We chose the point of suspension of the first pendulum as the origin and angles are measured from vertical line as it is shown in Fig.. The Lagrange function has the following form L = ( l µ θ + l m θ + l m θ ) + l l m cos(θ θ ) θ θ + l l m cos(θ θ ) θ θ (.) where µ = m + m + m. Let us remark that the dynamics of this system can be interpreted as a geodesic motion on torus T. The system has S symmetry. This is why the Lagrange function depends on differences of angles only. Hence it is reasonable to introduce new variables defined by γ = θ θ γ = θ θ γ = θ (.)
θ l y x m θ l θ m l m Figure : Geometry of the flail pendulum and corresponding momenta π = (l m + l m ) γ + l m γ + [l m + l m + l l m cos γ + l l m cos(γ + γ )] γ π = l m γ + l m γ + l m [l + l cos(γ + γ )] γ π = [l m + l m + l l m cos γ + l l m cos(γ + γ )] γ + l m [l + l cos(γ + γ )] γ + [l m + l m + l µ + l l m cos γ + l l m cos(γ + γ )] γ. Now γ is a cyclic variable and π is a first integral. Thus after this reduction the system has two degrees of freedom and its dynamics is determined by Hamiltonian H = [l µ(l m (π π ) + l m π ) + l l m m (π c) l (l l m π cos γ + l m cos(γ + γ )(l m π (c π ) + l l m (π π ) cos(γ + γ )) l l m m (π π ) cos γ (l (π c) + l π cos(γ + γ )))]/[l l l m m B] (.) where B := µ m cos γ m cos (γ + γ ) and c := π is a parameter. The aim of this paper is an analysis of integrability of this reduced system. To get an idea about the complexity of dynamics we made several Poincaré cross sections. In Figures and 4 we show such cross sections. We take γ = 0 and π > 0 as the cross section plane. The phase portraits are presented in (γ y ) plane. The cross sections are made for the following values of parameters: π = c = m = m = m = l = l = l =. (.4) Figures are ordered according to increasing energy of the system. They show that for chosen values of parameters the system is not integrable. The main problem considered in this paper is to prove that for a wide range of the parameters the system is in fact non-integrable. Our main result we formulate in the following theorem. Theorem.. If the reduced flail system governed by Hamiltonian (.) has parameters satisfying m = 0 c = 0 l = l and m l = m l ; or m = m l = l and c = 0 then it is nonintegrable in the class of meromorphic functions of coordinates and momenta.
Figure : The Poincare sections for E = 0.0 (on the left) and E = 0.0005 (on the right) Figure : The Poincare sections for E = 0.00 (on the left) and E = 0.0 (on the right) Figure 4: The Poincare sections for E = 0.0 (on the left) and E = 0.0 (on the right)
The plan of the paper is the following. Section contains the non-integrability proof and is divided into two subsections. In Section some remarks are given. Explicit forms of some complicated expressions are given in Appendix A. For convenience of readers in Appendix B basic information about linear second order differential equation with rational coefficients and Kovacic algorithm are given. Non-integrability proof In this section we prove Theorem. using the Morales-Ramis theory see [ ]. We formulate the main theorem of this theory in the following simple settings. Let us consider a complex Hamiltonian system. We assume that the underlying phase space is C n which is considered as a linear symplectic space equipped with the canonical symplectic form n ω = dq i dp i i= where x = (q... q n p... p n ) are the canonical coordinates. Let H : C n C be a holomorphic Hamiltonian and d dt x = v H(x) v H (x) = I n x H x C n t C (.) the associated Hamilton equations. Here I n denotes n n symplectic unit matrix. Let t ϕ(t) C n be a non-equilibrium solution of (.). The maximal analytic continuation of ϕ(t) defines a Riemann surface Γ with t as a local coordinate. Variational equations along ϕ(t) have the form d dt ξ = A(t)ξ A(t) = v H (ϕ(t)). (.) x We can attach to this equations the differential Galois group G. Roughly speaking G is an algebraic subgroup of GL(n C) which preserves polynomial relations between solutions of (.). As a linear algebraic group among other things it is a union of a finite number of disjoint connected components. One of them containing the identity is called the identity component of G and is denoted by G 0. For a precise definition of the differential Galois group and differential Galois theory see e.g. [5 0 4]. It appears that the integrability of the considered system manifests itself in properties of the differential Galois group of the variational equations. For Hamiltonian system and integrability in the Liouville sense this relation is particularly elegant. For these systems G is an algebraic subgroup of Sp(n C). In nineties of XX century Morales-Ruiz and Ramis showed that the integrability in the Liouville sense imposes a very restrictive condition on G 0 [ ]. Theorem. (Morales-Ruiz and Ramis). Assume that a Hamiltonian system is meromorphically integrable in the Liouville sense in a neigbourhood of a phase curve Γ. Then the identity component of the differential Galois group of NVEs associated with Γ is Abelian. In order to apply this theorem we need an effective method which allows to determine properties of the differential Galois group of linear equations. In the investigated system variational equations split into two subsystems of linear equations. Each of these subsystems can be transformed into equivalent second order equation with rational coefficients. For such equation there exists the Kovacic algorithm [7] see also Appendix B which allows to determine the differential Galois group. We apply this algorithm in our considerations. 4
. Case m = m and l = l At first we make the following non-canonical transformation 0 0 0 γ x γ π = A x y A = 0 0 0 0 l l. (.) π y 0 0 0 After this transformation equations of motion take the form where ẋ = J n x H x = [x x y y ] T (.4) J n := BIB T B = A H(x) := H(Ax). (.5) The explicit form of Hamiltonian in these variable is following H = [l m m (cl + l y l (y + y )) + l (l µ m y + m (l y l y ) ) ( ) l (m cos(x x ) l m y (cl + l y l (y + y )) + l m (l y l y ) cos(x x ) + l m m ( l y + l y ) ( l y + l ( c + y + y ) + l y cos(x x )) cos x )] + l l m y cos x /[l l l m m (µ m cos(x x ) m cos x )]. (.6) An explicit form of vector field v(x) := J n x H is given in Appendix A see equation (A.). These equations are still Hamiltonian but written in non-canonical variables. Let us assume that the following conditions are fulfilled. Then system (.4) has an invariant manifold m l = m l l = l (.7) N = {(x x y y ) C 4 x = y = 0 y = cl /(l l )}. In fact restricting the right hand sides of system (.4) to N and putting m = m l /l we obtain c ẋ = 0 ẋ = l m (l l ) ẏ = 0 ẏ = 0. (.8) Hence solving the above equations we obtain our particular solution C t ϕ(t). Let X := [X X Y Y ] T denote variations of [x x y y ] T. The variational equations along ϕ(t) take the form a 0 a a 4 d dt X = AX A = a 0 a a 4 a 0 a a 4. (.9) a 4 0 a 4 a 44 The explicit forms of entries a ij are given in equation (A.) in Appendix A. We note that equations on X Y and Y do not depend on X i.e. variable X decouples from remaining variables. Moreover one can observe that equations for Y and for Y are dependent see appropriate coefficients a ij. This fact is connected with the existence of a first integral of variational equations. The expansion of energy integral along the particular solution gives the first integral of variational equations h := grad H(ϕ) X = c l l m (l l ) [ l Y + (l + l )Y ]. (.0) 5
At the level h = 0 we have Y = l l + l Y. (.) Then substituting this into the sub-system of variational equations for X Y and Y and eliminating Y we obtain the following second order linear equation for X with coefficients Ẍ + aẋ + bx = 0 (.) a = [4c(l + l )(l l + l (l + l ) cos x )(l l m + l (l + l )m + l l m cos x ) sin x ] /[l (l l )(4l l m + l (l + l )(l m + (l + l )m ) + 4l l l (l + l )m cos x )(l m + (l + l )m (l + l )m cos(x ))] b = [c (l l + l (l + l ) cos x )(4l l l (l + l )m + (4l l m + l (l + l )(l m + (l + l )m )) cos x )]/[l l (l l ) m (4l l m + l (l + l )(l m + (l + l )m ) + 4l l l (l + l )m cos x )(l m + (l + l )m (l + l )m cos(x ))]. In order to transform this equation into an equation with rational coefficients we use the following transformation t z = cos x (t). (.) This change of variable together with transformations of derivatives d dt = ż d dz d d = (ż) dt dz + z d dz ż = sin x ẋ = c sin x l m (l l ) (ż) = c ( z ) l m (l l ) z = cos x ẋ = c z l m (l l ) convert equation (.) into X + p(z)x + q(z)x = 0 X = X (.4) where the prime denotes the derivative with respect to z. The explicit forms of coefficients p(z) and q(z) are the following p = ( z + ża) = z ż q = ḃ z = where we introduced notation σz (z )σ αγ 4σ 4 + σ(4z + σ + αγ) z + (σz + )(4σ + z(4 + σ(σ + αγ))) (z )[σ(z ) αγ)(4 + σ(4z + σ + αγ)] α = l l β = l l γ = m m σ := α + β. These parameters take non-negative values and we assume that they do not vanish simultaneously. Now we make the classical Tschirnhaus change of dependent variable [ X = w exp z ] p(ζ) dζ (.5) z 0 which transforms (.4) into its reduced form w = r(z)w r(z) = q(z) + p (z) + 4 p(z) (.6) 6
with coefficient [ ] [ ] r(z) = 6 (z ) + (z + ) 4 (z δ) + (z + δ) (z ρ) + σ δ ( 5z δ 4 ) + 4σ z ( 4 + δ z ) + 4σ ( 0 + δ 5z ) + 64z σ (z )(z δ )(z ρ) where δ := + αγ σ ρ := σ 4 δ σ. (.7) Let us notice that coefficient r(z) depends only on two parameters δ and σ. We assume that both are positive and δ. If δ = and σδ = then equation (.6) has six pairwise different regular singular points z = ± z 4 = ±δ z 5 = ρ z 6 =. The respective differences of exponents at these points are following = = = 4 = 0 5 = 6 =. (.8) Now we can prove the following result. Lemma.. For the differential Galois group of equation (.6) is SL( C). δ = and δσ = (.9) Proof. We prove this lemma by a contradiction. Thus let us assume that G = SL( C). Then according to Lemma B. G is either a finite subgroup of SL( C) or a subgroup of dihedral group or a triangular group of SL( C). At first we show that the first two possibilities do not occur. As 4 = 0 local solutions in a neighbourhood of z = z and z = z 4 contain a logarithmic term. More precisely two linearly independent solutions w and w of (.6) in a neighbourhood of z {z z 4 } have the following forms w (z) = (z z ) ρ f (z) w (z) = w (z) ln(z z ) + (z z ) ρ h(z) (.0) where f (z) and h(z) are holomorphic at z and f (z ) = 0. The local monodromy matrix corresponding to a continuation of solutions along a small loop encircling z counterclock-wise gives rise to a triangular monodromy matrix M = [ ] πi 0 for details see [9]. A subgroup of SL( C) generated by M is not finite and thus also the differential Galois group G of (.6) is not finite. As M is not diagonalisable G is not a subgroup of dihedral group. Thus G is a subgroup of the triangular group of SL( C). In order to check if it is really the case we apply the Kovacic algorithm for details see Appendix B. According to it if G is a subgroup of the triangular group of SL( C) the first case of this algorithm occurs i.e. the equation admits an exponential solution. To check this possibility we determine sets of exponents Thus we have E = E = E i := {( ± i )/} for i =... 6. { 4 } E = E = 4 { } { E 5 = } E 6 = {0 }. (.) 7
Next according to the algorithm we look for elements e = (e e e e 4 e 5 e 6 ) of Cartesian product E = 6 i= E i such that d(e) := e 6 5 i= e i N 0 (.) where N 0 is the set of non-negative integers. In our case there exists only one e E with this property namely ( e = 4 4 ) (.) for which d(e) = 0. Thus according to the algorithm we look for a polynomial P = 0 of degree d(e) such that it is a solution of the following equation where P + ωp + (ω + ω r)p = 0 (.4) ω = 5 i= e i z z i. In the considered case we have P = so equation (.4) gives equality [ ] (σz + ) (σ δ + 4)z + 4σ = 0 which cannot be fulfilled for σ > 0. This end the proof. Now let us analyse the cases excluded in Lemma.. Case σ = 0 is equivalent to α = l /l = 0 and β = l /l = 0. So it is equivalent to the case of zero length of the first pendulum. In effect the system consists of two simple independent pendula and is obviously integrable. In the case of confluence of singular point z and z 5 we prove the following. Lemma.. Let us assume that δ = σ. (.5) Then the differential Galois group of equation (.6) is SL( C). Proof. Then coefficient r in equation (.6) simplifies to the following form r = 6(z ) 6(z + ) 4(z z ) ρ(z z 4 ) + ρ 6(ρ + )(z + ) + + ρ 6(ρ )(z ) ρ(ρ )(z z ) (.6) where z = z 4 = δ := /σ. Now if δ > then equation (.6) has five singularities z = ± z and z 4 and z 5 =. One can easily check that the difference of exponents at z = z vanishes. Thus the local solutions around this point have the form (.0). One can repeat reasoning from the proof of Lemma.. In order to check whether the first case of the Kovacic algorithm can appear at first we calculate the sets of exponents E i for i =... 5 E = E = { 4 } E = 4 { } E 4 = {0 } E 5 = {0 }. (.7) One can easily check that for each e = (e e e e 4 e 5 ) 5 i= E i we have d(e) = e 5 4 i= 8 e i N 0. (.8)
Thus if δ > according to Kovacic algorithm the differential Galois group of (.6) is SL( C). For δ = i.e. for σ = coefficient r simplifies to r (z) = 7 + 4z. 4( + z ) (.9) Now equation (.6) has only three regular singularities i.e. it is the Riemann P equation. Using the well known Kimura theorem [6] one can check that its differential Galois group is SL( C). Remark.4 An equivalent form of condition σδ = is following l (l + l )(l m + l (m + m )) 4l l m = 0. For m = 0 and l 6= 0 using the Kovacic algorithm one can prove only that if the differential Galois group of variational equations along the mentioned particular solution is not SL( C) then it is a finite subgroup of SL( C). Calculations in first and second case of this algorithm are simple but in the third case becomes really complicated. Poincare sections given in Fig. 5 suggest that also in this case system is not integrable. Figure 5: Poincare sections for m = 0 m = m = l = l = l = π = c = and for E = 0.0095 (on the left) and E = 0.0 (on the right). Cross-section plane γ = 0 The proof of the first part of Theorem. is an immediate consequence of Lemma.. and Theorem... Case l = l and m = m If masses and lengths of two pendula attached to the first one are equal the one can expect that behaviour of such system is regular. But Poincare sections suggest that also in this case the system is not integrable. In order to prove that the system is not integrable we apply again the Morales-Ramis theory. For l = l and m = m the reduced Hamiltonian (.6) and the corresponding equations of motion (.4) simplify considerably. In this case for c = 0 the system possesses the following invariant manifold N = {( x x y y ) C4 x = y = 0} and its restriction to this manifold takes the form y x = 0 x = y = 0 l m 9 y = 0. (.0)
Figure 6: Poincare sections for m = m = m = l = l = l = π = c = (on the left) and E = 0.06 (on the right). Cross-section plane γ = 0 and E = 0.005 These equations determine our particular solution ϕ(t). Variational equations (.) along this particular solution are also simpler. Namely entries ai4 of matrix A vanish for i 4. Now equations for X and Y form a closed subsystem of variational equations. From them we obtain a second order equation for X. It takes form (.) with coefficients 4y (l + l cos x )(l (m + m ) + l m cos x ) sin x l (l m + l (m + m ) + 4l l m cos x )(m + m m cos(x )) y (l + l cos x )(4l l m + (l m + l (m + m )) cos x ). l l4 m (l m + l (m + m ) + 4l l m cos x )(m + m m cos(x )) a= b= Changing independent variable according to (.) we obtain equation of the form (.4) with coefficients 8z α + 6z4 α 4α( + γ) + 6z ( + α ( + γ)) z(6 + γ)( + α ( + γ)) (z )( + z γ)( + 4zα + α ( + γ)) ( + zα)(4α + z( + α ( + γ))) q=. (z )α( + z γ)( + 4zα + α ( + γ)) p= The coefficient r in reduced equation (.6) has the following expansion r= + 4( z z5 ) 6(z ) 6(z + ) 4( z δ ) 4( z + δ ) + 8 + α(5 + δ + α(6 + δ ( + α( + δ )))) 6α( + δ )( + α( + αδ ))(z ) + 8 + α( 5 + α) + α( + α( + α))δ α δ4 δ + α( 5 + δ(δ + α( + αδ( + δ )))) + 6α( + δ )( + α( + αδ ))(z + ) 4αδ( + αδ) ( + δ )(z δ) α(5 + α4 δ4 + α6 δ6 + α ( 6 + 7δ )) δ + α( 5 + δ(δ + α( + αδ( + δ )))) + 4αδ( + αδ) ( + δ )(z + δ) ( + α δ ) ( + α4 δ4 + α ( + δ ))(z z5 ) (.) where z5 = ( + α δ )/(α) and δ = + γ/. Equation (.4) with coefficient r given in (.) has singularities z = ± z4 = ±δ and z5 = ( + α δ )/(α) that are all poles of the second order and infinity has degree for all positive α and δ. At first we consider the generic case. 0
Lemma.5. If δ = /α and δ = then the differential Galois group of equation (.6) with r(z) given by (.) is SL( C). Proof. We apply the Kovacic algorithm. The sets of exponents E i in this case are { E = E = 4 } { } { E = E 4 4 = E 5 = } E 6 = {0 }. (.) Hence in a neighbourhood of z = z and z = z 4 local solutions have the form (.0). Thus we proceed in a similar way as in the proofs of previous lemmas. There exists only one e E = 5 i= E i such that d(e) N 0. This is ( e = 4 4 ) (.) and for it we have d(e) = 0. Then we construct rational functions ω = 5 i= e i z z i and check whether monic polynomial P = satisfies differential equation (.4) for some values of parameters α and δ. This equation is satisfied only when equalities α = 0 + α + α δ = 0 α + α δ = 0 hold but this system of algebraic equations is not consistent. Lemma.6. If δ = /α and δ = then the differential Galois group of equation (.6) with r(z) given by (.) is SL( C). Proof. In this case r given in (.) simplifies to the form r = 6(z ) 6(z + ) 4(z z ) + α 6(α )(z ) + + α 6(α + )(z + ) + α ( + α )(z z ) α (z z 4 ) (.4) where z = α and z 4 = α. Singularities z = ± and z are poles of the second order and z 4 is a pole of the first order. The infinity z 5 = is generically of degree i.e. it is a regular singular point. Its degree changes only for excluded values of α i.e. for α {0 }. The difference of exponents at z = z is zero. Thus local solutions around this point have the form (.0) and the differential Galois group is either a subgroup of the triangular group or it is SL( C). Thus we have to check only the first case of the Kovacic algorithm. The sets of exponents E i are now following E = E = { 4 4 } E = { } E 4 = E 5 = {0 }. Only for e = (/4 /4 / 0 ) 5 i= value of d(e) is an integer equal to d(e) = 0. Inserting rational function ω = 4 i= e i z z i and P = into differential equation (.4) we obtain α(z + α) = 0. This condition is fulfilled only for excluded value α = 0. Thus the differential Galois group is not a subgroup of the triangular group so it is SL( C).
Case when δ = corresponds to γ = 0 and now we prove the following lemma. Lemma.7. For γ = 0 differential Galois group of (.6) is SL(n C). Proof. Using particular solution ϕ defined by (.0) one can easy check that the first and the second case of the Kovacic cannot occur but an analysis of the third case is very onerous. Fortunately in this case we can find another particular solution for which analysis of the differential Galois group of the corresponding variational equations is simpler. Namely Hamilton equation governed by Hamiltonian (.) for m = 0 m = m and l = l have another particular solution lying on invariant manifold N = {(γ γ π π ) γ = 0 π = π } We consider only the level c = 0. It is convenient to make the following non-canonical change of variables γ x 0 0 0 γ π = A x y A = 0 0 0 0 0 0 π y 0 0. (.5) After this transformation equations of motion take the form ẋ = J x Ĥ J = BIB T B = A Ĥ = H(Ax) (.6) where x = [x x y y ] T. Hamilton function of the reduced system in these variables takes the form Ĥ = [4l y + 4l (y + y ) l (cos x ( 4l y (y + y ) + l (y y ) cos x ) + (y y )(l y + l (y + y ) cos x ) cos(x + x ) + l (y + y ) cos(x + x ) )]/[4l l m ( cos(x ) cos[(x + x )])]. In these variables our system has invariant manifold Equations of motion restricted to this manifold read N = {(x x y y ) x = y = 0}. (.7) ẋ = y (l + l + l l cos x ) l l m sin x ẏ = y ((l + l ) cos x + l l ( + cos(x ))) 4l l m sin x. (.8) These equations defines our particular solution and variational equations along it take the form Ẋ a 0 a 0 X Ẋ Ẏ = a a a a 4 X a 0 a 0 Y (.9) Ẏ a 4 a 4 0 a 44 Y
where a = y (l + l cos x ) l l m sin x a = l m ) a = y ((7l + 8l ) cos x + l (0l + 6l cos(x ) + l cos(x ))) 6l l m sin x a = y ((l + l ) cos x + l l ( + cos(x ))) l l m sin a = x l m a 4 = l + l + l l cos x l l m sin x a = y (0l l cos x + 4(l + l ) cos(x ) + l (6l cos(x ) + l ( + cos(4x )))) l l m sin 4 x a = y (l + l cos x ) l l m sin x a 4 = a 4 = y (l l cos x + 4(l + l )( + cos(x )) + l l cos(x )) 6l l m sin 4 x a 44 = y ((l + l ) cos x + l l ( + cos(x ))) l l m sin x Thus we have the closed subsystem of normal variational equations Ẋ = a X + a Y Ẏ = a X + a Y. (.40) This system can be rewritten as one second order linear equation Ẍ + b(t)x = 0 b(t) = y (l + l cos x ) To rationalise this equation we use transformation 4l l4 m sin4 x. (.4) t z = cos x (t). (.4) This transformation together with the following expressions on time derivatives d dt = ż d dz d d = (ż) dt dz + z d dz ż = sin x ẋ = y (l + l + l l cos x ) l l m (ż) = y (l + l + l l z) sin x 4l 4l4 m ( + z ) z = cos x ẋ sin x ẍ = y (l + l + l l z) 4l l m ( + z ) translates equation (.4) into the following one with rational coefficients p = ż z = (z z ) X + p(z)x + q(z)x = 0 X = X (.4) q = ḃ z = 4(z z ) + α ( + α) (z ) α (α ) (z + ) + α (α ) (z z ). Here z = ( + α )/(α) and α = l /l. Equation (.4) has four singularities z = ± z = + α α z 4 =.
Using the transformation (.5) to (.4) we obtain its reduced form w = r(z)w r(z) = q(z) + p (z) + 4 p(z) (.44) where the expansion of coefficient r is the following 7 r = 6(z z ) α (α + ) (z ) + α (α ) (z + ) α (α ) (z z ) (.45) and order of infinity is provided α = 0. At first let us assume that α =. To check the differential Galois group of equation (.44) we apply the Kovacic algorithm. In the first case sets E i are following E = E = {0 } E = E 4 = { i 4 + i }. (.46) 4 Next we select from the Cartesian product E = 4 i= E i those elements e = (e e e e 4 ) for which d(e) = e 4 i= e i N 0. (.47) We have two choices of e E namely ( e () = 0 0 i 4 ) ( i e () = 0 0 4 + i 4 ) + i (.48) 4 that give d(e (i) ) = 0 for i =. Then we construct rational functions ω = i= e i z z i and substitute into equation (.4) setting also P =. This equation is satisfied only for α = 0 for both choices e (i) i =. In the second case of the Kovacic algorithm sets E i are the following { E = E = {4} E = E 4 = i + i } Z = {}. Thus we see that there is no element e = (e e e e 4 ) for that ( d(e) = e 4 i= e i ) N 0. (.49) Similarly in the third case of this algorithm sets E i are following { E = E = {} E = E 4 = 6 + 6ki } k = 0 ± ±... ± n Z = {6} n where n = 4 6. But we see that also there is no e = (e e e e 4 ) with e i E i given above for that ( d(e) = n e 4 i= e i ) N 0. Thus differential Galois group is whole SL( C). For α = we have a confluence of singular points. In this case coefficient r simplifies to the following forms + 7z r = 6(z )( + z). We note that in these cases equation (.44) becomes a Riemman P equation. With a help of the Kimura theorem [6] one can show that the differential Galois group of this equation is SL( C). 4
The second part of our main Theorem. is an immediate consequence of Theorem. and Lemmas.5.6 and.7. Conclusions The integrability analysis for the flail pendulum with parameters satisfying m l = m l (.) was almost finished except the case when m = 0 m = m and l = l. The only integrable case was found for l = 0. There is an open question about integrability in the case when equality (.) does not hold. Poincaré sections suggest non-integrability for generic values of parameters. Figure 7: The Poincaré sections for m = m = m = l = l = l = π = c = and E = 0.04 (on the left) and m = m = m = l = l = l = π = c = and E = 0.0058 (on the right). Cross-section plane γ = 0 Acknowledgments The authors are very grateful to Andrzej J. Maciejewski for many helpful comments and suggestions concerning improvements and simplifications of some results. 5
A Explicit forms of reduced equations and variational equations The explicit form of vector field v = [v v v v 4 ] T in Section. defined by reduced equations of motion (.4) after linear transformation (.) is the following [ v = l l m m (cl + l y l (y + y )) l (m + m + m )(l m y + l m (l y l y )) +l (l l m (l y l y ) cos(x x ) + l m cos x (l m (l (c y y ) + l y ) +l l m y cos x ) + l m m cos(x x )(l (l (c y y ) + l y ) l ( l y ] +l (y + y )) cos x ))] / [l l l ( m m (m + m + m ) + m m (m cos(x x ) + m cos x )) [ v = l (m + m + m )(l y l y ) + l m (cl + l y l (y + y )) + l (l m (l y l y ) cos(x x ) + l m (cl + l y l (y + y )) cos x l m y cos(x x ) ] (l + l cos x ))] / [l l ( l m (m + m + m ) + l m (m cos(x x ) + m cos x )) v = [(l ( l y + l ( c + y + y )) + l l y cos(x x ) + l (l y l y ) cos x ) (l (l (l m + l (4m + m + 4m ))y (l l )m ( l y + l y ) cos x ) cos x sin x l (l + l )(l m y + m ( l y + l y ) cos x ) cos(x ) sin x + l m ( l y + l ( c + y +y )) cos(x ) sin(x ) + l (l (l m l (4m + m + 4m ))y cos x + (l y l y ) (l ((m + m ) + m ) + ( l + l )m cos(x ))) sin x + l (cl + l y l (y + y )) ( l m + l m cos(x )) sin(x ) + l (l + l )(l m y cos x + m ( l y + l y ) cos(x )) sin(x ))] / [l l l (m + m + m m cos[(x x )] m cos(x )) ] v 4 = [(l ( l y + l ( c + y + y )) + l l y cos(x x ) + l (l y l y ) cos x ) (m cos(x x )(l ( l y + l ( c + y + y )) + l (l y l y ) cos x ) + l l y (m [ +m + m m cos(x ))) sin(x x )] / l l l (m + m + m m cos[(x x )] m cos(x )) ] (A.) 6
Entries of matrix in variational equations (.9) are the following a = c(l l + l (l + l ) cos x ) sin x l l (l l )( l m l (m + m ) + (l + l )m cos x ) a = l l m + l (l m + (l + l )m ) + l l (l + l )m cos x l l m ( l m l (m + m ) + (l + l )m cos x ) a 4 = (l l )(l l + l (l + l ) cos x ) l l l ( l m l (m + m ) + (l + l )m cos x ) a = cl (l + l cos x ) sin x l l (l l )( l m l (m + m ) + (l + l )m cos x ) a = l l m l (l m + l (m + m )) + l l m cos x ( l + l cos x ) l l m ( l m l (m + m ) + (l + l )m cos x ) a 4 = (l + l )(l + l cos x ) l l m l l (l m + l (m + m ) (l + l )m cos x ) (A.) a = c (l + l ) sin x (l l ) ( l m (l + l )m + (l + l )m cos(x )) c(l a = + l )(l + l cos x ) sin x l l (l l )( l m (l + l )m + (l + l )m cos(x )) a 4 = a 4 = c(l + l ) sin x l l ( l m (l + l )m + (l + l )m cos(x )) c l sin x (l l ) ( l m (l + l )m + (l + l )m cos(x )) cl a 4 = (l + l cos x ) sin x l l (l l )( l m (l + l )m + (l + l )m cos(x )) c sin x a 44 = l ( l m (l + l )m + (l + l )m cos(x )). B Linear second order differential equation with rational coefficients and Kovacic algorithm Let us consider a linear second order differential equation with rational coefficients w + p(z)w + q(z)w = 0 p(z) q(z) C(z). (B.) A point z = c C is a singular point of this equation if it is a pole of p(z) or q(z). A singular point is a regular singular point if at this point p(z) = (z c)p(z) and q(z) = (z c) q(z) are holomorphic. An exponent of equation (B.) at point z = c is a solution of the indicial equation ρ(ρ ) + p 0 ρ + q 0 = 0 p 0 = p(c) q 0 = q(c). After a change of the dependent variable z /z equation (B.) reads w + P(z)w + Q(z)w = 0 P(z) = ( ) z p + z z Q(z) = ( ) z 4 q. z (B.) We say that the point z = is a singular point for equation (B.) if z = 0 is a singular point of equation (B.). Equation (B.) is called Fuchsian if all its singular points (including infinity) are regular see [5 4]. The Kovacic algorithm [7] allows to decide whether all solutions of equation (B.) with rational coefficients p(z) and q(z) are Liouvillian. Roughly speaking Liouvillian 7
functions are obtained from the rational functions by a finite sequence of admissible operations: solving algebraic equations integration and taking exponents of integrals. For a formal definition see e.g. [7]. If one (non-zero) solution w of equation (B.) is Liouvillian then all its solutions are Liouvillian. In fact the second solution w linearly independent from w is given by w = w w [ exp ] p. Putting [ w = y exp ] p into equation (B.) we obtain its reduced form y = r(z)y r(z) = q(z) + p (z) + 4 p(z). (B.) This change of variable does not affect the Liouvillian nature of the solutions. For equation (B.) its differential Galois group G is an algebraic subgroup of SL( C). The following lemma describes all possible types of G and relates these types to the forms of a solution of (B.) see [7 ]. Lemma B.. Let G be the differential Galois group of equation (B.). Then one of four cases can occur. Case I G is conjugate to a subgroup of the triangular group T = {[ ] a b 0 a a C b C} ; in this case equation (B.) has an exponential solution of the form y = P exp ω where P C[z] and ω C(z) Case II G is conjugate to a subgroup of D = {[ ] } c 0 0 c c C {[ ] } 0 c c c C ; 0 in this case equation (B.) has a solution of the form y = exp ω where ω is algebraic over C(z) of degree Case III G is primitive and finite; in this case all solutions of equation (B.) are algebraic thus y = exp ω where ω belongs to an algebraic extension of C(z) of degree n = 4 6 or. Case IV G = SL( C) and equation (B.) has no Liouvillian solution. Kovacic in paper [7] have formulated the necessary conditions for the respective cases from Lemma B. to hold. At first we introduce notation. We write r(z) C(z) in the form r(z) = s(z) t(z) s(z) t(z) C[z] where s(z) and t(z) are relatively prime polynomials and t(z) is monic. The roots of t(z) are poles of r(z). We denote Σ := {c C t(c) = 0} and Σ := Σ { }. The order ord(c) of c Σ is equal to the multiplicity of c as a root of t(z) the order of infinity is defined by ord( ) := deg t deg s. Lemma B.. The necessary conditions for the respective cases in Lemma B. are the following. 8
Case I. Every pole of r must have even order or else have order. The order of r at must be even or else be greater than. Case II. r must have at least one pole that either has odd order greater than or else has order. Case III. The order of a pole of r cannot exceed and the order of r at must be at least. If the partial fraction expansion of r is r(z) = i a i (z c i ) + j b j z d j (B.4) then i = + 4a i Q for each i j b j = 0 and if then + 4g Q. g = i a i + b j d j j In [7] Kovacic also formulated a procedure called now the Kovacic algorithm which allows to decide if an equation of the form (B.) possesses a Liouvillian solution and to find it in a constructive way. Applying it we also obtain information about the differential Galois group of the analysed equation. Beside the original formulation of this algorithm several its versions and improvements are known [ 8 ]. Now we describe the Kovacic algorithm for the respective cases from Lemma B.. Case I Step I. Let Γ is the set of poles of r. For each c Γ { } we define a rational function [ r] c and two complex numbers α + c α c as described below. (c ) If c Γ and ord(c) = then [ r] c = 0 α + c = α c =. (c ) If c Γ and ord(c) = then [ r] c = 0. If r has the expansion of the form (B.4) then α ± c = ( ± ) + 4a (B.5) (c ) If c Γ and ord c = ν 4 with k (only even orders are admissible in this case) then [ r ]c = a (z c) ν + + d (z c) ( r [ r ] c = b (z c) ν+ + O (z c) ν ) and α ± c = (± ba + ν ). ( ) If the order of r at infinity is > then [ r] = 0 α + = 0 α =. ( ) If the order of r at infinity is then [ r] = 0. 9
If the Laurent series expansion of r at takes the form r = a ( ) z + O z (B.6) then α ± = ( ± ) + 4a. (B.7) ( ) If the order of r at is ν 0 (necessarily even in this case) then r = az ν + + s (B.8) is the indicated part of the Laurent series expansion of r at and r [ r ] ( = bzν + O z ν ). Then α ± = (± ba ν ). (B.9) Step II. For each family s = (s(c) s( )) c Γ where s(c) and s( ) are either + either let we compute If d is a non-negative integer then d := e s( ) α s(c) c. c Γ ( ω(z) = s(c) [ r ] ) c Γ c + αs(c) c z c + s( ) [ r ] (B.0) is a candidate for ω. If there are no such elements equation (B.) does not have an exponential solution and the algorithm stops here. Step III. For each family from step II that gives d N 0 we search for a monic polynomial P = P(z) of degree d satisfying the following equation P + ω(z)p + (ω (z) + ω(z) r(z))p = 0. (B.) If such polynomial exists then equation (B.) possesses an exponential solution of the form y = P exp θ where θ = ω if not equation (B.) does not have an exponential solution. Case II Let Γ be the set of poles of r. Step I. For each c Γ we define E c as follows (c ) If ord(c) = then the set E c = {4} i (c ) If ord c = and r has the expansion of the form (B.4) then { ( E c := + ) ( + 4a )} + 4a Z (B.) (c ) If ord c = ν > then E c = {ν} ( ) If ord > then E = { 4} ( ) If ord = and the he Laurent series expansion of r at takes the form (B.6) then { ( E := + ) ( + 4a )} + 4a Z. (B.) ( ) If ord = ν < then E = {ν} 0
Step II. We consider all families (e c ) c Γ { } with e c E c and at least one of the coordinates is odd. We compute ( ) d := e e c. c Γ and select those families (e c ) c Γ { } for which d(e) is a non-negative integer. If there are no such elements Case II cannot occur and the algorithm stops here. Step III. For each family giving d N 0 we define ω = ω(z) = c Γ e c z c (B.4) and we search for a monic polynomial P = P(z) of degree d satisfying the following equation P + ωp + (ω + ω 4r)P + (ω + ωω + ω 4rω r )P = 0. (B.5) If such a polynomial exists then equation (B.) possesses a solution of the form y = exp θ where θ ψθ + ψ + ψ r = 0 ψ = ω + P P. If we do not find such polynomial then Case II in Lemma B. cannot occur. Case III Step I. For each c Γ { } we define a set E c of integers as follows (c ) If c is a pole of order then E = {}. (c ) If ord c = and r has the expansion of the form (B.4) then { E c = 6 + k } m + 4a k = 0 ± ±... ± Z. m Here and below in this case m {4 6 }. (B.6) ( ) For c = when the Laurent expansion of r at infinity is of the form (B.6) we define independently on order { E = 6 + k } m + 4a k = 0 ± ±... ± Z. (B.7) m Obviously it can appear that a = 0. Step II. For e E we calculate d(e) := m ( ) e e c. c Γ We select those elements e E for which d(e) is a non-negative integer. If there are no such elements Case II cannot occur and the algorithm stops here. Step III. We consider all families (e c ) c Γ { } with e c E c that give d(e) = n N 0. For each such family we define ω = ω(z) = m c Γ e c z c. Next we search for a monic polynomial P = P(z) of degree n satisfying a differential equation of degree m + defined in the following way. Put P m = P. Then calculate P i for i = m m... 0 according to the following formula P i = SP i + [(m i)s Sω]P i (m i)(i + )S rp i+ (B.8)
where S = c Σ (x c). Then P = 0 gives the desired equation for P. If such polynomial exists then equation (B.) possesses a solution of the form y = exp θ where θ is a solution of the equation n i=0 S i P i (m i)! θi = 0. If we do not find such polynomial then we repeat these calculations for the next element m {4 6 }. If for all m such polynomial does not exist then Case III in Lemma B. cannot occur. References [] A. Baider R. C. Churchill D. L. Rod and M. F. Singer. On the infinitesimal geometry of integrable systems. In Mechanics day (Waterloo ON 99) volume 7 of Fields Inst. Commun. pages 5 56. Amer. Math. Soc. Providence RI 996. [] F. Beukers. Differential Galois theory. In From number theory to physics (Les Houches 989) pages 4 49. Springer Berlin 99. [] A. Duval and M. Loday-Richaud. Kovaˇciˇc s algorithm and its application to some families of special functions. Appl. Algebra Engrg. Comm. Comput. (): 46 99. [4] E. L. Ince. Ordinary Differential Equations. Dover Publications New York 944. [5] I. Kaplansky. An introduction to differential algebra. Actualités Scientifiques et Industrielles No. 5 Publications de l Institut de Mathématique de l Université de Nancago No. V. Hermann Paris second edition 976. [6] T. Kimura. On Riemann s equations which are solvable by quadratures. Funkcial. Ekvac. :69 8 969/970. [7] J. J. Kovacic. An algorithm for solving second order linear homogeneous differential equations. J. Symbolic Comput. (): 4 986. [8] A. J. Maciejewski J.-M. Strelcyn and M. Szydłowski. Non-integrability of Bianchi VIII Hamiltonian System. J. Math. Phys. 4(4):78 74 00. [9] A. J. Maciejewski and M. Przybylska. Non-integrability of ABC flow. Phys. Lett. A 0(4):65 7 00. [0] A. R. Magid. Lectures on differential Galois theory volume 7 of University Lecture Series. American Mathematical Society Providence RI 994. [] J. J. Morales Ruiz. Differential Galois theory and non-integrability of Hamiltonian systems volume 79 of Progress in Mathematics. Birkhäuser Verlag Basel 999. [] J. J. Morales-Ruiz and J.-P. Ramis. Galoisian obstructions to integrability of Hamiltonian systems. I. Methods Appl. Anal. 8(): 95 00. [] V. N. Salnikov. On the dynamics of the triple pendulum: non-integrability topological properties of the phase space. 006. Lecture notes of conference Dynamical Integrability (CIRM) 006 pages. [4] M. van der Put and M. F. Singer. Galois theory of linear differential equations volume 8 of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag Berlin 00. [5] E. T. Whittaker and G. N. Watson. A Course of Modern Analysis. Cambridge University Press London 95.