Pulsating Different Curves of Zero Velocity around Triangular Equilibrium Points in Elliptical Restricted Three-Body Problem

The oblateness and the photogravitational effects of both the primaries on the location and the stability of the triangular equilibrium points in the elliptical restricted three-body problem have been discussed. The stability of the triangular points under the photogravitational and oblateness effects of both the primaries around the binary systems Achird, Lyeten, Alpha Cen-AB, Kruger 60, and Xi-Bootis, has been studied using simulation techniques by drawing different curves of zero velocity.


Introduction
The present paper is devoted to the analysis of the photogravitational and the oblateness effects of both primaries on the stability of triangular equilibrium points of the planar elliptical restricted three-body problem.The elliptical restricted three-body problem describes the dynamical system more accurately on account that the realistic assumptions of the motion of the primaries are subjected to move along the elliptical orbit.We have attempted to investigate the stability of triangular equilibrium points under the photogravitational and oblateness effects of both the primaries.The bodies of the elliptical restricted three-body problem are generally considered to be spherical in shape, but in actual situations, we have observed that several heavenly bodies are either oblate spheroid or triaxial rigid bodies.The Earth, Jupiter, and Saturn are examples of the oblate spheroid.The lack of sphericity in heavenly bodies causes large perturbation.In addition to the oblateness of heavenly bodies, the triaxiality, the radiation forces of the bodies, the atmospheric drag, and the solar wind are also causes of perturbation.
This motivates studies of stability of triangular equilibrium points under the influence of oblateness and radiation of the primaries in the elliptical restricted three-body problem.The stability of the infinitesimal around the triangular equilibrium points in the elliptical restricted three-body problem described in considerable details is due to [1] and the problem was also studied [2][3][4][5][6][7][8][9].The stability of motion of infinitesimal around one of the triangular equilibrium points ( 4 ) also depends on  and .Nonlinear stability of the triangular equilibrium points of the elliptical restricted three-body problem with or without radiation pressure was studied [10][11][12].Furthermore, the nonlinear stability of the infinitesimal in the orbits or the size of the stable region around  4 was studied numerically by [11] and the parametric resonance stability around  4 in the elliptical restricted threebody problem has been studied [10].
The influence of the eccentricity of the orbits of the oblate primary bodies with one of the photogravitational effects on the location of collinear and triangular equilibrium points and their stability has been investigated [27][28][29].The stability of triangular points in the elliptical restricted three-body problem under the radiating and oblate primaries has been recently discussed [30].A similar problem has been discussed in detail by applying different techniques to investigate the stability of the system using simulation technique [31].
The present study aims to examine the motion of the infinitesimal body in the elliptical restricted three-body problem, when the primaries are oblate spheroid and are also a source of radiation.We have obtained the coordinate of the triangular equilibrium points of the problem.For the circular problem, the primaries are fixed with respect to uniformly rotating axes and hence the Hamiltonian does not involve time explicitly.But when the primaries move on elliptical orbits, the introduction of nonuniformly rotating and pulsating coordinate system results again in fixed location of the primaries.The elliptical restricted three-body problem generalizes the original circular restricted three-body problems, while some useful problem of circular model still can be satisfied by the elliptical case.The Hamiltonian, however, does not depend explicitly on independent variable in this case.The dimensionless variables are introduced by using the distance  between primaries given by where  and  are the semimajor axis and the eccentricity of the elliptical orbit of the primaries, moving along each other and V is the true anomaly of  1 .A co-ordinate system which rotates with the variable angular velocity  is introduced.This angular velocity is given by where  * is dimensionless time.
The equation follows from the principal of the conservation of angular momentum in the problem of two bodies formed by the primaries of masses  1 and  2 .This principle is expressed by where  =  1 +  2 , and  1 and  2 are the products of the universal gravitational constants with the masses of primaries.
The force of radiation is given by where   is the gravitational attraction force;   is the radiation pressure;  is the mass radiation factor.For simplicity of calculation, we have considered  = 1 − ∈  ,  = 1, 2, where The present paper deals with the photogravitational and oblateness effects of both the primaries on the stability of infinitesimal by exploiting simulation techniques by drawing the different curves of zero velocity.
The present paper comprises three sections.The first section of the paper describes the variational equation of motion of the system.In the second section, we have described the triangular equilibrium points of the system, and in the third section the paper we have derived an expression which is suitable for tracing of different curves of zero velocity.
The curves of zero velocity around equilibrium points have been presented through simulation techniques, which shows the region of stability.The effect of the oblateness of primaries plays an important role in analyzing the stability of infinitesimal which is obvious from the curves of zero velocity traced.

