Next Article in Journal
Dirac, Majorana, Weyl in 4D
Previous Article in Journal
A Bayesian Approach to the Simulation Argument
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximate Analytical Periodic Solutions to the Restricted Three-Body Problem with Perturbation, Oblateness, Radiation and Varying Mass

1
School of Mathematical Science, Yangzhou University, Yangzhou 225002, China
2
Departament de Matemàtiques, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Catalonia, Spain
*
Author to whom correspondence should be addressed.
Universe 2020, 6(8), 110; https://doi.org/10.3390/universe6080110
Submission received: 17 June 2020 / Revised: 31 July 2020 / Accepted: 2 August 2020 / Published: 4 August 2020
(This article belongs to the Section Planetary Sciences)

Abstract

:
Against the background of a restricted three-body problem consisting of a supergiant eclipsing binary system, the two primaries are composed of a pair of bright oblate stars whose mass changes with time. The zero-velocity surface and curve of the problem are numerically studied to describe the third body’s motion area, and the corresponding five libration points are obtained. Moreover, the effect of small perturbations, Coriolis and centrifugal forces, radiative pressure, and the oblateness and mass parameters of the two primaries on the third body’s dynamic behavior is discussed through the bifurcation diagram. Furthermore, the second- and third-order approximate analytical periodic solutions around the collinear solution point L 3 in two-dimensional plane and three-dimensional spaces are presented by using the Lindstedt-Poincaré perturbation method.

1. Introduction

It is well known that Newton’s famous proposition “Principle” took a decisive step in our understanding of the universe. The origin of the universe is inseparable from the three-body dynamics, which has always existed since the universe’s birth. It can be said that there are no three-body dynamics; it is impossible to create a universe [1]. The three-body interactions in the early universe help people understand the black hole MACHO (abbreviation of MAssive Compact Halo Object) binary formation and the detectability of gravitational waves in black hole MACHO binaries coalesce [2,3]. In the theoretical framework of scalar-tensor gravitation, Zhou et al. [4] found that there is a collinear solution to the three-body problem in the presence of a scalar field, and studied the effect of the scalar field on this solution and the positions of Lagrange points through numerical examples. For the three-body problem in the spherical universe, their perturbation theory analysis showed that the rate of precession of two small and nearly circular solutions of identical particles is proportional to the square root of their initial distance and inversely proportional to the square of the radius of the universe [5]. In addition, the three-body problem is also widely used in the evolution of binary systems [6] and the dynamic analysis of binary asteroids [7,8], as well as other fields in the universe such as dark matter, galaxies, GW170817 (GW is short for gravitational wave), and Mukhanov-Sasaki Hamiltonian dynamics, and so forth (see References [9,10,11,12,13] for more information).
The mathematician Poincaré believed periodic solutions to be the unique avenue for solving the three-body problem. He stated that “what makes these solutions so precious to us, is that they are, so to say, the only opening through which we can try to enter a place.” There is too much literature in this area, mainly involving three aspects—qualitative analysis, approximate quantitative solution, and computer-aided numerical simulation (see References [14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32], and the references therein).
When considering multiple perturbation factors, Singh et al. did a great deal of work on the libration points and their stability of various restricted three-body models (see References [33,34,35,36,37,38,39,40,41]). These factors mainly include the oblateness of the primaries [35,36,37,38,40,41], the oblateness of the third body [37], and the small perturbation in centrifugal force and Coriolis force [35,36,37,38,39,40,41], the radiation pressure of the primaries [33,34,35,36,37,38,39,40,41], the changing mass of the primaries over time based on the Meshcherskii law [34,35,36], the changing mass of the third body governed by Jeans law [42], the triaxiality of the primaries [39]. Basic references on oblateness, Coriolis and centrifugal forces, and radiation pressure can be found in References [43,44,45,46,47].
In addition, some other researchers also discussed the periodic solutions and libration points of the restricted three-body problem (R3BP) under specific perturbations. For example, Abdullah et al. [48] investigated the locations and stability of libration points in the circular R3BP with Coriolis and centrifugal forces under the assumption that the two primaries and the infinitesimal body with variable masses. Under the assumption that the more massive primary has triaxiality, Abouelmagd et al. [49] obtained the approximate periodic solution near the collinear libration point of R3BP to the second-order. Since the star’s isotropic radiation or absorption is one of the factors that cause the mass of the celestial bodies to change with time, Bekov [50] considered the R3BP of variable masses, in which the motion of the two primaries was determined by the Gylden-Meshcherskii problem and found the collinear solutions, as well as the spatial solution and its domains of existence. When the more extensive primary is an oblate sphere, Zotos [51] numerically studied the dynamic behavior of the circular R3BP under the initial conditions of the solution are bounded, escaping and collisional, and found that the oblateness factor has a significant influence on the characteristics of the solutions. Based on the model of the binary system proposed in Reference [39], Gao and Wang [6] continued to study the analytical approximate periodic solutions around the collinear libration points. They examined the influence of the small perturbation in Coriolis and centrifugal forces, the triaxiality, and radiation pressure of the primaries on the third body’s dynamic behavior through the bifurcation diagram.
In this paper, the effects of small perturbations, Coriolis and centrifugal forces, radiation pressure, and the oblateness and variable mass of the primaries on the third body’s dynamic behavior will be discussed numerically through bifurcation diagram in the context of an R3BP consisting of a supergiant eclipsing binary system. Here we take into account the oblateness, Coriolis and centrifugal forces together because the oblateness of the primaries will cause both Coriolis force and centrifugal force to increase (see References [52,53]). By Lindstedt-Poincaré perturbation method, the approximate analytical solutions of the perturbed third body in two and three dimensions near the collinear libration points will be obtained. This paper is organized as follows—the motion equations in an autonomous system will be given in Section 2. The geometry structure of the zero velocity surfaces and curves will be illustrated in Section 3. Moreover, the bifurcation diagrams of the different perturbation parameters will be shown in Section 4. Furthermore, the periodic planar and spatial solutions near the collinear libration points will be obtained by means of Lindstedt-Poincaré in Section 5. In the last Section 6, conclusions will be stated.

2. Equations of Motion

