A Criterion for Non-Complete Integrability of Hamiltonian Systems (original) (raw)

A Criterion for Non-Complete Integrability of Hamiltonian Systems

Delphine Boucher

Universit� de Limoges (France)

Algorithms Seminar

January 15, 2001

[summary by Philippe Dumas Bruno Salvy]

A properly typeset version of this document is available in postscript and in pdf.

If some fonts do not look right on your screen, this might be fixed by configuring your browser (see the documentation here).

Abstract

Finding polynomial solutions of linear differential equations is a building block implemented in several algorithms of computer algebra systems. In particular, this is a necessary sub-step when looking for rational, algebraic or Liouvillian solutions of linear differential equations. When there are no parameters, several algorithms are available, but the general case with parameters is undecidable. However, special families can be handled by ad hoc methods. Such methods were developed by Boucher who applied them to the nice example of integrability of the 3-body problem. The key idea there is to rely on a recent result of Morales-Ruiz and Ramis who relate complete integrability and differential Galois group. It turns out that special properties of this group can be related to computable properties of an appropriate linear differential equation, which leads Boucher to a ``simple'' sufficient condition for non-complete integrability.

1 Polynomial Solutions of Linear Differential Equations

The classical method to find polynomial solutions of linear differential equations over K(x), where K is a field, starts by determining a bound on the degree of potential solutions. This is a bound on the integer solutions of the_indicial equation_ at infinity.

Once a bound on the degree has been found, one uses an indeterminate coefficients method. The linear system on these coefficients has a band-matrix structure which can be exploited to accelerate the computation [1]. This linear system is rectangular, with more equations than unknown coefficients, thus existence of solution is related to the vanishing of a determinant.

When parameters occur in the equation (K is a field of rational functions), there are two difficulties: the size of the matrix may depend on the parameters and even when it does not, the determinant which must vanish is a polynomial in the parameters. Using Matijasevich's result on the undecidability of Hilbert's 10th problem (Is there a finite process which determines if a polynomial equation is solvable in integers?), it is possible to show that this problem itself is undecidable. More precisely, Jacques-Arthur Weil observes that the equation

has rational solutions if and only if P(_a_1,...,a m)=0 has integral solutions.

There are still cases where all polynomial solutions can be found: this happens when either the size of the matrix is bounded independently of the parameters and the vanishing of the required determinant can be determined or when the structure of the matrix is sufficiently regular to make the decision possible. Examples of both cases are given in [4].

2.1 Hamiltonian Mechanics

In the Hamiltonian approach to classical mechanics, the state of a system is characterized by 2_n_ variables, q i (positions) and p i (momenta), _i_=1,...,n, living in an open subset U_of **R**2_n (the phase space). More generally, the phase space of a system is the cotangent fibre bundle T*M of an _n_-dimensional real manifold M. The formulae we give below are expressions in a chart of useful quantities. The state variables satisfy where a dot denotes a derivative with respect to time and H(p,q,t) is the Hamiltonian. Physically, the Hamiltonian often represents the energy of the system. The system (1) governs the evolution of the system (in the phase space U). Solutions g(t) of (1) are the trajectories of the system.

In a more abstract setting, **R**2_n_ is endowed with a non-degenerate 2-form

known as Liouville's symplectic 2-form. Since w is non-degenerate, it induces an isomorphism between **R**2_n_ and its dual under which -dH is the image of a vector field X H. In this language, the Hamiltonian system (1) reduces to
First integrals are functions F(p,q) that are constant along the solutions g(t). A necessary and sufficient condition is where {F,H} is known as the Poisson bracket of _F_and H. In particular, the Hamiltonian itself is a first integral.

Two first integrals are in involution if their Poisson bracket vanishes. A Hamiltonian system is completely integrable when it possesses a set of n first integrals in involution that are independent (i.e., their Jacobian matrix is regular in the open set U).

Informally, a completely integrable system can be ``solved'' in terms of its first integrals. Indeed, given a first integral, a process known as symplectic reduction makes it possible to reduce the number of degrees of freedom by 1, i.e., the dimension by 2 [2, p. 91].

2.2 Many-Body Problem

In the many-body problem, n particles obeying Newton's law are governed by the following Hamiltonian: Note that here each p i and q i has coordinates in **R**3, thus the phase space has dimension 6_n_.