Variational Equation of Motion
The differential equations of motion of infinitesimal in the elliptical restricted three-body problem under the oblate and radiating primaries in barycentric, pulsating, and nondimensional coordinates are represented as follows [31]: where where where Ω  denotes the partial differentiation of Ω with respect to , and Ω  denotes the differentiation of Ω partially with respect to, where  1 and  2 are the oblateness parameters of the primaries. 1 and  2 are the mass radiation factors due to the source of radiation of the bigger primary and the smaller primary.
The co-ordinates of the triangular equilibrium points  4 and  5 are determined as follows [31]: Journal of Mathematics  Thus by the coordinates of the triangular equilibrium points we obtain up to the first-order terms in the parameter  1 ,  2 ,  1 , and  2 , which is represented by (9), the location of equilibrium points is shown Figures 1, 2, 3, 4, and 5 for the infinitesimal moving around various binary systems.

Different Curves of Zero Velocity
In order to discuss the different curves of zero velocity of the infinitesimal in the elliptical restricted three-body problem, when both the primaries are oblate spheroid and radiating, multiplying the first equation of ( 5) by   and the second equation by   and adding, we get; we obtain Since Ω does not contain the time (true anomaly) explicitly, (10) can be integrated to give  Due to the presence of (1+ cos V) in the denominator of (11), the equation is not possible to integrate to any defined form.Hence, in elliptical restricted three-body problem, it does not adjust the Jacobi integral of the classical circular problem at least in its usual sense.
The elliptical restricted three-body problem is different from the classical restricted problem in the sense that the Jacobi integral does not exist [16], and energy along any orbit is a time-dependent quantity.As we know no exact, complete, and general solution to the elliptical restricted three-body problem, Ω can be obtained unlike in classical restricted three-body problem, but this mathematical inconvenience is overcome along investigation of certain special cases of the problem based on simplifying the mathematical model under consideration [2].Now, consider the potential function which is represented as follows: Hence, Ω(, ) depends not only on the position coordinate of the infinitesimal but also on an independent variable.We select the initial point V = 0 and we consider only a part of the trajectory V = 0 and V = , where  is arbitrary sufficiently small time interval, during which the primaries describe the small one.We may define a Jacobi constant in elliptical case as follows: Equation ( 13) describes different curves of zero velocity, at each given instant of time of elliptical restricted three-body problem.The zero velocity curves are now pulsating with frequency of the nominal elliptical motion.Therefore, in the planar elliptical restricted three-body problem, the zero velocity curves are obtained from the following equation: with the help of ( 14), the different curves of zero velocity have been traced using the software MATLAB 7.1 of the infinitesimal around the binary system Achird, Luyten, Alpha Cen AB, Kruger 60 and Xi Bootis, taking into account values of  oblateness parameters  1 and  2 and critical mass parameters  1 and  2 from Table 1.We have traced different curves of zero velocity of the infinitesimal about triangular equilibrium points.Figures 6, 7, and 8 represent the same around binary system Achird, and likewise Figures 9, 10, and 11 for binary system Luyten, Figures 12, 13, and 14 for binary system Alpha Cen-AB, Figures 15, 16, and 17 for binary system Kruger 60, and Figures 18, 19, and 20 for binary system Xi Bootis.
Hence, we observe typical behavior of the infinitesimal around the binary system, Achird, Luyten, Alpha Cen-AB, Kruger 60, and Xi Bootis.

Discussion and Conclusion
The oblateness and photogravitational effects of the primaries on the location and the stability of the triangular Lagrangian points in elliptical restricted three-body problem have been studied.The problem was studied under the assumption that the eccentricity of the orbit of the gravitating bodies is small.The oblateness of the more massive primary does not affect the motion of the smaller primary due to its large mass, whereas it affects the motion of infinitesimal body.
The differential equation governing the motion and stability of triangular equilibrium points of the elliptical restricted three-body problem under the oblate and radiating primaries has been analyzed, and configurations of the triangular equilibrium points are described.The stability of the triangular points under the photogravitational and oblateness effects of both the primaries around the binary systems Achird, Lyeten, Aipha Cen-AB, Kruger 60, and Xi Bootis have been studied using simulation technique by drawing different curves   of zero velocity around triangular equilibrium point.It is observed that the region within the curves, the infinitesimal will remain stable.

Table 1 :
[30]tion of triangular equilibrium points mentioned in the research paper of Jagdish and Aishetu[30].