Open access peer-reviewed chapter

# Topology and Integrability in Lagrangian Mechanics

By Leo T. Butler

Submitted: May 10th 2016Reviewed: October 4th 2016Published: May 3rd 2017

DOI: 10.5772/66147

## Abstract

This chapter reviews complete integrability in the setting of Lagrangian/Hamiltonian mechanics. It includes the construction of angle-action variables in illustrative examples, along with a proof of the Liouville-Arnol’d theorem. Results on the topology of the configuration space of a mechanical (or Tonelli) Hamiltonian are reviewed and several open problems are high-lighted.

### Keywords

• Hamiltonian mechanics
• Lagrangian mechanics
• integrability
• topological obstructions
• topological entropy

## 1. Introduction

Lagrangian mechanics employs the least-action principle to derive Newton’s equations from a scalar function, the action function L. In classical mechanics, Lis the difference of kinetic and potential energies and therefore appears as an artifice. It is somewhat mysterious, then, that the reformulation of Newtonian mechanics in terms of momentum and position, rather than velocity and position as in Lagrangian mechanics, leads immediately to the total energy function Hand a plethora of geometric structure that is hidden in the native setting.

Due to the advantages of the Hamiltonian perspective, this chapter studies Lagrangian systems from this dual point of view. The organization of the chapter is this: Section 2 recalls the classic construction of angle-action variables in 1degree of freedom via several examples, then states and proves the Liouville- Arnol’d theorem; Section 3 discusses the relationship between the topology of the configuration space and the existence of integrable mechanical systems; and it reviews several constructions of integrable systems whose configuration space is the sphere or torus. Section 3 provides a number of open problems that may stimulate interested researchers or students.

## 2. Integrability in Hamiltonian mechanics

### 2.1. Integrability in 1 degree of freedom

One of the central problems in classical mechanics is the integrability of the equations of motion. The classical notion of integrability is loosely related to exact solvability, and roughly corresponds to the ability to solve a system of differential equations by means of a finite number of integration steps.

2.1a. Example: Harmonic oscillator Let us take the simple harmonic oscillator, or an idealized Hookean spring-mass system, with mass mand spring constant k. If qis the displacement from equilibrium and pthe momentum, then the total energy is

H=12mp2+k2q2,andequationsofmotionare{q˙=p/m,p˙=kq}.E1

The change of variables (q,p)=(Q/λ,λP)transforms the system to, with λ=km4,

H=ω2(P2+Q2),andequationsofmotionare{Q˙=ωP,P˙=ωQ},E2

where ω=k/m. A second change of variables (Q,P)=(2Icosθ¸2Isinθ¸)transforms the system to

H=ωI,andequationsofmotionare{θ˙=ω,I˙=0}.E3

The differential equations in (3) are trivial to integrate since the right-hand sides are constants. Let us explain the sequence of transformations. The change of coordinates (q,p)(Q,P)is an area-preserving linear transformation that transforms the elliptical level sets of Hinto circles. The transformation (Q,P)(θ,I)is analogous to the introduction of polar coordinates–indeed the transformation (r,θ)=(2I,θ)is a transformation to polar coordinates. Because the area form dPdQ=rdrdθ, we see that the transformation dPdQ=dI=dθ.

Therefore, the change of coordinates (q,p)(θ¸I)not only reveals the exact solutions of the harmonic oscillator equations, it is area preserving.

Suppose that for some reason one did not know to introduce “polar” coordinates. One might still determine the change of coordinates using only that the transformation (Q,P)(θ¸I)preserves area. Indeed, since d(PdQIdθ)=0, there is a function ν=ν(Q,θ)such that PdQIdθ=dνor P=vQand I=vθ. Then, upon substituting the identity P=νQinto (2), one obtains

ν0Q2H/ωQ2dQ=12Q2H/ωQ2(H/ω)arccos(Q/2H/ω),E4

where indicates that νequals the right-hand side up to the addition of a 2π-periodic function of θ.

If (Q,P)make a complete circuit around the contour {H=c}then one obtains from (4) and the identity that P=νQthat

Δν={H=c}PdQ=(H/ω)2π.E5

On the other hand, since d2=0and Iis held constant on the contour, Green’s theorem implies that

Δν={H=c}dν+Idθ={H=c}Idθ=2πI.E6

Equating (5) and (6) shows that H=ωI.

These calculations show that one may determine Has a function of Iwithout explicit knowledge of the coordinate transformation (Q,P)(θ,I)–but one does need to solve the Hamilton-Jacobi equation

H(Q,νQ)=c,E7

for ν, as performed in Eq. (4). At this point, if one wants to derive the change of coordinates from ν, Eq. (4) shows that it is easier to write ν=ν(Q,I), in which case PdQ+θdI=dνor

θ=νI=arccos(Q/2I),E8

so Q=2Icos(θ)and P=νQ=2Isin(θ).

Let it be observed that if, in Eq. (4), one had chosen the anti-derivative to be arcsinrather than arccos, then Qwould be 2Isin(θ)and Pwould be ±2Icos(θ). However, because dPdQ=dIdθ, one would be obligated to choose the negative square root to define P; otherwise, dPdQ=dIdθ.

2.1b. Example: the planar pendulum. Let us take the idealized planar pendulum with a mass-less rigid rod of length lsuspended at a fixed end with a bob of mass mat the opposite end (Figure 1). The total energy is

H=12mp2+mlg(1cosq),and{q˙=p/m,p˙=mlgsinq}.E9

To simplify the exposition, assume that the mass m=1and let ω2=16lg, where ωis 4times the frequency of the linearized oscillations at q=p=0. The substitution q=2Q, p=P/2transforms the Hamiltonian to

8H=P2+ω2sin2Q,and{Q˙=P/4P˙=ω2sin(2Q)/8}.E10

