# Generalized Out-of-Plane Equilibrium Points in the Frame of Elliptic Restricted Three-Body Problem: Impact of Oblate Primary and Luminous-Triaxial Secondary

Authors: Aminu Abubakar Hussain; Aishetu Umar
Journal: Advances in Astronomy
 (2019)

---

## Abstract

This paper studies the motion of a third body near the1st family of the out-of-plane equilibrium points, L6,7, in the elliptic restricted problem of three bodies under an oblate primary and a radiating-triaxial secondary. It is seen that the pair of points (ξ0,0,±ζ0) which correspond to the positions of the 1st family of the out-of-plane equilibrium points, L6,7, are affected by the oblateness of the primary, radiation pressure and triaxiality of the secondary, semimajor axis, and eccentricity of the orbits of the principal bodies. But the point ±ζ0 is unaffected by the semimajor axis and eccentricity of the orbits of the principal bodies. The effects of the parameters involved in this problem are shown on the topologies of the zero-velocity curves for the binary systems PSR 1903+0327 and DP-Leonis. An investigation of the stability of the out-of-plane equilibrium points, L6,7 numerically, shows that they can be stable for 0.32≤μ≤0.5 and for very low eccentricity. L6,7 of PSR 1903+0327 and DP-Leonis are however linearly unstable.


---

## Body

## 1. Introduction

Generally, the motion of a body is influenced by the mutual gravitational attraction from other bodies. In a situation where we consider an isolated dynamical system of three gravitational interacting bodiesm1,m2, and m3, and where the third body m3 (called the infinitesimal mass) has a very small mass compared to the other bodies m1andm2 (called the principal bodies, where m1≥m2), the third body m3 has a very small effect (negligible) on the motion of the principal bodies m1andm2 (i.e., m3≪m1&m2). These principal bodies called the primaries (i.e., m1 is called bigger primary or primary and m2 is called smaller primary or secondary, respectively) act as a pair, which has their motion approximated by a two-body problem with their trajectories predetermined. The motion of the third body as affected by the principal bodies is the subject of the restricted three-body problem (R3BP). This is what made the R3BP a more tractable case than the general three-body problem. If the principal bodies are restricted to a circular or elliptic orbit while the third body of infinitesimal mass orbits in the plane of motion established by the principal bodies, then we, respectively, have circular or elliptic R3BP (CR3BP or ER3BP). See [1].Over decades, the study of the existence of some families of particular solutions in both circular and elliptic R3BP has received the attention of various researchers like [1–13] among others.Reference [14] was the first to point out the existence of the out-of-plane equilibrium points in the cases of Sun-Planet-Particle and Galaxy Kernel-Sun-Particle and he found the two equilibrium points L6,7 on the x-z plane symmetrically with respect to the orbital plane. The case where the Poynting-Robertson effect is taken into account was treated by [15].Reference [16] proved the existence of a second family of the out-of-plane equilibrium points L8,9 when the principal bodies are both luminous under certain condition (i.e., the relation between the mass and radiation pressure parameters). They also studied the stability for μ=0.5 and μ=0.2.The periodic solution about the out-of-plane equilibrium points in the photogravitational restricted three-body problem was carried out by [17]. They found that two such families exist: the first family originates and terminates on the same equilibrium point, while the second family terminates by flattening on the orbital plane.Afterwards, various studies of the out-of-plane equilibrium points have been carried out by other researchers like [18–23] and many more.Reference [5] studied the motion in the photogravitational elliptic restricted three-body problem under an oblate primary. In their study, they found that the out-of-plane equilibrium points are unstable for any combination of the parameters under their own consideration. Also, [24] studied the family of out-of-plane equilibrium points L6,7 in the elliptic restricted three-body problem with radiating and oblate primaries. They concluded that the positions and stability depend on the oblateness and radiation pressure of the primaries and the eccentricity of their orbits.Reference [25] examined the out-of-plane equilibrium points in the photogravitational CR3BP with oblateness and P-R drag. In their numerical study, they found that the out-of-plane equilibrium points are unstable in the sense of Lyapunov due to the presence of positive real root.Recently, [26] investigated the influence of Poynting-Robertson drag and oblateness on the existence and stability of the out-of-plane equilibrium points in the spatial elliptic restricted three-body problem. They also found L6,7 to be unstable. Reference [27] studied the Newton-Raphson basins of convergence of the out-of-plane equilibrium points in the Copenhagen problem with oblate primaries. He considered a multivariate Newton-Raphson iterative scheme for revealing the corresponding basins of convergence on x-y and y-z planes. He also demonstrated how the oblateness coefficient A influenced the positions of out-of-plane equilibrium points.Our aim in this paper is to study the existence and stability of the out-of-plane equilibrium points under an oblate primary and a luminous-triaxial secondary in the ER3BP with application to the binary systems DP-Leonis and PSR 1903-0327.This article is organized in 7 sections as follows: the first section is the introduction; the equations of motion are given in Section2; Section 3 contains the surface of zero-velocity curves (ZVC) in the (ξ-ζ) plane; the positions and stability of the out-of-plane equilibrium points are treated in Sections 4 and 5, respectively; Section 6 contains numerical application, while Section 7 is the discussion which concludes the study.


