## Keplerian orbits through the Conley-Zehnder index

Henry Kavle, Daniel Offin, Alessandro Portaluri ∗

January 10, 2023

## Abstract

It was discovered by Gordon [Gor77] that Keplerian ellipses in the plane are minimizers of the Lagrangian action and spectrally stable as periodic points of the associated Hamiltonian flow. The aim of this note is to give a direct proof of these results already proved by authors in [HS10, HLS14] through a self-contained and explicit computation of the Conley-Zehnder index through crossing forms in the Lagrangian setting.

The techniques developed in this paper can be used to investigate the higher dimensional case of Keplerian ellipses, where the classical variational proof no longer applies.

## In memory of our friend Florin Diacu

AMS Subject Classification: 70F05, 53D12, 70F15.

Keywords: Two body problem, Conley-Zehnder index, Linear and Spectral Stability.

## Introduction

In the remarkable paper [Gor77], Gordon was able to apply the Tonelli direct method in Calculus of Variations for the Lagrangian action functional of the gravitational central force problem, by proving that the infimum of the Lagrangian action functional on the space of loops in the plane avoiding the origin and having non-vanishing winding number about the origin is realized by the Keplerian orbits (ellipses), including the limiting case of the elliptic collision-ejection orbit which passes through the origin. Excluding the latter, these all have winding numbers 1 (for the direct orbits) or -1 (for the retrograde orbits).

Starting from the aforementioned seminal paper of Gordon dozens of papers using methods from Calculus of Variations for weakly singular Lagrangian problems were published in the last decades.

The main difficulties addressed by author in the aforementioned paper are due to the lack of compactness. The non-compactness arises because of the unboundedness of the configuration space as well as the presence of the singularity. The non-compactness due to the unboundedness of the configuration space can be cured by properly defining the class of loops for minimizing the action. This class of loops was defined by introducing a sort of tied condition in terms of the winding number. In more abstract terms, the infimum of the action functional in the two of the infinitely many (labelled by the winding number) path connected components of the loop space of the plane with the origin removed having winding number ± 1 , is attained precisely on the Keplerian ellipses (other than the elliptic ejection-collision solutions). We observe that all elliptic orbits with fixed period T have the same Lagrangian action.

In [HS10], authors were interested in studying the relation between the Morse index and the stability for the elliptic Lagrangian solutions of the three body problem. In particular, in the first half of this paper, starting from Gordon's theorem and by using an index theory of periodic solutions of Hamiltonian system and in particular the Bott-type iteration formula, the authors

∗ Authors are partially supported Prin 2015 'Variational methods, with applications to problems in mathematical physics and geometry' No. 2015KB9WPT \_ 001 .

were able to get a stability criterion for the elliptic Keplerian orbits as well as compute the Morse index of all of their iterations.

It is interesting to observe that through the approach developed in this paper, it is possible to explicitly compute the index properties of closed Keplerian orbits for surfaces of constant Gaussian curvature among others. We conclude by observing that, Gordon's theorem breaks down if the dimension of the configuration space is bigger than 2, due to the fact that the loop space of the Euclidean n -dimensional space (for n glyph[greaterorequalslant] 3 ) with the origin removed is path connected, contrary to what happens with our approach.

This paper is directed towards a twofold aim: to recover Gordon results through the use of the Conley-Zehnder index, and to recover the results on the Keplerian ellipses given by authors in [HS10] in a maybe more direct way through crossing forms in the Lagrangian setting.

The paper is organized as follows:

## Contents

| 1   | Variational and Geometrical framework         | Variational and Geometrical framework            |   3 |
|-----|-----------------------------------------------|--------------------------------------------------|-----|
|     | 1.1                                           | Variational setting for the Kepler problem . . . |   3 |
|     | 1.2                                           | Geometrical properties of Keplerian orbits . . . |   4 |
| 2   | Maslov index and Conley-Zehnder index         | Maslov index and Conley-Zehnder index            |   6 |
|     | 2.1                                           | On the Conley-Zehnder index . . . . . . . . . .  |   6 |
|     | 2.2                                           | The Maslov index in Λ(2) . . . . . . . . . . . . |   9 |
| 3   | Some explicit computations for paths in Sp(2) | Some explicit computations for paths in Sp(2)    |  12 |
| 4   | Indices and stability of Keplerian orbits     | Indices and stability of Keplerian orbits        |  17 |
|     | 4.1                                           | The Conley-Zehnder index in polar coordinates    |  18 |

## Acknowledgements

The third name author wishes to thank all faculties and staff at the Queen's University (Kingston) for providing excellent working conditions during his stay and especially his wife, Annalisa, that has been extremely supportive of him throughout this entire period and has made countless sacrifices to help him getting to this point.

## Notation

For the sake of the reader, we introduce some notation that we shall use throughout the paper.

- -We denote by R + = (0 , ∞ ) the positive real numbers. The symbol 〈· , ·〉 k (or just 〈· , ·〉 if no confusion can arise) denotes the Euclidean product in R k and |·| denotes the (Euclidean) norm. I k denotes the identity matrix in R k ; as shorthand, we use just the symbol I . By U ⊂ C we denote the unit circle (centered at 0 ) of the complex plane.
- -X = R 2 \{ } 0 denotes the configuration space, T X its tangent bundle, or state space, and T ∗ X its cotangent bundle, or phase space. Λ X denotes the Hilbert manifold of loops of length T on X having Sobolev regularity H 1 . T = R /T Z denotes the unit circle of length T .
- -P T (2 n ) is the set of continuous symplectic maps ψ : [0 , T ] → Sp(2 n ) such that ψ (0) = I , as defined in Equation (2.1).
- -B sym ( V ) is the set of symmetric bilinear forms on the vector space V . For B ∈ B sym ( V ) , we denote by σ B ( ) its spectrum and by n -( B ) , n + ( B ) and n 0 ( B ) , its index (total number of negative eigenvalues), its coindex (total number of positive eigenalues) and its nullity (dimension of the kernel), respectively.
- The signature of B is defined by sign ( B ) = n + ( B ) -n -( B ) .

## 1 Variational and Geometrical framework

The aim of this section is to briefly describe the problem in its appropriate variational setting and to discuss some basic properties of the bounded non-colliding motions that we shall use for proving our main results.

## 1.1 Variational setting for the Kepler problem

We consider the gravitational force interaction between two point particles (or bodies) in the Euclidean plane having masses m ,m 1 2 ∈ R + . It is well-known that the motion of the two bodies interacting through the gravitational potential is mathematically equivalent to the motion of a single body with a reduced mass equal to

$$\mu = \frac { m _ { 1 } m _ { 2 } } { m _ { 1 } + m _ { 2 } }$$

that is acted on by an (external) attracting gravitational central force (pointing toward the origin). So, we let X = R 2 \ { 0 } be the configuration space . The elements of the state space (namely the tangent bundle T X ∼ = X × R 2 ) are denoted by ( q, v ) where q ∈ X and v ∈ T q X . We denote by U : X → R the Keplerian potential (function) , which is defined as follows

$$U ( q ) \, = \frac { m } { | q | }$$

for m = Gm m 1 2 and G ∈ R + denoting the gravitational constant . We now consider the Lagrangian L : T X → R defined by

$$L ( q, v ) = K ( v ) + U ( q ) \quad \text{ where } K ( v ) \, = \frac { 1 } { 2 } \mu \, | v | ^ { 2 }$$

where µ is given in Equation (1.1) and U is given in Equation (1.2) and we observe that the potential U is a positively homogeneous function of degree -1 .

We denote by T ∗ X the phase space (i.e. the cotangent bundle of X ). Elements of T ∗ X are denoted by z = ( p, q ) where q ∈ X and p ∈ T ∗ q X . Since the Lagrangian function L defined in Equation (1.3) is fiberwise C 2 -convex (being quadratic in the velocity v ), then the Legendre transformation

↦

$$\mathcal { L } _ { L } \, \colon T \mathbb { X } \to T ^ { * } \mathbb { X }, \quad ( q, v ) \mapsto \left ( D _ { v } L ( q, v ), q \right ),$$

is a smooth (local) diffeomorphism. The Fenchel transform of L is the autonomous smooth Hamiltonian on T ∗ X defined by

$$H ( p, q ) \, = \, \max _ { v \in T _ { q } \mathbb { X } } \left ( p [ v ] - L ( q, v ) \right ) = p \left [ v ( p, q ) \right ] - L ( q, v ( p, q ) ),$$

where ( q, v ( p, q ) ) = L -1 L ( p, q ) .

Given T &gt; 0 , we denote by T = R /T Z the circle of length T and we denote by Λ X = H 1 ( T X , ) the Hilbert manifold of T -periodic loops in X having Sobolev regularity H 1 with respect to the scalar product

$$\langle \gamma _ { 1 }, \gamma _ { 2 } \rangle _ { H ^ { 1 } } = \int _ { 0 } ^ { T } \left [ \langle \gamma _ { 1 } ^ { \prime } ( t ), \gamma _ { 2 } ^ { \prime } ( t ) \rangle + \langle \gamma _ { 1 } ( t ), \gamma _ { 2 } ( t ) \rangle \right ] d t.$$

We consider the Lagrangian action functional A : Λ X → R given by

$$\mathbb { A } ( \gamma ) \, = \int _ { 0 } ^ { T } L \left ( \gamma ( t ), \dot { \gamma } ( t ) \right ) d t.$$

It is well known that the Lagrangian action A is of class C 2 (and further, is actually smooth). By a straightforward calculation of the first variation and up to some standard regularity arguments, it follows that the critical points of A are C 2 solutions of the Euler-Lagrange equation

$$\mu \ddot { \gamma } ( t ) & = - \frac { m } { | \gamma ( t ) | ^ { 3 } } \gamma ( t ), \quad t \in ( 0, T )$$

Given a classical solution γ of Equation (1.5), the second variation of A is represent by the index form given by

$$d ^ { 2 } \mathbb { A } ( \gamma ) [ \xi, \eta ] = \int _ { 0 } ^ { T } \left [ \langle P \dot { \xi } ( t ), \dot { \eta } ( t ) \rangle + \langle R ( t ) \xi ( t ), \eta ( t ) \rangle \right ] d t \quad \forall \xi, \eta \in \mathbb { A } \mathbb { X }$$

