# Effects of Radiation Pressure on the Elliptic Restricted Four-Body Problem

Authors: Sahar H. Younis; M. N. Ismail; Ghada F. Mohamdien; A. H. Ibrahiem
Journal: Journal of Applied Mathematics
 (2021)

---

## Abstract

In this paper, under the effects of the largest primary radiation pressure, the elliptic restricted four-body problem is formulated in Hamiltonian form. Moreover, the canonical equations are obtained which are considered as the equations of motion. The Lagrangian points within the frame of the elliptic restricted four-body problem are obtained. The true anomalies are considered as independent variables. An analytical and numerical approach had been used. A code of Mathematica version 12 is constructed to truncate these considerations and is applied on the Earth-Moon-Sun system. In addition, the stability and periodicity of the motion about the equilibrium points are studied by using the Poincare maps. The motion about the collinear point L2 is presented as an example for the obtained results, and some families of periodic orbits are presented.


---

## Body

## 1. Introduction

Several systems in space dynamics such as two-body, three-body, and four-body are considered. There are many attempts to study the problem of these objects using classical methods or Hamiltonian techniques which are good techniques for solving it (see Budrikis [1] and Stone and Leigh [2]). Many authors (Ershkov [3] and Ito [4]) have used these models with various perturbations such as eccentricity, solar radiation, the variable mass of primes, Coriolis and centrifugal forces to checkpoints of vibration, linear stability, and periodic orbits. Singh and Haruna [5] have studied the locations and linear stability of trigonometric points when two primaries radiate, and they studied their heterogeneity. The constrained two-circle four-body problem (4BP) is the simplest model used in the four-body problem (4BP) field Andreu [6]. In this model, the two primers move in circular orbits around the center of mass while this center of mass also moves in a circular orbit around the largest mass in the same plane. Kumari and Kushvah [7] and Ismail et al. [8] used Poincaré’s surface sections to examine the invariance of the constrained four-body problem and applied their study to the Sun-Earth-Moon system.El-Saftawy and Abd El-Salam [9] obtained the solution of the problem using the technique of Delva. Grebenikov et al. [10] studied the problem of the four constrained bodies, by building the Hamiltonian technique, and obtained equilibrium solutions and found six possible equilibrium configurations. Lakin [11] studied the elliptic restricted three-body problem using different coordinate systems to illustrate the stability. Assadian and Pourtakdoust [12] and Chakraborty and Narayan [13] used the problem of the four finite Belial bodies to examine the influence of the Sun on the points of vibration of the Earth-Moon system. Chakraborty and Narayan [14] investigated the equilibrium points, linear stability, zero velocity curves, and fractal trough of the four-body constrained elliptic problem. The model cited in Xu and Fu [15] is useful as a simple example of nonintegral dynamical systems and is successful in understanding many quasiperiodic phenomena in astronomy. Ibrahim [16] studied the effect of solar radiation pressure onto the motion of many-body problem, and applied the study on the Sun-Earth-Moon-spacecraft system.In this paper, we study the problem of the four bodies under the influence of the combined forces of gravity of the elementary bodies that revolve in an elliptical orbit around their center of mass and the pressure of radiation of the largest primary. The problem is formulated in the Hamiltonian form. Then, the canonical form is obtained.The locations of the equilibrium points were determined in terms of the mass ratio of the primariesμ and β the ratio of the radiation strength to the gravitational force of the largest mass. The stability of the motion about the points of vibration under the forces of gravity and radiation is studied. Finally, the canonical equations are applied to the Earth-Moon-Sun system model to calculate the corresponding solutions for both the trigonometric and linear equilibrium points. Finally, the Poincare maps are presented to show these stabilities.


## 2. Description of the Problem