In a barycentric coordinate O x y z , consider that two oblate variable-mass radiating primaries move at an angular velocity ω in a circular orbit around their center of mass. The plane of motion of the primaries (located on the x-axis) coincides with the x-y plane. The third body of infinitesimal mass that moves under the gravity of the primaries but does not affect their motion, then the dynamic governing equations (see Reference [35] or Reference [36] for more details) can be written as follows
d 2 x d t 2 2 ω d y d t α d ω d t α y = ω 2 x β μ 1 ( t ) q 1 x x 1 r 1 3 μ 2 ( t ) q 2 x x 2 r 2 3 3 μ 1 ( t ) q 1 A 1 ( t ) x x 1 2 r 1 5 3 μ 2 ( t ) q 2 A 2 ( t ) x x 2 2 r 2 5 , d 2 y d t 2 + 2 ω d x d t α + d ω d t α x = ω 2 y β μ 1 ( t ) q 1 y r 1 3 μ 2 ( t ) q 2 y r 2 3 3 μ 1 ( t ) q 1 A 1 ( t ) y 2 r 1 5 3 μ 2 ( t ) q 2 A 2 ( t ) y 2 r 2 5 , d 2 z d t 2 = μ 1 ( t ) q 1 z r 1 3 μ 2 ( t ) q 2 z r 2 3 3 μ 1 ( t ) q 1 A 1 ( t ) z 2 r 1 5 3 μ 2 ( t ) q 2 A 2 ( t ) z 2 r 2 5 ,
where r i 2 = x x i 2 + y 2 + z 2 ( i = 1 , 2 ) denotes the distance between the ith primary and the third body, and the two primaries are placed on x 1 , 0 , 0 and x 2 , 0 , 0 , respectively. For the distance r between the primaries, x 1 = μ 2 ( t ) r / μ ( t ) , x 2 = μ 1 ( t ) r / μ ( t ) , μ ( t ) = μ 1 ( t ) + μ 2 ( t ) , μ i ( t ) = f m i ( t ) , m i ( t ) represents the mass of the ith primary, f is the gravitational constant. In addition, the parameters q i and A i ( t ) stand for the radiation factor and the varying oblateness of the ith primary, respectively. The Coriolis force α and centrifugal force β satisfy α = 1 + ε , β = 1 + ε , here | ε | and | ε | are the corresponding small perturbations with values far less than 1. Furthermore, it can be followed from the reference [36] that the angular velocity holds the following expression
ω 2 = 3 μ ( t ) 2 κ r 3 ( t ) 2 3 + A 1 ( t ) + A 2 ( t ) r 2 ( t ) ,
where κ is an arbitrary dimensionless constant of a particular integral r μ ( t ) = κ C 2 in Gyldén-Meshcherskii problem (see References [54,55,56] for more information), it represents an arbitrary sum of the masses of the primaries, C is a constant of the area integral. Note that if the masses of the primaries are constant and the sum κ = 1 , then μ ( t ) = 1 by unitizing the gravitational constant. In addition, if the distance r between the primaries is 1, and without considering the oblateness of the second primary, that is, A 2 ( t ) = 0 . Then Equation (2) is simplified to ω 2 = 1 + 3 A 1 ( t ) / 2 , which is the same as the mean motion expression in Reference [47].
We now transform Equations (1) into the autonomous system u , v , w , τ by the Meshcherskii’s transformation (see Reference [35] or Reference [57])
x = u ( τ ) R ( t ) , y = v ( τ ) R ( t ) , z = w ( τ ) R ( t ) , d t / d τ = R 2 ( t ) , r i = d i R ( t ) ,
where d i is the distance between the third body and the ith primary, R ( t ) = ϕ t 2 + 2 φ t + γ , ϕ , φ and γ are constants, the unified Meshcherskii law
μ ( t ) = μ 0 / R ( t ) , μ 1 ( t ) = μ 10 / R ( t ) , μ 2 ( t ) = μ 20 / R ( t ) ,
where μ 0 , μ 10 , μ 20 are constants, and the particular solutions of Gyldén-Meshcherskii problem,
x 1 = u 1 R ( t ) , x 2 = u 2 R ( t ) , C = d 12 2 ω 0 , r = d 12 R ( t ) , ω = ω 0 / R 2 ( t ) ,
where d 12 is the distance between the two primaries at initial time. Applying transformation A i ( t ) = α i R 2 ( t ) , here α i = d E i 2 d P i 2 / 5 d 12 2 ( i = 1 , 2 ) refers to the oblateness of the ith primary at initial time, d E i and d P i are the equatorial radii and the polar radii of the ith primary, respectively. Similarly, following Singh and Leke [35], unify the mass and distance at the initial time and introduce the mass parameter μ m which satisfies μ 10 / μ 0 = 1 μ m , μ 20 / μ 0 = μ m , μ m ( 0 , 1 / 2 ] stands for the ratio of the smaller mass to the total mass of the two primaries, then μ 0 = κ , d 12 = 1 . Therefore, we obtain
ω 2 = 3 μ 0 2 κ d 12 3 R 4 ( t ) 2 3 + α 1 R 2 ( t ) + α 2 R 2 ( t ) d 12 2 R 2 ( t ) = 1 R 4 ( t ) 1 + 3 2 ( α 1 + α 2 ) .
It follows from Equation (5) that ω 0 2 = 1 + 3 ( α 1 + α 2 ) / 2 , α i 1 i = 1 , 2 . Then Equations (1) are reduced to the following form
u ¨ 2 α ω 0 v ˙ = Ω u , v ¨ + 2 α ω 0 u ˙ = Ω v , w ¨ = Ω w ,
with the potential function
Ω = ω 0 2 β + κ 1 2 u 2 + v 2 + ω 0 2 κ 1 w 2 2 + κ q 1 1 μ m d 1 + q 2 μ m d 2 + q 1 α 1 1 μ m 2 d 1 3 + q 2 α 2 μ m 2 d 2 3 ,
where d 1 = u + μ m 2 + v 2 + w 2 and d 2 = u + μ m 1 2 + v 2 + w 2 are the distances from the third body to the two primaries, respectively.

3. Zero Velocity Surfaces and Curves

The level surface of Equation (7) is an energy surface. The projection of this surface onto position space is named Hill’s region, and its boundary is the zero-velocity surface, on which the velocity of the third body is zero, so the third body can only move inside the zero-velocity surface. The intersection of this surface with the x-y plane is known as the zero-velocity curve. Note that Equation (7) admits a first integral
2 Ω u ˙ 2 + v ˙ 2 + w ˙ 2 = C ,
where C is the Jacobi constant and u ˙ 2 + v ˙ 2 + w ˙ 2 denotes the velocity of the relative motion. The geometric structure of 2 Ω = C in the two-dimensional plane and three-dimensional space is shown as a series of curves and surfaces for the different Jacobi constants C.
The zero-velocity surface (or curve) with perturbations helps us to analyze the third body’s dynamic behavior under the influence of a specific perturbation factor, which is different from that of classical R3BP. For example, when q 1 = q 2 = 1 , α 1 = α 2 = 0 and α = β = 1 , the zero-velocity surface is different from the classic R3BP due to the presence of κ . For the parameter values presented in Table 1, Figure 1 shows the four zero-velocity surfaces when κ takes 0.5, 1, 2 and 4, respectively. It is easy to find that the structure of the zero-velocity surface of Equation (9) is more open as the value of κ becomes smaller and smaller. Because the geometry of the zero-velocity surface (or curve) will vary with the value of the Jacobi constant C and κ , we can numerically analyze the zero-velocity surface (or curve) and then find the libration points of Equation (7). Without loss of generality, now we only study the surfaces and curves for the case κ = 1 .
For 2.82017 C 4.16993 , the forbidden region of the third body changes with the Jacobi constant C as shown in Figure 2 and Figure 3. With the decrease of C, five regular libration points L 1 , L 2 , , L 5 are discovered gradually and the forbidden region of motion of the third body is shrinking. More precisely, the third body can only fly in the two ellipsoidal areas around the two primaries when C = 4.16993 , and the two ellipsoidal motion areas are connected together through the “tangent point” (marked as librarian point L 1 ), and the third body can pass from one primary to the other through L 1 when C < 4.16993 , but it cannot fly out of these two ellipsoidal areas. When C = 3.59308 , the peanut-shaped flight area is connected to the outer space through a new “tangent point” L 2 , and when the value of C continues to decrease, the third body can only enter the outer space in this direction. When C drops to 3.57940, we will discover another new “tangent point” in the opposite direction, which is the libration point L 3 . At this time the third body can also fly from this direction to outer space when the value of C continues to decrease. For the Jacobi constant C = 2.82017 , the symmetrical libration points L 4 and L 5 appear at the same time. Now the third object can escape into outer space in any direction around the two primaries. Therefore, as a summary of Figure 2 and Figure 3, the corresponding five planar librations points can be identified as follows L 1 ( 0.018995 , 0 ) , L 2 ( 1.1989 , 0 ) , L 3 ( 1.19261 , 0 ) , L 4 ( 0.01424987 , 0.851953 ) , L 5 ( 0.01424987 , 0.851953 ) . Besides, there are several methods for determining the libration points. For example, the total number of collinear libration points for this problem can also be determined based on the theory of topological degree [58].

4. Bifurcation Analysis

