Next Article in Journal
The Effect of the LMC on the Milky Way System
Previous Article in Journal / Special Issue
The Simplest Parametrization of the Equation of State Parameter in the Scalar Field Universe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Periodic Orbits of the Perturbed Two- and Three-Body Problems

by
Elbaz I. Abouelmagd
1,*,
Juan Luis García Guirao
2 and
Jaume Llibre
3
1
Celestial Mechanics and Space Dynamics Research Group (CMSDRG), Astronomy Department, National Research Institute of Astronomy and Geophysics (NRIAG), Helwan 11421, Cairo, Egypt
2
Departamento de Matemática Aplicada y Estadística, Universidad Politécnica de Cartagena, Hospital de Marina, 30203 Cartagena, Spain
3
Departament de Matemàtiques, Universitat Autònoma de Barcelona, Bellaterra, 08193 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Galaxies 2023, 11(2), 58; https://doi.org/10.3390/galaxies11020058
Submission received: 19 March 2023 / Revised: 8 April 2023 / Accepted: 14 April 2023 / Published: 18 April 2023

Abstract

:
In this work, a perturbed system of the restricted three-body problem is derived when the perturbation forces are conservative alongside the corresponding mean motion of two primaries bodies. Thus, we have proved that the first and second types of periodic orbits of the rotating Kepler problem can persist for all perturbed two-body and circular restricted three-body problems when the perturbation forces are conservative or the perturbed motion has its own extended Jacobian integral.

1. Introduction

In celestial mechanics, the periodic orbits have a considerable importance due to the existence of a direct relationship among the periodic orbits and the motion of planetary systems, as well as the motion of most stellar systems. Thus, in the literature of celestial mechanics, there are many studies on periodic orbits. In 1897, Henri Poincaré studied periodic orbits, and he considered that the exploration of such orbits is important for the understanding of the dynamics of differential systems. Through his study on the restricted three-body problem, he distinguished three classes of periodic orbits [1]:
  • First, those emerging from the circular periodic orbits of the rotating Kepler problem;
  • Second, those coming from the elliptical periodic orbits of the rotating Kepler problem;
  • Third, those generated from the spatial rotating Kepler problem, where the plane of motion is inclined.
The analytical studies on these three types of periodic orbits are developed in [2,3].
In space science, the importance of periodic orbits is not only restricted to celestial mechanics but also to astrodynamics. There is a need to find periodic orbits to design different types of space missions. For example, in [4], the presence of periodic orbits of the first type is confirmed using the small-parameter technique when the period of a spacecraft’s motion and the undisturbed circular orbit of the primaries are equal. Furthermore, the obtained results can be applied to the design and ballistic analysis of electrically powered thrusting spacecraft. Additionally, polar periodic orbits with π / 2 inclination are employed to observe planet surfaces subjected to space missions [5,6].
Periodic solutions or periodic orbits play outstanding roles not only in physical, mathematical, engineering systems but also in biological biosphere systems, where a paramount problem is to estimate or explore if the automatic oscillatory activity can be continued when it is subjected to a small external effect. Many researchers have studied the perseverance and continuation problems of the periodic orbits. Although these studies are carried out from several standpoints all of them obtained the same basic outcome under an assumption regarding the period manifold, known as normal non-degeneracy, where this assumption is expressible as a condition on the unperturbed flow over the manifold [7,8,9,10,11,12,13].
Using Delaunay variables, the periodic orbits of the second type were investigated in [14]. The presence of second-type period orbits was developed under the perturbation secular potential for the perturbed three-body problem [15]. The existence of first-type orbits in the framework of the restricted photo-gravitational four-body problem was studied in [16]. Recently, in [17], the authors studied the resonance transition structures using the continuation method. Then, they described how a 1 : 1 resonance is progressively overlapped in the proximity of the first-order resonance through the increases in the perturbation. They also proposed a technique to find these orbits as well as the transition resonance with respect to the system of the Sun and Jupiter. Further considerable studies on periodic orbits that would serve as an excellent guide to readers were performed in [18,19,20,21,22,23,24,25,26,27]
The innovative part of this study is the construction of the general perturbed equations of an infinitesimal body under many different perturbations, which come from conservative forces, such as a lack of sphericity and the radiation pressure or asteroid belt effect. Thus, we will prove that the circular and elliptic periodic orbits of the planar rotating Kepler problem can be prolonged to the periodic orbits of first and second type for the perturbed two-body problem and the perturbed planar circular restricted three-body problem when the sources of perturbations are conservative forces.

2. Perturbation and Mean Motion

2.1. Perturbed Forces

