Existence of Resonance Stability of Triangular Equilibrium Points in Circular Case of the Planar Elliptical Restricted Three-Body Problem under the Oblate and Radiating Primaries around the Binary System

This paper analyzes the existence of resonance stability of the triangular equilibrium points of the planar elliptical restricted three-body problem when both the primaries are oblate spheroid as well as the source of radiation under the particular case, when e = 0. We have derived Hamiltonian function describing the motion of infinitesimal mass in the neighborhood of the triangular equilibrium solutions taken as a convergent series. Hamiltonian function for the system has been derived and also expanded in powers of the generalized components of momenta. We have used canonical transformation to make the Hamiltonian function independent of true anomaly. The most interesting and distinguishable results of this study are establishing the relation for determining the range of stability at and near the resonance ω 2 = 1/2 around the binary system.


Introduction
The elliptical restricted three-body problem describes the dynamical system more accurately on account of realistic assumptions of the motion of the primaries that are subjected to move along the elliptical orbit.If the primaries move along elliptical orbits, this can have significant effects because the orbits of the most of planets and stars are elliptical rather than circular, so this problem is turned as elliptical restricted threebody problem.The celestial bodies as we know are not perfectly spheres.They are either oblate spheroid or triaxial rigid body.The lack of sphericity of the heavenly bodies causes large perturbations from a two-body orbit.This motivates our studies of stability of triangular equilibrium points under the resonance, when both of the primaries are oblate spheroid and source of radiation.The stability of infinitesimal around the equilibrium points of the elliptical restricted three bodies has been studied [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20] and also the authors have investigated the different aspects of the same problem.The existence of the libration points and their stability in the radiational elliptical restricted three-body problem has been investigated [14].The stability of the motion of infinitesimal around the triangular equilibrium points is depending upon  and .The different aspects of the same problem in details have been investigated by [16,[21][22][23].The existence of libration points and their stability in the three-body photogravitational elliptical restricted problem have been studied [11].The analytical investigation concerning the structure of asymptotic perturbative approximation for small amplitude motions of the third-point mass in the neighborhood of a Lagrangian equilateral libration positions in the planar and elliptical restricted three bodies have been investigated [24].After a sequence of canonical transformations, they formulated the Hamiltonian function governing the motion of the negligible mass body using the eccentric anomaly of the primaries elliptical Keplerian orbit as the independent variable.They studied the liberalized system of differential equations of motion obtained from expanding the Hamiltonian around a Lagrangian solution.The approximated integrations of the elliptical restricted three-body problem by means of perturbation technique based on Lie series and development, which led to an approximated solution of the differential system of canonical equations of motion derived from the chosen Hamiltonian function, have been discussed [25].The present study aims to investigate the condition of existence of resonance stability of infinitesimal body around the triangular equilibrium points in the elliptical restricted three-body problem, when both of the primaries are oblate spheroid and also the source of radiation.We have adopted the method due to Markeev [21], in which the Hamiltonian function pertaining to the problem is made independent of time using several canonical transformations.The existence of resonance and the stability of infinitesimal near the resonance frequency has been analyzed.
The stability of the triangular equilibrium points in the circular restricted three-body problem considering the bigger primary as an oblate spheroid in linear cas, has been studied and established that range of the mass parameter given rise to stable triangular solutions is lowered [12].The values of critical mass ratio have been obtained for the various values of the oblateness parameters in this paper.
The nonlinear stability of the triangular equilibrium point  4 in the restricted problem has been studied under x 0 0.4563 0.5913 0.6053 0.9 1.22 1.18 the condition that the potential between the bodies is considered as the potential forces.In the first case they discussed the nonlinear stability of the same problem by taking into account of various conditions such that the bigger primary is an oblate spheroid and in the second case both the primaries are oblate spheroid, whose axes of symmetry is perpendicular to the plane of relative motion of the primaries.They also discussed the nonlinear stability of triangular equilibrium points in circular restricted three bodies as the primaries are spherical in shape and the bigger primary is a source of radiation [26].
The nonlinear stability zones of the triangular Lagrangian points are computed numerically, where the bigger primary is an oblate spheroid in the planer circular restricted threebody problem.It is found that the oblateness has a noticeable effect and this is identified to be related to the resonance case and the associated curves in the mass parameters  verses oblateness coefficient  1 configuration space and represented this measure in the parameter space.The location of triangular equilibrium points around the binary system is represented by Figures 1, 2, 3, 4, and 5.
The present paper described the effects of the oblateness and radiation of both the primaries on the existence of resonance and the stability of the triangular equilibrium points of the planar elliptical restricted three-body problem in particular case when  = 0. We have adopted the method due to Markeev [21,27] which made the Hamiltonian function independent of time by using several canonical transformations.The existence of resonance stability of the triangular equilibrium points near the resonance frequency has been analyzed using simulation technique.

Equation of Motion
The equation of motion of the planar elliptical restricted three-body problem under the photogravitational and oblateness of both the primaries in barycentric, pulsating, nondimensional coordinates is represented by [28]: where where and  = 1/(1 +  cos V).
Here, prime (  ) denotes the differentiation with respect to the true anomaly and Ω  and Ω  denote the partial differentiation of Ω with respect to  and , respectively, where  1 and  2 are the oblateness parameter of the primaries and  1 ,  2 are the source of radiation of both the primaries.There are two triangular equilibrium points in the plane of the finite bodies around the coordinates (, ) and then the three bodies form nearly equilateral triangle.Since the equilibrium points are symmetrical to each other, the nature of motion near the two triangular equilibrium points is the same Thus, the coordinate of the triangular equilibrium points has been obtained up to the first order terms in the parameters  1 ,  2 ,  1 , and  2 , which is represented by (4) referred from [22,29].
The system (1) described the motion of dynamical system with Lagrangian which is represented as We constructed the expression for the Hamiltonian function of the problem using the formula Hence, we have The perturbed Hamiltonian function of the given problem can be obtained using ( 5) and (8), which is reduced to the following form: Since, the two triangular equilibrium points are symmetrical, the nature of the oscillations of infinitesimal near the two Advances in Astronomy points will be the same.We investigated the motion of infinitesimal near the equilibrium point.In order to discuss the motion of infinitesimal near the equilibrium point  4 , shift the origin at  4 by using the change of variables given as follows: where, the displacement of infinitesimal at and near the equilibrium point  4 is represented as follows: The solution (11) in the new variables is given by the equilibrium position: Expanding the Hamiltonian function (9) in the powers of  1 and  1 , we obtained where  0 = (, ,   ,   ) = constant and  1 = 0. We evaluated  2 ,  3 , . .., using (10), and the terms of the Hamiltonian ( 9) are expanded one by one, the terms are not depending upon   and   , and those of order one are neglected.Thus, we obtain where where

Stability of Triangular Equilibrium Points of the Problem in Circular Case
We discussed the stability of infinitesimal under elliptical restricted three-body problem, when both of the primaries are oblate spheroid and the source of radiation under the particular case, when  = 0. We have constructed out a suitable Hamiltonian to test the stability of the perturbed system by substituting  = 0 in the expression (A.9) (see the appendix); we thus obtained the Hamiltonian for the circular problem up to the second order terms (see the appendix), which is represented as follows: The variational equation for the circular problem can be written as where the Hamiltonian of second order  2 is represented by (20).Thus, the variational equations for the circular problem are given as where  (2) −  (2) + 13 4  (1) ) , where  = 3 √ 3(1 − 2)/4.Now, the characteristic equation based on the problem can be obtained on putting The characteristic equation is obtained by substituting the results obtained from ( 24) in ( 22), and we get Solving the characteristic equation, we get where  * +  − 4 = −1 and The characteristic equation (26) has been reduced to the following form: We observed that  1 = 0 and  2 = 0, the characteristic equation ( 27) is representing the classical circular restricted three-body problem: Taking  1 =  1 and  2 =  2 , we obtained from ( 28) Hence, we have that is, ] . (31) The equilibrium position is stable, if  1,2 are purely imaginary.
(37) Also, we observed that in the case when The existence of resonance is possible in the neighborhood of the value of  for which the frequencies  1 and  2 given by (31) satisfy at least one of the following relations: where  is the natural number.
The dependence of  1 and  2 on  can be examined by using the simulation technique by drawing the curves by finding out correlation between  and  1,2 .It is clear that in the region (32), only the resonance  2 = 1/2 is possible.