where P = µI and R t ( ) = D qq L γ t , γ ( ( ) ˙ ( t ) ) . In particular, it is a continuous symmetric bilinear form on the Hilbert space W consisting of the T -periodic H 1 sections ξ of γ ∗ ( T X ) . Since the Lagrangian L is exactly quadratic with respect to v , it follows that the index form is a compact perturbation of the form

↦

$$( \xi, \eta ) \mapsto \int _ { 0 } ^ { T } \left [ \langle P \dot { \xi } ( t ), \dot { \eta } ( t ) \rangle + \langle \xi ( t ), \eta ( t ) \rangle \right ] d t$$

which is coercive on W and hence it is an essentially positive Fredholm quadratic form. In particular, its Morse index (i.e. the maximal dimension of the subspace of W such that the restriction of the index form is negative definite), is finite. Thus, we introduce the following definition:

Definition 1.1. Let γ be a critical point of A . We define the Morse index of γ as the Morse index of D 2 A ( γ ) , the second Frechét derivative of A at γ .

## 1.2 Geometrical properties of Keplerian orbits

↦

By changing to polar coordinates ( r, ϑ ) in the configuration space X , we observe that the Lagrangian L defined in Equation (1.3) along the smooth curve t → γ t ( ) = ( r t , ϑ ( ) ( ) t ) reduces to

$$L ( r, \vartheta, \dot { r }, \dot { \vartheta } ) & = \frac { 1 } { 2 } \mu ( \dot { r } ^ { 2 } + r ^ { 2 } \dot { \vartheta } ^ { 2 } ) + U ( r ).$$

Thus the Euler-Lagrange given in Equation (1.5) fits into the following

$$\begin{cases} \mu \ddot { r } - \mu r \dot { \vartheta } ^ { 2 } + \frac { m } { r ^ { 2 } } = 0 \quad \text{on } \ [ 0, T ] \\ \frac { d } { d t } ( \mu r ^ { 2 } \dot { \vartheta } ) = 0. \end{cases}$$

We refer to the first (respectively second) differential equation in Equation (1.8) as the radial Kepler (respectively transversal Kepler ) equation . By a direct integration in the second equation, we directly get

$$\dot { \vartheta } = \frac { k } { \mu r ^ { 2 } }, \quad k \in \mathbb { R } ^ { + }$$

where k denotes the modulus of the angular momentum (which is in fact a conserved quantity). A second constant of motion is given by the energy h . Energy is constant because there are no external forces acting on the reduced body; hence the Lagrangian is time independent. By substituting ˙ ϑ given in Equation (1.9) in the radial Kepler equation, we get

$$\mu \ddot { r } - \frac { k ^ { 2 } } { \mu r ^ { 3 } } + \frac { m } { r ^ { 2 } } = 0 \quad \text{ on } [ 0, T ].$$

So, we define the effective potential energy as

$$U _ { k } ( r ) \, = \frac { k ^ { 2 } } { 2 \mu r ^ { 2 } } - \frac { m } { r }$$

and we observe that Equation (1.10) is nothing but the equation of motion of a particle moving on the line in the force field generated by the potential function U k . By the energy conservation law we get that the energy level h is determined by

$$\frac { 1 } { 2 } \mu \dot { r } ^ { 2 } + \frac { k ^ { 2 } } { 2 \mu r ^ { 2 } } - \frac { m } { r } = h.$$

Let us now introduce the new time τ (usually called eccentric anomaly ) and defined by

$$d \tau = \frac { 1 } { r ( t ) } d t.$$

It is worth noticing that the scaling function is a priori unknown (as it is dependent on the unknown function r ). Thus, we get

$$\frac { d } { d t } = \frac { 1 } { r } \frac { d } { d \tau }, \quad \frac { d ^ { 2 } } { d t ^ { 2 } } = \frac { 1 } { r ^ { 2 } } \frac { d ^ { 2 } } { d \tau ^ { 2 } } - \frac { r ^ { \prime } } { r ^ { 3 } } \frac { d } { d \tau }$$

where we denote d dτ by a prime · ′ . The energy level h given in Equation (1.12) in the new time variable can be written as follows

$$\mu ^ { 2 } r ^ { \prime 2 } + k ^ { 2 } - 2 m \mu r = 2 \mu r ^ { 2 } h$$

and thus the radial Kepler equation given in Equation (1.10), reduces to the following linear second order (non-homogeneous ordinary differential) equation

$$\mu r ^ { \prime \prime } - 2 r h - m = 0 \quad \text{on } \, ( 0, \mathfrak { T } ).$$

We observe also that, in the new time variable τ , the (prime) period of the solution is given by

$$\mathfrak { T } & = \frac { 2 \pi } { \omega } \quad \text{where} \quad \omega \, = \sqrt { \frac { 2 \left | h \right | } { \mu } }.$$

Remark 1.2 . According to the value of the energy h and the angular momentum k , six cases can appear. However in this paper we are interested only in the bounded non-collision motions which correspond to the case of non-zero angular momentum and negative energy.

It is well-known, in fact, that all solutions can be written in terms of the orbital elements and in the particular case of non-zero angular momentum and negative energy, such solutions are ellipses given by

$$( 1. 1 8 ) \ r ( \tau ) = a [ 1 - e \cos \omega \tau ], \quad \text{where} \quad \sqrt { \frac { m } { \mu a ^ { 3 } } } \, t = \frac { 2 \pi } { T } t = \tau - e \sin \tau, \ \ a \text{ is the semi-major axis}$$

and where ε = [1 + 2 hk / µm 2 ( 2 )] 1 / 2 denotes the eccentricity. We observe that in the circular case ε = 0 and r τ ( ) = a is indeed constant. It is also well-known that in polar coordinates ( r, ϑ ) , the polar equation of the ellipses is given by

$$r ( \vartheta ) = \frac { r _ { 0 } } { 1 - e \cos \vartheta } \quad \text{where} \quad r _ { 0 } \, = \frac { k ^ { 2 } } { \mu m }.$$

r 0 is called semi-latus rectum and it is related to the eccentricity and the semi-major axis of the ellipses by a = r / 0 (1 -e 2 ) .

## 2 Maslov index and Conley-Zehnder index

This section is devoted to some classical definitions and basic properties of the Conley-Zehnder index. More precisely we need its generalization to paths with degenerate endpoints in the 2 and 4 dimensional symplectic space as well as a description of the Maslov index as an intersection index in the Lagrangian Grassmannian setting.

## 2.1 On the Conley-Zehnder index

We consider the standard symplectic space ( R 2 n , ω 0 ) where the (standard) symplectic form ω is defined by ω ( · , · ) = 〈 J , · ·〉 and where J is the 2 n × 2 n matrix given by J = [ 0 -I I 0 ] . We denote by Sp(2 n ) the symplectic group defined by Sp(2 n ) = { M ∈ GL(2 n ) ∣ ∣ M JM T = J } and for 0 glyph[lessorequalslant] k glyph[lessorequalslant] 2 n we set Sp (2 k n ) = { M ∈ Sp(2 n ) | dimker( M -I ) = k } . Thus, in particular, we get that

$$\text{Sp} ( 2 n ) = \bigcup _ { k = 0 } ^ { 2 n } \text{Sp} _ { k } ( 2 n ).$$

Given M ∈ Sp(2 n ) , we set I ( M ) = ( -1) n -1 det( M -I ) and we define the following hypersurface of

$$S p ^ { 0 } ( 2 n ) = \{ \, M \in S p ( 2 n ) \, | \, I ( M ) = 0 \, \} \subset S p ( 2 n ).$$

Now, setting Sp(2 n ) ∗ = Sp(2 n ) \ Sp (2 0 n ) and denoting by

$$\text{Sp} ^ { + } ( 2 n ) \, = \, \left \{ \, M \in \text{Sp} ( 2 n ) \, | \, I ( M ) < 0 \, \right \} \quad \text{and} \quad \text{Sp} ^ { - } ( 2 n ) \, = \, \left \{ \, M \in \text{Sp} ( 2 n ) \, | \, I ( M ) > 0 \, \right \}$$

it is possible to prove that these two are the only two path connected components of Sp(2 n ) ∗ and are simply connected in Sp(2 n ) .

↦

For any M ∈ Sp (2 0 n ) , we can define a transverse orientation at M through the positive direction d dt Me tJ | t =0 of the path t → Me tJ with t glyph[greaterorequalslant] 0 sufficiently small. We define the following set

$$\mathcal { P } _ { T } ( 2 n ) \, = \left \{ \, \psi \in \mathcal { C } ^ { 0 } \left ( [ 0, T ] ; \text{Sp} ( 2 n ) \right ) \, | \, \psi ( 0 ) = I \, \right \}.$$

In the following will be essentially interested to the cases P T (2) and P T (4) .

## The Conley-Zehnder index in Sp(4)

Since later on we shall work in Sp(4) , in this section we restrict to Sp(4) and Sp(2) in order to simplify the presentation.

Consider now the two square matrices 2 × 2 matrices M 1 and M 2 given by

$$M _ { k } \, = \, \begin{pmatrix} a _ { k } & b _ { k } \\ c _ { k } & d _ { k } \end{pmatrix} \in \text{Sp} ( 2 ) \quad \text{for} \quad k = 1, 2,$$

where a , b k k , c k , d k ∈ R . The symplectic sum of M 1 and M 2 is defined as the following 4 × 4 matrix below

$$M _ { 1 } \diamond M _ { 2 } \, = \, \begin{pmatrix} a _ { 1 } & 0 & b _ { 1 } & 0 \\ 0 & a _ { 2 } & 0 & b _ { 2 } \\ c _ { 1 } & 0 & d _ { 1 } & 0 \\ 0 & c _ { 2 } & 0 & d _ { 2 } \end{pmatrix} \in \text{Sp} ( 4 ).$$

The 2 -fold symplectic sum of M with itself is denoted by M glyph[diamondmath] 2 . The symplectic sum of two paths ψ j ∈ P T (2) , with j = 1 2 , , is defined by:

$$( \psi _ { 1 } \diamond \psi _ { 2 } ) ( t ) \, = \psi _ { 1 } ( t ) \diamond \psi _ { 2 } ( t ), \ \forall \, t \in [ 0, T ].$$

Given any two continuous paths φ , φ 1 2 : [0 , T ] → Sp(4) such that φ 1 ( T ) = φ 2 (0) , we denote by ∗ their concatenation . We also define a special continuous symplectic path ξ 2 : [0 , T ] → Sp(4) as follows:

