Abstract
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 . 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 around the binary system.
1. 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 three-body 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–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–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 in the restricted problem has been studied under 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 three-body 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 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 . 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.
2. 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 .
Here, prime () denotes the differentiation with respect to the true anomaly and and denote the partial differentiation of with respect to and , respectively, where and are the oblateness parameter of the primaries and , 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 , , , and , 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 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 , shift the origin at by using the change of variables given as follows: where, the displacement of infinitesimal at and near the equilibrium point 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 and , we obtained where = constant and .
We evaluated , 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
3. 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 . We have constructed out a suitable Hamiltonian to test the stability of the perturbed system by substituting 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 is represented by (20). Thus, the variational equations for the circular problem are given as where where .
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 and
The characteristic equation (26) has been reduced to the following form:
We observed that and , the characteristic equation (27) is representing the classical circular restricted three-body problem:
Taking and , we obtained from (28) Hence, we have that is, The equilibrium position is stable, if are purely imaginary. In order to be purely imaginary, the quantity must be negative and hence the discriminant of (29) is represented as follows: If the equality relation holds in (31), we get that is, that is, Since , the positive sign is inadmissible. Hence, the region of stability in first approximation can be written as Thus, the value of responsible for stable equilibrium points is given by 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 and given by (31) satisfy at least one of the following relations: where is the natural number.
The dependence of and on can be examined by using the simulation technique by drawing the curves by finding out correlation between and . It is clear that in the region (32), only the resonance is possible.
The corresponding value of for is given by
4. Conclusion
We discussed the stability of infinitesimal around the triangular equilibrium points under the elliptical restricted three-body problem, when both of the primaries are oblate spheroid and source of radiation for , 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 has been analyzed using simulation technique of the problem. The region of stability of the linear problem in and 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 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 three-body problem in the particular case when at and near the resonance frequency , 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 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.
Appendix
In order to expand in powers of and , we are required to obtain the expansion of , , , and with the help of Taylor’s series, and we get
In a similar manner, , , and are expanded in the same line of action as it was in the case of as mentioned above, applying Taylor’s theorem, whose coefficients are calculated as follows: Equation (9) can be expressed in the following forms: The Hamiltonian of the perturbed motion of infinitesimal around is obtained after using the values obtained from (14), (15), (16), and (17), substituting in (9), along with substituting the value of , , , and from (A.5), (A.6) in (17), and collecting only fourth degree terms of and , and we get Equating the coefficient of the second order form and the expression of Hamilton function (A.8), we get
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.