In dynamical astronomy or in stellar and solar system motion, perturbation is the complex motion of a celestial body that is moving not only under the effect of the main force of gravitation but is also subjected to other forces. These forces can involve the gravitational attraction of another body; two, three, or more bodies instead of the resistance of drag forces; the off-center oblate (or triaxial and misshapen body) attraction; or the effect from asteroid belts and Kirkwood gaps.
The simple motion of a body under a mutual gravitation effect of another body can be described in geometrical terms as conic sections, which are two unperturbed bodies in Kepler motion, but this description originates only from a hypothetical or theoretical point of view. The difference between the theoretical and real motion of celestial bodies considers the effect of perturbation forces or the additional gravitational forces of the surrounding and enclosure bodies.
If there is a third body (or multiple bodies) affecting the two-body motion, then the motion becomes perturbed, and it can be considered a three-body (or N-body) motion. However, the analytical solution with a closed form or a specified mathematical expression to identify the locations and velocities of the body at any time only exists for Kepler motion. Thus, the two-body problem becomes unsolvable if at least one of the two bodies has an irregular shape.
In the motion of most dynamical systems, which include many gravitational attractions (or other perturbed forces), there is one main body, which is predominant in its effect. The other gravitational effects of the remaining bodies (or other perturbed forces) can be treated as a perturbation for the hypothetical unperturbed motion of the planet (or satellite) around its star or planet. In this context, we will construct the perturbed mean motion of a two-body system in the next subsection.

2.2. Perturbed Mean Motion

We assume that the perturbed potential between two masses m i and m j can be written in the following form:
Γ i j = G m i m j 1 ρ i j + ϝ i j ( ρ i j , ε i j ) ,
where ρ i j is the distance between the two centers of the masses m i and m j . ϝ i j ( ρ i j , ε i j ) is the perturbation force such that it tends to zero (infinity) when ρ i j tends to infinity (zero); it can be also written as ϝ i j ( ρ i j , ε i j ) = ε i j Υ i j ( ρ i j ) , where Υ i j ( ρ i j ) is differentiable n −times with respect to ρ i j . In addition, ε i j is a small parameter, and its value depends on the kind of perturbation force.
Let ρ 1 and ρ 2 be the position vectors of masses m 1 and m 2 , respectively, in the inertial frame, while ρ 12 = ρ 2 ρ 1 is the relative position vector of m 2 with respect to m 1 . In addition, let F r ( r = 1 , 2 ) be the force between the masses m 2 and m 1 .
F 1 = F 2 = Γ 12 ρ 12 ρ ^ 12 ,
where ρ ^ 12 is the unit vector in the direction of ρ 12 and Γ 12 is the perturbed potential experienced by the second body given by Equation (1). Thus, the equations of motion of the two bodies are
m 1 ρ ¨ 1 = F 1 , m 2 ρ ¨ 2 = F 2 .
Utilizing Equations (1)–(3), the vectorial equation of motion of the mass m 2 with respect to the first mass m 1 is
ρ ¨ 12 = m 1 + m 2 m 1 m 2 Γ 12 .
In 2dimensions, the operator ∇ can be written in the polar coordinates ( ρ 12 , φ ) as
= ρ ^ 12 ρ 12 + φ ^ 1 ρ 12 φ ,
where ρ 12 is the radial distance and φ [ 0 , 2 π ] is the polar angle. Of course, ρ ^ 12 and φ ^ are two orthogonal unit vectors such that the first refers to the direction of the second body related to the first one, and φ ^ indicates the direction of the increasing polar angle φ . Hence, the vector acceleration in polar coordinates is identified by
ρ ¨ 12 = ρ ¨ 12 ρ 12 φ ˙ 2 ρ ^ 12 + 1 ρ 12 d d t ρ 12 2 φ ˙ φ ^ .
Substituting Equations (5) and (6) into (4), we obtain
ρ ¨ 12 ρ 12 φ ˙ 2 ρ ^ 12 + 1 ρ 12 d d t ρ 12 2 φ ˙ φ ^ = ( m 1 + m 2 ) m 1 m 2 Γ 12 ρ 12 ρ ^ 12 .
From Equation (1), we obtain
Γ 12 ρ 12 = G m 1 m 1 1 ρ 12 2 ε 12 d Υ 12 d ρ i j .
Now, utilizing Equations (7) and (8), we obtain
ρ ¨ 12 ρ 12 φ ˙ 2 ρ ^ 12 + 1 ρ 12 d d t ρ 12 2 φ ˙ φ ^ = G ( m 1 + m 1 ) 1 ρ 12 2 ε 12 d Υ 12 d ρ i j ρ ^ 12 .
Equation (9) has particular solutions as ρ 12 = constant and φ ˙ = ϖ , where ϖ is the mean motion, which can be written as
ϖ 2 = G ( m 1 + m 1 ) ρ 12 1 ρ 12 2 ε 12 d Υ 12 d ρ i j .
Expression (10) represents the general formula for the perturbed mean motion when the two bodies are moving under the effect of some conservative forces, such as a lack of sphericity, oblateness, and triaxial and radiation pressure effects [28,29,30,31].
Kepler’s laws can be employed with Newton’s laws in the case where the force is oriented toward a central fixed point and proportional to the inverse of distance square. Although the sun is considered the source of the main force in solar system, it is not fixed; rather, it experiences small accelerations due to the motions of the planets. In addition, the planets attract each other, and their motions are affected in accordance with Newton’s second and third laws. Thus, the Sun is not only the source of effective force on the motions of the planets. However, the other planets disturb the elliptical motion of the planet that would have happened if the planet was the only body moving around the Sun. Thereby, Kepler’s laws are considered as the only approximation for real motion, which motivated us to find the mean motion in the perturbed sense.