$$( 2. 3 ) \quad \xi _ { 2 } ( t ) \, = \, \begin{bmatrix} 2 - \frac { t } { T } & 0 \\ 0 & \left ( 2 - \frac { t } { T } \right ) ^ { - 1 } \end{bmatrix} ^ { \otimes 2 } \, = \, \begin{bmatrix} \left ( 2 - \frac { t } { T } \right ) I & 0 \\ 0 & \left ( 2 - \frac { t } { T } \right ) ^ { - 1 } I \end{bmatrix} \,, \quad \forall \, t \in [ 0, T ]$$

where in the righthand side I denotes the identity 2 × 2 . Setting D a ( ) = diag[ a, a -1 ] , we define the two matrices M + 2 = D (2) glyph[diamondmath] D (2) ∈ Sp + (4) and M -2 = D ( -2) glyph[diamondmath] D (2) ∈ Sp -(4) . Given ψ ∈ P T (2) and m ∈ N \ { 0 } , the m -th iteration of ψ is ψ m : [0 , mT ] → Sp(2) defined as

$$\psi ^ { m } ( t ) \, = \psi ( t - j T ) ( \psi ( T ) ) ^ { j }, \quad \text{for } j T \leqslant t \leqslant ( j + 1 ) T, \ \ j = 0, \dots, m - 1.$$

Remark 2.1 . All definitions given in this subsection in dimension 2 or 4 can be carried over in any even dimension.

Wenow introduce the definition of ι CZ -index which is an intersection index between a symplectic path starting from identity and the singular hypersurface defined above.

Definition 2.2. We consider a continuous path ψ : [0 , T ] → Sp(4) such that ψ (0) = I . We define the Conley-Zehnder index of the path ψ as the integer given by

$$\iota ^ { C z } ( \psi ( t ), t \in [ 0, T ] ) \, = \left [ e ^ { - \varepsilon J } \psi * \xi _ { 2 } \, \colon \text{Sp} ^ { 0 } ( 4 ) \right ]$$

where the righthand side in Equation (2.4) is the usual homotopy intersection number and ε is a positive real sufficiently small number .

↦

Remark 2.3 . It is worth to observe that the advantage to perturb the path t → ( ψ ∗ ξ 2 )( ) t is in order to get a new path having nondegenerate endpoints (and in particular of nondegenerate endpoint which is the original assumption made for defining the Conley-Zehnder index). We also notice that the advantage to concatenate the original path ψ with ξ 2 is in order to simplify some relations in the specific case in which the symplectic path ψ is the fundamental solution of a Hamiltonian system whose Hamiltonian is the Fenchel transform of a Lagrangian one.

## Properties of the ι CZ -index

We list some of the the basic properties of the index that we need to use later on.

- (i) ( glyph[diamondmath] -additivity) Let ψ 1 : [0 , T ] → Sp(2) and ψ 2 : [0 , T ] → Sp(2) be two symplectic paths. Then we have

$$\iota ^ { \text{CZ} } ( ( \psi _ { 1 } \diamond \psi _ { 2 } ) ( t ), t \in [ 0, T ] ) = \iota ^ { \text{CZ} } ( \psi _ { 1 } ( t ), t \in [ 0, T ] ) + \iota ^ { \text{CZ} } ( \psi _ { 2 } ( t ), t \in [ 0, T ] ).$$

- (ii) (Homotopy invariance) For any two paths ψ 1 and ψ 2 , if ψ 1 ∼ ψ 2 in Sp(2) with either fixed or always non-degenerate endpoints, there holds

$$\iota ^ { \text{cz} } ( \psi _ { 1 } ( t ), t \in [ 0, T ] ) = \iota ^ { \text{cz} } ( \psi _ { 2 } ( t ), t \in [ 0, T ] ).$$

- (iii) (Affine scale invariance) For all k &gt; 0 and ψ ∈ P kT (4) , we have

$$\iota ^ { \text{cz} } ( \psi ( k t ), t \in [ 0, T ] ) = \iota ^ { \text{cz} } ( \psi ( t ), t \in [ 0, k T ] ).$$

Figure 1: The singular surface Sp(2) 0 1 . The representation is in cylindrical coordinates ( x, y, z ) = ( r cos θ, r sin θ, z ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00075v2-figure-12122e07-7eeb-4181-8e32-8d4c61e93660.jpg)

## The geometric structure of Sp(2)

The symplectic group Sp(2) captured the attention of I. Gelfand and V. Lidskii first, who in 1958 described a toric representation of it. The R 3 -cylindrical coordinate representation of Sp(2) , that we shall use throughout, was introduced by Y. Long in 1991. Every real invertible matrix A can be decomposed in polar form

$$A = P O,$$

where P = ( AA T ) 1 / 2 is symmetric and positive definite and O is a proper rotation:

$$O = \begin{pmatrix} \cos \theta & - \sin \theta \\ \sin \theta & \cos \theta \end{pmatrix}.$$

In particular, the matrix P can be written in the following form (cfr. [BJP16, Appendix A ] and references therein for further details )

$$P = \begin{pmatrix} r & z \\ z & \frac { 1 + z ^ { 2 } } { r } \end{pmatrix}$$

and then every symplectic 2 × 2 matrix M can be written as the product

$$M & = \begin{pmatrix} r & z \\ z & \frac { 1 + z ^ { 2 } } { r } \end{pmatrix} \begin{pmatrix} \cos \theta & - \sin \theta \\ \sin \theta & \cos \theta \end{pmatrix},$$

where ( r, θ, z ) ∈ (0 , + ∞ × ) [0 , 2 π ) × R . Viewing ( r, θ, z ) as cylindrical coordinates in R 3 \ { z -axis } we obtain a smooth global diffeomorphism Ψ : Sp(2) → R 3 \{ z -axis } . We shall henceforth identify elements in Sp(2) with their image under Ψ .

The eigenvalues of a symplectic matrix M written as in Equation (2.5) are

$$\lambda _ { \pm } \, = \, \frac { 1 } { 2 r } \left [ ( 1 + r ^ { 2 } + z ^ { 2 } ) \cos \theta \pm \sqrt { ( 1 + r ^ { 2 } + z ^ { 2 } ) ^ { 2 } \cos ^ { 2 } \theta - 4 r ^ { 2 } } \right ].$$

Thus, we get

$$I ( M ) = 2 - \left ( r + \frac { 1 + z ^ { 2 } } { r } \right ) \cos \theta$$

and define

$$S p ^ { \pm } ( 2 ) = \left \{ \, ( r, \theta, z ) \in ( 0, + \infty ) \times [ 0, 2 \pi ) \times \mathbb { R } \, \left | \, \pm ( 1 + r ^ { 2 } + z ^ { 2 } ) \cos \theta > 2 r \, \right \},$$

$$\text{Sp} ^ { 0 } ( 2 ) = \left \{ \, ( r, \theta, z ) \in ( 0, + \infty ) \times [ 0, 2 \pi ) \times \mathbb { R } \, \Big | \, \pm ( 1 + r ^ { 2 } + z ^ { 2 } ) \cos \theta = 2 r \, \right \}.$$

$$& \times \mathbb { R } \, \Big | \, \pm ( 1 + r ^ { 2 } + z ^ { 2 } ) \cos \theta > 2 r \, \Big \} \,, \\ & \times \mathbb { R } \, \Big | \, \pm ( 1 + r ^ { 2 } + z ^ { 2 } ) \cos \theta = 2 r \, \Big \} \,.$$

The set Sp (2) = Sp ∗ + (2) ∪ Sp -(2) is named the regular part of Sp(2) , while Sp (2) 0 is its singular part ; the former corresponds to the subset of 2 × 2 symplectic matrices which do not have 1 as an eigenvalue, whereas those matrices admitting 1 in their spectrum belong to the latter. We are particularly interested in Sp (2) 0 , the singular part of Sp(2) associated with the eigenvalue 1 , a representation of which is depicted in Figure 1. The 'pinched' point is the identity matrix, and it is the only element satisfying dimker( M -I ) = 2 . If we denote by

$$\text{Sp} ( 2 ) _ { \pm } ^ { 0 } \, = \left \{ \, ( r, \theta, z ) \in \text{Sp} ^ { 0 } ( 2 ) \, \Big | \, \pm \sin \theta > 0 \, \right \},$$

we see that Sp (2) 0 \ { I } = Sp(2) 0 + ∪ Sp(2) 0 -, and each subset is a path-connected component diffeomorphic to R 2 \ { 0 } . In Figure 2, it is represented as a horizontal section (specifically the

Figure 2: Intersection of Sp (2) 0 with the plane z = 0 . The representation is in Cartesian coordinates ( x, y ) = ( r cos θ, r sin θ ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00075v2-figure-926fa93c-6703-4d41-913b-bdc401297b9c.jpg)

intersection with the plane z = 0 ) of the hypersurface Sp (2) 0 .

What the ι CZ index of a symplectic path in Sp(2) starting from the identity actually counts are the algebraic (signed) intersections of the path with the surface Sp (2) 0 . If we imagine projecting the symplectic path into the horizontal plane { z = 0 } a simple way to think about this index is an algebraic count of intersection with the curve depicted in Figure 2of a continuous path (which is the projection of the original one) on the ̂ xOy -plane. (For further details, we refer the interested reader to [BJP16] and references therein). We observe that, in these coordinate system the (unbounded) path connected component Sp + (2) depicted in Figure 2, corresponds to the region of the plane containing the portion of the x -axis emanating from the identity.

## 2.2 The Maslov index in Λ(2)

In the 4 -dimensional standard symplectic space ( R 4 , ω 0 ) a Lagrangian subspace is an 2 -dimensional subspace L ⊂ R 4 on which ω 0 vanishes identically. It is well-known [Dui76] that the set of all Lagrangian subspaces Λ(2) , usually called the Lagrangian Grassmannian of the 2 -dimensional Lagrangian subspaces has the structure of a three dimensional compact, real-analytic submanifold of the Grassmannian of all 2 -dimensional subspaces of R 4 . The real-analytic atlas can be described as follows. For all L ∈ Λ(2) and k ∈ 0 1 , , let Λ ( 0 L ) denote the dense and open set of all Lagrangian subspaces L ′ such that L ′ ∩ L 0 = { 0 } . Given L 1 ∈ Λ(2) , let L 0 ∈ Λ( L 1 ) and let us consider the map