If one tries to solve for a generating function ν=ν(Q, I)of a coordinate change (Q,P)(θ,I)such that H=H(I), then one obtains from P=νQthat

ν0Q8Hω2sinQdQ=ωkE(Q,k)E11

where indicates equality up to a 2π-periodic function of θ, H=ω2/(8k2)and Eis the elliptic integral of the second kind defined by E(x,k)=0x1k2sin2x dx.

If (Q,P)make a complete circuit around the contour {H=c}, then one obtains from Eq. (11) that

Δν=4ωkK(k)E12

where K(k)=E(Q+(k),k)and Q+(k)=arcsin(1/k)if k>1and π/2if k<1(in which case, Kis the complete elliptic integral of the second kind). The area of the shaded region Kin Figure 2 shows the geometric meaning of K(k)for k>1. Along with the identity (6), one obtains

I=2πωkK(k),E13

which determines H=H(I)implicitly.

Figure 3 graphs Has a function of Iusing the definition of Iin (13) with ω=1, along with the graph for the harmonic oscillator. Although Happears to be a smooth function of Ion the interval depicted, this is a numerical artifact. Indeed, there are two distinct proofs that Hcannot be differentiable in Iover the interval [0,1]. Without loss of generality, it is assumed that ω=1.

The first, calculus-based, proof is this: as k1+(H1/8),I/k. If His a differentiable function of I, then H/I=0at I=2/π (H=1/8). But then the entire level set consists of fixed points, which is false.

The second, topological, proof is this: each level set {H=c}, c<1/8, is connected; each level set for c>1/8has exactly two connected components (c.f. Figure 2). If the generating function v were differentiable in (Q,I)on any rectangle containing R/πZ×{2/π}, then Eq. (13) would determine a homeomorphism H=H(I), and so the level sets of Hwould remain connected on either side of the critical level at height 1/8. Absurd.

To derive the change of coordinates (Q,P)(θ,I)from the generating function ν, one uses the identity θ=νIand properties of the elliptic integrals to deduce

θ=π2F(Q,k)(Q+(k),k)Q=amk(2F+πθ)E14

where F(x,k)=0xdx/1k2sin2xis the elliptic integral of the first kind, F+=F(Q+(k), k)and amk(u)is the Jacobian amplitude function, a local inverse to F([1], Chapter 2). Along with P=νQ, (14) implies that

P=ωkdnk(2F+πθ),E15

where dnk(u)is the Jacobian elliptic function.

2.1c. Example: a mechanical system. Let V=V(Q)be a smooth potential function of a 1-degree-of-freedom Hamiltonian system with

H=12P2+V(Q).E16

If one attempts to find the generating function ν=ν(Q,I)of an area-preserving transformation (Q,P)(θ,I)that transforms H=H(I), then one deduces that

νQ0Q2(HV(Q))dQ,E17

up to a function depending only on I. Then, in a complete circuit around the connected contour {H=c}, one has 2π=Δθ=ΔνIidentically, so

2πI={H=c}PdQ.E18

and, upon solving (18) for H=H(I), one inverts

θ=12Q0QHIHV(Q)dQ,E19

to obtain Q=Q(θ,I), and finally P=νQyields P=P(θ, I). Since the change of coordinates is area-preserving, the Hamiltonian form of the equations of motion are preserved, so the resulting equations are

H=H(I)and{θ˙=H/II˙=0}.E20

### 2.2. The generating function

The above three examples use a generating function ν=ν(Q,I)of a mixed system of coordinates in order to create an area-preserving change of coordinates to angle-action variables (θ,I).

2.2a. Question: why do the angle-action variables exist? In order to understand the generating function, it is necessary to clarify the existence of the coordinates (θ,I), which are commonly called angle-action variables. Let H:XRbe a smooth function from an oriented surface Xto the reals. If it is assumed that AXis an open, connected, saturated (H1(H (A))A=A)subset of the domain of H, H|Ahas no critical points and H|Ais proper, then H|Ais a submersion onto the interval B=H(A)R. Since H|Ais proper, for each bB, the level set (H|A)1(b)is a compact one-manifold and hence its components are circles. Since A is connected and H|Ais critical-point free, the level set must be connected, so it is a circle. Therefore, the submersion theorem implies that Ais diffeomorphic to A= S1×B.[1] - To make this system of coordinates concrete, note that there is a complete vector field Uon Asuch that dH(U)1. Let γAbe a segment of an integral curve of Uwhich is maximal (i.e. an integral curve that strictly contains γintersects XA). For each aA, let t=t(a)be the time along the flow line of the Hamiltonian vector field XHbeginning at the initial condition γH1(H(a)). The function tis multi-valued, since the flow line is closed, so it should be considered as a function on the universal cover of A.

Since the tangent space at aAis spanned by XHand U,Ωis determined by Ω(XH,U). But Q(XH,U)=dH(U)=1, so Ω=dH dt.

Let Tbe such that 2πTis the least period of the function t(i.e. the first return time to γ). Then T=T(H)is a function of Halone. Define θby

θ= t/T(H)(mod 2π)and I bydI=T(H)dH.E21

The function θis the normalized time along the flow lines of Hamiltonian vector field XH, while dH/dI=1/T(I)is the frequency. One computes that the oriented area form Ω=dH dt=dI dθ. Moreover, in the coordinates (θ,I), the Hamiltonian vector field

XH={θ˙=1/T(H)=dHdI,I˙=0}.E22

This proves the existence of an area-preserving diffeomorphism ϕ:D×S1A, where DRis an open interval, such that the Hamiltonian His transformed to a function of Ialone; and ϕis as smooth as Hand the area form Ωare (e.g. if both are real-analytic, then ϕis real-analytic).