Apart from the Hamiltonian itself, known first integrals for this system are the momentum of the centre of mass and the angular momentum

�q i_� p i. Thus, the number of degrees of freedom can be reduced from 3_n to 3_n_-6 (or from 2_n_ to 2_n_-4 in the planar case).

For the 3-body problem, Poincar� proved that there are no other _complex analytic_first integrals. Bruns proved a similar result for complex algebraic first integrals.

2.3 Theorem of Morales-Ruiz and Ramis

We now present a simple version of a result of Morales-Ruiz and Ramis in [10, 11, 12, 13] (see also [3]) on non-complete integrability in terms of meromorphic first integrals. The Hamiltonian is analytic over an open set of **C**2_n_ and t (time) is a complex variable. Given a non-stationary trajectory G(t), following an idea of Poincar�, one considers the linear differential equation that must satisfy a ``small'' variation h, such that G(t)+h(t) is solution of the Hamiltonian system. This equation is called the variational equation along G. A theorem of Morales-Ruiz and Ramis relates complete integrability and Galois group of this equation. (For an introduction to differential Galois theory, see [14] or the summary of Ulmer's talk in this seminar in 1994.) However, since the Galois group is often very difficult to compute, it is useful to consider a differential equation of lower order. This is achieved by the following result.

Theorem 1 [Morales-Ruiz and Ramis] If the system possesses n meromorphic first integrals in the neighbourhood of G, independent and in involution, then the connected component of identity in the differential Galois group of the normal variational equation along G is abelian.

Similar earlier results of Ziglin based on the monodromy group and of Churchill, Singer et alii based on the Galois group did not extend to the case where the variational equation has an irregular singular point. In this theorem, the normal variational equation is an equation obtained from the variational equation through symplectic reduction. Indeed, dH(G(t))�h is a first integral of the variational equation, as can be seen by a first-order expansion.

3 Boucher's Criterion and its Application

It is not necessary to compute the Galois group of a linear differential equation in order to detect that it is not abelian. Thanks to a sufficient criterion [5, 6], Boucher has proved that the planar 3-body problem is not completely integrable in terms of meromorphic first integrals. Unfortunately, the formulae involved in this derivation are much too large to be reproduced here. Thus we content ourselves with a sketch of the steps and a description of the tools used in the calculations.

3.1 Criterion

Theorem 2 Assume that the linear differential operator L can be factored as KM, with M=lcm(L 1 ,...,L m ) where the L i , i=1,...,m, are irreducible (and lcm denotes the least common left multiple). Assume moreover that M(y)=0 has a formal solution with a logarithm. Then the connected component of the differential Galois group of L(y)=0 is not abelian.

Given a linear differential equation, this theorem reduces the task to factoring and finding formal solutions. Factoring can be done by an algorithm of van Hoeij [19, 20], and formal solutions can be computed at any singularity, including infinity [15, 20].

3.2 Application to the 3-Body Problem

Tsygvintsev and Boucher have proved independently that the planar 3-body problem is not completely integrable in terms of meromorphic first integrals. Their approaches [5, 17] follow the same initial steps till the normal variational equation. Then [17] uses Ziglin's result. We now outline Boucher's approach.

Reduced Hamiltonian

Using the first integrals obtained in Section 2.2, the problem is reduced to a Hamiltonian with three degrees of freedom, given in [17]. The parameters in this equation are the three masses _m_1, _m_2, _m_3and the value c of the angular momentum (which reduces to a scalar in this dimension). By homogeneity, we can freely assume _m_3=1. (Note that these transformations make the resulting expressions asymmetric with respect to the bodies.)

In order to apply Theorem 1, we need a particular solution of the system. This is provided by the celebrated Lagrange solutions. In these solutions, the three particles have orbits on similar conics with a common focus located at their centre of mass (see [8, p. 400]). Since any particular solution can be chosen, Tsygvintsev and Boucher concentrate on the parabolic orbit (for angular momentum c

�0).

Variational Equation

The variational equation (2) is a linear system of order _n_=6. The normal variational equation is obtained via a linear change of symplectic basis as follows. We observe that X H itself is a solution of the variational equation. It will be the first vector _e_1 of the new basis. Next, we compute a basis (_e_1=X,_e_2,...,e n,e n+2,...,e_2_n) of the kernel of dH(G(t)) satisfying w(e i,e n+i)=1 for 1<_i_� n and w(e i,e j)=0 otherwise. Finally, we compute a vector e n+1=Y such that w(e i,Y)=0 for _i_�1 and w(X,Y)=1. In the new basis (_e_1,...,e_2_n), the first column of the matrix of the variational equation is 0, since X H is a solution. Now, for any vector field h,w(X,h)=-dH(G(t))�h, therefore for any solution h of the variational equation, the value of this first integral is the coordinate of h on the vector Y in the new basis. The normal variational equation is obtained by setting this coordinate to 0 and considering the induced matrix A on the subspace with basis (_e_2,...,e n,e n+2,...,e_2_n).

Cyclic Vector

The criterion of Theorem 1 applies to equations rather than systems. A classical method to convert a system of order m into an equation L(u)=0 is to start from a random vector u and find a linear dependency between the m+1 vectors u,u',...,u(m) where the derivatives are computed using the matrix A. Unfortunately, this process generically introduces spurious singularities that are roots of the determinant of the change of basis (u,u',...,u(_m_-1)). Boucher therefore selects cyclic vectors in such a way that no new singularity occurs and this requires distinguishing two cases depending on the value of the mass _m_1.