We note that the three-body problem has a high degree of complexity due to its dynamic equations’ extremely strong nonlinearity. Any subtle changes in different parameters may substantially change the third body’s dynamic behavior. Thus, we will analyze the effect of different perturbation parameters on the third body’s dynamic behavior through the bifurcation diagram. After a large number of numerical simulations, we select u 0 = ( 0.08 , 0.001 , 0.2 , 0.002 , 1.6 , 0.001 ) as the initial value of the iteration. Combining the conditions in Section 2, we restrict the Coriolis force α and centrifugal force β to [1, 1.1] in the α - β -u frame, Figure 4a,b show the influence of the perturbation of α and β on the variables u-v, respectively, and reflect the corresponding chaotic characteristics. However, the joint perturbation of α and β shows that the amplitude of the third body in u direction no longer changes when 1.055 < α < 1.1 (see Figure 4c).
Consider the perturbation range of the oblateness coefficients of the two primaries is [ 0 , 0.1 ] . As shown in Figure 5a,b, as the oblateness coefficient α 1 of the first primary increases, the effect of this perturbation factor on the third body’s amplitude in u direction gradually increases, but the effect on the amplitude in v direction gradually decreases. When the oblateness coefficient α 2 of the second primary changes, the third body exhibits apparent chaotic behavior. However, the third body keeps moving periodically in u direction under the joint action of α 1 and α 2 (see Figure 5c). This indicates that the oblateness coefficients play a certain degree of “correction” to the two primaries’ dynamic behavior so that they can maintain regular motion as much as possible.
For the parameter κ , which represents the arbitrary sum of the primaries’ masses, we only show the bifurcation diagrams whose value is between 0 and 1. Figure 6a–c clearly show that no matter how the value of κ changes, it has a considerable influence on the third body’s dynamic behavior. The third body reflects the intricate, chaotic properties in the frames κ -u-v, κ -u-w, and κ -v-w. It is consistent with our common sense because the masses of the primaries are closely related to their gravitation, and this plays a significant role in determining the behavior of the third body in the sphere of gravity.
Next, we show the bifurcation diagrams of the mass parameter μ m in the interval ( 0 , 0.5 ] in Figure 7. It is easy to find that when μ m ’s value is small, the mass of one primary is enormous compared to the other. It has a significant impact on the third body’s amplitude in the three directions of u, v, and w. When μ m = 0.07 , the third body’s amplitude in u direction seems to be stable from its projection (see Figure 7a,b), but it is still apparent in the three-dimensional space. However, when the value of μ m gradually increases, that is, the difference between the masses of the two primaries is getting smaller and smaller, its influence on the third body’s amplitude in the v direction is gradually more considerable than that in u and w directions (see Figure 7b,c).
Similarly, we examine the effect of the radiation factors q 1 and q 2 on the third body’s dynamic behavior at the interval [ 0 , 1 ] . Figure 8a,b demonstrate the influence of the radiation factors of the first and second primaries, respectively, and Figure 8c reflects the impact of the two radiation factors acting simultaneously. As shown in Figure 8a,b, the radiation factor q 1 or q 2 has almost the same influence on the third body’s movement in u and v directions. However, the influence of q 1 on the third body’s amplitude is greater than that of q 2 , especially when q 1 < 0.57 or 0.67 < q 1 < 0.85 . In contrast, the impact of q 2 is much more moderate. Figure 8c shows that when q 1 0.0202 , regardless of how q 2 changes, the amplitude of the third body changes little in the u direction. In other words, the two radiations from the two primaries in the current state have almost little effect on the dynamic behavior of the third body in the u direction, which may be due to the weak radiation of the first primary. However, when q 1 [ 0.0202 , 0.3434 ] , the third body’s amplitude has a significant change, even if q 2 is very small. When the value of q 1 is greater than 0.3434, the third body exhibits approximately symmetrical motion in a larger area under the joint action of q 1 and q 2 .

5. Periodic Solutions near the Collinear Libration Points

5.1. Expansion of Two-Dimensional Dynamic Equations

The third body’s equations of motion in the u-v plane are
u ¨ 2 α ω 0 v ˙ = U u , v ¨ + 2 α ω 0 u ˙ = U v ,
where the potential function
U = ω 0 2 β + κ 1 2 u 2 + v 2 + κ q 1 1 μ m d a + q 2 μ m d b + q 1 α 1 1 μ m 2 d a 3 + q 2 α 2 μ m 2 d b 3 ,
where d a = u + μ m 2 + v 2 and d b = u + μ m 1 2 + v 2 are the distances from the third body to two primaries respectively. Then the right hand side of Equation (10) can be written as
U u = ω 0 2 β + κ 1 u q 1 κ 1 μ m u + μ m d a 3 q 2 κ μ m u + μ m 1 d b 3 3 q 1 κ α 1 1 μ m u + μ m 2 d a 5 3 q 2 κ α 2 μ m u + μ m 1 2 d b 5 , U v = ω 0 2 β + κ 1 v q 1 κ 1 μ m v d a 3 q 2 κ μ m v d b 3 3 q 1 κ α 1 1 μ m v 2 d a 5 3 q 2 κ α 2 μ m v 2 d b 5 .
Now we study the periodic solutions of the third body around the collinear libration points L i i = 1 , 2 , 3 , and u L i is the specific coordinate corresponding to L i . By letting u = u L i + ξ and v = η , then Equation (10) are reduced to the following form
ξ ¨ 2 α ω 0 η ˙ = U ξ , η ¨ + 2 α ω 0 ξ ˙ = U η .
Expanding the right hand side of Equation (13) to the second-order yields
ξ ¨ 2 α ω 0 η ˙ = K 1 ξ + K 2 ξ 2 + K 3 η 2 , η ¨ + 2 α ω 0 ξ ˙ = P 1 η + P 2 ξ η ,
and to the third-order in the same way, then we obtain
ξ ¨ 2 α ω 0 η ˙ = K 1 ξ + K 2 ξ 2 + K 3 η 2 + K 4 ξ 3 + K 5 ξ η 2 , η ¨ + 2 α ω 0 ξ ˙ = P 1 η + P 2 ξ η + P 3 η 3 + P 4 ξ 2 η ,
where corresponding coefficients can be found in Appendix A.

Two-Dimensional Periodic Solutions

Assuming that the periodic solutions for Equation (14) with respect to orbital parameter e e 1 have the following form
ξ = ξ 1 e + ξ 2 e 2 , η = η 1 e + η 2 e 2 .
Substituting Equation (16) into Equation (14) and comparing the terms with the same power of e, we have
ξ ¨ 1 2 α ω 0 η ˙ 1 K 1 ξ 1 = 0 , η ¨ 1 + 2 α ω 0 ξ ˙ 1 P 1 η 1 = 0 ,
and
ξ ¨ 2 2 α ω 0 η ˙ 2 K 3 η 1 2 K 1 ξ 2 K 2 ξ 1 2 = 0 , η ¨ 2 + 2 α ω 0 ξ ˙ 2 P 1 η 2 P 2 ξ 1 η 1 = 0 .
By using the Lindstedt-Poincaré perturbation method (see Chapter 5 in [59] for more details), the periodic solution of Equation (17) can be written as
ξ 1 = a 1 cos ϖ t + b 1 sin ϖ t , η 1 = a 2 cos ϖ t + b 2 sin ϖ t ,
with period T = 2 π / ϖ .
Substituting Equation (19) into Equation (17) and comparing the coefficients of sin ϖ t and cos ϖ t , then we get
K 1 + ϖ 2 a 1 + 2 α ω 0 ϖ b 2 = 0 , 2 α ω 0 ϖ a 1 + P 1 + ϖ 2 b 2 = 0 ,
and
K 1 + ϖ 2 b 1 2 α ω 0 ϖ a 2 = 0 , 2 α ω 0 ϖ b 1 + P 1 + ϖ 2 a 2 = 0 .
According to the homogeneous linear Equations (20) and (21), in order to get a nonzero solution, the corresponding determinant satisfies
ϖ 4 + K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + K 1 P 1 = 0 ,
then we set a 1 = 1 , b 1 = 0 , and substitute them into Equations (20) and (21) respectively, such that a 2 = 0 and b 2 = 2 α ω 0 ϖ 2 α ω 0 ϖ P 1 + ϖ 2 P 1 + ϖ 2 , so the periodic solution of Equation (17) is
ξ 1 = cos ϖ t , η 1 = b 2 sin ϖ t .
Similarly, we write the periodic solution of Equation (18) as
ξ 2 = a 3 + a 4 cos ϖ t + a 5 cos 2 ϖ t + b 3 sin ϖ t + b 4 sin 2 ϖ t , η 2 = a 6 cos ϖ t + a 7 cos 2 ϖ t + b 5 sin ϖ t + b 6 sin 2 ϖ t .
By using the same method as above, we obtain
ξ 2 = a 3 + a 5 cos 2 ϖ t , η 2 = b 6 sin 2 ϖ t .
Accordingly, the periodic solution of the Equation (14) with respect to e writes
ξ = ξ 0 + cos ϖ t e + a 3 + a 5 cos 2 ϖ t e 2 , η = b 2 sin ϖ t e + b 6 sin 2 ϖ t e 2 ,
where ξ 0 , 0 is the coordinate of the collinear libration point.
For Equation (15), suppose that the periodic solution in powers of parameter e takes the form
ξ = ξ 1 e + ξ 2 e 2 + ξ 3 e 3 , η = η 1 e + η 2 e 2 + η 3 e 3 .
Substituting Equation (27) into Equation (15) and comparing the terms with the same power of e. We obtain three sets of equations, of which the first two sets are identical to Equations (17) and (18) respectively, and the third set of equations has the form
ξ ¨ 3 2 α ω 0 η ˙ 3 = K 1 ξ 3 + 2 K 2 ξ 1 ξ 2 + 2 K 3 η 1 η 2 + K 4 ξ 1 3 + K 5 ξ 1 η 1 2 , η ¨ 3 + 2 α ω 0 ξ ˙ 3 = P 1 η 3 + P 2 ξ 1 η 2 + P 2 ξ 2 η 1 + P 3 η 1 3 + P 4 ξ 1 2 η 1 .
Let the periodic solution of Equations (28) be
ξ 3 = a 8 cos ϖ t + a 9 cos 2 ϖ t + a 10 cos 3 ϖ t + b 7 sin ϖ t + b 8 sin 2 ϖ t + b 9 sin 3 ϖ t , η 3 = a 11 cos ϖ t + a 12 cos 2 ϖ t + a 13 cos 3 ϖ t + b 10 sin ϖ t + b 11 sin 2 ϖ t + b 12 sin 3 ϖ t .
Substituting Equations (23), (25) and (29) into Equation (28), and comparing the coefficients of sin j ϖ t and cos j ϖ t ( j = 1 , 2 , 3 ) respectively, then we have
ξ 3 = a 10 cos 3 ϖ t , η 3 = b 12 sin 3 ϖ t ,
where the expressions of a 10 and b 12 can be found from Appendix B.
Therefore, the periodic solution of Equation (15) in powers of parameter e becomes
ξ = ξ 0 + cos ϖ t e + a 3 + a 5 cos 2 ϖ t e 2 + a 10 cos 3 ϖ t e 3 , η = b 2 sin ϖ t e + b 6 sin 2 ϖ t e 2 + b 12 sin 3 ϖ t e 3 ,
where the expression of coefficients, please refer to Appendix B again.
Note that different values of κ not only determine the positions of the libration points, but also affect the shape of the second- and third-order periodic solutions through Equation (22). This is because the coefficients of Equation (22) are related to κ (see Appendix A for details). In order to better illustrate the approximate periodic solution obtained above, we show the phase portraits of the above approximate analytical periodic solutions near the libration point L 3 under the four cases of κ = 0.5 , 1, 2, and 4 respectively (see Figure 9a,d). When the value of κ is small, it is easy to find that the difference between the second- and the third-order approximate periodic solution is insignificant (see Figure 9a,b). Especially for the periodic solutions in Figure 9a, the second- and third-order approximate analytical periodic solutions on the two-dimensional plane around L 3 agree very well, and they are almost indistinguishable visually. In this state, the second-order approximate periodic solution can be regarded as the alternative to the third-order case. If the values of other parameters remain unchanged, we can easily find from Figure 9c,d that the difference between the second- and third-order approximate periodic solutions become more and more significant as κ increases. Therefore, we can conclude according to Figure 9 that the smaller the value of κ , the better the degree of agreement between the second- and third-order approximate periodic solutions, and vice versa.

