## 1 Introduction

The three-body problem has a special relevance, particularly in astrophysics and astrodynamics. In general, the three-body problem is classified into two types:

- (i) The general three-body problem, which describes the motion of three bodies of arbitrary masses under their mutual attraction due to the gravitational field. The motion of the bodies take place in three dimensions and there are no restrictions on their masses nor on the initial conditions.
- (ii) The restricted three-body problem in which the mass of the third body is very small in comparison with the masses of the primaries, and it does not affect their motion. In this case, the primaries move around their center of mass either along circular or elliptical orbits.

Thus, the general problem has some applications in celestial mechanics such as the dynamics of triple star systems (and only a very few in space dynamics and solar system dynamics), whereas the restricted problem plays an important role in studying the motion of artificial satellites. It can be used also to evaluate the motion of the planets, minor planets and comets. The restricted problem gives an accurate description not only regarding the motion of the Moon but also with respect to the motion of other natural satellites. Furthermore, the restricted problem has many applications not only in celestial mechanics but also in physics, mathematics and quantum mechanics, to name a few. In quantum mechanics, a general form of the restricted problem is formed to solve the Schrödinger equation of helium-like ions. Furthermore, in modern solid state physics, the restricted problem can be used to discuss the motion of an infinitesimal mass affected not only by the gravitational field but also by light pressure from one (or both) of the primaries, which is called the photogravitational problem.

In the last decades, the restricted three-body problem has been enhanced by a great number of researches. Many of them deal with the effects of the perturbed forces such as the lack of sphericity, photogravitational force, Coriolis and centrifugal forces, variation of masses, the Pointing-Robertson and Yarkovsky effects and the atmospheric drag and solar wind. The Kirkwood gaps in the ring of the asteroid’s orbits lying between the orbits of Mars and Jupiter are examples of the perturbation produced by Jupiter on an asteroid. This motivates many researchers to study the restricted three-body problem under the effects of small perturbations in Coriolis and centrifugal forces when the three involved objects are oblate, . . . , etc.

The significance of the restricted problem when describing actual physical situations can be judged by the results obtained when these are compared with observations. It is worth mentioning that the utility might be prejudged by order of magnitude evaluations regarding the masses and the distances of the participating bodies. In this way, a classic example in space dynamics is the Sun-Earth-Moon system.

Some important contributions related to the libration points in the restricted three-body problem with one or both primaries are oblate spheroids when the equatorial plane is coincident with the plane of motion were studied by Subbarao and Sharma [28], Sharma and Subbarao [23], and Markellos et al. [20]. Abouelmagd [9] also studied the effects of oblateness *J*_{2} and *J*_{4} for the more massive primary in the planar restricted three-body problem on the locations of the triangular points and their linear stability. He found that these locations are affected by the coefficients of oblateness. Furthermore, he investigated that the triangular points are stable for 0 < *μ* < *μ _{c}* and unstable when

*μ*≤

_{c}*μ*≤ 1/2, where

*μ*is the critical mass parameter which depends on the coefficients of oblateness.

_{c}The existence of libration points and their linear stability as well as periodic orbits around these points when the more massive primary is radiating and the smaller is an oblate spheroid were studied by Abouelmagd and Sharaf [3]. Their study also includes the effects of oblateness of

The restricted problem when the three participating bodies are oblate spheroids was also studied by Elipe and Ferrer [15], and El-Shaboury and El-Tantawy [18]. When one or two of the primaries are triaxial bodies this problem was introduced by El-Shaboury et al. [17], Khanna and Bhatnagar [19], and Sharma et al. [24].

Various researchers made studies in the restricted three-body problem under the effects of small perturbations in centrifugal and Coriolis forces such as in Szebehely [29], Bhatnagar and Hallan [12], Devi and Singh [14], and Shu and Lu [25], to quote some of them.

The effect of small perturbations *ε*,*ε*′ in Coriolis and centrifugal forces with variable mass in the restricted three-body problem has been studied by Singh [26]. He found that in the nonlinear sense the triangular points are stable for all mass ratios in the range of linear stability except for three mass ratios which depend on *ε*,*ε*′ and the constant *β* due to the variation of mass governed by Jeans’ law.

Mittal et al. [21] have studied periodic orbits generated by Lagrangian solutions of the restricted three-body problem when the bigger body is a source of radiation and the smaller is an oblate spheroid. They used the definition of Karimov and Sokolsky for mobile coordinates to determine these orbits, they also used the predictor method to draw them.

Singh and Begha [27] have studied the existence of periodic orbits around the triangular points in the restricted three-body problem when the bigger primary is a triaxial and the smaller primary is considered as an oblate spheroid in the range of linear stability with the perturbed Coriolis and centrifugal forces. They deduced that long and short periodic orbits exist around these points and their periods, orientation and eccentricities are affected by the non sphericity and the perturbation in Coriolis and centrifugal forces.

The existence of libration points and their linear stability when the three participating bodies are oblate spheroids and the primaries are radiation source as well were studied by Abouelmagd and El-Shaboury [2]. They found that the collinear points are still unstable while the triangular points are stable for 0 < *μ* < *μ _{c}*, and unstable for

*μ*≤

_{c}*μ*≤ 1/2, where