Right Factors

In the simplest case of Boucher's criterion, the operator L has an irreducible right factor M whose formal solutions exhibit logarithms. This requires M to have order at least 2. Factors of order k are found by constructing an auxiliary equation _L_� k of order () whose solutions are Wronskians of k independent solutions of L [7]. (Note that this can be computed directly from L.) Indeed, a monic right factor of order k has for coefficient of order _k_-1 the logarithmic derivative w'/w of some particular Wronskian of its solutions. Finding right factors then amounts to looking for so-called exponential solutions of_L_� k (i.e., those with logarithmic derivative that is rational). From a basis of such solutions, corresponding to linear combinations of Wronskians, Pl�cker's relations help select those that are indeed Wronskians [16]. From there, the complete factor can be reconstructed. Exponential solutions are found by looking at formal solutions at all singularities of the equation [19]. This requires a discussion in the parametric case. If a factor is found, the next step is to check whether this factor is irreducible, or to find conditions on the parameters that make it irreducible. This is done again by searching for factors of the factor. It turns out that in this application, in all cases an irreducible right factor of order 2 is found.

Logarithms

Logarithms in formal solutions occur when the indicial equation at a singularity has roots that differ by an integer. A necessary and sufficient condition has been given by Frobenius [9, p. 404--406]. Again, in all generality nothing can be said when parameters are present but Boucher manages to show that logarithms are present in all cases for this application.

4 Conclusion

This application is a very good showcase for many of the algorithms that have been developed in computer algebra for linear differential equations: formal solutions, factorization, polynomial solutions, ...

What Boucher has shown is that, even in the presence of parameters, these algorithms can be exploited to provide useful information by concentrating on those points where specific quantities such as the indicial equation or its solutions do not depend ``too much'' on the parameters.

A recent trend in computer algebra is to revisit all these algorithms that have been designed for equations and extend them to deal with systems, without using the cyclic vector. It would be a natural step to try and adapt Boucher's criterion so that the symplectic structure is not lost. (Work on this has been started by Boucher and Weil.)

Remark 1 _A new result of Tsygvintsev [_18] shows the stronger result that there is no additional meromorphic first integral. Also, Theorem 2 has been extended to the case when L is a product of irreducible factors one of which has a solution with logarithms.

References

[1]

Abramov (Sergei A.), Bronstein (Manuel), and Petkovsek (Marko). -- On polynomial solutions of linear operator equations. In Levelt (A. H. M.) (editor), Symbolic and Algebraic Computation. pp. 290--296. -- ACM Press, New York, 1995. Proceedings of ISSAC'95, July 1995, Montreal, Canada.

[2]

Arnol'd (V. I.), Kozlov (V. V.), and Neishtadt (A. I.). --Dynamical systems. III. -- Springer-Verlag, Berlin, 1988, xiv+291p. Mathematical Aspects of Classical and Celestial Mechanics. Translated from the Russian.

[3]