2.2b. Question: what kind of “function” is ν? In the first instance, νis not single-valued. Indeed, one postulates the area-preserving change of coordinates ϕ:(Q,P)(θ,I)to deduce that

d(PdQ+θdI)=0,E23

so that locally there is a function νsuch that

PdQ+θdI=dν.E24

But since θis an angle variable, this equation can only hold globally modulo 2πZdI. So, in this formulation of the generating function, νcan only be defined globally modulo 2πZI. Or, equivalently, νIis a function with values in the circle R/2πZ.

The way to resolve these ambiguities or difficulties is simple: the domain of the change of coordinates ϕmust be non-simply connected (a disjoint union of open annuli, in fact, as can be deduced from the discussion above) and so one should view (24) as holding globally on the universal cover of this annulus where θis a single-valued real function (c.f. 21). In this case, the lift of a closed contour {H=c}is a path that projects to the contour and whose endpoints differ by a deck transformation–which in the angle-action coordinates is (θ,I)(θ+2π,I). Since Iis constant along this path, the path integral of P dQequals the path integral of dν, i.e. Δν, the change in νfrom one preimage to its translate. With this understanding, Eq. (18) is correct. And, indeed, one sees that the integral in Eq. (17) is defined not on the domain of the coordinate change ϕbut on its universal cover; the same is true for the integral in Eq. (19), but the marvellous fact about that integral is that it is 2π-periodic: this follows from the observation that Δθ=2πidentically around a closed connected contour in {H=c}.

So to answer the question that started the section, the generating function νis a function defined on the universal cover of the union of regular compact levels of Hwhich implicitly defines a 2π-periodic change of coordinates to “angle-action” variables (θ,I).

### 2.3. Integrability in 2 or more degrees of freedom and Tonelli Hamiltonians

Integrability in 2 or more degrees of freedom is substantially more involved than the case of 1 degrees of freedom. Of course, a sum of n distinct, non-interacting 1-degree-of-freedom Hamiltonians is a simple case; and upon reflection, a not-so- simple case, because this condition is not coordinate independent. Indeed, a necessary and sufficient condition is that the Hamiltonian vector field be Hamiltonian with respect to two distinct non-degenerate Poisson brackets {,}ithat are compatible in the sense that the linear space spanned by the brackets is a space of Poisson brackets, and maximal in the sense that a “recursion” operator naturally defined from the two brackets has a maximal number of functionally independent eigenvalue fields [2].

Let us turn now to a definition which generalizes mechanical Hamiltonians.

Definition 2.1 (Tonelli Hamiltonian). Let Σbe a smooth n-manifold and T*Σits cotangent bundle. A smooth function H:T*ΣRwhich satisfies (T1) H|Tx*Σis strictly convex for each xΣ; and (T2)H(x,tp)/tuniformly as t, is called a Tonelli Hamiltonian.

As noted, Tonelli Hamiltonians are natural generalizations of mechanical systems. For this reason, Σwill be referred as the configuration space of the Hamiltonian H.

If Qiare coordinates on Σand Θ=iPidQiare the coordinates of the 1-form Θ, then the canonical symplectic structure Ω=dΘ=dPidQion T*Σ. The symplectic form Ωequips the space of smooth functions on T*Σwith a Poisson bracket denoted {,} that satisfies

{Pi,Qj}={Qj,Pi}=δij{Qi,Qj}={Pi,Pi}=0E25

for all i,j. The Poisson bracket is fundamental to Hamiltonian mechanics. For each smooth function H, one has a smooth vector field XH={H,}, and the skew symmetry of the bracket implies that His preserved by the flow. One says that H1and  H2Poisson commute if {H1,H2}0.

A fundamental result in Hamiltonian mechanics is the Liouville-Arnol’d theorem, which provides a semi-local description of a completely integrable Hamiltonian and the Poisson bracket.

Theorem 2.1 (Liouville-Arnol’d). Let H:T*ΣRbe a smooth Hamiltonian. Assume there exists n functionally independent, Poisson commuting conserved quantities F=(F1=H,,Fn):T*ΣRn. If LF1(c)is a compact component of a regular level set, then there is a neighbourhood Wof Land a diffeomorphism ϕ=(θ,I):Tn×BnWsuch that

F=F(I){Ii,θj}=δij,{Ii,Ij}={θi,θj}=0,XFi=Fi(I)Ijθj,

that maps L to Tn×{0}.

In such a situation, it is said that His Liouville, or completely, integrable. The torus Tn×{I0}is a Liouville torus, the neighbourhood Tn×Bnis a toroidal ball and the conserved quantities are first integrals. Systems with k first integrals, of which l<kPoisson commute with all kfirst integrals, where k+l=2nare called non-commutatively integrable; when k=2n1, the system is also called super-integrable c.f. [3, 4].

There are several proofs of the Liouville-Arnol’d theorem in the literature. The basic ideas are already captured in the one-dimensional case discussed in Section 2.2.

It can be assumed, without loss, that L=F1(c). Since cRnis a regular value of F, the submersion theorem implies that there is an open neighbourhood Cof cconsisting of regular values of Fand the open set F1(C)is diffeomorphic to L×C. Therefore, there is a smooth n-dimensional submanifold MF1(C)such that Mtransversely intersects each level set Lf=F1(f),fC. Possibly by shrinking the open set C, it can be assumed that Mis Lagrangian: Ω|M0.[1] -

Because the functions Fi,,FnPoisson commute and are functionally independent, the Hamiltonian vector fields XF1,,XFnspan the tangent space TxLf, for each xLf, fC. Because Lfis compact, each vector field is complete, so there is a well-defined flow map ϕFi:R×F1(C)F1(C). Because Fi,,FnPoisson commute, the respective flow maps commute, so there is an action of Rnon F1(C)defined by