$$\ast \quad \varphi _ { L _ { 0 }, L _ { 1 } } \colon \Lambda _ { 0 } ( L _ { 1 } ) \to \text{B} _ { \text{sym} } ( L _ { 0 } ) \text{ defined by } \varphi _ { L _ { 0 }, L _ { 1 } } ( L ) = \omega _ { 0 } ( \cdot, T \cdot ) | _ { L _ { 0 } \times L _ { 0 } }$$

where T : L 0 → L 1 denotes the unique linear map whose graph is the Lagrangian (subspace) L . 1 Let L 0 ∈ Λ(2) and for j ∈ { 0 1 2 , , } , we set

$$\Lambda _ { j } ( L _ { 0 } ) & = \left \{ \, L \in \Lambda ( 2 ) \, | \, \dim ( L \cap L _ { 0 } ) = j \, \right \}.$$

It is easy to see that Λ ( j L 0 ) is a connected j ( j +1) 2 / -codimensional submanifold of Λ(2) and in particular Λ ( 2 L 0 ) = L 0 . Moreover, the set

$$\Sigma ( L _ { 0 } ) & = \Lambda _ { 1 } ( L _ { 0 } ) \cup \Lambda _ { 2 } ( L _ { 0 } )$$

↦

is the (topological) closure of the top stratum Λ ( 1 L 0 ) usually called the Maslov cycle . Λ ( 1 L 0 ) has a canonical transverse orientation, meaning that there exists δ &gt; 0 for each L ∈ Λ ( 1 L 0 ) , the path of Lagrangian subspaces t → e tJ L for t ∈ ( -δ, δ ) crosses transversally Λ ( 1 L 0 ) and as t increases the path is pointing towards the transverse positive direction. Thus this cycle is two-sidedly embedded in Λ(2) . Following authors in [CLM03] we introduce the following definition.

Definition 2.4. Let L 0 ∈ Λ(2) and let glyph[lscript] : [0 , T ] → Λ(2) be a continuous path. We term ι CLM -index the integer defined by

$$\iota ^ { C L M } ( L _ { 0 }, \ell ( t ), t \in [ 0, T ] ) \, = \, \left [ e ^ { - \epsilon J } \, \ell ( t ) \, \colon \Sigma ( L _ { 0 } ) \right ]$$

where the right-hand side is the intersection number and glyph[epsilon1] ∈ (0 1) , is sufficiently small.

Remark 2.5 . We observe that the ι CLM -index given in Definition 2.4 could be defined by using the Seifert-Van Kampen theorem for groupoids. More precisely, we denote throughout by π (Λ(2)) the fundamental groupoid of Λ(2) , namely the set of fixed-endpoints homotopy classes [ glyph[lscript] ] of continuous paths γ in Λ(2) , endowed with the partial operation of concatenation ∗ . For all L 0 ∈ Λ(2) , there exists a unique Z / 2 -valued groupoid homomorphism ι RS : π (Λ(2)) → Z / 2 such that