The elliptic restricted four-body problem models represent the motion of the infinitesimal mass,m, under the influence of the gravitational field of the three primaries m1, m2, and m3 and radiation pressure of the third primary m3. The three primaries and the infinitesimal mass m4 are moving in the same plane, i.e., coplanar motion. The two smallest primaries m1 and m2 are moving in an elliptic orbit with e2,f2 about the barycenter CM1 orbiting the largest primary m3 in an elliptic orbit with e1,f1, such that 0<e1, and e2<1. CM is the barycenter between CM1 and m3. The distance between CM1 and m3 is referred to by D, and the distance between m1 and m2 is referred to by d as shown in Figure 1.Figure 1
Present the motion ofm1, m2, m3, and m4.Then, we assume the following two hypotheses: (1)m3≫m1+m2≫m2 and (2) D≫d. Then, the distance between m3 and CM1 from CM is D1 and D2, respectively, and that of m1 and m2 from CM1 is d1 and d2, respectively. We have to keep in our minds that the above hypotheses are in case that the movement of the infinitesimal mass m4 takes place in space as a spatial elliptic restricted four-body problem [17].The main idea in this research is to study the effect of radiation pressure from the more massive bodym3 and to identify the disturbance which will affect the smallest mass m4, so that m3 is considered as perturbing body and the formation of the Hamiltonian system is needed to obtain the canonical form which is considered as the equations of motion.To introduce the dimensionless system, let the distance betweenm1 and m2 equals the unity [18], that means
(1)Gm1+m2=1,μ=m2m1+m2,1−μ=m1m1+m2,μ3=m3m1+m2.Additionally, let the value of the orbital angular momentum with respect to the motion of the primariesm1 and m2 as unity (see Zebehely [17] and Ismail et al. [18, 19]).The coordinates ofm1, m2, m3, and m4 are given as follows:
(2)m1:m1ξ1,η1,ζ1:D2cosf1−d1cosf2,D2sinf1−d1sinf2,0,m2:m2ξ2,η2,ζ2:D2cosf1+d1cosf2,D2sinf1+d1sinf2,0,m4:m4ξ,η,ζ:ξ,η,ζ,m3:m3ξ3,η3,ζ3:−D1cosf1,−D1sinf1,0.Relative to the approximation involved in the assumption of elliptic motion, the distancesD and d are defined by
(3)D=a11−e121+e1cosf1,d=a21−e221+e2cosf2,where a1 and a2 are the semimajor axis of the elliptic orbits; thus, D1= a1m3,D2=d1m1+m2, and D=D1+D2. Similarly, d1=a2m1, d2=a2m2, and d=d1+d2 and f1,f2 are the true anomalies of the bodies m3 and m1, respectively.Now, the system (ξ,η,ζ) will be transformed to the system (x,y,z). The new coordinates of m1, m2, m3, and m4 become
(4)m1:m1x1,y1,z1:−d1,0,0,m2:m2x2,y2,z2:d2,0,0,m3m3x3,y3,z3:−Dscosf2−f1,−Dssinf2−f1,0,where Ds=Rs1−e12/1+e1cosf1.


## 3. The Hamiltonian of the Problem

The equations of motion of the small particlem4 are formulated in Hamiltonian form under the effects of the gravitational potential of small primaries m1 and m2 taking into consideration the perturbing gravitational and radiation of the third large primary m3. The system has a Lagrangian function represented by
(5)Lx,y,z,x2̇,y2̇,ż2=T−V,where T is the total kinetic energy and V is the total potential energy of the system. The potential is defined as
(6)V=1−μr1+μr2+1−βμ3r3,where
(7)r1=x+μ2++y2+z2,r2=x−1+μ2+y2+z2,r3=x−Dscosf1−f22+y−Dssinf1−f22+z2.β is the radiation pressure coefficient of the third body. Then,
(8)L=12ẋ−y2+ẏ+x2+ż2+1−μr1+μr2+1−βμ3r3.The Hamiltonian function defined as(9)H=PTẊ−L,(10)P=∂LẊ.p
x
,
p
y
,
p
z are the components of the momentum P and (ẋ,ẏ,ż) are the components of the velocity Ẋ̇;since the motion is in the plane x−y, then ż=0 and pz=0; from Equations (8), (9), and (10), the Hamiltonian is given by
(11)H=12px2+py2−pxy−pyx−1−μr1+μr2+1−βμ3r3


## 4. The Equations of Motion and Libration Points

The equations of motions are obtained by considering the canonical form of the Hamiltonian which are(12)Ẋ=∂H∂P,(13)Ṗ=−∂H∂X.Now, the libration points are obtained by applying the following conditions:(14)ẋ=0,(15)ẏ=0,(16)pẋ=0,(17)pẏ=0.For collinear libration points, Equations (12), (13), and (14) and (15), (16), and (17) are solved together for x and y=0, while for the nonlinear libration points, Equations (12), (13), and (14) and (15), (16), and (17) are solved together for x,y.


## 5. Results and Discussion

### 5.1. Stability about Equilibrium Points