5.2. Expansion of Three-Dimensional Dynamic Equations

Applying the translations u = u L i + ξ , v = η and w = ζ to Equation (7), then we have
ξ ¨ 2 α ω 0 η ˙ = Ω ξ , η ¨ + 2 α ω 0 ξ ˙ = Ω η , ζ ¨ = Ω ζ .
Expanding the right hand side of Equation (32) to the second- and to the third-order respectively, then we obtain
ξ ¨ 2 α ω 0 η ˙ = K 1 ξ + K 2 ξ 2 + K 3 η 2 + K 6 ζ 2 , η ¨ + 2 α ω 0 ξ ˙ = P 1 η + P 2 ξ η , ζ ¨ = Q 1 ζ + Q 2 ξ ζ ,
and
ξ ¨ 2 α ω 0 η ˙ = K 1 ξ + K 2 ξ 2 + K 3 η 2 + K 6 ζ 2 + K 7 ξ 3 + K 8 ξ η 2 + K 9 ξ ζ 2 , η ¨ + 2 α ω 0 ξ ˙ = P 1 η + P 2 ξ η + P 5 η 3 + P 6 ξ 2 η + P 7 η ζ 2 , ζ ¨ = Q 1 ζ + Q 2 ξ ζ + Q 3 ζ 3 + Q 4 ξ 2 ζ + Q 5 η 2 ζ ,
where the corresponding coefficients are displayed in the Appendix A.

Three-Dimensional Periodic Solutions

By means of successive approximations, we assume that the periodic solution of Equation (33) admits the following form
ξ = ξ 1 e + ξ 2 e 2 , η = η 1 e + η 2 e 2 , ζ = ζ 1 e + ζ 2 e 2 .
Substituting Equation (35) into Equation (33) and comparing these terms with the same power e, then we obtain
ξ ¨ 1 2 α ω 0 η ˙ 1 K 1 ξ 1 = 0 , η ¨ 1 + 2 α ω 0 ξ ˙ 1 P 1 η 1 = 0 , ζ ¨ 1 Q 1 ζ 1 = 0 ,
and
ξ ¨ 2 2 α ω 0 η ˙ 2 K 1 ξ 2 K 2 ξ 1 2 K 3 η 1 2 K 6 ζ 1 2 = 0 , η ¨ 2 + 2 α ω 0 ξ ˙ 2 P 1 η 2 P 2 ξ 1 η 1 = 0 , ζ ¨ 2 Q 1 ζ 2 Q 2 ξ 1 ζ 1 = 0 .
For the third equation of Equation (36), the corresponding periodic solution takes the form
ζ 1 = a cos ϖ t + b sin ϖ t ,
and substitute Equation (38) into the third equation of Equation (36) yields
Q 1 + ϖ 2 ζ 1 = 0 .
Note that if there is a non-zero solution to Equation (39), then Q 1 = ϖ 2 < 0 . So we choose a = 0 and b = 1 , and one of the solutions is ζ 1 = sin ϖ t .
Let the periodic solution of Equation (36) admits the following form
ξ 1 = a 14 cos ϖ t + b 13 sin ϖ t , η 1 = a 15 cos ϖ t + b 14 sin ϖ t .
Substituting Equation (40) into Equation (36) and comparing the coefficients of sin ϖ t and cos ϖ t , then we obtain
K 1 + ϖ 2 a 14 + 2 α ω 0 ϖ b 14 = 0 , 2 α ω 0 ϖ a 14 + P 1 + ϖ 2 b 14 = 0 ,
and
K 1 + ϖ 2 b 13 2 α ω 0 ϖ a 15 = 0 , 2 α ω 0 ϖ b 13 + P 1 + ϖ 2 a 15 = 0 .
Based on the results of Equations (20) and (21), then the periodic solution of Equation (36) writes
ξ 1 = cos ϖ t , η 1 = b 14 sin ϖ t , ζ 1 = sin ϖ t .
Similarly, let the periodic solution of Equation (37) has the following form
ξ 2 = a 16 + a 17 cos ϖ t + a 18 cos 2 ϖ t + b 15 sin ϖ t + b 16 sin 2 ϖ t , η 2 = a 19 cos ϖ t + a 20 cos 2 ϖ t + b 17 sin ϖ t + b 18 sin 2 ϖ t , ξ 2 = a 21 cos ϖ t + a 22 cos 2 ϖ t + b 19 sin ϖ t + b 20 sin 2 ϖ t .
By using the same method as mentioned above, then we get
ξ 2 = a 16 + a 18 cos 2 ϖ t , η 2 = b 18 sin 2 ϖ t , ζ 2 = b 20 sin 2 ϖ t .
Therefore, the periodic solution of Equation (33) with respect to e becomes
ξ = ξ 0 + cos ϖ t e + a 16 + a 18 cos 2 ϖ t e 2 , η = b 14 sin ϖ t e + b 18 sin 2 ϖ t e 2 , ζ = sin ϖ t e + b 20 sin 2 ϖ t e 2 ,
where the expressions of these coefficients can be obtained from Appendix B.
As for the third-order approximate periodic solution of Equation (34), suppose that the periodic solution in the form of e can be written as
ξ = ξ 1 e + ξ 2 e 2 + ξ 3 e 3 , η = η 1 e + η 2 e 2 + η 3 e 3 , ζ = ζ 1 e + ζ 2 e 2 + ζ 3 e 3 .
Substituting Equation (47) into Equation (34), and comparing the terms containing the same power of e, then we obtain three set of equations, the equations corresponding to e and e 2 are the same as Equations (36) and (37), respectively, and the equations corresponding to e 3 are as follows
ξ ¨ 3 2 α ω 0 η ˙ 3 = K 1 ξ 3 + 2 K 2 ξ 1 ξ 2 + 2 K 3 η 1 η 2 + 2 K 6 ζ 1 ζ 2 + K 7 ξ 1 3 + K 8 ξ 1 η 1 2 + K 9 ξ 1 ζ 1 2 , η ¨ 3 + 2 α ω 0 ξ ˙ 3 = P 1 η 3 + P 2 ξ 1 η 2 + P 2 ξ 2 η 1 + P 5 η 1 3 + P 6 ξ 1 2 η 1 + P 7 η 1 ζ 1 2 , ζ ¨ 3 = Q 1 ζ 3 + Q 2 ξ 1 ζ 2 + Q 2 ξ 2 ζ 1 + Q 3 ζ 1 3 + Q 4 ξ 1 2 ζ 1 + Q 5 η 1 2 ζ 1 .
Writing the periodic solution of Equations (48) as
ξ 3 = a 23 cos ϖ t + a 24 cos 2 ϖ t + a 25 cos 3 ϖ t + b 21 sin ϖ t + b 22 sin 2 ϖ t + b 23 sin 3 ϖ t , η 3 = a 26 cos ϖ t + a 27 cos 2 ϖ t + a 28 cos 3 ϖ t + b 24 sin ϖ t + b 25 sin 2 ϖ t + b 26 sin 3 ϖ t , ζ 3 = a 29 cos ϖ t + a 30 cos 2 ϖ t + a 31 cos 3 ϖ t + b 27 sin ϖ t + b 28 sin 2 ϖ t + b 29 sin 3 ϖ t .
Substituting Equations (43), (45) and (49) into Equation (48), and comparing the coefficients of sin j ϖ t and cos j ϖ t ( j = 1 , 2 , 3 ) , then we obtain
ξ 3 = a 25 cos 3 ϖ t , η 3 = b 26 sin 3 ϖ t , ζ 3 = b 29 sin 3 ϖ t .
Therefore, the third-order approximate analytical periodic solution of Equation (7) in three-dimensional space can be written as
ξ = ξ 0 + cos ϖ t e + a 16 + a 18 cos 2 ϖ t e 2 + a 25 cos 3 ϖ t e 3 , η = b 14 sin ϖ t e + b 18 sin 2 ϖ t e 2 + b 26 sin 3 ϖ t e 3 , ζ = sin ϖ t e + b 20 sin 2 ϖ t e 2 + b 29 sin 3 ϖ t e 3 .
Note that one can find the specific expressions of the above coefficients from Appendix B.
Based on the results of the above analysis, we show the three-dimensional periodic solution near the collinear libration point L 3 for four different κ values of 0.5, 1, 2, and 4 in Figure 10. When κ = 0.5 and κ = 1 , the second- and third-order three-dimensional periodic solutions agree very well, as shown in Figure 10a,b. When κ = 2, although Figure 10c shows that the second- and third-order periodic solutions are somewhat different, the degree of agreement is also good. However, Figure 10d (corresponding to κ = 4) indicates that the second- and third-order periodic solutions differ significantly. Therefore, it is easy to find that the smaller the value of κ , the better the approximation effect.