Audin (Mich�le). -- Int�grabilit� et non-int�grabilit� de syst�mes hamiltoniens [d'apr�s S. Ziglin, J. Morales-Ruiz, J.-P. Ramis, ...]. -- To appear in Ast�risque. Bourbaki seminar 884, March 2001. Preliminary version available athttp://irma.u-strasbg.fr/~maudin/publications.html.

[4]

Boucher (Delphine). -- About the polynomial solutions of homogeneous linear differential equations depending on parameters. In Dooley (Sam) (editor), Proceedings of the 1999 International Symposium on Symbolic and Algebraic Computation (Vancouver, BC). pp. 261--268. -- ACM, New York, 1999.

[5]

Boucher (Delphine). -- Sur la non-int�grabilit� du probl�me plan des trois corps de masses �gales. Comptes Rendus de l'Acad�mie des Sciences. S�rie I. Math�matique, vol. 331, n°5, 2000, pp. 391--394.

[6]

Boucher (Delphine). --Sur les �quations diff�rentielles lin�aires param�tr�es, une application aux syst�mes hamiltoniens. -- PhD Thesis, Universit� de Limoges, October 2000.

[7]

Bronshtein (M.) and Petkovshek (M.). -- Ore rings, linear operators and factorization. Rossiiskaya Akademiya Nauk. Programmirovanie, vol. 1994, n°1, 1994, pp. 27--44.

[8]

Hestenes (David). --New foundations for classical mechanics. -- D. Reidel Publishing Co., Dordrecht, 1986, xii+644p.

[9]

Ince (E. L.). --Ordinary differential equation. -- Dover Publications, New York, 1956, viii+558p. Reprint of the 1926 edition.

[10]

Morales-Ruiz (Juan J.) and Ramis (Jean Pierre). -- Galoisian obstructions to integrability of Hamiltonian systems: statements and examples. In Hamiltonian systems with three or more degrees of freedom (S'Agar�, 1995), pp. 509--513. -- Kluwer Acad. Publ., Dordrecht, 1999.

[11]

Morales-Ruiz (Juan J.) and Ramis (Jean-Pierre). -- Galoisian obstructions to integrability of Hamiltonian systems.Methods and Applications of Analysis, vol. 8, n°1, 2001.

[12]

Morales-Ruiz (Juan J.) and Ramis (Jean-Pierre). -- Galoisian obstructions to integrability of Hamiltonian systems. II. Methods and Applications of Analysis, vol. 8, n°1, 2001.

[13]

Morales-Ruiz (Juan J.) and Ramis (Jean-Pierre). -- A note on the non-integrability of some Hamiltonian systems with a homogeneous potential. Methods and Applications of Analysis, vol. 8, n°1, 2001.

[14]

Singer (Michael F.). -- An outline of differential Galois theory. In Computer algebra and differential equations, pp. 3--57. -- Academic Press, London, 1990.

[15]

Tournier (�velyne). --Solutions formelles d'�quations diff�rentielles. -- Doctorat d'�tat, Universit� scientifique, technologique et m�dicale de Grenoble, 1987.

[16]

Tsar�v (S. P.). -- Some problems that arise in the factorization of linear ordinary differential operators. Rossiiskaya Akademiya Nauk. Programmirovanie, n°1, 1994, pp. 45--48. -- (Russian). Translation in Programming and Computer Software, vol. 20, n°1, 1994, pp. 27--29.

[17]

Tsygvintsev (Alexei). -- La non-int�grabilit� m�romorphe du probl�me plan des trois corps. Comptes Rendus de l'Acad�mie des Sciences. S�rie I. Math�matique, vol. 331, n°3, 2000, pp. 241--244.

[18]

Tsygvintsev (Alexei). -- Sur l'absence d'une int�grale premi�re m�romorphe suppl�mentaire dans le probl�me plan des trois corps. Comptes Rendus de l'Acad�mie des Sciences. S�rie I. Math�matique, vol. 333, n°2, 2001, pp. 125--128.

[19]

van Hoeij (Mark). -- Factorization of differential operators with rational functions coefficients. Journal of Symbolic Computation, vol. 24, n°5, 1997, pp. 537--561.

[20]

van Hoeij (Mark). -- Formal solutions and factorization of differential operators with power series coefficients. Journal of Symbolic Computation, vol. 24, n°1, 1997, pp. 1--30.

This document was translated from LATEX byHEVEA.