*μ*∊ (0,1/2). They also deduced that for these points the range of stability will decrease. In addition they studied the periodic orbits around these points in the range 0 <

_{c}*μ*<

*μ*.

_{c}From a physical point of view, it is unreasonable to consider all objects as being point masses with no physical dimensions. This is conflict with the real cases for the celestial bodies. On the other hand the effect of rotation causes deformation in the shape of the objects at the poles as might be expected. For this reason most objects may be treated to a good approximation as oblate spheroids, in this context there are a number of contributions established by: Abouelmagd [5, 6], Abouelmagd et al. [7, 8], and Abouelmagd et al. [9–11]. Therefore in this paper we will generalize the restricted three-body, in which the three participating bodies are oblate spheroids as well as the existence of small perturbations in the Coriolis and centrifugal forces. The existence of libration points and their linear stability and the periodic orbits around these points will be studied. This model could be applicable in astrodynamics or astrophysics.

Now we assume that the three participating bodies are oblate spheroids. Therefore we refer to oblateness parameters of the bigger, smaller primaries and infinitesimal body by *A*_{1},*A*_{2}, and *A*, respectively. We also denote the small perturbations in Coriolis and centrifugal forces by *ε _{i}* :

*i*= 1,2, where 0 <

*A*

_{1}≪ 1, 0 <

*A*

_{2}≪ 1, 0 <

*A*≪ 1, and

*ε*≪ 1.

_{i}## 2 The equations of motion

It is well known that there are five exact solutions in three-body problem which the three masses maintain a constant configuration which revolves with constant angular velocity. An important specialization of the three-body problem is the restricted three-body problem in which *m* is infinitesimal body, and *m*_{1} and *m*_{2} are the bigger and smaller primaries, respectively. They move in circular orbits around their barycenter. The smallness of *m *means that it does not influence the motion of *m*_{1} and *m*_{2}. For many purposes, it is convenient to describe the motion of *m* in a coordinate system which is attached to both *m*_{1} and *m*_{2}. In this rotating coordinate system, the five Lagrangian solutions show up as five fixed points at which *m* would be stationary if placed there with zero velocity.

If the object has axial symmetry, it can be shown from potential theory that the gravitational potential experienced by the satellite can be written as in (Murray and Dermott, [22]), namely,

where

*G*is the universal constant.*m*_{0}and*R*_{0}are the mass and the mean radius of the object.*p*(sin_{n}*δ*) denote the Legendre polynomials of degree*n*.*r*is the distance from the center of the object to the satellite.*J*is the dimensionless coefficient that characterizes the size of non-spherical components of the potential._{n}

If the motion of the satellite and the object are in the same plane (*δ* = 0), then Eq. (1) becomes

where

Now we assume that *m*_{1},* m*_{2} and *m* are the masses of the bigger, smaller and infinitesimal bodies, *r*_{1} and *r*_{2} are the distances from *m*_{1} and *m*_{2} to *m*, respectively.

Let *m*_{1} and *m*_{2} have a circular orbit about their common center of mass and *m* is moving under their gravitational field in the same plane but it is so small that it does not affect their motion. Furthermore let:

=*R**ui*+*ν j*be the position vector of*m*_{2}with respect to*m*_{1}.*R*is the separation distance between*m*_{1}and*m*_{2}.

Using Eq. (3), then the potential *V*_{12} between *m*_{1} and *m*_{2} can be written as

Thus, the equation of motion for *m*_{2} is

where

Hence Eq. (5) has the following particular solutions:

- (i)
*R*= constant. - (ii)
*u*=*R*cos*θ*. - (iii)
*ν*=*R*sin*θ*. - (iv)
*θ*=*nt*.

Therefore, the perturbed mean motion is given by the following relation

Now we assume that *XY* is the sidereal frame while *xy* are the synodic frames such that they have the same origin at the center of mass for the primaries. We also chose the synodic frame to rotate with angular velocity *n* in positive direction (see Fig. (1)). Let the coordinates of *m*_{1},*m*_{2}, and *m* in the sidereal frame be (*X*_{1},* Y*_{1}), (*X*_{2},* Y*_{2}), and (*X, Y* ), respectively, but in the synodic frame are (*x*_{1},* y*_{1}),(*x*_{2},* y*_{2}), and (*x, y*), respectively, too.

Based on the above discussion, the sidereal coordinates are related to the synodic coordinates via the following relation:

Using Newton’s law, the equations of motion for infinitesimal body in a sidereal frame are as follows:

From Eq. (3), the potential *V*_{1} between *m* and *m*_{1} is given by

and the potential *V* between *m, m*_{1} and *m*_{2} is also given by

while the total potential *V* between *m, m*_{1}, and *m*_{2} is controlled by

once the notation and terminology of Szebehely [30] has been adopted. As a consequence, the distance between the primaries does not change and is equal to one, the sum of masses of the primaries is also taken as one, the unit time is chosen as to make the unperturbed mean motion and the gravitational constant unity, then it follows that *m*_{1} = 1 − *μ*, whose coordinate is (*μ*, 0),* m*_{2} = *μ* ≤ 1/2 and located at (*μ* − 1,0).

Substituting Eqs. (7-8) and (11-13) into (9-138) that gives the equations of motion in a synodic coordinate system (*x* − *y*) for infinitesimal body with dimensionless variables as

- (i)
*ẍ*− 2*nẏ*=*U*._{x} - (ii)
*ÿ*+ 2*nẋ*=*U*,_{y}

where

Eqs. (i)-(ii) and (14) represent the equation of motion for the infinitesimal body in a rotating coordinate system when the three participating bodies are oblate spheroids.

In the previous equations the perturbed effects of Coriolis and centrifugal forces are switched off. In physics the Coriolis effect is a deflection of moving objects when they are viewed in a rotating reference frame. In a reference frame with clockwise rotation the deflection is to the left of the motion of the object, in one with counter-clockwise rotation the deflection is to the right. The Newton’s laws of motion govern the motion of an object in a non-accelerating inertial frame of reference. When Newton’s laws are transformed to a uniformly rotating frame of reference, the Coriolis and centrifugal forces appear.

Both forces are proportional to the mass of the object. The Coriolis force is proportional to the rotation rate and the centrifugal force is proportional to its square. The Coriolis force acts in a direction perpendicular to the rotation axis and to the velocity of the body in the rotating frame and is proportional to the object’s speed in the rotating frame. The centrifugal force acts outwards in the radial direction and is proportional to the distance of the body from the axis of the rotating frame.

It is important to note that these forces do not arise from any physical agency: they arise solely as a result of rotation for the coordinate system. Thus, if the angular velocity is reduced to zero, both the Coriolis and the centrifugal forces will vanish. For this reason, the Coriolis and centrifugal forces are referred to as “fictitious” or “inertial” forces. These additional forces allow the application of Newton’s laws to a rotating system. They are correcting factors that do not exist in a non-accelerating or inertial reference frame.

Anyway, if the perturbed effects of Coriolis and centrifugal forces are considered, let *φ*, *ψ* denote these forces, respectively. Assume that *ε*_{1} represents the perturbation in Coriolis force, while *ε*_{2} is the accompanying perturbation in the centrifugal force such that

Hence, Eqs. (i)-(ii) and (14) can be written as follows:

where

It is worth mentioning that the subscripts on the right handside denote the partial derivation with respect to *x* and *y*, respectively.

## 3 The locations of the libration points

The libration points are singular points of the differential equations of motions in the restricted problem of three bodies. They are also equilibrium points since the gravitational force on a mass placed in such a point are balanced by the centrifugal force. The three libration points are called the collinear points which are found on the line connecting with the primaries. The other two are called the triangular points which are symmetrical with respect to this line and form triangles, see Fig. (1).

From Eqs. (i)-(91), we can obtain the Jacobi integral as

where *c* is a constant of integration.

The locations of these points will be determined by solving the following expressions:

where the first partial derivatives of Ω will be governed by

From Eqs. (25-26), two cases yield. They will be investigated in the following subsection.

### 3.1 The collinear points

The collinear points *L _{i}* :

*i*= 1, 2, 3 are the solutions of Eqs. (24) and (25-26) when

*y*= 0. Therefore,

where *i* indicates the value of *ν* at *L _{i}* :

*i*= 1, 2, 3. Further,

Lagrange developed a useful method for inverting series expansions which could be applicable in the present work. He investigated that if a variable *χ* can be expressed as a function of *𝒴* as in the below relation (Murray and Dermott, [22]):

Therefore, ϒ can be also expressed as a function of *χ* as follows:

### 3.1.1 Location of *L*_{1}

The solution of *L*_{1} lies beyond the smaller mass as in Fig. (2). Therefore, *r*_{1} − *r*_{2} = 1, and hence, *r*_{2} = *μ* − *x* − 1,* r*_{1} = *μ* − *x*, which implies that

provided that we take *r*_{2} = *r*, *r*_{1} = 1 + *r*, and substitute them into Eq. (29). After some simple calculations the value of *ν* will be written as

where the values *a*_{4},* a*_{5},* a*_{6} and *a*_{7} are given as in Appendix I.

Now, we use Lagrange’s inversion method to invert the above series and hence to express *r* as a function of *ν* which is written as

Consequently, the location of *L*_{1} will be given by

### 3.1.2 Location of *L*_{2}

The solution of *L*_{2} lies between the two finite bodies as in Fig. (2), where *r*_{1} + *r*_{2} = 1. Hence, *r*_{2} = *x* − *μ* + 1,* r*_{1} = *μ* − *x*. This implies that

Thus, if we take *r*_{2} = *ρ*, *r*_{1} = 1 − *ρ*, and substitute them into Eq. (29), then the value of *ν* will be given by

where the values *b*_{4},* b*_{5},* b*_{6}, and *b*_{7} are governed as in Appendix I.

As mentioned previously, we use Lagrange’s inversion method in order to invert the above series and hence, express *ρ* as a function of *ν*, which can be written as

Therefore, the location of *L*_{2} could be written as

### 3.1.3 Location of *L*_{3}

The solution of *L*_{3} lies beyond the large mass as in Fig. (2), where *r*_{2} − *r*_{1} = 1. Hence, *r*_{2} = 1 − *μ* + *x*, *r*_{1} = *x* − *μ*. This implies that