Since the primaries are moving in elliptical orbits, then to study the motion about any of the libration points and its stability, it is more convenient to introduce the independent variablef2 instead of the independent variable t, so that let
(18)∈=1−e221+e2cosf2,x=∈x¯,y=∈y¯,z=∈z¯,px=∈px¯,py=∈py¯,pz=∈pz¯,H=∈H¯.The new Hamiltonian is given by(19)H¯=ϵ12px¯2+py¯2−px¯y¯−py¯x¯−1ϵ21−μr¯1+μr¯2+1−βμ3r¯3.Then, the equations of motion in the new coordinates are given by(20)x´=σpx¯−y¯−σ1x¯,y´=σpy¯+x¯−σ1y¯,px¯´=−σpy¯+x¯−σ1px¯−σ∈3μ−1+x¯+μy¯2+z¯2+−1+x¯+μ23/2+1−μx¯+μy¯2+z¯2+x¯+μ23/2+−Dscosf1−f2+xμ31−βx¯−Dscosf1−f22+y¯−Dssinf1−f22+z¯23/2,py¯´=−σpx¯−y¯−σ1py¯−σ∈3y¯μy¯2+z¯2+−1+x¯+μ23/2+y¯1−μy¯2+z¯2+x¯+μ23/2+y¯−Dssinf1−f2μ31−βx¯−Dscosf1−f22+y¯−Dssinf1−f22+z¯23/2,where
(21)ḟ=1+e2cosf221−e223/2,∈´=d∈df=e2sinf21−e221+e2cosf22,σ=1ḟ,σ1=∈´∈.It is well known that the characteristic equation for the linearized system is given by(22)λ4+4−Vxx−Vyyλ2+VxxVyy=0,where Vxx is the second derivative of V with respect to x and Vyy is the second derivative of V with respect to y. The four roots of Equation (22) illustrate the stability of motion which depends on the kind of the root if it is real or imaginary or complex numbers. This was treated in more details in many pervious works. Now, in this work, the Poincare map is used to illustrate the stability.


### 5.2. Poincare Map

To apply the obtained equations of motion, some assumptions must be considered: (i) the relation betweenf1 and f2 and (ii) the method used to solve the equations of motion. In this study, the numerical method is used using the explicit Runge-Kutta method. A cod was constructed to do this, and the Poincare map is shown to illustrate the region of the stability of motion about any libration. So that if the Sun-Earth-Moon system is chosen, then it is known that for Sun Ds=Rs1−e12/1+e1cosf1, Rs=328900.48, e=0.016, the true anomaly of orbit Earth-Moon about Sun F takes the values 0 to 2π, μs=328900.48, μ=0.012, and e1=0.0549, the true anomaly of the orbit of Earth-Moon f1 takes the values 0 to 2π. Notice that when the Earth-Moon complete one revolution about the Sun, the Earth-Moon system completes 12 revolutions about their center of mass. That means the true anomaly f2 is considered as a faster variable, so that the change in the Sun’s true anomaly f1 is so slow in comparison with the change in the Earth-Moon true anomaly f2. These notes are taken into consideration when the cod is truncated. β is the radiation pressure coefficient of the Sun and is obtained from β=Frad/Fgr, where Fgr is the gravitational force of the Sun and Frad the radiation force of the Sun. Now, let f1=10°, f2 truncated from 0 to 100π, and β=10−6; these values are used in the cod and are applied on the canonical system (16) and (17), so Figure 2 illustrates the Poincare maps about the libration point Lx2=1.06. In Figure 3, the Poincare map is obtained with f1=170°, to compare the results with those of f1=10° obtained in Figure 2, while Figure 4 illustrates the three dimensions of the Poincare map which appears as elliptical trajectories centered at the libration point Lx, and many periodic trajectories about it are included. In Figure 5, the Poincare map with F=30° is compared with that in Figure 6f1=150°. It is noticed that the change in the value of the true anomaly f1 does not affect the stability of motion about the libration point just changing the phase of map by2π. The same was done for the zero velocity curves which are shown in Figures 7–9 with f1=30°, f1=30°, and f1=150°, respectively. It is noted that the regions about Lx are still unchanged while the energy levels out of these regions are changed.Figure 2
Poincare map about Lx2 withf1=10°.Figure 3
Poincare map about Lx2 withf1=170°.Figure 4
Poincare map about Lx2 withf1=10° in three dimensions.Figure 5
Poincare map about Lx2 withF=30 .Figure 6
Poincare map about Lx2 withF=150 .Figure 7
ZVC about Lx2 withf1=10°.Figure 8
ZVC about Lx2 withf1=30°.Figure 9
ZVC about Lx2 withf1=150°.


## 5.1. Stability about Equilibrium Points