ϕt=ϕFntnϕF1t1E26

for all tRn. Define a map

Φ(t,m)=ϕt(m),tRn,mM.E27

This is a smooth map which is a local diffeomorphism of Rn× Mwith F1(C). Indeed, ϕtcarries each level Lfinto itself and carries Minto a submanifold ϕt(M)transverse to Lfat ϕt(m); on the other hand, the derivative of ϕtwith respect to tis a surjective linear map onto Tϕt(m)Lf. Therefore, dΦis surjective, so injective, hence Φis a local diffeomorphism onto its image. Compactness and connectedness of the levels Lfimply that the image of Φis F1(C).

For each mM, let P(m)Rnbe the set of t such that Φ(t,m)=m. Since each level set is compact, P(m)is a discrete subgroup of Rnisomorphic to Zn. This is the “period lattice” of the action ϕ. If one selects a basis of P(m), one obtains a map MGL(n;R),m2πT(m). The implicit function theorem implies that there is a smooth map amongst these maps. Moreover, since F|Mis a bijection onto its image, one can take the components of Fas coordinates on M, or in other words, T=T(F).

Define functions θ=(θ1,,θn)by

θ=T(F)1t(mod2π),θ:Rn×CRn/2πZn.E28

The flow map Φtherefore induces a diffeomorphism F1(C)Tn×C:x(θ(x),F(x)).

To complete the proof, one might show that each vector field /θiis Hamiltonian with Hamiltonian function Iiand that Fis functionally dependent on Iso that (θ,I)is a canonical system of coordinates on F1(C). This is performed indirectly. Define the functions Ii=Ii(F)by

2πIi=Γi(F)ξ,E29

where ξ=PdQis the primitive of the symplectic form Ωand Γi(F)is the cycle on LFon which θiincreases from 0 to 2πand the other angle variables are held equal to 0. To show that (θ,I)is a system of coordinates on F1(C), one computes the Jacobian [Ii/Fj]:

2πIiFj=lims01sCj(F,s)Ω,E30

where, in the (t,F)coordinate system,

Cj(F,s)={(uT(F+υej)i,F+υej)|u[0,2π],υ[0,s]}

is the “cylinder” obtained by sweeping out the cycles Γj(F+vej)as the j-th component of Fincreases from Fjto Fj+s, and Tiis the i-th column of the period matrix T. Since

Ω(Fj,tk)=FkFj=δjk,E31

which implies

IiFj=Tji.E32

Since the period matrix Tis non-singular, the transformation (θ,F)(θ,I)is a diffeomorphism.

Finally, the functions I1,,InPoisson commute and since Mis Lagrangian, the functions ti,,tnPoisson commute, which implies θ1,,θnPoisson commute. And, since {Fi,tj}=δij, this implies that {Ii,θj}=δij.

The remainder of the theorem follows from the fact that the angle-action coordinates (θ,I)are canonical and F=F(I).

## 3. Topology of configuration spaces

The central problem in the theory of completely integrable Tonelli Hamiltonian systems is to

Problem 3.1. Determine necessary conditions on the configuration space Σ for the existence of a completely integrable Tonelli Hamiltonian H.

This is a broad, overarching problem which has motivated research by many authors over an almost 40-year period, including many of the author’s publications. It is helpful to pose several sub-problems which address aspects of this problem and that appear to be amenable to solution. The remainder of this section is devoted to an elaboration of this problem, along with known results. We start with two-dimensional configuration spaces.

### 3.1. Surfaces of genus more than one

As a rule, completely integrable Tonelli Hamiltonians are quite rare, as are the configuration spaces Σwhich support such Hamiltonians. Indeed, in two dimensions, the compact surfaces that are known to support a completely integrable Tonelli Hamiltonian are the 2-sphere, S2, the 2-torus T2and their non-orientable counterparts. With some quite mild restrictions on the singular set–called condition , and assuming that the Hamiltonian is Riemannian, Bialy has proven these are the only compact examples [5]. This extended an earlier result of V. V. Kozlov [6]; the author has obtained a similar result for super-integrable Tonelli Hamiltonians [7].

V. Bangert has suggested to the author that Bialy’s argument should extend to prove the non-existence of a C2integral that is independent of the Hamiltonian when Σis a compact surface of negative Euler characteristic (c.f. [8]). The idea of such a proof would be the following (assuming that His Riemannian): Suppose that Henjoys a C2integral Fthat is independent on a dense set, hence that the union of Liouville tori is dense. Let ΓH1(12)be the union of orbits which project to minimizing geodesics. It is known, due to results of Manning and Katok [9, 10], that Γ¯contains a hyperbolic invariant set Λon which the flow is conjugate to a horseshoe. Let λΛbe a closed orbit of the geodesic flow of period T. Since the union of Liouville tori is dense, for each ϵ>0, there is a Liouville torus Lλ,ϵthat contains an orbit of the geodesic flow that remains within a distance ϵof λover the interval [0,T]. Hence, π1(Lλ,ϵ)has a homotopy class mapping onto λ. Since λis minimizing, it has no conjugate points and so for ϵsufficiently small, the same is true for the orbit on Lλ,ϵover the time interval [0,T]. This implies that the image of π1(Lλ,ϵ)is (free) cyclic and the kernel is generated by a cycle that bounds a disc–in classical terminology, this means that Lλ,ϵis compressible. It follows that Lλ,ϵbounds a solid torus TλT1×B2that is invariant for the geodesic flow. The integral F|Tλinduces a singular fibration of the solid torus by invariant 2-tori.

Thus, for each closed orbit λin the hyperbolic invariant set Λ, we have produced an invariant solid torus Tλthat shadows λ–at least in some rough, homotopic sense. This fact alone should suffice to achieve a contradiction.