## 2. Equations of Motion

The equations of motion of a test particle in the ER3BP with a bigger oblate primary and a radiating-triaxial secondary in dimensionless-pulsating (rotating) coordinate systemξ,η,ζ are given in [5, 6] as follows:(1)ξ′′-2η′=Ωξη′′+2ξ′=Ωηζ′′=Ωζwhere the prime in equations (1) denotes differentiation with respect to the eccentric anomaly E (E is an angular parameter that defines the position of a body that is moving along the elliptic Keplerian orbit; see [5] for more details) and Ω is the potential-like function defined by(2)Ω=1-e2-1/2ξ2+η22+1n21-μr1+μq2r2+1-μA2r13-31-μζ2A2r15+μq22σ1-σ22r23-3μq22σ1-σ2ζ22r25.With r12=ξ+μ2+η2+ζ2 and r22=ξ+μ-12+η2+ζ2 as the distances of the third body m3 from the primary body m1 and secondary body m2, respectively. Also e and a are the eccentricity and semimajor axis of the orbit of the primaries, respectively, n is the perturbed mean motion of the primaries, while q2 is the radiation pressure factor of the secondary body which is given by Fg-Fp=Fg1-Fp/Fg=q2Fg, where q2=1-Fp/Fg such that 0<1-q2=δ≪1 (Radzievskii [14]), and Fg and Fp are, respectively, the gravitational and radiation pressure forces. Additionally, A=a12-a32/5R2≪1, σ1=b12-b32/5R2≪1, and σ2=b22-b32/5R2≪1, where a1 and a3 are, respectively, the equatorial and polar radii of the primary body, b1,b2, and b3 are the lengths of the semiaxes of the secondary body, and R is the dimensional distance between the primaries.


### 2.1. Jacobi Constant

Equations (1) admit the Jacobi constant:(3)2Ω-ξ′2+η′2+ζ′2=C,where C is the constant of integration known as Jacobi constant.


### 2.2. Derivation of the Perturbed Mean Motion

In the case of elliptic orbit, the distance between the primaries isr=a1-e2/1+ecos⁡ν and the mean distance between them is given by(4)12π∫02πrdν=a1-e21+e21/2.The orbits of m1 and m2 with respect to the centre of mass, with semimajor axes a1=am2 and a2=am1, respectively, have the same eccentricity with (see Szebehely [1]):Their equations are given by(5)n2am11-e21+e21/2=k2m21+3A2+32σ1-σ22,n2am21-e21+e21/2=k2m11+3A2+32σ1-σ22,where k is the Gaussian constant of gravitation.Adding equations (5) together we have(6)n2=k21+e21/2a1-e21+3A2+32σ1-σ22.Here, we adopt that the sum of the masses of the primaries and the distances between them as the units of mass and length. We choose the unit of time such as making k2=1. Hence, (7)n2=1+e21/2a1-e21+3A2+32σ1-σ22.Making the simplification of the above n2, by rejecting the second and higher order terms of e2, we get(8)n2=1a1+3e22+3A2+32σ1-σ22.


## 2.1. Jacobi Constant

Equations (1) admit the Jacobi constant:(3)2Ω-ξ′2+η′2+ζ′2=C,where C is the constant of integration known as Jacobi constant.


## 2.2. Derivation of the Perturbed Mean Motion

In the case of elliptic orbit, the distance between the primaries isr=a1-e2/1+ecos⁡ν and the mean distance between them is given by(4)12π∫02πrdν=a1-e21+e21/2.The orbits of m1 and m2 with respect to the centre of mass, with semimajor axes a1=am2 and a2=am1, respectively, have the same eccentricity with (see Szebehely [1]):Their equations are given by(5)n2am11-e21+e21/2=k2m21+3A2+32σ1-σ22,n2am21-e21+e21/2=k2m11+3A2+32σ1-σ22,where k is the Gaussian constant of gravitation.Adding equations (5) together we have(6)n2=k21+e21/2a1-e21+3A2+32σ1-σ22.Here, we adopt that the sum of the masses of the primaries and the distances between them as the units of mass and length. We choose the unit of time such as making k2=1. Hence, (7)n2=1+e21/2a1-e21+3A2+32σ1-σ22.Making the simplification of the above n2, by rejecting the second and higher order terms of e2, we get(8)n2=1a1+3e22+3A2+32σ1-σ22.


## 3. Surface of Zero-Velocity Curves (ZVC) in the (ξ-ζ) Plane

We introduce the ZVC as our quantitative method of understanding information about the motion near out-of-plane equilibrium points without actually solving the differential equations. The potentials of the principal bodies under consideration in this study have significant effects on the existence and the locations of the out-of-plane equilibrium points and, also, on the structure of the boundary between region of possible motion and forbidden region. So, we present the ZVC in the (ξ-ζ) plane in this section for the binary systems under consideration (i.e., PSR 1903+0327 and DP-Leonis. See Section 6 for more details of the binaries). In Figures 1–4, we demonstrate the possible topologies of the curves for PSR 1903+0327, while Figures 5–10 demonstrate the possible topologies of the curves for DP-Leonis. For clarity, in each case we show the 3D (i.e., the right frame in each row) representation of the curves corresponding to the left frame in each row. In the case of PSR 1903+0327, Figures 1 and 2 and Figures 3 and 4 show the topologies of the ZVC for varying the values of A and (σ1&σ2), respectively, when keeping remaining parameters constant. For the system DP-Leonis, Figures 5 and 6, Figures 7 and 8, and Figures 9 and 10 show the topologies of the ZVC for varying values of A, (σ1&σ2), and e, respectively, when keeping remaining parameters constant.Figure 1
Zero-velocity curves in the (ξ-ζ) plane of PSR 1903+0327 for A=0.01,a=0.95,σ1=0.03andσ2=0.015.
(a) (b)Figure 2
Zero-velocity curves in the (ξ-ζ) plane of PSR 1903+0327 for A=0.1,a=0.95,σ1=0.03andσ2=0.015.
(a) (b)Figure 3
Zero-velocity curves in the (ξ-ζ) plane of PSR 1903+0327 for A=0.02,a=0.95,σ1=0.01andσ2=0.005.
(a) (b)Figure 4
Zero-velocity curves in the (ξ-ζ) plane of PSR 1903+0327 for A=0.02,a=0.95,σ1=0.01andσ2=0.005.
(a) (b)Figure 5
Zero-velocity curves in the (ξ-ζ) plane of DP-Leonis for A=0.01,e=0.2,a=0.95,σ1=0.03andσ2=0.015.
(a) (b)Figure 6
Zero-velocity curves in the (ξ-ζ) plane of DP-Leonis for A=0.1,e=0.2,a=0.95,σ1=0.03andσ2=0.015.
(a) (b)Figure 7
Zero-velocity curves in the (ξ-ζ) plane of DP-Leonis for A=0.02,e=0.2,a=0.95,σ1=0.01andσ2=0.005.
(a) (b)Figure 8
Zero-velocity curves in the (ξ-ζ) plane of DP-Leonis for A=0.02,e=0.2,a=0.95,σ1=0.04andσ2=0.02.
(a) (b)Figure 9
Zero-velocity curves in the (ξ-ζ) plane of DP-Leonis for A=0.02,e=0.1,a=0.95,σ1=0.03andσ2=0.015.
(a) (b)Figure 10
Zero-velocity curves in the (ξ-ζ) plane of DP-Leonis for A=0.02,e=0.5,a=0.95,σ1=0.03andσ2=0.015.
(a) (b)


## 4. Positions of the Out-of-Plane Equilibrium Points

The solutions of the systemΩξ=Ωη=Ωζ=0 when (ξ≠0,η=0andζ≠0) represent the positions of the out-of-plane equilibrium points (ξ0,0,ζ0); that is,(9)Ωξ=1-e2-1/2ξ-1n21-μξ+μr13+μq2ξ+μ-1r23+31-μξ+μA2r15+3μq2ξ+μ-12σ1-σ22r25-151-μξ+μζ2A2r17-15μq2ξ+μ-12σ1-σ2ζ22r27=0,(10)Ωη=1-e2-1/2η1-1n21-μr13+μq2r23+31-μA2r15+3μq24σ1-3σ22r25-151-μAζ22r17-15μq22σ1-σ2ζ22r27=0,(11)Ωζ=1-e2-1/2-ζn21-μr13+μq2r23+91-μA2r15+9μq22σ1-σ22r25-151-μAζ22r17-15μq22σ1-σ2ζ22r27=0.From (10) with η=0 we have(12)n2-1-μr13-μq2r23-31-μA2r15-3μq24σ1-3σ22r25+151-μAζ22r17+15μq22σ1-σ2ζ22r27≠0.Also, from (11) with ζ≠0 we have(13)1-μr13+μq2r23+91-μA2r15+9μq22σ1-σ22r25-151-μAζ22r17-15μq22σ1-σ2ζ22r27=0.Here, we let Q1=1-μ and Q2=μq2; then (13) becomes(14)Q1r13+Q2r23+9Q1A2r15+9Q22σ1-σ22r25-15Q1Aζ22r17-15Q22σ1-σ2ζ22r27=0,while (9) becomes(15)n2ξ-Q1ξ+μr13-Q2ξ+μ-1r23-3Q1ξ+μA2r15-3Q2ξ+μ-12σ1-σ22r25+15Q1ξ+μζ2A2r17+15Q2ξ+μ-12σ1-σ2ζ22r27=0.This implies that(16)
ξ
1
-
1
n
2
Q
1
r
1
3
+
Q
2
r
2
3
+
3
Q
1
A
2
r
1
5
-
15
Q
1
ζ
2
A
2
r
1
7
+
3
Q
2
2
σ
1
-
σ
2
2
r
2
5
-
15
Q
2
2
σ
1
-
σ
2
ζ
2
2
r
2
7
-
μ
n
2
Q
1
r
1
3
+
Q
2
r
2
3
+
3
Q
1
A
2
r
1
5
-
15
Q
1
ζ
2
A
2
r
1
7
+
3
Q
2
2
σ
1
-
σ
2
2
r
2
5
-
15
Q
2
2
σ
1
-
σ
2
ζ
2
2
r
2
7
+
1
n
2
Q
2
r
2
3
+
3
Q
2
2
σ
1
-
σ
2
2
r
2
5
-
15
Q
2
2
σ
1
-
σ
2
ζ
2
2
r
2
7
=
0
.Making use of (14) in (16), we have(17)ξ=Q1/r13+9Q1A/2r15+3Q2Q12σ1-σ2/r25-15Q1Aζ2/2r17-3μQ1A/r15n2+3Q1A/2r15+3Q22σ1-σ2/2r25.From (14), we have(18)ζ2=2r17r2715Q1Ar27+15Q22σ1-σ2r17Q1r13+9Q1A2r15+Q2r23+9Q22σ1-σ22r25.Here, we use ξ0=1-μ and ζ0=3A as initial approximation. Using Wolfam Mathematica 10.3 (software package) we are able to obtain the positions of a pair of points ξ0,0,±ζ0 corresponding to out-of-plane equilibrium points L6,7 which can be approximated in the forms of power series to the third order terms in A from (17) and (18) as follows:(19)ξ0=1-μ-331-μ21-a+3e2+32σ1-σ2A5/22aμq22σ1-σ2+0A7/2ζ0=35A5+5A3/252σ1-σ2-5A5/2302σ1-σ22+315μ-1A35μq22σ1-σ22+0A7/2The effects of the parameters involved in the positions of the 1st family of out-of-plane equilibrium points are shown graphically for the binary systems PSR 1903+0327 and DP-Leonis, respectively, in Figures 11–14. Figures 11 and 13 show the positions of L6,7 in the ξ-ζ plane as a function of oblateness A in the interval [0,0.1], keeping the remaining parameters constant for the binary systems PSR 1903+0327 and DP-Leonis, respectively, while Figures 12 and 14 show the positions of L6,7 in the ξ-ζ plane as a function of triaxiality σ1&σ2 in the intervals [0.01, 0.02] and [0.005, 0.01], respectively, keeping remaining parameters constant for the binary systems PSR 1903+0327 and DP-Leonis, respectively.Figure 11
Positions ofL6,7 in the (ξ-ζ) plane as a function of oblateness A in the interval [0,0.1] of PSR 1903+0327 for a=0.95,σ1=0.03andσ2=0.015.Figure 12
Positions ofL6,7 in the (ξ-ζ) plane as a function of triaxiality σ1andσ2 in the interval [0.01,0.02]and[0.005,0.01], respectively, of PSR 1903+0327 for a=0.95,A=0.1.Figure 13
Positions ofL6,7 in the (ξ-ζ) plane as a function of oblateness A in the interval [0,0.1] of DP-Leonis for a=0.95,e=0.2,σ1=0.03andσ2=0.015.Figure 14
Positions ofL6,7 in the (ξ-ζ) plane as a function of triaxiality σ1andσ2 in the interval [0.01,0.02]and[0.005,0.01], respectively, of DP-Leonis for a=0.95,e=0.2andA=0.1.


## 5. Linear Stability of the Out-of-Plane Equilibrium Points

In order to study the stability of motion near the out-of-plane equilibrium point, we adopt the characteristic equation of the system as in [4, 21], given as follows:(20)λ6+4-Ωξξ0-Ωηη0-Ωζζ0λ4+Ωηη0Ωζζ0+Ωξξ0Ωζζ0+Ωξξ0Ωηη0-4Ωζζ0-Ωξζ02λ2-Ωξξ0Ωηη0Ωζζ0-Ωξζ02Ωηη0=0.where the superscript 0 denotes that the partial derivatives are evaluated at the out-of-plane point ξ0,0,ζ0.So, at the out-of-plane equilibrium points we have(21)
Ω
ξ
ξ
0
=
1
-
e
2
-
1
/
2
1
+
1
n
2
3
Q
1
ξ
0
+
μ
2
r
10
5
-
Q
1
r
10
3
+
3
Q
2
ξ
0
+
μ
-
1
2
r
20
5
-
Q
2
r
20
3
+
15
Q
1
ξ
0
+
μ
2
A
2
r
10
7
-
3
Q
1
A
2
r
10
5
-
105
Q
1
ξ
0
+
μ
2
ζ
0
2
A
2
r
10
9
+
15
Q
1
ζ
0
2
A
2
r
10
7
+
15
Q
2
ξ
0
+
μ
-
1
2
2
σ
1
-
σ
2
2
r
20
7
-
3
Q
2
2
σ
1
-
σ
2
2
r
20
5
-
105
Q
2
ξ
0
+
μ
-
1
2
2
σ
1
-
σ
2
ζ
0
2
2
r
20
9
+
15
Q
2
2
σ
1
-
σ
2
ζ
0
2
2
r
20
7
,
(22)
Ω
η
η
0
=
1
-
e
2
-
1
/
2
1
-
1
n
2
Q
1
r
10
3
+
Q
2
r
20
3
+
3
Q
1
A
2
r
10
5
-
15
Q
1
ζ
0
2
A
2
r
10
7
+
3
Q
2
4
σ
1
-
3
σ
2
2
r
20
5
-
15
Q
2
2
σ
1
-
σ
2
ζ
0
2
2
r
20
7
,
(23)
Ω
ζ
ζ
0
=
1
-
e
2
-
1
/
2
1
n
2
3
Q
1
ζ
0
2
r
10
5
-
Q
1
r
10
3
+
3
Q
2
ζ
0
2
r
20
5
-
Q
2
r
20
3
+
45
Q
1
ζ
0
2
A
2
r
10
7
-
9
Q
1
A
2
r
10
5
-
105
Q
1
ζ
0
4
A
2
r
10
9
+
45
Q
2
2
σ
1
-
σ
2
ζ
0
2
2
r
20
7
-
9
Q
2
2
σ
1
-
σ
2
2
r
20
5
-
105
Q
2
2
σ
1
-
σ
2
ζ
0
4
2
r
20
9
,
(24)
Ω
ξ
ζ
0
=
1
-
e
2
-
1
/
2
3
ζ
0
n
2
3
Q
1
ξ
0
+
μ
r
10
5
+
Q
2
ξ
0
+
μ
-
1
r
20
5
+
15
Q
1
ξ
0
+
μ
A
2
r
10
7
-
35
Q
1
ξ
0
+
μ
ζ
0
2
A
2
r
10
9
+
15
Q
2
ξ
0
+
μ
-
1
2
σ
1
-
σ
2
2
r
20
7
-
35
Q
2
ξ
0
+
μ
-
1
2
σ
1
-
σ
2
ζ
0
2
2
r
20
9
.


## 6. Numerical Application

In this section, we apply the data (barrowed from NASA ADS and [7]) of the binaries PSR 1903+0327 and DP-Leonis (Table 1) to the analytical results obtained in the previous sections under the assumption that the primary body is an oblate spheroid while the secondary body is a luminous-triaxial body both moving in elliptic Keplerian orbits. The positions of the out-of-plane equilibrium points (Equation (19)) together with the roots of the characteristic equation (Equation (20)) are computed (for some assumed values of the potential of the system under consideration) in Tables 2–4 for PSR 1903+0327 and in Tables 5–7 for DP-Leonis with the help of the Wolfram Mathematica 10.3.Table 1
Relevant numerical data.
Binary System Masses (M⊙) Luminosity Mass ratio Radiation pressure Eccentricity m 1 m 2 (L⊙) (μ) (q2) (e) PSR 1903+0327 1.6700 1.1000 2.5000 0.397112 0.997581 0.437  DP-Leonis 0.600 0.090 0.00449 0.130435 0.999947 -Table 2
Effects of oblateness of PSR 1903+0327 on the locations and stability of out-of-plane points fora=0.95,e=0.437,σ1=0.03andσ2=0.015.
A ξ ± ζ λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.01 0.602141 0.143655 ±39.12118 ± 98.95687 i ± 35.94876 unstable  0.02 0.598661 0.215135 - 12.78173 ± 0.6731746 i ± 36.07696 i 12.78173 ± 0.6731746 i unstable  0.03 0.591240 0.276159 - 6.249728 ± 0.8922171 i ± 19.32210 i 6.249728 ± 0.8922171 i unstable  0.04 0.578977 0.331025 - 3.497018 ± 0.9379422 i ± 12.24203 i 3.497018 ± 0.9379422 i unstable  0.05 0.561117 0.380711 - 2.031369 ± 0.9503968 i ± 8.546200 i 2.031369 ± 0.9503968 i unstableTable 3
Effects of eccentricity of PSR 1903+0327 on the locations and stability of out-of-plane points fora=0.95,A=0.02,σ1=0.03andσ2=0.015.
e ξ ± ζ λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.1 0.602643 0.143655 ±37.44600 ± 94.78687 i ±34.45998 unstable  0.2 0.602560 0.143655 ±37.72920 ± 95.49115 i ±34.71129 unstable  0.3 0.602421 0.143655 ±38.19624 ± 96.65320 i ±35.12607 unstable  0.4 0.602227 0.143655 ±38.84000 ± 98.25624 i ±35.69849 unstable  0.5 0.601977 0.143655 ±39.65119 ± 100.2783 i ±36.42090 unstableTable 4
Effects of triaxiality of PSR 1903+0327 on the locations and stability of out-of-plane points fora=0.95,e=0.437andA=0.02.
σ 1 σ 2 ξ ± ζ λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.01 0.005 0.591620 0.253439 - 2.546035 ± 0.771621 i ± 14.71489 i 2.546035 ± 0.771621 i Unstable  0.02 0.010 0.596901 0.226272 - 8.552855 ± 0.7998829 i ± 26.47024 i 8.552855 ± 0.7998829 i Unstable  0.03 0.015 0.598669 0.215135 - 12.78173 ± 0.6731746 i ± 36.07696 i 12.78173 ± 0.6731746 i Unstable  0.04 0.020 0.599541 0.209176 - 16.19443 ± 0.4955597 i ± 44.06055 i 16.19443 ± 0.4955597 i Unstable  0.05 0.025 0.600069 0.205475 - 19.07191 ± 0.183883 i ± 50.88468 i 19.07191 ± 0.183883 i UnstableTable 5
Effects of oblateness of DP-Leonis on the locations and stability of out-of-plane points fora=0.95,e=0.2,σ1=0.03andσ2=0.015.
A ξ ± ζ λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.01 0.868127 0.143390 ±21.23742 ± 20.52800 ± 55.03901 i Unstable  0.02 0.861428 0.213009 - 7.258800 ± 0.8860236 i ± 20.43774 i 7.258800 ± 0.8860236 i Unstable  0.03 0.847143 0.268985 - 3.688747 ± 0.9403290 i ± 11.31166 i 3.688747 ± 0.9403290 i Unstable  0.04 0.823537 0.314020 - 2.172492 ± 0.9468649 i ± 7.487427 i 2.172492 ± 0.9468649 i Unstable  0.05 0.789157 0.347500 - 1.312145 ± 0.9376641 i ± 5.479935 i 1.312145 ± 0.9376641 i UnstableTable 6
Effects of eccentricity of DP-Leonis on the locations and stability of out-of-plane points fora=0.95,A=0.02,σ1=0.03andσ2=0.015.
e ξ ± ζ λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.1 0.863491 0.213009 - 7.219975 ± 0.8852772 i ± 20.31470 i 7.219975 ± 0.8852772 i Unstable  0.2 0.861428 0.213009 - 7.258800 ± 0.8860236 i ± 20.43774 i 7.258800 ± 0.8860236 i Unstable  0.3 0.857990 0.213009 - 7.311630 ± 0.8892137 i ± 20.62122 i 7.311630 ± 0.8892137 i Unstable  0.4 0.853177 0.213009 - 7.360223 ± 0.8975893 i ± 20.83215 i 7.360223 ± 0.8975893 i Unstable  0.5 0.846988 0.213009 - 7.378476 ± 0.9140821 i ± 21.02381 i 7.378476 ± 0.9140821 i UnstableTable 7
Effects of triaxiality of DP-Leonis on the locations and stability of out-of-plane points fora=0.95,e=0.2andA=0.02.
σ 1 σ 2 ξ ± ζ λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.01 0.005 0.851343 0.247062 - 1.61690 ± 0.8535960 i ± 8.628061 i 1.61690 ± 0.8535960 i Unstable  0.02 0.010 0.858907 0.223084 - 4.92729 ± 0.9.117207 i ± 15.14238 i 4.92729 ± 0.9.117207 i Unstable  0.03 0.015 0.861428 0.213009 - 7.25880 ± 0.8860236 i ± 20.43774 i 7.25880 ± 0.8860236 i Unstable  0.04 0.020 0.862689 0.207582 - 9.136514 ± 0.852032 i ± 24.82868 i 9.136514 ± 0.852032 i Unstable  0.05 0.025 0.863445 0.204200 - 10.71929 ± 0.8148062 i ± 28.57866 i 10.71929 ± 0.8148062 i Unstable


## 7. Discussion

The equations of motion of a third body under the influence of an oblate primary and a radiating-triaxial secondary moving in an elliptic Keplerian orbits are described by the governing equations (1).It is seen in Figures1–4 that the regions above and below the principal bodies are allowed to motion, while the regions between the principal bodies, the regions at the left of the primary body, the regions at the right of the secondary body, and the small regions nearer to the principal bodies are not allowed to motion. Also, a close study of Figures 5–9 reveals that all the regions are allowed to motion except the regions between the principal bodies and the small regions nearer to the principal bodies for each figure (i.e., Figures 5–9), while in Figure 10 all the regions are allowed to motion except the small regions nearer to the principal bodies due to the high value of eccentricity e. It is also seen that the topologies of the ZVC for the binary systems PSR 1903+0327 and DP-Leonis differ. This happens due to the effects of their mass ratio μ. We also observed that the topology of the ZVC may change depending on the energy of the third body.Equations (19) give the pair of points (ξ0,0,±ζ0) which correspond to the 1st family of the out-of-plane equilibrium points, L6,7. It can be seen from (19) that the point ξ0 is affected by the oblateness of the primary body, radiation and triaxiality of the secondary body, semimajor axis, and eccentricity of the orbits of the principal bodies. Meanwhile, the point ±ζ0 is independent of the semimajor axis and eccentricity of the orbits of the principal bodies; it is only affected by the oblateness of the primary body and radiation pressure as well as triaxiality of the secondary body. This is similar to [5, 21] in which the point ±ζ0 is independent of the semimajor axis and eccentricity of the orbits of principal bodies.It is seen in Figures12 and 14 that there exist values of triaxiality corresponding to the out-of-plane equilibrium points which coincides with the values corresponding to the collinear equilibrium points in particular, the inner collinear equilibrium points L2 for the binary systems PSR 1903+0327, and DP-Leonis, respectively.The numerical computations of the roots of characteristic Equation (20), given in Tables 2–7, reveal the existence of at least a positive real root and complex root with positive part for each particular set of values. This means that the motion of an infinitesimal mass near L6,7 is unbounded. Hence, the topologies of the ZVC and the nature of the computed roots of the characteristic equation (20) for the binary systems PSR 1903+0327 and DP-Leonis confirm the instability of the out-of-plane equilibrium points which agrees with [5, 14, 17, 20–22], while the characteristic roots given in Table 8 reveal that all the roots are purely imaginary. This means that the motion of an infinitesimal mass near L6,7 is stable. The stability of the out-of-plane equilibrium points resulted from the value of mass ratio μ, 0.32≤μ≤0.5, and the small value of eccentricity.Table 8
Effects of mass ratioμ on the locations and stability of out-of-plane points for a=0.95,e=0.2,q2=0.9999,A=0.08,σ1=0.03andσ2=0.015.
μ ξ 0 ± ζ 0 λ 1 ,  2 λ 3 ,  4 λ 5 ,  6 Remark 0.32 0.597001 0.481528 ± 0.5279522 i ± 1.326220 i ± 4.025670 i Stable  0.33 0.590700 0.484032 ± 0.5278752 i ± 1.329252 i ± 4.058162 i Stable  0.34 0.584181 0.486388 ± 0.5259979 i ± 1.333737 i ± 4.089879 i Stable  0.35 0.577463 0.488610 ± 0.5226098 i ± 1.339435 i ± 4.120928 i Stable  0.36 0.570563 0.490709 ± 0.5179634 i ± 1.346133 i ± 4.151399 i Stable  0.37 0.563495 0.492694 ± 0.5122757 i ±1.353646i ±4.181367i Stable  0.38 0.556273 0.494575 ± 0.5057316 i ± 1.361814 i ± 4.210892 i Stable  0.39 0.548909 0.496359 ± 0.4984878 i ± 1.370503 i ± 4.240027 i Stable  0.40 0.541413 0.498054 ± 0.4906756 i ± 1.379600 i ± 4.268812 i Stable  0.41 0.533794 0.499666 ± 0.4824047 i ± 1.389008 i ± 4.297285 i Stable  0.42 0.526062 0.501202 ± 0.4737663 i ± 1.398650 i ± 4.325475 i Stable  0.43 0.518225 0.502666 ± 0.4648359 i ± 1.408461 i ± 4.353408 i Stable  0.44 0.510290 0.504064 ± 0.4556760 i ± 1.418386 i ± 4.381104 i Stable  0.45 0.502262 0.505399 ± 0.4463377 i ± 1.428382 i ± 4.408583 i Stable  0.46 0.494149 0.506676 ± 0.4368632 i ± 1.438414 i ± 4.435859 i Stable  0.47 0.485956 0.507900 ± 0.4272867 i ± 1.448451 i ± 4.462947 i Stable  0.48 0.477687 0.509072 ± 0.4176362 i ± 1.458471 i ± 4.489856 i Stable  0.49 0.469348 0.510196 ± 0.4079342 i ± 1.468455 i ± 4.516598 i Stable  0.50 0.460942 0.511275 ± 0.3981986 i ± 1.478387 i ± 4.543180 i Stable



---

*Source: 3278946-2019-02-26.xml*