Thus, if we take *r*_{1} = 1 *−* +*δ*, *r*_{2} = 2 + *δ* , and substitute them into Eq. (29), then *v* is as follows:

where *c*_{0},* c*_{1} and *c*_{2} are given as in Appendix I.

We also use Lagrange’s inversion method to invert the above series and to express *δ* as a function of *v*. Hence, the following expression yields:

Thus, the location of *L*_{3} will be controlled by

Note that Eqs. (36), (38) and (41) represent semi-closed forms for the positions of collinear points. Furthermore, the small perturbation in the Coriolis force has no effect on these positions.

### 3.1.4 The triangular points

The solutions of Eqs. (24) and (25-26) when *y* ≠ 0 will give the triangular points *L*_{4} and *L*_{5}. These solutions will lead to

Therefore, we get

To study the effect of oblateness in the perturbed problem on the locations of the triangular points, we will assume that the three participating bodies are not oblate spheroids, namely, *A*_{1} = *A*_{2} = *A* = 0.

Hence Eqs. (43-44) will admit

where *π _{i}* is the additive displacement in

*r*as a result of oblateness effects. Substituting Eq. (45) into (43-44) and restricting ourselves to the linear terms in

_{i}*π*

_{1},

*π*

_{2},

*A*

_{1},

*A*

_{2}, and

*A*, we obtain an appropriate approach to

*π*as follows:

_{i}From Eqs. (16-17), the exact solutions of the triangular points could be written as

Substituting Eqs. (46-47) and (45) into both (48-49) and limiting ourselves to the previous assumptions, it holds that

Eqs. (50-51) state that the presence of parameters that represent the oblateness of the infinitesimal body, the small perturbations in Coriolis and centrifugal forces have no effect on the *x*−coordinate of the triangular points. While the *y*−coordinate of the triangular points may be affected by all the parameters of the perturbed forces except for the small perturbations in the Coriolis force.

## 4 The stability of the libration points

The models which describe the restricted three-body problem are developed so far they provide a concise description regarding the third body motion. These models depend on the state variables that define the initial conditions and the variety of parameters that either affects the measurement process or dynamical motions, as well. There is a difficulty to solve these models directly for any choice of parameters from a given set of measurements, due to the great complexity underlying these models. This is the main reason for which our attention is focused on linearizing the dynamical system to obtain simplified expressions that could be handled more easily. To understand the possible motion of the third body in the vicinity of equilibrium points, the equations of motion should be linearized along the initial state vector. In other words, we will expand the right handside of motion equations around the equilibrium points. Consequently, the linear part in these equations is called the *variational equations*. We also assume that the variation vector is related to the initial state vector by the following expression:

where *r* ≡ (*x*, *y*),__r___{0} ≡ (*x*_{0},* y*_{0}), and Δ* r* ≡ (

*ξ*,

*η*), provided that

*ξ*and

*η*are small enough. Thus, substituting Eq. (52) into Eq. (20-21), and expanding them via a Taylor series, one obtain that

where *𝒪*[3] stands from the third and higher order terms in *ξ* and *η*. When these terms and the second order terms are omitted, then Eqs. (53-54) give us the linear variational equations and *ξ* and *η* are called the *variations*. Indeed,

where the subscripts *x* and *y* denote the second partial derivatives of Ω and superscript 0 means that these derivatives are calculated with respect to one of the five equilibrium points (*x*_{0},* y*_{0}).

Eqs. (55-56) represent the linear differential system of second order in ℝ^{2}. Hence, we transform such a system in a first order system in ℝ^{4}. To deal with, let us assume that *ξ*_{1} = *ξ*, *ξ*_{2} = *η*, *ξ*_{3} = *ξ*_{1},* ξ*_{4} = *ξ*_{2}, and let also substitute them into Eqs. (55-56). Then the system can be written as follows:

The stability of libration points is a well-known problem explored alongn the classical literature. Therefore, in this section and also in forthcoming sections, we will follow the analysis based on Szebehely computations [29].

### 4.1 Characteristic equation

Let us consider a dynamical system within four degree of freedom which is described through the following set of differential equations:

where * χ* ∊ ℝ

^{4}, and

*= (*

__g__*g*

_{1},

*g*

_{2},

*g*

_{3},

*g*

_{4}) is a vector function. This system is non-autonomous whether the vector function

*explicitly depends on time. Otherwise, it is autonomous and can be describes in the following form:*

__g__Thus, the solution of the set of equations * g*(

__χ___{0}) = 0 for the system at Eq. (62) gives us an equilibrium point

*=*

__χ__

__χ___{0}. Hence, the solution of the system in Eq. (62) at time

*t*can be written as

__χ___{0}. Further, let

*= (*

__χ__*x*

_{1},

*x*

_{2},

*x*

_{3},

*x*

_{4})

*, where*

^{S}*S*denotes the transposed vector. Hence, if

*M*is the matrix of the dynamical system, then the linearized system can be written as

It is worth mentioning that this system has a non-trivial solution if and only if

where *λ* is an eigenvalue of matrix *M* or a root of the characteristic polynomial appeared in Eq. (64). Moreover, if * e_{i}* is an eigenvector associated with

*λ*, then the relation

_{i}*M*=

__e___{i}*λ*will be achieved. Consequently, the general solution of Eq. (63) can be written as

__e___{i}The properties of stability for *t* ≥ *t*_{0} of the linearized system depend on the nature regarding the eigenvalues of the matrix *M*. According to that, next we distinguish three main cases:

- (i) If all the eigenvalues are negative real numbers, then the solution becomes stable. However, if any one of them is positive, then the solution is unstable.
- (ii) If the eigenvalues are complex numbers having negative real parts, then the equilibrium point becomes asymptotically stable. Otherwise, if some of the eigenvalues or all of them have positive real parts, then the equilibrium point becomes unstable.
- (iii) If the eigenvalues are pure imaginary numbers, then the solution is stable. Nevertheless, if there are multiple roots, then the solution contains periodic and secular terms, so the equilibrium point becomes unstable (see Szebehely [29] and Celletti (2010)).

From the previous discussion, we can write the system in Eq. (57)-(60) in the form

where

Thus, *M* is the matrix of the linear system such that the derivatives of the potential functions are be governed by the following expressions:

where

Further, the characteristic equation corresponding to Eq. (39) is as follows:

The most fundamental questions regarding the motion near libration points are those about the stability of such points. Thus, in the next subsections we will survey the results obtained regarding the stability.

### 4.2 Stability of collinear points

To investigate the stability of collinear points, we have to study the motion in the vicinity of such points. For this purpose, Eqs. (67-69) could be written at the collinear points as follows:

It is worth mentioning that Eqs. (73)-(76) show that *L _{i}* :

*i*= 1, 2, 3. Thus, as a result of this,

Let *λ*^{2} = Λ. Therefore, Eq. (72) could be written as

where

Hence, the roots of Eq. (77) are

Eqs. (78-79) and (80) show that the roots of the characteristic Eqs. (72) could be written as *λ*_{1,2} = ±σ, *λ*_{3,4} = ±*iτ*, where *σ* and *τ* are real numbers. Furthermore, it holds that

In addition, the general solution of Eqs. (53-54) can be written as

Therefore,

Hence, we can establish that the motion around the collinear points is unbounded. However, it holds that *λ*_{3,4} are pure imaginary, and *λ*_{1,2} are real. Accordingly, the solutions of these points become unstable.

### 4.3 Stability of triangular points

In this case, Eqs. (67-69) could be written as follows:

Let *λ*^{2} = *ω*. Then Eq. (72) leads to

where

where the superscript 0 refers to the second partial derivatives which have to be evaluated at the triangular points. Hence,

where

where *α*, *β* , and *γ* have to be evaluated as in Appendix 2.

The variations *ξ* and *η* will represent stable solutions in the proximity of *L*_{4} and *L*_{5}, provided that the four roots of the characteristic Eq. (72) are purely imaginary numbers. Otherwise, if any of such roots are real or complex, then these solutions will increase with time, and therefore, the solutions become unstable. Moroever, from Eq. (92), it is easy to show that

Eq. (93) establishes that *D* is a decreasing function on the interval (0, 1/2). Consequently, there is only one value of *μ*, say *μ _{c}*, on that interval for which

*D*equals zero and

*D*is positive provided that 0 <

*μ*<

*μ*. In that case, the solutions are stable.

_{c}### 4.4 Critical mass

From the discussion above, it holds that the value of *μ* for which *D* equals zero is called the *critical mass value* (see Eq. (92)). This value will be governed by

Thus, if we substitute *α, β* , and *γ* into Eq. (94) and restrict ourselves to the linear terms of *A _{i}*,

*A*,

*ε*, as well as to the coupling terms in

_{i}*A*and

_{i}ε_{i}*Aε*:

_{i}*i*= 1,2, then the critical mass value could be written as

where *p* represents these effects when the three participating bodies are oblate spheroids. Accordingly,

where *μ*_{00}, *μ*_{10}, *μ*_{i1}, *μ*_{i2}, and *μ*_{i3} are as in Appendix 3.

It is worth mentioning that Eq. (96) consists of three main parts. The first part gives the effects of small perturbations *ε*_{1} and *ε*_{2} in Coriolis and centrifugal forces, respectively. The second part represents oblateness influence of the three participating bodies, namely, *A*_{1},* A*_{2}, and *A*_{3}. Finally, the third part represents the mixed effects for only one of the parameters *A*_{1},* A*_{2},* A*, together with *ε*_{1} and *ε*_{2}. Thus, that part will disappear provided that one of the parameters is ignored.

## 5 Periodic orbits

The periodic orbits have great significance generally in the celestial mechanics and specially in space dynamics. That significance makes necessity for studying the periodic orbits. There are some important reasons that motivated this study. First of all, note that the non-integrable dynamical systems prevent us to obtain a complete solution for any orbits unless they are asymptotic or periodic orbits. For any particular solution of the restricted three-body problem, one can always find out a periodic solution. We can also get periodic orbits via linearized solutions for the motion of infinitesimal body around the libration points. Periodic orbits play a crucial role in order to classify the different classes of orbits and also to reduce the dimensionality of the problem in the phase space. Additionally, they can be used as reference orbits.

### 5.1 Periodic orbits around collinear points

It is easy to find out periodic orbits about the collinear points *L*_{1,2,3}. However, these points are unstable, i.e., if a body in any of these points is disturbed, then it will move away. This is quite reasonable from a physical point of view. In fact, for these points there are saddle points in the potential function. Thus, substituting Eq. (83) into (55-56), simple calculations give the relations between the coefficients *ρ _{i}* and

*δ*. In fact, it holds that

_{i}where

These relations show that the coefficients *δ _{i}* and

*ρ*:

_{i}*i*= 1, 2, 3, 4 are dependent. Therefore, the four initial conditions

*ξ*

_{0},

*η*

_{0},

*δ*and

_{i}*ρ*, where the subscript 0 refers to these quantities have to be evaluated at the initial time (

_{i}*t*=

*t*

_{0}). Hence, substituting Eq. (97) into both Eqs. (83) and (84), the following expressions hold:

Since the determinant (Δ) of the system appeared in Eqs. (99)-(102) is not zero, then

The coefficients in that system can be written as functions of the initial conditions. Thus, if the initial conditions *ξ*_{0} and *η*_{0} are properly chosen for *δ*_{1} = *δ*_{2} = 0, then Eq. (83 could be written in the form

Thus, from Eqs. (104-106), the following holds:

This means that once the components of the initial conditions *ξ*_{0} and *η*_{0} are chosen, then we cannot select the associated initial velocities

Eq. (109) establishes that the trajectory of the infinitesimal body around the collinear points becomes an ellipse whose center is located at these points. Further, the semi-major (*a _{c}*) and the semi-minor

*b*axes are parallel to the

_{c}*y*−axis and the

*x*−axis, respectively. Accordingly, these axes, the eccentricity (

*e*), and the period (

_{c}*T*) can be written in the following form:

_{c}Since *ξ*_{0} ≠ 0, and *η*_{0} = 0, then we conclude that the motion along the orbits is retrograde.

### 5.2 Periodic orbits around triangular points

Since the triangular points are linearly stable provided that 0 < *μ* < *μ _{c}*, and the characteristic equation presents four purely imaginary roots, then we have a bounded motion around the triangular points which consists of two harmonic motions. Consequently, this motion will be governed by

where *s*_{1} and *s*_{2} are the angular frequencies with respect to long and short periodic orbits, respectively. Moreover, the terms within the coefficients *s*_{1} and *s*_{2} will be given by

Eqs. (116-117) give the frequencies for the orbits of the long and short periodic motion when the influence of oblateness, small perturbations of Coriolis (*ε*_{1}), and centrifugal (*ε*_{2}) forces are considered. We would like also to show that these equations are valid in the range 0 < *μ* < *μ*_{c}, where *μ*_{c} is the critical mass value.

Substituting Eqs. (114-115) into Eqs. (55-56), and equaling the coefficients of sine and cosine terms, respectively, then the following relation among the coefficients of the long and short periodic terms yields:

where

and *t* = 0, i.e., the short periodic terms are removed. Thus, substituting these quantities into Eqs. (114-115) and (118-119), then the following expressions hold:

Let us assume that the triangular points represent the origin of the coordinate system and also that the infinitesimal body starts its motion at the origin. Thus, from Eqs. (50-51) and (52), it holds that the initial conditions will be controlled by (*ξ*_{0},* η*_{0}) = (−*x*_{0}, −*y*_{0}), where

where the ∓ sign means that the origin is at *L*_{4} (*L*_{5}).

### 5.3 Elliptic orbits

After we have removed the long or the short periodic terms, it holds that the path of the infinitesimal body will be an ellipse. We can justify this fact by rewriting Eqs. (114-115) with respect to the long periodic terms in the following form:

On the other hand, if we replace Eqs. (121-124) into Eqs. (129-130), and remove the time effect, then these expressions could be reduced to

This equation represents an ellipse centred at the origin of the coordinate system *ξ* and *η*, since

where

### 5.4 The orientation of principal axes of the ellipse

Since Eq. (131) includes a bilinear term of the form *ξ η*, then the principal axes of the ellipse are rotated by an angle *θ* corresponding to the coordinate system (*ξ*, *η*). This motivates the introduction of a new coordinate

Hence, let us substitute Eqs. (136-137) into Eq. (131), and equate the coefficient of

where the + sign refers to the center of the ellipse at *L*_{4}, while the − sign corresponds to the *L*_{5} case. Furthermore, it holds that the lengths of the semi-major (*a _{t}*) and the semi-minor (

*b*) axes, the eccentricity (

_{t}*e*), as well as the period of the motion (

_{t}*T*), are controlled by

_{t}Next, we will introduce an algorithm which allows to find out the elements in the periodic orbits around the equilibrium points. This algorithm will be formulated via the following steps:

- (i) Determine the parameters
*μ*,*A*_{1},*A*_{2},*A, ε*_{1}, and*ε*_{2}, for any given system. - (ii) Evaluate
, and${\Omega}_{xx}^{0},{\Omega}_{xy}^{0}$ .${\Omega}_{yy}^{0}$ - (iii) Evaluate
*τ*. - (iv) Evaluate
*s*_{1}and*s*_{2}. - (v) Evaluate
*m*_{3}. - (vi) Evaluate
*ξ*_{0}and*η*_{0}. - (vii) Evaluate
*α*_{1},*β*_{1}, and*γ*_{1}. - (viii) Evaluate sin 2
*θ*and cos 2*θ*. - (ix) Apply steps (v) and (vi) to obtain the lengths of the semi-major (
*a*) and the semi-minor (_{c}*b*) axes as well as the value of the eccentricity (_{c}*e*)._{c} - (x) Apply step (v) to find the period of motion
*T*._{c} - (xi) Use steps (vii) and (viii) to get the lengths of the semi-major (
*a*) and the semi-minor (_{t}*b*) axes as well as the value of the eccentricity (_{t}*e*)._{t} - (xii) Use step (iv) to find the period of motion
*T*._{t} - (xiii) Use steps (ix) and (x) for the periodic orbits around the collinear points.
- (xiv) Use steps (xi) and (xii) for the periodic orbits around the triangular points.

## 6 Conclusions

The restricted three-body problem is investigated in the framework the three participating are oblate spheroids, as well as the effect of small perturbations in the Coriolis and the centrifugal forces. The existence of libration points and their linear stability have been studied, under the parameter effects consisting of the oblateness and the small perturbations in the Coriolis and the centrifugal forces. The Lagrangian method is applied to invert expansion series in order to determined the locations of the collinear points by a semi-closed formula. We state that under the aforementioned effects, the collinear points remain unstable and the triangular points are stable for 0 < *μ* < *μ _{c}*, and unstable for

*μ*

_{c}≤

*μ*≤ 1/2.

We prove that the trajectories of the infinitesimal body around the five libration points are ellipses. In this context the elements of these ellipses (the frequencies, the semi-major,the semi-minor axes, the eccentricity and the period of motion) are also found. In addition, the directions of the principal axes for the orbits around the triangular points and the coefficients of long and short periodic terms have been evaluated. Furthermore, an algorithm has been provided in order to calculate the elements of the periodic orbits around such points.

^{}

Communicated by Jaume Llibre

- The series of coefficients at
*L*_{1}in Eqs. (20-21) are given as follows:$$\begin{array}{c}{a}_{4}={\displaystyle \frac{1-a-p}{b}},\text{\hspace{1em}\hspace{1em}}{a}_{5}={\displaystyle \frac{1+4a+2p}{b}}.\\ {a}_{6}=-{\displaystyle \frac{10ab+p-ap+3bp+{p}^{2}}{{b}^{2}}},\text{\hspace{1em}\hspace{1em}}{a}_{7}={\displaystyle \frac{20ab-p-4ap+4bp-2{p}^{2}}{{b}^{2}}}.\end{array}$$ - The series of coefficients at
*L*_{2}in Eq. (24) are given by:$$\begin{array}{c}{b}_{4}={\displaystyle \frac{a+p-1}{b}},\text{\hspace{1em}\hspace{1em}}{b}_{5}={\displaystyle \frac{1+4a+2p}{b}}.\\ {b}_{6}={\displaystyle \frac{10ab+p-ap+3bp+{p}^{2}}{{b}^{2}}},\text{\hspace{1em}\hspace{1em}}{b}_{7}={\displaystyle \frac{20ab-p-4ap+4bp-2{p}^{2}}{{b}^{2}}}.\end{array}$$ - The series of coefficients at
*L*_{3}in Eq. (29) are as follows:$${c}_{0}={\displaystyle \frac{16(a+p-1)}{32-b-4p}},\text{\hspace{1em}\hspace{1em}}{c}_{1}=-{\displaystyle \frac{16(16+144a-3b-2ab-12ap+72p-4{p}^{2})}{{(32-b-4p)}^{2}}}.$$ $$\begin{array}{c}{c}_{2}={\displaystyle \frac{8}{{(32-b-4p)}^{3}}}(512+25088a-192b-608ab-7{b}^{2}+7a{b}^{2}+8384p-4160ap+66bp\\ \phantom{\rule[-0.0ex]{5em}{0.0ex}}+70abp-40{p}^{2}+200a{p}^{2}+{b}^{2}p-832{p}^{2}+6b{p}^{2}+40{p}^{3}).\end{array}$$

## References

- [1]
Abouelmagd, E.I. (2012), Existence and stability of triangular points in the restricted three-body problem with numerical applications, Astrophys. Space Sci. 342: 45–53.

- [2]↑
Abouelmagd, E.I., El-Shaboury, S.M. (2012), Periodic orbits under combined effects of oblateness and radiation in the restricted problem of three bodies, Astrophys. Space Sci. 341: 331–341.

- [3]↑
Abouelmagd, E.I., Sharaf, M.A. (2013), The motion around the libration points in the restricted three-body problem with the effect of radiation and oblateness, Astrophys. Space Sci. 344: 321–332.

- [4]
Abouelmagd, E.I., Asiri, H.M., Sharaf, M.A. (2013), The effect of oblateness in the perturbed restricted three-body problem, Meccanica 48: 2479–2490.

- [5]↑
Abouelmagd, E.I. (2013), Stability of the triangular points under combined effects of radiation and oblateness in the restricted three-body problem, Earth Moon and Planets 110: 143–155.

- [6]↑
Abouelmagd, E.I. (2013), The effect of photogravitational force and oblateness in the perturbed restricted three-body problem, Astrophysics Space Science 346: 51–69.

- [7]↑
Abouelmagd, E.I., Guirao, Juan L.G., Mostafa, A. (2014), Numerical integration of the restricted thee-body problem with Lie series, Astrophysics Space Science 354: 369–378.

- [8]↑
Abouelmagd, E.I., Awad, M.E., Elzayat, E.M.A., Abbas, I.A. (2014), Reduction the secular solution to periodic solution in the generalized restricted three-body problem, Astrophysics Space Science 350: 495–505.

- [9]↑
Abouelmagd, E.I., Alhothuali, M.S., Guirao, Juan L.G., Malaikah, H.M. (2015), The effect of zonal harmonic coefficients in the framework of the restricted three body problem, Advances in Space Research 55: 1660–1672.

- [10]
Abouelmagd, E.I., Alhothuali, M.S., Guirao, Juan L.G., Malaikah, H.M. (2015), Periodic and secular solutions in the restricted three-body problem under the effect of zonal harmonic parameters, Applied Mathematics & Information Science 9 (4): 1–11.

- [11]↑
Abouelmagd, E.I., Guirao, Juan L.G., Vera, Juan A. (2015), Dynamics of a dumbbell satellite under the zonal harmonic effect of an oblate body, Communications in Nonlinear Science and Numerical Simulation 20: 1057–1069.

- [12]↑
Bhatnagar, K.B., Hallan, P.P. (1978), Effect of perturbations in Coriolis and centrifugal forces on the stability of libration points in the restricted problem, Celest. Mech. 18: 105–112.

- [13]
Bhatnagar, K.B., Hallan, P.P. (1979), Effect of perturbed potentials on the stability of libration points in the restricted problem, Celest. Mech. 20: 95–103.

- [14]↑
Devi, G.S., Singh, R. (1994), Location of equilibrium points in the perturbed photogravitational circular restricted problem of three bodies, Bull Astr Soc India 22: 433–437.

- [15]↑
Elipe, A., Ferrer, S. (1985), On the equilibrium solution in the circular planar restricted three rigid bodies problem, Celest. Mech. 37: 59–70.

- [16]
El-Shaboury, S.M. (1989), The equilibrium solutions of restricted problem of three axisymmetric rigid bodies, Earth, Moon and Planets 45: 205–211.

- [17]↑
El-Shaboury, S.M., Shaker, M.O., El-Dessoky, A.E., El Tantawy, M.A. (1991), The Libration points of axisym-metric satellite in the gravitational field of two triaxial rigid body, Earth, Moon and Planets 52: 69–81.

- [18]↑
El-Shaboury, S.M., El-Tantawy, M.A. (1989), Eulerian libration points of restricted problem of three oblate spheroids, Earth, Moon and Planets 63: 23–28.

- [19]↑
Khanna, M., Bhatnagar, K.B. (1999), Existence and stability of libration points in the restricted three body problem when the smaller primary is a triaxial rigid body and the bigger one an oblate spheroid, Indian J. Pure Appl. Math. 30 (7): 721–733.

- [20]↑
Markellos, V.V., Papadakis, K.E., Perdios, E.A. (1996), Non-linear stability zones around triangular equilibria in the plane circular restricted three-body problem with oblateness, Astrophys Space Sci. 245: 157–164.

- [21]↑
Mittal, A., Ahmad, I., Bhatnagar, K.B. (2009), Periodic orbits generated by Lagrangian solution of the restricted three body problem when one of the primaries is an oblate body, Astrophys Space Sci. 319: 63–73.

- [22]↑
Murray, C.D., Dermott, S.F. (1999), Solar system dynamics, Cambridge University Press, Cambridge.

- [23]↑
Sharma, R.K., Subbarao, P.V. (1978), A case of commensurability induced by oblateness, Celest. Mech. 18: 185–194.

- [24]↑
Sharma, R.K., Taqvi, Z.A., Bhatnagar, K.B. (2001), Existence of libration Points in the restricted three body problem when both primaries are triaxial rigid bodies, Indian J. Pure Appl. Math. 32 (1): 125–141.

- [25]↑
Shu, Sh., Lu, Bk. (2005), Effect of perturbation of Coriolis and Centrifugal forces on the location and linear stability of the libration points in the Robe problem, ChA&A 29: 421–429.

- [26]↑
Singh, J. (2009), Effect of perturbations on the nonlinear stability of triangular points in the restricted three-body problem with variable mass, Astrophys. Space Sci. 321: 127–135.

- [27]↑
Singh, J., Begha, J.M. (2011), Periodic orbits in the generalized perturbed restricted three-body problem, Astrophys. Space Sci. 332: 319–324.

- [28]↑
Subbarao, P.V., Sharma, R.K., (1975), A note on the stability of the triangular points of equilibrium in the restricted three-body problem, Astron. & Astrophys 43: 381–383.

- [29]↑
Szebehely, V. (1967), Stability of the points of equilibrium in the restricted problem, Astron. J. 72: 7–9.