Conclusion
We discussed the stability of infinitesimal around the triangular equilibrium points under the elliptical restricted threebody problem, when both of the primaries are oblate spheroid and source of radiation for  = 0, in circular case near the resonance frequency.We have constructed a suitable normalized convergent Hamiltonian function and investigated the stability of infinitesimal around the triangular equilibrium points and the perturbed system analytically and numerically due to oblateness and radiation of primaries in circular case up to the second order terms.The region of stability and instability of the linear problem for the value of eccentricity "" = 0 has been analyzed using simulation technique of the problem.The region of stability of the linear problem in  and  1,2 in the plane has been clearly marked.In the region outside the boundary, the triangular points are unstable while, in the region inside the boundary curve, the necessary condition of linear stability is satisfied in  and  1,2 plane.We concluded that the effect of the oblateness and radiation of primaries affects the location and resonance stability of triangular equilibrium points of the elliptical restricted threebody problem in the particular case when  = 0 at and near the resonance frequency  2 = 1/2, which is analyzed from the graphical behavior of the triangular equilibrium points around the binary system.Hence, we arrived at the conclusion that the shifting of the locations and the resonance stability of triangular equilibrium points in the particular case when  = 0 around the well-known binary system would be possible by changing the oblateness and radiation parameters which are observed from Figures 6,7,8,10,11,12,14,15, and 16 and the combined Figures 9,13,17,18,19,20,and 21.Nonlinear behavior of infinitesimal around the triangular equilibrium points under the above mentioned perturbing forces is yet to be investigated.

Figure 3 :
Figure 3: Location of triangular points for  Cen-AB.
Figure 14: Correlation between  and  1,2 for binary system Achird-I.Figure 15: Correlation between  and  1,2 for binary system Achird-II.Figure 16: Correlation between  and  1,2 for binary system Achird-III.Figure 17: Combined graph for binary systems Achird-I, II, and III.Figure 18: Combined graph for binary systems Luyten-I, II, and III.
Figure 21: Combined graph for binary systems Xi-Booties-I, II, and III.