Since the primaries are moving in elliptical orbits, then to study the motion about any of the libration points and its stability, it is more convenient to introduce the independent variablef2 instead of the independent variable t, so that let
(18)∈=1−e221+e2cosf2,x=∈x¯,y=∈y¯,z=∈z¯,px=∈px¯,py=∈py¯,pz=∈pz¯,H=∈H¯.The new Hamiltonian is given by(19)H¯=ϵ12px¯2+py¯2−px¯y¯−py¯x¯−1ϵ21−μr¯1+μr¯2+1−βμ3r¯3.Then, the equations of motion in the new coordinates are given by(20)x´=σpx¯−y¯−σ1x¯,y´=σpy¯+x¯−σ1y¯,px¯´=−σpy¯+x¯−σ1px¯−σ∈3μ−1+x¯+μy¯2+z¯2+−1+x¯+μ23/2+1−μx¯+μy¯2+z¯2+x¯+μ23/2+−Dscosf1−f2+xμ31−βx¯−Dscosf1−f22+y¯−Dssinf1−f22+z¯23/2,py¯´=−σpx¯−y¯−σ1py¯−σ∈3y¯μy¯2+z¯2+−1+x¯+μ23/2+y¯1−μy¯2+z¯2+x¯+μ23/2+y¯−Dssinf1−f2μ31−βx¯−Dscosf1−f22+y¯−Dssinf1−f22+z¯23/2,where
(21)ḟ=1+e2cosf221−e223/2,∈´=d∈df=e2sinf21−e221+e2cosf22,σ=1ḟ,σ1=∈´∈.It is well known that the characteristic equation for the linearized system is given by(22)λ4+4−Vxx−Vyyλ2+VxxVyy=0,where Vxx is the second derivative of V with respect to x and Vyy is the second derivative of V with respect to y. The four roots of Equation (22) illustrate the stability of motion which depends on the kind of the root if it is real or imaginary or complex numbers. This was treated in more details in many pervious works. Now, in this work, the Poincare map is used to illustrate the stability.


## 5.2. Poincare Map

To apply the obtained equations of motion, some assumptions must be considered: (i) the relation betweenf1 and f2 and (ii) the method used to solve the equations of motion. In this study, the numerical method is used using the explicit Runge-Kutta method. A cod was constructed to do this, and the Poincare map is shown to illustrate the region of the stability of motion about any libration. So that if the Sun-Earth-Moon system is chosen, then it is known that for Sun Ds=Rs1−e12/1+e1cosf1, Rs=328900.48, e=0.016, the true anomaly of orbit Earth-Moon about Sun F takes the values 0 to 2π, μs=328900.48, μ=0.012, and e1=0.0549, the true anomaly of the orbit of Earth-Moon f1 takes the values 0 to 2π. Notice that when the Earth-Moon complete one revolution about the Sun, the Earth-Moon system completes 12 revolutions about their center of mass. That means the true anomaly f2 is considered as a faster variable, so that the change in the Sun’s true anomaly f1 is so slow in comparison with the change in the Earth-Moon true anomaly f2. These notes are taken into consideration when the cod is truncated. β is the radiation pressure coefficient of the Sun and is obtained from β=Frad/Fgr, where Fgr is the gravitational force of the Sun and Frad the radiation force of the Sun. Now, let f1=10°, f2 truncated from 0 to 100π, and β=10−6; these values are used in the cod and are applied on the canonical system (16) and (17), so Figure 2 illustrates the Poincare maps about the libration point Lx2=1.06. In Figure 3, the Poincare map is obtained with f1=170°, to compare the results with those of f1=10° obtained in Figure 2, while Figure 4 illustrates the three dimensions of the Poincare map which appears as elliptical trajectories centered at the libration point Lx, and many periodic trajectories about it are included. In Figure 5, the Poincare map with F=30° is compared with that in Figure 6f1=150°. It is noticed that the change in the value of the true anomaly f1 does not affect the stability of motion about the libration point just changing the phase of map by2π. The same was done for the zero velocity curves which are shown in Figures 7–9 with f1=30°, f1=30°, and f1=150°, respectively. It is noted that the regions about Lx are still unchanged while the energy levels out of these regions are changed.Figure 2
Poincare map about Lx2 withf1=10°.Figure 3
Poincare map about Lx2 withf1=170°.Figure 4
Poincare map about Lx2 withf1=10° in three dimensions.Figure 5
Poincare map about Lx2 withF=30 .Figure 6
Poincare map about Lx2 withF=150 .Figure 7
ZVC about Lx2 withf1=10°.Figure 8
ZVC about Lx2 withf1=30°.Figure 9
ZVC about Lx2 withf1=150°.


## 6. Conclusion

In this work, the Hamiltonian of the elliptic four-body problem is studied. The new coordinates which depend on the true anomaly as independent variable are used. The effects of the true anomaly on the stability of motion have a great rule through this study. The sizes of gaps between regions of motion about the libration point are affected by the value of the true anomaly, which means that the position of the primaries in their orbits must be taken into consideration when treating these problems. Also, the energy levels are affected with the values of the true anomaly.



---

*Source: 5842193-2021-11-30.xml*