Problem 3.2. Let Σbe a compact surface of negative Euler characteristic. Extend the above argument to prove the non-existence of a smooth Tonelli Hamiltonian H:T*ΣRwith a second C2integral Fthat is independent on a dense set; or give an example of a completely integrable Tonelli Hamiltonian H:T*ΣR.

V. Bangert proposes similar problems in his contribution in ([8], Problems 1.1, 1.2).

There is a similar, but possibly more accessible, problem for twist maps. Recall that if we discretize time, the notion of a Tonelli Hamiltonian is replaced by that of a twist map f:T*Σ→T*Σ which is a symplectomorphism that satisfies a condition analogous to T1. If fenjoys nindependent, Poisson commuting first integrals, then the Liouville-Arnol’d theorem implies that some power of facts a translation on the Liouville tori. We noted above that the Hamiltonian flow of a Tonelli Hamiltonian has a horseshoe on an energy level.

Problem 3.3. Let f:T*T1T*T1be a twist map. If f has a horseshoe and a C1first integral F, is Fnecessarily constant on an open set?

### 3.2. The 2-torus

Let us turn now to the torus. The 2-torus T2admits a family of completely integrable Riemannian Hamiltonians which are called Liouville. These are of the form

H=px2+py22[f(x)+g(y)]F=g(y)px2f(x)py2f(x)+g(y)E33

where f,g:T1Rare smooth positive functions and (x,y,px,py)is a canonical system of coordinates on T*T2. The degenerations of the Liouville family include the rotationally symmetric (fconst.)and flat (f,gconst.).

The Liouville family is obtained from two uncoupled mechanical oscillators with periodic potentials,

G=12(px2+py2)+a(x)+b(y),E34

on an energy level E=α+β>maxa+maxbsuch that f=αa, g=βb. The Maupertuis principle states that orbits of the Hamiltonian flow of Gon the energy level {G=E}are orbits of the Hamiltonian flow of Hup to a change in time along the orbit. The complete integrability of Gis explained in Sections 2.1c and 2.3.

It is a remarkable fact that the Liouville family exhausts the list of known completely integrable Riemannian Hamiltonians whose configuration space is T2. Indeed, in 1989, Fomenko conjectured that these are the only examples possible when the second integral in polynomial-in-momenta [11]. Most recently, in 2012, Kozlov, Denisova and Treschëv reiterate Fomenko’s conjecture ([12], p. 908).

Let us note that it is a well-known fact that, if the first integral Fis real- analytic, then F=N0FNwhere each term FNis polynomial-in-momenta with real-analytic coefficients, homogeneous and of degree Nand since {H, FN}is polynomial-in-momenta, homogeneous and of degree N+1, each graded piece of Fis a first integral. So, there is no loss in generality in restricting attention to polynomial-in-momenta first integrals–and, indeed, a slight increase in generality because the coefficients of the polynomial-in-momenta first integral are not assumed to be real-analytic.

In [13, 14], Kozlov and Denisova prove that if, when (x,y)are isothermal coordinates, and

H=12Λ(px2+py2),E35

with the conformal factor Λa trigonometric polynomial, then the existence of a second independent first integral that is polynomial-in-momenta implies that His Liouville.

In [12], Denisova, Kozlov and Treschëv prove that, if one only assumes Λis smooth, then Hhas no irreducible polynomial-in-momenta first integral Fthat is of degree 3 or 4 that is independent of H. Mironov separately proves the non-existence of Fof degree 5, but as noted in ([12], p. 909), Λsatisfies an extra unstated hypothesis [15]. The line of attack used in these papers is pioneered in [16], where Kozlov and Treschëv introduce the notion of the spectrum S⊂2nZ2 of the function Λas the support of the Fourier transform of Λ. This spectrum is finite iff Λis a trigonometric polynomial; Denisova and Kozlov prove that, in this case, any first integral of His dependent on Hunless the spectrum S is contained in a pair of orthogonal lines through (0, 0), in which case His Liouville and has a second independent first integral that is quadratic-in-momenta. Without the hypothesis that S is finite, the problem becomes significantly more delicate. The bulk of [12], for example, is devoted to a study of solutions to a PDE that characterizes the first integral Fby means of Fourier analysis.

An alternative approach, due to Bialy and Mironov, is to observe that the equation {H,F}=0coupled with the hypothesis that Fis polynomial-in-momenta of degree Nimplies that when we write Fas

F=j=0Naj(x,y)pxNjpyjE36

then the coefficients a0,,ansatisfy a semi-linear PDE [17, 18]. Indeed, there is a system of coordinates (τ,υ)on T2such that, when Fis written in the adapted canonical coordinates as F=j=0Nuj(τ,υ)(pυ/g)jpr(Nj)then this equation is of the form

uυ+T(u)ur=0E37

where u0=1, u1=g, u=(u1,u2,,uN)and

T(u)ij={ui+1ifj=i+1,(i+1)ui+1(N1i)ui1ifj=1,0otherwise,E38

where we adopt the convention that u_1= uN+10.

A standard technique to solve a quasi-linear PDE like (37) is to diagonalize it, that is, to find Riemann invariants, so that it is equivalent to

rυ+Δ(r)rr=0whereΔ(r)=diag(δ1(r),,δN(r)),r=(r1,,rN).E39

To find Riemann invariants, Bialy and Mironov employ the following trick: let pυ=gcos(θ),pr=sin(θ)parameterize cotangent fibres of H1(12). The invariance condition {H, F}=0translates to Fvg1cos(θ) +Fτsin(θ)=0along the locus where Fθ=0, i.e. where dFand dHare co-linear. If one supposes that θi=θi(τ,υ),i=1,,N, is a smooth parameterization of the critical-point set, then the critical values ri=F(r,υ,θi(r,υ))are Riemann invariants with δi=g(τ,v)×tan(θi). Of course, the main problem is to determine the relationship between the Liouville foliation–the singular foliation of T*T2by the Liouville tori and their degenerations–and the system 39.

In ([18], Theorems 1 and 2), Bialy and Mironov prove that if N4, then in any region where a multiplier δiis non-real, the metric is Liouville. One can view the result of Bialy and Mironov as a partial confirmation of Fomenko’s conjecture and an important step toward resolving that conjecture.

The key step in Bialy and Mironov’s proof is to show that, in any region where δiis non-real, the imaginary part of the Riemann invariant risatisfies an elliptic PDE. It appears that the properties of this PDE are key to proving stronger results.

Problem 3.4. Extend Bialy and Mironov’s work to show that there are no regions where any multiplier δiis non-real on T2, i.e. show that (39) is a hyperbolic system.

There is good reason to believe that the multipliers δiare always real. When δi1(C\R)T2, Bialy and Mironov prove that the Riemann invariant riis real and constant, say ri=si. This implies that the common level set F1(si)H1(12), a subset of the complexified cotangent bundle TC*T2, has a tangent with the fibres of TC*T2on an open set. That picture is dramatically at odds with the real picture, where the tangency can occur along a one-cycle at most. Because of this, it seems likely that there is a geometric proof of Problem 3.4.

Hyperbolicity of Eq. (39) has additional meaning. As the previous paragraph alluded to, the points where Fθ=0are the critical points of the canonical projection map π:T*T2T2restricted to a common level F1(r)H1(12). Such tori necessarily bound a solid torus in H1(12)and are not minimizing. Based on Fomenko’s conjecture, it is expected that these solid tori must be quite rigid in a well-defined sense: in homology, they should generate at most two transverse subgroups of H1(T*T2).

There is an alternative approach to Fomenko’s conjecture that is based on topological entropy. In a series of papers based on Glasmachers dissertation results, Glasmachers and Knieper study Riemannian Hamiltonians on T*T2with zero topological entropy [19, 20]. They prove the closure of one of the above-mentioned solid tori is a union of one or two closed, minimizing geodesic orbits and their stable and unstable manifolds ([20], Theorem 3.7c).[1] - The picture that emerges from their work is that there is a family of minimizing closed geodesics of the same homology class, and their stable and unstable manifolds, which bound a family of invariant solid tori. Bialy [5] describes the boundary of this set as a separatrix chain. The projection of the separatrix chain covers T2. A neighbourhood of the separatrix chain in the complement is fibred by invariant Lagrangian tori that are graphs, i.e. that are a union of minimizing orbits. The multipliers δi, or rather the angles θimentioned above, define sections of the unit cotangent bundle trapped within a separatrix chain.

Let us reformulate this as:

Problem 3.5. Prove the vanishing of the topological entropy of the geodesic flow of a Riemannian Hamiltonian on T*T2 that is completely integrable with a polynomial-in-momenta first integral F.

In various special cases, such as when Fis real-analytic or Morse-Bott, it is known that the topological entropy vanishes [21].

Finally, since topological entropy is an important invariant in the study of these systems, let us state a number of problems that are directly relevant to the preceding discussion. If one assumes Fomenko's conjecture is true and that the Liouville family of Riemannian Hamiltonians equals the set of completely integrable Riemannian Hamiltonians on T2, then it should be true that

Problem 3.6. The topological entropy of a non-Liouville Riemannian Hamiltonian on T*T2is positive.

Glasmachers and Knieper [20, 19] have studied the structure of geodesic flows with zero topological entropy on T*T2. The picture that emerges is the phase portrait looks remarkably like that of an integrable system. It seems likely that their results admit a strengthening: in particular, they are unable to determine the number of primitive homology classes represented by non-minimizing geodesics (for Liouville metrics, this is at most 4).

On the other hand, it is known, from results of Contreras, Contreras and Paternain and Knieper and Weiss that an open and dense set of Riemannian Hamiltonians have positive topological entropy [2224]. In the case of this particular problem, the natural point of departure is to look at Riemannian Hamiltonians that are close to Liouville, i.e. where the conformal factor in (35) is of the form

Λϵ=Λ0+ϵΛ1+O(ϵ2)E40

where Λ0is Liouville-and has no T1symmetry–and Λϵis not Liouville for all ϵ0. Based on the study in [25, 26] of the phase portrait of such systems, it should be possible to prove that the perturbed flow develops transverse homoclinic points.

### 3.3. The 2-sphere

The unit two-dimensional sphere S2R3admits a completely integrable geodesic flow. Indeed, the geodesic flow of an ellipsoid is also completely integrable with the second integral of motion that is, in general, a quadratic form in the momenta.

The fundamental problem is to describe the moduli space of completely integrable Hamiltonians on T*S2. The sub-problem of describing the integrable Riemannian (resp. natural or mechanical) Hamiltonians Hhas received wide-spread attention. When His Riemannian, the most common approach is to assume the second integral Fis polynomial-in-momenta, and without loss of generality, homogeneous. If the degree of Fis fixed, then the problem of determining H&Fis reducible to a non-linear PDE in the coefficients of F. When the degree is 1, the first integral Fis a momentum map of a T1isometry group (see below). When the degree is 2, then the Hamiltonian is Liouville, a classical result due to Darboux c.f. [27]. In degree 3, there is the well-known case due to Goryachev-Chaplygin, and more recent cases due to Selivanova, Dullin and Matveev and Dullin, Matveev and Topalov and Valent [2833]. In degree 4, Selivanova and Hadeler & Selivanova have produced a family of examples using the results of Kolokol’tsov [34, 27]. Beyond degree 4, Kiyohara has provided a construction of a smooth Riemannian metric Hwith an independent first integral Fof degree kfor any k1. In this construction, the metric Hdepends on a functional modulus, and so for each k, the set is infinite dimensional [35].

### 3.4. Super-integrable systems with a linear-in-momenta first integral

Let us review the work of Matveev and Shevchishin in more detail [36]. These authors impose an additional formal constraint that the metric possess one first integral that is linear-in-momenta. In conformal coordinates (x,y)where H=12c(x)(px2+py2), the existence of a cubic integral is reduced to a second-order ODE involving c.

From a geometric perspective, it is more natural to introduce coordinates adapted to the isometry group. That is, the existence of a linear-in-momenta first integral is equivalent to the existence of an isometry group containing T1. The action of T1on S1induces a cohomogeneity-1 structure. The fixed set of the T1action is a set of points {p,p+}which are equidistant along any minimal geodesic; and the principal T1-orbits are orthogonal to these geodesics. If γ:[T,T]S2is a minimal geodesic such that γ(±T)=p±, then we can let (r,θ)be ‘polar’ coordinates adapted to this structure. The Hamiltonian H and polynomial-in-momenta integral F can be written in the adapted coordinates as

H=12(pr2+s(r)pθ2),F=eivθ×j=0aj(r)prjpθNj,E41

where vZ, 3Nis a positive integer and the coefficients ajare to be determined. The equation {H,F}0is equivalent to a differential system that couples the coefficients a0,,aN, sand an anti-derivative Sof vs:

dS=vsdr,E42a
daj=12(N+2j)aj2dsaj1dS;(j=0,,N),E42b
ds=2vaN/aN1drE42c

where a2=a1=0. It is clear that the general solution of (42b), without the compatibility condition (42c), is obtained via repeated quadratures of products of s and S. The compatibility condition distinguishes those solutions which may arise from (41). The behaviour of s at r=±Tultimately determines whether the solution obtained arises from a T1-invariant Riemannian Hamiltonian Hand an independent first integral Fon T*S2.

In case N=3, the differential system reduces to a third-order nonlinear ODE similar to that studied by Chazy, in his generalization of the Painlêvé classification ([37], Eq. (6)). Based on the work of Matveev and Shevchishin [36], we know the solutions to this equation are real-analytic and define a parameterized family of super-integrable Riemannian metrics with cubic-in-momenta first integral. The latter authors do not solve the ODE explicitly.

Problem 3.7. Solve the N=3case of the differential system (42).

It appears to the author that this differential system may be soluble via hypergeometric functions. A successful resolution to the N=3case will naturally lead to the higher degree cases, which appear to be somewhat more involved.

Problem 3.8. Solve the higher degree cases of the differential system (42).

### 3.5. Super-integrable systems with a higher degree first integral

The author believes that the differential system 42 provides the key to understanding the subspace of super-integrable Riemannian Hamiltonians which admit a cohomogeneity-1 structure. Super-integrability alone does not imply the existence of such a cohomogeneity-1 structure. Without this additional hypothesis, there is very little known. Indeed, the extremely valuable construction of Kiyohara is the only construction that provides a smooth Riemannian Hamiltonian with a polynomial-in-momenta first integral of degree N>3super-integrable or not [35, 38].

Let us explain Kiyohara’s construction in some detail. Let H0be the Riemannian Hamiltonian of the standard unit sphere in R3. Let F0, F1be linear-in-momenta first integrals of H0that are linearly independent and let lk1be integers such that N=k+l3. Define a polynomial-in-momenta first integral G0=F0kF1l. For almost all qS2, the functions G0, H0|Tq*S2are dependent along two distinct lines through 0; this defines a pair of mutually transverse line bundles L±*over S2\{p0±,p1±}. The excluded, singular set consists of the anti-podal points pj±where Fjvanishes identically on the fibre (equivalently, the corresponding Killing field vanishes). This pair of line bundles provides a branched double covering

Φ:T2=R2/2πZ2S2E43

with simple branch points at {p0±,p1±}=Φ(πZ2). The line bundles L±*pullback to the line bundles Rdxjon T2={(x1,x2)mod2πZ}. Kiyohara shows that in these coordinates, the pullback of the function rwhich measures the time along the unique geodesic γthrough {p0±,p1±}(see Figure 4) satisfies the second-order PDE

2rx1x2+1B1+B2B1x2rx1+1B1+B2B2x1rx2=0E44

where B1and B2are functions that describe the line bundles L±in terms of the basis {dr,sin(r)dθ}.

Kiyohara writes a function R=r0+rwhere r0is the solution to (44) given by Φ*rand ris a solution of (44) with C2small boundary conditions satisfying

r(s,0)=u1(s),r(0,s)=u2(s),E45a
whereui(s)=ui(s)=ui(πs),for alli,s,andE45b
ui([ϵ,ϵ])=0.E45c

Then, by means of this perturbed function R, Kiyohara writes down an explicit formula for the perturbed Riemannian Hamiltonian Hand polynomial-in-momenta first integral F. The condition for the Poisson bracket {H,F}to vanish is shown to reduce to the satisfaction of Eq. (44) by R for the given values of B1and B2(this legerdemain is the real trick that makes the construction work).

Condition (45b) ensures that Rfactors through Φto a function on S2, while the condition (45c) ensures that Ris Con S2and coincides with ron a neighbourhood of the branch set {p0±,p1±}(hence that Hand Fcoincide with H0and F0, respectively, on a neighbourhood of the cotangent fibres of the branch set).

Let us now state several problems related to Kiyohara’s construction. First, Kiyohara’s vanishing condition on the boundary values (45c) is used to deduce the Riemannian Hamiltonians are not real-analytic. Since all the remaining constructions involve real-analytic data, this serves to show his examples are genuinely different.

Problem 3.9. Does Kiyohara’s construction extend to real-analytic boundary conditions u1, u2that satisfy (45b)? Do these real-analytic metrics include other known cases?

In particular, the obtained metrics are unlikely to have a T1isometry group, so the question is really whether the known examples in degree 3 and 4 are obtainable via this construction [1214, 2830, 34, 39, 40].

Second, Kiyohara’s construction produces a polynomial-in-momenta first integral Ffactors as A0lA1mwhere Aiare linear-in-momenta functions. It is clear that the reducibility of the first integral Fis forced by the desire to use a very simple branched covering.

Problem 3.10. Is reducibility of the first integral F necessary?

It ought to be fruitful to ask three related questions. The reducibility of Fis very special, with just two distinct factors.

Problem 3.11. Is it possible to extend Kiyohara’s construction so that the polynomial-in-momenta first integral F has more than 2 distinct linear factors?

It would be natural to try to extend the construction to the case where the zeros all lie on the same geodesic γ. More generally, one might attempt to mirror Kiyohara’s construction but in a more abstract way: start with a simple ramified covering Φ:S2with a branch set YS2. Let F0be a product of linear first integrals of H0that vanishes identically on TY*S2and not elsewhere. The stumbling block is that we need to clarify the intrinsic geometric meaning of the PDE that governs the perturbed systems (44).

Problem 3.12. Describe in explicit terms the third, independent first integral of H that is of least degree.

Kiyohara proves in his paper that His super-integrable (he proves the geodesic flow is 2π-periodic, in fact), but that proof does not proceed by finding this third first integral.

### 3.6. Three-dimensional configuration spaces

In comparison to the wealth of results and examples for surfaces that were surveyed above, comparatively little is known about the three-dimensional analogues. Tăĭmanov tells us that if the Tonelli Hamiltonian is completely integrable with real-analytic first integrals, then the three-dimensional configuration space has a finite covering p^:^such that the fundamental group π1(^)is abelian and of rank at most 3 [4143]. Based on the resolution of the Poincaré conjecture, this result implies that, up to finite covering the only such configuration spaces are

S3,S2×T1orT3.E46

The author generalized Kozlov’s result on surfaces to three-manifolds. In this result, if the Tonelli Hamiltonian is completely integrable and the singular set is topologically tame, then Tăĭamanov’s list extends to include those three-manifolds such that π1()is almost solvable (equivalently, due to the resolution of the geometrization conjecture, admits either a Nil or Sol geometry) [44]. Both results are sharp, like Kozlov’s, in the sense that all such admissible configuration spaces admit a geometric structure and the Riemannian Hamiltonian of such a structure is completely integrable with first integrals of the requisite type [45, 46].

There are a large number of questions that this strand of research has opened. Let us sketch a few.

### 3.7. The 3-sphere

The case of S3is perhaps best understood. It has been known since Jacobi proved the complete integrability of the geodesic flow of an ellipsoid via separation of variables, that the Liouville family of metrics on S3is completely integrable. These systems possess three independent quadratic-in-momenta first integrals.

Based on the analogous problem for the two-sphere,

Problem 3.13. Describe the structure of the super-integrable Riemannian Hamiltonians on S3.

Researchers who specialize in super-integrable classical and quantum systems have developed tools for constructing and classifying super-integrable systems c.f. [4749]. Unfortunately, some key ingredients in these constructions lead to systems with singularities.

The first method is based on the cohomogeneity-1 structure of S3with the group G=SO(3)acting as the linear isometry group of R3R4. If one represents

S3={(x,r)|xR3,rR,|x|2+|r|2=1},E47

then we see that Gacts freely on T*S3\TF*S3where F={(0,±1)}is the fixed-point set of the G-action on S3. This is enough to see that any G-invariant Hamiltonian on T*S3is non-commutatively integrable (analogous to the same fact for S2). If K:so(3)*Ris a positive-definite quadratic form, and Ψ:T*S3so(3)*is the momentum map of the SO(3)-action, then an invariant Riemannian Hamiltonian can be written as

H=12pr2+12s(r)Ψ*K,E48

for some function s>0such that s×(1±r)2const.0as r1.

If one employs the ansatz of Matveev & Shevchishin (c.f. Section 3.3), one would like to find first integrals that are polynomial-in-momenta of the form

F=j=0Nbj(x,r)prjΨ*ηNjE49

where ηNj:so(3)*Ris a homogeneous polynomial of degree Nj. In (41), the pre-factor exp(ivθ)appears to ensure that the coefficients of the first integral Fare common eigenfunctions of the Casimir ΔS1=2θ2parameterized by r. In the current case, the ansatz suggests that the coefficients bjshould factor as ϕλ(θ)aj(r)where ϕλis an eigenfunction of the Casimir ΔS2with eigenvalue λand θ=x/|x|.

Problem 3.14. Extend the construction sketched, above to higher dimensional spheres.

## Notes

• If one prefers a purely “elementary” proof, one might apply the inverse function theorem at this point.
• The existence of M is a consequence of Darboux’s theorem. Of course, a less elementary proof would appeal to Weinstein’s theorem and Moser’s isotopy lemma.
• Although the minimizing orbits have stable and unstable manifolds, it is not suggested that they are hyperbolic.

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Leo T. Butler (May 3rd 2017). Topology and Integrability in Lagrangian Mechanics, Lagrangian Mechanics, Hüseyin Canbolat, IntechOpen, DOI: 10.5772/66147. Available from:

### chapter statistics

1Crossref citations

### More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

### Related Content

#### Lagrangian Mechanics

Edited by Hüseyin Canbolat

Next chapter

By Yury Yanilkin

#### Electrostatics

Edited by Hüseyin Canbolat

First chapter

#### Electrostatic Interactions in Dense DNA Phases and Protein-DNA Complexes

By Andrey G. Cherstvy

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books