$$\iota ^ { \text{RS} } ( \ell ( t ), L _ {$$

for all continuous curve glyph[lscript] : [0 , 1] → Λ ( 0 L 1 ) and for all L 1 ∈ Λ ( 0 L 0 ) . 2 By [LZ00, Equation (3.7)] we get that

$$\mathcal { v$$

where h i ( ) = dim ( L 0 ∩ glyph[lscript] ( ) i ) for i = 0 1 , . Thus locally the ι CLM -index with respect to the fixed Lagrangian L 0 could be defined equally well as the unique Z -valued groupoid homomorphism ι RS : π (Λ(2)) → Z / 2 such that

$$\iota ^ { \text{RS} } ( \ell ( t ), L _ {$$

A different choice has been considered by authors in [GPP04].

Remark 2.6 . It is well-known that the Lagrangian Grassmannian can be realized as homogeneous space, through an action of the unitary group. In fact it is easy to show that

$$\Lambda ( 2 ) = U ( 2 ) / O ( 2 ).$$

As proved by Arnol'd in [Arn86, Section 3], Λ(2) is the nonoriented total space of the nontrivial bundle with fiber S 2 and with base the circle. The proof provided by Arnol'd is based on the short exact sequences of fibrations of the unitary and orthogonal Lie groups.

1 The symmetry of the restriction of the bilinear map ω 0 ( T , · · ) onto L 0 × L 0 is consequence of the fact that L is Lagrangian.

2 This index was defined in a slightly different manner by authors in [RS93].

## Computing ι CLM -index through crossing forms

The ι CLM -index defined above, is in general, quite hard to compute. However one efficient way to do so is via crossing forms as introduced by authors in [RS93]. Let glyph[lscript] : [0 , T ] → Λ(2) be a C 1 -curve of Lagrangian subspaces and let L 0 ∈ Λ(2) . Fix t ∈ [0 , T ] and let W be a fixed Lagrangian complement of glyph[lscript] ( ) t . If s belongs to a suitable small neighborhood of t for every v ∈ l ( ) t we can find a unique vector w s ( ) ∈ W in such a way that v + w s ( ) ∈ l ( s ) .

Definition 2.7. The crossing form Γ( glyph[lscript], L , t 0 ∗ ) at t ∗ is the quadratic form Γ( glyph[lscript], L , t 0 ∗ ) : l ( t ∗ ) ∩ L 0 → R defined by

$$\Gamma ( l, L _ { 0 }, t _ { * } ) [ v ]$$

glyph[negationslash]

The number t ∗ is said to be a crossing instant for l with respect to L 0 if l ( t ∗ ) ∩ L 0 = 0 { } and it is called regular if the crossing form is non-degenerate.

Let us remark that regular crossings are isolated and hence on a compact interval there are finitely many. The following result is well-known.

Proposition 2.8. ([LZ00, Theorem 3.1]) Let L 0 ∈ Λ(2) and glyph[lscript] ∈ C 1 ( [0 , T ] , Λ(2) ) having only regular crossings. Then the ι CLM -index of glyph[lscript] with respect to L 0 is given by

$$\iota ^ { \text{CLM} } ( L _ { 0 }, \ell, [ 0, T ] ) \, = \, \mathfrak { n } _ { + } \, \Gamma ( \ell, L _ { 0 }, 0 ) + \sum _ { t _ { * } \in ( 0, T ) } \text{sgn} \, \Gamma ( \ell, L _ { 0 }, t _ { * } ) - \mathfrak { n } _ { - } \, \Gamma ( \ell, L _ { 0 }, T )$$

where the summation runs over all crossings instants t ∗ ∈ (0 , T ) .

On the Euclidean space V = R 2 × R 2 , we introduce the symplectic form ω J : V × V → R defined by

$$\psi \psi _ { \vec { \lambda } ^ { \prime } } \quad \omega _ { \vec { \vartheta } } ( v, w ) \, = \, \langle \vec { \vartheta } v, w \rangle, \quad \forall \, v, w \in V \quad \text{where} \quad \vec { \vartheta } \, = \, \begin{pmatrix} - J & 0 \\ 0 & J \end{pmatrix}.$$

↦

↦

The next result (in the general setting) was proved by authors in [LZ00, Corollary 2.1] (cfr. [HS09, Lemma 4.6]) and in particular put on evidence the relation between the ι CZ -index associated to a path of symplectic matrices and the ι CLM -index of the corresponding path of Lagrangian subspaces with respect to the diagonal ∆ = Gr( ) I .

By a direct calculation it follows that, if M ∈ Sp(2) then its graph Gr( M ) = { ( x, Mx ) T | x ∈ R 2 } is a Lagrangian subspace of the symplectic space ( V, ω J ) . Thus, a path of symplectic matrices t → ψ t ( ) induces a path of Lagrangian subspaces of ( V, ω J ) defined through its graph by t → Gr( ψ t ( )) .

Proposition 2.9. For any continuous symplectic path ψ ∈ P T (2 n ) (thus starting at the identity), we get the following equality

$$\iota ^ { \text{CZ} } ( \psi ( t ), t \in [ 0, T ] ) + n = \iota ^ { \text{CLM} } ( \Delta, \text{Gr} ( \psi ( t ) ), t \in [ 0, T ] )$$

In the special case n = 1 , if φ ∈ P T (2) , it holds that

$$\iota ^ { \text{CZ} } ( \phi ( t ), t \in [ 0, T ] ) + 1 = \iota ^ { \text{CLM} } ( \Delta, \text{Gr} ( \phi ( t ) ), t \in [ 0, T ] ).$$

We conclude this section with an index theorem which relates the Morse index ι ( γ ) of a periodic solution γ of Equation (1.5) seen as critical point of the Lagrangian action functional and the ι CZ -index of the fundamental solution ψ of the linearized Hamiltonian system at γ .

Proposition 2.10 (Morse Index Theorem, [Lon02, page 172]) . Under the above notation, we have

$$\iota ( \gamma ) = \iota ^ { \text{cz} } ( \psi ( t ), t \in [ 0, T ] ).$$

## 3 Some explicit computations for paths in Sp(2)

The aim of this subsection is to explicitly compute the Conley-Zehnder index as well as the ι CLM -index introduced in Subsection 1.2 using crossing forms.

Lemma 3.1. Let φ : [0 , T ] → Sp(2) be the path

$$\phi ( t ) \, = \, \begin{pmatrix} a ( t ) & b ( t ) \\ c ( t ) & d ( t ) \end{pmatrix}, \quad t \in [ 0, T ]$$

glyph[negationslash]

with a, b, c, d ∈ C 1 ([0 , T ] , R ) and let glyph[lscript] ( ) t = Gr ( φ t ( ) ) ∈ Λ(2) . We assume that t ∗ ∈ [0 , T ] is a crossing instant for glyph[lscript] (with respect to ∆ ) such that glyph[lscript] ( t ∗ ) ∩ ∆ = { 0 } . Assuming that d t ( ∗ ) = 0 , then the crossing form at t = t ∗ is given by glyph[negationslash]

$$\begin{array} { c c c c c c c c c } \psi \, \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi & \psi \\ ( 3. 1 ) & \xi ^ { \prime } ( t _ { * } ) = a ^ { \prime } ( t _ { * } ) x _ { 0 } + b ^ { \prime } ( t _ { * } ) y _ { 0 } - \frac { b ( t _ { * } ) } { d ( t _ { * } ) } \left [ c ^ { \prime } ( t _ { * } ) x _ { 0 } + d ^ { \prime } ( t _ { * } ) y _ { 0 } \right ], & \quad \forall \, v \, = ( x _ { 0 }, y _ { 0 }, x _ { 0 }, y _ { 0 } ) \in \Delta. \\ - & \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

Proof. In order to compute the crossing form given in Equation (2.11), we first consider the Lagrangian subspace

$$W \, = \{ 0 \} \times \mathbb { R } \times \mathbb { R } \times \{ 0 \}$$

and we observe that this gives a Lagrangian decomposition of R 4 , specifically R 4 = ∆ ⊕ W . Now, for any v = ( x , y 0 0 , x 0 , y 0 ) ∈ ∆ let us choose w t ( ) = (0 , η ( ) t , ξ ( ) t , 0) ∈ W in order that v + w t ( ) ∈ l ( ) t . This means that η t ( ) and ξ t ( ) solve the equations

$$( 3. 2 ) \quad \ \ x _ { 0 } + \xi ( t ) = a ( t ) x _ { 0 } + b ( t ) \left ( y _ { 0 } + \eta ( t ) \right ), \quad y _ { 0 } = c ( t ) x _ { 0 } + d ( t ) \left ( y _ { 0 } + \eta ( t ) \right ).$$

Since in a crossing instant t ∗ we have ξ t ( ∗ ) = η t ( ∗ ) = 0 , differentiating the above identities gives

$$\xi ^ { \prime } ( t _ { * } ) & = a ^ { \prime } ( t _ { * } ) x _ { 0 } + b ^ { \prime } ( t _ { * } ) y _ { 0 } - \frac { b ( t _ { * } ) } { d ( t _ { * } ) } \left [ c ^ { \prime } ( t _ { * } ) x _ { 0 } + d ^ { \prime } ( t _ { * } ) y _ { 0 } \right ], \\ \eta ^ { \prime } ( t _ { * } ) & = - \frac { 1 } { d ( t _ { * } ) } \left [ c ^ { \prime } ( t _ { * } ) x _ { 0 } + d ^ { \prime } ( t _ { * } ) y _ { 0 } \right ].$$

By a direct computation we obtain

$$( 3. 4 ) \ \omega _ { \hat { g } } \left ( v, w ( t ) \right ) = \left \langle \hat { y } v, w ( t ) \right \rangle = - \left \langle J \left ( _ { y _ { 0 } } ^ { x _ { 0 } } \right ), \left ( _ { \tau ( t ) } ^ { 0 } \right ) \right \rangle + \left \langle J \left ( _ { y _ { 0 } } ^ { x _ { 0 } } \right ), \left ( _ { 0 } ^ { \xi ( t ) } \right ) \right \rangle = - x _ { 0 } \eta ( t ) - y _ { 0 } \xi ( t ).$$

Hence the crossing form at the crossing instant t = t ∗ is given by

$$\Gamma ( \ell, \Delta, t _ { * } ) ( v ) = \frac { d } { d t } \omega _ { \mathcal { J } } \left ( v, w ( t ) \right ) \Big | _ { t = t _ { * } } = - x _ { 0 } \eta ^ { \prime } ( t _ { * } ) - y _ { 0 } \xi ^ { \prime } ( t _ { * } ).$$

$$\Gamma ( \ell, \Delta, t _ { * } ) ( v ) = \frac { d } { d t } \omega _ { \beta } ( v, w ( t ) ) \Big | _ { t = t _ { * } } = - x _ { 0 } \eta ^ { \prime } ( t _ { * } ) - y _ { 0 } \xi ^ { \prime } ( t _ { * } ).$$

↦

Remark 3.2 . If d t ( ∗ ) = 0 , then it is enough to replace the path t → glyph[lscript] ( ) t by the path t → glyph[lscript] ε ( ) t = Gr( e -εJt φ t ( )) with ε &gt; 0 sufficiently small. By the well-definedness of the ι CLM -index (cfr. [CLM03, pag. 138]) the result follows.

↦

We are going to apply the computation provided in Lemma 3.1 in some specific cases that we shall need later.

Example 3.3. We let a , a 1 2 ∈ R + , β = √ a a 1 2 and we consider the path R β : [0 , T ] → Sp(2) defined by

$$R _ { \beta } ( t ) = \begin{bmatrix} \cos ( \beta t ) & - \frac { \beta } { a _ { 1 } } \sin ( \beta t ) \\ \frac { \beta } { a _ { 2 } } \sin ( \beta t ) & \cos ( \beta t ) \end{bmatrix}.$$

↦

We aim to compute the crossing form with respect to ∆ of the Lagrangian path t → glyph[lscript] β ( ) t = Gr ( R β ( ) t ) . Bearing in mind previous notation, we get

$$a ( t ) = d ( t ) = \cos ( \beta t ), \ \ b ( t ) = - \frac { \beta } { a _ { 1 } } \sin ( \beta t ), \quad c ( t ) = \frac { \beta } { a _ { 2 } } \sin ( \beta t ).$$

We observe that t ∗ is a crossing instant if and only if t ∗ ∈ 2 π β Z . Moreover by a direct calculation we get

$$( 3. 6 ) \\ a ( t _ { * } ) = d ( t _ { * } ) = 1, \quad b ( t _ { * } ) = c ( t _ { * } ) = 0.$$

## Moreover,

$$a ^ { \prime } ( t _ { * } ) & = - \beta \sin ( \beta t _ { * } ) = d ^ { \prime } ( t _ { * } ), \quad b ^ { \prime } ( t _ { * } ) = - a _ { 2 } \cos ( \beta t _ { * } ), \quad c ^ { \prime } ( t _ { * } ) = a _ { 1 } \cos ( \beta t _ { * } ) \\ a ^ { \prime } ( t _ { * } ) & = 0 = d ^ { \prime } ( t _ { * } ) \quad b ^ { \prime } ( t _ { * } ) = - a _ { 2 }, \quad c ^ { \prime } ( t _ { * } ) = a _ { 1 }.$$

Thus ξ ′ ( t ∗ ) = -a y 2 0 and η ′ ( t ∗ ) = -a x 1 0 . Using Lemma (3.1), we directly get

$$\Gamma ( \ell _ { \beta }, \Delta, t _ { * } ) [ v ] = - x _ { 0 } \eta ^ { \prime } ( t _ { * } ) - y _ { 0 } \xi ^ { \prime } ( t _ { * } ) = a _ { 2 } x _ { 0 } ^ { 2 } + a _ { 1 } y _ { 0 } ^ { 2 }.$$

Since Γ( glyph[lscript] β , ∆ , t ∗ ) is a positive definite quadratic form on a 2 dimensional vector space, it follows that it is non-degenerate and its signature (which coincides with the coindex) is 2 . Summing up all these computations we obtain

$$\ell ^ { C L M } ( \Delta, \ell _ { \beta } ( t ), t \in [ 0, T ] ) = \begin{cases} 2 \left \lfloor \frac { T \beta } { 2 \pi } \right \rfloor & \text{if $T \in \frac{2\pi}{\beta} \mathbb{Z}$} \\ 2 \left \{ \left \lfloor \frac { T \beta } { 2 \pi } \right \rfloor + 1 \right \} & \text{if $T \notin \frac{2\pi}{\beta} \mathbb{Z}$}, \end{cases}$$

where glyph[floorleft]·glyph[floorright] denotes the greatest integer less than or equal to its argument.

Remark 3.4 . Let us consider the path S β : [0 , T ] → Sp(2) defined by

$$S _ { \beta } ( t ) = \begin{bmatrix} \cos ( \beta t ) & \frac { \beta } { a _ { 1 } } \sin ( \beta t ) \\ - \frac { \beta } { a _ { 2 } } \sin ( \beta t ) & \cos ( \beta t ) \end{bmatrix}$$

↦

and let us define the Lagrangian path t → m β ( ) = Gr t ( S β ( ) t ) . By the very same calculations as before, we get

$$\iota ^ { \text{CLM} } ( \Delta, m _ { \beta } ( t ), t \in [ 0, T ] ) = \begin{cases} - 2 \left \lfloor \frac { T \beta } { 2 \pi } \right \rfloor & \text{ if } T \notin \frac { 2 \pi } { \beta } \mathbb { Z } \\ - 2 \left \{ \left \lfloor \frac { T \beta } { 2 \pi } \right \rfloor + 1 \right \} & \text{ if } T \in \frac { 2 \pi } { \beta } \mathbb { Z }, \end{cases}$$

Summing up we have the following result.

Lemma 3.5. Let us consider the path

↦

$$\psi \colon [ 0, T ] \longrightarrow \text{Sp} ( 2 ) \, \colon t \longmapsto e ^ { t J S }$$

where S is either symmetric positive or negative definite and for every t ∈ [0 , T ] , we let n t ( ) = Gr ( ψ t ( ) ) . Thus, we get (3.10) √

$$\tau ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ 0, T ] ) = \begin{cases} 2 \left \lfloor \frac { \sqrt { a _ { 1 } a _ { 2 } } T } { 2 \pi } \right \rfloor & \text{if $S$ is positive definite and $T\in\frac{2\pi}{\sqrt{a_{1}a_{2}}Z$} } \\ - 2 \left \lfloor \frac { \sqrt { a _ { 1 } a _ { 2 } } T } { 2 \pi } \right \rfloor & \text{if $S$ is negative definite and $T\notin\frac{2\pi}{\sqrt{a_{1}a_{2}}Z$} } \\ 2 \left \{ \left \lfloor \frac { \sqrt { a _ { 1 } a _ { 2 } } T } { 2 \pi } \right \rfloor + 1 \right \} & \text{if $S$ is positive definite and $T\notin\frac{2\pi}{\sqrt{a_{1}a_{2}}Z$} } \\ - 2 \left \{ \left \lfloor \frac { \sqrt { a _ { 1 } a _ { 2 } } T } { 2 \pi } \right \rfloor + 1 \right \} & \text{if $S$ is negative definite and $T\in\frac{2\pi}{\sqrt{a_{1}a_{2}}Z$} } \end{cases} \\ \text{where} \, a _ { 1 } \text{ and } a _ { 2 } \text{ are the eigenvalues of $S$}.$$

where a 1 and a 2 are the eigenvalues of S .

Proof. Since S is symmetric, we diagonalize it in the orthogonal group and we get a symplectic basis in R 2 of the folloiwing matrices:

$$S & = \begin{bmatrix} a _ { 1 } & 0 \\ 0 & a _ { 2 } \end{bmatrix} \quad \text{ and } J S = \begin{bmatrix} 0 & - a _ { 2 } \\ a _ { 1 } & 0 \end{bmatrix}.$$

We start by setting

$$\beta \, = \sqrt { a _ { 1 } a _ { 2 } } \quad \text{if} \quad a _ { 1 } a _ { 2 } > 0$$

and we observe that by the definition of the matrix exponential we have the following expression for ψ :

## 1. If the sign S &gt; 0 ( a 1 &gt; 0 and a 2 &gt; 0 ), then we get

$$\psi ( t ) = \begin{bmatrix} \cos ( \beta t ) & - \frac { \beta } { a _ { 1 } } \sin ( \beta t ) \\ \frac { \beta } { a _ { 2 } } \sin ( \beta t ) & \cos ( \beta t ) \end{bmatrix}$$

## 2. If the sign S &lt; 0

$$\psi ( t ) = \begin{bmatrix} \cos ( \beta t ) & \frac { \beta } { a _ { 1 } } \sin ( \beta t ) \\ - \frac { \beta } { a _ { 2 } } \sin ( \beta t ) & \cos ( \beta t ) \end{bmatrix}$$

The proof now follows by invoking Equation (2.12), Lemma 3.1 and Example 3.3. This concludes the proof.

As a consequence of Lemma 3.5 and the homotopy invariance of the ι CLM -index we get the following result.

Proposition 3.6. Let us consider the path pointwise defined by ψ ε ( ) t = e -εJ ψ t ( ) , where ψ is given in Lemma 3.5 and let n ε ( ) t = Gr ( ψ ε ( ) t ) . Then, for ε &gt; 0 sufficiently small, we get

$$( 3. 1 2 ) & \iota ^ { \text{CLM} } ( \Delta, n _ { \varepsilon } ( t ), t \in [ 0, T ] ) = \iota ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ 0, T ] ).$$

Proof. This result is a direct consequence of the well definedness of the ι CLM -index. (Cfr.[CLM03, pag. 138]).

Another situation which occurs often in the applications when there is the presence of a conservation law, it is described in the following result.

Lemma 3.7. Let N : [0 , T ] → Sp(2) be the path pointwise defined by

$$N ( t ) \, = \, \begin{pmatrix} 1 & f ( t ) \\ 0 & 1 \end{pmatrix}, \quad t \in [ 0, T ]$$

with either f ( ) = t t or f ( ) = t -t . We denote by n the induced path of Lagrangian subspaces in R 4 pointwise defined by n t ( ) = Gr ( N t ( ) ) . Then we have

$$\iota ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ 0, T ] ) = \begin{cases} 0 & \text{if $f(t)=t$} \\ 1 & \text{if $f(t)=-t$.} \end{cases}$$

Moreover

$$\iota ^ { \text{cz} } ( N ( t ), t \in [ 0, T ] ) = \begin{cases} - 1 & \text{if $f(t)=t$} \\ 0 & \text{if $f(t)=-t$.} \end{cases}$$

Figure 3: c

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00075v2-figure-c5cc25bf-8952-4415-9734-6dfe8a4a01ad.jpg)

↦

Proof. By the well definedness of the ι CLM -index, (cfr.[CLM03, pag. 138]), it is enough to compute the ι CLM -index for the path t → n ε ( ) t where n ε ( ) = Gr( t N ε ( )) t for

$$N _ { \varepsilon } ( t ) \coloneqq \begin{pmatrix} \cos \varepsilon & f ( t ) \cos \varepsilon + \sin \varepsilon \\ - \sin \varepsilon & - f ( t ) \sin \varepsilon + \cos \varepsilon \end{pmatrix}$$

The crossing instants are the zeros of the equation:

$$( 3. 1 4 ) & & 2 - 2 \cos \varepsilon + f ( t ) \sin \varepsilon = 0.$$

We let g ε ( ) t = 2 -2 cos ε + f ( ) sin t ε and we observe that g ε (0) = 2 -2 cos ε &gt; 0 for a positive sufficiently small ε . Moreover g ε ( τ ) = 0 iff

$$t _ { \varepsilon } \, = \, & \begin{cases} \frac { 2 ( 1 - \cos \varepsilon ) } { \sin \varepsilon } > 0 & \text{if $f(t)=-t$} \\ \frac { 2 ( \cos \varepsilon - 1 ) } { \sin \varepsilon } < 0 & \text{if $f(t)=t$} \end{cases} \\ \dots \, \frac { \varepsilon + \varepsilon + \varepsilon + \dots + \varepsilon + \varepsilon + \dots + \varepsilon + \varepsilon + \dots } \end{cases}.$$

↦

By this it directly follows that in the case f ( ) = t t , the (perturbed) path t → n ε ( ) t has no crossing instants in the interval [0 , T ] . 3 By this, immediately follows that

$$\iota ^ { C L M } ( \Delta, n _ { \varepsilon } ( \tau ), \tau \in [ 0, T ] ) = 0.$$

3 More generally if f ( ) t &gt; 0 for every t ∈ [ α, β ] ⊆ [0 , T ] , then the perturbed path has no crossing instants in [ α, β ] . In fact the equation:

$$f ( t _ { \varepsilon } ) = \frac { 2 ( 1 - \cos \varepsilon ) } { \sin \varepsilon }, \quad \forall \, t \in [ \alpha, \beta ] \\ \text{id to non-with-but-but-hand-oids-in-out}$$

has no solution as the left hand side is negative whilst the right hand side is positive.

↦

We let us now compute the ι CLM -index for f ( ) = t -t . We observe that we are in a very degenerate situation, since for every t ∈ [0 , T ] the path t → n ε ( ) t is contained in Sp (2) 0 , although it is not entirely contained in a fixed stratum, as the starting point is N (0) = I . Thus by the stratum homotopy invariance of the ι CLM -index, we immediately get that

$$\iota ^ { \text{CLM} } ( \Delta, n _ { \varepsilon } ( t ), t \in [ \delta, T ] ) = 0, \quad \forall \ \delta \in ( 0, T ]$$

and hence

$$( 3. 1 5 ) & \iota ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ \delta, T ] ) = \mathbf n _ { + } \, \Gamma \left ( n ( t ), \Delta, 0 \right ).$$

By Lemma 3.1, we get that

$$\eta ^ { \prime } ( t _ { \varepsilon } ) = - \frac { 1 } { d ( t _ { \varepsilon } ) } d ^ { \prime } ( t _ { \varepsilon } ) y _ { 0 }, \quad \xi ^ { \prime } ( t _ { \varepsilon } ) = b ^ { \prime } ( t _ { \varepsilon } ) y _ { 0 } - \frac { b ( t _ { \varepsilon } ) } { d ( t _ { \varepsilon } ) } d ^ { \prime } ( t _ { \varepsilon } ) y _ { 0 }.$$

Since v := ( x , y 0 0 ) T ∈ ker( N t glyph[epsilon1] ( ε ) -I ) , then we have x 0 = b t ( ε ) 1 -cos glyph[epsilon1] y 0 . By using Lemma 3.1 once more, we infer that

$$\Gamma ( n _ { \varepsilon }, \Delta, t _ { \varepsilon } ) ( v ) & = - x _ { 0 } \eta ^ { \prime } ( t _ { \varepsilon } ) - y _ { 0 } \xi ^ { \prime } ( t _ { \varepsilon } ) \\ & = \frac { b ( t _ { \varepsilon } ) } { d ( t _ { \varepsilon } ) } \frac { d ^ { \prime } ( t _ { \varepsilon } ) } { 1 - \cos \epsilon } y _ { 0 } ^ { 2 } - b ^ { \prime } ( t _ { \varepsilon } ) y _ { 0 } ^ { 2 } + \frac { b ( t _ { \varepsilon } ) } { d ( t _ { \varepsilon } ) } d ^ { \prime } ( t _ { \varepsilon } ) y _ { 0 } ^ { 2 } \\ & = \left [ \frac { b ( t _ { \varepsilon } ) } { d ( t _ { \varepsilon } ) } \left ( \frac { d ^ { \prime } ( t _ { \varepsilon } ) } { 1 - \cos \epsilon } + d ^ { \prime } ( t _ { \varepsilon } ) \right ) - b ^ { \prime } ( t _ { \varepsilon } ) \right ] y _ { 0 } ^ { 2 }. \\ \mu \quad \cdot \quad e \quad \nu ( t _ { \varepsilon } ) \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{cases}$$

The crossing form Γ( l, ∆ , t ε ) is a non-degenerate quadratic form on a one-dimensional vector space. In order to determine its inertia indexes, it is enough to know the sign of

$$C ( \Gamma ) \, = \frac { b ( t _ { \varepsilon } ) } { d ( t _ { \varepsilon } ) } \left ( \frac { d ^ { \prime } ( t _ { \varepsilon } ) } { 1 - \cos \varepsilon } + d ^ { \prime } ( t _ { \varepsilon } ) \right ) - b ^ { \prime } ( t _ { \varepsilon } ).$$

Since t ε = 2(1 -cos ε ) sin ε then we get b ′ ( t ε ) = -cos ε, d ′ ( t ε ) = sin ε, b ( t ε ) = -t ε cos ε +sin ε, d ( t ε ) = t ε sin ε +cos glyph[epsilon1].

Thus

$$( 3. 1 8 )$$

$$C ( \Gamma ) & = \frac { - t _ { \varepsilon } \cos \varepsilon + \sin \varepsilon } { t _ { \varepsilon } \sin \varepsilon + \cos \varepsilon } ( \frac { \sin \varepsilon } { 1 - \cos \varepsilon } + \sin \varepsilon ) + \cos \varepsilon \\ & = \frac { ( 1 - \cos \varepsilon ) ^ { 2 } } { 2 - \cos \varepsilon \sin \varepsilon } ( \frac { \sin \varepsilon } { 1 - \cos \varepsilon } + \sin \varepsilon ) + \cos \varepsilon \\ & = \frac { ( 1 - \cos \varepsilon ) \sin \varepsilon } { 2 - \cos \varepsilon \sin \varepsilon } + \frac { ( 1 - \cos \varepsilon ) ^ { 2 } \sin \varepsilon } { 2 - \cos \varepsilon \sin \varepsilon } + \cos \varepsilon. \\. \. \. \. \. \. \ \ ^ { \wedge \ + \quad. \ ( 1 - \cos \varepsilon ) \sin \varepsilon \quad. \ ( 1 - \cos \varepsilon ) ^ { 2 } \sin \varepsilon$$

It's easy to check that, for ε → 0 + , we get (1 -cos ε ) sin ε 2 -cos ε sin ε → 0 and (1 -cos ε ) 2 sin ε 2 -cos ε sin ε → 0 . Therefore C (Γ) ∼ 1 holds for sufficient small ε &gt; 0 . Hence the crossing form C (Γ) is positive definite. Thus by the previous computation and by Equation (3.15), we get

$$( 3. 1 9 )$$

$$\iota ^ { \text{CLM} } ( n _ { \epsilon }, \Delta,, t \in [ 0, T ] ) = 1.$$

This concludes the proof.

Remark 3.8 . A different proof of Lemma 3.7 could be as follows. By the localization axiom of the Maslov index [Gut14, Lemma 5.2 (Shear property) ], it follows that

$$\iota ^ { \text{RS} } ( n ( t ), \Delta, t \in [ 0, T ] ) = \begin{cases} - \frac { 1 } { 2 } & \text{if $f(t)=t$} \\ \frac { 1 } { 2 } & \text{if $f(t)=-t$}. \end{cases}$$

Now, by [LZ00, Equation (3.7) pag.97], it follows that

$$\iota ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ 0, T ] ) = \iota ^ { \text{RS} } ( n ( t ), \Delta, t \in [ 0, T ] ) + \frac { 1 } { 2 },$$

we can conclude that

$$\iota ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ 0, T ] ) = \begin{cases} 0 & \text{if $f(t)=t$} \\ 1 & \text{if $f(t)=-t$.} \end{cases}$$

The second claim follows by Equation (3.20) and Proposition 2.9. This concludes the proof.

↦

Remark 3.9 . A different proof of Lemma 3.7, without using any perturbation could be conceived even by using crossing forms. However, the reader should be aware on the fact that the case of symplectic shear is degenerate and a priori the formula for computing the ι CLM index through crossing forms, is not available in that form. A different proof of Lemma 3.7 goes as follows. We observe that we are in a very degenerate situation, since for every t ∈ [0 , T ] the path t → N t ( ) is contained in Sp (2) 0 even though it is not entirely contained in a fixed stratum, because of the starting point being in fact N (0) = I . Thus by the stratum homotopy invariant of the ι CLM -index, we immediately get that

$$\iota ^ { C L M } ( \Delta, n ( t ), t \in [ \varepsilon, T ] ) = 0, \quad \forall \ \varepsilon \in ( 0, T ]$$

and hence

$$\iota ^ { \text{CLM} } ( \Delta, n ( t ), t \in [ \varepsilon, T ] ) = \mathbf n _ { + } \, \Gamma ( n ( t ), \Delta, 0 ).$$

By a direct computation of the crossing form, we get

$$& a ( t ) = d ( t ) = 1, \quad b ( t ) = f ( t ), \quad c ( t ) = 0. \quad \text{Moreover}, \\ & a ^ { \prime } ( t ) = d ^ { \prime } ( t ) = c ^ { \prime } ( t ) = 0, \quad \quad b ^ { \prime } ( t ) = f ^ { \prime } ( t ).$$

Thus ξ ′ (0) = f ′ (0) y 0 and η ′ (0) = 0 . Using Equations (3.3) we get

$$\Gamma ( n, \Delta, 0 ) [ v ] = - x _ { 0 } \eta ^ { \prime } ( t _ { * } ) - y _ { 0 } \xi ^ { \prime } ( t _ { * } ) = - f ^ { \prime } ( 0 ) y _ { 0 } ^ { 2 }.$$

In particular for such a crossing form,

- · if f ( ) = t t , we get that n + Γ ( n t , ( ) ∆ 0 , ) = 0
- · if f ( ) = t -t , we get n + Γ ( n t , ( ) ∆ 0 , ) = 1

As the crossing form is degenerate, the conclusion follows by using the generalization through partial signatures as given by authors in [GPP04, Proposition 2.1] .

The second claim follows by the first one and Proposition 2.9. This concludes the proof.

By Lemma 3.7 and by the well-posedness of the ι CLM -index, the following is clear.

Corollary 3.10. Let N : [0 , T ] → Sp(2) be the symplectic path defined in Lemma 3.7 and N ε pointwise defined by N ε ( ) t = e -εJ N t ( ) . Thus we have:

$$\iota ^ { \text{CZ} } ( N _ { \varepsilon } ( t ), t \in [ 0, T ] ) = \begin{cases} - 1 & \text{if $f(t)=t$} \\ 0 & \text{if $f(t)=-t$.} \end{cases}$$

## 4 Indices and stability of Keplerian orbits

The aim of this section is to compute the Conley-Zehnder index of a Keplerian ellipses with eccentricity e ∈ [0 , 1) . In Subsection 4.1 we explicitly compute the ι CZ -index in the case of circular Keplerian orbit and finally we conclude the general case of Keplerian ellipses having eccentricity e ∈ [0 , 1) .

## 4.1 The Conley-Zehnder index in polar coordinates

Let us now come back to the Lagrangian function given at Equation 1.7 whose induced Hamiltonian function is given by

$$H ( r, \vartheta, \dot { r }, \dot { \vartheta } ) = \frac { 1 } { 2 } \mu ( \dot { r } ^ { 2 } + r ^ { 2 } \dot { \vartheta } ^ { 2 } ) - U ( r ) = \frac { 1 } { 2 } \left [ \frac { p _ { r } ^ { 2 } } { \mu } + \frac { p _ { \vartheta } ^ { 2 } } { r ^ { 2 } \mu } \right ] - \frac { m } { r }$$

where ( p , p r ϑ ) = ( µr, µr ˙ 2 ˙ ϑ ) . The induced Hamiltonian system is given by

$$( 4. 2 ) & \quad \begin{cases} \dot { p } _ { r } = - \partial _ { r } H ( p _ { r }, p _ { \emptyset }, r, \emptyset ) = - \frac { m } { r ^ { 2 } } + \frac { p _ { \emptyset } ^ { 2 } } { \mu r ^ { 3 } } \\ \dot { p } _ { \emptyset } = - \partial _ { \emptyset } H ( p _ { r }, p _ { \emptyset }, r, \emptyset ) = 0 \\ \dot { r } = \partial _ { p _ { r } } H ( p _ { r }, p _ { \emptyset }, r, \emptyset ) = \frac { p _ { r } } { r ^ { 2 } \mu } \\ \dot { \emptyset } = \partial _ { p _ { \emptyset } } H ( p _ { r }, p _ { \emptyset }, r, \emptyset ) = \frac { p _ { \emptyset } } { r ^ { 2 } \mu } \end{cases}$$

By linearizing the Hamiltonian system given in Equation (4.2) at the circular solution r t ( ) = re iωt for t ∈ [0 , 2 π/ω ] , then we get

$$\text{on } \iota \in \{ \mathfrak { v }, \varepsilon \, \overline { n } / \omega \}, \, \text{even we gev} \\ ( 4. 3 ) & \quad \begin{cases} \dot { y } _ { r } = \left [ \frac { 2 m } { r ^ { 3 } } - \frac { 3 p _ { \vartheta } ^ { 2 } } { \mu r ^ { 4 } } \right ] x _ { r } + \frac { 2 p _ { \vartheta } } { \mu r ^ { 3 } } y _ { \vartheta } \\ \dot { y } _ { \vartheta } = 0 \\ \dot { x } _ { r } = \frac { y _ { r } } { \mu } \\ \dot { x } _ { \vartheta } = - \frac { 2 p _ { \vartheta } } { r ^ { 3 } \mu } x _ { r } + \frac { y _ { \vartheta } } { r ^ { 2 } \mu }. \\ \varepsilon \, \cdots \, \cdots \, \quad \, \int \, \mathfrak { T } \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdot \end{cases}$$

Setting w = ( y , y r ϑ , x r , x ϑ ) T , then the linearized Hamiltonian system given at Equation (4.3) can be written as ˙ w = Lw where L is the four by four matrix given by

$$L \, = \, & \begin{bmatrix} 0 & \frac { 2 p _ { \vartheta } } { \mu r ^ { 3 } } & \left [ \frac { 2 m } { r ^ { 3 } } - \frac { 3 p _ { \vartheta } ^ { 2 } } { \mu r ^ { 4 } } \right ] & 0 \\ 0 & 0 & 0 & 0 \\ \frac { 1 } { \mu } & 0 & 0 & 0 \\ 0 & \frac { 1 } { r ^ { 2 } \mu } & - \frac { 2 p _ { \vartheta } } { r ^ { 3 } \mu } & 0 \end{bmatrix} \, = \, & \begin{bmatrix} 0 & A & C & 0 \\ 0 & 0 & 0 & 0 \\ D & 0 & 0 & 0 \\ 0 & B & - A & 0 \end{bmatrix} \\ \text{where}$$

where

$$A \, = \frac { 2 p _ { \vartheta } } { \mu r ^ { 3 } }, \ \ B \, = \frac { 1 } { r ^ { 2 } \mu }, \ \ C \, = \left [ \frac { 2 m } { r ^ { 3 } } - \frac { 3 p _ { \vartheta } ^ { 2 } } { \mu r ^ { 4 } } \right ] \text{ and finally } \ D \, = \frac { 1 } { \mu }.$$

It is worthwhile to observe that, the matrix L is a time independent Hamiltonian matrix. Thus, by setting, -CD = ω 2 , the fundamental (matrix) solution is given by φ 0 ( ) = t e Lt where t ∈ [ 0 , 2 π ω ] . Now, since the determinant of L -λI 4 is λ 2 ( λ 2 + ω 2 ) , there exists a symplectic matrix P ∈ Sp(4)

$$L = P \left ( \begin{bmatrix} 0 & s ( r ) \\ 0 & 0 \end{bmatrix} \diamond \begin{bmatrix} 0 & - \omega \\ \omega & 0 \end{bmatrix} \right ) P ^ { - 1 }$$

such that

↦

where r → s r ( ) is a positive function (to be determined) and by using the direct sum property of the Conley-Zehnder index (cf. Section 2), then we get that

$$\iota ^ { \text{cz} } \Big ( \phi _ { 0 } ( t ), t \in [ 0, T ] \Big ) = \iota ^ { \text{cz} } \Big ( \overline { \phi } _ { 0 } ( t ), t \in [ 0, T ] \Big ) + 1$$

where φ t ( ) = [ 1 s r t ( ) 0 1 ] , and t ∈ [0 , T ] and T = 2 π/ω . In order to compute the function s r ( ) (actually, by using Lemma 3.7, we only need to compute the sign of this function), we proceed

as follows. Denoting by E = { e , e 1 2 , e 3 , e 4 } the canonical basis of R 4 we start to observe that Le 4 = 0 . Moreover, setting v = (0 , -1 , A/C, 0) T , then we get that

$$L v = \left ( \frac { A ^ { 2 } } { C } - B \right ) e _ { 4 } = \left ( \frac { A ^ { 2 } D } { \omega ^ { 2 } } - B \right ) e _ { 4 } = - \frac { 1 } { \mu r ^ { 2 } } \left [ 1 + \frac { 4 p _ { \vartheta } ^ { 2 } } { 2 m \mu r - 3 p _ { \vartheta } ^ { 2 } } \right ] e _ { 4 }.$$

Since for circular motions, by the discussion performed at the end of Section 1 and more precisely at Equation (1.19), we get that r = r 0 = k / µm 2 ( ) = p 2 ϑ / µm ( ) were k denotes the angular momentum, we finally get that

$$L v = \frac { 3 } { \mu r ^ { 2 } } e _ { 4 }.$$

↦

We observe that, being ω 0 ( e , v 4 ) = 1 , it readily follows that { e , v 4 } is a symplectic basis of that invariant subspace and the function r → s r ( ) = 3 µr 2 is positive. In particular, by using Lemma 3.7, we get that

$$\iota ^ { \text{cz} } ( \overline { \phi } _ { 0 } ( t ), t \in [ 0, T ] ) = - 1.$$

Summing up the previous computation we finally get the following result.

Proposition 4.1. The Conley-Zehnder index of a planar circular solution of Equation (1.5) on a prime period, vanishes.

Remark 4.2 . In the case of homogeneous singular potentials an analogous result was already proved by authors in [BJP16]. A similar result in the case of Keplerian orbits on constant curvature surfaces has been proved in [DDZ19].

Proposition 4.3. Let φ e be the monodromy matrix of the Keplerian ellipse having period T and eccentricity e . Then there exists P ∈ Sp(4) such that

$$M = P ^ { - 1 } \left ( N _ { 1 } ( 1, 1 ) \diamond I _ { 2 } \right ) P,$$

$$\ w h e r e \ N _ { 1 } ( 1, 1 ) \, = \, \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}.$$

Proof. By [HS10, Lemma 3.3] we get that there exists P ∈ Sp(4) such that

$$M = P ^ { - 1 } \left ( N _ { 1 } ( 1, 1 ) \diamond M _ { 1 } \right ) P$$

where the 2 × 2 Jordan block N 1 (1 , 1) corresponds to the energy integral. By the integrability of the Kepler problem, we know that the angular momentum is a first integral. This implies that σ M ( 1 ) = { 1 } . By the basic normal form of a 2 × 2 symplectic matrix, M 1 has to be symplectically similar to a matrix of the form

[

]

N

1

b

1

(1

, b

)

=

0

1

where

b

∈ {

0 1

,

,

- }

1

.

Now, we observe that in the negative energy 3 -dimensional hypersurface of the phase space, every solution is an elliptic orbit having prime period T . So the time T fundamental solution restricted to the fixed negative energy hypersurface Σ h corresponding to the energy level h (which is a φ e ( ) t -invariant manifold) has to be the identity map (otherwise Σ h would not be invariant). By this argument, we directly conclude that

$$\dim \ker ( M - I _ { 4 } ) = 3$$

and so M is symplectically similar to

$$N _ { 1 } ( 1, 1 ) \diamond I _ { 2 }.$$

This concludes the proof.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00075v2-figure-9631f1f6-a242-42ee-badf-59b59ef63aa7.jpg)

Theorem 4.4. Let γ be a Keplerian ellipse (i.e. a solution of Equation (1.5) with prime period T ) and let γ k its k -th iteration. Then, we have

$$\iota ( \gamma ^ { k } ) = 2 ( k - 1 ) \quad \forall \, k \in \mathbb { N } ^ { * }$$

where, as before, we denoted by ι ( γ ) the Morse index of γ . In particular

$$\iota ( \gamma ) = 0 \quad a n d \ \lim _ { k \to + \infty } \frac { \iota ( \gamma ^ { k } ) } { k } = 2.$$

Proof. Since the time rescaling doesn't change the Morse index, by Proposition 2.10, we get

$$\iota ( \gamma ^ { k } ) = \iota ^ { \text{cz} } ( \phi _ { e } ( t ), t \in [ 0, T k ] ).$$

By invoking Equation (4.5) and as direct consequence of Lemma 3.7 and Lemma 3.5 and by using the additivity property of the Conley-Zehnder index under concatenation of paths, we finally get that

$$\iota ^ { C Z } ( \phi _ { e } ( t ), t \in [ 0, T k ] ) = 2 k - 2$$

To conclude, we observe that since the Conley-Zehnder index of the fundamental solution φ e only depends on the monodromy matrix, the thesis follows by previous computation and by Equation (4.5). This concludes the proof of the first part. The second follows straightforward.

Remark 4.5 . It is worth noticing that the contribution given by the conservation law of the energy to the ι CZ -index is -1 and the symplectic normal form is given by the 2 × 2 Jordan block relative to the eigenvalue 1 whereas the symplectic normal form corresponding to the conservation law of the angular momentum is given by the 2 × 2 identity matrix. As expected we are in a very degenerate situation.

We conclude the section by summarizing the stability properties of the Keplerian ellipses.

Theorem 4.6. Let γ be a Keplerian ellipse. Then it is elliptic, meaning that all the eigenvalues belongs to U . Moreover it is spectrally stable and not linearly stable.

Proof. Denoting by M the monodromy matrix, it follows that there exists P ∈ Sp(4) such that

$$M & = P ^ { - 1 } \left ( \begin{bmatrix} 1 & 1 \\ 0 & 1 \end{bmatrix} \diamond \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} \right ) P.$$

In particular σ M ( ) = { 1 } ∈ U and the algebraic multiplicity of its (unique) Floquet multiplier is 4 . This proves the first claim. The second follows by observing that M is not diagonalizable (having a non-trivial Jordan block); thus in particular γ is spectrally and not linearly stable. We also observe that its nullity is 3 = dim ker( M -I ) .

## References

[Arn86] Arnol'd, V. I. Sturm theorems and symplectic geometry. (Russian) Funktsional. Anal. i Prilozhen. 19 (1985), no. 4, 1-10, 95.

[BJP14] Barutello, Vivina L.; Jadanza, Riccardo D.; Portaluri, Alessandro Linear instability of relative equilibria for n-body problems in the plane. J. Differential Equations 257 (2014), no. 6, 1773-1813.

[BJP16] Barutello, Vivina; Jadanza, Riccardo D.; Portaluri, Alessandro Morse index and linear stability of the Lagrangian circular orbit in a three-body-type problem via index theory. Arch. Ration. Mech. Anal. 219 (2016), no. 1, 387-444.

[CLM03] Cappell, Sylvain E.; Lee, Ronnie; Miller, Edward Y. On the Maslov index. Comm. Pure Appl. Math. 47 (1994), no. 2, 121-186.

[DDZ19] Deng, Yanxia; Diacu, Florin; Zhu, Shuqiang Variational property of Keplerian orbits by Maslov-type index To appear in J. Differential Equations

[Dui76] Duistermaat, J. J. On the Morse index in Variational Calculus. Adv. Math. 21 (1976), 173-195.

[GPP04] Giambò, Roberto; Piccione, Paolo; Portaluri, Alessandro Computation of the Maslov index and the spectral flow via partial signatures. C. R. Math. Acad. Sci. Paris 338 (2004), no. 5, 397-402.

[Gor77] Gordon, William B. A minimizing property of Keplerian orbits. Amer. J. Math. 99 (1977), no. 5, 961-971.

[Gut14] Gutt, Jean Normal forms for symplectic matrices. Port. Math. 71 (2014), no. 2, 109139.

[HLS14] Hu, Xijun; Long, Yiming; Sun, Shanzhong Linear stability of elliptic Lagrangian solutions of the planar three-body problem via index theory. Arch. Ration. Mech. Anal. 213 (2014), no. 3, 993-1045.

[HP17] Hu, Xijun; Portaluri, Alessandro Index theory for heteroclinic orbits of Hamiltonian systems. Calc. Var. Partial differential Equations 56 (2017), no. 6, Art. 167, 24 pp.

[HS10] Hu, Xijun; Sun, Shanzhong Morse index and stability of elliptic Lagrangian solutions in the planar three-body problem. Adv. Math. 223 (2010), no. 1, 98-119.

[HS09] Hu, Xijun; Sun, Shanzhong Index and stability of symmetric periodic orbits in Hamiltonian systems with application to figure-eight orbit. Comm. Math. Phys. 290 (2009), no. 2, 737-777.

[Lon02] Long, Yiming Index theory for symplectic paths with applications. Birkhäuser Verlag, Basel, 2002.

[LZ00] Long, Yiming; Zhu, Chaofeng Maslov-type index theory for symplectic paths and spectral flow. II. Chinese Ann. Math. Ser. B 21 (2000), no. 1, 89-108.

[MS05] Meyer, Kenneth R.; Schmidt, Dieter S. Elliptic relative equilibria in the N-body problem. J. Differential Equations 214 (2005), no. 2, 256-298.

[RS93] Robbin, Joel; Salamon, Dietmar The Maslov index for paths. Topology 32 (1993), no. 4, 827-844.

Henry Kavle Department of Mathematics and Statistics Queen's University, Kingston (ON) K7K 3N6 Kingston (Ontario) Canada E-mail: kavle.h@queensu.ca

Prof. Daniel Offin Department of Mathematics and Statistics Queen's University, Kingston (ON)

K7K 3N6 Kingston (Ontario) Canada E-mail: offind@queensu.ca

Prof. Alessandro Portaluri

DISAFA

Università degli Studi di Torino

Largo Paolo Braccini, 2

10095 Grugliasco, Torino

Italy

Website: https://sites.google.com/view/alessandro-portaluri/home

E-mail:

alessandro.portaluri@unito.it
Henry Kavle Daniel Offin Alessandro Portaluri
Keplerian orbits through the Conley-Zehnder index
It was discovered by Gordon in 1977 that Keplerian ellipses in the plane are
minimizers of the Lagrangian action and spectrally stable as periodic points of
the associated Hamiltonian flow. The aim of this paper is to give a homotopy
theoretical proof of these results through a self-contained, explicit and
simple computation of the Conley-Zehnder index.
  The techniques developed in this paper can be used to investigate the higher
dimensional case of Keplerian ellipses, where the classical variational proof
no longer applies.