6. Conclusions

For a restricted three-body problem consisting of a supergiant eclipsing binary system, we investigated the geometric structure of the corresponding zero-velocity surfaces and curves of the coplanar libration points and discussed the third body’s dynamic behavior under the influence of eight perturbation parameters through bifurcation diagrams. These perturbations include: Coriolis force α , centrifugal force β , oblateness coefficients α 1 and α 2 , random sums κ of the masses of the two primaries, mass ratio parameter μ m , and the primaries’ radiation factors q 1 and q 2 . We found that any perturbation of Coriolis force and centrifugal force has a greater impact on the third body’s dynamic behavior, but their combined effect has a smaller effect on the third body. Similarly, the oblateness coefficient α 2 has a more significant influence on the third body’s dynamic behavior than that of α 1 . However, the degree of their combined effect on the third body’s dynamic behavior is minimal. Additionally, the parameter κ has a more extensive influence on the dynamic behavior of the third body, any tiny changes can cause apparent chaos phenomena. Furthermore, the smaller the mass ratio parameter, the more considerable the dynamic behavior of the third body, and vice versa. Besides, the influence of the radiation factor q 1 on the three-body dynamics is more significant than that of q 2 . The combined effect of q 1 and q 2 has a notable effect on the dynamics of the third body when the radiation factor q 1 > 0.3434 .
In addition, we also applied the Lindstedt-Poincaré perturbation method to present the second- and third-order approximate analytical periodic solutions near the colinear libration points in the plane and space respectively, and illustrated the corresponding phase portraits with respect to κ = 0.5 , 1, 2, and 4 around the libration point L 3 . We found that the smaller the value of κ , the better the coincidence of the corresponding second-order and third-order periodic solutions. Besides, for the same κ value, the coincidence between the three-dimensional periodic solutions of different orders is better than that of the two-dimensional case.

Author Contributions

Formal analysis, Y.W.; software, F.G.; writing-original draft, Y.W.; writing—review and editing, F.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) though grant No. 11672259, the China Scholarship Council through grant No. 201908320086.

Acknowledgments

We are very grateful to the anonymous reviewers whose comments and suggestions helped improve and clarify this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

K 1 = ω 0 2 β + κ 1 + 2 q 1 κ 1 μ m 1 d 01 3 + 3 α 1 d 01 5 + 2 q 2 κ μ m 1 d 02 3 + 3 α 2 d 02 5 ,
K 2 = 3 q 1 κ 1 μ m 1 d 01 4 + 5 α 1 d 01 6 V 1 3 q 2 κ μ m 1 d 02 4 + 5 α 2 d 02 6 V 2 ,
K 3 = 3 2 q 1 κ 1 μ m 1 d 01 4 + 5 α 1 2 d 01 6 V 1 + 3 2 q 2 κ μ m 1 d 02 4 + 5 α 2 2 d 02 6 V 2 ,
K 4 = 2 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 d 01 7 + 2 q 2 κ μ m 2 d 02 5 + 15 α 2 d 02 7 ,
K 5 = 3 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 2 d 01 7 3 q 2 κ μ m 2 d 02 5 + 15 α 2 2 d 02 7 ,
K 6 = 3 2 q 1 κ 1 μ m 1 d 01 4 + 5 α 1 2 d 01 6 V 1 + 3 2 q 2 κ μ m 1 d 02 4 + 5 α 2 2 d 02 6 V 2 ,
K 7 = 2 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 d 01 7 + 2 q 2 κ μ m 2 d 02 5 + 15 α 2 d 02 7 ,
K 8 = 3 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 2 d 01 7 3 q 2 κ μ m 2 d 02 5 + 15 α 2 2 d 02 7 ,
K 9 = 3 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 2 d 01 7 3 q 2 κ μ m 2 d 02 5 + 15 α 2 2 d 02 7 ,
P 1 = ω 0 2 β + κ 1 q 1 κ 1 μ m 1 d 01 3 + 3 α 1 2 d 01 5 q 2 κ μ m 1 d 02 3 + 3 α 2 2 d 02 5 ,
P 2 = 3 q 1 κ 1 μ m 1 d 01 4 + 5 α 1 2 d 01 6 V 1 + 3 q 2 κ μ m 1 d 02 4 + 5 α 2 2 d 02 6 V 2 ,
P 3 = 3 2 q 1 κ 1 μ m 1 d 01 5 + 5 α 1 2 d 01 7 + 3 2 q 2 κ μ m 1 d 02 5 + 5 α 2 2 d 02 7 ,
P 4 = 3 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 2 d 01 7 3 q 2 κ μ m 2 d 02 5 + 15 α 2 2 d 02 7 ,
P 5 = 3 2 q 1 κ 1 μ m 1 d 01 5 + 5 α 1 2 d 01 7 + 3 2 q 2 κ μ m 1 d 02 5 + 5 α 2 2 d 02 7 ,
P 6 = 3 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 2 d 01 7 3 q 2 κ μ m 2 d 02 5 + 15 α 2 2 d 02 7 ,
P 7 = 3 2 q 1 κ 1 μ m 1 d 01 5 + 5 α 1 2 d 01 7 + 3 2 q 2 κ μ m 1 d 02 5 + 5 α 2 2 d 02 7 ,
V 1 = s g n u L i + μ m ,
V 2 = s g n u L i + μ m 1 ,
d 01 = u L i + μ m ,
d 02 = u L i + μ m 1 ,
Q 1 = ω 0 2 κ 1 q 1 κ 1 μ m 1 d 01 3 + 3 α 1 2 d 01 5 q 2 κ μ m 1 d 02 3 + 3 α 2 2 d 02 5 ,
Q 2 = 3 q 1 κ 1 μ m 1 d 01 4 + 5 α 1 2 d 01 6 V 1 + 3 q 2 κ μ m 1 d 02 4 + 5 α 2 2 d 02 6 V 2 ,
Q 3 = 3 2 q 1 κ 1 μ m 1 d 01 5 + 5 α 1 2 d 01 7 + 3 2 q 2 κ μ m 1 d 02 5 + 5 α 2 2 d 02 7 ,
Q 4 = 3 q 1 κ 1 μ m 2 d 01 5 + 15 α 1 2 d 01 7 3 q 2 κ μ m 2 d 02 5 + 15 α 2 2 d 02 7 ,
Q 5 = 3 2 q 1 κ 1 μ m 1 d 01 5 + 5 α 1 2 d 01 7 + 3 2 q 2 κ μ m 1 d 02 5 + 5 α 2 2 d 02 7 .

