Next Article in Journal
Optimal Stopping and Loading Rules Considering Multiple Attempts and Task Success Criteria
Previous Article in Journal
Remarks on the “Onsager Singularity Theorem” for Leray–Hopf Weak Solutions: The Hölder Continuous Case
Previous Article in Special Issue
Darboux Transformation and Soliton Solution of the Nonlocal Generalized Sasa–Satsuma Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrable Systems: In the Footprints of the Greats

by
Velimir Jurdjevic
Department of Mathematics, University of Toronto, Toronto, ON M5S 3G3, Canada
Mathematics 2023, 11(4), 1063; https://doi.org/10.3390/math11041063
Submission received: 6 January 2023 / Revised: 4 February 2023 / Accepted: 16 February 2023 / Published: 20 February 2023
(This article belongs to the Special Issue Completely Integrable Equations: Algebraic Aspects and Applications)

Abstract

:
In his 1842 lectures on dynamics C.G. Jacobi summarized difficulties with differential equations by saying that the main problem in the integration of differential equations appears in the choice of right variables. Since there is no general rule for finding the right choice, it is better to introduce special variables first, and then investigate the problems that naturally lend themselves to these variables. This paper follows Jacobi’s prophetic observations by introducing certain “meta” variational problems on semi-simple reductive groups G having a compact subgroup K. We then use the Maximum Principle of optimal control to generate the Hamiltonians whose solutions project onto the extremal curves of these problems. We show that there is a particular sub-class of these Hamiltonians that admit a spectral representation on the Lie algebra of G. As a consequence, the spectral invariants associated with the spectral curve produce a large number of integrals of motion, all in involution with each other, that often meet the Liouville complete integrability criteria. We then show that the classical integrals of motion associated, with the Kowalewski top, the two-body problem of Kepler, and Jacobi’s geodesic problem on the ellipsoid can be all derived from the aforementioned Hamiltonian systems. We also introduce a rolling geodesic problem that admits a spectral representation on symmetric Riemannian spaces and we then show the relevance of the corresponding integrals on the nature of the curves whose elastic energy is minimal.
MSC:
53C17; 53C22; 53B21; 53C25; 30C80; 26D05; 49J15; 58E40

1. Introduction

The theory of integrable systems begins with W.R. Hamilton who in 1835 pronounced that the equations of motion of an n body system conform to the principle of least action, and consequently can be represented as
d q i d t = H p i , d p i d t = H q i , i = 1 , , n ,
under the transformation
p i = T q ˙ i ( q 1 , , q n , q ˙ 1 , , q ˙ n ) , q ˙ i = d q i d t , i = 1 , , n ,
where H = T + V  is the total energy, with T the kinetic and V the potential energy of the system. He then observed that H is conserved along the solutions of the system. Hamilton’s discovery gave rise to a new class of differential equations of the form
d x d t = H y ( x ( t ) , y ( t ) ) d y d t = H x ( x ( t ) , y ( t ) ) .
associated with any function H of 2 n  variables x = x 1 , , x n  and y = y 1 , , y n . Such equations became known as the canonical equations. Then the transformations ( x , y ) ( x , y )  that preserved the canonical form of these equations were also called canonical, and the functions whose values were conserved by canonical systems became known as integrals.
Hamilton’s discovery had an immediate impact on the scientific community of the nineteenth century. Canonical equations became the central object of study in the mathematics of that period with the contributions of J. Liouville, S.D. Poisson, C.G. Jacobi and H. Poincaré leading the way towards a new branch in mathematics known today as the theory of integrable systems. This theory was principally driven by a lasting interest in the existence of extra integrals of motion and the symmetries that are accountable for the existence of these integrals. Its defining moment may be attributed to S.D. Poisson who in 1809 [1] introduced his bracket (known since as the Poisson bracket)
{ f , g } = i = 1 n f x i g y i f y i g x i
for functions f and g in the canonical variables x 1 , , x n , y 1 , , y n .
The introduction of the Poisson bracket greatly facilitated the emerging theory of that period. It provided an alternative definition of canonical systems as differential systems that satisfy
d x i d t = { x i , H } , and d y i d t = { y i , H } , i = 1 , , n
and it also redefined integrals of motion associated with H as functions F that satisfy { F , H } = 0 . It was Jacobi, however, who noticed the fundamental property of the Poisson bracket
{ f , { g , h } } + { h , { f , g } } + { g , { h , f } } = 0 ,
that has been known ever since as the Jacobi’s identity. It is then an easy consequence of Jacobi’s identity that F 3 = { F 1 , F 2 }  is a third integral of motion for H for any two integrals F 1  and F 2  (known as Poisson’s theorem [2]). Alternatively integrals of motion were detected through a suitable change of canonical coordinates. Jacobi characterized such changes of coordinates through a generating function S ( x , y ) . According to Jacobi ( x , y ) ( x , y )  is canonical if and only if y i = S x i , x i = S y .  Poincaré characterized canonical change of coordinates in terms of differential forms: ( x , y ) ( x , y )  is canonical if and only if i = 1 n x i d y i x i d y i = d S  for some function S.
From contemporary perspectives the theory of integrable systems begins with C.G. Jacobi and his seminal book Lectures in Dynamics [3]. Jacobi demostrated that the canonical Equation (3) can be integrated with the aid of a partial differential equation
H ( x 1 , , x n , d S d x 1 , , d S d x n ) = c ,
in terms of an unknown function S. He showed that if a particular solution of (7) can be found in terms of n arbitrary constants of motion h 1 , , h n  then c = ϕ ( h 1 , , h n )  for some function ϕ , and the transformation
y i = S x i , h = S h i
transforms the canonical coordinates ( x , y )  into new canonical coordinates ( h , h )  relative to which the canonical Equation (3) are transformed into the equations
d h d t = d ϕ d h , d h d t = d ϕ d h = 0 ,
whose solutions are given by
h ( t ) = c 2 t + c 3 , h ( t ) = c 1 , c 2 = d ϕ d h .
Canonical coordinates whose solutions are given by (9) are called action-angle coordinates [4].
Equation H ( x 1 , , x n , S x 1 , , S x n ) = c  is known as Jacobi’s equation. Poincaré referred to the above result as the first theorem of Jacobi in his treatise of celestial mechanics [2]. Jacobi’s solution of the above partial differential equation in terms of the elliptic coordinates stands out as the most original and, perhaps, the most enigmatic contribution to the theory of canonical systems. Jacobi’s use of elliptic coordinates suggested the existence of a special class of variational problems whose solutions can be described by Abelian integrals in some privileged system of coordinates, exemplified by the geodesic problem on the ellipsoid. In the absence of any apparent symmetries on the ellipsoid that account for the integrability of the geodesic problem, this result of Jacobi seemed particularly mysterious.
In Jacoby summary, the main problem in the integration of differential equations appears in the choice of right variables. Given no general rule for finding the right choice, it is better to introduce special variables first, and then investigate the problems that naturally lend themselves to these variables [3]. Jacobi, however, does not comment on another exceptional aspect of his discovery, namely the mysterious presence of partial differential equations for the problems of variational calculus, an issue that remained open for a long time.
Almost a hundred years later, C. Carathéodory in the introduction to his famous book on the calculus of variations [5] remarks that “ neither Jacobi, nor his students, nor the many other prominent men who so brilliantly represented and advanced this discipline during the nineteenth century, thought in any way of the relationship between the calculus of variations and partial differential equation”. H. Poincaré also sidestepped this issue by treating canonical systems as the solutions of a dynamical system
d d t k = 1 2 n x i d y i d α k d d α k k = 1 2 n x i d y i d t = d F d α k , k = 1 , , 2 n ,
where α 1 , , α 2 n  denote the constants x i ( t 0 ) = α i , y i ( t 0 ) = α i + n , i = 1 , , n . Since d F d α k = i = 1 n F x i x i α k + F y i y i α k  the above differential equation can be reformulated as
i = 1 n ( d x i d t F y i ) d y i d α k ( d y i d t + F x i ) d x i d α k = 0 ,
which shows that Equations (3) and (10) have the same solutions. Poincaré equation used Equation (10) to show that a transformation ( x , y ) ( x , y )  is canonical if and only if the differential form i = 1 n y i d x i  satisfies i = 1 n y i d x i y i d x i = d S  for some function S ( x , x ) .
Among many other stellar advancements of that epoch, the following result of J. Liouville, reported in 1855 [6], seemed particularly influential for the present mathematics [4]. Liouville considered a differential system
d x d t = y F ( t , x ( t ) , y ( t ) ) , d y d t = x F ( t , x ( t ) , y ( t ) )
associated with a function F ( t , x 1 , , x n , y 1 , , y n ) . He then assumed the existence of n integrals of motion h 1 ( t , x , y ) , , h n ( t , x , y )  such that the system of equations
h 1 = h 1 ( t , x , y ) , h 2 = h 2 ( t , x , y ) , , h n = h n ( t , x , y )
can be solved for y 1 , , y n  in the variables t , x 1 , , x n , h 1 , h n . He also imposed the condition that h 1 , , h n  are in involution, that is,
{ h i , h j } = k = 1 n h i x k h j y h h i y k h j x k = 0 , 1 i , j n .
Liouville interpreted d x i d t = y i F ( t , x ( t ) , y ( t ) )  as the exactness condition for the differential form i = 1 n y i ( t , x , h ) d x i F ( t , x , p ( t , x , h ) ) d t  and concluded that there is a function S ( t , x , h )  such that
i = 1 n y i d x i F ( t , x , p ( t , x , h ) ) d t = S x i d x i + d S d t d t ,
that is, y i = S x , d S d t + F ( t , x , S x ) = 0 . But then S can be used as the generating function for the canonical transformation ( x , y ) ( h , h )  where h = S h . Liouville refers to
y i = S x i , h i = S h i , i = 1 , , n
as a complete system. Indeed, in the new coordinates h 1 , , h n  remain constants of motion and therefore d h i d t = 0 , i = 1 , , n . Since d h i d t = F h , F is independent of h , that is, F is a function of t and h. But then F ˜ h  is a given function of time, and h ( t )  is given by its integral. When F is a function of x and y, and not explicitly dependent on time, then F ˜  is only a function of h. Therefore, the general solution is given by
h ( t ) = h ( 0 ) , h ( t ) = ω t + h ( 0 ) , ω = F ˜ h .
This heritage from 19-th century mathematics forms a core of knowledge indispensable for problems of mathematical physics, symplectic geometry, calculus of variations and optimal control theory, and its unanswered questions still motivate much of the current research in integrable systems.
This paper will address the “hidden” symmetries that account for the existence of extra integrals of motion. We will show that the canonical integrable systems, such as Jacobi’s geodesic problem on the ellipsoid, Neumann’s mechanical problem on the sphere, Euler’s top, and the associated heavy tops, all derive their constants of motion from certain “meta” systems on Lie groups that admit isospectral representations of the form
d L λ d t ( t ) = [ M λ ( t ) , L λ ( t ) ]
on the Lie algebra g  of G.
We will confine our attention to semi-simple Lie groups G having a compact subgroup K, for then the Lie algebra g  admits a decomposition g = k + p , where k  the Lie algebra of K and p  is the orthogonal complement of k  relative to the Killing form K l ( A , B ) = T r ( a d A a d ( B ) . But then [ p , k ] p ]  and therefore g  as a vector space also carries the semi-direct Lie algebra g s  associated with the semi-direct product G s = p K . We will then single out a class of left-invariant variational problems on G that admit an isospectral representation with
L λ , s = L p λ L k + ( λ 2 s ) A ,
where s = 0  in the semi-direct case and s = 1  in the semi-simple case, L = L k + Ł p , L k k , L p p , and where A is a fixed element in p . It is then known that the spectral invariants ϕ λ , s k ( L ) = T r ( L λ , s k )  are in involution relative to the canonical Poisson bracket on g , respectively on g s . We will show that these invariants shed light on the hidden symmetries that surround many of the aforementioned integrable systems. In the process we will be able to demonstrate that the quest for the geometric origins behind the “mysterious” integrals of motions also leads to new and unexpected encounters with problems of Riemannian and sub-Riemannian geometry in which geometric control theory plays a major role.

2. Symplectic Background, Hamiltonian Systems

The theoretic framework upon which above claims are made is rooted in symplectic geometry. Below is a brief summary of the theoretical ingredients required for our main results.
Recall that a manifold M together with a non-degenerate and closed 2-form ω  is called symplectic. The symplectic form yields a correspondence between functions and vector fields: to every function f there is a vector field f  defined by ω ( f , X ) = d f ( X )  for all vector fields X on M. Then f  is called the Hamiltonian vector field generated by f. Every symplectic manifold is even dimensional, and at each point of M there is a neighbourhood with coordinates ( x 1 , , x n , p 1 , , p n )  on which Hamiltonian vector fields are given by
f = i = 1 n f p i x i f x i p i .
This choice of coordinates in which f  is given by (17) is called symplectic, or canonical in the terminology of the 19-th century.
Every cotangent bundle T * M  is a symplectic manifold with its canonical symplectic form, ω = d p d x  in terms of the symplectic coordinates ( x 1 , , x n , p 1 , , p n ) . As a symplectic manifold the cotangent bundle is special, in the sense that it is also a vector bundle. Hence every vector field X on M can be lifted to a unique Hamiltonian vector field f X  in T * M  via the function f X ( ξ ) = ξ ( X ( x ) ) , ξ T x * M . Vector field f X  is called the Hamiltonian lift of X. The same procedure is applicable to any time varying vector field, and by extension to any differential system on M. Thus any differential system in M can be lifted to a Hamiltonian system in T * M . This fact is also important for problems of optimal control where the Maximum Principle singles out the appropriate Hamiltonian lifts that govern the optimal solutions [7].
When the base manifold is a Lie group G, and when the underlying differential system is either left or right invariant, then there is a special system of coordinates based on the representation of T * G  as G × g * , with g *  the dual of g . This coordinate system preserves the left invariant symmetries and elucidates the conserved quantities of the associated Hamiltonian systems. The passage to these coordinates and the associated formalism was amply documented in my earlier publications [7,8,9]. Below we will highlight the main points in this theory required for our results.

2.1. Left-Invariant Trivializations and the Symplectic Form

Having in mind applications that involve left-invariant variational systems the cotangent bundle T * G  and the tangent bundle T G  will be viewed as the products G × g *  and G × g  via the left-translations. More explicitly, tangent vectors v T g G  will be identified with pairs ( g , X ) G × g  via the relation v = L g * X , where L g *  denotes the tangent map associated with the left translation L g ( h ) = g h . Similarly, points ξ T g * G  will be identified with pairs ( g , ) G × g *  via ξ = · L g * 1 . Then  T ( T * G ) , the tangent bundle of the cotangent bundle  T * G , will be identified with  ( G × g * ) × ( g × g * ) , with the understanding that an element  ( ( g , ) , ( A , a ) ) ( G × g * ) × ( g × g * )  denotes the tangent vector  ( A , a )  at the base point  ( g , ) .
We will make use of the fact that  G × g *  is a Lie group in its own right since  g * , as a vector space, is an abelian Lie group. Then left-invariant vector fields V in  G × g *  will be denoted by  V ( g , ) = ( g A , a ) ( g , )  in  G × g * . In this setting the canonical symplectic form on  T * G  is given by
ω ( g , ) ( V 1 , V 2 ) = a 2 ( A 1 ) a 1 ( A 2 ) ( [ A 1 , A 2 ] )
for any left-invariant vector fields  V 1 = ( g A 1 , a 1 )  and  V 2 = ( g A 2 , a 2 )  [7]. The above form is invariant under the left-translations in  G × g * , and is especially revealing for the Hamiltonian vector fields generated by left-invariant functions on  G × g * .
A function H on  G × g *  is left-invariant if  H ( h g , ) = H ( g , )  for all  g , h G  and all  g * . That is, left-invariant functions coincide with functions of  g * . Each left-invariant vector field  X ( g ) = g A  on G lifts to a linear function  ( A )  on  g *  because
h X ( ξ ) = ξ ( X ( g ) ) = L g * 1 ( L g ) * ( A ) = ( A ) , ξ T g * G .
Functions H on  g *  generate Hamiltonian vector fields  H  on  G × g *  whose integral curves are the solutions of
d g d t ( t ) = g ( t ) d H ( t ) , d d t ( t ) = ad * d H ( t ) ( ( t ) ) .
In a more general case, where H depends on both  g G  and  g * , the integral curves of  H  are the solutions of
d g d t ( t ) = g ( t ) d H ( t ) , d d t ( t ) = ad * d H ( t ) ( ( t ) ) d H g L g * ,
that can be easily shown through the relations
b ( d H ) + d H g L g * B = b ( A ) a ( B ) [ A , B ] .
This situation occurs in problems of mechanics in the presence of potential functions. For example, the movements of a three-dimensional rigid body with a potential function  V : S O ( 3 ) R  are described by the Hamiltonian
H ( R , ) = H 0 ( ) + V ( α 1 , α 2 , α 3 )
on the cotangent bundle of  S O ( 3 ) , where  α 1 , α 2 , α 3  denote the columns of the matrix transpose of the rotation R in  S O ( 3 ) . For then the directional derivative of V in the direction  R X  is given by
d V ( R X ) = i = 1 3 V α i α i , X
where  ,  denotes the standard inner product  1 2 T r ( X Y )  in  so ( 3 ) . Thus  d H g d L g = i = 1 3 V α i α i  and the equations of motion for H are given by
d g d t ( t ) = g ( t ) d H 0 ( ( t ) ) , d d t ( t ) = ad * d H 0 ( ( t ) ) ( ( t ) ) + i = 1 3 α i V α i .
These equations extend to an “n-dimensional rigid body” with the Hamiltonian  H ( R , ) = H 0 ( ) + V ( α 1 , , α n )  where
d R d t = R ( t ) Ω ( t ) , d M d t = [ Ω ( t ) , M ( t ) ] + i = 1 n α i V α i P ( Ω ( t ) ) = M ( t ) , α i ( t ) = R T ( t ) e i , i = 1 , , n .
In this context,  M ( t )  is the generalization of the angular momentum,  Ω ( t )  is the generalization of the angular velocity,  P  is the generalized inertia tensor, and  i = 1 n α i V α i  is the external torque.

2.2. Poisson Manifolds, Coadjoint Orbits

Equation (19) lend themselves to an insightful description in terms of the Poisson structure on  g *  inherited from the symplectic form  ω . Recall that a manifold M together with a bilinear, skew-symmetric form
{ , } : C ( M ) × C ( M ) C ( M )
that satisfies
{ f g , h } = f { g , h } + g { f , h } ,   ( L e i b n i z s   r u l e ) ,   a n d { f , { g , h } } + { h , { f , g } } + { g , { h , f } } = 0 ,   ( J a c o b i s   i d e n t i t y ) ,
for all functions  f , g , h  on M, is called a Poisson manifold.
Every symplectic manifold is also a Poisson manifold with the Poisson bracket given by  { f , g } ( p ) = ω p ( f ( p ) , g ( p ) ) , p M . However, the converse may not be true due to the fact that the Poisson bracket may be degenerate at some points of M. Nevertheless, each function f on M induces a Poisson vector field  f  through the formula  f ( g ) = { f , g }  as in the symplectic case. Poisson vector fields clarify the relation with symplectic manifolds through the following fundamental fact: every Poisson manifold is foliated by the orbits of its family of Poisson vector fields and each orbit is a symplectic submanifold of M with its symplectic form  ω p ( f , h ) = { f , h } ( p )  [7].
The dual  g *  of a Lie algebra  g  is a Poisson manifold with the Poisson bracket
{ f , h } ( ) = ( [ d h , d f ] )
for any functions f and h on  g * . In the literature on integrable systems the bracket  { f , h } ( ) = ( [ d f , d h ] )  is known as the Lie-Poisson bracket [10]. We have taken its negative to be compatible with the projections of left-invariant Hamiltonian vector fields on  g *  (and also to agree with the sign conventions in [7]).
It follows that each function H on  g *  defines a Poisson vector field  H  on  g *  via the formula  H ( f ) ( ) = { H , f } ( ) = ( [ d H , d f ] )  in which case the integral curves of  H  are the solutions of
d d t ( t ) = ad * d H ( t ) ( ( t ) ) .
Thus, as we already mentioned above, each function H on  g *  may be simultaneously viewed as a Hamiltonian on  T * G , and a function on the Poisson space  g * . Of course, Poisson equations coincide with the projections of the Hamiltonian equations on  g *  .
Solutions of Equation (24) are intimately linked with the coadjoint orbits of G through the following proposition. due to of A.A. Kirillov [11] (the proof is also given in [7]).
Proposition 1. 
Let  F  denote the family of Poisson vector fields on  g *  and let  M = O F ( 0 )  denote the orbit of  F  through a point  0 g * . Then M is equal to the connected component of the coadjoint orbit of G that contains  0 . Consequently each coadjoint orbit is a symplectic submanifold of  g * .
Recall that the coadjoint orbit of G through a point  g *  is given by  Ad g * ( ) = { Ad g 1 , g G } .
The fact that the Poisson equations can be naturally restricted to coadjoint orbits implies useful reductions in the theory of Hamiltonian systems.

2.3. Representation of Coadjoint Orbits on Lie Algebras

On semi-simple Lie groups Poisson Equation (24) can be expressed on  g  as
d L d t = [ d H , L ] ,
because the Killing form, or any scalar multiple of it  ,  is non-degenerate, and invariant, in the sense that,  X , [ Y , Z ] = [ X , Y ] , Z , X , Y , Z g , and can be used to identify  g  with  g *  via the formula
L , X = ( X ) , g * , X g .
Then coadjoint orbits are identified with the adjoint orbits and the Poisson vector fields  f X ( ) = ad * X ( )  are identified with vector fields  X ( L ) = [ X , L ] . Each vector field  [ X , L ]  is tangent to an orbit at L, and  ω L ( [ X , L ] , [ Y , L ] ) = L , [ Y , X ] X , Y  in  g  is the symplectic form on each orbit  O ( L 0 ) .
In a reductive semi-simple Lie group G there is also the semi-direct product  G 0 = p K  described earlier which generates its own coadjoint orbits on the dual of the Lie algebra  g 0  of  G 0 . Recall that the Lie algebra  g 0  of  G 0  consists of pairs  ( A , B ) , A p , B k  together with the Lie bracket
[ ( A 1 , B 1 ) , ( A 2 , B 2 ) ] = ( [ A 1 , B 2 ] [ A 2 , B 1 ] , [ B 1 , B 2 ] ) .
When the elements  ( A , B ) g 0  are identified with the sums the sums  A + B  in  g g  as a vector space, carries a double Lie algebra; the semi-direct product Lie algebras  g 0 , and the semi-simple Lie algebra  g 1 = g . We then have
[ A + B , C + D ] = [ A , B ] s + [ A , D + [ B , C ] + [ A , D ] , s = 0 , 1 ,
for any  A , C  in  p  and any  B , D  in  k , with  s = 0  in the semi-direct case, and  s = 1  in the semi-simple case.
Since both  g  and  g 0  Lie algebras over the same vector space, the Poisson equations on  g 0 *  can be also represented on  g 0  via the quadratic form  , , but the resulting expression takes a slightly different form. To see the difference, let  d H = d H p + d H k  and  L = L p + L k  denote the decompositions of  d H  and L onto the factors  p  and  k . On the semi-direct product Poisson equations reduce to
d L k d t = [ d H k , L k ] + [ d H p , L p ] , d L p d t = [ d H k , L p ] .
This equation can be combined with the equations for the semi-simple case in terms of the parameter s as
d L k d t = [ d H k , L k ] + [ d H p , L p ] , d L p d t = [ d H k , L p ] + s [ d H p , L k ] , s = 0 , 1 .
One can show that
P = A d h ( P 0 ) , Q = [ A d h ( P 0 ) , X ] + A d h ( Q 0 ) , ( X , h ) G 0
is the coadjoint orbit through  P 0 p , Q 0 k  under the action of  G 0 = p K  when  0 g s *  is identified with  L 0 = P 0 + Q 0  in  g 0 , and when  = Ad ( X , h ) * ( 0 )  is identified with  L = P + Q  [7].
The adjoint orbits of a non-compact semi-simple Lie groups G can be realized as the cotangent bundles of flag manifolds [12], and the same has been shown recently for the coadjoint orbits under the action of the semi-direct products [13,14]. We will make use of that fact later on in the paper.

3. Affine-Quadratic Problems

As stated earlier, we will restrict our attention to semi-simple Lie groups G and compact subgroups K with zero centre. We refer to  ( G , K )  as a reductive pair. Then  g  and  k  will denote their Lie algebras, and  p  will denote the orthogonal complement of  k  in  g  relative to the Killing  K l ( A , B ) = T r ( a d A a d ( B ) )  in  g . Recall that  K l  is non-degenerate and satisfies
K l ( A , [ B , C ] ) = K l ( [ A , B ] , C ) , A , B , C in g .
Hence  p  is well defined and satisfies  [ p , k ] p  ( in fact,  [ p , k ] = p  because  g  is semi-simple). We will also assume that  [ p , p ] = k . Note that the Killing form is negative-definite on  k  because K has zero centre [15], hence any negative scalar multiple  ,  of it is positive definite on  k . We shall assume that such a scalar product is fixed.
An affine quadratic problem is defined through a positive definite quadratic form Q on  k , and a regular element A in the Cartan space  p . An element A in  p  is called regular if  { X p : [ A , X ] = 0 }  is an abelian subalgebra in  p . The corresponding affine-quadratic problem consists of finding the solutions  g ( t )  in G of the affine control system
d g d t = g ( t ) ( A + U ( t ) ) ,
generated by a square-integrable control  U ( t )  in  k  that transfers a given state  g 0  in G to a given terminal state  g 1  in time T with a minimal energy  1 2 0 T Q ( U ( s ) ) d s . Any positive definite quadratic form Q is of the form  Q ( U ) = 1 2 P ( U ) , U  for some self-adjoint and positive linear operator  P  on  k . Then there exists an orthonormal basis  U 1 , , U k  in  k  such that  P  is diagonal relative to it. That is, if  U ( t ) = i = 1 k u i ( t ) U i  then  P ( U ( t ) ) = i = 1 k c i u i ( t ) U i  for some constants  c 1 , , c k . Then (29) can be rewritten as as
d g d t = X 0 ( g ) + i = 1 k u i ( t ) X i ( g ) ,
where  X 0 , , X k  are the left-invariant vector fields with  X 0 ( g ) = g A  and  X i ( g ) = g U i , i = 1 , , k , with  1 2 0 T i = 1 k c i u i 2 ( t ) d t .  the energy associated with each solution. The most natural case occurs when  P = I , that is, when  c i = 1 , i = 1 , , k . We will refer to this case as the canonical affine-quadratic problem.
When A is regular, then (29) is controllable, a consequence of our assumption  [ p , p ] = k , that is, any terminal state  g 1  can be reached in some finite time  T > 0  from any initial state  g 0 . But then there is an optimal solution  ( g ¯ ( t ) , U ¯ ( t ) )  on the interval  [ 0 , T ]  for which the energy of transfer  0 T Q ( U ¯ ( s ) ) d s  is minimal (see [7] for the proof). Therefore the above optimal control problem is well-posed.
To each affine-quadratic problem there is an analogous “shadow problem” defined on the semi-direct product  G o = p K  defined by the same data as in the original problem. It follows that every affine space  Γ = { A + U : U k }  that defines an affine left-invariant system on G also defines a corresponding left-invariant affine system on the semi-direct product  G 0 . Thus behind every affine quadratic optimal problem on G there is a corresponding affine-quadratic “shadow” problem on the semi-direct product  G s . The shadow problem is also well defined in the sense that optimal solutions exist on some interval  [ 0 , T ]  for each pair of boundary points  g ( 0 ) = g 0  and  g ( T ) = g 1 .
According to Pontriyagin’s Maximum Principle every optimal trajectory generated by a bounded and measurable control is the projection of an extremal curve, and each extremal curve is an integral curve of a suitable Hamiltonian system on the cotangent bundle of the ambient space. The Maximum Principle is also valid for optimal problems with  L 2  controls over affine systems with quadratic costs ([16]).
Let now  g ( t )  be an optimal trajectory generated by a control  u ( t ) . According to the Maximum Principle,  g ( t )  is the projection of an extremal curve  ξ ( t )  in  T * G  along which the cost extended Hamiltonian
λ 2 i = 1 k c i u i 2 ( t ) + H 0 ( ξ ) + i = 1 k u i ( t ) H i ( ξ ( t ) ) , λ = 0 , 1
is maximal at  u ( t )  relative to all competing controls  u ( t ) . In this notation, each  H i  is the Hamiltonian lift of  X i , i.e.,  H i ( ξ ( t ) ) = ξ ( t ) ( X i ( g ( t ) ) ) . In the abnormal case, which we will not treat here,  λ = 0 , and the Maximum principle results in the constraints  H i ( ξ ( t ) ) = 0 , i = 1 , , k .  In the normal case,  λ = 1 , the maximality condition implies that the optimal controls are of the form  u i ( t ) = 1 c i H i ( ξ ( t ) ) , i = 1 , , k . Consequently, optimal solutions are the projections of solution curves of a single Hamiltonian vector field  H  generated by the Hamiltonian
H ( ξ ) = 1 2 i = 1 k 1 c i H i 2 ( ξ ) + H 0 ( ξ ) = 1 2 i = 1 k 1 c i ( ( U i ) ) 2 + ( A ) .
Recall that each lift  H i ( ξ )  is a linear function on  g *  given by  H i ( ξ ) = ( U i )  with  H 0 ( ξ ) = ( A ) . Thus H is left-invariant, hence its Hamiltonian equations are given by
d g d t = X 0 ( g ) + i = 1 n 1 c i H i ( ( t ) ) X i ( g ( t ) ) , d d t = a d * d H ( ( t ) ) ( ( t ) ) .
The associated Poisson equations can be now written in  g  as
d L k d t = [ P 1 ( L k ) , L k ] + [ A , L p ] , d L p d t = [ P 1 ( L k ) , L p ] + s [ A , L k ] , s = 0 , 1 ,
after the identification of  g *  with  L g  via the scalar product  , , and the decomposition  L = L k + L p , L k k , L p p  (Equation (27)). In the canonical case ( P = I ) the preceding equations reduce to
d L k d t = [ A , L p ] , d L p d t = [ L k , L p ] + s [ A , L k ] , s = 0 , 1 .
Note that  s L k , L k + L p , L p  is an integral for (32). This integral is a universal integral of motion in the sense that it remains constant for any left-invariant Hamiltonian on  g s .

3.1. Isospectral Representations

We now single out a remarkable class of affine-quadratic Hamiltonians that plays a prominent role in the theory of integrable systems. It consists of Hamiltonians  H = 1 2 P 1 L k , L k + L p , A  that admit a spectral representation of the form
d L λ d t = [ M λ , L λ ] , with   M λ = P 1 ( L k ) λ A , and L λ , s ( L ) = L p λ L k + ( λ 2 s ) B ,
for some element  B p  that comutes with A, where  L p  and  L k  are the solutions of the Poisson Equation (32). Such a class is called isospectral and  L λ ( s )  is called the associated spectral curve. This terminology has origins in J. Zimmerman’s PhD thesis in 2002, in which he showed that the rolling sphere problem is isospectral [17]. We will return to Zimmerman’s problem and relate its results to the canonical affine-quadratic problem [18].
For Hamiltonian systems that admit an isospectral representation, the discrete spectral invariants of L are replaced by the functional invariants  ϕ λ , s ( k ) ( L ) = T r a c e ( L λ , s k ( L ) ) . Remarkably, the functional invariants  ϕ λ , s k  are in involution with each other, both with respect to the semi-simple and the semi-direct product Lie bracket, and in some instances generate a sufficient number of integrals of motion to ensure complete integrability ([7], 9.2). For instance, the family of functions
F 0 = { ϕ λ , 0 k , k 1 , λ R } { h X : [ X , B ] = 0 , X k }
is completely integrable on each coadjoint orbit in  p K  [19]. This means that H is completely integrable on each coadjoint orbit in  p k  whenever H is in involution with the Hamiltonian lifts  h X ( L ) = L , X , X k , [ X , B ] = 0 . This implies that the canonical affine Hamiltonian is completely integrable on coadjoint orbits since each left-invariant vector field with values in the isotropy group of A is a symmetry for the canonical system. It is reasonable to expect that the analogous family of functions is also completely integrable on coadjoint orbits of G, but, to the best of my knowledge, the proofs have not yet appeared in the literature.
The focus on the affine-quadratic problem and the associated Hamiltonians allows for the following characterization of isospectral Hamiltonians (proved in [7]).
Theorem 1. 
An affine Hamiltonian  H = 1 2 P 1 L k , L k + L p , A  is isospectral if and only  [ P 1 ( L k ) , B ] = [ L k , A ]  for some element  B p  that commutes with A. In the isospectral case,  L p = s B  is an invariant set for equations (32). On this set (32) are given by
d L k d t = [ P 1 ( L k ) , L k ] ,
and admit the reduced spectral representation
d d t ( L k λ B ) = [ P 1 ( L k ) λ A , L k λ B ] .
This theorem shows that the fundamental results A.T. Fomenko, A. S. Mischenko, and V.V. Trofimov on integrable left-invariant Riemannian metrics on compact Lie groups [20,21] based on Manakov’s seminal work on the n-dimensional Euler’s top [22] are subordinate to the isospectral properties of the affine Hamiltonian system, in the sense that the spectral invariants of  L k λ B  on  k  are always in involution with a larger family of functions generated by the spectral invariants of  L λ = L p + λ L k + ( λ 2 s ) B  on  g s  associated with an affine Hamiltonian H.

3.2. Affine Hamiltonians and Mechanical Tops

Let us now draw comparisons between the semi-direct Poisson equations
d L k d t = [ P 1 ( L k ( t ) ) , L k ( t ) ] + [ A , L p ( t ) ] , d L p d t = [ P 1 ( L k ( t ) ) , L p ( t ) ]
and the “top-like” equations:
d R d t = R ( t ) ( P 1 ( M ( t ) ) ) , d M d t = [ P 1 ( M ( t ) ) , M ( t ) ] + i = 1 n α i ( t ) V α i ,
associated with the Hamiltonian  H = 1 2 P 1 ( M ) , M + V ( α 1 , , α n ) . We will consider two cases- tops with linear potentials and tops with quadratic potentials.
Linear potentials. Equation (38) will be referred to heavy top-like equations when the potential energy V is generated by a linear Newtonian field, that is, when  V = i = 1 n c i ( α i , a ) , where a is a vector in  R n , and  c 1 , , c n  are constants. When  a = 0 , the external torque  i = 1 n α i ( t ) V α i  is equal to zero, and Equation (38) reduces to the Hamiltonian equation associated with a left-invariant Riemannian metric induced by the operator  P  (called the n-dimensional Euler’s top in some Russian literature [20]).
Heavy top-like equations can be written more compactly as
d R d t = R ( t ) Ω ( t ) , d M d t = [ Ω ( t ) , M ( t ) ] + a p ( t ) ,
where  Ω ( t ) = P 1 M ( t ) , and  p ( t ) = i = 1 n c i α i ( t ) . Since  α i ( t ) = R ( t ) T e i p ( t )  is a solution of  d p d t = Ω ( t ) p ( t ) . Hence each solution resides on the sphere  { p R n : | | p ( t ) | | = | | p ( 0 ) | | } .
Our theorems below relate Equation (39) to the Poisson Equation (37) on the reductive Lie algebras  so ( n + 1 )  and  so ( 1 , n )  associated with reductive pairs  ( S O ϵ , K )  where  S O ϵ  is  S O ( n + 1 )  when  ϵ = 1  and  S O ( 1 , n )  when  ϵ = 1 , and  K = { 1 } × S O ( n ) .
We will tackle both cases simultaneously but first we will need to introduce additional notation and terminology. We will use  so ϵ  to denote the Lie algebra of  S O ϵ  endowed with the trace form  A , B = 1 2 T r ( A B ) . Relative to  S O ϵ  we define its invariant bilinear form  ( x , y ) ϵ = x 0 y 0 + ϵ i = 1 n x i y i  in the ambient space  R n + 1 .
Then  a ϵ b , a R n + 1 , b R n + 1 , will denote the matrix defined by
( a ϵ b ) x = ( a , x ) ϵ b , x R n + 1 .
and  a ϵ b  denotes the matrix  a ϵ b b ϵ a . Since
( ( a ϵ b ) x , y ) ϵ + ( x , ( a ϵ b ) y ) ϵ = 0 ,
a ϵ b  belongs to  so ϵ ( n + 1 )  for any  a , b  in  R n + 1 . We then have
Theorem 2. 
Heavy top-like Equation (39) are isomorphic to the Poisson Equation (37) on the coadjoint orbit through  P 0 = p ( 0 ) ϵ e 0 , Q 0 = 0 0 0 M ( 0 )  under the coadjoint action of  p ϵ S O ( n ) . The passage to the affine Hamiltonian is via the following correspondences
A = ϵ a ϵ e 0 , L p = p ϵ e 0 , L k = 0 0 0 M , P 1 ( L k ) = 0 0 0 P 1 ( M ) .
For a proof see [14]. The preceding theorem clarifies the presence of heavy tops in the Hamiltonian equations on Lie algebras [10]. It also clarifies the relation between the tops and elastic rods initiated by G. Kirchhoff known as the “kinetic analogues” [23,24]. It also proves that the classification of completely integrable elastic rods in [7,8] carries over to the heavy tops.
Quadratic potentials. We will now show that the tops with quadratic potential V are also present in the equations of affine Hamiltonians, but this time on the tangent bundle of  S L ( n ) , or more precisely on the tangent bundle of the semi-direct product  s y m 0 ( n ) S O ( n )  where  s y m 0 ( n )  denotes the space of symmetric  n × n  matrices with zero trace. For that purpose let
H ( R , M ) = 1 2 ( P 1 ( M ) , M ) + 1 2 i = 1 n a i S α i , α i ,
with  R S O ( n ) , M s o ( n ) R T e i = α i , and S a symmetric  n × n . In accordance with (38) the Hamiltonian equations of  H  are given by
d R d t = R ( t ) Ω ( t ) , d M d t = [ Ω ( t ) , M ( t ) ] + i = 1 n a i α i ( t ) S α i ( t ) ,
where  Ω ( t ) = P 1 ( M ( t ) ) .
Theorem 3. 
Top-like Equation (41) are isomorphic with the Poisson equations generated by the affine Hamiltonian  H = 1 2 P 1 ( L k ) , L k + L p , S  on the coadjoint orbit through  P 0 = i = 1 n a i ( e i e i ) ( 1 n i = 1 n a i ) I  and  Q 0 = M ( 0 )  under the action of the semi-direct product  s y m 0 ( n ) S O ( n ) .
Proof. 
Every solution  ( M ( t ) , R ( t ) )  of (41) generates symmetric matrices  L p ( t )  and  X ( t )  given by
L p ( t ) = A d h ( t ) P 0 = i = 1 n a i ( α i ( t ) α i ( t ) ) 1 n i = 1 n a i I , X ( t ) = Ad h ( t ) Y ( t ) , Y ( t ) = 0 t Ad h 1 ( s ) S d s ,
with  h ( t ) = R T ( t ) . Then,
d L p d t = i = 1 n a i ( Ω α i α i + α i Ω α i ) = [ Ω , L p ] , d X d t = [ Ω ( t ) , X ( t ) ] + Ad h ( t ) Y ˙ = [ Ω ( t ) , X ( t ) ] S .
Additionally,
[ S , L p ( t ) ] = i = 1 n ( a i ( α i α i ) S S a i ( α i α i ) ) = i = 1 n a i α i S α i a i S α i α i = i = 1 n a i ( α i S α i ) ,
which in turn implies that (41) can be written as
d R d t = R ( t ) Ω ( t ) , d M d t = [ Ω ( t ) , M ( t ) ] + [ S , L p ( t ) ] .
Let now  Q ( t ) = [ Ad h ( t ) ( P 0 ) , X ( t ) ] + Ad h ( t ) Q 0 = [ L p ( t ) , X ( t ) ] + Ad h ( t ) Q 0 . Note first that
[ [ Ω , L p ] , X ] = [ [ X , Ω ] , L p ] [ [ L p , X ] , Ω ] = [ [ X , Ω ] , L p ] + [ Ω , Q ] [ Ω , Ad h Q 0 ] .
Then,
d Q d t = [ [ Ω ( t ) , L p ( t ) ] , X ( t ) ] + [ L p ( t ) , d X d t ( t ) ] + [ Ω ( t ) , Ad h ( t ) ( Q 0 ] = [ Ω ( t ) , Q ( t ) ] + [ L p , [ X , Ω ( t ) ] ] + [ L p , d X d t ] = [ Ω ( t ) , Q ( t ) ] + [ S , L p ] .
Therefore  Q ( t )  and  M ( t )  satisfy the same differential equation. Hence  Q ( t ) = M ( t )  whenever  Q 0 = M ( 0 ) . If we now rename  Q ( t )  as  L k ( t )  we get the Poisson equations for the shadow Hamiltonian  H = 1 2 P 1 ( L k ) , L k + S , L p .   □
The preceding theorem links isospectral Hamiltonians to the equations of the top under quadratic potentials and paves a way to the n-dimensional generalization of O. Bogoyavlensky’s famous result on integrability of three-dimensional mechanical tops in the presence of a quadratic potential [25]. The path to isospectral Hamiltonians is provided by Manakov’s observation that the inertia tensor  P ( U ) , U  for a rigid body is confined to the transformations  P ( U ) = S U + U S , for some positive definite matrix S. For then,  [ P 1 ( M ) , S 2 ] = [ M , S ] . Indeed, in this situation  P ( U ) = S U + U S = M , and
[ P 1 M , S 2 ] = [ U , S 2 ] = [ S U + U S , S ] = [ M , S ] .
Hence the corresponding affine Hamiltonian  H ^ = 1 2 P 1 L k , L k + S , L p  is isospectral on  sl ( n )  (Theorem 1). Since the equations of the Hamiltonian
H = 1 2 P 1 M , M + i = 1 n a i ( α i , S α i )
corresponding to the top with quadratic potential  V = 1 2 i = 1 n a i ( S α i , α i )  can be identified with the Poisson equations of  H ^  on the coadjoint orbit through  L p = P 0 , L k = M ( 0 ) , the isospectral invariants of
L λ = i = 1 n a i α i λ M + λ 2 S
are integrals of motion for the top. (Theorem 3). Since  { X so ( n ) : [ X , S ] = 0 } = 0  for each non-singular symmetric matrix S, the spectral invariants of
L λ ( s ) = i = 1 n a i ( α i α i λ M + ( λ 2 s ) S
form a completely integrable family of functions on each coadjoint orbit in  sl ( n )  (semi-simple and semi-direct). the top with a quadratic potential is completely integrable in all dimensions.

3.3. Three-Dimensional Tops- Kirchhoff-Kowalewski Type

We will now turn our attention to the class of affine-quadratic systems of Kirchhoff-Kowalewski type on complex Lie algebras with a particular interest on the symmetries that account for the existence of Kowalewski’s integral reported in her seminal paper on the motions of a rigid body around a fixed point under the influence of gravity [26]. We will follow our recent paper [27] and show that there is a natural Hamiltonian on  sp ( 4 , C )  that answers the fundamental questions raised by Kowalewki’s paper, namely, what is the geometric rational behind her approach in which all the variables were treated as complex quantities, and secondly. what are the symmetries that account for the existence of not only her integral of motion, but also of similar integrals, known as Kowalewski type integrals, that subsequently appeared in the literature on integrable systems [8,28,29,30,31].
Theorem 2 suggests that the search for the answers to the above questions should begin with the Poisson equations associated with an affine-quadratic Hamiltonian on  so ( 4 , C )  since both  so ( 1 , 3 )  and  so ( 4 )  are real forms for  so ( 4 , C )  (see also [24]). We will show that Kowalewski’s “mysterious” change of variables appear naturally in the passage from  so ( 4 , C )  to  sl ( 2 , C ) × sl ( 2 , C )  an important intermediate step towards the right Hamiltonian on  sp ( 4 , C ) . The journey from  so ( 4 , C )  to  sp ( 4 , C )  to this remarkable Hamiltonian begins with
H = 1 2 ( m 1 λ 1 + m 2 λ 2 + m 3 λ 3 ) + b 1 p 1 + b 2 p 2 + p 3 b 3 ,
where  L = m 1 A 1 + m 2 A 2 + m 3 A 3 + p 1 B 1 + p 2 B 2 + p 3 B 3  is the coordinate representation of a point L in  so ( 4 , C )  relative to an orthonormal basis  A 1 , A 2 , A 3 , B 1 , B 2 , B 3  that conforms to the following Lie bracket Table 1:
Then
d L k d t = [ d H k , L k ] + [ B , L p ] , d L p d t = [ d H k , L p ] + s [ B , L k ] , s = 0 , 1 ,
are the Poisson equations generated by H, where  B = b 1 B 1 + b 2 B 2 + b 3 B 3  denote the drift element in  p d H k = i = 1 3 m i λ i A i L k = i = 1 3 m i A i , and  L p = i = 1 3 p i B i . The same equations can be also expressed as
d m 1 d t = m 2 m 3 ( 1 λ 3 1 λ 2 ) + p 2 b 3 p 3 b 2 , d m 2 d t = m 1 m 3 ( 1 λ 1 1 λ 3 ) + p 3 b 1 p 1 b 3 , d m 3 d t = m 1 m 2 ( 1 λ 2 1 λ 1 ) + p 1 b 2 p 2 b 1 , d p 1 d t = 1 λ 3 p 2 m 3 1 λ 2 p 3 m 2 + s ( m 2 b 3 m 3 b 2 ) , d p 2 d t = 1 λ 1 p 3 m 1 1 λ 3 p 1 m 3 + s ( m 3 b 1 m 1 b 3 ) , d p 3 d t = 1 λ 2 p 1 m 2 1 λ 1 p 2 m 1 + s ( m 1 b 2 m 2 b 1 ) .
When  s = 0  the above equations formally coincide with the equations of the top:  d M d t = [ Ω ( t ) , M ( t ) ] + b p , d p d t = Ω ( t ) p ( t )  (Equation (21)).
On  g s  there are two Casimirs:
I 1 = L p , L p + s L k , L k = ( p , p ) + s ( m , m ) , I 2 = L k , L p = ( m , p ) ,
Hence generic coadjoint orbits in  g s  are four-dimensional. Since each coadjoint orbit is symplectic, integrable cases occur whenever there is an extra integral of motion functionally independent of H I 1 , and  I 2 . Since the motion of the top is subordinate to the Poisson system of H on  se ( 3 , C ) , the search for integrable tops reduces to the search for an additional integral of motion functionally independent from  I 1 , I 2  and H.
Let us now come to the conditions of Kowalewski
λ = λ 1 = λ 2 = 2 λ 3 , b 3 = 0 .
and her “mysterious” variables
z 1 = m 1 + i m 2 , z 2 = m 1 i m 2 , w 1 = p 1 + i p 2 , w 2 = p 1 i p 2 , z 3 = i m 3 , w 3 = i p 3 , b = b 1 + i b 2 , b ¯ = b 1 i b 2 .
After the substitutions, Equation (46) become
d z 1 d t = 1 λ z 1 z 3 + b w 3 , d z 2 d t = 1 λ z 2 z 3 b ¯ w 3 d z 3 d t = 1 2 ( b w 2 b ¯ w 1 ) d w 1 d t = 1 λ z 1 w 3 2 λ z 3 w 1 + s b z 3 , d w 2 d t = 1 λ z 2 w 3 + 2 λ z 3 w 2 s b ¯ z 3 d w 3 d t = 1 2 λ ( z 1 w 2 z 2 w 1 ) + s 2 ( b z 2 b ¯ z 1 ) ,
from which it can be easily extracted that
I = ( z 1 2 2 λ b w 1 + 1 2 s λ b 2 ) ( z 2 2 2 λ b ¯ w 2 + 1 2 s λ b ¯ 2 ) ,
is an integral of motion. Following the terminology in [7] we will refer refer to this integral as the Kirchhoff-Kowalewski integral. It is only in the special case  s = 0  and  λ = 2  that this integral coincides with the integral of motion found by Kowalewski. The real versions of the Kirchhoff-Kowalewski integral were originally discovered by V. Kuznetsov and I.V. Komarov in their studies of the hydrogen atom [28,29].
Let us reveal the geometric rational behind Kowalewski’s change of variables. The explanations are most naturally articulated through the root system in  so ( 4 , C ) . Recall that any maximal commutative sub-algebra  h  of a Lie algebra  g  is called a Cartan subalgebra. All Cartan subalgebras in a semi-simple Lie algebra are conjugate, and hence all have the same dimension. The dimension of any Cartan algebra is the rank of  g . The rank of  so 4 ( C )  is two. Evidently each pair  ( A i , B i ) , i = 1 , 2 , 3 , in Table 1 generates a Cartan algebra in  so ( 4 , C ) . Since these algebras are conjugate, there is no preferential choice. However, in regard to the equations of the top, there is a preferential choice when two moments of inertia are equal. In the case that  λ 1 = λ 2  the natural choice is the Cartan algebra generated by the pair  { A 3 , B 3 } .
An element  α  in the dual  h *  of a Cartan algebra  h  is called a root if for some  v g [ h , v ] = α ( h ) v  for all  h h . An easy calculation shows that there are four roots  ± α 1 ± α 2  given by
α 1 ( x A 3 + y B 3 ) = i ( x + y ) , α 2 ( x A 3 + y B 3 ) = i ( x y ) , x , y C .
The corresponding root spaces are one dimensional, and are generated by
C 1 = 1 2 ( A 1 i A 2 ) + 1 2 ( B 1 i B 2 ) , α = α 1 , C 2 = 1 2 ( A 1 + i A 2 ) + 1 2 ( B 1 + i B 2 ) , α = α 1 , D 1 = 1 2 ( A 1 i A 2 ) 1 2 ( B 1 i B 2 ) , α = α 2 , D 2 = 1 2 ( A 1 + i A 2 ) 1 2 ( B 1 + i B 2 ) , α = α 2 .
Together with  C 3 = i 2 ( A 3 + B 3 )  and  D 3 = i 2 ( A 3 B 3 )  these matrices form a basis for  so ( 4 , C ) . A simple calculation shows that
α 1 ( C 3 ) = 1 , α 2 ( D 3 ) = 1 , α 1 ( D 3 ) = α 2 ( C 3 ) = 0 , hence , [ C 3 , C 1 ] = C 1 , [ C 3 , C 2 ] = C 2 , [ D 3 , D 1 ] = D 1 , [ D 3 , D 2 ] = D 2 .
Furthermore,  [ C 1 , C 2 ] = 2 C 3 [ D 1 , D 2 ] = 2 D 3 , and  [ C i , D j ] = 0 , for all i and j.
The Lie algebras  g 1  and  g 2  spanned by  C 1 , C 2 , C 3 , and  D 1 , D 2 , D 3  satisfy  [ g 1 , g 2 ] = 0 . and each is isomorphic to  sl ( 2 , C )  under the identification
C 1 , D 1 0 1 0 0 , C 2 , D 2 0 0 1 0 , C 3 , D 3 1 2 1 0 0 1 .
An easy calculation shows that the coordinates  a 1 , a 2 , a 3 , b 1 , b 2 , b 3  of an arbitrary point  X so ( 4 , C )  relative to the basis  A 1 , A 2 , A 3 , B 1 , B 2 , B 3  are transformed to the coordinates  c 1 , c 2 , c 3 , d 1 , d 2 , d 3  relative to the basis  C 1 , C 2 , C 3 , D 1 , D 2 , D 3  according to the following formulas:
c 1 = 1 2 ( a 1 + i a 2 ) + 1 2 ( b 1 + i b 2 ) , d 1 = 1 2 ( a 1 + i a 2 ) 1 2 ( b 1 + i b 2 ) , c 2 = 1 2 ( a 1 i a 2 ) + 1 2 ( b 1 i b 2 ) , d 2 = 1 2 ( a 1 i a 2 ) 1 2 ( b 1 i b 2 ) . c 3 = i ( a 3 + b 3 ) , d 3 = i ( a 3 b 3 ) .
Let us now  Φ : so ( 4 , C ) sp ( 4 , C )  be given by
Φ ( i = 1 3 ( c i C i + d i D i ) ) = 1 0 0 0 c 3 2 c 1 c 2 c 3 2 + 0 0 0 1 d 3 2 d 1 d 2 d 3 2 = I 1 2 1 2 ( c 3 + d 3 ) c 1 + d 1 ( c 2 + d 2 ) 1 2 ( c 3 + d 3 ) + 1 0 0 1 1 2 1 2 ( c 3 d 3 ) c 1 d 1 ( c 2 d 2 ) 1 2 ( c 3 d 3 ) = I 1 2 i a 3 a 1 + i a 2 a 1 + i a 2 i a 3 + 1 0 0 1 1 2 i b 3 b 1 + i b 2 b 1 + i b 2 i b 3 .
where  A B  denotes the Kronecker product of matrices A and B.
To see that  Φ ( so ( 4 , C ) ) sp ( 4 , C )  recall first that  sp ( 4 , C )  consists of matrices M that satisfy  J M J 1 = M T , where J is the matrix that defines the symplectic form  ( z , J w )  on  C 4 , i.e.,  J 2 = I . It is easy to check that both  I 1 2 i a 3 a 1 + i a 2 a 1 + i a 2 i a 3  and  1 0 0 1 1 2 i b 3 b 1 + i b 2 b 1 + i b 2 i b 3  satisfy  J M J 1 = M T  with  J = 1 0 0 1 0 1 1 0 . Since  J 2 = I 2 I 2 = I , our claim follows.
We will identify  sl ( 2 , C )  with pure complex quaternions  Q  via the correspondence
q = q 1 i + q 2 j + q 3 k Q = q 1 E 1 + q 2 E 2 + q 3 E 3 = i q 3 q 1 + i q 2 q 1 + i q 2 i q 3 .
Then the standard basis  i , j , k  in  Q  is identified with the matrices
E 1 = 0 1 1 0 , E 2 = 0 i i 0 , E 3 = i 0 0 i ,
and any element  X = a b c a  in  sl ( 2 , C )  is represented by the quaternion  X = 1 2 i a 3 a 1 + i a 2 a 1 + i a 2 i a 3 a 1 = b c a 2 = ( b + c ) a 3 = 2 i a .
Let now  A i = 1 2 E i A i = I A i , B i = E 3 A i , so that  Φ ( A i ) = A i  and  Φ ( B i ) = B i  for each i i = 1 , 2 , 3 , and let  so 4 = Φ ( so ( 4 , C ) ) . Matrices  A i , B i , i = 1 , 2 , 3 .  form an orthonormal basis in  so 4  relative to the inner product  X , Y = T r ( X Y )  on  sp ( 4 , C ) .
It is easy to verify that  Φ : so ( 4 , C ) so 4  is a Poisson map. Therefore
H ˜ ( ˜ ) = H ( Φ * ( ˜ ) ) = H ( ) , Φ * ( ˜ ) = , ˜ so 4 ,
for any function H on  so * ( 4 , C ) , where  Φ *  denotes the dual map of  Φ . After the identification of  so 4 *  with  so 4  via the trace form, the Poisson equations of  H ˜  associated with H in (44) become
d d t ( I Z ) = [ I Ω , I Z ] + [ E 3 B , E 3 W ] = I ( [ Ω , Z ] + [ B , W ] ) , d d t ( E 3 W ) = [ I Z , E 3 W ] + s [ E 3 B , I Z ] = E 3 ( [ Ω , W ] + s [ B , Z ] ) ,
or, in simpler form,
d Z d t = [ Ω , Z ] + [ B , W ] , d W d t = [ Ω , W ] + s [ B , Z ] ,
where  Z = 1 2 i m 3 m 1 + i m 2 m 1 + i m 2 i m 3 , W = 1 2 i p 3 p 1 + i p 2 p 1 + i p 2 i p 3 ,
Ω = 1 2 1 λ 3 z 3 1 λ 1 m 1 + i λ 2 m 2 1 λ 1 m 1 + i λ 2 m 2 1 λ 3 z 3 , and  B = 1 2 i b 3 b 1 + i b 2 b 1 + i b 2 i b 3 .
Now we see Kowalewski variables  z 1 = m 1 + i m 2 , z 2 = m 1 i m 2 , w 1 = p 1 + i p 2 , w 2 = p 1 i p 2  as the natural coordinates in this Poisson representation. Under Kowalewski’s conditions  λ = λ 1 = λ 2 = 2 λ 3 b 3 = 0 , Equation (55) reduce to Equation (49). The passage from H to  H ˜  reveals the geometric rational behind the ad-hoc change of variables in (48) and serves as a natural segue to our ultimate Hamiltonian on  sp ( 4 , C ) .

3.4. Kowalewski’s Conditions and Isospectral Representations

We now address the origins of the “enigmatic” conditions (47) through an extended affine-quadratic Hamiltonian
H = i = 1 3 m i 2 2 λ i + b i p i + c i q i
on  sp * ( 4 , C )  defined by complex numbers  λ 1 , λ 2 , λ 3 b 1 , b 2 , b 3 , and  c 1 , c 2 , c 3  and an extended basis  A i , B i , i = 1 , 2 , 3 A 4 = 1 2 E 2 I ,  and  C 1 = E 1 A 1 , C 2 = E 1 A 2 , C 3 = E 1 A 3 , where  E 1 = 0 1 1 0 , E 2 = 0 1 1 0 ,  and  E 3 = 1 0 0 1 .
The reader can easily verify that  g = sp ( 4 , C )  has the following decomposition
g = k p , p = p 1 p 2 , k = k 0 C A 4
where  k 0  is the Lie algebra spanned by  A 1 , A 2 , A 3  and where  p 1  and  p 2  are respectively the linear spans of  B i , i = 1 , 2 , 3  and  C i , i = 1 , 2 , 3 . These spaces conform to the following Lie algebraic relations:
[ A 4 , k 0 ] = 0 , [ A 4 , p 1 ] = p 2 , [ A 4 , p 2 ] = p 1 , [ k o , p 1 ] = p 1 , [ k 0 , p 2 ] = p 2 , [ p 1 , p 1 ] = k 0 , [ p 2 , p 2 ] = k 0 , [ p 1 , p 2 ] = C A 4 .
After  g *  is identified with  g  via the scalar product  X , Y = 1 2 T r ( X Y )  the above Hamiltonian can be written as
H = 1 2 P ( L k 0 ) , L k 0 + L p , A ,
where  L k 0 = m 1 A 1 + m 2 A 2 + m 3 A 3 , L p = i = 1 3 p i B i q i C i  and  A = B + C B = i = 1 3 b i B i C = i = 1 3 c i C i  (note that  ,  is negative on  p 2  which accounts for the negative signs in the expression for  L p ).
Since  g = sp ( 4 , C )  is semi-simple the Poisson equations for  H  are given by
d L k d t = [ d H k , L k ] + [ A , L p ] , d L p d t = [ d H k , L p ] + s [ A , L k ] , s = 0 , 1 .
These equations can be written in a more succinct form as
d Z d t = [ Ω , Z ] + [ B , W ] + [ C , S ] , d m 4 d t = ( C · W + B · S ) , d W d t = [ Ω , W ] + s ( [ B , Z ] m 4 C ) , d S d t = [ Ω , S ] s ( [ C , Z ] m 4 B ) .
in terms of the following notations:
L k = I Z + m 4 A 4 , L p = E 3 W E 1 S , d H = I Ω + E 3 B + E 1 C , Ω = i = 1 3 m i ˘ i A i , B = i = 1 3 b i A i , C = i = 1 3 c i A i .
Z = 1 2 z 3 z 1 z 2 z 3 , W = 1 2 w 3 w 1 w 2 w 3 , z 1 , 2 = m 1 ± i m 2 , z 3 = i m 3 , w 1 , 2 = p 1 ± i p 2 , w 3 = i p 3 , as in the previous section, and  S = 1 2 s 3 s 1 s 2 s 3 ,  with  s 1 , 2 = q 1 ± i q 2 , s 3 = i q 3 .
We now come to the crux of the matter, the existence of integrals of motion for the above system. The intermediate question is the existence of an integral I of the form  I = α 1 m 1 + α 2 m 2 + α 3 m 3 + β m 4  for some constants  α 1 , α 2 , α 3 ,  and  β .
Proposition 2. 
I = α 1 m 1 + α 2 m 2 + α 3 m 3 + β m 4  is an integral of motion for  H  in exactly two cases: when  λ 1 = λ 2  and  b 1 = b 2 = c 1 = c 2 = 0 , then  I = m 3 , and in the second case, when  λ 1 = λ 2 , b 3 = c 3 = 0 , b 1 = ± i c 2 , b 2 = i c 1 , then  I = i m 3 + m 4 . (for the proof see [27]).
This first condition singles out the top of Lagrange, while the second condition is a precursor to Kowalewski’s top as will be demonstrated below. Note that the second condition  b 1 = i c 2  and  b 2 = i c 1  can be also written as  b = c , and  b ¯ = c ¯  where  b = b 1 + i b 2 b ¯ = b 1 i b 2 c = c 1 + i c 2  and  c ¯ = c 1 i c 2 . Then  C = 1 2 0 b b ¯ 0 , and since it is orthogonal to  B = 1 2 0 b b ¯ 0 , it will be denoted by  B .
We will say that (59) satisfies the preliminary condition of Kowalewski whenever  λ 1 = λ 2  and  C = B . It follows that the preliminary condition of Kowalewski is synonymous with the integral of motion  I = i m 3 + m 4 = z 3 + m 4 . (This integral of motion was also discovered earlier by A.M Savu in [32]).
We will now assume that the preliminary condition holds and we will pursue conditions on the ratio  δ = λ λ 3 , where  λ = λ 1 = λ 2 , that guarantee extra integrals of motion for system (59). Note that in this situation  A = B + C = E 3 B + E 1 B . Systems that satisfy the preliminary condition of Kowalewski and also satisfy  δ = 2  will be said to satisfy the Kowalewsky conditions. The following proposition provides an important characterization of Kowalewski’s conditions for both  s = 0  and  s = 1 .
Proposition 3. 
Assume that (59) satisfies the preliminary condition of Kowalewski and is restricted to the manifold  z 3 + m 4 = 0 . Then  H  satisfies the isospectrality condition  [ d H k , B ] = [ L k , A ] (as in Theorem 1) for some matrix  B p , with  [ A , B ] = 0  if only if (59) satisfies the Kowalewski conditions. In fact,  B = λ ( E 1 B + E 3 B ) = λ A  [27].
Indeed, under Kowalewski’s conditions
Ω = 1 λ Z + 1 2 λ z 3 E 3 and L k = I Z z 3 A 4 .
Therefore,
[ d H k , B ] = [ I × Ω , B ] = [ I × Z , A ] + z 3 2 [ I × E 3 , A ] , and [ L k , A ] = [ I Z , A ] z 3 [ A 4 , A ] .
Since  [ z 3 2 I E 3 , A ] = z 3 [ A 4 , A ] [ I Ω , B ] = [ L k , A ] .
It follows from Theorem 1 that Kowalewski’s condition is necessary and sufficient for the existence of isospectral representation
d L μ d t = [ M μ , L μ ] , L μ = L p μ ( L k 0 z 3 A 4 ) + ( μ 2 s ) λ A
on the invariant manifold  z 3 + m 4 = 0 . Consequently,  ϕ k = T r ( L μ 2 k ) = L μ k , L μ k  are integrals of motion for (59), in involution with each other for each  s = 0 , or  s = 1 . Remarkably, the prototype of Kowalewski’s integrals of motion is found among the above spectral invariants. (see also [33,34] for other spectral representations).
We will show the existence of Kowalewki’s integral of motion directly from the equations
d Z d t = [ Ω , Z ] + [ B , U ] + [ B , V ] , d U d t = [ Ω , U ] , d V d t = [ Ω , V ] ,
obtained from (59) under the change of variables  U = W s λ B , V = S + s λ B .  Equation (60) may be seen as a semisimple extension of the Kowalewski-type gyrostat in two constant fields introduced in [35].
Equations (60) may be also expressed in terms of the coordinates as
d z 1 d t = 1 λ z 1 z 3 + b ( u 3 + v 3 ) , d z 2 d t = 1 λ z 2 z 3 b ¯ ( u 3 v 3 ) , d z 3 d t = 1 2 ( b u 2 b ¯ u 1 + b v 2 + b ¯ v 1 ) , d u 1 d t = u 3 z 1 λ 2 u 1 z 3 λ , d u 2 d t = 2 u 2 z 3 λ z 2 u 3 λ , d u 3 d t = 1 2 λ ( z 1 u 2 u 1 z 2 ) d v 1 d t = v 3 z 1 λ 2 v 1 z 3 λ , d v 2 d t = 2 v 2 z 3 λ z 2 v 3 λ , d v 3 d t = 1 2 λ ( z 1 v 2 v 1 z 2 ) .
One readily obtains the following fundamental equalities
d d t ( u 1 + v 1 ) = z 1 λ ( u 3 + v 3 ) 2 z 3 λ ( u 1 + v 1 ) , d d t ( u 2 v 2 ) = 2 z 3 λ ( u 2 v 2 ) z 2 λ ( u 3 v 3 )
Let now
e 1 = z 1 2 2 λ b ( u 1 + v 1 ) , and e 2 = z 2 2 2 λ b ¯ ( u 2 v 2 ) .
Then
d e 1 d t = d d t ( z 1 2 2 λ b ( u 1 + v 1 ) ) = z 1 λ ( 1 λ z 1 z 3 + b ( u 3 + v 3 ) ) b ( z 1 λ ( u 3 + v 3 ) 2 z 3 λ ( u 1 + v 1 ) ) = 2 z 3 λ ( z 1 2 2 λ b ( u 1 + v 1 ) ) = 2 z 3 λ e 1 ,
and
d e 2 d t = d d t ( z 2 2 2 λ b ¯ ( u 2 v 2 ) ) = z 2 λ ( 1 λ z 2 z 3 b ¯ ( u 3 v 3 ) ) b ¯ ( 2 z 3 λ ( u 2 v 2 ) z 2 λ ( u 3 v 3 ) = 2 z 3 λ ( z 2 2 2 λ b ¯ ( u 2 v 2 ) ) = 2 z 3 λ e 2 .
Hence,  c = e 1 e 2  is an integral of motion for (61) since
d d t c = d d t e 1 e 2 = d e 1 d t e 2 + e 1 d e 2 d t = 0 .
An interested reader may want to show that the following are also integrals of motion
c 0 = | | V 2 | | , c 1 = | | U | | 2 , c 2 = U , V , c 3 = [ U , V ] · ( λ ( [ B , V ] + [ B , U ] ) z 3 Z ) + 1 2 ( ( V · Z ) 2 ( U · Z ) 2 ) .
The preceding calculation also draws attention to the following general fact:
Proposition 4. 
c = ( z 1 2 2 λ b ( u 1 + v 1 ) ) ( z 2 2 2 λ b ¯ ( u 2 v 2 ) )  is a constant of motion for any differential system in the variables  z i , u i , v i , i = 1 , 2 , 3  that satisfy
d z 1 d t = 1 λ z 1 z 3 + b ( u 3 + v 3 ) , d z 2 d t = 1 λ z 2 z 3 b ¯ u 3 v 3 ) , d d t ( u 1 + v 1 ) = z 1 λ ( u 3 + v 3 ) 2 z 3 λ ( u 1 + v 1 ) , d d t ( u 2 v 2 ) = 2 z 3 λ ( u 2 v 2 ) z 2 λ ( u 3 v 3 ) ,
independently of the equations that govern the evolution of  u 3  and  v 3 .
To come back to the top of Kowalewski, note that  V = 0  is an invariant subsystem for (60). On this set,  S = λ s B , and c reduces to
c = ( z 1 2 2 λ b u 1 ) ( z 2 2 2 λ b ¯ u 2 ) = ( z 1 2 2 λ b w 1 + s λ b 2 ) ( z 2 2 2 λ a ¯ b w 2 + s λ b ¯ 2 ) ,
and remains an integral of motion for the reduced system
d Z d t = [ Ω , Z ] + [ B , U ] , d U d t = [ Ω , U ] ) ,
with its fundamental relations (63)
d z 1 d t = 1 λ z 1 z 3 + b u 3 , d z 2 d t = 1 λ z 2 z 3 b ¯ u 3 , d d t u 1 = z 1 λ u 3 2 z 3 λ u 1 , d d t u 2 = 2 z 3 λ u 2 z 2 λ u 3 .
This reduced system coincides the Kirchhoff-Kowalewski system on  se ( 3 , C )  (Equation (49),  s = 0 , after u is replaced by w). Then
c = ( z 1 2 2 λ b u 1 ) ( z 2 2 2 λ b ¯ u 2 )
coincides with the integral of motion discovered by Kowalewski. The remaining isospectral integrals of motion  c 1 = | | U | | 2  and  c 3 = ( U · Z )  coincide with the Casimirs on  se ( 3 , C ) .
To recover the semi-simple form of the Kirchhoff-Kowalewski integral, let  Y = U + s 2 λ B . In terms of  Z  and  Y  the preceding system becomes
d Z d t = [ Ω , Z ] + [ B , Y ] , d Y d t = [ Ω , Y ] s λ 2 [ Ω , B ] .
This system satisfies the same equations as the Kirchhoff-Kowalewski system except for  d y 3 d t . Indeed,
s λ 2 [ Ω , B ] = s 2 ( [ Z , B ] + z 3 [ A 3 , B ] ) = s [ Z , B ] + s 2 ( [ Z , B ] + z 3 B ) = s [ B , Z ] + s 2 ( b ¯ z 1 b z 2 ) A 3 .
The remaining equations given by
d z 1 d t = 1 λ z 1 z 3 + b y 3 , d z 2 d t = 1 λ z 2 z 3 b ¯ y 3 , d y 1 d t = z 1 λ y 3 2 z 3 λ y 1 + s b z 3 , d y 2 d t = 2 z 3 λ y 2 z 2 λ y 3 s b ¯ z 3 ,
are the same as (66), and consequently yield
c = ( z 1 2 2 λ b y 1 + s 2 λ b 2 ) ( z 2 2 2 λ b ¯ y 2 + s 2 λ b ¯ 2 )
as an integral of motion for the preceding system, as well as for the Kirchhoff-Kowalewski system (Equation (49)) when  Y  is replaced by  W .
The papers of V. Dragovi c ´  and K. Kuki c ´  [30] and V. V. Sokolov [31] produce differential systems which admit Kowalewski type integrals different from the ones in this paper and yet follow the same integration procedure used by S. Kowalewski in her original paper. Remarkably, all these systems satisfy the fundamental relations (63) from which the existence of their extra integrals of motion could be easily ascertained.

4. Kepler, Jacobi, Neumann and Moser

Let us now return to  G = S L ( n + 1 )  and its Lie algebra  sl ( n + 1 )  endowed with the trace form  X , Y = 1 2 T r ( X Y ) . As a vector space V, the set of  ( n + 1 ) × ( n + 1 )  matrices with zero trace admits several kinds of Lie algebras and each of these Lie algebras induces its own Poisson structure on V. The most common Lie algebra is  sl ( n + 1 )  itself. Then  K = S O ( n + 1 )  induces the orthogonal decomposition  sl ( n + 1 ) = s y m 0 so ( n + 1 )  where  s y m 0  denotes the vector space of symmetric matrices in V. But then V also carries the semi-direct product structure  s y m 0 so ( n + 1 ) .
However,  K = S O ( p , q ) , p + q = n + 1 ,  is also a closed subgroup of G and hence the pair  ( S L ( n + 1 ) , S O ( p , q ) )  induces its own Cartan decomposition  V = p k , where  p  is the orthogonal complement to  k = so ( p , q ) . In fact K is the set of points in G fixed by the automorphism  σ ( g ) = D g T 1 D 1  where D denotes diagonal matrix with its first p diagonal entries equal to 1 and the remaining q diagonal entries equal to  1 .  The set of points  g G  such that  σ ( g ) = g  satisfies  D = g D g T , that is,  g S O ( p , q ) . It follows that its tangent map  σ *  induces the above decomposition with
k = { X V : D X D T = X } , p = { X V : D X T D = X } .
Consequently, matrices in  p  are symmetric relative to the scalar product  ( x , y ) p , q = ( x , D y ) x , y  in  R n + 1 .
We will now return to the canonical Hamiltonians
H ( L ) = 1 2 L k , L k + A , L p
and their Poisson Equation (33) restricted to the coadjoint orbits through rank one matrices  X 0  in  sl ( n + 1 ) . We will consider two cases: the coadjoint orbit through a symmetric rank-one matrix  X 0  of unit length under the action of  G 1 = s y m 0 S O ( n + 1 ) , and the second case, the coadjoint orbit through rank-one matrix  X 0  of unit length, symmetric relative to the Lorentzian inner product in  R n + 1  under the action of  p S O ( 1 , n ) . The above matrices can be naturally expressed in terms of the notations introduced earlier, the scalar product  ( x , y ) ϵ = x 0 y 0 + ϵ i = 1 n x i y i , ϵ = ± 1 ,  in the ambient space  R n + 1 , and matrices  a ϵ b  and  a ϵ b = a ϵ b b ϵ a . For then
X 0 = x 0 ϵ x 0 ( x 0 , x 0 ) ϵ n + 1 I .
If  ( x 0 , x 0 ) ϵ > 0  then let  S ϵ n = { x R n + 1 : ( x , x ) ϵ = ( x 0 , x 0 ) ϵ , x 0 . 0 } . It follows that  S ϵ n  is the Euclidean sphere of radius  | | x 0 | |  when  ϵ = 1  and a hyperboloid of two sheets when  ϵ = 1 . We have chosen  S 1 n  to be the sheet defined by  x 0 > 0 .
Proposition 5. 
The coadjoint orbit through  X 0 = x 0 ϵ x 0 ( x 0 , x 0 ) ϵ n + 1 I  is symplectomorphic to the cotangent bundle of the real projective space  P n + 1  in the semi-simple case, and it is symplectomorphic to the cotangent bundle of  S ϵ n  in the semi-direct case.
For the proof see [36]. Here it is implicitly understood that the cotangent bundles are identified with the tangent bundles via the ambient inner product  ( , ) ϵ . Then each tangent vector  ( x , y ) , x S ϵ n , ( x , y ) ϵ = 0  is identified with  L p = x ϵ x ( x 0 , x 0 ) ϵ n + 1 I  in  p ϵ  and  L k = x ϵ y  in  k ϵ .
On the orbit through  X 0 H = 1 2 ( x , x ) ϵ ( y , y ) ϵ 1 2 ( A x , x ) ϵ , and the associated Poisson equations are of the form
d d t ( x ϵ y ) = [ A , x ϵ x ] , d d t ( x ϵ x ) = [ x ϵ y , x ϵ x ]
A simple calculation show that
x ˙ = ( x , x ) ϵ y , y ˙ = A x ( A x , x ) ϵ ( x , x ) ϵ + ( y , y ) ϵ ) x .
On the unit sphere, Equation (71) after A is replaced by  A  coincide with the equations for the mechanical problem of C. Neumann for a particle on the sphere moving under a quadratic potential [37]. The preceding equations for  ϵ = 1  could be analogously interpreted as the equations on the hyperboloid for a particle moving under quadratic potential [7].
The canonical affine-quadratic problem illuminates deep and beautiful connections between Kepler’s gravitational problem, Jacobi’s geodesic problem on the ellipsoid, and Neumann’s mechanical problems.
Let us first examine the isospectral integrals associated with the spectral curve  L λ = L p λ L k + λ 2 A  on the coadjoint orbit through rank-one matrices. The zero trace requirement is inessential for the calculations below and will be disregarded. Additionally, A will be replaced by  A  and  L λ  will be rescaled by dividing by  λ 2  to read
L λ = 1 λ 2 L p + 1 λ L k + A . = 1 λ 2 x ϵ x + 1 λ x ϵ y + A .
The spectrum of  L λ  is then given by
0 = D e t ( z I L λ ) = D e t ( z I A ) D e t ( I ( z I A ) 1 ( 1 λ 2 L p + 1 λ L k ) ) ,
Matrix  M = I ( z I A ) 1 ( 1 λ 2 L p + 1 λ L k )  is of the form.
M = I + 1 λ 2 R z x ϵ x 1 λ ( R z x ϵ y R z y ϵ x ) ,
where  R z = ( z I A ) 1 .  We then have the following proposition
Lemma 1. 
D e t ( M ) = 1 λ 2 ( ( R z x , x ) ϵ + ( R z x , x ) ϵ ( R z y , y ) ϵ ( R z x , y ) ϵ 2 ) + 1 .
For the proof see [7] (p. 200).
Corollary 1. 
Function  F ( z ) = ( R z x , x ) ϵ + ( R z x , x ) ϵ ( R z y , y ) ϵ ( R z x , y ) ϵ 2 , z R  is an integral of motion for H.
Function F is a rational function with poles at the eigenvalues of the matrix  A .  Hence,  F ( z )  is an integral of motion for H if and only if the residues of F are constants of motion for H.
In the Euclidean case the eigenvalues of  A  are real and distinct since  A   is symmetric and regular. Hence there is no loss in generality in assuming that A is diagonal. Let  a 0 , , a n  denote its diagonal entries Then
F ( z ) = k = 0 n F k z a k ,
where  F 0 , , F n  denote the residues of  F .  It follows that
F ( z ) = k = 0 n x k 2 z a k + k = 0 n j = 0 n x k 2 y j 2 ( z a k ) ( z a j ) ( k = 0 n x k y k z a k ) 2 = k = 0 n x k 2 z a k + k = 0 n j = 0 , j k n x k 2 y j 2 ( z a k ) ( z a j ) 2 k = 0 n j = 0 , j k n x k y k x j y j ( z a k ) ( z a j ) .
Hence,
F k = lim z a k ( z a k ) F ( z ) = x k 2 + j = 0 , j k n x j 2 y k + x k 2 y j 2 ( a k a j ) 2 j = 0 , j k n x k y k x j y j ( a k a j ) = x k 2 + j = 0 , j k n ( x j y k x k y j ) 2 ( a k a j ) , k = 0 , , n .
The preceding calculation yields the following proposition.
Proposition 6. 
Each residue  F k = x k 2 + j = 0 , j k n ( x j y k x k y j ) 2 ( a k a j ) , k = 0 , , n  is an integral of motion for Newmann’s spherical system, and functions  F 0 , , F n  are in involution.
These results coincide with the ones reported in [38,39,40], but the connection with the affine-quadratic problem shows that similar integrals of motion exist for the hyperbolic Neumann problem [7] (p. 191).
In the literature on integrable systems the integrals of motion for Neumann’s problem are related to the integrals of motion for Jacobi’s problem on the ellipsoid through the transformation of H. Knörrer that transforms the Neumann’s equations on energy level  H = 0  onto the equations of Jacobi on the ellipsoid [39,41]. Our exposition takes another route: we will instead show that an “elliptic” problem on the sphere is completely integrable with its integrals of motion as in Neumann’s problem, and then we will show that the Hamiltonian equations for the elliptic problem on the sphere and Jacobi’s problem on the ellipsoid are symplectomorphic. We will then use this symplectomorphism to show the existence of Jacobi’s integrals of motion on the ellipsoid.
Let now  H = 1 2 D 1 ( L k ) D 1 , L k + D 1 , L p  denote an affine-quadratic Hamiltonian on  G = S L ( n + 1 )  defined by a diagonal matrix D with positive diagonal entries. As before, we will dispense with zero-trace requirements since they are inessential. The above Hamiltonian is generated by a positive definite operator  P ( X ) = D X D , X so ( n + 1 )  and the drift  A = D 1 . We will call this Hamiltonian elliptic for reasons that will be made clear later on. Let  B = D . Then
[ P 1 ( L k ) , B ] = [ P 1 ( L k ) , D ] = L k D 1 + D 1 L k = [ L k , D 1 ] = [ L k , A ] .
Therefore  H  is isospectral (Proposition 1) and its Hamiltonian equations admit a representation
d L λ d t = [ M λ , L λ ] , L λ = L p λ L k ( λ 2 s ) A .
Since  L λ = L p λ L k ( λ 2 s ) A  is a spectral curve for the canonical affine Hamiltonian  H = 1 2 L k , L k A , L p  we have the following corollary.
Corollary 2. 
The spectral invariants of  L p λ L k ( λ 2 s ) A  are common integrals of motion for both the canonical Hamiltonian  H = 1 2 L k , L k A , L p  and the elliptic Hamiltonian  H = 1 2 A 1 L k A 1 , L k + A 1 , L p .
As before, on coadjoint orbit through  X 0 = x 0 ϵ x 0 , ( x 0 , x 0 ) ϵ = 1 , the Poisson equations of  H  ( the semi-direct version) are given by
d d t ( x ϵ y ) = [ A 1 ( x ϵ y ) A 1 , x ϵ y ] + [ A 1 , x ϵ x ] , d d t ( x ϵ x ) = [ A 1 ( x ϵ y ) A 1 ) , x ϵ x ] .
We then have
A 1 L k A 1 , L k = A 1 x y A 1 , x y = ( A 1 x · x ) ( A 1 y · y ) ( A 1 x · y ) 2 , A 1 , L p = 1 2 ( x · A 1 x ) .
which shows that the Hamiltonian  H = 1 2 A 1 L k A 1 , L k + A 1 , L p  is given by
H = 1 2 ( ( A 1 y · y ) ( A 1 x · y ) 2 ( A 1 x · x ) 1 ) ( A 1 x · x ) .
The correspondence  ( x , y ) x ϵ y + x ϵ x  defines a symplectomorphism between the cotangent bundle of  S ϵ n  with its canonical Poisson bracket and the coadjoint orbit through  X 0  (Proposition 5).
Proposition 7. 
On energy level  H = 0  Equation (72) correspond to
d x d t = ( A 1 x · x ) ( A 1 y ( A 1 x · y ) ( A 1 x · x ) A 1 x ) d y d t = ( A 1 x · x ) ( ( A 1 x · y ) ( A 1 x · x ) A 1 y ( A 1 x · y ) 2 ( A 1 x · x ) 2 A 1 x x . )
under the correspondence  ( x x , x y ) ( x , y ) .
These equations can be reparametrized by a parameter  s = ( A 1 x ( t ) · x ( t ) ) d t  to read
d x d s = d x d t d t d s = A 1 y ( A 1 x · y ) ( A 1 x · x ) A 1 x d y d s = d y d t d t d s = ( A 1 x · y ) ( A 1 x · x ) ( A 1 y ( A 1 x · y ) ( A 1 x · x ) A 1 x ) x .
We will presently show that Equation (73) are Hamiltonian equations that correspond to the geodesic problem on the sphere relative to the elliptic metric  1 2 ( A x ˙ , x ˙ ) .
As an intermediate step we will now derive the Hamiltonian equations associated with the geodesic problem on the quadric surface  ( A 1 x , x ) = 1  induced by the scalar product  1 2 ( D x ˙ , x ˙ ) . We will follow the procedure based on the version of the Maximum Principle for variational problems with constraints outlined in [7] (p. 218) and identify the quadric surface with the submanifold  N = { x R n + 1 : ( x , A 1 x ) = 1 } . Then its cotangent bundle will be defined in terms of the constraints
G 1 = { ( x , A 1 x ) 1 = 0 } , and G 2 = { ( x , A 1 p ) = 0 } .
The Hamiltonian lift of a curve  x ˙ = u ( t )  that belongs to  T x ( t ) N  is given by
h u ( t ) ( x , p ) = 1 2 ( D u , u ) + ( p , u ) + λ 1 G 1 + λ 2 G 2
for the multipliers  λ 1  and  λ 2  that satisfy  { h u ( t ) , G 1 } = { h u ( t ) , G 2 } = 0 . If  h u 0 = 1 2 ( D u , u ) + ( p , u ) , then
{ h u , G 1 } = { h u 0 , G 1 } + λ 2 { G 2 , G 1 } , { h u , G 2 } = { h u 0 , G 2 } + λ 1 { G 1 , G 2 } .
It follows that
{ h u 0 , G 1 } = 2 u · A 1 x , { h u 0 , G 2 } = u · A 1 p , { G 1 , G 2 } = 2 A 1 x · A 1 x .
Hence,
λ 1 = 1 { G 1 , G 2 } { h u 0 , G 2 } = 1 2 ( u , A 1 p ) ( A 1 x , A 1 x ) , λ 2 = 1 { G 2 , G 1 } { h u 0 , G 1 } = ( u , A 1 x ) ( A 1 x , A 1 x ) .
According to the Maximum Principle an extremal control  u ( t )  must optimize
h u 0 = 1 2 ( u · D u ) + p · u
on  G 1 = G 2 = 0  over all controls that satisfy  u · A 1 x = 0 . Hence extremal controls are the critical points of  1 2 ( u , D u ) + ( p , u ) α 0 ( u , A 1 x )  for some multiplier  α 0 , that is, they are solutions of  D u + p α 0 A 1 x = 0 . It follows that the extremal controls are of the form  u = D 1 ( p α 0 A 1 x ) , But then  ( u , A 1 x ) = 0  implies that  α 0 = ( A 1 x , D 1 p ) ( D 1 A 1 x , A 1 x ) .  For this choice of controls
h u = H + λ 1 G 1 + λ 2 G 2 ,
where  H = 1 2 ( D 1 ( p α 0 A 1 x ) , p ) ) . An easy calculation shows that
H = 1 2 ( ( D 1 p , p ) ( D 1 p , A 1 x ) 2 ( A 1 x , A 1 x ) ) .
Then the extremal curves are the solutions of the following differential equation:
d x d t = H p = D 1 ( p α 0 A 1 x ) , d p d t = H x 2 λ 1 G 1 x = α 0 ( A 1 ( D 1 ( p α 0 A 1 x ) ) 2 λ 1 A 1 x ,
which emanate from  H = 1 2 , that is, satisfy
( D 1 p , p ) ( D 1 A 1 x , A 1 x ) ( D 1 p , A 1 x ) 2 = ( D 1 A 1 x , A 1 x ) .
We will now single out the cases relevant for our earlier claims.
The geodesic problem on the ellipsoid. In this classic case initiated by C. Jacobi  D = I  and  ( A 1 x , p ) = 0 . Hence
α 0 = ( A 1 x , p ) ( A 1 x , A 1 x ) = 0 , λ 1 = 1 2 ( p , A 1 p ) ( A 1 x , A 1 x ) .
Then Equation (74) reduce to
d x d t = p , d p d t = ( p , A 1 p ) ( A 1 x , A 1 x ) A 1 x .
The preceding equation agree with the equations in J. Moser [39].
The elliptic problem on the sphere. Here the ambient metric is defined by a positive-definite matrix D and  A = I . In such a case Equation (74) are given by
( D 1 p , p ) ( D 1 p , x ) 2 ( D 1 x , x ) = 1 .
Furthermore,
λ 1 = 1 2 D 1 ( p α 0 x ) , p ) = 1 2 ( D 1 p , p ) ( x , D 1 p ) 2 ( D 1 x , x ) ) = 1 2 , α 0 = ( D 1 x , p ) ( D 1 x , x ) .
The Hamiltonian equations are then given by
d x d t = D 1 p ( D 1 x , p ) ( D 1 x , x ) D 1 x , d p d t = ( D 1 x , p ) ( D 1 x , x ) ( D 1 p ( D 1 x , p ) ( D 1 x , x ) D 1 x ) ) x ,
which agrees with Equation (73) when  D = A .
Proposition 8. 
The Hamiltonian systems that correspond to the elliptic problem on the sphere and the geodesic problem on the ellipsoid are symplectomorphic.
Proof. 
Let  ( x , y )  denote the coordinates on the tangent bundle of the sphere and let  ( q , p )  denote the coordinates on the tangent bundle of the ellipsoid  E = { q R n + 1 : ( q , A 1 q ) 1 = 0 } . In these coordinates the systems in question are given by
d x d t = u , d y d t = α u x , and d q d t = p , d p d t = ( A 1 p , p ) ( A 1 q , A 1 q ) A 1 q ,
where  u = A 1 ( y α x )  and  α = ( A 1 x , y ) ( A 1 x , A 1 x ) .
Let  Φ  denote the mapping from the cotangent bundle of the sphere to the cotangent bundle of E defined by
q = A 1 2 x , p = A 1 2 ( y α x ) = A 1 2 u .
Let  θ = i = 0 n p i d q i = ( p , d q ) , ( d q , A 1 q ) = 0 ,  denote the Liouville-Poincaré canonical form on  T * E . Then
Φ * θ = ( A 1 2 ( y α x ) , A 1 2 d x ) = ( y , d x ) α ( x , d x ) = ( y , d x ) ,
because  0 = d q · A 1 q = A 1 2 d x · A 1 A 1 2 x = d x · x . Since  Φ *  takes the Liouville form on  T * E  to the Liouville form on  T * S n , it also takes the canonical symplectic form on  T * E  to the canonical symplectic for on  T * S n  and hence is a symplectomorphism.
It now follows from (78) that  d u d t = ( 1 + d α d t ) A 1 x  and that  1 + d α d t = ( u , u ) ( A 1 x , x ) . Then,
d q d t = A 1 2 d x d t = A 1 2 u = p , d p d t = A 1 2 d u d t = ( ( u , u ) ( A 1 x , x ) ) A 1 q = ( A 1 p , p ) ( A 1 q , A 1 q ) A 1 q ,
and thus  Φ *  takes the Hamiltonian flow on the sphere onto the Hamiltonian flow on the ellipsoid. □
Proposition 9. 
Jacobi’s problem on the ellipsoid is completely integrable. Functions
G k = p k 2 + j = 1 , j k n + 1 ( q j p k q k p j ) 2 ( a k a j ) , k = 1 , , ( n + 1 )
are constants of motion, all in involution with each other, for the Hamiltonian system  d q d t = p , d p d t = ( A 1 p , p ) ( A 1 q , A 1 q ) A 1 q  on the cotangent bundle of the ellipsoid.
Proof. 
We have shown that
F k = x k 2 + j = 0 , j k n ( x j y k x k y j ) 2 ( α k α j ) , k = 0 , , n
are an involutive family of integrals of motion for the elliptic-geodesic problem on the sphere. We have also shown that the above integrals of motion are the residues of the function
F ( z ) = ( R z x , x ) + ( R z x , x ) ( R z y , y ) ( R z x , y ) 2 , R z = ( z I A ) 1 .
We will now show that functions (79) are the residues of the pull-back of F under the symplectomorphism  Φ . First note that F remains unchanged if the variable y is replaced by  y + α x  with  α  an arbitrary number. Since  A 1 2 p = y ( A 1 x , y ) ( A 1 x , x ) x , we may replace y by  A 1 2 p  and x by  A 1 2 q . Also note that  ( p , p ) = ( A 1 y , y ) ( A 1 x , y ) 2 ( A 1 x , x ) = 1  (use Equation (75)). Then,
1 + ( R z y , y ) = 1 + ( R z A p , p ) = 1 + k = 0 n a k p k 2 z a k = k = 0 n p k 2 + a k p k 2 z a k = z k = 0 n p k 2 z a k = z ( R z p , p ) , ( R z x , x ) = ( R z A 1 q , q ) = k = 0 n q k 2 a k ( z a k ) = 1 z k = 0 n q k a k + q k 2 z a k = 1 z ( 1 + ( R z q , q ) , and ( R z x , y ) = ( R z q , p ) .
It follows that
F ( z ) = ( R z p , p ) ( 1 + R z q , q ) ( R z q , p ) 2
is constant along the solutions of Jacobi’s equations.. A calculation identical to the one used for Neumann’s system shows that
G k = p k 2 + j = 0 , j k n ( q j p k q k p j ) 2 ( a k a j ) , k = 1 , , ( n + 1 )
are the residues of F, and hence are integrals of motion for Jacobi’s equations. □

Degenerate Case  A = 0  and Kepler’s Problem

Let us now return to the Hamiltonian equations generated by the canonical affine Hamiltonian  H = 1 2 ( x , x ) ϵ ( y , y ) ϵ 1 2 ( A x , x ) ϵ  on the coadjoint orbit through  X 0 = a ϵ a  for some  a R n + 1  with
d d t ( x ϵ y ) = [ A , x ϵ x ] , d d t ( x ϵ x ) = [ x ϵ y , x ϵ x ] ,
and their equivalent formulation on the tangent bundle of  S ϵ n :
x ˙ = ( x , x ) ϵ y , y ˙ = A x ( A x , x ) ϵ ( x , x ) ϵ + ( y , y ) ϵ ) x .
When  A = 0  the Hamiltonian H reduces to  H = 1 2 ( x , x ) ϵ ( y , y ) ϵ  and the corresponding equations reduce to
d d t ( x ϵ y ) = 0 , d d t ( x ϵ x ) = [ x ϵ y , x ϵ x ] .
Then Equation (82) yield an integral of motion  x ϵ y = c o n s t , and Equation (81) reduce to
x ˙ = | | x | | ϵ 2 y , y ˙ = | | y | | ϵ 2 x , where | | x | | ϵ 2 = ( x , x ) ϵ , | | y | | ϵ 2 = ( y , y ) ϵ .
Upon differentiating we get
x ¨ + | | x | | ϵ 2 | | y | | ϵ 2 x = 0 .
We will now assume that  ( a , a ) ϵ = h 2  so that  S 1 n  is the hyperboloid  x 0 2 = h 2 + i = 1 n x i 2 , x 0 0 . On energy level  H = ϵ 2 h 2 | | x | | ϵ 2 | | y | | ϵ 2 = ϵ h 2  the solutions of (83) are given by  x ( t ) = c 1 cos t h ϵ + c 2 sin t h ϵ  where  c 1  and  c 2  are constant vectors (complex when  ϵ = 1 ) that satisfy  | | c | | ϵ 2 = | | c | | ϵ 2 = h 2 , ( c 1 , c 2 ) ϵ = 0 .
For  ϵ = 1  the above curves trace great circles on the sphere  | | x | | 2 = | | x 0 | | 2  and for  ϵ = 1  the solutions trace great hyperbolas on the hyperboloid  | | x | | 1 2 = | | x 0 | | 1 2  (an immediate consequence of the fact that  x ( t ) ϵ x ˙ ( t ) = a ϵ b ). That is, solutions are the geodesics on spaces of constant non-zero curvature. The zero curvature case may be obtained by considering  ϵ  as a continuous parameter and then letting it tend to zero (as will be explained below).
We will now show that there exists a canonical change of coordinates  { ( x 0 , , x n , y 0 , , y n ) , ( x , x ) ϵ = h 2 , ( x , y ) ϵ = 0 } ( p 1 , , p n , q 1 , , q n )  in which p is the stereographic projection through the point  x 0 = h e 0  given by
λ ( x h e 0 ) + h e 0 = ( 0 , p ) with λ = h h x 0 , x S ϵ n
such that in the new coordinates the preceding geodesic differential system is transformed into the n-dimensional Kepler’s system, an n-dimensional generalization of the Hamiltonian equations that describe the motion of a planet around an immovable planet in the presence of the gravitational force.
Equation (84) yields  p = h h x 0 x ¯ , where  x ¯ = ( x 1 , , x n ) . Therefore the inverse map  x = Φ ϵ ( p )  is given by
x 0 = h ( | | p | | 2 ε h 2 ) | | p | | 2 + ε h 2 , x ¯ = 2 ϵ h 2 | | p | | 2 + ε h 2 p .
Assume that the cotangent bundle of  R n  is identified with its tangent bundle  R n × R n  via the Euclidean inner product  ( · ) , and let  ( p , q )  denote the points of  R n × R n . We will next find  q = Ψ ( x , y )  such that  ( d x , y ) ϵ = ( d p · Ψ ( x , y ) ) ,  for all  ( x , y )  with  x S ϵ n ( h )  and  ( x , y ) ϵ = 0 . For then the transformation  ( p , q ) R n × R n ( x , y ) T S ϵ n ( h )  is a symplectomorphism since it pulls back the Liouville form  ( d x , y ) ϵ  on  T S ϵ n ( h )  onto the Liouville form  ( d p · q )  in  R n × R n  (The symplectic form is the exterior derivative of the Liouville-Poincaré form). It follows that
( y , d x ) ϵ = j = 1 n y 0 x 0 p j d p j + ϵ i = 1 n j = 1 n y i x i p j d p j = j = 1 n ( y 0 x 0 p j + ϵ i = 1 n y i x i p j ) d p j = j = 1 n q j d p j = ( d p · q ) .
Therefore,
q j = x 0 p j y 0 + ϵ i = 1 n y i x i p j , j = 1 , , n .
After the appropriate differentiations in (85) we get
q = 2 h 2 | | p | | 2 + ϵ h 2 ( 2 ϵ h y 0 | | p | | 2 + ϵ h 2 p + y ¯ 2 ( y ¯ · p ) | | p | | 2 + ϵ h 2 p ) , y ¯ = ( y 1 , , y n ) .
Hence,
| | p | | 2 + ϵ h 2 2 h 2 ( q · p ) = 2 ϵ h y 0 | | p | | 2 | | + ϵ h 2 | | p | | 2 | | p | | 2 ϵ h 2 | | p | | 2 + ϵ h 2 ( y ¯ · p ) .
Since y is orthogonal to x ( y ¯ · p ) = y 0 2 h ( | | p | | 2 ϵ h 2 ) . Therefore,
y 0 = 1 h q · p , y ¯ = | | p | | 2 + ε h 2 2 h 2 q q · p h 2 p .
After the substitutions  | | p | | 2 = h 2 h x 0 | | x ¯ | | 2  into the preceding equation we get
q = 2 ( h x 0 ) | | x ¯ | | 2 + ( h x 0 ) ( ( h x 0 ) 2 y ¯ + y 0 x ¯ ) .
To pass to the problem of Kepler, write the Hamiltonian  H = 1 2 | | x | | ε 2 | | y | | ε 2  in the variables  ( p , q ) . An easy calculation in (86) yields  ( y , y ) ϵ = ε ( | | p | | 2 + ε h 2 ) 2 4 h 4 | | q | | 2 . Therefore,
H = 1 2 h 2 ε ( | | p | | 2 + ε h 2 ) 2 4 h 4 | | q | | 2 = 1 2 ε ( | | p | | 2 + ε h 2 ) 2 4 h 2 | | q | | 2 .
The corresponding flow is given by
d p d s = H q = ε ( | | p | | 2 + ε h 2 ) 2 4 h 2 q , d q d s = H p = ε | | p | | 2 + ε h 2 2 h 2 | | q | | 2 p
On energy level  H = ε 2 h 2 , ( | | p | | 2 + ε h 2 ) 2 4 | | q | | 2 = 1 ,  and the preceding equations reduce to
d p d s = ε q h 2 | | q | | 2 , d q d s = ε | | q | | h 2 p .
After the reparametrization  t = ε h 2 0 s | | q ( τ ) | | d τ  Equation (23) become
d p d t = d p d s d s d t = q | | q | | 3 , d q d t = d q d s d s d t = p .
On  H = ϵ h 2 ( | | p | | 2 + ε h 2 ) 2 4 | | q | | 2 = 1  and
E = 1 2 | | p | | 2 1 | | q | | = 1 2 | | q | | ( | | p | | 2 | | q | | 2 ) = 1 2 | | q | | ( 2 ε h 2 | | q | | 2 ) = 1 2 ε h 2 .
So  E < 0  in the spherical case and  E > 0  in the hyperbolic case.
The Euclidean case  E = 0  can be obtained by a limiting argument in which  ε  is regarded as a continuous parameter which tends to zero.
To explain in more detail, let
w 0 = lim ϵ 0 = h , w = lim ϵ 0 1 2 ϵ h 2 x ¯ = 1 | | p | | 2 .
The transformation  p w  with  w = 1 | | p | | 2 p  is the inversion about the circle  | | p | | 2 = 1  in the affine hyperplane  w 0 = h , and  | | d w | | 2 = 1 | | p | | 4 | | d p | | 2  is the corresponding transformation of the Euclidean metric  | | d p | | 2 . The Hamiltonian  H 0  associated with this metric is equal to  1 2 | | p | | 4 4 | | q | | 2 .  This Hamiltonian can be also obtained as the limit of  ( h 2 ϵ ) 1 2 ( | | p | | 2 + ϵ h 2 ) 2 4 h 2 | | q | | 2  when  ϵ 0 . On energy level  H = 1 2 , | | p | | 2 | | q | | = 2  and therefore,  E = 0 . Of course, the solutions of (12) tend to the Euclidean geodesics as  ϵ  tends to zero. Consequently,  w ( t ) = lim ϵ 0 1 2 h 2 ε ( x ¯ ( t ) )  is a solution of  d 2 w d t 2 = 0 ,  and hence, is a geodesic corresponding to the standard Euclidean metric.
Let us also note that the angular momentum  L = q p  and the Laplace-Runge-Lenz vector  F = L p q | | q | |  for Kepler’s problem have simple geometric interpretation on the coadjoint orbits according to the following proposition.
Proposition 10. 
Let  x = x 0 e 0 + x ¯  and  y = y 0 e 0 + y ¯ .  On energy level  H = ϵ 2 h 2 ,
L = ( y ¯ ϵ x ¯ ) and F = h ( y 0 ( e 0 x ¯ ) ε x 0 ( e 0 y ¯ ) ε ) e 0 .
For a proof see [7].
This remarkable discovery that the solutions of Kepler’s problem are intimately related to the geometry of spaces of constant curvature goes back to A.V. Fock’s paper of 1935 [42] in which he reported that the symmetry group for the motions of the hydrogen atom is  O 4 ( R )  for negative energy,  E 3 O 3 ( R )  for zero energy and  O ( 1 , 3 )  for positive energy. It is then not altogether surprising that similar results apply to the problem of Kepler since the energy function for Kepler’s problem is formally the same as the energy function for the hydrogen atom.
This connection between the problem of Kepler and the geodesics on the sphere was reported by J. Moser in 1970 [43], while Y. Osipov [44] reported similar results later for geodesics on spaces of negative constant curvature. In spite of their brilliance, these papers did not attempt any explanations in regard to this enigmatic connection between planetary motions and geodesics on space forms. This issue later inspired V. Guillemin and S. Sternberg to take up the problem of Kepler in a larger geometric context, with Moser’s observation as the background, in a paper titled Variations on a theme by Kepler [45]. The introduction of Kepler’s problem through the canonical affine-quadratic problem exemplifies, once again, this fascinating and recurrent interplay between mathematical physics, geometry and integrable systems.

5. Homogeneous Riemannian Manifolds and Rolling Geodesics

Our overview of integrable systems raises a natural question: what is the geometric origin behind the affine-quadratic problem that accounts for its ubiquitous presence in the theory of integrable systems? A partial answer to this question comes, somewhat unexpectedly, from a new class of variational problems, called rolling problems. We will take up this issue next. Since the underlying variational problems require new concepts and terminology, we will be obliged to make a slight detour into an earlier paper [9] in order to introduce the necessary ingredients.
The general setting is defined by a reductive pair  ( G , K )  with G semi-simple and K compact. We assume that the Lie algebra decomposition  g = p k , with  p  the orthogonal complement of  k  relative to the Killing form on  g  satisfies the strong Cartan conditions
[ p , k ] = p , [ p , p ] = k , [ k , k ] k .
We will also assume the Killing form is of definite sign on  p  in which case  ,  will denote a scalar multiple of the Killing form that is positive on  p . We recall that the Killing form is invariant under any linear automorphism of  g  and hence the quadratic form  ,  is  A d G  invariant [15].
We consider G a semi-Riemannian manifold relative to the left-invariant metric  g X , g Y g = X , Y , X , Y g  induced by  ,  (the Killing form is not necessarily positive on  g , hence the metric is in general of indefinite sign, i.e., it is semi-Riemannian [46]). The left-invariant distributions  D ( g ) = { g X : X p }  and  V ( g ) = { g X ) : X k }  are called horizontal and vertical respectively. Then curves that are tangent to  D , i.e., satisfy  d g d t D ( g ( t ) )  are called horizontal. Likewise curves that are tangent to  V  are called vertical. It follows that
D ( g ) V ( g ) = T g G , g G .
We will assume that  M = G / K  consisting of the left coset  g K  is endowed with a manifold structure so that the natural projection  π ( g ) = g K  is a smooth surjection [46]. A curve  g ( t )  in G is called a lift of a curve  p ( t ) M  if  π ( g ( t ) ) = p ( t ) . A lift is called horizontal when  g ( t )  is a horizontal curve. Every curve  p ( t )  in M is the projection of a horizontal curve  g ( t ) . If a curve  g ( t )  is a solution of  d g d t = g ( t ) U ( t )  for some curve  U ( t ) p  then  d g ( t ) π ( g ( t ) U ( t ) ) = d p d t .  The correspondence  D ( g ) T π ( g ) M  given by  d g π ( g U ) = d p d t  is an isomorphism and induces a metric on M
( d g π ( g V ) , d g π ( g W ) π ( g ) = g V , g W g = V , W , V , W p .
Let now  { τ g : g G }  denote the group of diffeomorphisms defined by the left action
π ( L g ( h ) ) = τ g ( π ( h ) ) , h G , L g ( h ) = g h .
We then have
Proposition 11. 
The metric (89) is invariant under  { τ g : g G } , that is,
( d o τ g ( V ( p ) , d o τ g ( W ( p ) ) τ g ( p ) = ( V ( p ) , W ( p ) ) p ,
for any  g G  and any tangent vectors  V ( p )  and  W ( p )  in  T p M .
For a proof see [47].
It follows that each  τ g  is an isometry. Since G acts transitively on M, M can be represented by the orbit  { τ g ( o ) : g G }  where  o = π ( e )  and e is the group identity in G. It follows that  π ( ( exp t U ) g ) = τ exp t U ( π ( g ) )  for any  U g . Note that  g ( exp t U ) g  is the flow generated by a right-invariant vector field  U r ( g ) = U g . Therefore the flow of  U r  is  π -related to the flow  { τ exp t U : t R }  in M. We will let  U  denote the infinitesimal generator of the flow  { τ exp t U : t R } .
It follows that each  U  is a Killing vector field on M. A vector field whose flow acts on M by isometries is called a Killing vector field (see [46] for additional details). The correspondence  U r ( g ) U ( π ( g ) )  is one to one and onto  T π ( g ) M . Since the Lie brackets of vector fields related by a mapping F are also F-related, the Lie brackets  [ U r , V r ]  are  d π -related to  [ U , V ] . Therefore the correspondence  U r ( g ) U ( π ( g ) )  is a Lie algebra homomorphism, and hence  F = { U : U g }  is a finite dimensional Lie algebra of Killing vector fields that satisfies  F ( p ) = T p M  for each  p M .
Note that  π ( exp t U ) = τ e t U ( o ) = exp t U ( o ) . So if  U k  then  π ( exp t U ) = o  and therefore  U ( o ) = 0 . It then follows that  d e π ( U ) = U ( o )  is an isometric isomorphism from  p  onto  T o M . More generally if  g ( t )  is any horizontal curve then  p ( t ) = π ( g ( t ) ) = τ g ( t ) π ( e )  implies that
d p d t = d g ( t ) ( π ( g ( t ) U ( t ) ) ) = d o τ g ( t ) d e π ( U ( t ) ) = d o τ g ( t ) U ( t ) ( o ) ,
and
( d o τ g ( t ) U ( t ) ( o ) , d o τ g ( t ) V ( t ) ( o ) ) p ( t ) = ( U ( o ) , V ( o ) ) o .
Therefore  d o τ g ( t )  is an isometry that maps  T o M  onto  T p ( t ) M .
A homogeneous manifold  M = G / K  with a G-invariant metric defined by a reductive pair (G,K) with G semi-simple and K compact, will be referred to as semi-simple (it is defined by a semi-simple Lie group G, a compact subgroup K, and the metric induced by the Killing form). It can be shown that any symmetric Riemannian space with no Euclidean factors can be reduced to a semi-simple manifold (so that  [ p , p ] = k  holds). Conversely, every semi-simple manifold is locally symmetric. It is symmetric when G is simply connected (see [48], Proposition 6.27). We will not pursue further proximities with symmetric spaces since the present exposition makes no use of geodesic symmetries.
We now come to the main topic of this section, rolling of semi-simple manifolds on their tangent spaces. We begin by recalling the basic definition.
Definition 1. 
A curve  α ( t )  on a Riemannian manifold M rolls on a curve  α ^ ( t )  on another Riemannian manifold  M ^  if there exists an isometry  A ( t ) : T α ( t ) M T α ^ ( t ) M ^  that satisfies:
d α ^ d t = A ( t ) d α d t ,
and also satisfies the condition that  A ( t ) v ( t )  is a parallel vector field in  M ^  along  α ^ ( t )  for each parallel vector field  v ( t )  along  α ( t )  in M.
This intrinsic definition of rolling was introduced in [49], and later used in [50,51]. In this context the triple  ( α ( t ) , α ^ ( t ) , A ( t ) )  is called a rolling curve. It is clear that rolling is reflexive in the sense that if  α ( t )  is rolled on  α ^ ( t )  by an isometry  A ( t )  then  α ^ ( t )  is rolled on  α ( t )  by the isometry  A 1 ( t ) , and therefore  ( α ^ ( t ) , α ( t ) , A 1 ( t ) )  is also a rolling curve.
We will now examine rollings of semi-simple manifolds on their tangent planes. It comes as a pleasant surprise that such rollings are essentially described by Equation (91) reinterpreted in terms of rolling. So the passage to rolling becomes largely a question of semantics, as demonstrated in the text below.
We will consider rollings of M on  M ^ = T o M  with its metric  ( u , v ) o  defined by (89). The rollings on other tangent spaces are conjugate to the rollings on  T o M  [47]. Let  α ( t )  be an arbitrary curve in M and let  α ^ ( t )  be a curve in  M ^  that  α ( t )  is rolled on. It follows that  α ( t ) = π ( g ( t ) ) = τ g ( t ) ( o )  for some horizontal curve  g ( t ) . If  g ( t )  is a solution of  d g d t = g ( t ) U ( t ) , U ( t ) p  then according to (91)
d α ( t ) d t = d g ( t ) π ( g ( t ) U ( t ) ) = d o τ g ( t ) U ( t ) ( o ) ,
If we now let  α ^ ( t )  be any solution in  M ^  of  d α ^ ( t ) d t = U ( t ) ( o )  then  A ( t ) = d o τ g ( t )  is an isometry that rolls  α ^ ( t )  on  α ( t )  since the parallel transport condition is satisfied (for proofs see [47]). Of course, then  A 1 ( t )  rolls  α ( t )  on  α ^ ( t ) .
It follows that each horizontal curve  g ( t )  in G defines a family of curves  α ^ ( t )  in  M ^ , each a solution of  d α ^ d t = U ( t ) ( o ) , with  U ( t )  induced by  U ( t ) = g 1 ( t ) d g d t , that roll on  α ( t ) = π ( g ( t ) ) . The converse is also true: every solution  ( g ( t ) , α ^ ( t ) )  of the differential system
d g d t = g ( t ) U ( t ) , d α ^ ( t ) d t = U ( t ) ( o ) , U ( t ) p
defines a curve  α ( t ) = π ( g ( t ) )  in M on which  α ^ ( t )  in  M ^  is rolled by the isometry  d o τ g ( t ) .
We will regard (93) as the fundamental object in rolling defined on  G = G × M ^ , M ^ = T o M , a Lie group with its group operation
g h = ( g , p ) ( h , q ) = ( g h , p + q ) , g = ( g , p ) , h = ( h , q ) .
Then  G = g × T o M  will denote the Lie algebra of  G  with the Lie bracket  [ ( X , U ( o ) ) , ( Y , V ( o ) ] = ( [ X , Y ] , 0 ) .
Let now  H ( g , p ) = { ( g U , U ( o ) ) : U p } , ( g , p ) G . We will view  H  as a left-invariant distribution on  G  defined by the left-translates of vector space  Γ = { ( U , U ( o ) ) : U p }  in  G . The distribution  H  is called the rolling distribution and its integral curves are called rolling motions. Any rolling motion  g ( t ) = ( g ( t ) , p ( t ) )  is a solution of
d g d t = g ( t ) U ( t ) , d p d t = U ( t ) ( o ) ,
and can be associated with the rolling curve  ( α ^ ( t ) , α ( t ) ) , d o τ g ( t ) ) , where  α ( t ) = τ g ( t ) ( o )  and  α ^ ( t ) = p ( t ) .
Since  p  and  k  satisfy strong Cartan conditions  [ p , k ] = p  and  [ p , p ] = k Γ  satisfies  [ Γ , Γ ] = ( k , 0 ) ,  and  [ Γ , [ Γ ] , Γ ] ] = ( p , 0 ) . Therefore,
Γ + [ Γ , Γ ] + [ Γ , [ Γ , Γ ] ] = G ,
Hence the Lie algebra generated by the left-invariant vector fields tangent to  H  is equal to  G , therefore any two points in  G  can be connected by a rolling motion, and each rolling motion inherits a natural length  0 T U ( t ) , U ( t ) d t  from G. It is then known that any pair of points in  G  can be connected by an integral curve of  H  of minimal length because vector fields in  H  are complete [50]. The above shows that  G  with the above metric is a sub-Riemannian manifold. We will refer to the associated sub-Riemannian geodesics as the rolling geodesics.
We will now turn to the Maximum principle to find the necessary conditions that the rolling geodesics must satisfy. To put the matter in the control theoretic context, let  A 1 , , A m  be an orthonormal basis in  p  so that  ( A i , A i ( o ) )  becomes an orthonormal basis in  Γ . Then an absolutely continuous curve  g ( t ) = ( g ( t ) , p ( t ) )  is a rolling motion if and only if
d g d t = i = 1 m u i ( t ) g ( t ) A i , d p d t = i = 1 m u i ( t ) A i ( o ) ,
for some bounded and measurable control functions  u 1 ( t ) , , u m ( t ) , in which case the length of  g ( t )  is given by  0 T u 1 2 ( t ) + + u m 2 ( t ) d t . The rolling problem is an optimal control problem and consists of finding the solutions  g ( t ) = ( g ( t ) , p ( t ) )  on a fixed time interval  [ 0 , T ]  that satisfy the given boundary conditions  g ( 0 ) = g 0  and  g ( T ) = g 1  along which the energy of transfer  1 2 0 T i = 1 m u i 2 ( t ) d t  is minimal. It is known that each rolling geodesic is locally optimal and hence is a solution to the above control problem [7,50].

5.1. Rolling Hamiltonians

To emphasize the invariant properties of the problem we will rewrite (96) as
d g d t = i = 1 m u i ( t ) X i ( g ) ,
where each  X i  a left-invariant vector field  X i ( g ) = ( g A i , A i ( o ) ) g = ( g , p ) . If  g ( t )  is an optimal trajectory then, according to the Maximum Principle,  g ( t )  is the projection of an extremal curve  ξ ( t )  in  T * G  along which the cost extended Hamiltonian
λ 2 i = 1 m u i 2 ( t ) + i = 1 m u i ( t ) H i ( ξ ( t ) ) , λ = 0 , 1
is maximal relative to all other control functions. Here  H i  is the Hamiltonian lift of  X i , i.e.,  H i ( ξ ( t ) ) = ξ ( t ) ( X i ( g ( t ) ) .
There are two kinds of extremal curves depending whether  λ = 0  (abnormal case) or  λ = 1  (normal case). In the abnormal case the Maximum principle results in the constraints
H i ( ξ ( t ) ) = 0 , i = 1 , , m ,
and beyond that gives no further information about the optimal control in question. In the normal case, however, the above maximum yields  u i ( t ) = H i ( ξ ( t ) ) , i = 1 , , m , where  ξ ( t )  is a solution curve of a single Hamiltonian vector field corresponding to the Hamiltonian
H ( ξ ) = 1 2 i = 1 m H i 2 ( ξ ) .
Each optimal solution  g ( t )  is either the projection of an abnormal or a normal extremal curve. If  g ( t )  is the projection of a normal extremal curve  ξ ( t )  then  ξ ( t )  is an integral curve of  H  and the control  u ( t )  that generates  g ( t )  is of the form  u i ( t ) = H i ( ξ ( t ) ) , i = 1 , , m .
We will not concern ourselves with the abnormal extremals. It is very likely that every optimal trajectory is the projection of a normal extremal curve, as in [52], in which case the abnormal extremals could be ignored. Instead we will turn to the normal Hamiltonian H and its Hamiltonian equations
d g d t = i = 1 n H i ( ( t ) ) X i ( g ( t ) ) , d d t = a d * d H ( ( t ) ) ( ( t ) ) .
Let us first consider the solutions of the associated Poisson equation
d d t = a d * d H ( ( t ) ) ( ( t ) )
and the structure of the coadjoint orbits.
Since  M ^  is a Euclidean vector space, its tangent space at the origin can be identified with  M ^ . Then the Lie algebra  G  will be identified with  g × M ^ , and its dual with  G * = g * M ^ * , where
g * = { G * : ( p ˙ ) = 0 , p ˙ M ^ } , M ^ * = { G * : ( g ) = 0 } .
It then follows that every  G *  can be written as  = 1 + 2  with  1 g *  and  2 M ^ * . Since  M ^  is an abelian algebra the projection  2  on  M ^ *  is constant on each coadjoint orbit of  G . The argument is straightforward:
A d g * ( ) ( X + p ˙ ) = ( A d g 1 ( X + p ˙ ) ) = ( A d g 1 ( X ) + p ˙ ) = 1 ( A d g 1 ( X ) ) + 2 ( p ˙ ) ,
for any  g = ( g , p ) G . It follows that the coadjoint orbits in  G  are of the form
{ A d g * ( 1 ) : g G } + 2 , for any = 1 + 2 .
This fact can be also verified directly from Equation (100): we have
d d t V = [ d H , V ] , for any V = X + p ˙ in G ,
where  d H = i = 1 m H i ( ) ( A i + A i ( o ) )  and  H i ( ) = 1 ( A i ) + 2 ( A i ( o ) ) . Therefore,
d 1 d t ( X ) + d 2 d t ( p ˙ ) = ( 1 + 2 ) ( [ d H , X + x ˙ ] ) = i = 1 m H i ( i ) [ A i , X ] .
from which follows that
d 1 d t ( X ) = i = 1 n H i ( i ) [ A i , X ] , X g , d 2 d t ( p ˙ ) = 0 .
Since  p ˙  is arbitrary  d 2 d t = 0 .
To uncover other constants of motion identify  G *  with  G  via the natural quadratic forms on each of the factors, and then recast the preceding equations on  G . More precisely, identify each  2  in  M ^ *  with a tangent vector  l = i = 1 m l i A i ( o )  via the formula  2 ( p ˙ ) = ( l , p ˙ ) , p ˙ M ^ . Similarly, identify  1 g *  with  L g  via the formula  1 ( X ) = L , X , X g . Then decompose  L g  into the sum  L = L p + L k L p p  and  L k k . Relative to the basis  A 1 , , A m  in  p L p = i = 1 m P i A i  where  P i = 1 ( A i ) = L , A i . It follows that
H i ( ξ ) = ( A i + A i ( o ) ) = 1 ( A i ) + 2 ( A i ( o ) ) = P i + l i ,
and
d 1 d t ( X ) = d L d t , X = L , [ i = 1 m ( l i + P i ) A i , X ] = [ L , i = 1 m ( l i + P i ) A i ] , X , ( d l d t , p ˙ ) = d 2 d t ( t ) ( p ˙ ) = 0
Since X and  p ˙  are arbitrary,
d L d t = [ i = 1 m ( l i + P i ) A i , L ] = [ A + L p , L ] , A = i = 1 m l i A i , d l d t = 0 .
Coupled with
d g d t = g ( t ) ( A + L p ) , d p d t = i = 1 n ( l i + P i ) A i ( o ) ,
Equation (101) constitute the Hamiltonian equations on  G × G  generated by the Hamiltonian  H = 1 2 i = 1 m H i 2 = 1 2 i = 1 m ( l i + P i ) 2 .
Each extremal curve projects onto a geodesic  g ( t ) = ( g ( t ) , p ( t ) ) , and each geodesic further projects onto the pair of curves  α ( t ) = τ g ( t ) ( o )  in M and  β ( t ) = p ( t )  in  M ^  that are rolled upon each other by  g ( t ) . Note that in this identification of the Lie algebras with their duals, coadjoint orbits  { A d g * ( 1 ) + 2 : g G }  are identified with the affine sets  { A d g ( L ) + l : g G } .
Recall now the Hamiltonian equations associated with the canonical affine-quadratic problem (Equation (33)):
d g d t = g ( t ) ( A + L k ( t ) ) , d L k d t = [ A , L p ] , d L p d t = [ L k , L p ] + s [ A , L k ] , s = 0 , 1 .
The propositions below reveal a remarkable fact that the Poisson equations of a canonical affine-quadratic Hamiltonian are subordinate to the Poisson equations associated with a rolling Hamiltonian. This connection identifies the drift term in the affine-quadratic system with a coadjoint invariant of the rolling Poisson system. To keep the systems apart we will use bold letters when referring to the variables in the rolling Hamiltonian in contrast to the variables in the affine-quadratic Hamiltonian which will remain the same.
Proposition 12. 
Let  ( g ( t ) , p ( t ) ) , L p ( t ) , L k ( t )  be an integral curve of the rolling Hamiltonian  H = 1 2 | | A + L p | | 2 , that is,
d g d t = g ( t ) ( A + L p ( t ) ) , d p d t = i = 1 m ( l i + P i ) A i ( o ) , d L k d t = [ A , L p ] , d L p d t = [ A + L p , L k ] , A = i = 1 m l i A i
Then
g ˜ ( t ) = g ( t ) h ( t ) , L p ( t ) = A d h 1 ( t ) ( L p ( t ) ) , L k = A d h 1 ( t ) ( L k ( t ) )
is an integral curve of the affine Hamiltonian  H = 1 2 L k , L k + A , L p , where  A = A d h 1 ( t ) ( A + L p ( t ) ) , and  h ( t )  is the solution of  d h d t = L k ( t ) h ( t )  with  h ( 0 ) = I .
Moreover, if  x ( t )  a solution of  d x d t = A + L p ( t )  then  g ˜ ( t ) = ( x ( t ) , h ( t ) )  in  p K  is the projection of an extremal curve
L k ( t ) = A d h 1 ( t ) L k ( t ) , L p ( t ) = A d h 1 ( t ) ( L p ( t ) ) A , A = A d h 1 ( t ) ( A + L p ( t ) )
associated with the shadow Hamiltonian  H = 1 2 L k , L k + A , L p .
The converse also holds according to the following proposition.
Proposition 13. 
Suppose that  ( g ˜ ( t ) , L p ( t ) , L k ( t ) )  is an extremal curve of the affine Hamiltonian  H = 1 2 L k , L k + A , L p . Let
g ( t ) = g ˜ ( t ) h 1 ( t ) , L p ( t ) = A d h ( t ) ( L p ( t ) ) , L k ( t ) = A d h ( t ) ( L k ( t ) ) , A = A d h ( t ) ( A L p ( t ) )
where  h ( t )  is a solution of  d h d t = h ( t ) ( L k ( t ) )  and let  p ( t )  be a solution of  d p d t = A ( o ) + L p ( t ) ( o ) .  Then  ( g ( t ) , p ( t ) )  together with
L p ( t ) = A d h ( t ) ( L p ( t ) ) , L k ( t ) = A d h ( t ) ( L k ( t ) ) , A = A d h ( t ) ( A L p ( t ) )
is an extremal curve of the rolling Hamiltonian  H = 1 2 A + L p , A + L p .
However, if  g ˜ ( t ) = ( x ( t ) , R ( t ) ) L p ( t ) + L k ( t ) )  is an extremal curve of the shadow Hamiltonian H, then  ( g ( t ) , p ( t ) ) , solutions of
d g d t = g ( t ) A d R ( t ) ( A ) ) , d p d t = d x d t ( o ) ,
together with
L p ( t ) = A d R ( t ) ( A + L p ( t ) ) , L k ( t ) = A d R ( t ) ( L k ( t ) ) , A = A d R L p
define an extremal curve of the rolling Hamiltonian  H = 1 2 A + L p , A + L p .
The proofs follow by straightforward calculations (also done in [9]).
Let us now come back to isospectral representations and Zimmerman’s method [17,52]. For that purpose let   X 0 ( t ) = A + L p ( t ) , X 1 ( t ) = L k ( t ) , X 2 ( t ) = A , X 3 = 0 . Then Poisson’s equations for the rolling problem can be written as
d X i d t = [ X 0 ( t ) , X i + 1 ( t ) ] , i = 0 , 1 , 2 .
These equations are invariant under a dilational change of variables  X i λ i 1 X i . It then follows that
L λ = i = 0 3 λ i X i = L p ( t ) + λ L k ( t ) + ( 1 λ 2 ) A
satisfies the equation
d L λ d t = [ M λ ( t ) , L λ ( t ) ] , M λ ( t ) = 1 λ ( A + L p ( t ) ) .
Therefore  L λ  is the spectral curve for  H . But then the Poisson system associated with the affine-quadratic Hamiltonian also admits an isospectral representation after the substitutions
A = A d h ( t ) ( A L p ) , L k = A d h ( t ) ( L k ) , L p = A d h ( t ) ( L p ) , d h d t = h ( t ) L k ( t ) .
For then  ( L p ( t ) , L k ( t ) )  are the extremal curves for the Poisson system associated with the affine-quadratic system (Proposition 13) and further satisfy
L λ = A d h ( t ) ( L p ) + λ A d h ( t ) ( L k ) + ( 1 λ 2 ( A d h ( t ) ( A L p ) = A d h ( t ) ( λ 2 L p + λ L k + ( 1 λ 2 ) A ) = A d h ( t ) L λ .
But then
A d h ( t ) [ 1 λ A , L λ ] = [ 1 λ ( A + L p ) , L λ ] = d L λ d t = d d t ( A d h ( t ) ( L λ ) = A d h ( t ) [ L λ , L k ] + A d h ( t ) d L λ d t
implies
d L λ d t = [ L k , L λ ] + [ 1 λ A , L λ ] = [ 1 λ A + L k , L λ ] .
To be consistent with my earlier publications, replace  λ  by  1 λ  to get
d L λ d t = [ M λ , L λ ] ,
where  M λ = L k λ A , and  L λ = L p λ L k + ( λ 2 1 ) A . Equation (107) agrees with the isospectral representation (34).
To get the spectral curve  L λ  for the shadow Hamiltonian, use relations  L k = A d h ( L k ) L p = A d h ( L p + A )  and  A = A d h L p  from Proposition 13. Then
L λ = L p + λ L k + ( 1 λ 2 ) A = A d h L λ , L λ = λ 2 L p + λ L k + A .
Then a calculation analogous to the one above gives  d L λ d t = [ 1 λ A + L k , L λ ] . After the rescaling  λ 1 λ  we get a modified Lax pair
d L λ d t = [ M λ , L λ ] , M λ = L k λ A , L λ = L p λ L k + λ 2 A .

5.2. Rolling Problem on Spaces of Constant Curvature

We will now introduce another optimal problem intertwined with the rolling problem. It consists of finding a continuously differentiable curve  p ( t )  in M in an interval  [ 0 , T ] , with its tangent vector  p ˙ ( t )  of unit length and its covariant derivative bounded and measurable in  [ 0 , T ]  that satisfies fixed tangential directions  p ˙ ( 0 ) = v 0 , v 0 T p ( 0 ) M  and  p ˙ ( T ) = v 1 , v 1 T p ( T ) M  along which the integral  1 2 0 T κ 2 ( s ) d s  minimal among all other curves that satisfy the same boundary conditions. Here  κ ( t ) = | | d D p ( t ) d t ( p ˙ ( t ) ) | | , where  d D p ( t ) d t  denotes the covariant derivative along  p ( t ) . The integral  1 2 0 T κ 2 ( s ) d s  is known as the elastic energy of the curve  p ( t )  [24]. Curves  p ( t )  defined on some interval  [ 0 , T ]  are called elastic if for each  t ( 0 , T )  there exits an interval  [ t 0 , t 1 ] [ 0 , T ]  over which the elastic energy of  p ( t )  is minimal relative to the boundary conditions  p ˙ ( t 0 )  and  p ˙ ( t 1 )  [7].
On semi-simple manifolds the curvature problem can be lifted to the unit tangent bundle of G, and it is this lifted version of the problem that will be of interest for this paper. In this formulation of the problem the tangent bundle of G is realized as the product  G × g  with  ( g , X ) G × g  identified with  g X T g G . Then each tangent vector  v T p M  is the projection of a manifold  V = { ( g h , A d h ( U ) ) , h K }  in  G × g  where  p = π ( g )  and  v = d g π ( g ) U , U p . The lifted curvature problem consists of finding a curve  ( g ( t ) , Λ ( t ) )  in  G × S p S p = { Λ p : Λ , Λ = 1 } , a solution of
d g d t = g ( t ) Λ ( t ) , d Λ d t = U ( t ) , U ( t ) , Λ ( t ) = 0 ,
that originates in the manifold  V 0 = { ( g 0 h , A d h 1 Λ 0 ) , h K , Λ 0 p }  at  t = 0  and terminates at the manifold  V 1 = { ( g 1 h , A d h 1 Λ 1 ) : h K , Λ 1 p }  at  t = T  for which the energy of transfer  1 2 0 T | | U ( s ) | | 2 d s  is minimal. If  p ( t ) = π ( g ( t ) ) = τ g ( t ) ( o )  is the projected curve, then  p ( t )  is the solution of
p ˙ ( t ) = d g ( t ) π ( g ( t ) ) Λ ( t ) = d o τ g ( t ) Λ ( t ) ( o ) .
that satisfies  | | p ˙ ( t ) | | = 1  and the boundary conditions
p ( 0 ) = π ( g ( 0 ) ) , p ˙ ( 0 ) = d g ( 0 ) π ( V 0 ) = d g ( 0 ) π ( g ( 0 ) Λ 0 ) , p ( T ) = π ( g ( T ) ) , p ˙ ( T ) = d g ( T ) π ( V 1 ) = d g ( T ) π ( g ( T ) Λ 1 ) .
It is a simple exercise to show that a curve  p ( t )  is elastic if and only if it is the projection of a solution of the lifted curvature problem on a fixed interval  [ 0 , T ] .
The Hamiltonian system for the curvature problem (Equation (109)) can also be obtained through the Maximum principle properly modified to account for the constraints, as outlined in ([7], Chapter 11). To go into these details would take us away from the central theme of the paper, so instead, we will just quote the relevant equations from [7] (pp. 354–355).
The curvature Hamiltonian H is given by  H = 1 2 | | X | | 2 + Λ , P  together with the associated Hamiltonian equations
d g d t = g Λ ( t ) , d P d t = [ Λ , Q ] , d Q d t = [ Λ , P ] , d Λ d t = X ( t ) , d X d t = P ( | | X | | 2 P , Λ ) Λ ,
subject to the transversality condition  Q ( t ) + [ Λ ( t ) , X ( t ) ] = 0 . The transversality condition can be incorporated into the above equations to yield an equivalent system
d g d t = g Λ ( t ) , d Λ d t = X ( t ) , d X d t = P ( | | X | | 2 P , Λ ) Λ , d P d t = [ Λ , [ Λ , X ] ] , d Q d t = [ Λ , P ] .
We will now confine our attention to spaces of constant curvature, with a particular interest on the connections between the rolling problems and the elastic curves reported in [52]. For those reasons let us return to the “spheres”  S ϵ n ( ρ ) = { x R n + 1 : ( x , x ) ϵ = ρ 2 , x 0 > 0 when ϵ = 1 }  and their rollings on the isometry groups groups  S O ϵ ϵ = ± 1  endowed with the quadratic form  A , B ϵ = 1 2 ϵ ρ 2 T r ( A B ) .  The rolling equations associated with the rollings of  S ϵ n ( ρ )  on the tangent plane  M ^ = T ρ e 0 S ϵ n ( ρ )  are given by
d g d t = g ( t ) ( u ( t ) ϵ e 0 ) ) , d p d t ( t ) = ρ u ( t ) .
In what follows we will make use of the following isospectral integrals of motion associated with the preceding rolling problem extracted from the functions  f 2 , λ = T r ( L λ 2 )  and  f 4 , λ = T r ( L λ 4 )
I 0 = 2 H = | | A + L p | | 2 , I 1 = | | L p | | 2 + ϵ | | L k | | 2 I 2 = | k | | | L k | | 2 | | L p | | 2 | | [ L p , L k ] | | 2 + ϵ 2 ( k | | L k | | 4 | | L k 2 | | 2 ) , I 4 = | k | | | L k | | 2 | | A + L p | | 2 | | [ A + L p , L k ] | | 2 .
These integrals of motion are rescaled variants of the integrals of motion in [52] after the metric is replaced by  A , B ρ = ρ A , B ϵ  (the metric in this paper is a scalar multiple of the metric used in [52]).
Recall that on spaces of constant Riemannian curvature the curvature k is defined by
[ V , [ V , X ] ] = k X ,
for any V and X in  p  that satisfy  | | V | | = 1  and  V , X = 0 . In particular k is equal to  ϵ ρ 2  on  S ϵ n ( ρ ) . Note that  1 2 T r ( A B ) = k A , B .
Proposition 14. 
Rolling geodesics that are the projections of the extremal curves on  H = 1 2  and  I 4 = 0  project on the elastic curves in  S ϵ n ( ρ ) . Conversely each elastic curve in  S ϵ n ( ρ )  is the projection of such an extremal curve.
Proof. 
Each elastic curve on  S ϵ n ( ρ )  is the projection of an extremal curve corresponding to the curvature problem (Equation (110)). On spaces of constant Riemannian curvature
[ Λ , [ Λ , X ] ] = k X .
Therefore, Equation (110) can be written as
d g d t = g Λ ( t ) , d Λ d t = X ( t ) , d X d t = P ( | | X | | 2 P , Λ ) Λ , d P d t = k X , d Q d t = [ Λ , P ] .
It follows that  k d Λ d t d P d t = 0 , and therefore,  k Λ P = k A  for some constant element A in  p . The transversality condition  Q + [ Λ , X ] = 0  can be recast as  0 = [ Λ , Q ] + [ Λ , [ Λ , X ] ] = [ Λ , Q ] k X . These observations can be incorporated in the preceding equations to get
d g d t = g ( t ) Λ ( t ) = g ( t ) 1 k ( k A + P ) , d P d t = k X = [ Λ , Q ] = 1 k [ k A + P , Q ] , d Q d t = [ Λ , P ] = 1 k [ k A + P , P ] = 1 k [ k A , P ] .
If we now identify  1 k P  with  L p , and  1 k Q  with  L k , then the preceding equations reduce to the rolling Hamiltonian system. Moreover,
Λ = A + 1 k P = A + L p , and L k = 1 k Q = 1 k [ A + L p , X ] .
Hence  | | A + L p | | = 1  so the first constraint is satisfied. To verify the second constraint note that  L k = 1 k [ A + L p , X ] , and therefore
| | L k | | 2 = 1 k 2 | | [ A + L p , X ] | | 2 = 1 k 2 [ [ A + L p , X ] , A + L p ] , X = 1 | k | | | X | | 2 ,
and  | | [ A + L p , L k ] | | 2 = 1 k 2 | | [ A + L p , [ A + L p , X ] ] | 2 = | | X | | 2 . Therefore,
I 4 = | | L k | | 2 | k | | | A + L p | | 2 | | [ A + L p , L k ] | | 2 = | | X | | 2 | | X | | 2 = 0 .
To prove the converse assume that  g ( t ) , p ( t ) , A , L k ( t ) , L p ( t )  is a rolling extremal curve on  I 4 = 0 . As a geodesic it satisfies  H = 1 2 , or  | | A + L p | | = 1 . We need to show that  L k ( t ) = [ A + L p ( t ) , X ( t ) ]  for some  X ( t ) p  such that  X ( t ) , A + L p ( t ) = 0 .
Let
Λ ( t ) = A + L p ( t ) , p Λ ( t ) = { X ( t ) p : X ( t ) , Λ ( t ) = 0 } , k Λ ( t ) = { Q ( t ) k : [ Q ( t ) , Λ ( t ) ] = 0 } , k Λ ( t ) = { Q k ; : Q , k Λ = 0 } .
Then  Λ ( t ) = λ ( t ) ϵ e 0 , ( λ ( t ) , e 0 ) ϵ = 0  for some vector  λ ( t ) R n + 1 . It then follows that  p Λ ( t ) = { u ( t ) ϵ e 0 : ( u ( t ) , e 0 ) ϵ = ( λ ( t ) , u ( t ) ) ϵ = 0 } , and  k Λ ( t ) = { λ ( t ) ϵ u ( t ) : ( u ( t ) , λ ( t ) ) ϵ = ( e 0 , u ( t ) ) ϵ = 0 } .
Hence, dim ( p Λ ( t ) ) =  dim( k Λ ) . The mapping  F ( t ) X = a d Λ ( t ) ( X ) , X p Λ ( t )  satisfies  F ( p Λ ( t ) ) k Λ ( t )  because  [ Λ , X ] , k Λ = 0  . On spaces of non-zero constant curvature, the kernel of this mapping is zero because  a d Λ ( t ) X = 0  implies that  0 = a d 2 Λ ( t ) ( X ( t ) ) = ϵ ρ 2 X ( t ) . Since  p Λ ( t )  and  k Λ ( t )  have the same dimension, F maps  p Λ ( t )  onto  k Λ ( t ) . So every curve  L ( t ) k Λ ( t )  is of the form  L ( t ) = [ Λ ( t ) , X ( t ) ]  for some  X ( t ) p  perpendicular to  Λ ( t ) .
It remains to show that  L k ( t )  belongs to  k Λ ( t )  when the rolling geodesic is on  I 4 = 0 , that is, when  | | L k | | 2 = 1 k | | [ Λ ( t ) , L k ( t ) | | 2 . Now assume that  L k ( t ) = U 1 ( t ) + U 2 ( t ) U 1 ( t ) k Λ ( t )  and  U 2 ( t ) k Λ ( t ) . It follows from above that  U 2 ( t ) = [ Λ ( t ) , X ( t ) ] , and therefore
| | U 2 ( t ) | | 2 = | | [ Λ ( t ) , X ( t ) ] , [ Λ ( t ) , X ( t ) ] | | 2 = | a d 2 Λ ( t ) ( X ) , X ( t ) | = | k | | | X | | 2
Hence,
1 | k | | | [ Λ ( t ) , U 1 ( t ) + U 2 ( t ) ] | | 2 = 1 | k | | | [ Λ ( t ) , U 2 ( t ) ] | | 2 = | k | | | X | | 2 = | | L k | | 2 .
But  | | L k ( t ) | | 2 = | | U 1 | | 2 + | | U 2 ( t ) | | 2 = | | U 1 ( t ) | | 2 + | k | | | X | | 2 , and therefore  U 1 ( t ) = 0 . □
The following proposition characterizes elastic curves [7].
Proposition 15. 
Let  κ ( t )  and  τ ( t )  denote the geodesic curvature and the torsion of the projection curve  p ( t )  associated with an extremal curve of the curvature problem. Then  ξ ( t ) = κ 2 ( t )  is the solution of the following equation
( d ξ d t ) 2 = ξ 3 + 4 ( H ϵ ) ξ 2 + 4 ( I 1 H 2 ) ξ 4 I 2 ,
and  ( κ 2 ( t ) τ ( t ) ) 2 = k I 2 . All other curvatures in the Serret-Frenet frame along  p ( t )  are zero.
I believe that the proof given below is more to the point than similar proofs given elsewhere [7,52].
Proof. 
We leave it to the reader to verify that  k | | L k | | 4 | | L k 2 | | 2 = 0  when  L k = 1 k [ A + L p , X ] . Let  P ( t ) = k L p ( t )  and let  Q ( t ) = k L k ( t ) . Then
I 2 k 2 = k 2 ( | k | | | L k | | 2 | | L p | | 2 | | [ L p , L k ] | | 2 ) = | | P | | 2 | | X | | 2 | | | | [ L p , Q ] | | 2 = | | P | | 2 | | X | | 2 | | [ L p , [ A + L p , X ] ] | | 2 = | | P | | 2 | | X | | 2 A + L p , L p 2 | | X | | 2 P , X 2 .
Also
I 1 k 2 = k 2 ( | | L p | | 2 + ϵ | | L k | | 2 ) = | | P | | 2 + ϵ k 2 | | Q | | 2 = | | P | | 2 + k | | X | | 2 .
Since  κ 2 ( t ) = | | X | | 2 d ξ d t = 2 X , X ˙ = 2 X , P . Therefore,
( d ξ d t ) 2 = 4 P , X 2 = 4 ( | | P | | 2 | | X | | 2 A + L p , P 2 | | X | | 2 ) 4 k 2 I 2 = 4 ( ( I 1 k 2 k | | X | | 2 ) | | X | | 2 ( H 1 2 | | X | | 2 ) 2 | | X | | 2 ) 4 k 2 I 2 = 4 ( I 1 k 2 k ξ ) ξ 4 ( H 1 2 ξ ) 2 ξ 4 k 2 I 2 = ξ 3 + 4 ( H k ) ξ 2 + ( I 1 k 2 H 2 ) ξ 4 k 2 I 2 .
As to the second part, let  T = A + L p ( t ) . Since  | | A + L p ( t ) | | = 1 T ( t )  is a unit vector that projects onto the tangent vector  p ˙ ( t ) . Then
d T d t = [ A + L p ( t ) , L k ( t ) ] = [ A + L p ( t ) , [ A + L p ( t ) , 1 k X ( t ) ] ] = X ( t ) .
Therefore  d T d t = κ ( t ) N ( t )  where  N ( t ) = 1 | | X ( t ) | | X ( t )  is a unit vector in  p  that projects onto the unit normal  n ( t )  along  p ( t ) . Continuing,
d N d t = 1 | | X ( t ) | | ( P ( | | X | | 2 Λ , P ( A + L p ) ) 1 | | X | | 2 X , X ˙ X ) = | | X | | ( A + L p ) + 1 | | X | | ( P + A + L p , P ( A + L p ) ) + 1 | | X | | 2 P , X X ) = κ ( t ) T ( t ) + Y ( t ) ,
where
Y ( t ) = 1 | | X | | ( P ( t ) + T ( t ) , P ( t ) T ( t ) ) + 1 | | X | | 2 P ( t ) , X ( t ) X .
Since  Y ( t )  is is orthogonal to  A + L p  and X, it is in the direction of the binormal vector  B ( t ) . So if we define  τ ( t ) = | | Y ( t ) | |  and  B ( t ) = 1 | | Y | | Y  then  d N d t = κ ( t ) T ( t ) + τ B ( t )  and  B ( t )  projects onto the binormal vector  b ( t )  along  p ( t ) . Hence,
| | X | | 2 τ 2 = | | P | | 2 A + L p , P 2 1 | | X | | 2 P , X 2 ,
or
| ( κ 2 τ ) 2 = | | X | | 4 τ 2 = | | P | | 2 | | X | | 2 A + L p , P 2 P , X 2 = k 2 I 2 ,
Evidently  d B d t  is in the linear span of  T ( t ) , N ( t ) , B ( t ) , hence the Serret-Frenet frame along  p ( t )  terminates. □
Corollary 3. 
Elastic curves in  M ϵ = S ϵ n ( ρ )  are rolled on the elastic curves in the tangent space  M ^ = T e M .
Proof. 
Since the geodesic curvature is preserved under rolling, the elastic curves in  S ϵ n ( ρ )  are rolled on the elastic curves in  M ^  relative to the Euclidean metric inherited from the metric on  p . So the statement follows from the rolling definition. □
This remarkable relation between the elastic curves and the rolling geodesics breaks down on spaces of non-constant Riemannian curvature, as it becomes evident when one compares Equation (110) for the curvature problem to the Equation (101) for the rolling problem. It is interesting to note that the solutions of either of these two Equations (101) and (110) are not known beyond the spaces of constant curvature. While the curvature equation seems particularly challenging beyond the spaces of constant curvature, the rolling geodesic equations remain integrable on all semi-simple spaces and should be “solvable” according to the general theory of integrable systems.
Apart from the above remarks, there is another spectacular property of elastic curves that makes them special: elastic curves appear as soliton solutions in the non-linear Schroedinger equation [53]. More generally it was shown in [53] that the space of periodic horizontal curves of fixed length L in the isometry group G over a three dimensional space of constant curvature can be given a structure of an infinite dimensional Poisson manifold relative to which some famous equations of mathematical physics appear as Poisson equations associated with geometric invariants of curves on the base space. In particular, Heisenberg’s magnetic equation and Schroedinger’s non-linear equation appear as Poisson equation associated with  f 0 ( g ( s ) ) = 1 2 0 L | | d Λ d s ( s ) | | 2 d s  where  Λ ( s ) = g 1 ( s ) d g d s ( s ) | | Λ ( s ) | | = 1 . Since this function can be also expressed as  1 2 0 L κ 2 ( s ) d s  where  κ ( s )  is the geodesic curvature of the projected curve in the underlying symmetric space, elastic curves appear naturally in this setting (see also [54,55,56,57] for related results). This leap to infinite dimensional Hamiltonians and related hierarchies of commuting Hamiltonians further illustrates the relevance of Lie algebraic methods in the theory of integrable systems.

Funding

The research received no external funding.

Data Availability Statement

No additional data.

Conflicts of Interest

There is no conflict of interest.

References

  1. Poisson, S.D. Sur les inégalités séculaires des moyens mouvemens des planétes. J. L’École Polytech. 1809, 8, 15. [Google Scholar]
  2. Poincaré, H. Les Méthodes Nouvelles de la Mécanique Célecte; Tome I; Gauther-Villars: Paris, France, 1892. [Google Scholar]
  3. Jacobi, C.G.J. Vorlersungen Über Dynamic; Druck und Verlag: Berlin, Germany, 1884. [Google Scholar]
  4. Arnold, V.I. Mathematical Methods of Classical Mechanics; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
  5. Carathéodory, C. Calculus of Variations and Partial Differential Equations of the First Order; Second (revised) English Translation; Chelsea Publishing Co.: White River Junction, VT, USA, 1982. [Google Scholar]
  6. Liouville, J. Note sur l’intégration des équations différentielles de la Dynamique, présenté au Bureau des Longitudes le 29 juin 1853. JMPA 1855, Tome XX, 137–138. [Google Scholar]
  7. Jurdjevic, V. Optimal Control and Geometry: Integrable Systems; Cambridge Studies in Advanced Mathematics 154; Cambridge University Press: Cambridge, UK, 2016. [Google Scholar]
  8. Jurdjevic, V. Integrable Hamiltonian Systems on Lie groups: Kowalewski type. Ann. Math. 1999, 150, 605–644. [Google Scholar] [CrossRef]
  9. Jurdjevic, V. Rolling geodesics, mechanical systems and elastic curves. Mathematics 2022, 10, 4827. [Google Scholar] [CrossRef]
  10. Reyman, A.G.; Semenov Tian-Shansky, M.A. Group Theoretic Methods in the Theory of Finite Dimensional Inrtegrable Systems; Dynamical Systems VII, Chapter 2, Encyclopaedia of Mathematical Sciences; Springer: Berlin/Heidelberg, Germany, 1994; Volume 16. [Google Scholar]
  11. Kirillov, A.A. Elements of the Theory of Representations; Springer: Berlin/Heidelberg, Germany, 1976. [Google Scholar]
  12. Gasparim, E.; Grama, L.; San Martin, L. Adjoint orbits of semi-simple Lie groups and Lagrangian submanifolds. Proc. Edinb. Math. Soc. 2017, 60, 361–385. [Google Scholar] [CrossRef]
  13. Báez, J.; San Martin, L.A.B. Deformations of adjoint orbits for semisimple Lie algebras and Lagrangian submanifolds. Differ. Geom. Its Appl. 2021, 75, 101719. [Google Scholar] [CrossRef]
  14. Jurdjevic, V. Affine-Quadratic problems on Lie Groups: Tops and Integrable Systems. J. Lie Theory 2020, 30, 425–444. [Google Scholar]
  15. Helgason, S. Differential Geometry, Lie Groups and Symmetric Spaces; Academic Press: New York, NY, USA, 1978. [Google Scholar]
  16. Trélat, E. Some properties of the value function and its level sets for affine control systems with quadratic cost. J. Dyn. Control Syst. 2000, 6, 511–541. [Google Scholar] [CrossRef]
  17. Zimmerman, J. The Rolling Sphere Problem. Ph.D. Thesis, University of Toronto, Toronto, ON, Canada, 2002; 94p. [Google Scholar]
  18. Zimmerman, J.A. Optimal control of the sphere Sn rolling on En. Math. Control Signals Syst. 2005, 17, 14–37. [Google Scholar] [CrossRef]
  19. Bolsinov, A. A completeness criterion for a family of functions in involution obtained by the shift method. Soviet Math. Dokl. 1989, 38, 161–165. [Google Scholar]
  20. Fomenko, A.T.; Mischenko, A.S. Euler equation on finite-dimensional Lie groups. Math. USSR Izv. 1978, 12, 371–389. [Google Scholar]
  21. Fomenko, A.T.; Trofimov, V.V. Integrability in the sense of Liouville of Hamiltonian systems on Lie algebras. Uspekhi Mat. Nauk 1984, 2, 3–56. (In Russian) [Google Scholar]
  22. Manakov, S.V. Note on the integration of Euler’s equations of the dynamics of an n dimensional rigid body. Funct. Anal. Appl. 1976, 10, 328–329. [Google Scholar] [CrossRef]
  23. Love, A.E. A Treatise on the Mathematical Theory of Elasticity, 4th ed.; Dover: New York, NY, USA, 1927. [Google Scholar]
  24. Jurdjevic, V. Integrable Hamiltonian Systems on Complex Lie Groups. Mem. Am. Math. Soc. 2005, 178, 838. [Google Scholar] [CrossRef]
  25. Bogoyavlenski, O. New Integrable Problem of Classical Mechanics. Commun. Math. Phys. 1984, 94, 255–269. [Google Scholar] [CrossRef]
  26. Kowalewski, S. Sur le problème de la rotation d’un corps solide autor d’un point fixé. Acta Math. 1889, 12, 177–232. [Google Scholar] [CrossRef]
  27. Jurdjevic, V. Kowalewski top and complex Lie algebras. Anal. Math. Phys. 2021, 11, 38. [Google Scholar] [CrossRef]
  28. Komarov, I.V.; Kuznetsov, V.B. Kowalewski top on the Lie algebras o(4), e(3) and o(3, 1). J. Phys. A 1990, 23, 841–846. [Google Scholar] [CrossRef]
  29. Komarov, I.V. Kowalewski top for the hydrogen atom. Theor. Math. Phys. 1981, 47, 67–72. [Google Scholar] [CrossRef]
  30. Dragović, V.; Kukić, K. Systems of Kowalevski Type and Discriminantly Separable Polynomials. Regul. Chaotic Dyn. 2014, 19, 162–184. [Google Scholar]
  31. Sokolov, V.V. A New Integrable Case for the Kirchhoff Equation. Theor. Math. Phys. 2001, 129, 1335–1340. [Google Scholar] [CrossRef]
  32. Ivanescu, C.; Savu, A. The Kowalewski top as a reduction of a Hamiltonian system on Sp(4,ℝ). Proc. Am. Math. Soc. 2003, 131, 607–618. [Google Scholar] [CrossRef]
  33. Haine, L.; Horozov, E. A Lax pair for Kowalewski top. Phys. D 1987, 29, 173–180. [Google Scholar] [CrossRef]
  34. Horozov, E.; van Moerbeke, P. The full geometry of Kowalewski’s top and (1,2)-abelian surfaces. Commun. Pure Appl. Math. 1989, 42, 357–407. [Google Scholar] [CrossRef]
  35. Bobenko, A.I.; Reyman, A.G.; Semenov-Tian Shansky, M.A. The Kowalewski top 99 years later; a Lax pair, generalizations and explicit solutions. Commun. Math. Phys. 1989, 122, 321–354. [Google Scholar] [CrossRef]
  36. Jurdjevic, V. Affine-quadratic problems on Lie groups. Math. Control Rel. Fields 2013, 3, 347–374. [Google Scholar] [CrossRef]
  37. Neumann, C. De problemate quodam mechanico, quod ad primam integralium ultraellipticorum classem revocatum. J. Reine. Angew. Math. 1859, 56, 46–63. [Google Scholar]
  38. Ratiu, T. The C. Newmann problem as a completely integrable system on a coadjoint orbit. Trans. Am. Mat. Soc. 1981, 264, 321–329. [Google Scholar] [CrossRef] [Green Version]
  39. Moser, J. Integrable Hamiltonian Systems and Spectral Theory; Lezioni Fermiane, Academia Nazionale dei Lincei, Scuola Normale Superiore: Pisa, Italy, 1981. [Google Scholar]
  40. Perelomov, A.M. Integrable Systems of Classical Mechanics and Lie Algebras; Birkhauser Verlag: Basel, Switzerland, 1990; Volume 1. [Google Scholar]
  41. Knörrer, H. Geodesics on quadrics and a mechanical problem of C. Newmann. J. Riene Angew. Math. 1982, 334, 69–78. [Google Scholar]
  42. Fock, V.A. The hydrogen atom and non-Euclidean geometry. Izv. Akad. Nauk SSSR Ser. Fizika 1935, 8. [Google Scholar]
  43. Moser, J. Regularization of Kepler’s problem and the averaging method on a manifold. Commun. Pure Appl. Math. 1970, 23, 609–623. [Google Scholar] [CrossRef]
  44. Osipov, Y. The Kepler problem and geodesic flows in spaces of constant curvature. Celest. Mech. 1977, 16, 191–208. [Google Scholar] [CrossRef]
  45. Guillemin, V.; Sternberg, S. Variations on a Theme by Kepler; American Mathematical Society: Providence, RI, USA, 1990; Volume 42. [Google Scholar]
  46. O’Neill, B. Semi-Riemannian Geometry; Academic Press: Cambridge, MA, USA; Elsevier: Amsterdam, The Netherlands, 1983. [Google Scholar]
  47. Jurdjevic, V.; Markina, I.; Silva Leite, F. Symmetric spaces rolling on flat spaces. J. Geom. Anal. 2023, 33, 94. [Google Scholar] [CrossRef]
  48. Ziller, W. Lie Groups. Representation Theory and Symmetric Spaces; University of Pennsylvania: Philadelphia, PA, USA, 2010. [Google Scholar]
  49. Bryant, R.; Hsu, L. Rigidity of integral curves of rank 2 distributions. Invent. Math. 1993, 114, 435–461. [Google Scholar] [CrossRef]
  50. Agrachev, A.; Sachkov, Y. Control Theory from the Geometric Viewpoint; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  51. Chitour, Y.; Godoy-Molina, M.; Kokkonen, P. The Rolling Problem: Overview and Challenges. Geometric Control Theory and Sub-Riemannian Geometry; Springer INdAM Ser. 5; Springer: Berlin/Heidelberg, Germany, 2014; pp. 103–122. [Google Scholar]
  52. Jurdjevic, V.; Zimmerman, J. Rolling sphere problems on spaces of constant curvature. Math. Proc. Camb. Phil. Soc. 2008, 144, 729–747. [Google Scholar] [CrossRef]
  53. Jurdjevic, V. The symplectic structure of curves in three dimensional spaces of constant curvature and the equations of mathematical physics. Ann. I. H. Poincaré 2009, 26, 1483–1515. [Google Scholar] [CrossRef] [Green Version]
  54. Langer, J.; Perline, R. Poisson Geometry of the Filament Equation. J. Nonlinear Sci. 1978, 1, 71–93. [Google Scholar] [CrossRef]
  55. Chabat, C.; Zakharov, V. Exact theory of two dimensional self-focusing and one dimensional self-modulation of waves in non-linear media. Sov. Phys. JETP 1972, 34, 62–69. [Google Scholar]
  56. Hasimoto, H. A soliton on a vortex element. J. Fluid Mech. 1972, 51, 477–485. [Google Scholar] [CrossRef]
  57. Millson, J.; Zombro, B.A. A Kähler structure on the moduli spaaces of isometric maps of a circle into Euclidean spaces. Invent. Math. 1995, 123, 35–59. [Google Scholar] [CrossRef]
Table 1. Lie brackets for  s = 0 , 1 .
Table 1. Lie brackets for  s = 0 , 1 .
[ , ]   A 1   A 2   A 3   B 1   B 2   B 3
  A 1 0   A 3   A 2 0   B 3   B 2
  A 2   A 3 0   A 1   B 3 0   B 1
  A 3   A 2   A 1 0   B 2   B 1 0
  B 1 0   B 3   B 2 0   s A 3   s A 2
  B 2   B 3 0   B 1   s A 3 0   s A 1
  B 3   B 2   B 1 0   s A 2   s A 1 0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jurdjevic, V. Integrable Systems: In the Footprints of the Greats. Mathematics 2023, 11, 1063. https://doi.org/10.3390/math11041063

AMA Style

Jurdjevic V. Integrable Systems: In the Footprints of the Greats. Mathematics. 2023; 11(4):1063. https://doi.org/10.3390/math11041063

Chicago/Turabian Style

Jurdjevic, Velimir. 2023. "Integrable Systems: In the Footprints of the Greats" Mathematics 11, no. 4: 1063. https://doi.org/10.3390/math11041063

APA Style

Jurdjevic, V. (2023). Integrable Systems: In the Footprints of the Greats. Mathematics, 11(4), 1063. https://doi.org/10.3390/math11041063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop