Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem with Oblateness

Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem... Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem with Oblateness Hindawi Publishing Corporation Home Journals About Us About this Journal Submit a Manuscript Table of Contents Journal Menu Abstracting and Indexing Aims and Scope Article Processing Charges Articles in Press Author Guidelines Bibliographic Information Contact Information Editorial Board Editorial Workflow Free eTOC Alerts Reviewers Acknowledgment Subscription Information Open Special Issues Published Special Issues Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Full-Text ePUB Linked References How to Cite this Article Complete Special Issue Advances in Mathematical Physics Volume 2012 (2012), Article ID 679063, 18 pages doi:10.1155/2012/679063 Research Article Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem with Oblateness Jagadish Singh 1 and Abubakar Umar Sandah 2 1 Department of Mathematics, Faculty of Science, Ahmadu Bello University, Zaria, Nigeria 2 Department of Mathematics Statistics and Computer Science, College of Science and Technology, Kaduna Polytechnic, Kaduna, Nigeria Received 26 March 2012; Accepted 2 July 2012 Academic Editor: Burak Polat Copyright © 2012 Jagadish Singh and Abubakar Umar Sandah. This is an open access article distributed under the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract This paper investigates the positions and linear stability of an infinitesimal body around the equilibrium points in the framework of the Robe’s circular restricted three-body problem, with assumptions that the hydrostatic equilibrium figure of the first primary is an oblate spheroid and the second primary is an oblate body as well. It is found that equilibrium point exists near the centre of the first primary. Further, there can be one more equilibrium point on the line joining the centers of both primaries. Points on the circle within the first primary are also equilibrium points under certain conditions and the existence of two out-of-plane points is also observed. The linear stability of this configuration is examined and it is found that points near the center of the first primary are conditionally stable, while the circular and out of plane equilibrium points are unstable. 1. Introduction Robe [ 1 ] considered a new kind of restricted three-body problem in which, one of the primaries of mass 𝑚 1 is a rigid spherical shell, filled with homogenous, incompressible fluid of density 𝜌 1 ; the second one is a point mass 𝑚 2 located outside the shell and moving around the mass 𝑚 1 in a Keplerian orbit; the infinitesimal mass 𝑚 3 is a small sphere of density 𝜌 3 , moving inside the shell and is subject to the attraction of 𝑚 2 and the buoyancy force due to the fluid of the first primary. Further, he discussed the linear stability of an equilibrium point obtained in two cases. In the first case, the orbit of 𝑚 2 around 𝑚 1 is circular and in the second case, the orbit is elliptic, but the shell is empty (there is no fluid inside it) or densities of 𝑚 1 and 𝑚 3 are equal. Since then various studies (e.g., [ 2 – 4 ]) under different assumptions have been carried out. In his study, Robe [ 1 ] assumed that the pressure field of the fluid 𝜌 1 has a spherical symmetry around the center of the shell and he took into account only one out of the three components of the pressure field which is due to the own gravitational field of the fluid 𝜌 1 . He did not consider the other two components arising from the attraction of 𝑚 2 and the centrifugal force. Taking care of all these three components of the pressure field, A. R. Plastino and A. Plastino [ 5 ] reanalyzed the Robe’s. But in their study, they assumed the hydrostatic equilibrium figure of the first primary as Roche’s ellipsoid (see Figure 1 ). They found that when the density parameter 𝐷 is taken as zero, every point inside the fluid is an equilibrium point; otherwise the center of the ellipsoid is the only equilibrium point and it is linearly stable. Figure 1: The Robe’s CRTBP with oblate primaries. Hallan and Rana [ 3 ] investigated the existence of all equilibrium point and their stability in the Robe’s [ 1 ] restricted three-body problem. It was seen that the Robe’s elliptic restricted three-body problem has only one equilibrium point for all values of the density parameter 𝐾 and the mass parameter 𝜇 , while the Robe’s circular restricted three-body problem can have two, three, or infinite numbers of equilibrium points. As regards to the stability of these equilibria, they confirmed the stability result given by Robe [ 1 ] of the equilibrium point ( − 𝜇 , 0 , 0 ) , whereas triangular and circular points are always unstable. The equilibrium point collinear with the center of the shell and the second primary was found to be stable under some conditions. Hallan and Mangang [ 4 ] studied the Robe’s [ 1 ] restricted three-body problem by considering the full buoyancy force as in A. R. Plastino and A. Plastino [ 5 ] and assuming the hydrostatic equilibrium figure of the first primary as an oblate spheroid. They derived the pertinent equations of motion and discussed the existence of equilibrium point and their linear stability. The participating bodies in the classical restricted three-body problem are strictly spherical in shape, but in actual situations several heavenly bodies, such as Saturn and Jupiter, are sufficiently oblate. The minor planets and meteoroids have irregular shape. The lack of sphericity, or the oblateness, of the planet causes large perturbations from a two-body orbit. The motions of artificial Earth satellites are examples of this. Global studies of problems with oblateness have been carried out by many researchers (e.g., [ 6 – 9 ]). Therefore, our effort in this paper aims at investigating the equilibrium points and their stability in the Robe’s circular restricted three-body problem when the hydrostatic equilibrium figure of the fluid of the first primary is an oblate spheroid and the second one is an oblate spheroid as well. The model of this study can be used to study the small oscillation of the Earth’s inner core taking into account the Moon’s attraction. This paper is organized as follows; Section 2 represents the equations of motion; the existence of the equilibrium points is mentioned in Section 3 , while Section 4 investigates their linear stability; Section 5 discusses the results obtained; the conclusion is drawn in Section 6 . 2. Equation of Motion Let the first primary 𝑚 1 be a fluid of density 𝜌 1 in the shape of an oblate spheroid as assumed by Hallan and Mangang [ 4 ]; let the second primary 𝑚 2 be an oblate body too as Sharma and Subba Rao [ 6 ] assumed, which describes a circular orbit around 𝑚 1 . We adopt a uniformly rotating coordinate system 𝑂 𝑥 1 𝑥 2 𝑥 3 with origin at the center of mass 𝑚 1 , 𝑂 𝑥 1 pointing towards 𝑚 2 , with 𝑂 𝑥 1 𝑥 2 being the orbital plane of 𝑚 2 coinciding with the equatorial plane of 𝑚 1 . Then, the equations of motion of the infinitesimal body of density 𝜌 3 in the coordinate system take the form [ 4 , 6 ]: ̈ 𝑥 1 − 2 𝑛 ̇ 𝑥 2 = 𝜕 𝑈 𝜕 𝑥 1 , ̈ 𝑥 2 + 2 𝑛 ̇ 𝑥 1 = 𝜕 𝑈 𝜕 𝑥 2 , ̈ 𝑥 3 = 𝜕 𝑈 𝜕 𝑥 3 , ( 2 . 1 ) where 𝑛 𝑈 = 𝑉 + 2   𝑥 1 −  𝑚 2 /  𝑚 1 + 𝑚 2 𝑅    2 + 𝑥 2 2  2 , 𝑉 = 𝐵 + 𝐵  − 𝜌 1 𝜌 3 ⎡ ⎢ ⎢ ⎣ 𝐵 + 𝐵  + 𝑛 2   𝑥 1 −  𝑚 2 /  𝑚 1 + 𝑚 2 𝑅    2 + 𝑥 2 2  2 ⎤ ⎥ ⎥ ⎦ , 𝐵 = 𝜋 𝐺 𝜌 1  𝐼 − 𝐴 1 𝑥 2 1 − 𝐴 1 𝑥 2 2 − 𝐴 2 𝑥 2 3  , 𝐵  = 𝐺 𝑚 2   𝑅 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  1 / 2 + 𝐺 𝑚 2 𝛼 2 2   𝑅 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝐺 𝑚 2 𝛼 2 𝑥 2 3 2   𝑅 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 , 𝐼 = 2 𝑎 2 1 𝐴 1 + 𝑎 2 2 𝐴 2 , 𝐴 1 = 𝑎 2 1 𝑎 2  ∞ 0 𝑑 𝑢 Δ  𝑎 2 1  + 𝑢 , 𝐴 2 = 𝑎 2 1 𝑎 2  ∞ 0 𝑑 𝑢 Δ  𝑎 2 2  , Δ + 𝑢 2 =  𝑎 2 1  + 𝑢 2  𝑎 2 2  , 𝑛 + 𝑢 2 = 𝐺  𝑚 1 + 𝑚 2  𝑅 2  3 1 + 2 𝛼 1 + 3 2 𝛼 2  ; 𝛼 1 = 𝑎 2 1 − 𝑎 2 3 5 𝑅 2 , 𝛼 2 = 𝑎 2 2 − 𝑎 2 4 5 𝑅 2 , 𝜌 𝐷 = 1 − 1 𝜌 3 . ( 2 . 2 ) Here 𝑉 is the potential that explains the combined forces upon the infinitesimal mass, 𝐵 denotes the potential due to the fluid mass of the first primary, 𝐵 ′ stands for the potential due to the second primary, 𝑅 is the distance between the primaries, and 𝐺 is the gravitational constant. 𝑛 is the mean motion. 𝑎 1 , 𝑎 2 and 𝑎 3 , 𝑎 4 are the equatorial and polar radii of the first and second primary, respectively. 𝐼 stands for the polar moment of inertia, while 𝐴 𝑖 ( 𝑖 = 1 , 2 ) are the index symbols. 𝛼 1 and 𝛼 2 are the oblateness coefficients of the first and second primaries, respectively. We choose the unit of mass such that the sum of the masses of the primaries is taken as unity, thus we take 𝑚 2 = 𝜇 , 0 < 𝜇 = 𝑚 2 / ( 𝑚 1 + 𝑚 2 ) < 1 . For the unit of length, we take the distance between the primaries as unity, that is, 𝑅 = 1 and the unit of time is also selected such that 𝐺 = 1 . With these units and substituting the expression for the potential 𝐵 due to the fluid in the first primary and the potential 𝐵 ′ due to the second oblate primary, the equations of motion ( 2.1 ) are recast to the form: ̈ 𝑥 1 − 2 𝑛 ̇ 𝑥 2 = 𝜕 𝑈 𝜕 𝑥 1 , ̈ 𝑥 2 + 2 𝑛 ̇ 𝑥 1 = 𝜕 𝑈 𝜕 𝑥 2 , ̈ 𝑥 3 = 𝜕 𝑈 𝜕 𝑥 3 , ( 2 . 3 ) where ⎡ ⎢ ⎢ ⎢ ⎣ 𝑈 = 𝐷 𝜋 𝜌 1  𝐼 − 𝐴 1  𝑥 2 1 + 𝑥 2 2  − 𝐴 2 𝑥 2 3  + 𝜇   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  1 / 2 + 𝜇 𝛼 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 + 𝑛 2   𝑥 1  − 𝜇 2 + 𝑥 2 2  2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑛 2 =  3 1 + 2 𝛼 1 + 3 2 𝛼 2  . ( 2 . 4 ) These above equations of motion of the infinitesimal mass 𝑚 3 under the framework of the Robe’s circular restricted three-body problem have been obtained by taking into account the shapes of the primaries, the full buoyancy force, the forces due to the gravitational attraction of the second primary, and the gravitational force exerted by the fluid of density 𝜌 1 . In the case when the second primary is not an oblate spheroid (i.e., 𝛼 2 = 0 ), the equations are the same as those of Hallan and Mangang [ 4 ]. 3. Position of Equilibrium Points The equilibrium points are the solutions of the equations: 𝑈 𝑥 1 = 𝑈 𝑥 2 = 𝑈 𝑥 3 = 0 . ( 3 . 1 ) That is, 𝑈 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝑥 1 𝐴 1 + 𝜇  1 − 𝑥 1    1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 + 3 𝜇 𝛼 2  1 − 𝑥 1  2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 − 1 5 𝜇 𝛼 2  1 − 𝑥 1  𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  7 / 2 + 𝑛 2  𝑥 1  ⎤ ⎥ ⎥ ⎥ ⎦ 𝑈 − 𝜇 = 0 , ( 3 . 2 ) 𝑥 2 = 𝐷 𝑥 2 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  7 / 2 + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ 𝑈 = 0 , ( 3 . 3 ) 𝑥 3 = 𝐷 𝑥 3 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 2 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 − 3 𝜇 𝛼 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  7 / 2 ⎤ ⎥ ⎥ ⎥ ⎦ = 0 . ( 3 . 4 ) 3.1. Equilibrium Points Near the Centre of the First Primary The positions of the equilibrium points near the first primary are the solutions of ( 3.2 ) when 𝑈 𝑥 1 = 0 , 𝑥 1 ≠ 0 , 𝑥 2 = 𝑥 3 = 0 , 𝐷 ≠ 0 , and 𝑛 2 = 1 + ( 3 / 2 ) ( 𝛼 1 + 𝛼 2 ) . The 𝑥 1 coordinate of the equilibrium points are then the roots of the equation: − 2 𝜋 𝜌 1 𝑥 1 𝐴 1 + 𝜇  1 − 𝑥 1  | | 1 − 𝑥 1 | | 3 + 3 𝜇 𝛼 2  1 − 𝑥 1  2 | | 1 − 𝑥 1 | | 5 +  3 1 + 2 𝛼 1 + 3 2 𝛼 2   𝑥 1  − 𝜇 = 0 . ( 3 . 5 ) We first determine the roots of ( 3.5 ) in the absence of oblateness, that is, the case when the primaries are spherical. In this case, the roots are [ 4 ] 𝑥 1 1 √ = 1 + 𝜇 + 𝜇 2 + 8 𝜇 𝜋 𝜌 1 𝐴 1 − 4 𝜇 2  1 − 2 𝜋 𝜌 1 𝐴 1  , 𝑥 1 2 √ = 1 + 𝜇 − 𝜇 2 + 8 𝜇 𝜋 𝜌 1 𝐴 1 − 4 𝜇 2  1 − 2 𝜋 𝜌 1 𝐴 1  . ( 3 . 6 ) The term 𝐴 1 which appears in ( 3.2 ) and is due to the fluid mass affects these roots. Therefore, these roots will be real if the discriminant is nonnegative, that is if 𝜇 + 8 𝜋 𝜌 1 𝐴 1 − 4 ≥ 0 . ( 3 . 7 ) When ( 1 / 4 ) 𝜇 ≥ 1 − 2 𝜋 𝜌 1 𝐴 1 > 0 , both roots are greater than unity and we reject them because they lie outside the first primary. Now, if 1 − 2 𝜋 𝜌 1 𝐴 1 < 0 , we have 𝑥 1 2 > 1 and 𝑥 1 1 < 1 . Further, we see that 𝑥 1 1 > − 1 when 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 . Thus, in the case when 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 , the point ( 𝑥 1 1 , 0 , 0 ) lies within the first primary if | 𝑥 1 1 | < 𝑎 1 . When 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 , | 𝑥 1 1 | < 𝑎 1 ; 𝑥 1 1 is a root of ( 3.5 ). Hence for 1 − 2 𝜋 𝜌 1 𝐴 1 = 0 , the only root is 𝑥 1 1 = 2 which lies outside the first primary and we neglect it. Hence, for 𝛼 1 = 0 , 𝛼 2 = 0 , 𝑥 1 1 = 0 is always a root of ( 3.5 ) and 𝑥 1 = 𝑥 1 1 is also a root provided 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 , | 𝑥 1 1 | < 𝑎 1 . Now, we find the roots of ( 3.5 ) when oblateness of both primaries is considered (i.e., 𝛼 1 ≠ 0 , 𝛼 2 ≠ 0 ). Let the roots be such that 𝑥 1 = 0 + 𝑝 1 , | | 𝑝 1 | | 𝑥 ≪ 1 , 1 = 𝑥 1 1 + 𝑝 2 , | | 𝑝 2 | | ≪ 1 . ( 3 . 8 ) Putting these values in ( 3.5 ), multiplying throughout by ( 1 − 𝑝 1 ) 4 , expanding and neglecting second and higher powers of 𝑝 1 , 𝛼 1 , 𝛼 2 , as they are very small quantities, we have 𝑝 1 ≅ − 3 𝛼 1 2 𝜇 2 𝜋 𝜌 1 𝐴 1 . − ( 1 + 2 𝜇 ) ( 3 . 9 ) Similarly, putting 𝑥 1 = 𝑥 1 1 + 𝑝 2 in ( 3.5 ) and then simplifying it, we get  1 − 𝑥 1 1  4   𝑥 1 1   3 − 𝜇 1 + 2 𝛼 1 + 3 2 𝛼 2  − 2 𝜋 𝜌 1 𝐴 1 𝑥 1 1  + 𝑝 2  1 − 𝑥 1 1  3   1 − 3 𝑥 1 1   1 − 2 𝜋 𝜌 1 𝐴 1    + 2 𝜇 + 𝜇 1 − 𝑥 1 1  2 − 2 𝑝 2  1 − 𝑥 1 1    1 − 𝑥 1 1  2  𝑥 1 1 − 2 𝜋 𝜌 1 𝐴 1 𝑥 1 1   − 𝜇 + 𝜇 = − 3 𝜇 𝛼 2 2 . ( 3 . 1 0 ) Multiplying ( 3.8 ) by ( 1 − 𝑥 1 1 ) 2 , simplifying and then using it in ( 3.10 ), we get 𝑝 2  ≅ − 1 − 𝑥 1 1 𝑥    1 1 − 𝜇   ( 3 / 2 ) 𝛼 1 + ( 3 / 2 ) 𝛼 2     1 − 3 𝑥 1 1   1 − 2 𝜋 𝜌 1 𝐴 1   − + 2 𝜇 3 𝜇 𝛼 2 2  1 − 𝑥 1 1  3   1 − 3 𝑥 1 1   1 − 2 𝜋 𝜌 1 𝐴 1   . + 2 𝜇 ( 3 . 1 1 ) A substitution of ( 3.11 ) in the second equation of ( 3.8 ) at once gives the position of the other equilibrium point near the center of the first primary. 3.1.1. Positions of Circular Points The positions of the circular points are sought using the first two equations of system ( 3.1 ) with the conditions 𝑥 1 ≠ 0 , 𝑥 2 ≠ 0 , 𝑥 3 = 0 ; that is, they are the solutions of 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2 3 + 1 + 2 𝛼 1 + 3 2 𝛼 2 ⎤ ⎥ ⎥ ⎥ ⎦ + 𝜇   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 + 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2  3 − 𝜇 1 + 2 𝛼 1 + 3 2 𝛼 2  = 0 , ( 3 . 1 2 ) − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2 3 + 1 + 2 𝛼 1 + 3 2 𝛼 2 = 0 . ( 3 . 1 3 ) Solving the above equations and knowing that 𝜇 ≠ 0 , we get 1   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 + 3 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2 − 𝑛 2 = 0 . ( 3 . 1 4 ) We let  1 − 𝑥 1  2 + 𝑥 2 2 = 𝑟 2 . ( 3 . 1 5 ) Substituting ( 3.15 ) in ( 3.14 ), and simplifying, we get 𝑛 2 𝑟 5 − 𝑟 2 − 3 2 𝛼 2 = 0 . ( 3 . 1 6 ) Now, we let 𝑟 = 1 + 𝜀 , 𝜀 ≪ 1 . ( 3 . 1 7 ) Substituting ( 3.16 ) in ( 3.15 ), neglecting second and higher powers of 𝜀 , we get 1 𝜀 = − 2 𝛼 1 . ( 3 . 1 8 ) Therefore, ( 3.17 ) is now expressed as 1 𝑟 ≅ 1 − 2 𝛼 1 . ( 3 . 1 9 ) A substitution of ( 3.14 ) in ( 3.13 ) yields 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) . ( 3 . 2 0 ) Therefore, when 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) , the points on the circle given by ( 3.15 ) with 𝑥 3 = 0 and 𝑟 = 1 − ( 1 / 2 ) 𝛼 1 lying within the first primary are also equilibrium points. The general coordinates of these circular points are given by ( 1 + 𝑟 c o s 𝜃 , 𝑟 s i n 𝜃 , 0 ) , where 𝜃 is a parameter. When 𝑦 = 0 , the circular points coalesce to those lying on the line joining the primaries. 3.1.2. Positions of Out-of-Plane Equilibrium Points The out-of-plane points have no analogy in the classical restricted three-body problem. However the investigation concerning these points in the photogravitational restricted three-body problem was first carried out by Radzievskii [ 10 ]. Afterwards, other researchers, for instance Douskos and Markellos [ 8 ], Singh and Leke [ 11 ], and so forth, have worked on the out-of-plane points. In this section, we locate these points for our study, as it has remained an open problem to date. The positions of the out-of-plane equilibrium points of the Robe’s problem with oblate primaries are the solutions of the first and last equations of ( 3.1 ) with 𝑥 2 = 0 , 𝐷 ≠ 0 ; that is, 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ + 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 + 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 − 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 − 𝑛 2 𝑥 𝜇 = 0 , ( 3 . 2 1 ) 3 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 2 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 − 9 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎤ ⎥ ⎥ ⎥ ⎦ = 0 . ( 3 . 2 2 ) From ( 3.22 ), since 𝑥 3 ≠ 0 , we have − 2 𝜋 𝜌 1 𝐴 2 = 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 + 9 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 − 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 . ( 3 . 2 3 ) Let 𝑙 2 =  1 − 𝑥 1  2 + 𝑥 2 3 . ( 3 . 2 4 ) Then, ( 3.23 ) and ( 3.21 ) may be written respectively: 1 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 = 2 𝜋 𝜌 1 𝐴 2 + 𝜇 𝑙 3 + 9 𝜇 𝛼 2 2 𝑙 5 , 𝑥 1  − 2 𝜋 𝜌 1 𝐴 1 − 𝜇 𝑙 3 − 3 𝜇 𝛼 2 2 𝑙 5 + 1 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 + 𝑛 2  + 𝜇 𝑙 3 + 3 𝜇 𝛼 2 2 𝑙 5 − 1 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 − 𝑛 2 𝜇 = 0 . ( 3 . 2 5 ) Now, from first equations ( 3.25 ), we get 𝑥 3 𝑙 = ± √ 1 5 𝜇 𝛼 2  𝜇  9 𝛼 2 + 2 𝑙 2  + 4 𝜋 𝜌 1 𝐴 2 𝑙 5  1 / 2 . ( 3 . 2 6 ) The use of ( 3.24 ) in second equation of ( 3.25 ) yields 𝑥 1 =  2 𝜋 𝜌 1 𝐴 2 + 𝑛 2 𝜇  𝑙 5 + 3 𝜇 𝛼 2  2 𝜋 𝜌 1  𝐴 2 − 𝐴 1  + 𝑛 2  𝑙 5 + 3 𝜇 𝛼 2 . ( 3 . 2 7 ) We use the software package Mathematica (Wolfram 2004) to compute the coordinates of the out-of-plane equilibrium points denoted by 𝐿 6 and 𝐿 7 starting with the initial values 𝑥 1 = 1 − 𝜇 and 𝑥 3 = √ 3 √ 𝛼 2 in the case where we have kept up to first order terms in both the numerator and the denominator; we then get 𝑥 1 = 2 𝜇 + 3 𝜇 𝛼 1 + 4 𝐴 2 𝜋 𝜌 1 2 + 3 𝛼 1 − 4 𝐴 1 𝜋 𝜌 1 + 4 𝐴 2 𝜋 𝜌 1 − 3  𝜇 𝐴 1 + ( 1 − 𝜇 ) 𝐴 2  𝛼 2 𝜋 𝜌 1  1 + 3 𝛼 1  1 + 2 𝐴 2 𝜋 𝜌 1 − 2 𝐴 1 𝜋 𝜌 1  + 4 𝜋 𝜌 1  𝐴 2 − 𝐴 1 + 𝐴 2 1 𝜋 𝜌 1 + 𝐴 2 2 𝜋 𝜌 1 − 2 𝐴 1 𝐴 2 𝜋 𝜌 1 , 𝑥   3 = √  2 / 1 5 𝜇 𝜇 3  1 + 2 𝜋 𝐴 2 𝜇 2 𝜌 1  √ 𝜇 𝛼 2 + 7 √  3 𝜇 / 1 0 𝜇 + 2 𝜋 𝐴 2 𝜇 3 𝜌 1  √ 𝛼 2 2  𝜇 3  1 + 2 𝜋 𝐴 2 𝜇 2 𝜌 1  . ( 3 . 2 8 ) The location of the out-of-plane equilibrium points can be obtained by solving numerically equations ( 3.26 ) and ( 3.27 ) using ( 3.24 ). Now, from the expression for the density parameter  𝜌 𝐷 = 1 − 1 𝜌 3  . ( 3 . 2 9 ) We assume that 𝜌 1 ≠ 𝜌 3 , then 𝐷 > 𝑜 𝑟 < 0 . In the case when the density parameter is positive, we have 𝜌 1 < 𝜌 3 . ( 3 . 3 0 ) Hence, numerically we choose 𝑎 2 1 = 0 . 9 4 , 𝑎 2 2 = 0 . 9 , 𝑎 2 3 = 0 . 8 2 , 𝑎 2 4 = 0 . 8 , 𝜇 = 0 . 0 1 , 𝜋 = 3 . 1 4 . ( 3 . 3 1 ) Then, 𝛼 1 = 0 . 0 2 4 , 𝛼 2 = 0 . 0 2 , 𝜌 1 = 0 . 2 3 6 . ( 3 . 3 2 ) Now, we perform a numerical exploration of computing the out-of-plane points in the case of the Earth-Moon system. To do this, we arbitrarily choose values for the 𝐴 𝑖 ( 𝑖 = 1 , 2 ) . We found that when 𝐴 1 = 2 . 5 0 7 6 and 𝐴 2 = 2 . 5 5 5 , the positions of the out-of-plane points ( 𝑥 1 , 0 , ± 𝑥 3 ) : 𝑥 1 = 3 . 3 5 2 7 , 𝑥 3 = 0 . 2 7 1 4 1 8 . ( 3 . 3 3 ) The abscissae of the out-of-plane point is outside the possible region of motion of the infinitesimal mass and so we neglect it. However, in the case when the 𝐴 𝑖 ( 𝑖 = 1 , 2 ) are chosen such that | | 𝐴 1 − 𝐴 2 | | ≪ 1 , 𝐴 1 = 0 . 7 , 𝐴 2 = 0 . 6 8 ( s a y ) , ( 3 . 3 4 ) 𝐴 𝑖 ∈ ( 0 , 0 . 7 ] , the point 𝐿 6 , and 𝐿 7 are, respectively, 𝑥 1 = 0 . 9 8 6 1 1 , 𝑥 3 = 0 . 2 7 1 3 8 1 ( 3 . 3 5 ) and lies within the fluid. 4. Linear Stability of the Equilibrium Points In order to study the linear stability of any equilibrium point ( 𝑥 1 0 , 𝑥 2 0 , 𝑥 3 0 ) of an infinitesimal body, we displace it to the position ( 𝑥 1 , 𝑥 2 , 𝑥 3 ) such that  𝑥 1 0 + 𝜉 , 𝑥 2 0 + 𝜂 , 𝑥 3 0  + 𝜁 , ( 4 . 1 ) where 𝜉 , 𝜂 , 𝜁 are small displacements, and then linearize equation ( 2.3 ) to obtain the equations: ̈  𝑈 𝜉 − 2 𝑛 ̇ 𝜂 = 0 𝑥 1 𝑥 1   𝑈 𝜉 + 0 𝑥 1 𝑥 2   𝑈 𝜂 + 0 𝑥 1 𝑥 3  ̇  𝑈 𝜁 , ̈ 𝜂 + 2 𝑛 𝜉 = 0 𝑥 1 𝑥 2   𝑈 𝜉 + 0 𝑥 2 𝑥 2   𝑈 𝜂 + 0 𝑥 2 𝑥 3  ̈  𝑈 𝜁 , 𝜁 = 0 𝑥 1 𝑥 3   𝑈 𝜉 + 0 𝑥 2 𝑥 3   𝑈 𝜂 + 0 𝑥 3 𝑥 3  𝜁 , ( 4 . 2 ) where the partial derivatives are evaluated at the equilibrium points. 4.1. Equilibrium Points Near the Center of the First Primary In order to consider the motion near any equilibrium point in the 𝑥 1 𝑥 2 -plane, we let solutions of the first two equations of ( 4.2 ) be 𝜉 = 𝐴 e x p ( 𝜆 𝑡 ) , 𝜂 = 𝐵 e x p ( 𝜆 𝑡 ) , ( 4 . 3 ) where 𝐴 , 𝐵 , and 𝜆 are constants. Taking first and second derivatives of the above, substituting them into the first two equations of system ( 4.2 ) and has a non-zero solution when | | | | | | |  𝜆 2 − 𝑈 0 𝑥 1 𝑥 1   2 𝑛 𝜆 + 𝑈 0 𝑥 1 𝑥 2   2 𝑛 𝜆 − 𝑈 0 𝑥 1 𝑥 2 𝜆   2 − 𝑈 0 𝑥 2 𝑥 2  | | | | | | | = 0 . ( 4 . 4 ) Expanding the determinant, we have 𝜆 4 −  𝑈 0 𝑥 1 𝑥 1 + 𝑈 0 𝑥 2 𝑥 2 − 4 𝑛 2  𝜆 2 + 𝑈 0 𝑥 1 𝑥 1 𝑈 0 𝑥 2 𝑥 2 −  𝑈 0 𝑥 1 𝑥 2  2 = 0 . ( 4 . 5 ) Equation ( 4.5 ) is the characteristic equation corresponding to the variational equations ( 4.2 ) in the case when motion is considered in the 𝑥 1 , 𝑥 2 -plane. Now, the values of the second-order partial derivatives of the equilibrium point ( 𝑥 𝐿 , 0 , 0 ) , where 𝑥 𝐿 = 𝑝 1 stands for the first equilibrium and 𝑥 𝐿 = 𝑥 1 1 + 𝑝 2 for the second one, are given as 𝑈 0 𝑥 1 𝑥 1  −  = 𝐷 𝜇 1 − 𝑥 𝐿  3 −  3 𝛼 2 / 2   1 − 𝑥 𝐿  + 2 𝑥 𝐿  1 − 𝑥 𝐿  2 + 6 𝑥 𝐿 𝛼 2 + 𝑛 2  1 − 𝑥 𝐿  5 𝑥 𝐿  1 − 𝑥 𝐿  5  , 𝑈 0 𝑥 2 𝑥 2  −  = 𝐷 𝜇 1 − 𝑥 𝐿  3 − ( 3 / 2 ) 𝛼 2  1 − 𝑥 𝐿  − 𝑥 𝐿  1 − 𝑥 𝐿  2 − ( 3 / 2 ) 𝛼 2 𝑥 𝐿 + 𝑛 2  1 − 𝑥 𝐿  5 𝑥 𝐿  1 − 𝑥 𝐿  5  , 𝑈 0 𝑥 3 𝑥 3  = − 𝐷 2 𝜋 𝜌 1 𝐴 2 + 𝜇  1 − 𝑥 𝐿  3 + 9 𝜇 𝛼 2 2  1 − 𝑥 𝐿  5  , 𝑈 0 𝑥 1 𝑥 2 = 0 = 𝑈 0 𝑥 2 𝑥 3 = 𝑈 0 𝑥 1 𝑥 3 . ( 4 . 6 ) Substituting these in ( 4.2 ), we at once have the variational equations: ̈ 𝜉 − 2 𝑛 ̇ 𝜂 = 𝑈 0 𝑥 1 𝑥 1 ̇ 𝜉 , ̈ 𝜂 + 2 𝑛 𝜉 = 𝑈 0 𝑥 2 𝑥 2 ̈  𝜇 𝜂 , ( 4 . 7 ) 𝜁 = − 𝐷  1 − 𝑥 𝐿  3 + 9 𝜇 𝛼 2 2  1 − 𝑥 𝐿  5 + 2 𝜋 𝜌 1 𝐴 2  𝜁 , ( 4 . 8 ) where the partial derivatives have been computed at each equilibrium point 𝑥 𝐿 . Now, ( 4.8 ) is independent of ( 4.7 ), the solution being a periodic function is bounded and therefore, the motion of the infinitesimal body in the 𝑥 3 direction is stable. Now, the characteristic equation of the equilibrium points ( 𝑥 𝐿 , 0 , 0 ) corresponding to the system ( 4.7 ) is 𝜆 4 −  𝑈 0 𝑥 1 𝑥 1 + 𝑈 0 𝑥 2 𝑥 2 − 4 𝑛 2  𝜆 2 + 𝑈 0 𝑥 1 𝑥 1 𝑈 0 𝑥 2 𝑥 2 = 0 , ( 4 . 9 ) where 𝑈 0 𝑥 1 𝑥 1  𝛼 = 3 𝐷 𝜇 1 2 𝑥 𝐿  𝑈 , ( 4 . 1 0 ) 0 𝑥 2 𝑥 2  𝛼 = 3 𝐷 𝜇 − 1 + 1 2 𝑥 𝐿  . ( 4 . 1 1 ) These equations have been obtained using binomial expansion and ignoring terms with second and higher power in 𝑝 1 , 𝑝 2 , 𝛼 2 , and their product. Now, let 𝜆 2 1 and 𝜆 2 2 be the roots of ( 4.9 ), then, the equilibrium point is stable if both the roots are real and negative. This means that their sum must be negative and their product must be positive. Hence, the points ( 𝑥 𝐿 , 0 , 0 ) will be stable if the following two conditions hold: 𝜆 2 1 + 𝜆 2 2 = 𝑈 0 𝑥 1 𝑥 1 + 𝑈 0 𝑥 2 𝑥 2 − 4 𝑛 2 𝜆 < 0 , ( 4 . 1 2 ) 2 1 𝜆 2 2 = 𝑈 0 𝑥 1 𝑥 1 𝑈 0 𝑥 2 𝑥 2 > 0 . ( 4 . 1 3 ) Now, in the case of the first equilibrium point 𝑥 𝐿 = 𝑝 1 , if we suppose in ( 4.10 ) that 𝑝 1 < 0 then, 𝑈 0 𝑥 1 𝑥 1 < 0 since 0 < 𝜇 < 1 , 𝐷 > 0 , 0 < 𝛼 1 ≪ 1 and when 𝑝 1 > 0 , we have 𝑈 0 𝑥 1 𝑥 1 > 0 . Similarly, in ( 4.11 ), if we suppose 𝑝 1 < 0 then 𝑈 0 𝑥 2 𝑥 2 < 0 . For the case 𝑝 1 > 0 , we will have 𝑈 0 𝑥 2 𝑥 2 > 0 when 𝛼 1 > 2 | 𝑝 1 | which is not possible, hence 𝑈 0 𝑥 2 𝑥 2 < 0 . In the case 0 < 𝑝 1 < 𝛼 1 / 2 , 𝑈 0 𝑥 1 𝑥 1 > 0 , and 𝑈 0 𝑥 2 𝑥 2 > 0 . Also, if 0 < 𝛼 1 / 2 < 𝑝 1 , we see that 𝑈 0 𝑥 1 𝑥 1 > 0 and 𝑈 0 𝑥 2 𝑥 2 < 0 . Thus, for the case 𝑝 1 < 0 , the equilibrium point is stable. For 0 < 𝑝 1 < 𝛼 1 / 2 , the equilibrium point is stable if the condition ( 4.12 ) holds. When 0 < 𝛼 1 / 2 < 𝑝 1 , the equilibrium point is unstable. Next, for the other equilibrium point positioned at 𝑥 𝐿 = 𝑥 1 1 + 𝑝 2 , when 𝑥 1 1 > 0 , then 𝑥  1 1 > 0 since | 𝑝 2 | ≪ 1 and the equilibrium point is stable if the conditions ( 4.12 ) and ( 4.13 ) are satisfied. If 𝑥 1 1 < 0 then 𝑥  1 1 < 0 ; it makes 𝑈 0 𝑥 1 𝑥 1 < 0 , 𝑈 0 𝑥 2 𝑥 2 < 0 . Therefore, when 𝑥 1 1 < 0 , both the conditions ( 4.12 ) and ( 4.13 ) are fulfilled and the equilibrium point is stable. 4.2. Circular Points At circular points ( 1 + 𝑟 c o s 𝜃 , 𝑟 s i n 𝜃 , 0 ) , the values of the second partial derivatives with the use of ( 3.14 ) and neglecting the product 𝛼 1 𝛼 2 are 𝑈 0 𝑥 1 𝑥 1 = 3 𝐷 𝜇 c o s 2 𝜃  𝑛 2 + 𝛼 2  , 𝑈 0 𝑥 1 𝑥 2  𝑛 = 3 𝐷 𝜇 c o s 𝜃 s i n 𝜃 2 + 𝛼 2  , 𝑈 0 𝑥 2 𝑥 2 = 3 𝐷 𝜇 s i n 2 𝜃  𝑛 2 + 𝛼 2  , 𝑈 0 𝑥 3 𝑥 3  = − 𝐷 2 𝜋 𝜌 1 𝐴 2  𝑛 + 𝜇 2 + 3 𝛼 2 .   ( 4 . 1 4 ) Substituting these values in the variational equations ( 4.2 ), we get ̈ 𝜉 − 2 𝑛 ̇ 𝜂 = 3 𝐷 𝜇 c o s 2 𝜃  𝑛 2 + 𝛼 2   𝑛 𝜉 + 3 𝐷 𝜇 c o s 𝜃 s i n 𝜃 2 + 𝛼 2  ̇  𝑛 𝜂 + ( 0 ) 𝜁 , ̈ 𝜂 + 2 𝑛 𝜉 = 3 𝐷 𝜇 c o s 𝜃 s i n 𝜃 2 + 𝛼 2  𝜉 + 3 𝐷 𝜇 s i n 2 𝜃  𝑛 2 + 𝛼 2  ̈  𝜂 + ( 0 ) 𝜁 , ( 4 . 1 5 ) 𝜁 = − 𝐷 2 𝜋 𝜌 1 𝐴 2  𝑛 + 𝜇 2 + 3 𝛼 2   𝜁 . ( 4 . 1 6 ) Equation ( 4.16 ) is independent of ( 4.15 ), it shows that the motion of the infinitesimal mass along the 𝑥 3 -direction is stable. Now, a substitution of these partial derivatives in the characteristic equation ( 4.5 ) yields 𝜆 4 −   𝑛 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2  𝜆 2 = 0 . ( 4 . 1 7 ) Let 𝜆 2 = Λ in ( 4.17 ) then, we have Λ    𝑛 Λ − 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2   = 0 . ( 4 . 1 8 ) Hence, either  𝑛 Λ = 0 o r Λ = 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2 , ( 4 . 1 9 ) which implies that   𝑛 𝜆 = 0 t w i c e , o r 𝜆 = ± 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2  1 / 2 . ( 4 . 2 0 ) Therefore, ( 4.20 ) gives the roots of the characteristic equation ( 4.17 ). Hence we conclude that the circular points are unstable due to the presence of multiple roots. 4.3. Out-of-Plane Points To determine the stability of the out-of-plane equilibrium points, we consider the following partial derivatives: 𝑈 𝑥 1 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 ⎧ ⎪ ⎨ ⎪ ⎩  2  + 𝜇 1 − 𝑥 1  2 − 𝑥 2 3    1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 3 2 𝜇 𝛼 2 ⎧ ⎪ ⎨ ⎪ ⎩  4  1 − 𝑥 1  2 − 𝑥 2 3    1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ − 1 5 2 𝜇 𝛼 2 𝑥 2 3 ⎧ ⎪ ⎨ ⎪ ⎩  6  1 − 𝑥 1  2 − 𝑥 2 3    1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑈 𝑥 1 𝑥 3 ⎡ ⎢ ⎢ ⎢ ⎣  = 𝐷 − 3 𝜇 1 − 𝑥 1  𝑥 3   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 − 1 5 2 𝜇 𝛼 2  1 − 𝑥 1  𝑥 3   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 − 1 5 𝜇 𝛼 2  1 − 𝑥 1  2 2 𝑥 3   1 − 𝑥 1  2 + 𝑥 2 3  − 7 𝑥 3 3   1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑈 𝑥 2 𝑥 2 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 ⎧ ⎪ ⎨ ⎪ ⎩  − 𝜇 1 − 𝑥 1  2 + 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ − 3 𝜇 𝛼 2 2 ⎧ ⎪ ⎨ ⎪ ⎩  1 − 𝑥 1  2 + 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 1 5 𝜇 𝛼 2 𝑥 2 3 2 ⎧ ⎪ ⎨ ⎪ ⎩  1 − 𝑥 1  2 + 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑈 𝑥 3 𝑥 3 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝐴 2 ⎧ ⎪ ⎨ ⎪ ⎩  − 𝜇 1 − 𝑥 1  2 − 2 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ − 9 𝜇 𝛼 2 2 ⎧ ⎪ ⎨ ⎪ ⎩  1 − 𝑥 1  2 − 4 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 1 5 𝜇 𝛼 2 2 ⎧ ⎪ ⎨ ⎪ ⎩ 3 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  − 7 𝑥 4 3   1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ ⎤ ⎥ ⎥ ⎥ ⎦ . ( 4 . 2 1 ) Since 𝑥 2 = 0 , therefore the partial derivatives to be computed at the out-of-plane equilibrium points are 𝑈 0 𝑥 1 𝑥 2 = 0 = 𝑈 0 𝑥 2 𝑥 1 = 𝑈 0 𝑥 2 𝑥 3 = 𝑈 0 𝑥 3 𝑥 2 , 𝑈 0 𝑥 1 𝑥 1  = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 − 2 𝜋 𝜌 1 𝐴 2 + 𝜇 𝑙 3 − 3 𝜇 𝑥 2 3 𝑙 5 + 3 𝜇 𝛼 2 2 𝑙 5 − 9 𝜇 𝛼 2 2 𝑙 5 − 4 5 𝜇 𝛼 2 𝑥 2 3 𝑙 7 + 1 0 5 𝜇 𝛼 2 𝑥 4 3 2 𝑙 9 + 𝑛 2  , 𝑈 0 𝑥 2 𝑥 2  = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 + 2 𝜋 𝜌 1 𝐴 2 + 3 𝜇 𝛼 2 𝑙 5 + 𝑛 2  , 𝑈 0 𝑥 1 𝑥 3  = − 3 1 − 𝑥 1  𝑥 3 𝐷  𝜇 𝑙 5 + 1 5 𝜇 𝛼 2 2 𝑙 7 − 3 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 9  , 𝑈 0 𝑥 3 𝑥 3  = 𝐷 3 𝜇 𝑥 2 3 𝑙 5 + 7 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 − 1 0 5 𝜇 𝛼 2 𝑥 4 3 2 𝑙 9  . ( 4 . 2 2 ) Now, we let, 𝑈 0 𝑥 1 𝑥 1 = 𝑈 1 1 , 𝑈 0 𝑥 2 𝑥 2 = 𝑈 2 2 , 𝑈 0 𝑥 1 𝑥 3 = 𝑈 1 3 , 𝑈 0 𝑥 3 𝑥 3 = 𝑈 3 3 , 𝑈 0 𝑥 1 𝑥 2 = 𝑈 1 2 , 𝑈 0 𝑥 2 𝑥 3 = 𝑈 2 3 . ( 4 . 2 3 ) Using ( 4.23 ), the variational equation can be recast in the form: ̈ 𝜉 − 2 𝑛 ̇ 𝜂 = 𝑈 1 1 𝜉 + 𝑈 1 3 ̇ 𝜁 , ̈ 𝜂 + 2 𝑛 𝜉 = 𝑈 2 2 𝜂 + 𝑈 2 3 ̈ 𝜁 , 𝜁 = 𝑈 1 3 𝜉 + 𝑈 3 3 𝜁 . ( 4 . 2 4 ) In order to consider the motion of the out-of-plane points, we let solution of the system ( 4.24 ) be 𝜉 = 𝐴 e x p ( 𝜆 𝑡 ) , 𝜂 = 𝐵 e x p ( 𝜆 𝑡 ) , 𝜁 = 𝐶 e x p ( 𝜆 𝑡 ) , ( 4 . 2 5 ) where 𝐴 , 𝐵 , 𝐶 , and 𝜆 are constants. 𝜉 , 𝜂 , and 𝜁 are the small displacements in the coordinates of the infinitesimal body. Now, the characteristic equation corresponding to the variational equations ( 4.24 ) in the case of the out-of-plane point may be expressed as 𝜆 6 − 𝑎 1 𝜆 4 + 𝑎 2 𝜆 2 + 𝑎 3 = 0 , ( 4 . 2 6 ) where the coefficients of the characteristic equation ( 4.26 ) are such that 𝑎 1  = 𝐷 − 4 𝜋 𝜌 1 𝐴 1 − 2 𝜋 𝜌 1 𝐴 2 + 2 𝑛 2  , 𝑎 2 = 1 4 𝑙 1 6  𝐷  − 9 𝜇 2 𝐷 𝑥 2 3  − 1 + 𝑥 1  2  2 𝑙 4 + 5 𝛼 2  3 𝑙 2 − 7 𝑥 2 3   + 3 𝑙 4 𝜇 𝑥 2 3  2 𝑙 4 + 5 𝛼 2  5 𝑙 2 − 7 𝜇 𝑥 2 3 ×    3 𝐷 𝜇 𝛼 2 + 2 𝑙 5  ( − 4 + 𝐷 ) 𝑛 2 − 2 𝐷 𝐴 1 𝜋 𝜌 1 + 2 𝐷 𝐴 2 𝜋 𝜌 1    + 𝐷 6 𝑙 4 𝜇 𝑥 2 3 + 3 𝛼 2  𝑙 4 𝜇 − 2 5 𝑙 2 − 3 5 𝜇 𝑥 4 3  + 2 𝑙 9  𝑛 2 − 2 𝐴 1 𝜋 𝜌 1 + 2 𝐴 2 𝜋 𝜌 1 ×    3 𝛼 2  𝑙 4 𝜇 − 2 𝑙 2  1 5 + 𝑙 2  𝑥 2 3 + 3 5 𝜇 𝑥 4 3  + 2 𝑙 6  𝑙 3 𝑛 2 + 𝜇 − 2 𝑙 3 𝐴 1 𝜋 𝜌 1 − 2 𝑙 3 𝐴 2 𝜋 𝜌 1 𝑎     3 = 3 𝐷 3 4 𝑙 1 8 𝜇 𝑥 2 3  2 𝑙 4 + 5 𝛼 2  5 𝑙 2 − 7 𝜇 𝑥 2 3  𝑛   2 + 3 𝜇 𝛼 2 2 𝑙 5 − 2 𝐴 1 𝜋 𝜌 1 + 2 𝐴 2 𝜋 𝜌 1  ×  3 𝛼 2  𝑙 4 𝜇 − 2 𝑙 2  1 5 + 𝑙 2  𝜇 𝑥 2 3 + 3 5 𝜇 𝑥 4 3  + 2 𝑙 6  𝑙 3 𝑛 2 + 𝜇 − 2 𝑙 3 𝐴 1 𝜋 𝜌 1 − 2 𝑙 3 𝐴 2 𝜋 𝜌 1 ,   ( 4 . 2 7 ) where 𝑙 2 = ( 1 − 𝑥 1 ) 2 + 𝑥 2 3 . These computations have been done using the software package Mathematica. For the stability analysis of the out-of-plane equilibrium point, we compute numerically the partial derivatives calculated at the out-of-plane points with the use of ( 3.28 ) and the following numerical values: 𝜇 = 0 . 0 1 , 𝜋 = 3 . 1 4 , 𝛼 1 = 0 . 0 2 4 , 𝛼 2 𝐴 = 0 . 0 2 1 = 0 . 7 , 𝐴 2 = 0 . 6 8 , 𝜌 1 = 0 . 2 3 6 , 𝐷 = 0 . 2 1 3 3 . ( 4 . 2 8 ) Now, substituting the above values in the characteristic equation ( 4.27 ), we get 𝜆 6 − 0 . 1 9 2 4 3 7 𝜆 4 + 0 . 2 4 8 0 3 4 𝜆 2 − 0 . 0 0 0 0 1 9 7 6 2 4 = 0 . ( 4 . 2 9 ) Its roots are: 𝜆 1 , 2 𝜆 = − 0 . 5 4 5 0 6 6 ± 0 . 4 4 8 2 3 9 𝑖 , 3 , 4 𝜆 = ± 0 . 0 0 8 9 2 6 4 2 , 5 , 6 = 0 . 5 4 5 0 6 6 ± 0 . 4 4 8 2 3 9 𝑖 . ( 4 . 3 0 ) The positive root and the positive real part of the complex roots induce instability at the out-of-plane point. Hence, the motion of the infinitesimal mass around the out-of-plane equilibrium points is unstable for the specific numerical example given here. However, fuller discussion of their stability remains a theme for future research. 5. Discussion The equation of motion ( 2.3 ) is different from those of Hallan and Mangang [ 4 ] due to oblateness of the second primary. If we assume that the second primary is not oblate (i.e., 𝛼 2 = 0 ), then these equations will fully coincide with those of Hallan and Mangang [ 4 ]. Equation ( 3.9 ) gives the equilibrium position of the point ( 𝑃 1 , 0 , 0 ) near the center of the first primary and fully coincides with that of Hallan and Mangang [ 4 ]. It shows that the position of this equilibrium point does not depend on oblateness of the second primary, while the other equilibrium point ( 𝑥 1 1 + 𝑃 2 , 0 , 0 ) given by ( 3.11 ) is different from that of Hallan and Mangang [ 4 ] due to the appearance of oblateness of the second primary. When 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) , points on the circle ( 1 − 𝑥 1 ) 2 + 𝑥 2 2 = 𝑟 2 , 𝑥 3 = 0 lying within the first primary are also equilibrium points. These points are affected by oblateness of both primaries. Equations ( 3.28 ) give the positions of the out-of-plane points when only linear terms in oblateness of the second primary are retained. We have been able to show that the oblateness of the primaries allows the existence of the out-of-plane equilibrium points in the 𝑥 1 𝑥 3 -plane within the first primary. These points have no analogy in the previous studies of the Robe’s restricted three-body problem. The linear stability analysis of the equilibrium solutions of the problem is investigated with the help of characteristic roots. The characteristic equation ( 4.9 ) in the case of the equilibrium point 𝑥 𝐿 = 𝑝 1 near the center of the first primary is the same as that of Hallan and Mangang [ 4 ], while that of the other point 𝑥 𝐿 = 𝑥 1 1 + 𝑝 2 near the center differs from that of Hallan and Mangang [ 4 ] due to oblateness of the second primary. The characteristic equation of the circular case ( 4.17 ) also differs from that of Hallan and Mangang [ 4 ] due to oblateness of the second primary. The stability in the first approximation of this configuration shows that points near the centre of the first primary are conditionally stable; the circular points are unstable. This confirms the earlier results of Hallan and Rana [ 3 ], Hallan and Mangang [ 4 ]. A numerical exploration shows that the out-of-plane equilibrium points are also unstable. This outcome validates the earlier results of Douskos and Markellos [ 8 ] and Singh and Leke [ 11 ] that the points are unstable. 6. Conclusion We have derived the equations of motion and established the positions of the equilibrium points of the infinitesimal body in the Robe’s [ 1 ] restricted three-body problem with oblateness. The term “oblateness” is used in the sense that both primaries are considered as oblate spheroids under the effects of the full buoyancy force exerted by the fluid on the infinitesimal mass. We have obtained one equilibrium point ( 𝑃 1 , 0 , 0 ) near the centre of the first primary which will be on the left or right of the centre of the first primary accordingly as 2 𝜋 𝜌 1 𝐴 1 − 2 𝜇 > < 1 . This point is the same as that of Hallan and Mangang [ 4 ]. In addition to this, another equilibrium point ( 𝑥 1 1 + 𝑃 2 , 0 , 0 ) is found within the first primary on the line joining the center of the primaries when 1 − 2 𝜋 𝜌 1 𝐴 1 < − 3 𝜇 / 4 and | 𝑥 1 1 | < 𝑎 1 . When 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) , points on the circle ( 1 − 𝑥 1 ) 2 + 𝑥 2 2 = 𝑟 2 , 𝑥 3 = 0 lying within the first primary are also equilibrium points. We call them circular points. Finally, we have been able to show that the oblateness of the primaries allows the existence of the out-of-plane equilibrium points in the 𝑥 1 𝑥 3 -plane within the first primary. The result of this paper can be summarized as follows. The restricted three-body problem under the framework of the Robe’s [ 1 ] problem with oblate primaries has the equilibrium points of the type: points near the center of the first primary, points on the circle (circular points), and two out-of-plane points 𝐿 6 , 7 . It is seen that points near the first primary are conditionally stable, the circular points are unstable, while the out-of-plane equilibrium points are unstable for the specific numerical example given here. The effect of drag forces as considered by Giordano et al. [ 2 ] under the present context, particularly as regards the analysis of the properties of the equilibrium points located inside the first primary, will be interesting. References H. A. G. Robe, “A new kind of three body problem,” Celestial Mechanics and Dynamical Astronomy , vol. 16, pp. 343–351, 1977. C. M. Giordano, A. R. Plastino, and A. Plastino, “Robe's restricted three-body problem with drag,” Celestial Mechanics and Dynamical Astronomy , vol. 66, no. 2, pp. 229–242, 1996. View at Publisher · View at Google Scholar P. P. Hallan and N. Rana, “The existence and stability of equilibrium points in the Robe's restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy , vol. 79, no. 2, pp. 145–155, 2001. View at Publisher · View at Google Scholar P. P. Hallan and K. B. Mangang, “Existence and linear stability of equilibrium points in the Robe's restricted three body problem when the first primary is an oblate spheroid,” Planetary and Space Science , vol. 55, no. 4, pp. 512–516, 2007. View at Publisher · View at Google Scholar · View at Scopus A. R. Plastino and A. Plastino, “Robe's restricted three-body problem revisited,” Celestial Mechanics and Dynamical Astronomy , vol. 61, no. 2, pp. 197–206, 1995. View at Publisher · View at Google Scholar · View at Scopus R. K. Sharma and P. V. Subba Rao, “Stationary solutions and their characteristic exponents in the restricted three-body problem when the more massive primary is an oblate spheroid,” Celestial Mechanics , vol. 13, no. 2, pp. 137–149, 1976. J. Singh and B. Ishwar, “Stability of triangular points in the generalized photogravitational restricted three-body problem,” Bulletin of the Astronomical Society of India , vol. 27, pp. 415–424, 1999. C. N. Douskos and V. V. Markellos, “Out-of-plane equilibrium points in the restricted three-body problem with oblateness,” Astronomy and Astrophysics , vol. 446, no. 1, pp. 357–360, 2006. View at Publisher · View at Google Scholar · View at Scopus J. Singh and O. Leke, “Stability of the photogravitational restricted three-body problem with variable masses,” Astrophysics and Space Science , vol. 326, pp. 305–314, 2010. V. V. Radzievskii, “The restricted problem of three bodies taking account of light pressure,” Astronomicheskii Zhurnal , vol. 27, pp. 250–256, 1950. J. Singh, “Combined effects of perturbations, radiation and oblateness on the nonlinear stability of triangular points in the restricted three-body problem,” Astrophysics and SpaceScience , vol. 332, pp. 331–339, 2011. // http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Advances in Mathematical Physics Hindawi Publishing Corporation

Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem with Oblateness

Loading next page...
 
/lp/hindawi-publishing-corporation/existence-and-linear-stability-of-equilibrium-points-in-the-robe-s-WsYZnnk36m

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

Publisher
Hindawi Publishing Corporation
Copyright
Copyright © 2012 Jagadish Singh and Abubakar Umar Sandah.
ISSN
1687-9120
eISSN
1687-9139
Publisher site
See Article on Publisher Site

Abstract

Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem with Oblateness Hindawi Publishing Corporation Home Journals About Us About this Journal Submit a Manuscript Table of Contents Journal Menu Abstracting and Indexing Aims and Scope Article Processing Charges Articles in Press Author Guidelines Bibliographic Information Contact Information Editorial Board Editorial Workflow Free eTOC Alerts Reviewers Acknowledgment Subscription Information Open Special Issues Published Special Issues Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Full-Text ePUB Linked References How to Cite this Article Complete Special Issue Advances in Mathematical Physics Volume 2012 (2012), Article ID 679063, 18 pages doi:10.1155/2012/679063 Research Article Existence and Linear Stability of Equilibrium Points in the Robe’s Restricted Three-Body Problem with Oblateness Jagadish Singh 1 and Abubakar Umar Sandah 2 1 Department of Mathematics, Faculty of Science, Ahmadu Bello University, Zaria, Nigeria 2 Department of Mathematics Statistics and Computer Science, College of Science and Technology, Kaduna Polytechnic, Kaduna, Nigeria Received 26 March 2012; Accepted 2 July 2012 Academic Editor: Burak Polat Copyright © 2012 Jagadish Singh and Abubakar Umar Sandah. This is an open access article distributed under the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract This paper investigates the positions and linear stability of an infinitesimal body around the equilibrium points in the framework of the Robe’s circular restricted three-body problem, with assumptions that the hydrostatic equilibrium figure of the first primary is an oblate spheroid and the second primary is an oblate body as well. It is found that equilibrium point exists near the centre of the first primary. Further, there can be one more equilibrium point on the line joining the centers of both primaries. Points on the circle within the first primary are also equilibrium points under certain conditions and the existence of two out-of-plane points is also observed. The linear stability of this configuration is examined and it is found that points near the center of the first primary are conditionally stable, while the circular and out of plane equilibrium points are unstable. 1. Introduction Robe [ 1 ] considered a new kind of restricted three-body problem in which, one of the primaries of mass 𝑚 1 is a rigid spherical shell, filled with homogenous, incompressible fluid of density 𝜌 1 ; the second one is a point mass 𝑚 2 located outside the shell and moving around the mass 𝑚 1 in a Keplerian orbit; the infinitesimal mass 𝑚 3 is a small sphere of density 𝜌 3 , moving inside the shell and is subject to the attraction of 𝑚 2 and the buoyancy force due to the fluid of the first primary. Further, he discussed the linear stability of an equilibrium point obtained in two cases. In the first case, the orbit of 𝑚 2 around 𝑚 1 is circular and in the second case, the orbit is elliptic, but the shell is empty (there is no fluid inside it) or densities of 𝑚 1 and 𝑚 3 are equal. Since then various studies (e.g., [ 2 – 4 ]) under different assumptions have been carried out. In his study, Robe [ 1 ] assumed that the pressure field of the fluid 𝜌 1 has a spherical symmetry around the center of the shell and he took into account only one out of the three components of the pressure field which is due to the own gravitational field of the fluid 𝜌 1 . He did not consider the other two components arising from the attraction of 𝑚 2 and the centrifugal force. Taking care of all these three components of the pressure field, A. R. Plastino and A. Plastino [ 5 ] reanalyzed the Robe’s. But in their study, they assumed the hydrostatic equilibrium figure of the first primary as Roche’s ellipsoid (see Figure 1 ). They found that when the density parameter 𝐷 is taken as zero, every point inside the fluid is an equilibrium point; otherwise the center of the ellipsoid is the only equilibrium point and it is linearly stable. Figure 1: The Robe’s CRTBP with oblate primaries. Hallan and Rana [ 3 ] investigated the existence of all equilibrium point and their stability in the Robe’s [ 1 ] restricted three-body problem. It was seen that the Robe’s elliptic restricted three-body problem has only one equilibrium point for all values of the density parameter 𝐾 and the mass parameter 𝜇 , while the Robe’s circular restricted three-body problem can have two, three, or infinite numbers of equilibrium points. As regards to the stability of these equilibria, they confirmed the stability result given by Robe [ 1 ] of the equilibrium point ( − 𝜇 , 0 , 0 ) , whereas triangular and circular points are always unstable. The equilibrium point collinear with the center of the shell and the second primary was found to be stable under some conditions. Hallan and Mangang [ 4 ] studied the Robe’s [ 1 ] restricted three-body problem by considering the full buoyancy force as in A. R. Plastino and A. Plastino [ 5 ] and assuming the hydrostatic equilibrium figure of the first primary as an oblate spheroid. They derived the pertinent equations of motion and discussed the existence of equilibrium point and their linear stability. The participating bodies in the classical restricted three-body problem are strictly spherical in shape, but in actual situations several heavenly bodies, such as Saturn and Jupiter, are sufficiently oblate. The minor planets and meteoroids have irregular shape. The lack of sphericity, or the oblateness, of the planet causes large perturbations from a two-body orbit. The motions of artificial Earth satellites are examples of this. Global studies of problems with oblateness have been carried out by many researchers (e.g., [ 6 – 9 ]). Therefore, our effort in this paper aims at investigating the equilibrium points and their stability in the Robe’s circular restricted three-body problem when the hydrostatic equilibrium figure of the fluid of the first primary is an oblate spheroid and the second one is an oblate spheroid as well. The model of this study can be used to study the small oscillation of the Earth’s inner core taking into account the Moon’s attraction. This paper is organized as follows; Section 2 represents the equations of motion; the existence of the equilibrium points is mentioned in Section 3 , while Section 4 investigates their linear stability; Section 5 discusses the results obtained; the conclusion is drawn in Section 6 . 2. Equation of Motion Let the first primary 𝑚 1 be a fluid of density 𝜌 1 in the shape of an oblate spheroid as assumed by Hallan and Mangang [ 4 ]; let the second primary 𝑚 2 be an oblate body too as Sharma and Subba Rao [ 6 ] assumed, which describes a circular orbit around 𝑚 1 . We adopt a uniformly rotating coordinate system 𝑂 𝑥 1 𝑥 2 𝑥 3 with origin at the center of mass 𝑚 1 , 𝑂 𝑥 1 pointing towards 𝑚 2 , with 𝑂 𝑥 1 𝑥 2 being the orbital plane of 𝑚 2 coinciding with the equatorial plane of 𝑚 1 . Then, the equations of motion of the infinitesimal body of density 𝜌 3 in the coordinate system take the form [ 4 , 6 ]: ̈ 𝑥 1 − 2 𝑛 ̇ 𝑥 2 = 𝜕 𝑈 𝜕 𝑥 1 , ̈ 𝑥 2 + 2 𝑛 ̇ 𝑥 1 = 𝜕 𝑈 𝜕 𝑥 2 , ̈ 𝑥 3 = 𝜕 𝑈 𝜕 𝑥 3 , ( 2 . 1 ) where 𝑛 𝑈 = 𝑉 + 2   𝑥 1 −  𝑚 2 /  𝑚 1 + 𝑚 2 𝑅    2 + 𝑥 2 2  2 , 𝑉 = 𝐵 + 𝐵  − 𝜌 1 𝜌 3 ⎡ ⎢ ⎢ ⎣ 𝐵 + 𝐵  + 𝑛 2   𝑥 1 −  𝑚 2 /  𝑚 1 + 𝑚 2 𝑅    2 + 𝑥 2 2  2 ⎤ ⎥ ⎥ ⎦ , 𝐵 = 𝜋 𝐺 𝜌 1  𝐼 − 𝐴 1 𝑥 2 1 − 𝐴 1 𝑥 2 2 − 𝐴 2 𝑥 2 3  , 𝐵  = 𝐺 𝑚 2   𝑅 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  1 / 2 + 𝐺 𝑚 2 𝛼 2 2   𝑅 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝐺 𝑚 2 𝛼 2 𝑥 2 3 2   𝑅 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 , 𝐼 = 2 𝑎 2 1 𝐴 1 + 𝑎 2 2 𝐴 2 , 𝐴 1 = 𝑎 2 1 𝑎 2  ∞ 0 𝑑 𝑢 Δ  𝑎 2 1  + 𝑢 , 𝐴 2 = 𝑎 2 1 𝑎 2  ∞ 0 𝑑 𝑢 Δ  𝑎 2 2  , Δ + 𝑢 2 =  𝑎 2 1  + 𝑢 2  𝑎 2 2  , 𝑛 + 𝑢 2 = 𝐺  𝑚 1 + 𝑚 2  𝑅 2  3 1 + 2 𝛼 1 + 3 2 𝛼 2  ; 𝛼 1 = 𝑎 2 1 − 𝑎 2 3 5 𝑅 2 , 𝛼 2 = 𝑎 2 2 − 𝑎 2 4 5 𝑅 2 , 𝜌 𝐷 = 1 − 1 𝜌 3 . ( 2 . 2 ) Here 𝑉 is the potential that explains the combined forces upon the infinitesimal mass, 𝐵 denotes the potential due to the fluid mass of the first primary, 𝐵 ′ stands for the potential due to the second primary, 𝑅 is the distance between the primaries, and 𝐺 is the gravitational constant. 𝑛 is the mean motion. 𝑎 1 , 𝑎 2 and 𝑎 3 , 𝑎 4 are the equatorial and polar radii of the first and second primary, respectively. 𝐼 stands for the polar moment of inertia, while 𝐴 𝑖 ( 𝑖 = 1 , 2 ) are the index symbols. 𝛼 1 and 𝛼 2 are the oblateness coefficients of the first and second primaries, respectively. We choose the unit of mass such that the sum of the masses of the primaries is taken as unity, thus we take 𝑚 2 = 𝜇 , 0 < 𝜇 = 𝑚 2 / ( 𝑚 1 + 𝑚 2 ) < 1 . For the unit of length, we take the distance between the primaries as unity, that is, 𝑅 = 1 and the unit of time is also selected such that 𝐺 = 1 . With these units and substituting the expression for the potential 𝐵 due to the fluid in the first primary and the potential 𝐵 ′ due to the second oblate primary, the equations of motion ( 2.1 ) are recast to the form: ̈ 𝑥 1 − 2 𝑛 ̇ 𝑥 2 = 𝜕 𝑈 𝜕 𝑥 1 , ̈ 𝑥 2 + 2 𝑛 ̇ 𝑥 1 = 𝜕 𝑈 𝜕 𝑥 2 , ̈ 𝑥 3 = 𝜕 𝑈 𝜕 𝑥 3 , ( 2 . 3 ) where ⎡ ⎢ ⎢ ⎢ ⎣ 𝑈 = 𝐷 𝜋 𝜌 1  𝐼 − 𝐴 1  𝑥 2 1 + 𝑥 2 2  − 𝐴 2 𝑥 2 3  + 𝜇   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  1 / 2 + 𝜇 𝛼 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 + 𝑛 2   𝑥 1  − 𝜇 2 + 𝑥 2 2  2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑛 2 =  3 1 + 2 𝛼 1 + 3 2 𝛼 2  . ( 2 . 4 ) These above equations of motion of the infinitesimal mass 𝑚 3 under the framework of the Robe’s circular restricted three-body problem have been obtained by taking into account the shapes of the primaries, the full buoyancy force, the forces due to the gravitational attraction of the second primary, and the gravitational force exerted by the fluid of density 𝜌 1 . In the case when the second primary is not an oblate spheroid (i.e., 𝛼 2 = 0 ), the equations are the same as those of Hallan and Mangang [ 4 ]. 3. Position of Equilibrium Points The equilibrium points are the solutions of the equations: 𝑈 𝑥 1 = 𝑈 𝑥 2 = 𝑈 𝑥 3 = 0 . ( 3 . 1 ) That is, 𝑈 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝑥 1 𝐴 1 + 𝜇  1 − 𝑥 1    1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 + 3 𝜇 𝛼 2  1 − 𝑥 1  2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 − 1 5 𝜇 𝛼 2  1 − 𝑥 1  𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  7 / 2 + 𝑛 2  𝑥 1  ⎤ ⎥ ⎥ ⎥ ⎦ 𝑈 − 𝜇 = 0 , ( 3 . 2 ) 𝑥 2 = 𝐷 𝑥 2 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  7 / 2 + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ 𝑈 = 0 , ( 3 . 3 ) 𝑥 3 = 𝐷 𝑥 3 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 2 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 − 3 𝜇 𝛼 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 2 + 𝑥 2 3  7 / 2 ⎤ ⎥ ⎥ ⎥ ⎦ = 0 . ( 3 . 4 ) 3.1. Equilibrium Points Near the Centre of the First Primary The positions of the equilibrium points near the first primary are the solutions of ( 3.2 ) when 𝑈 𝑥 1 = 0 , 𝑥 1 ≠ 0 , 𝑥 2 = 𝑥 3 = 0 , 𝐷 ≠ 0 , and 𝑛 2 = 1 + ( 3 / 2 ) ( 𝛼 1 + 𝛼 2 ) . The 𝑥 1 coordinate of the equilibrium points are then the roots of the equation: − 2 𝜋 𝜌 1 𝑥 1 𝐴 1 + 𝜇  1 − 𝑥 1  | | 1 − 𝑥 1 | | 3 + 3 𝜇 𝛼 2  1 − 𝑥 1  2 | | 1 − 𝑥 1 | | 5 +  3 1 + 2 𝛼 1 + 3 2 𝛼 2   𝑥 1  − 𝜇 = 0 . ( 3 . 5 ) We first determine the roots of ( 3.5 ) in the absence of oblateness, that is, the case when the primaries are spherical. In this case, the roots are [ 4 ] 𝑥 1 1 √ = 1 + 𝜇 + 𝜇 2 + 8 𝜇 𝜋 𝜌 1 𝐴 1 − 4 𝜇 2  1 − 2 𝜋 𝜌 1 𝐴 1  , 𝑥 1 2 √ = 1 + 𝜇 − 𝜇 2 + 8 𝜇 𝜋 𝜌 1 𝐴 1 − 4 𝜇 2  1 − 2 𝜋 𝜌 1 𝐴 1  . ( 3 . 6 ) The term 𝐴 1 which appears in ( 3.2 ) and is due to the fluid mass affects these roots. Therefore, these roots will be real if the discriminant is nonnegative, that is if 𝜇 + 8 𝜋 𝜌 1 𝐴 1 − 4 ≥ 0 . ( 3 . 7 ) When ( 1 / 4 ) 𝜇 ≥ 1 − 2 𝜋 𝜌 1 𝐴 1 > 0 , both roots are greater than unity and we reject them because they lie outside the first primary. Now, if 1 − 2 𝜋 𝜌 1 𝐴 1 < 0 , we have 𝑥 1 2 > 1 and 𝑥 1 1 < 1 . Further, we see that 𝑥 1 1 > − 1 when 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 . Thus, in the case when 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 , the point ( 𝑥 1 1 , 0 , 0 ) lies within the first primary if | 𝑥 1 1 | < 𝑎 1 . When 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 , | 𝑥 1 1 | < 𝑎 1 ; 𝑥 1 1 is a root of ( 3.5 ). Hence for 1 − 2 𝜋 𝜌 1 𝐴 1 = 0 , the only root is 𝑥 1 1 = 2 which lies outside the first primary and we neglect it. Hence, for 𝛼 1 = 0 , 𝛼 2 = 0 , 𝑥 1 1 = 0 is always a root of ( 3.5 ) and 𝑥 1 = 𝑥 1 1 is also a root provided 1 − 2 𝜋 𝜌 1 𝐴 1 < − ( 3 / 4 ) 𝜇 , | 𝑥 1 1 | < 𝑎 1 . Now, we find the roots of ( 3.5 ) when oblateness of both primaries is considered (i.e., 𝛼 1 ≠ 0 , 𝛼 2 ≠ 0 ). Let the roots be such that 𝑥 1 = 0 + 𝑝 1 , | | 𝑝 1 | | 𝑥 ≪ 1 , 1 = 𝑥 1 1 + 𝑝 2 , | | 𝑝 2 | | ≪ 1 . ( 3 . 8 ) Putting these values in ( 3.5 ), multiplying throughout by ( 1 − 𝑝 1 ) 4 , expanding and neglecting second and higher powers of 𝑝 1 , 𝛼 1 , 𝛼 2 , as they are very small quantities, we have 𝑝 1 ≅ − 3 𝛼 1 2 𝜇 2 𝜋 𝜌 1 𝐴 1 . − ( 1 + 2 𝜇 ) ( 3 . 9 ) Similarly, putting 𝑥 1 = 𝑥 1 1 + 𝑝 2 in ( 3.5 ) and then simplifying it, we get  1 − 𝑥 1 1  4   𝑥 1 1   3 − 𝜇 1 + 2 𝛼 1 + 3 2 𝛼 2  − 2 𝜋 𝜌 1 𝐴 1 𝑥 1 1  + 𝑝 2  1 − 𝑥 1 1  3   1 − 3 𝑥 1 1   1 − 2 𝜋 𝜌 1 𝐴 1    + 2 𝜇 + 𝜇 1 − 𝑥 1 1  2 − 2 𝑝 2  1 − 𝑥 1 1    1 − 𝑥 1 1  2  𝑥 1 1 − 2 𝜋 𝜌 1 𝐴 1 𝑥 1 1   − 𝜇 + 𝜇 = − 3 𝜇 𝛼 2 2 . ( 3 . 1 0 ) Multiplying ( 3.8 ) by ( 1 − 𝑥 1 1 ) 2 , simplifying and then using it in ( 3.10 ), we get 𝑝 2  ≅ − 1 − 𝑥 1 1 𝑥    1 1 − 𝜇   ( 3 / 2 ) 𝛼 1 + ( 3 / 2 ) 𝛼 2     1 − 3 𝑥 1 1   1 − 2 𝜋 𝜌 1 𝐴 1   − + 2 𝜇 3 𝜇 𝛼 2 2  1 − 𝑥 1 1  3   1 − 3 𝑥 1 1   1 − 2 𝜋 𝜌 1 𝐴 1   . + 2 𝜇 ( 3 . 1 1 ) A substitution of ( 3.11 ) in the second equation of ( 3.8 ) at once gives the position of the other equilibrium point near the center of the first primary. 3.1.1. Positions of Circular Points The positions of the circular points are sought using the first two equations of system ( 3.1 ) with the conditions 𝑥 1 ≠ 0 , 𝑥 2 ≠ 0 , 𝑥 3 = 0 ; that is, they are the solutions of 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2 3 + 1 + 2 𝛼 1 + 3 2 𝛼 2 ⎤ ⎥ ⎥ ⎥ ⎦ + 𝜇   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 + 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2  3 − 𝜇 1 + 2 𝛼 1 + 3 2 𝛼 2  = 0 , ( 3 . 1 2 ) − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2 3 + 1 + 2 𝛼 1 + 3 2 𝛼 2 = 0 . ( 3 . 1 3 ) Solving the above equations and knowing that 𝜇 ≠ 0 , we get 1   1 − 𝑥 1  2 + 𝑥 2 2  3 / 2 + 3 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 2  5 / 2 − 𝑛 2 = 0 . ( 3 . 1 4 ) We let  1 − 𝑥 1  2 + 𝑥 2 2 = 𝑟 2 . ( 3 . 1 5 ) Substituting ( 3.15 ) in ( 3.14 ), and simplifying, we get 𝑛 2 𝑟 5 − 𝑟 2 − 3 2 𝛼 2 = 0 . ( 3 . 1 6 ) Now, we let 𝑟 = 1 + 𝜀 , 𝜀 ≪ 1 . ( 3 . 1 7 ) Substituting ( 3.16 ) in ( 3.15 ), neglecting second and higher powers of 𝜀 , we get 1 𝜀 = − 2 𝛼 1 . ( 3 . 1 8 ) Therefore, ( 3.17 ) is now expressed as 1 𝑟 ≅ 1 − 2 𝛼 1 . ( 3 . 1 9 ) A substitution of ( 3.14 ) in ( 3.13 ) yields 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) . ( 3 . 2 0 ) Therefore, when 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) , the points on the circle given by ( 3.15 ) with 𝑥 3 = 0 and 𝑟 = 1 − ( 1 / 2 ) 𝛼 1 lying within the first primary are also equilibrium points. The general coordinates of these circular points are given by ( 1 + 𝑟 c o s 𝜃 , 𝑟 s i n 𝜃 , 0 ) , where 𝜃 is a parameter. When 𝑦 = 0 , the circular points coalesce to those lying on the line joining the primaries. 3.1.2. Positions of Out-of-Plane Equilibrium Points The out-of-plane points have no analogy in the classical restricted three-body problem. However the investigation concerning these points in the photogravitational restricted three-body problem was first carried out by Radzievskii [ 10 ]. Afterwards, other researchers, for instance Douskos and Markellos [ 8 ], Singh and Leke [ 11 ], and so forth, have worked on the out-of-plane points. In this section, we locate these points for our study, as it has remained an open problem to date. The positions of the out-of-plane equilibrium points of the Robe’s problem with oblate primaries are the solutions of the first and last equations of ( 3.1 ) with 𝑥 2 = 0 , 𝐷 ≠ 0 ; that is, 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 1 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 − 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ + 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 + 3 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 − 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 − 𝑛 2 𝑥 𝜇 = 0 , ( 3 . 2 1 ) 3 ⎡ ⎢ ⎢ ⎢ ⎣ − 2 𝜋 𝜌 1 𝐴 2 − 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 − 9 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 + 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎤ ⎥ ⎥ ⎥ ⎦ = 0 . ( 3 . 2 2 ) From ( 3.22 ), since 𝑥 3 ≠ 0 , we have − 2 𝜋 𝜌 1 𝐴 2 = 𝜇   1 − 𝑥 1  2 + 𝑥 2 3  3 / 2 + 9 𝜇 𝛼 2 2   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 − 1 5 𝜇 𝛼 2 𝑥 2 3 2   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 . ( 3 . 2 3 ) Let 𝑙 2 =  1 − 𝑥 1  2 + 𝑥 2 3 . ( 3 . 2 4 ) Then, ( 3.23 ) and ( 3.21 ) may be written respectively: 1 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 = 2 𝜋 𝜌 1 𝐴 2 + 𝜇 𝑙 3 + 9 𝜇 𝛼 2 2 𝑙 5 , 𝑥 1  − 2 𝜋 𝜌 1 𝐴 1 − 𝜇 𝑙 3 − 3 𝜇 𝛼 2 2 𝑙 5 + 1 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 + 𝑛 2  + 𝜇 𝑙 3 + 3 𝜇 𝛼 2 2 𝑙 5 − 1 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 − 𝑛 2 𝜇 = 0 . ( 3 . 2 5 ) Now, from first equations ( 3.25 ), we get 𝑥 3 𝑙 = ± √ 1 5 𝜇 𝛼 2  𝜇  9 𝛼 2 + 2 𝑙 2  + 4 𝜋 𝜌 1 𝐴 2 𝑙 5  1 / 2 . ( 3 . 2 6 ) The use of ( 3.24 ) in second equation of ( 3.25 ) yields 𝑥 1 =  2 𝜋 𝜌 1 𝐴 2 + 𝑛 2 𝜇  𝑙 5 + 3 𝜇 𝛼 2  2 𝜋 𝜌 1  𝐴 2 − 𝐴 1  + 𝑛 2  𝑙 5 + 3 𝜇 𝛼 2 . ( 3 . 2 7 ) We use the software package Mathematica (Wolfram 2004) to compute the coordinates of the out-of-plane equilibrium points denoted by 𝐿 6 and 𝐿 7 starting with the initial values 𝑥 1 = 1 − 𝜇 and 𝑥 3 = √ 3 √ 𝛼 2 in the case where we have kept up to first order terms in both the numerator and the denominator; we then get 𝑥 1 = 2 𝜇 + 3 𝜇 𝛼 1 + 4 𝐴 2 𝜋 𝜌 1 2 + 3 𝛼 1 − 4 𝐴 1 𝜋 𝜌 1 + 4 𝐴 2 𝜋 𝜌 1 − 3  𝜇 𝐴 1 + ( 1 − 𝜇 ) 𝐴 2  𝛼 2 𝜋 𝜌 1  1 + 3 𝛼 1  1 + 2 𝐴 2 𝜋 𝜌 1 − 2 𝐴 1 𝜋 𝜌 1  + 4 𝜋 𝜌 1  𝐴 2 − 𝐴 1 + 𝐴 2 1 𝜋 𝜌 1 + 𝐴 2 2 𝜋 𝜌 1 − 2 𝐴 1 𝐴 2 𝜋 𝜌 1 , 𝑥   3 = √  2 / 1 5 𝜇 𝜇 3  1 + 2 𝜋 𝐴 2 𝜇 2 𝜌 1  √ 𝜇 𝛼 2 + 7 √  3 𝜇 / 1 0 𝜇 + 2 𝜋 𝐴 2 𝜇 3 𝜌 1  √ 𝛼 2 2  𝜇 3  1 + 2 𝜋 𝐴 2 𝜇 2 𝜌 1  . ( 3 . 2 8 ) The location of the out-of-plane equilibrium points can be obtained by solving numerically equations ( 3.26 ) and ( 3.27 ) using ( 3.24 ). Now, from the expression for the density parameter  𝜌 𝐷 = 1 − 1 𝜌 3  . ( 3 . 2 9 ) We assume that 𝜌 1 ≠ 𝜌 3 , then 𝐷 > 𝑜 𝑟 < 0 . In the case when the density parameter is positive, we have 𝜌 1 < 𝜌 3 . ( 3 . 3 0 ) Hence, numerically we choose 𝑎 2 1 = 0 . 9 4 , 𝑎 2 2 = 0 . 9 , 𝑎 2 3 = 0 . 8 2 , 𝑎 2 4 = 0 . 8 , 𝜇 = 0 . 0 1 , 𝜋 = 3 . 1 4 . ( 3 . 3 1 ) Then, 𝛼 1 = 0 . 0 2 4 , 𝛼 2 = 0 . 0 2 , 𝜌 1 = 0 . 2 3 6 . ( 3 . 3 2 ) Now, we perform a numerical exploration of computing the out-of-plane points in the case of the Earth-Moon system. To do this, we arbitrarily choose values for the 𝐴 𝑖 ( 𝑖 = 1 , 2 ) . We found that when 𝐴 1 = 2 . 5 0 7 6 and 𝐴 2 = 2 . 5 5 5 , the positions of the out-of-plane points ( 𝑥 1 , 0 , ± 𝑥 3 ) : 𝑥 1 = 3 . 3 5 2 7 , 𝑥 3 = 0 . 2 7 1 4 1 8 . ( 3 . 3 3 ) The abscissae of the out-of-plane point is outside the possible region of motion of the infinitesimal mass and so we neglect it. However, in the case when the 𝐴 𝑖 ( 𝑖 = 1 , 2 ) are chosen such that | | 𝐴 1 − 𝐴 2 | | ≪ 1 , 𝐴 1 = 0 . 7 , 𝐴 2 = 0 . 6 8 ( s a y ) , ( 3 . 3 4 ) 𝐴 𝑖 ∈ ( 0 , 0 . 7 ] , the point 𝐿 6 , and 𝐿 7 are, respectively, 𝑥 1 = 0 . 9 8 6 1 1 , 𝑥 3 = 0 . 2 7 1 3 8 1 ( 3 . 3 5 ) and lies within the fluid. 4. Linear Stability of the Equilibrium Points In order to study the linear stability of any equilibrium point ( 𝑥 1 0 , 𝑥 2 0 , 𝑥 3 0 ) of an infinitesimal body, we displace it to the position ( 𝑥 1 , 𝑥 2 , 𝑥 3 ) such that  𝑥 1 0 + 𝜉 , 𝑥 2 0 + 𝜂 , 𝑥 3 0  + 𝜁 , ( 4 . 1 ) where 𝜉 , 𝜂 , 𝜁 are small displacements, and then linearize equation ( 2.3 ) to obtain the equations: ̈  𝑈 𝜉 − 2 𝑛 ̇ 𝜂 = 0 𝑥 1 𝑥 1   𝑈 𝜉 + 0 𝑥 1 𝑥 2   𝑈 𝜂 + 0 𝑥 1 𝑥 3  ̇  𝑈 𝜁 , ̈ 𝜂 + 2 𝑛 𝜉 = 0 𝑥 1 𝑥 2   𝑈 𝜉 + 0 𝑥 2 𝑥 2   𝑈 𝜂 + 0 𝑥 2 𝑥 3  ̈  𝑈 𝜁 , 𝜁 = 0 𝑥 1 𝑥 3   𝑈 𝜉 + 0 𝑥 2 𝑥 3   𝑈 𝜂 + 0 𝑥 3 𝑥 3  𝜁 , ( 4 . 2 ) where the partial derivatives are evaluated at the equilibrium points. 4.1. Equilibrium Points Near the Center of the First Primary In order to consider the motion near any equilibrium point in the 𝑥 1 𝑥 2 -plane, we let solutions of the first two equations of ( 4.2 ) be 𝜉 = 𝐴 e x p ( 𝜆 𝑡 ) , 𝜂 = 𝐵 e x p ( 𝜆 𝑡 ) , ( 4 . 3 ) where 𝐴 , 𝐵 , and 𝜆 are constants. Taking first and second derivatives of the above, substituting them into the first two equations of system ( 4.2 ) and has a non-zero solution when | | | | | | |  𝜆 2 − 𝑈 0 𝑥 1 𝑥 1   2 𝑛 𝜆 + 𝑈 0 𝑥 1 𝑥 2   2 𝑛 𝜆 − 𝑈 0 𝑥 1 𝑥 2 𝜆   2 − 𝑈 0 𝑥 2 𝑥 2  | | | | | | | = 0 . ( 4 . 4 ) Expanding the determinant, we have 𝜆 4 −  𝑈 0 𝑥 1 𝑥 1 + 𝑈 0 𝑥 2 𝑥 2 − 4 𝑛 2  𝜆 2 + 𝑈 0 𝑥 1 𝑥 1 𝑈 0 𝑥 2 𝑥 2 −  𝑈 0 𝑥 1 𝑥 2  2 = 0 . ( 4 . 5 ) Equation ( 4.5 ) is the characteristic equation corresponding to the variational equations ( 4.2 ) in the case when motion is considered in the 𝑥 1 , 𝑥 2 -plane. Now, the values of the second-order partial derivatives of the equilibrium point ( 𝑥 𝐿 , 0 , 0 ) , where 𝑥 𝐿 = 𝑝 1 stands for the first equilibrium and 𝑥 𝐿 = 𝑥 1 1 + 𝑝 2 for the second one, are given as 𝑈 0 𝑥 1 𝑥 1  −  = 𝐷 𝜇 1 − 𝑥 𝐿  3 −  3 𝛼 2 / 2   1 − 𝑥 𝐿  + 2 𝑥 𝐿  1 − 𝑥 𝐿  2 + 6 𝑥 𝐿 𝛼 2 + 𝑛 2  1 − 𝑥 𝐿  5 𝑥 𝐿  1 − 𝑥 𝐿  5  , 𝑈 0 𝑥 2 𝑥 2  −  = 𝐷 𝜇 1 − 𝑥 𝐿  3 − ( 3 / 2 ) 𝛼 2  1 − 𝑥 𝐿  − 𝑥 𝐿  1 − 𝑥 𝐿  2 − ( 3 / 2 ) 𝛼 2 𝑥 𝐿 + 𝑛 2  1 − 𝑥 𝐿  5 𝑥 𝐿  1 − 𝑥 𝐿  5  , 𝑈 0 𝑥 3 𝑥 3  = − 𝐷 2 𝜋 𝜌 1 𝐴 2 + 𝜇  1 − 𝑥 𝐿  3 + 9 𝜇 𝛼 2 2  1 − 𝑥 𝐿  5  , 𝑈 0 𝑥 1 𝑥 2 = 0 = 𝑈 0 𝑥 2 𝑥 3 = 𝑈 0 𝑥 1 𝑥 3 . ( 4 . 6 ) Substituting these in ( 4.2 ), we at once have the variational equations: ̈ 𝜉 − 2 𝑛 ̇ 𝜂 = 𝑈 0 𝑥 1 𝑥 1 ̇ 𝜉 , ̈ 𝜂 + 2 𝑛 𝜉 = 𝑈 0 𝑥 2 𝑥 2 ̈  𝜇 𝜂 , ( 4 . 7 ) 𝜁 = − 𝐷  1 − 𝑥 𝐿  3 + 9 𝜇 𝛼 2 2  1 − 𝑥 𝐿  5 + 2 𝜋 𝜌 1 𝐴 2  𝜁 , ( 4 . 8 ) where the partial derivatives have been computed at each equilibrium point 𝑥 𝐿 . Now, ( 4.8 ) is independent of ( 4.7 ), the solution being a periodic function is bounded and therefore, the motion of the infinitesimal body in the 𝑥 3 direction is stable. Now, the characteristic equation of the equilibrium points ( 𝑥 𝐿 , 0 , 0 ) corresponding to the system ( 4.7 ) is 𝜆 4 −  𝑈 0 𝑥 1 𝑥 1 + 𝑈 0 𝑥 2 𝑥 2 − 4 𝑛 2  𝜆 2 + 𝑈 0 𝑥 1 𝑥 1 𝑈 0 𝑥 2 𝑥 2 = 0 , ( 4 . 9 ) where 𝑈 0 𝑥 1 𝑥 1  𝛼 = 3 𝐷 𝜇 1 2 𝑥 𝐿  𝑈 , ( 4 . 1 0 ) 0 𝑥 2 𝑥 2  𝛼 = 3 𝐷 𝜇 − 1 + 1 2 𝑥 𝐿  . ( 4 . 1 1 ) These equations have been obtained using binomial expansion and ignoring terms with second and higher power in 𝑝 1 , 𝑝 2 , 𝛼 2 , and their product. Now, let 𝜆 2 1 and 𝜆 2 2 be the roots of ( 4.9 ), then, the equilibrium point is stable if both the roots are real and negative. This means that their sum must be negative and their product must be positive. Hence, the points ( 𝑥 𝐿 , 0 , 0 ) will be stable if the following two conditions hold: 𝜆 2 1 + 𝜆 2 2 = 𝑈 0 𝑥 1 𝑥 1 + 𝑈 0 𝑥 2 𝑥 2 − 4 𝑛 2 𝜆 < 0 , ( 4 . 1 2 ) 2 1 𝜆 2 2 = 𝑈 0 𝑥 1 𝑥 1 𝑈 0 𝑥 2 𝑥 2 > 0 . ( 4 . 1 3 ) Now, in the case of the first equilibrium point 𝑥 𝐿 = 𝑝 1 , if we suppose in ( 4.10 ) that 𝑝 1 < 0 then, 𝑈 0 𝑥 1 𝑥 1 < 0 since 0 < 𝜇 < 1 , 𝐷 > 0 , 0 < 𝛼 1 ≪ 1 and when 𝑝 1 > 0 , we have 𝑈 0 𝑥 1 𝑥 1 > 0 . Similarly, in ( 4.11 ), if we suppose 𝑝 1 < 0 then 𝑈 0 𝑥 2 𝑥 2 < 0 . For the case 𝑝 1 > 0 , we will have 𝑈 0 𝑥 2 𝑥 2 > 0 when 𝛼 1 > 2 | 𝑝 1 | which is not possible, hence 𝑈 0 𝑥 2 𝑥 2 < 0 . In the case 0 < 𝑝 1 < 𝛼 1 / 2 , 𝑈 0 𝑥 1 𝑥 1 > 0 , and 𝑈 0 𝑥 2 𝑥 2 > 0 . Also, if 0 < 𝛼 1 / 2 < 𝑝 1 , we see that 𝑈 0 𝑥 1 𝑥 1 > 0 and 𝑈 0 𝑥 2 𝑥 2 < 0 . Thus, for the case 𝑝 1 < 0 , the equilibrium point is stable. For 0 < 𝑝 1 < 𝛼 1 / 2 , the equilibrium point is stable if the condition ( 4.12 ) holds. When 0 < 𝛼 1 / 2 < 𝑝 1 , the equilibrium point is unstable. Next, for the other equilibrium point positioned at 𝑥 𝐿 = 𝑥 1 1 + 𝑝 2 , when 𝑥 1 1 > 0 , then 𝑥  1 1 > 0 since | 𝑝 2 | ≪ 1 and the equilibrium point is stable if the conditions ( 4.12 ) and ( 4.13 ) are satisfied. If 𝑥 1 1 < 0 then 𝑥  1 1 < 0 ; it makes 𝑈 0 𝑥 1 𝑥 1 < 0 , 𝑈 0 𝑥 2 𝑥 2 < 0 . Therefore, when 𝑥 1 1 < 0 , both the conditions ( 4.12 ) and ( 4.13 ) are fulfilled and the equilibrium point is stable. 4.2. Circular Points At circular points ( 1 + 𝑟 c o s 𝜃 , 𝑟 s i n 𝜃 , 0 ) , the values of the second partial derivatives with the use of ( 3.14 ) and neglecting the product 𝛼 1 𝛼 2 are 𝑈 0 𝑥 1 𝑥 1 = 3 𝐷 𝜇 c o s 2 𝜃  𝑛 2 + 𝛼 2  , 𝑈 0 𝑥 1 𝑥 2  𝑛 = 3 𝐷 𝜇 c o s 𝜃 s i n 𝜃 2 + 𝛼 2  , 𝑈 0 𝑥 2 𝑥 2 = 3 𝐷 𝜇 s i n 2 𝜃  𝑛 2 + 𝛼 2  , 𝑈 0 𝑥 3 𝑥 3  = − 𝐷 2 𝜋 𝜌 1 𝐴 2  𝑛 + 𝜇 2 + 3 𝛼 2 .   ( 4 . 1 4 ) Substituting these values in the variational equations ( 4.2 ), we get ̈ 𝜉 − 2 𝑛 ̇ 𝜂 = 3 𝐷 𝜇 c o s 2 𝜃  𝑛 2 + 𝛼 2   𝑛 𝜉 + 3 𝐷 𝜇 c o s 𝜃 s i n 𝜃 2 + 𝛼 2  ̇  𝑛 𝜂 + ( 0 ) 𝜁 , ̈ 𝜂 + 2 𝑛 𝜉 = 3 𝐷 𝜇 c o s 𝜃 s i n 𝜃 2 + 𝛼 2  𝜉 + 3 𝐷 𝜇 s i n 2 𝜃  𝑛 2 + 𝛼 2  ̈  𝜂 + ( 0 ) 𝜁 , ( 4 . 1 5 ) 𝜁 = − 𝐷 2 𝜋 𝜌 1 𝐴 2  𝑛 + 𝜇 2 + 3 𝛼 2   𝜁 . ( 4 . 1 6 ) Equation ( 4.16 ) is independent of ( 4.15 ), it shows that the motion of the infinitesimal mass along the 𝑥 3 -direction is stable. Now, a substitution of these partial derivatives in the characteristic equation ( 4.5 ) yields 𝜆 4 −   𝑛 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2  𝜆 2 = 0 . ( 4 . 1 7 ) Let 𝜆 2 = Λ in ( 4.17 ) then, we have Λ    𝑛 Λ − 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2   = 0 . ( 4 . 1 8 ) Hence, either  𝑛 Λ = 0 o r Λ = 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2 , ( 4 . 1 9 ) which implies that   𝑛 𝜆 = 0 t w i c e , o r 𝜆 = ± 3 𝐷 𝜇 2 + 𝛼 2  − 4 𝑛 2  1 / 2 . ( 4 . 2 0 ) Therefore, ( 4.20 ) gives the roots of the characteristic equation ( 4.17 ). Hence we conclude that the circular points are unstable due to the presence of multiple roots. 4.3. Out-of-Plane Points To determine the stability of the out-of-plane equilibrium points, we consider the following partial derivatives: 𝑈 𝑥 1 𝑥 1 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 ⎧ ⎪ ⎨ ⎪ ⎩  2  + 𝜇 1 − 𝑥 1  2 − 𝑥 2 3    1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 3 2 𝜇 𝛼 2 ⎧ ⎪ ⎨ ⎪ ⎩  4  1 − 𝑥 1  2 − 𝑥 2 3    1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ − 1 5 2 𝜇 𝛼 2 𝑥 2 3 ⎧ ⎪ ⎨ ⎪ ⎩  6  1 − 𝑥 1  2 − 𝑥 2 3    1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑈 𝑥 1 𝑥 3 ⎡ ⎢ ⎢ ⎢ ⎣  = 𝐷 − 3 𝜇 1 − 𝑥 1  𝑥 3   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 − 1 5 2 𝜇 𝛼 2  1 − 𝑥 1  𝑥 3   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 − 1 5 𝜇 𝛼 2  1 − 𝑥 1  2 2 𝑥 3   1 − 𝑥 1  2 + 𝑥 2 3  − 7 𝑥 3 3   1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑈 𝑥 2 𝑥 2 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 ⎧ ⎪ ⎨ ⎪ ⎩  − 𝜇 1 − 𝑥 1  2 + 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ − 3 𝜇 𝛼 2 2 ⎧ ⎪ ⎨ ⎪ ⎩  1 − 𝑥 1  2 + 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 1 5 𝜇 𝛼 2 𝑥 2 3 2 ⎧ ⎪ ⎨ ⎪ ⎩  1 − 𝑥 1  2 + 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 𝑛 2 ⎤ ⎥ ⎥ ⎥ ⎦ , 𝑈 𝑥 3 𝑥 3 ⎡ ⎢ ⎢ ⎢ ⎣ = 𝐷 − 2 𝜋 𝜌 1 𝐴 2 ⎧ ⎪ ⎨ ⎪ ⎩  − 𝜇 1 − 𝑥 1  2 − 2 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  5 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ − 9 𝜇 𝛼 2 2 ⎧ ⎪ ⎨ ⎪ ⎩  1 − 𝑥 1  2 − 4 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  7 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ + 1 5 𝜇 𝛼 2 2 ⎧ ⎪ ⎨ ⎪ ⎩ 3 𝑥 2 3   1 − 𝑥 1  2 + 𝑥 2 3  − 7 𝑥 4 3   1 − 𝑥 1  2 + 𝑥 2 3  9 / 2 ⎫ ⎪ ⎬ ⎪ ⎭ ⎤ ⎥ ⎥ ⎥ ⎦ . ( 4 . 2 1 ) Since 𝑥 2 = 0 , therefore the partial derivatives to be computed at the out-of-plane equilibrium points are 𝑈 0 𝑥 1 𝑥 2 = 0 = 𝑈 0 𝑥 2 𝑥 1 = 𝑈 0 𝑥 2 𝑥 3 = 𝑈 0 𝑥 3 𝑥 2 , 𝑈 0 𝑥 1 𝑥 1  = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 − 2 𝜋 𝜌 1 𝐴 2 + 𝜇 𝑙 3 − 3 𝜇 𝑥 2 3 𝑙 5 + 3 𝜇 𝛼 2 2 𝑙 5 − 9 𝜇 𝛼 2 2 𝑙 5 − 4 5 𝜇 𝛼 2 𝑥 2 3 𝑙 7 + 1 0 5 𝜇 𝛼 2 𝑥 4 3 2 𝑙 9 + 𝑛 2  , 𝑈 0 𝑥 2 𝑥 2  = 𝐷 − 2 𝜋 𝜌 1 𝐴 1 + 2 𝜋 𝜌 1 𝐴 2 + 3 𝜇 𝛼 2 𝑙 5 + 𝑛 2  , 𝑈 0 𝑥 1 𝑥 3  = − 3 1 − 𝑥 1  𝑥 3 𝐷  𝜇 𝑙 5 + 1 5 𝜇 𝛼 2 2 𝑙 7 − 3 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 9  , 𝑈 0 𝑥 3 𝑥 3  = 𝐷 3 𝜇 𝑥 2 3 𝑙 5 + 7 5 𝜇 𝛼 2 𝑥 2 3 2 𝑙 7 − 1 0 5 𝜇 𝛼 2 𝑥 4 3 2 𝑙 9  . ( 4 . 2 2 ) Now, we let, 𝑈 0 𝑥 1 𝑥 1 = 𝑈 1 1 , 𝑈 0 𝑥 2 𝑥 2 = 𝑈 2 2 , 𝑈 0 𝑥 1 𝑥 3 = 𝑈 1 3 , 𝑈 0 𝑥 3 𝑥 3 = 𝑈 3 3 , 𝑈 0 𝑥 1 𝑥 2 = 𝑈 1 2 , 𝑈 0 𝑥 2 𝑥 3 = 𝑈 2 3 . ( 4 . 2 3 ) Using ( 4.23 ), the variational equation can be recast in the form: ̈ 𝜉 − 2 𝑛 ̇ 𝜂 = 𝑈 1 1 𝜉 + 𝑈 1 3 ̇ 𝜁 , ̈ 𝜂 + 2 𝑛 𝜉 = 𝑈 2 2 𝜂 + 𝑈 2 3 ̈ 𝜁 , 𝜁 = 𝑈 1 3 𝜉 + 𝑈 3 3 𝜁 . ( 4 . 2 4 ) In order to consider the motion of the out-of-plane points, we let solution of the system ( 4.24 ) be 𝜉 = 𝐴 e x p ( 𝜆 𝑡 ) , 𝜂 = 𝐵 e x p ( 𝜆 𝑡 ) , 𝜁 = 𝐶 e x p ( 𝜆 𝑡 ) , ( 4 . 2 5 ) where 𝐴 , 𝐵 , 𝐶 , and 𝜆 are constants. 𝜉 , 𝜂 , and 𝜁 are the small displacements in the coordinates of the infinitesimal body. Now, the characteristic equation corresponding to the variational equations ( 4.24 ) in the case of the out-of-plane point may be expressed as 𝜆 6 − 𝑎 1 𝜆 4 + 𝑎 2 𝜆 2 + 𝑎 3 = 0 , ( 4 . 2 6 ) where the coefficients of the characteristic equation ( 4.26 ) are such that 𝑎 1  = 𝐷 − 4 𝜋 𝜌 1 𝐴 1 − 2 𝜋 𝜌 1 𝐴 2 + 2 𝑛 2  , 𝑎 2 = 1 4 𝑙 1 6  𝐷  − 9 𝜇 2 𝐷 𝑥 2 3  − 1 + 𝑥 1  2  2 𝑙 4 + 5 𝛼 2  3 𝑙 2 − 7 𝑥 2 3   + 3 𝑙 4 𝜇 𝑥 2 3  2 𝑙 4 + 5 𝛼 2  5 𝑙 2 − 7 𝜇 𝑥 2 3 ×    3 𝐷 𝜇 𝛼 2 + 2 𝑙 5  ( − 4 + 𝐷 ) 𝑛 2 − 2 𝐷 𝐴 1 𝜋 𝜌 1 + 2 𝐷 𝐴 2 𝜋 𝜌 1    + 𝐷 6 𝑙 4 𝜇 𝑥 2 3 + 3 𝛼 2  𝑙 4 𝜇 − 2 5 𝑙 2 − 3 5 𝜇 𝑥 4 3  + 2 𝑙 9  𝑛 2 − 2 𝐴 1 𝜋 𝜌 1 + 2 𝐴 2 𝜋 𝜌 1 ×    3 𝛼 2  𝑙 4 𝜇 − 2 𝑙 2  1 5 + 𝑙 2  𝑥 2 3 + 3 5 𝜇 𝑥 4 3  + 2 𝑙 6  𝑙 3 𝑛 2 + 𝜇 − 2 𝑙 3 𝐴 1 𝜋 𝜌 1 − 2 𝑙 3 𝐴 2 𝜋 𝜌 1 𝑎     3 = 3 𝐷 3 4 𝑙 1 8 𝜇 𝑥 2 3  2 𝑙 4 + 5 𝛼 2  5 𝑙 2 − 7 𝜇 𝑥 2 3  𝑛   2 + 3 𝜇 𝛼 2 2 𝑙 5 − 2 𝐴 1 𝜋 𝜌 1 + 2 𝐴 2 𝜋 𝜌 1  ×  3 𝛼 2  𝑙 4 𝜇 − 2 𝑙 2  1 5 + 𝑙 2  𝜇 𝑥 2 3 + 3 5 𝜇 𝑥 4 3  + 2 𝑙 6  𝑙 3 𝑛 2 + 𝜇 − 2 𝑙 3 𝐴 1 𝜋 𝜌 1 − 2 𝑙 3 𝐴 2 𝜋 𝜌 1 ,   ( 4 . 2 7 ) where 𝑙 2 = ( 1 − 𝑥 1 ) 2 + 𝑥 2 3 . These computations have been done using the software package Mathematica. For the stability analysis of the out-of-plane equilibrium point, we compute numerically the partial derivatives calculated at the out-of-plane points with the use of ( 3.28 ) and the following numerical values: 𝜇 = 0 . 0 1 , 𝜋 = 3 . 1 4 , 𝛼 1 = 0 . 0 2 4 , 𝛼 2 𝐴 = 0 . 0 2 1 = 0 . 7 , 𝐴 2 = 0 . 6 8 , 𝜌 1 = 0 . 2 3 6 , 𝐷 = 0 . 2 1 3 3 . ( 4 . 2 8 ) Now, substituting the above values in the characteristic equation ( 4.27 ), we get 𝜆 6 − 0 . 1 9 2 4 3 7 𝜆 4 + 0 . 2 4 8 0 3 4 𝜆 2 − 0 . 0 0 0 0 1 9 7 6 2 4 = 0 . ( 4 . 2 9 ) Its roots are: 𝜆 1 , 2 𝜆 = − 0 . 5 4 5 0 6 6 ± 0 . 4 4 8 2 3 9 𝑖 , 3 , 4 𝜆 = ± 0 . 0 0 8 9 2 6 4 2 , 5 , 6 = 0 . 5 4 5 0 6 6 ± 0 . 4 4 8 2 3 9 𝑖 . ( 4 . 3 0 ) The positive root and the positive real part of the complex roots induce instability at the out-of-plane point. Hence, the motion of the infinitesimal mass around the out-of-plane equilibrium points is unstable for the specific numerical example given here. However, fuller discussion of their stability remains a theme for future research. 5. Discussion The equation of motion ( 2.3 ) is different from those of Hallan and Mangang [ 4 ] due to oblateness of the second primary. If we assume that the second primary is not oblate (i.e., 𝛼 2 = 0 ), then these equations will fully coincide with those of Hallan and Mangang [ 4 ]. Equation ( 3.9 ) gives the equilibrium position of the point ( 𝑃 1 , 0 , 0 ) near the center of the first primary and fully coincides with that of Hallan and Mangang [ 4 ]. It shows that the position of this equilibrium point does not depend on oblateness of the second primary, while the other equilibrium point ( 𝑥 1 1 + 𝑃 2 , 0 , 0 ) given by ( 3.11 ) is different from that of Hallan and Mangang [ 4 ] due to the appearance of oblateness of the second primary. When 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) , points on the circle ( 1 − 𝑥 1 ) 2 + 𝑥 2 2 = 𝑟 2 , 𝑥 3 = 0 lying within the first primary are also equilibrium points. These points are affected by oblateness of both primaries. Equations ( 3.28 ) give the positions of the out-of-plane points when only linear terms in oblateness of the second primary are retained. We have been able to show that the oblateness of the primaries allows the existence of the out-of-plane equilibrium points in the 𝑥 1 𝑥 3 -plane within the first primary. These points have no analogy in the previous studies of the Robe’s restricted three-body problem. The linear stability analysis of the equilibrium solutions of the problem is investigated with the help of characteristic roots. The characteristic equation ( 4.9 ) in the case of the equilibrium point 𝑥 𝐿 = 𝑝 1 near the center of the first primary is the same as that of Hallan and Mangang [ 4 ], while that of the other point 𝑥 𝐿 = 𝑥 1 1 + 𝑝 2 near the center differs from that of Hallan and Mangang [ 4 ] due to oblateness of the second primary. The characteristic equation of the circular case ( 4.17 ) also differs from that of Hallan and Mangang [ 4 ] due to oblateness of the second primary. The stability in the first approximation of this configuration shows that points near the centre of the first primary are conditionally stable; the circular points are unstable. This confirms the earlier results of Hallan and Rana [ 3 ], Hallan and Mangang [ 4 ]. A numerical exploration shows that the out-of-plane equilibrium points are also unstable. This outcome validates the earlier results of Douskos and Markellos [ 8 ] and Singh and Leke [ 11 ] that the points are unstable. 6. Conclusion We have derived the equations of motion and established the positions of the equilibrium points of the infinitesimal body in the Robe’s [ 1 ] restricted three-body problem with oblateness. The term “oblateness” is used in the sense that both primaries are considered as oblate spheroids under the effects of the full buoyancy force exerted by the fluid on the infinitesimal mass. We have obtained one equilibrium point ( 𝑃 1 , 0 , 0 ) near the centre of the first primary which will be on the left or right of the centre of the first primary accordingly as 2 𝜋 𝜌 1 𝐴 1 − 2 𝜇 > < 1 . This point is the same as that of Hallan and Mangang [ 4 ]. In addition to this, another equilibrium point ( 𝑥 1 1 + 𝑃 2 , 0 , 0 ) is found within the first primary on the line joining the center of the primaries when 1 − 2 𝜋 𝜌 1 𝐴 1 < − 3 𝜇 / 4 and | 𝑥 1 1 | < 𝑎 1 . When 2 𝜋 𝜌 1 𝐴 1 = 𝑛 2 ( 1 − 𝜇 ) , points on the circle ( 1 − 𝑥 1 ) 2 + 𝑥 2 2 = 𝑟 2 , 𝑥 3 = 0 lying within the first primary are also equilibrium points. We call them circular points. Finally, we have been able to show that the oblateness of the primaries allows the existence of the out-of-plane equilibrium points in the 𝑥 1 𝑥 3 -plane within the first primary. The result of this paper can be summarized as follows. The restricted three-body problem under the framework of the Robe’s [ 1 ] problem with oblate primaries has the equilibrium points of the type: points near the center of the first primary, points on the circle (circular points), and two out-of-plane points 𝐿 6 , 7 . It is seen that points near the first primary are conditionally stable, the circular points are unstable, while the out-of-plane equilibrium points are unstable for the specific numerical example given here. The effect of drag forces as considered by Giordano et al. [ 2 ] under the present context, particularly as regards the analysis of the properties of the equilibrium points located inside the first primary, will be interesting. References H. A. G. Robe, “A new kind of three body problem,” Celestial Mechanics and Dynamical Astronomy , vol. 16, pp. 343–351, 1977. C. M. Giordano, A. R. Plastino, and A. Plastino, “Robe's restricted three-body problem with drag,” Celestial Mechanics and Dynamical Astronomy , vol. 66, no. 2, pp. 229–242, 1996. View at Publisher · View at Google Scholar P. P. Hallan and N. Rana, “The existence and stability of equilibrium points in the Robe's restricted three-body problem,” Celestial Mechanics & Dynamical Astronomy , vol. 79, no. 2, pp. 145–155, 2001. View at Publisher · View at Google Scholar P. P. Hallan and K. B. Mangang, “Existence and linear stability of equilibrium points in the Robe's restricted three body problem when the first primary is an oblate spheroid,” Planetary and Space Science , vol. 55, no. 4, pp. 512–516, 2007. View at Publisher · View at Google Scholar · View at Scopus A. R. Plastino and A. Plastino, “Robe's restricted three-body problem revisited,” Celestial Mechanics and Dynamical Astronomy , vol. 61, no. 2, pp. 197–206, 1995. View at Publisher · View at Google Scholar · View at Scopus R. K. Sharma and P. V. Subba Rao, “Stationary solutions and their characteristic exponents in the restricted three-body problem when the more massive primary is an oblate spheroid,” Celestial Mechanics , vol. 13, no. 2, pp. 137–149, 1976. J. Singh and B. Ishwar, “Stability of triangular points in the generalized photogravitational restricted three-body problem,” Bulletin of the Astronomical Society of India , vol. 27, pp. 415–424, 1999. C. N. Douskos and V. V. Markellos, “Out-of-plane equilibrium points in the restricted three-body problem with oblateness,” Astronomy and Astrophysics , vol. 446, no. 1, pp. 357–360, 2006. View at Publisher · View at Google Scholar · View at Scopus J. Singh and O. Leke, “Stability of the photogravitational restricted three-body problem with variable masses,” Astrophysics and Space Science , vol. 326, pp. 305–314, 2010. V. V. Radzievskii, “The restricted problem of three bodies taking account of light pressure,” Astronomicheskii Zhurnal , vol. 27, pp. 250–256, 1950. J. Singh, “Combined effects of perturbations, radiation and oblateness on the nonlinear stability of triangular points in the restricted three-body problem,” Astrophysics and SpaceScience , vol. 332, pp. 331–339, 2011. //

Journal

Advances in Mathematical PhysicsHindawi Publishing Corporation

Published: Aug 16, 2012

There are no references for this article.