Appendix B

a 3 = K 2 + K 3 b 2 2 2 K 1 ,
a 5 = P 1 + 4 ϖ 2 K 3 b 2 2 K 2 + 4 α ω 0 ϖ P 2 b 2 32 ϖ 4 + 8 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + 2 K 1 P 1 ,
a 10 = P 1 + 9 ϖ 2 K 2 a 5 + K 3 b 2 b 6 1 4 K 4 + 1 4 K 5 b 2 2 + 3 α ω 0 ϖ P 2 b 6 + P 2 a 5 b 2 1 2 P 3 b 2 3 + 1 2 P 4 b 2 81 ϖ 4 + 9 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + K 1 P 1 ,
a 16 = K 2 + K 3 b 14 2 + K 6 2 K 1 ,
a 18 = P 1 + 4 ϖ 2 K 2 + K 3 b 14 2 + K 6 + 4 α ω 0 ϖ P 2 b 14 32 ϖ 4 + 8 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + 2 K 1 P 1 ,
a 25 = 1 81 ϖ 4 + 9 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + K 1 P 1 P 1 + 9 ϖ 2 K 2 a 18 + K 3 b 14 b 18 + K 6 b 20 1 4 K 7
+ 1 4 K 8 b 14 2 + 1 4 K 9 6 α ω 0 ϖ 1 2 P 2 b 18 1 2 P 2 a 18 b 14 + 1 4 P 5 b 14 3 1 4 P 6 b 14 + 1 4 P 7 b 14 ,
b 2 = 2 α ω 0 ϖ P 1 + ϖ 2 ,
b 6 = P 2 b 2 K 1 + 4 ϖ 2 4 α ω 0 ϖ K 3 b 2 2 K 2 32 ϖ 4 + 8 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + 2 K 1 P 1 ,
b 12 = K 1 + 9 ϖ 2 1 2 P 2 b 6 1 2 P 2 a 5 b 2 + 1 4 P 3 b 2 3 1 4 P 4 b 2 6 α ω 0 ϖ K 2 a 5 + K 3 b 2 b 6 1 4 K 4 + 1 4 K 5 b 2 2 81 ϖ 4 + 9 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + K 1 P 1 ,
b 14 = 2 α ω 0 ϖ P 1 + ϖ 2 ,
b 18 = P 2 b 14 K 1 + 4 ϖ 2 4 α ω 0 ϖ K 2 + K 3 b 14 2 + K 6 32 ϖ 4 + 8 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + 2 K 1 P 1 ,
b 20 = Q 2 6 ϖ 2 ,
b 26 = 1 81 ϖ 4 + 9 K 1 + P 1 4 α 2 ω 0 2 ϖ 2 + K 1 P 1 K 1 + 9 ϖ 2 1 2 P 2 b 18 1 2 P 2 a 18 b 14 + 1 4 P 5 b 14 3
1 4 P 6 b 14 + 1 4 P 7 b 14 6 α ω 0 ϖ K 2 a 18 + K 3 b 14 b 18 + K 6 b 20 1 4 K 7 + 1 4 K 8 b 14 2 + 1 4 K 9 ,
b 29 = 1 2 Q 2 b 20 1 2 Q 2 a 18 + 1 4 Q 3 1 4 Q 4 + 1 4 Q 5 b 14 2 8 ϖ 2 .