3. Restricted Three-Body Motion

3.1. Unperturbed Restricted Three–Body Motion

We assume that m 1 and m 2 are the two primary masses and m is the mass of the test particle, where its positions vectors with respect to the masses m 1 and m 2 are ρ 1 = ( x 1 , y 1 , z 1 ) , ρ 2 = ( x 2 , y 2 , z 2 ) while its position with respect to the origin of the inertial reference frame O X Y Z is ρ = ( x , y , z ) . We also suppose that the test particle does not affect the motion of the primaries and moves under their mutual gravitational.
Now, the motion of the test particle can be written in the following vectorial form:
m ρ ¨ = V I
where V I is the potential function with the following form:
V I ( ρ 1 , ρ 2 ) = m m 1 ρ 1 + m 2 ρ 2
The relative distances of the test particle are given by
ρ 1 2 = ( x x 1 ) 2 + ( y y 1 ) 2 + ( z z 1 ) 2 ρ 2 2 = ( x x 2 ) 2 + ( y y 2 ) 2 + ( z z 2 ) 2 ρ 2 = ξ 2 + y 2 + z 2
We intend to study the periodic orbits of the perturbed motion for the test particle, which has negligible mass with respect to the primary masses. Thus, we consider a synodic frame O ξ η ζ , which is moving around the ζ −axis with an angular velocity ϖ where ζ and Z axes are congruent.
To normalise the variable and use dimensionless coordinates, it is assumed that the universal constant of gravitation, the distance between the two primaries, the sum of their masses, and the angular velocity on the circular orbit are equal to one; see for details [32]. Thus, the mass of the massive big primary is m 1 = 1 μ , and the mass of the small primary is m 2 = μ , with μ = m 2 / ( m 1 + m 2 ) ( 0 , 1 / 2 ) . In synodical coordinates, the position of the massive primary is ( ξ 1 , η 1 , ζ 1 ) = ( μ , 0 , 0 ) and the position of the small one is ( ξ 2 , η 2 , ζ 2 ) = ( 1 μ , 0 ) , while the position of the test particle is ( ξ , η , ζ ) . Furthermore, we denote the magnitudes of the vectors ρ , ρ 1 , and ρ by ρ , ρ 2 , and ρ 2 , respectively. Then, the following vectorial equation of motion of the test particle in the framework of the classical restricted three-body problem in normalized synodic coordinates is given in [32]:
ρ ¨ + 2 ϖ ρ ˙ = V ,
where
V = 1 2 | ϖ ρ | 2 + 1 μ ρ 1 + μ ρ 2 .
The angular velocity and the position vectors in synodic coordinates are defined as follows:
ϖ = [ 0 0 ϖ ] T , ρ 1 = [ ( ξ + μ ) η ζ ] T , ρ 2 = [ ( ξ + μ 1 ) η ζ ] T , ρ = [ ξ η ζ ] T , ρ ˙ = [ ξ ˙ η ˙ ζ ˙ ] T ,
where ϖ equals one. In Equation (16), the distances among the infinitesimal body and the primaries are
ρ 1 2 = ( ξ + μ ) 2 + η 2 + ζ 2 , ρ 2 2 = [ ( ξ + μ 1 ) 2 + η 2 + ζ 2 , ρ 2 = ξ 2 + η 2 + ζ 2 .

3.2. Perturbed Restricted Three–Body Motion

In the case that the motion of the infinitesimal body is perturbed by very small forces with respect to the main gravitational forces of the primary bodies, the infinitesimal body will be affected by an additional acceleration in a synodic frame, which is called the perturbing acceleration A . Using Equations (8), (14) and (15), the perturbed motion of the infinitesimal body in synodic frame is given by
ρ ¨ + 2 ϖ ρ ˙ = Ψ ,
where Ψ = ( V + A ) , and A may be written as A = A 1 + A 2 + A 3 . Here, the magnitudes of the perturbing accelerations are A 1 = ε 1 Υ ( ρ 1 ) and A 2 = ε 2 Υ ( ρ 2 ) due to the massive and smaller primaries, respectively, while A 3 = ε 3 Υ ( ρ ) , may represent some additional perturbed forces, for example coming from asteroids belt. We mean that the first two perturbed forces are regarding to the two primaries, such that they may be are radiating bodies or have non–sphericity in their shapes (i.e. they are irregular bodies, these mean that they are oblate or triaxial bodies). In small word, the sources of all perturbations are conservative forces, which lead to total energy of the perturbed system is always has a constant amount. Hence, Ψ will be given by
Ψ = 1 2 ϖ 2 ( ξ 2 + η 2 ) + 1 μ ρ 1 + μ ρ 2 + A ,
where ρ 1 and ρ 2 are given in (17). Using Equation (10), the perturbed mean motion ϖ can be calculated by
ϖ 2 = 1 ε 12 Υ 12 ,
where Υ 12 = d Υ 12 / d ρ 12 when ρ 12 = 1 .
Utilizing Equations (18) and (19), the three dimensional equations of motion in the Cartesian synodic frame are controlled by
ξ ¨ 2 ϖ η ˙ = Ψ ξ , η ¨ + 2 ϖ ξ ˙ = Ψ η , ζ ¨ = Ψ ζ ,
where Ψ ξ = Ψ / ξ , Ψ η = Ψ / η , Ψ ζ = Ψ / ζ .
Equations (21) represents the dynamical system of the perturbed motion of the infinitesimal body in the framework of the spatial-restricted three-body problem when the perturbing forces are conservative. Thus, Equations (19) and (20) admit a Jacobian integral of the following form:
ξ ˙ 2 + η ˙ 2 + ζ ˙ 2 = 2 Ψ C ,
where C is the Jacobian constant, which can be used to identify the real motion as well as the possible regions of motion.

4. Hamiltonian of Perturbed Motion

In classical mechanics, the Hamiltonian function is a mathematical expression constructed in 1835 by Sir William Rowan Hamilton to characterize the evolution of a physical system that is composed of a set of differential equations. The features of this expression are that it provides significant insight into the properties of dynamical motion even if the physical system has neither an analytical solution nor is solvable. One of the famous examples is three-body systems in planetary and solar systems or stellar systems of motions. These systems have neither analytical nor closed-form solutions to their own systems, where Poincaré investigated that these systems revelation deterministic chaos in motion.
The total energy of a physical system (i.e., the sum of its kinetic energy and its potential energy) is identified by the Hamiltonian function. If this function does not depend on time, then the total energy of the system is conserved. One of the important aspects of this feature is that an infinitesimal phase–space volume is preserved [33]. A direct result to this property is Liouville’s theorem, which states that in the Hamiltonian system, the phase–space volume of a closed surface is preserved under time evolution.
In general, Hamiltonian systems are also called canonical systems, and if it has the properties of the autonomous system (we mean that the Hamilton function is not an explicit function of time), it may be indicated as a conservative system because the total energy of system will be a constant along time of motion and the Hamiltonian functions will also be called the first integral of motion. Thus, we will use the Hamiltonian function to describe the perturbed motion of the restricted three-body problem to prove the main result of this study, which is included in Theorem 1.
From [34] and using Equations (14) and (15), the Hamiltonian ( H 0 ) of the spatial unperturbed circular restricted 3 −body problem in synodical coordinates is
H 0 = 1 2 ( P ξ 2 + P η 2 + P ζ 2 ) + η P ξ ξ P η 1 μ ρ 1 μ ρ 2 ,
where H 0 H 0 ( ξ , η , ζ , P ξ , P η , P ζ ) , and ( P ξ , P η , P ζ ) is the conjugate momentum of the infinitesimal body at the position ( ξ , η , ζ ) .
In the case that the infinitesimal body is moving under some conservative perturbing force or its motion is governed by System (21) and admits the Jacobian integral (22), the Hamiltonian of the perturbed motion ( H H ( ξ , η , ζ , P ξ , P η , P ζ ) ) is given by
H = H 0 + ε 12 Υ 12 η P ξ ξ P η + ε 1 Υ ( ρ 1 ) + ε 2 Υ ( ρ 2 ) .
It is clear that the Hamiltonians (23) and (24) are equivalent in the absence of the perturbing forces.
Since the quantities ε 12 , ε 1 , ε 2 , and ε 3 are small, they can be written with a factor λ 1 , i.e., ε 12 = λ ε ¯ 12 , ε 1 = λ ε ¯ 1 , ε 2 = λ ε ¯ 2 and ε 3 = λ ε ¯ 3 . In this context, the Hamiltonian (24) becomes
H = 1 2 ( P ξ 2 + P η 2 + P ζ 2 ) + η P ξ ξ P η 1 μ ρ 1 μ ρ 2 + O ( λ ) .
We remark that the unperturbed motion of the infinitesimal body in the framework of the restricted three–body problem can be represented by the Hamiltonian (25) in the case of the whole perturbed forces are ignored. Moreover, if the mass ratio μ is very small, such that μ = λ μ ¯ , we have that the Hamiltonian (25) becomes
H = 1 2 ( P ξ 2 + P η 2 + P ζ 2 ) + η P ξ ξ P η 1 ξ 2 + η 2 + ζ 2 + O ( λ ) .
We must take care with O ( λ ) because it has terms that go to infinity near the primaries; therefore, we must exclude a neighborhood of the primaries. When the infinitesimal body is moving in the primaries’ plane, we have that ζ = 0 and P ζ = 0 ; hence, the Hamiltonian (26) becomes
H ( ξ , η , P ξ , P η ) = 1 2 ( P ξ 2 + P η 2 ) + η P ξ ξ P η 1 ξ 2 + η 2 + O ( λ ) .
Then, equations of motion related to the Hamiltonian (27) will take the following form:
ξ ¨ = 2 n η ˙ + n 2 ξ ξ ξ 2 + η 2 3 / 2 + O ( λ ) , η ¨ = 2 n ξ ˙ + n 2 η η ξ 2 + η 2 3 / 2 + O ( λ ) .
Note that the Hamiltonian (27) and System (28) for λ = 0 coincide with the planar Kepler problem in a rotating reference frame. Hence, we will conclude the main result of this study in the following theorem:
Theorem 1.
The circular and elliptic periodic orbits of the planar rotating Kepler problem can be prolonged to the periodic orbits of first and second type for the perturbed two-body problem and the perturbed planar circular restricted three-body problem when the perturbation is made with conservative forces.

5. Periodic Orbits of Perturbed Motion

In celestial mechanics, changes in two-body orbits occur rather slowly with respect to the perturbations, and it is considered a good first approximation for the perturbed motion. The perturbations can be applied if the perturbed forces are very small with respect to the main or fundamental gravitational forces. General perturbation techniques are preferred for some types of perturbed motions as the source of a certain predestined motion is easily found. This is not necessary with respect to particular perturbations where the motion would be predicted with similar reliability. However, there are no data on the configuration of the perturbing bodies.
In the solar system, planets are subjected to periodic perturbations from other planets, which pass in their orbit; these perturbations are composed of small impulses each time a planet crosses or passes the orbit of the other planet. This leads to the motion of bodies being periodic or quasi–periodic, such as the Moon in its strongly perturbed orbit. Thus, in the next subsection, we study the existence of first- and second-type periodic orbits for the perturbed problem.
Proof of Theorem 1. 
It is clear that System (28) represents Kepler motion in the rotating reference frame, also known as the so-called rotation Kepler problem when the parameter-to-mass ratio ( μ equals zero. This means that the third body or the infinitesimal body will be moving under the effect of the gravitational field of the massive body only. In this case, the massive body will take its place at the origin of the reference frame.
We also remark that if λ = 0 , the periodic solutions of System (28) are circular or elliptical orbits, namely ( ξ ( t ) , ( η ( t ) ) when the motion of the two-body system is bounded. In general, these orbits are called the first (second) type when the unperturbed motion is represented by circular (elliptical) orbits.
We emphasize that System (28) is the same or similar to the system in page 4 Equation (6) of [2]. Therefore in the case of λ = 0 (or λ 0 , but small, which means that the mass ratio μ is also very small, and μ 0 ), the obtained results in [2] can be applied to our system of perturbed motion (28). Hence, the results of Theorems 3.1, 3.2 and 5.1 in [2] can be applied to perturbed circular restricted three-body motion. These theorems will be stated and indexed in the present work as Theorems 2–4: □
Theorem 2.
Assume that the restricted three-body problem is perturbed through a very small perturbed force and that the perturbed system has a first integral. Additionally, we assume that 1 / ( 1 n j 3 ) Z , where j 3 1 / n and Z is the set of integer numbers. Then, the periodic orbits of first type with angular momentum j can be continued to the perturbed circular restricted three-body problem.
Theorem 3.
Assume that the restricted three-body problem is perturbed through a small perturbing force and that the perturbed system has a first integral. If τ 0 ( l + 1 / 2 ) π / n with l Z + , where Z + is the set of positive integer numbers and ϵ > 0 is very small, then the perturbed restricted problem has periodic orbits, which tend to the periodic orbits of the first type during the τ 0 period when ϵ tends to zero.
Theorem 4.
Assume that I and J are prime integers, K Z , where Z is the set of integer numbers, and τ = 2 π I / | K | . Then, the elliptical solution with a τ period for the rotating Kepler problem satisfying
g ( 0 ) = π , l ( 0 ) = π , L 3 ( 0 ) = I / J ,
and such that it does not cross the location of the smaller primary can be continued to the perturbed circular restricted three-body problem when ϵ > 0 is small and its motion Equations (9) (in Reference [2]) come from Equations (6) (in Reference [2]), satisfying the following invariant symmetry by ( x , y , t ) ( x , y , t ) .
Furthermore, the proof of Theorem 1 comes directly from the proofs of Theorems 3.1 and 5.1 in [2], here Theorems 2 and 4. Thereby, the circular (elliptical) periodic orbits of the planar rotating Kepler problem can be prolonged in the periodic orbits of first (second) types for the perturbed two-body problem and for the perturbed restricted three-body problem by conservative perturbing forces. This means that both theorems are satisfied when the perturbed model has a first integral or the perturbed motion has the so-called Jacobian constant.
Now, we show that Theorem 1 can be applied to many models, for example, to the perturbed two–body problem, by different perturbed conservative forces, some of these models are listed in what follows.
  • Perturbed two-body problem where both bodies are oblate [35].
  • An anisotropic perturbed Kepler problem [36].
  • Relativistic effect is a source of perturbation [37].
  • Perturbed two-body problem when the main body has a continuation fraction potential. [38].
We can also list some related works for the perturbed restricted three-body problem, as follows:
  • The oblateness of a massive primary body is a source of perturbed force [28].
  • The oblateness of a smaller primary and the radiation of a bigger primary are sources of perturbed forces [29].
  • A lack of sphericity of the three participants is a source of perturbed forces as well as the primaries’ radiations [30].
  • The effect of zonal harmonic coefficients ( J 2 and J 4 ) of the massive body are considered [39].
  • The first primary is oblate, while the second is a source of radiation [40].
  • Both primaries are oblate and radiating bodies, as well as subject to the asteroid belt effect [41].
  • The effects of zonal harmonic coefficients ( J 2 and J 4 ) of both primary bodies are considered [42].
  • The triaxial of a massive primary body and the oblateness of a smaller primary [43].
  • Both primaries are taken as triaxial rigid bodies [31].
  • The asteroid belt effect is considered a source of a perturbation force [44].
  • Quantum corrections are considered a source of perturbation forces [45].
The perturbed models for two and three bodies with their own first integral are not limited to those stated in the previous two lists, which are mentioned as examples. Furthermore, we demonstrate that this study is an extension and generalization of the obtained results in [2,46]. Theorem 1 can be applied in the case that the system is subjected to conservative forces, by which we mean that all forces can be derivable from the total energy of system, which must be a constant. However, the principle of conservation of energy is an expression relating the velocity of the mass and its own position at any time it has a fixed value for all time.

6. Conclusions

In this work, an overview of periodic orbits with their importance is stated in terms of both celestial mechanics and astrodynamics. The general perturbed mean motion of the two-body problem is derived if the perturbation forces are conservative, such as radiation pressure oblateness and the asteroid belt effects. Thereby, the unperturbed and perturbed motion of the test particle in the framework of the restricted three-body problem is also derived under the conservative perturbed forces. Thus, we have proved that the first and second types of periodic orbits of the rotating Kepler problem can persist for all perturbed two-body and circular restricted three-body problems when the perturbation forces are conservative or the perturbed motion has its own extended Jacobian integral. We emphasize that the terminologies of Poincaré and canonical transformation using Delaunay’s variables can be applied to prove the obtained results.

Author Contributions

Formal analysis: E.I.A., J.L.G.G. and J.L.; investigation: E.I.A., J.L.G.G.; methodology: E.I.A., J.L.G.G. and J.L.; project administration: J.L.; software: E.I.A., J.L.G.G. and J.L.; validation: E.I.A., J.L.G.G. and J.L.; visualization: E.I.A., J.L.G.G. and J.L.; writing—original draft: E.I.A. and J.L.; writing—review & editing: E.I.A. and J.L. All authors have read and agreed to the published version of the manuscript.

Funding

The third author is partially supported by the Agencia Estatal de Investigación grant PID2019-104658GB-I00, the H2020 European Research Council grant MSCA-RISE-2017-777911, the AGAUR (Generalitat de Catalunya) grant 2022-SGR 00113, and by the Acadèmia de Ciències i Arts de Barcelona.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study does not report any data.

Acknowledgments

This work was funded by the National Research Institute of Astronomy and Geophysics (NRIAG), Helwan 11421, Cairo, Egypt. The first author, therefore, acknowledges with thanks NRIAG technical and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Darwin, G. Periodic Orbits. Acta Math. 1897, 21, 99–242. [Google Scholar] [CrossRef]
  2. Abouelmagd, E.I.; Guirao, J.L.G.; Llibre, J. Periodic orbits for the perturbed planar circular restricted 3–body problem. Discret. Contin. Dyn. Syst.–B 2019, 24, 1007–1020. [Google Scholar] [CrossRef]
  3. de Bustos, M.T.; Fernández, A.; López, M.A.; Martínez, R.; Vera, J.A. Periodic Orbits of Third Kind in the Zonal J2 + J3 Problem. Symmetry 2019, 11, 111. [Google Scholar] [CrossRef]
  4. Krasilnikov, P.; Saraeva, A. Poincaré periodic orbits of the first kind in the planar circular restricted problem of three bodies with small acceleration. Cosmic Res. 2015, 53, 469–475. [Google Scholar] [CrossRef]
  5. Bloom, H. The National Polar-orbiting Operational Environmental Satellite System. In Proceedings of the 35th COSPAR Scientific Assembly, Paris, France, 18–25 July 2004; Volume 35, p. 2230. [Google Scholar]
  6. Lee, H.; Grogan, P.T. Measuring Strategic Risk Dominance Using the Multi-actor Value Model: A Study of the National Polar-orbiting Operational Environmental Satellite System. In Proceedings of the AIAA Scitech 2019 Forum, San Diego, CA, USA, 7–11 January 2019; p. 1031. [Google Scholar]
  7. Farkas, M. Determination of controllably periodic perturbed solutions by Poincaré’s method. Stud. Sci. Math. Hung. 1972, 7, 257–266. [Google Scholar]
  8. Hale, J.K.; Táboas, P.Z. Interaction of damping and forcing in a second order equation. Nonlinear Anal. Theory Methods Appl. 1978, 2, 77–84. [Google Scholar] [CrossRef]
  9. Ermentrout, G.B. n:m Phase-locking of weakly coupled oscillators. J. Math. Biol. 1981, 12, 327–342. [Google Scholar] [CrossRef]
  10. Chicone, C.; Rhouma, M. On the Continuation of Periodic Orbits. Methods Appl. Anal. 2000, 7, 85–104. [Google Scholar]
  11. Chicone, C. Lyapunov-Schmidt reduction and Melnikov integrals for bifurcation of periodic solutions in coupled oscillators. J. Differ. Equ. 1994, 112, 407–447. [Google Scholar] [CrossRef]
  12. Chicone, C. Periodic solutions of a system of coupled oscillators near resonance. SIAM J. Math. Anal. 1995, 26, 1257–1283. [Google Scholar] [CrossRef]
  13. Chicone, C. A geometric approach to regular perturbation theory with an application to hydrodynamics. Trans. Am. Math. Soc. 1995, 347, 4559–4598. [Google Scholar] [CrossRef]
  14. Barrar, R. Existence of periodic orbits of the second kind in the restricted problems of three bodies. Astron. J. 1965, 70, 3–4. [Google Scholar] [CrossRef]
  15. Sharma, R.K. Periodic orbits of the second kind in the restricted three–body problem when the more massive primary is an oblate spheroid. Astrophys. Space Sci. 1981, 76, 255–258. [Google Scholar] [CrossRef]
  16. Ahmad, A. Existence of periodic orbits offirst kind in the photogravitational circular restricted problem of four bodies. Bull. Astr. Soc. India 1995, 23, 165–176. [Google Scholar]
  17. Pan, S.; Hou, X. Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem. Appl. Sci. 2022, 12, 8952. [Google Scholar] [CrossRef]
  18. Koon, W.S.; Lo, M.W.; Marsden, J.E.; Ross, S.D. Heteroclinic connections between periodic orbits and resonance transitions in celestial mechanics. Chaos Interdiscip. J. Nonlinear Sci. 2000, 10, 427–469. [Google Scholar] [CrossRef]
  19. Henon, M. Families of asymmetric periodic orbits in Hill’s problem of three bodies. Celest. Mech. Dyn. Astron. 2005, 93, 87–100. [Google Scholar] [CrossRef]
  20. Kalantonis, V.; Perdios, E.; Perdiou, A. The Sitnikov family and the associated families of 3D periodic orbits in the photogravitational RTBP with oblateness. Astrophys. Space Sci. 2008, 315, 323–334. [Google Scholar] [CrossRef]
  21. Dutt, P.; Sharma, R. Analysis of periodic and quasi–periodic orbits in the Earth–Moon system. J. Guid. Control Dyn. 2010, 33, 1010–1017. [Google Scholar] [CrossRef]
  22. Perdiou, A.; Perdios, E.; Kalantonis, V. Periodic orbits of the Hill problem with radiation and oblateness. Astrophys. Space Sci. 2012, 342, 19–30. [Google Scholar] [CrossRef]
  23. Gong, S.; Li, J. Solar sail periodic orbits in the elliptic restricted three-body problem. Celest. Mech. Dyn. Astron. 2015, 121, 121–137. [Google Scholar] [CrossRef]
  24. Kalantonis, V.S. Numerical investigation for periodic orbits in the Hill three-body problem. Universe 2020, 6, 72. [Google Scholar] [CrossRef]
  25. Kotoulas, T.; Voyatzis, G. Planar retrograde periodic orbits of the asteroids trapped in two–body mean motion resonances with Jupiter. Planet. Space Sci. 2020, 182, 104846. [Google Scholar] [CrossRef]
  26. Alberti, A.; Vidal, C. First kind symmetric periodic solutions and their stability for the Kepler problem and anisotropic Kepler problem plus generalized anisotropic perturbation. Nonlinear Anal. Real World Appl. 2021, 58, 103238. [Google Scholar] [CrossRef]
  27. Yousuf, S.; Kishor, R. Families of periodic orbits about Lagrangian points L1, L2 and L3 with continuation method. Planet. Space Sci. 2022, 217, 105491. [Google Scholar] [CrossRef]
  28. Subbarao, P.; Sharma, R.K. A note on the stability of the triangular points of equilibrium in the restricted three-body problem. Astron. Astrophys. 1975, 43, 381–383. [Google Scholar]
  29. Sharma, R.K. The linear stability of libration points of the photogravitational restricted three–body problem when the smaller primary is an oblate spheroid. Astrophys. Space Sci. 1987, 135, 271–281. [Google Scholar] [CrossRef]
  30. Abouelmagd, E.I.; El-Shaboury, S. Periodic orbits under combined effects of oblateness and radiation in the restricted problem of three bodies. Astrophys. Space Sci. 2012, 341, 331–341. [Google Scholar] [CrossRef]
  31. Alzahrani, F.; Abouelmagd, E.I.; Guirao, J.L.; Hobiny, A. On the libration collinear points in the restricted three–body problem. Open Phys. 2017, 15, 58–67. [Google Scholar] [CrossRef]
  32. Murray, C.D.; Dermott, S.F. Solar System Dynamics; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  33. De Almeida, A.M.O. Hamiltonian Systems: Chaos and Quantization; Cambridge University Press: Cambridge, UK, 1988. [Google Scholar]
  34. Meyer, K.R.; Offin, D.C. Introduction to Hamiltonian Dynamical Systems and the N-Body Problem; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  35. Abouelmagd, E.I.; Mortari, D.; Selim, H.H. Analytical study of periodic solutions on perturbed equatorial two–body problem. Int. J. Bifurc. Chaos 2015, 25, 1540040. [Google Scholar] [CrossRef]
  36. Abouelmagd, E.I.; Llibre, J.; Guirao, J.L.G. Periodic orbits of the planar anisotropic Kepler problem. Int. J. Bifurc. Chaos 2017, 27, 1750039. [Google Scholar] [CrossRef]
  37. Abouelmagd, E.I.; Elshaboury, S.; Selim, H. Numerical integration of a relativistic two–body problem via a multiple scales method. Astrophys. Space Sci. 2016, 361, 38. [Google Scholar] [CrossRef]
  38. Abouelmagd, E.I. Periodic solution of the two–body problem by KB averaging method within frame of the modified newtonian potential. J. Astronaut. Sci. 2018, 65, 291–306. [Google Scholar] [CrossRef]
  39. Abouelmagd, E.I. Existence and stability of triangular points in the restricted three–body problem with numerical applications. Astrophys. Space Sci. 2012, 342, 45–53. [Google Scholar] [CrossRef]
  40. Abouelmagd, E.I.; Sharaf, M. The motion around the libration points in the restricted three–body problem with the effect of radiation and oblateness. Astrophys. Space Sci. 2013, 344, 321–332. [Google Scholar] [CrossRef]
  41. Singh, J.; Taura, J.J. Motion in the generalized restricted three-body problem. Astrophys. Space Sci. 2013, 343, 95–106. [Google Scholar] [CrossRef]
  42. Abouelmagd, E.I.; Alhothuali, M.; Guirao, J.L.; Malaikah, H. The effect of zonal harmonic coefficients in the framework of the restricted three–body problem. Adv. Space Res. 2015, 55, 1660–1672. [Google Scholar] [CrossRef]
  43. Singh, J.; Kalantonis, V.; Gyegwe, J.M.; Perdiou, A. Periodic motions around the collinear equilibrium points of the R3BP where the primary is a triaxial rigid body and the secondary is an oblate spheroid. Astrophys. J. Suppl. Ser. 2016, 227, 13. [Google Scholar] [CrossRef]
  44. Abozaid, A.A.; Selim, H.; Gadallah, K.A.; Hassan, I.; Abouelmagd, E.I. Periodic orbit in the frame work of restricted three bodies under the asteroids belt effect. Appl. Math. Nonlinear Sci. 2020, 5, 157–176. [Google Scholar] [CrossRef]
  45. Alshaery, A.; Abouelmagd, E.I. Analysis of the spatial quantized three-body problem. Results Phys. 2020, 17, 103067. [Google Scholar] [CrossRef]
  46. Abouelmagd, E.I.; García Guirao, J.L.; Llibre, J. Periodic Orbits of Quantised Restricted Three-Body Problem. Universe 2023, 9, 149. [Google Scholar] [CrossRef]
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

Abouelmagd, E.I.; Guirao, J.L.G.; Llibre, J. On the Periodic Orbits of the Perturbed Two- and Three-Body Problems. Galaxies 2023, 11, 58. https://doi.org/10.3390/galaxies11020058

AMA Style

Abouelmagd EI, Guirao JLG, Llibre J. On the Periodic Orbits of the Perturbed Two- and Three-Body Problems. Galaxies. 2023; 11(2):58. https://doi.org/10.3390/galaxies11020058

Chicago/Turabian Style

Abouelmagd, Elbaz I., Juan Luis García Guirao, and Jaume Llibre. 2023. "On the Periodic Orbits of the Perturbed Two- and Three-Body Problems" Galaxies 11, no. 2: 58. https://doi.org/10.3390/galaxies11020058

APA Style

Abouelmagd, E. I., Guirao, J. L. G., & Llibre, J. (2023). On the Periodic Orbits of the Perturbed Two- and Three-Body Problems. Galaxies, 11(2), 58. https://doi.org/10.3390/galaxies11020058

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