The dynamics of the relativistic Kepler problem
Elbaz I. Abouelmagda, Juan Luis Garc´ıaGuiraob, Jaume Llibrec
aCelestial Mechanics and Space Dynamics Research Group (CMSDRG), Astronomy Department, National Research Institute of Astronomy and Geophysics (NRIAG), Helwan – 11421, Cairo, Egypt. bDepartamento de Matem´atica Aplicada y Estad´ıstica. Universidad Polit´ecnica de Cartagena, Hospital de Marina, 30203-Cartagena, Regi´onde Murcia, Spain. cDepartament de Matem`atiques.Universitat Aut`onomade Barcelona, Bellaterra, 08193-Barcelona, Catalonia, Spain.
Abstract We deal with the Hamiltonian system (HS) associated to the Hamiltonian in polar coordinates ! 1 p2 1 H = p2 + φ − − , 2 r r2 r 2r2 where is a small parameter. This Hamiltonain comes from the correction given by the special relativity to the motion of the two–body problem, or by the first order correction to the two–body problem coming from the general relativity. This Hamiltonian system is completely integrable with the angular momentum C and the Hamiltonian H. We have two objectives. First we describe the global dynamics of the Hamiltonian system (HS) in the following sense. Let Sh and Sc are the subset of the phase space where H = h and C = c, respectively. Since C and H are first integrals, the sets Sc, Sh and Shc = Sh ∩ Sc are invariant by the action of the flow of the Hamiltonian system (HS). We determine the global dynamics on those sets when the values of h and c vary. Second recently Tudoran in [33] provided a criterion which detects when a non–degenerate equilibrium point of a completely integrable system is Lyapunov stable. Every equilibrium point q of the completely integrable Hamiltonian system (HS) is degenerate and has zero angular momentum, so the mentioned criterion cannot be applied to it. But we will show that this criterion is also satisfied when it is applied to the Hamiltonian system (HS) restricted to zero angular momentum. Keywords: Kepler problem, Global dynamics, Lyapunov, Stability.
1. Introduction The Kepler problem appears in various fields, some of these are beyond the physics, which have been studied In celestial mechanics, the Kepler problem is a special by Kepler himself. This problem is important in celestial dynamical system coming from the two–body problem. In mechanics because the Newtonian gravity obeys the law this model, two objects move under their mutual Newto- of inverse square distance between two bodies. Thus for nian gravitational force [1]. This attractive force varies in example, the motion of two stars around each other, the size as the inverse square of the separation distance be- planets moving surrounding the Sun, a satellite orbiting a tween them and it is proportional to the product of the planet, and many other examples of orbital motion. The masses of the two bodies. This system is used to find Kepler problem has also a significant relevance in the study the position and the velocity vectors of the two bodies at of the motion of two charged particles, because Coulomb’s specified time. Using the laws of classical mechanics, the law of electrostatics follows the law of inverse square dis- solution can be provided as a Kepler orbit through finding tance too. For Example the hydrogen atom, muonium and the six orbital elements. This problem is called the Kepler positronium, which have played serious roles as models of problem in honor of the German astronomer Johannes Ke- dynamical systems to test physical theories and measuring pler, after he proposed Kepler’s laws of the motion of the constants of nature, see [3–5]. planets and illustrated the kinds of forces which can pro- The importance of the Kepler problem is not only re- vide orbits obeying those laws [2]. lated to its varied applications in different fields of sci- ence, but rather to its use in developing some mathemati- Email addresses: [email protected] or cal methods. Thus, this problem has been used to develop [email protected] (Elbaz I. Abouelmagd), new methods in classical mechanics, like Hamiltonian me- [email protected] (Juan Luis Garc´ıaGuirao), [email protected] (Jaume Llibre) chanics, the Hamilton–Jacobi equation, Lagrangian me-
Preprint submitted to Results in Physics August 23, 2020 chanics and action-angle coordinates. Furthermore, the observation. Kepler and simple harmonic oscillator problems are two There were several pre– and post–relativistic attempts of the most fundamental problems in classical mechanics. to get such modified models. This is the case of the Manev They are also the only integrable dynamical systems which Hamiltonian introduced by Manev in [15–18]. See for more have closed orbits for open sets of possible initial condi- information on the two–body problem, for instance, [19– tions 23]. The Kepler problem is an unperturbed version form of The correction given by the special relativity to the the two–body problem, which must be extended or gener- motion of the two–body problem, or by the first order cor- alized to include to additional forces as perturbations, that rection to this problem coming from the general relativity represent realistic natural phenomena or to obtain precise is ! and accuracy results. Within the frame of anisotropic per- 1 p2 1 H = p2 + φ − − , (1) turbation in [6], the authors have showed that the anisotropic 2 r r2 r 2r2 Kepler problem is equivalent to two massless particles which move in a plane or on perpendicular lines and which in- where || 1 is a small parameter, for details see [24– teract according to the Newton’s law of gravitation. They 28]. Note that, when = 0, we have the rotating Kepler have also proved that the generalized issue of anisotropic problem, for more information on this last problem see Kepler problem and anisotropic two centres problem are [29]. non–integrable systems. For an analysis of the Kepler Note that the Hamiltonian H is a particular case of the problem, specific to radial periodic perturbation of a cen- Manev Hamiltonian H when a = −1 and b = −/2. tral force field, it has been proved the existence of rotating The Hamiltonian system defined by the Hamiltonian periodic solutions nearly circular orbits [7]. Implicit Func- H is r˙ = p , tion Theorem is used to find such solutions. r p The solution of the Kepler problem allowed researchers φ˙ = φ , to investigate the unperturbed and the perturbed Lagrange’s r2 (2) planetary motion. Thus, these motions were clarified en- p2 − 1 p˙ = φ − , tirely by the classical mechanics and the Newton’s law of r r3 r2 universal gravitation [8, 9]. More recently, the perturbed p˙ = 0. Kepler problem have been investigated by [10, 11]. The φ perturbed models are not limited to Kepler models but For a fixed < 0 small, this Hamiltonian system has the also include the restricted three–body problem, which can circle of equilibria be reduced to the perturbed Kepler problem in some cases 1 [12–14]. {q(φ) = (r, φ, pφ, pr) = (−, φ, 0, 0) : φ ∈ S }. The Manev Hamiltonian is The Hamiltonian dynamical system (2) is perfectly in- ! 1 p2 a b tegrable within frame the of Liouville–Arnold because it H = p2 + φ + + , 2 r r2 r r2 has two independent first integrals the Hamiltonian H and the angular momentum pφ in involution. For more details where a and b are arbitrary constants. This Hamiltonian on perfectly integrable Hamiltonian systems see [30, 31]. describes the motion of a two–body problem, which is gov- The objective of this paper is double. Our first objec- erned by the potential a/r+b/r2, where r is the separation tive is to describe the global dynamics of the Hamiltonian system (2) in the following sense. Let Sh, Sc are the subset distance between the two bodies and (pr, pφ) are the mo- menta in polar coordinates. of the phase space, H = h and pφ = c where h and c are The analysis of the motion in the two–body problems the constants of integration. Due to the fact that H and pφ has a long history. Just after the Newton’s work on the are first integrals, the sets Sh, Sc and Shc = Sh ∩Sc are in- two–body gravitational problem, some discrepancies ap- variant under the flow of the Hamiltonian system (2). We peared between the theoretical motions of the pericenters determine the global dynamics on those sets when h and of the planets and the observed ones. Consequently, some c vary. Moreover, we describe the foliation of the phase doubts on the accuracy of the inverse square law of Newton space by the invariant sets Sh, and the foliation of Sh by for gravitation. Motivated scientists to construct alterna- the invariant sets Shc. See section 2. tive gravitational models and corrections trying to recon- In fact, this first objective is a particular case of the cile these discrepancies. In fact, Newton was the first to general cases studied inside the Manev Hamiltonian H consider what we now call the Manev systems, see the with arbitrary values of a and b, see [32]. However, from book I, section IX, proposition XLIV, theorem XIV and that work, it is not easy to obtain the global dynamics for corollary 2 of the Principia. the particular case a = −1 and b = −/2, corresponding Several authors tried to find or construct an appro- to our Hamiltonian (1). We cover this case in this work. priate model with the features of the Newtonian one but Second, we shall show that every equilibrium point q(φ) with the convenient corrections, to make closer theory and of the Hamiltonian system (2) restricted to zero angular 2 momentum is Lyapunov stable. However, our interest in If h = 0 then Rh is homeomorphic to this Lyapunov stability rises in showing that the criterion h provided recently by Tudoran in [33] which detects when − , ∞ × 1 if < 0, 2 S (4) a non–degenerate equilibrium point of a completely inte- + 1 R × S if ≥ 0. grable system is Lyapunov stable, also can be extended to the degenerate equilibrium points of the completely in- If h > 0 then Rh is homeomorphic to tegrable system (2) restricted to zero angular momentum, √ where the equilibrium points of system (2) live, see section 1 − 2h − 1 1 , ∞ × S if < 0, 3. h (5) R+ × S1 if ≥ 0. 2. On the global dynamics 2.3. The sets S 2.1. Critical points and critical values h We will consider the notation and terminology of [32], Now we determine the topology of the invariant energy levels Sh by using the fact: as consequence, we suppose that Sh and Sc are the subset of the phase space where H = h and pφ. Science H and Sh = {(r, φ, pr, pφ) ∈ ℵ : g(r, pr, pφ) = h} pφ are first integrals, the sets Sh, Sc and Shc = Sh ∩Sc are −1 1 invariant under the flow of the Hamiltonian system (2), ≈ {g (h)} × S , i.e. if an orbit of the Hamiltonian system has a point on where the set Shc then the whole orbit is contained in this set. + 1 2 ! Let r ∈ R = (0, ∞), φ ∈ S ,(pr, pφ) ∈ R and ℵ = 1 p2 1 + 1 2 2 φ × × , where g(r, pr, pφ) = p + − − , R S R 2 r r2 r 2r2 Sh = {(r, φ, pr, pφ) ∈ ℵ : H(r, φ, pr, pφ) = h}, If h is a regular momentum value of the map g : R+ × Sc = {(r, φ, pr, pφ) ∈ ℵ : pφ = c}. R2 −→ R and the set of points {g−1(h)}= 6 ∅, then the set In this case the set of critical points of H is {g−1(h)} is a surface in R+ × R2. Hence the intersection −1 C = {(r, φ, p , p ) ∈ ℵ : r + = 0, φ ∈ 1} of {g (h)} with {r = r0 = constant}, is r φ S 1 either an ellipse if + + h > 0, since r > 0, then 2 r0 2r0 ( ∅ if ≥ 0, ) 1 C = or a point if + 2 + h = 0, − if < 0. r0 2r0 1 hence the critical value of H is 1/2 if < 0 or the empty set if + 2 + h < 0. r0 2r0 2.2. Hill Regions From the definition of the energy levels Sh we obtain + 1 Let ℵ and R = R × S be the phase space and the [ configuration space of the Hamiltonian system (2), and let Sh = ℵ(r,φ) (6) Γ: ℵ → R be the projection from ℵ to R. Then for each (r,φ)∈Rh h belongs to real set R, the regions of motion Rh (Hill where regions) of Sh are defined by Γ(Sh) = Rh, for more details p2 see [32, 34], hence 2 φ 1 2 ℵ(r,φ) = {(r, φ, pr, pφ) ∈ ℵ : pr+ = 2hr + 2r + }. 1 r2 4r2 R = {(r, φ) ∈ R : − − ≤ h} h r 2r2 For every point (r, φ) ∈ Rh the set ℵ(r,φ) is an ellipse, a = {r ∈ R+ : 2hr2 + 2r + ≥ 0} × S1. point or the empty set if the point (r, φ) belongs to the interior of the Hill region Rh, to the boundary of Rh, or Therefore, if h < 0 then Rh is homeomorphic to does not belong to Rh, respectively. Hence every energy 1 ∅ if < , level Sh of the planar relativistic Kepler problem is home- 2h omorphic, from (3) and (6), to 1 1 − × S1 if = , 1 h 2h ∅ if < , 2h √ √ 1 − 1 − 2h 1 + 1 − 2h 1 1 1 − , − × S if < < 0, 1 if = , h h 2h S 2h √ 1 + 1 − 2h 1 0, − × 1 if ≥ 0. S2 × S1 if < < 0, h S 2h (3) (S2 × S1) \ S1 if 0 ≤ , 3 if h < 0; and from (4), (5) and (6) to The manifold Sh is homeomorphic to a solid torus with- out its boundary and without its central circular axis. 2 1 1 (S × S ) \ S if < 0, Hence we can find the dynamics on Sh by rotating Fig. + 1 1 2(b) around the e axis. From this figure we obtain one R × S × S if ≥ 0, cylinder Shc for every c ∈ R \{0}, and two cylinders Shc if h ≥ 0. for every c = 0, which foliate Sh. Case 3: h ≥ 0 and < 0. Now g−1(h) is homeomor- 2 2.4. The sets Sc phic to a plane R see Fig. 3(a), and the curves γhc are homeomorphic to . Since Sc = {(r, φ, pr, pφ) ∈ ℵ : pφ = c}, we get that Sc R + 1 The manifold S is homeomorphic to a solid torus with- is homeomorphic to R × S × R for all c ∈ R. h out its boundary. The dynamics of Sh can be obtained by rotating Fig. 3(b) around the e axis. From this figure 2.5. The foliation of Sh by Shc the foliation of Sh is done by the cylinders Shc for every We can evaluate the invariant set Shc from knowing −1 c ∈ R \{0}. Figures (1, 2, 3) already have been appeared the set {g (h)} and in [32] in a more general context.
Shc = Sh ∩ Sc 3. On the Lyapunov stability of the equilibria = Sh ∩ {pφ = c}