References

  1. Marchal, C. The Three-Body Problem; Elsevier Science: New York, NY, USA, 1990. [Google Scholar]
  2. Ioka, K.; Chiba, T.; Tanaka, T.; Nakamura, T. Black hole binary formation in the expanding universe: Three body problem approximation. Phys. Rev. D 1998, 58, 063003. [Google Scholar] [CrossRef] [Green Version]
  3. Nakamura, T.; Sasaki, M.; Tanaka, T.; Thorne, K.S. Gravitational waves from coalescing black hole MACHO binaries. Astrophys. J. Lett. 1997, 487, L139. [Google Scholar] [CrossRef] [Green Version]
  4. Zhou, T.Y.; Cao, W.G.; Xie, Y. Collinear solution to the three-body problem under a scalar-tensor gravity. Phys. Rev. D 2016, 93, 064065. [Google Scholar] [CrossRef]
  5. Lindner, J.F.; Roseberry, M.I.; Shai, D.E.; Harmon, N.J.; Olaksen, K.D. Precession and chaos in the classical two-body problem in a spherical universe. Int. J. Bifurc. Chaos 2008, 18, 455–464. [Google Scholar] [CrossRef]
  6. Gao, F.B.; Wang, R.F. Bifurcation analysis and periodic solutions of the HD 191408 system with triaxial and radiative perturbations. Universe 2020, 6, 35. [Google Scholar] [CrossRef] [Green Version]
  7. Beletsky, V.V. Generalized restricted circular three-body problem as a model for dynamics of binary asteroids. Cosm. Res. 2007, 45, 408–416. [Google Scholar] [CrossRef]
  8. Hou, X.Y.; Xin, X.S.; Feng, J.L. Forced motions around triangular libration points by solar radiation pressure in a binary asteroid system. Astrodynamics 2020, 4, 17–30. [Google Scholar] [CrossRef]
  9. Braaten, E.; Kang, D.; Laha, R. Production of dark-matter bound states in the early universe by three-body recombination. J. High Energy Phys. 2018, 2018, 84. [Google Scholar] [CrossRef] [Green Version]
  10. Trani, A.A.; Fujii, M.S.; Spera, M. The Keplerian three-body encounter. I. Insights on the origin of the S-stars and the G-objects in the galactic center. Astrophys. J. 2019, 875, 42. [Google Scholar] [CrossRef]
  11. Muzzio, J.C.; Wachlin, F.C.; Carpintero, D.D. Regular and chaotic motion in a restricted three-body problem of astrophysical interest. Int. Astron. Union Colloq. 2000, 174, 281–285. [Google Scholar] [CrossRef] [Green Version]
  12. Logoteta, D.; Bombaci, I. Constraints on microscopic and phenomenological equations of state of dense matter from GW170817. Universe 2019, 5, 204. [Google Scholar] [CrossRef] [Green Version]
  13. Fahn, M.J.; Giesel, K.; Kobler, M. Dynamical properties of the Mukhanov-Sasaki Hamiltonian in the context of adiabatic vacua and the Lewis-Riesenfeld invariant. Universe 2019, 5, 170. [Google Scholar] [CrossRef] [Green Version]
  14. Gao, F.B.; Zhang, W. A study on periodic solutions for the circular restricted three-body problem. Astron. J. 2014, 148, 116. [Google Scholar] [CrossRef] [Green Version]
  15. Alshaery, A.A.; Abouelmagd, E.I. Analysis of the spatial quantized three-body problem. Results Phys. 2020, 17, 103067. [Google Scholar] [CrossRef]
  16. Zotos, E.E.; Chen, W.; Abouelmagd, E.I.; Han, H.T. Basins of convergence of equilibrium points in the restricted three-body problem with modified gravitational potential. Chaos Solitons Fractals 2020, 134, 109704. [Google Scholar] [CrossRef]
  17. Alzahrani, F.; Abouelmagd, E.I.; Guirao, J.L.G.; Hobiny, A. On the libration collinear points in the restricted three-body problem. Open Phys. 2017, 15, 58–67. [Google Scholar] [CrossRef]
  18. Abouelmagd, E.I.; Guirao, J.L.G.; Llibre, J. Periodic orbits for the perturbed planar circular restricted 3-body problem. Discret. Contin. Dyn. Syst.-Ser. B 2019, 24, 1007–1020. [Google Scholar] [CrossRef] [Green Version]
  19. Selim, H.H.; Guirao, J.L.G.; Abouelmagd, E.I. Libration points in the restricted three-body problem: Euler angles, existence and stability. Discret. Contin. Dyn. Syst.-Ser. B 2019, 12, 703–710. [Google Scholar] [CrossRef] [Green Version]
  20. Abouelmagd, E.I.; Mostafa, A.; Guirao, J.L.G. A first order automated lie transform. Int. J. Bifurc. Chaos 2015, 25, 1540026. [Google Scholar] [CrossRef] [Green Version]
  21. Abouelmagd, E.I.; Asiri, H.M.; Sharaf, M.A. The effect of oblateness in the perturbed restricted three-body problem. Meccanica 2013, 48, 2479–2490. [Google Scholar] [CrossRef]
  22. Pathak, N.; Thomas, V.O.; Abouelmagd, E.I. The perturbed photogravitational restricted three-body problem: Analysis of resonant periodic orbits. Discret. Contin. Dyn. Syst.-Ser. B 2019, 12, 849–875. [Google Scholar] [CrossRef] [Green Version]
  23. Niederman, L.; Pousse, A.; Robutel, P. On the co-orbital motion in the three-body problem: Existence of quasi-periodic horseshoe-shaped orbits. Commun. Math. Phys. 2020, 377, 551–612. [Google Scholar] [CrossRef] [Green Version]
  24. Koon, W.S.; Lo, M.W.; Marsden, J.E.; Ross, S.D. Dynamical Systems, the Three-Body Problem and Space Mission Design; Marsden Books: Wellington, New Zealand, 2011; ISBN 978-0-615-24095-4. [Google Scholar]
  25. Qian, Y.J.; Yang, X.D.; Yang, L.Y.; Zhang, W. Approximate analytical methodology for the restricted three-body and four-body models based on polynomial series. Int. J. Aerosp. Eng. 2016, 2016, 9747289. [Google Scholar] [CrossRef] [Green Version]
  26. Pathak, N.; Abouelmagd, E.I.; Thomas, V.O. On higher order resonant periodic orbits in the photo-gravitational planar restricted three-body problem with oblateness. J. Astronaut. Sci. 2019, 66, 475–505. [Google Scholar] [CrossRef]
  27. Hénon, M. Generating Families in the Restricted Three-Body Problem; Springer: New York, NY, USA, 1997. [Google Scholar]
  28. Chenciner, A.; Montgomery, R. A remarkable periodic solution of the three-body problem in the case of equal masses. Ann. Math. 2000, 152, 881–901. [Google Scholar] [CrossRef] [Green Version]
  29. Šuvakov, M.; Dmitrašinović, V. Three classes of Newtonian three-body planar periodic solutions. Phys. Rev. Lett. 2013, 110, 114301. [Google Scholar] [CrossRef] [Green Version]
  30. Janković, M.R.; Dmitrašinović, V.; Šuvakov, M. A guide to hunting periodic three-body orbits with non-vanishing angular momentum. Comput. Phys. Commun. 2020, 250, 107052. [Google Scholar] [CrossRef]
  31. Li, X.M.; Liao, S.J. More than six hundred new families of Newtonian periodic planar collisionless three-body orbits. Sci. China Phys. Mech. Astron. 2017, 60, 129511. [Google Scholar] [CrossRef] [Green Version]
  32. Kalantonis, V.S. Numerical investigation for periodic solutions in the Hill three-body problem. Universe 2020, 6, 72. [Google Scholar] [CrossRef]
  33. Singh, J. Photogravitational restricted three-body problem with variable mass. Indian Indian J. Pure Appl. Math. 2003, 34, 335–341. [Google Scholar]
  34. Singh, J.; Leke, O. Stability of the photogravitational restricted three-body problem with variable masses. Astrophys. Space Sci. 2010, 326, 305–314. [Google Scholar] [CrossRef]
  35. Singh, J.; Leke, O. Equilibrium points and stability in the restricted three-body problem with oblateness and variable masses. Astrophys. Space Sci. 2012, 340, 27–41. [Google Scholar] [CrossRef]
  36. Singh, J.; Leke, O. Effect of oblateness, perturbations, radiation and varying masses on the stability of equilibrium points in the restricted three-body problem. Astrophys. Space Sci. 2013, 344, 51–61. [Google Scholar] [CrossRef]
  37. Singh, J.; Haruna, S. Equilibrium points and stability under effect of radiation and perturbing forces in the restricted problem of three oblate bodies. Astrophys. Space Sci. 2014, 349, 107–116. [Google Scholar] [CrossRef]
  38. Singh, J.; Perdiou, A.E.; Gyegwe, J.M.; Kalantonis, V.S. Periodic orbits around the collinear equilibrium points for binary Sirius, Procyon, Luhman 16, α-Centuari and Luyten 726-8 systems: The spatial case. J. Phys. Commun. 2017, 1, 025008. [Google Scholar] [CrossRef]
  39. Singh, J.; Perdiou, A.E.; Gyegwe, J.M.; Perdios, E.A. Periodic solutions around the collinear equilibrium points in the perturbed restricted three-body problem with triaxial and radiating primaries for binary HD 191408, Kruger 60 and HD 155876 systems. Appl. Math. Comput. 2018, 325, 358–374. [Google Scholar] [CrossRef]
  40. Abdulraheem, A.R.; Singh, J. Combined effects of perturbations, radiation and oblateness on the periodic solutions in the restricted three-body problem. Astrophys. Space Sci. 2008, 317, 9–13. [Google Scholar] [CrossRef] [Green Version]
  41. Abdulraheem, A.R.; Singh, J. Combined effects of perturbations, radiation, and oblateness on the stability of equilibrium points in the restricted three-body problem. Astron. J. 2006, 131, 1880. [Google Scholar] [CrossRef]
  42. Singh, J.; Ishwar, B. Effect of perturbations on the location of equilibrium points in the restricted problem of three bodies with variable mass. Celestial Mech. 1984, 32, 297–305. [Google Scholar] [CrossRef]
  43. Bhatnagar, K.B.; Chawla, J.M. A study of the Lagrangian points in the photogravitational restricted three-body problem. Indian J. Pure Appl. Math. 1979, 10, 1443–1451. [Google Scholar]
  44. Szebehely, V. Stability of the points of equilibrium in the restricted problem. Astron. J. 1967, 72, 7–9. [Google Scholar] [CrossRef]
  45. Subba Rao, P.V.; Sharma, R.K. Perturbations of the critical mass in the restricted three-body problem. Astrophys. Space Sci. 1976, 42, L17–L18. [Google Scholar]
  46. Sharma, R.K. Perturbations of Lagrangian points in the restricted three-body problem. Indian J. Pure Appl. Math. 1975, 6, 1099–1102. [Google Scholar]
  47. Sharma, R.K.; Subba Rao, P.V. Stationary solutions and their characteristic exponents in the restricted three-body problem when the more massive primary is an oblate spheroid. Celestial Mech. 1976, 13, 137–149. [Google Scholar] [CrossRef]
  48. Abdullah, A.A.; Alhussain, Z.A.; Kellil, R. Locations and stability of the libration points in the CR3BP with perturbations. J. Math. Anal. 2017, 8, 131–144. [Google Scholar]
  49. Abouelmagd, E.I.; Alzahrani, F.; Hobiny, A.; Guirao, J.L.G.; Alhothuali, M. Periodic solutions around the collinear libration points. J. Nonlinear Sci. Appl. 2016, 9, 1716–1727. [Google Scholar] [CrossRef]
  50. Bekov, A.A. Particular solutions in the restricted collinear three body problem with variable masses. Soviet Astron. 1991, 35, 103–106. [Google Scholar]
  51. Zotos, E.E. How does the oblateness coefficient influence the nature of solutions in the restricted three-body problem? Astrophys. Space Sci. 2015, 358, 33. [Google Scholar] [CrossRef] [Green Version]
  52. Bhatnagar, K.B.; Hallan, P.P. Effect of perturbations in Coriolis and centrifugal forces on the stability of libration points in the restricted problem. Celestial Mech. 1978, 18, 105–112. [Google Scholar] [CrossRef]
  53. Subbarao, P.V.; Sharma, R.K. A note on the stability of the triangular points of equilibrium in the restricted three-body problem. Astron. Astrophys. 1975, 43, 381–383. [Google Scholar]
  54. Gyldén, H. Die bahnbewegungen in einem systeme von zwei Körpern in dem falle, dass die massen veränderungen unterworfen sind. Astron. Nachr. 1884, 109, 1–6. [Google Scholar] [CrossRef] [Green Version]
  55. Bekov, A.A. Periodic solutions of the Gylden-Meshcherskii problem. Astron. Rep. 1993, 37, 651–654. [Google Scholar]
  56. Bekov, A.A.; Momynov, S.B. Parametric solutions of the Gylden-Meshchersky problem. Int. J. Non-Linear Mech. 2019, 116, 195–199. [Google Scholar] [CrossRef]
  57. Meshcherskii, I.V. Works on the Mechanics of Bodies of Variable Mass; GITTL: Moscow, Russia, 1952; p. 205. (In Russian) [Google Scholar]
  58. Kalantonis, V.S.; Perdios, E.A.; Perdiou, A.E.; Vrahatis, M.N. Computing with certainty individual members of families of periodic orbits of a given period. Celestial Mech. Dyn. Astron. 2001, 80, 81–96. [Google Scholar] [CrossRef]
  59. Hu, H.Y. Applied Nonlinear Dynamics; China Aeronautical Industry Press: Beijing, China, 2000. (In Chinese) [Google Scholar]
Figure 1. Zero velocity surfaces for κ = 0.5 , 1 , 2 and 4, respectively.
Figure 1. Zero velocity surfaces for κ = 0.5 , 1 , 2 and 4, respectively.
Universe 06 00110 g001
Figure 2. Zero-velocity surfaces of the third body for κ = 1 : (a) C = 4.16993 , (b) C = 3.59308 , (c) C = 3.57940 , (d) C = 2.82017 .
Figure 2. Zero-velocity surfaces of the third body for κ = 1 : (a) C = 4.16993 , (b) C = 3.59308 , (c) C = 3.57940 , (d) C = 2.82017 .
Universe 06 00110 g002
Figure 3. Zero-velocity curves of the third body: (a) C = 4.16993 , (b) C = 3.59308 , (c) C = 3.57940 , (d) C = 2.82017 .
Figure 3. Zero-velocity curves of the third body: (a) C = 4.16993 , (b) C = 3.59308 , (c) C = 3.57940 , (d) C = 2.82017 .
Universe 06 00110 g003
Figure 4. (ac) Bifurcation diagrams of Coriolis force α and centrifugal force β in the frames α u v , β u v and α β u , respectively, when κ = 1 , μ m = 0.48785 , q 1 = 0.9988 , q 2 = 0.9985 , α 1 = 0.024 , α 2 = 0.02 .
Figure 4. (ac) Bifurcation diagrams of Coriolis force α and centrifugal force β in the frames α u v , β u v and α β u , respectively, when κ = 1 , μ m = 0.48785 , q 1 = 0.9988 , q 2 = 0.9985 , α 1 = 0.024 , α 2 = 0.02 .
Universe 06 00110 g004
Figure 5. (ac) Bifurcation diagrams of oblateness coefficient of the ith ( i = 1 , 2 ) primary in the frames α 1 u v , α 2 u v and α 1 α 2 u , respectively, when κ = 1 , μ m = 0.48785 , q 1 = 0.9988 , q 2 = 0.9985 , α = 1.001 , β = 1.002 .
Figure 5. (ac) Bifurcation diagrams of oblateness coefficient of the ith ( i = 1 , 2 ) primary in the frames α 1 u v , α 2 u v and α 1 α 2 u , respectively, when κ = 1 , μ m = 0.48785 , q 1 = 0.9988 , q 2 = 0.9985 , α = 1.001 , β = 1.002 .
Universe 06 00110 g005
Figure 6. (ac) Bifurcation diagrams of the arbitrary sum κ of the masses of the primaries in the frames κ u v , κ u w and κ v w , respectively, when μ m = 0.48785 , q 1 = 0.9988 , q 2 = 0.9985 , α 1 = 0.024 , α 2 = 0.02 , α = 1.001 , β = 1.002 .
Figure 6. (ac) Bifurcation diagrams of the arbitrary sum κ of the masses of the primaries in the frames κ u v , κ u w and κ v w , respectively, when μ m = 0.48785 , q 1 = 0.9988 , q 2 = 0.9985 , α 1 = 0.024 , α 2 = 0.02 , α = 1.001 , β = 1.002 .
Universe 06 00110 g006
Figure 7. (ac) Bifurcation diagrams of the mass parameter μ m in the frames μ m u v , μ m u w and μ m v w , respectively, when κ = 1 , q 1 = 0.9988 , q 2 = 0.9985 , α 1 = 0.024 , α 2 = 0.02 , α = 1.001 , β = 1.002 .
Figure 7. (ac) Bifurcation diagrams of the mass parameter μ m in the frames μ m u v , μ m u w and μ m v w , respectively, when κ = 1 , q 1 = 0.9988 , q 2 = 0.9985 , α 1 = 0.024 , α 2 = 0.02 , α = 1.001 , β = 1.002 .
Universe 06 00110 g007
Figure 8. (ac) Bifurcation diagrams of radiation factor of the ith ( i = 1 , 2 ) primary in the frames q 1 u v , q 2 u v and q 1 q 2 u , respectively, when κ = 1 , μ m = 0.48785 , α 1 = 0.024 , α 2 = 0.02 , α = 1.001 , β = 1.002 .
Figure 8. (ac) Bifurcation diagrams of radiation factor of the ith ( i = 1 , 2 ) primary in the frames q 1 u v , q 2 u v and q 1 q 2 u , respectively, when κ = 1 , μ m = 0.48785 , α 1 = 0.024 , α 2 = 0.02 , α = 1.001 , β = 1.002 .
Universe 06 00110 g008
Figure 9. The approximate periodic solutions in two-dimensional plane for (a) κ = 0.5 , (b) κ = 1 , (c) κ = 2 , and (d) κ = 4 , respectively.
Figure 9. The approximate periodic solutions in two-dimensional plane for (a) κ = 0.5 , (b) κ = 1 , (c) κ = 2 , and (d) κ = 4 , respectively.
Universe 06 00110 g009
Figure 10. The approximate periodic solutions in three-dimensional space for (a) κ = 0.5 , (b) κ = 1 , (c) κ = 2 , and (d) κ = 4 , respectively.
Figure 10. The approximate periodic solutions in three-dimensional space for (a) κ = 0.5 , (b) κ = 1 , (c) κ = 2 , and (d) κ = 4 , respectively.
Universe 06 00110 g010
Table 1. Physical parameters in Equation (7).
Table 1. Physical parameters in Equation (7).
Parameters μ m q 1 q 2 α 1 α 2 α β
Values0.487850.99880.99850.0240.021.0011.002

Share and Cite

MDPI and ACS Style

Gao, F.; Wang, Y. Approximate Analytical Periodic Solutions to the Restricted Three-Body Problem with Perturbation, Oblateness, Radiation and Varying Mass. Universe 2020, 6, 110. https://doi.org/10.3390/universe6080110

AMA Style

Gao F, Wang Y. Approximate Analytical Periodic Solutions to the Restricted Three-Body Problem with Perturbation, Oblateness, Radiation and Varying Mass. Universe. 2020; 6(8):110. https://doi.org/10.3390/universe6080110

Chicago/Turabian Style

Gao, Fabao, and Yongqing Wang. 2020. "Approximate Analytical Periodic Solutions to the Restricted Three-Body Problem with Perturbation, Oblateness, Radiation and Varying Mass" Universe 6, no. 8: 110. https://doi.org/10.3390/universe6080110

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop