We consider a version of the relativistic restricted three-body problem (R3BP) which includes the effects of oblateness of the secondary and radiation of the primary. We determine the positions and analyze the stability of the triangular points. We find that these positions are affected by relativistic, oblateness, and radiation factors. It is also seen that both oblateness of the secondary and radiation of the primary reduce the size of stability region. Further, a numerical exploration computing the positions of the triangular points and the critical mass ratio of some binaries systems consisting of the Sun and its planets is given in the tables.
1. Introduction
The problem of three bodies moving under their mutual gravitation attraction has been a subject of keen interest to interested mathematicians since Newton [1] who first postulated the system as three bodies whose forces decrease in a duplicate ratio of the distances, attracting each other mutually; and the accelerative attractions of any two towards the third are between themselves reciprocally as the squares of the distances. Poincare [2] proved that no analytic integrals of motion exist for the general three-body problem other than the energy and angular momentum. This implies that the problem may not be solved in terms of algebraic functions and integrals, but only through numerical integration or infinite series representation.
Regarding the general three-body problem, we may cite Chenciner [3] and Valtonen and Karttunen [4]. There are also various forms of three-body problems in general relativity, for which we may cite Renzetti [5], Nordtvedt [6], and Iorio [7]. If the mass of any one of the three bodies is negligible as compared to the mass of the other two bodies, then it imposes no gravitational influence on the motion of the other two. This then introduces a special form of the three-body problem referred to as the restricted three-body problem (R3BP). The larger bodies are referred to as the primaries and the third body as spacecraft or test particle.
In 1765, Euler [8] found a collinear solution for the restricted three-body problem, where one of the three bodies is a test mass. Soon later, his solutions were extended for a general three-body problem by Lagrange [9] who also found an equilateral triangular solutions in 1772. Now, the solutions for the restricted three-body problem are called Lagrange points L1, L2, L3, L4, and L5. The collinear points L1, L2, and L3 are unstable, they lie on the line joining the primaries, and the triangular points L4 and L5 are stable for the mass ratio μ<0.03852… [10].
The solar and heliospheric observations (SOHO) and WMAP launched by NASA are in operation at L1 and L2 of the Sun-Earth system, respectively. The laser interferometer space antenna (LISA) pathfinder is planned to go to L1. Lagrange points have recently attracted renewed interests for relativistic astrophysics [11, 12], where they have discussed gravitational radiation reaction on L4 and L5 by numerical methods. As a pioneering work, Nordtvedt [6] pointed out that the location of the triangular points is very sensitive to the ratio of the gravitational mass to the inertial one.
The Lagrangian points, by considering one or both primaries are oblate spheroids, whose equatorial planes coincide with the plane of motion are discussed by Subbarao and Sharma [13] and Markellos et al. [14]. Subbarao and Sharma [13] studied the problem in a synodic coordinates system when the massive primary is an oblate spheroid and found that the range of stability decreases with the increase of oblateness factor.
As we know stars (including the Sun) have not only gravitational interaction but also radiation pressure on celestial bodies moving around them. These two actions can be expressed as an equivalent force, called a photogravitation. It is practical and important to study the motion of a small body under the photogravitational restricted three-body problem. Radzievskii [15, 16] formulated the photogravitational restricted three-body problem. Sharma [17] studied the linear stability of triangular libration points of the restricted three-body problem when the more massive primary is a source of radiation and oblate spheroid as well. He found that the eccentricity of the conditional retrograde elliptic periodic orbits around the triangular points at the critical mass μc increases with the increase in the oblateness coefficient and the radiation force and becomes unity when μc is zero.
Simmons et al. [18] gave a complete solution of the restricted three-body. They discussed the existence and linear stability of the equilibrium points for all values of radiation pressures of both luminous bodies and all values of mass ratios.
Sharma [19] studied the stationary solution of the planner restricted three-body problem when the smaller primary is an oblate spheroid with its equatorial plane coinciding with the plane of motion and the bigger primary is a source of radiation as well. He found that the collinear points have conditional retrograde elliptical periodic orbits around them in the linear sense, while the triangular points have long or short periodic retrograde elliptic orbits when parameter of mass is in the range 0≤μ<μcrit, and he deduced that the value of critical mass decreases with the increase of oblateness and radiation force.
AbdulRaheem and Singh [20] studied the stability of equilibrium points under the influence of small perturbations in the Coriolis and the centrifugal forces, together with the effects of oblateness and radiation pressures of primaries, and they found that the Coriolis force has a stabilizing tendency, while the centrifugal force, radiation, and oblateness have destabilizing effects; therefore, the overall effect is that the range of stability of the triangular point decreases.
The existence of libration points and their linear stability as well as periodic orbits around these points when the more massive primary is radiating and smaller is an oblate spheroid studied by Abouelmagd and Sharaf [21]. Their study includes the effects of oblateness of J-2i(i=1,2) with respect to the smaller primary in the restricted three-body problem.
In contrast to the Newtonian theory of gravitation, the general relativity equations of motion of the gravitating masses as the source of the field are intimately related to the field equation as was demonstrated by Einstein and Grommer [22]. This is due to the nonlinear form of the field equations and the existence of the four Bianchi identities.
In general relativity, even writing down the equations of motion in the simplest case N=2 is difficult. Unlike in Newton’s theory, it is impossible to express the acceleration by means of the positions and velocities, in a way which would be valid within the “Exact” theory. Therefore, the approximation method is needed.
Historically, the equations of motion of the problem of N bodies considered to be point masses were first obtained by generalizing the geodesic principle. By the use of this method de Sitter [23] first derived the relativistic equations of motion of N-body problem. Some arithmetic errors occurred in these equations and are reproduced in the encyclopedic paper of Kottler [24] and the treatise by Chazy [25, 26] but were corrected by Eddington and Clark [27].
It is important to highlight some important articles in this field.
Krefetz [28] computed the post-Newtonian deviations of the triangular Lagrangian points from their classical positions in a fixed frame of reference for the first time, but without explicitly stating the equations of motion.
Brumberg [29, 30] studied the problem more in detail and collected most of the important results on relativistic celestial mechanics. He did not obtain only the equations of motion for the general problem of the three bodies but also deduced the equations of motion for the restricted problem of three bodies.
Bhatnagar and Hallan [31] studied the existence and stability of the triangular points L4,5 in the relativistic R3BP, and they concluded that L4,5 are always unstable in the whole range 0≤μ≤1/2 in contrast to the previous results of the classical restricted three-body problem where they are stable for 0≤μ<μ0, where μ is the mass ratio and μ0=0.03852… is Routh’s value.
Douskos and Perdios [32] examined the stability of the triangular points in the relativistic R3BP and, contrary to the results of Bhatnagar and Hallan [31], they obtained a region of linear stability in the parameter space 0≤μ<μ0-1769/486c2 where μ0=0.03852… is Routh’s value.
The locations of libration points in the relativistic R3BP when including one or more additional effects are included in the potential due to radiation pressure and oblateness of the primaries was studied by [33, 34]. The locations of triangular points and their linear stability when the bigger primary is radiating in the relativistic R3BP were examined by Singh and Bello [35]. In all the studies previously mentioned in the relativistic R3BP, no work is performed in the direction of linear stability of triangular points in the presence of both radiation and oblateness.
Our interest in the present paper is not only to obtain the locations of the triangular points but also to study their linear stability when the bigger primary is radiating and smaller one is oblate.
This paper is organized as follows: In Section 2, the equations governing the motion are presented; Section 3 describes the positions of triangular points, while their linear stability is analyzed in Section 4. A discussion of these results and their numerical applications are given in Sections 5 and 6, respectively; finally Section 7 summarizes the conclusions of our paper.
2. Equations of Motion
The pertinent equations of motion of an infinitesimal mass in the relativistic R3BP in barycentric synodic coordinate system ξ,η and dimensionless variables can be written as Brumberg [29] and Bhatnagar and Hallan [31]:(1)ξ¨-2nη˙=∂W∂ξ-ddt∂W∂ξ˙,η¨+2nξ˙=∂W∂η-ddt∂W∂η˙with(2)W=12ξ2+η2+1-μρ1+μρ2+1c2-321-13μ1-μξ2+η2+18ξ˙2+η˙2+2ξη˙-ηξ˙+ξ2+η22+321-μρ1+μρ2ξ˙2+η˙2+2ξη˙-ηξ˙+ξ2+η2-121-μ2ρ12+μ2ρ22+μ1-μ4η˙+72ξ1ρ1-1ρ2-η22μρ13+1-μρ23+-1ρ1ρ2+3μ-22ρ1+1-3μ2ρ2,(3)n=1-32c21-13μ1-μ,(4)ρ12=ξ+μ2+η2,ρ22=ξ+μ-12+η2,where 0<μ≤1/2 is the ratio of the mass of the smaller primary to the total mass of the primaries; ρ1 and ρ2 are distances of the infinitesimal mass from the bigger and smaller primary, respectively; n is the mean motion of the primaries; c is the velocity of light.
We now introduce the oblateness effect of the smaller primary and in the radiation effect of the bigger primary with the help of the parameters A2(0≤A2≪1) and q1(0<1-q1≪1), respectively. Ignoring the second and higher power of A2, we take equations of motion as(5)ξ¨-2nη˙=∂W∂ξ-ddt∂W∂ξ˙,η¨+2nξ˙=∂W∂η-ddt∂W∂η˙,where(6)W=121+3A22ξ2+η2+q11-μρ1+μρ21+A22ρ22+1c2-321-13μ1-μ1+3A24ξ2+η2+18ξ˙2+η˙2+2+3A22ξη˙-ηξ˙+1+3A22ξ2+η22+32q11-μρ1+μρ21+A22ρ22ξ˙2+η˙2+2+3A22ξη˙-ηξ˙+1+3A22ξ2+η2-12q121-μ2ρ12+μ2ρ221+A22ρ222+q1μ1-μ4+3A2η˙+72+214A2ξ1ρ1-1ρ21+A22ρ22-2+3A24η2μρ13+q11-μρ23+1+3A12·-1ρ1ρ2+μ2ρ1-q11-μρ1+q11-μ2ρ2-μρ21+A22ρ22,with the perturbed mean motion(7)n=1+34A2-32c21-13μ1-μ,where A2=AE2-AP2/5R2 [36], AE and AP are, respectively, the equatorial and polar radii and the radiation factor q1 is given by FP1=Fg11-q1 such that 0<1-q1≪1 Radzievskii [15] where Fg1 and FP1 are, respectively, the gravitation and radiation pressure force.
3. Location of Triangular Points
The libration points can be obtained from (5) after putting ξ˙=η˙=ξ¨=η¨=0.
These points are the solutions of the equations:(8)∂W∂ξ=0=∂W∂ηwithξ˙=η˙=0.For simplicity, putting q1=1-1-q1 and neglecting second and higher power of 1-q1 as 0<1-q1=δ≪1 and also the product of A2 and δ, we obtain(9)ξ-1-μξ+μρ13-μξ-1+μρ23+32A2ξ-μξ-1+μρ25+δ1-μξ+μρ13+1c2-3ξ1-μ1-μ3+12ξξ2+η2-32ξ2+η21-μξ+μρ13+μξ-1+μρ23+31-μρ1+μρ2ξ+1-μ2ρ14+μ2ξ-1+μρ24+μ1-μ721ρ1-1ρ2+72ξ-ξ+μρ13+ξ-1+μρ23+32η2μξ+μρ15+1-μξ-1+μρ25+ξ-1+μρ1ρ23+ξ+μρ13ρ2-3μ-2ξ+μ2ρ13-1-3μξ-1+μ2ρ23+32ξ2+η21-μξ+μρ13-31-μξρ1-21-μ2ξ+μρ14δ-μ1-μδ721ρ1-1ρ2+72ξ-ξ+μρ13+ξ-1+μρ23+3η2μξ+μ2ρ15+1-μξ-1+μρ25+ξ-1+μρ1ρ23+ξ+μρ13ρ2+4-5μξ+μ2ρ13+ξ-1+μ-2+4μ2ρ23+A2-941-μ1-μ3ξ+32ξ2+η2ξ+-9μξ-1+μ4ρ25-91-μξ+μ4ρ13-9μξ-1+μ4ρ23ξ2+η2+3μ2ρ23+9μ2ρ2+91-μ2ρ1ξ+2μ2ξ-1+μρ26+μ1-μA2-74ρ23+2141ρ1-1ρ2+72·ξ3ξ-1+μ2ρ25+214ξ-ξ+μρ13+ξ-1+μρ23+94η2μξ+μρ15+1-μξ-1+μρ25+3μξ-1+μ2ρ25+3ξ+μ2ρ13ρ2+3ξ-1+μ2ρ1ρ23-33μ-2ξ+μ4ρ13-31-μξ-1+μ4ρ23+3μξ-1+μ2ρ23=0,ηF=0,with(10)F=1-1-μρ13-μρ23+32A2-3μ2ρ25A2+δ1-μρ13+1c2-31-μ1-μ3+12ξ2+η2+31-μρ1+μρ2-32ξ2+η21-μρ13+μρ23+1-μ2ρ14+μ2ρ24+μ1-μ72·ξ-1ρ13+1ρ23-μρ13+1-μρ23+32η2μρ15+1-μρ25+1ρ13ρ2+1ρ1ρ23-3μ-22ρ13-1-3μ2ρ23+-31-μρ1+32ξ2+η21-μρ13-21-μ2ρ14δ-μ1-μδ72ξ-1ρ13+1ρ23-μρ13-21-μρ23+3η2μ2ρ15+1-μρ251ρ13ρ2+1ρ1ρ23+4-5μ2ρ13+4μ-22ρ23+A2-941-μ1-μ3+32ξ2+η2+-91-μ4ρ13-9μ4ρ23-9μ4ρ25ξ2+η2+3μ2ρ23+91-μ2ρ1+9μ2ρ2+2μ2ρ26+μ1-μA2214ξ1ρ25+1ρ23-1ρ13-32μρ13+1-μρ23+94η2μρ15+1-μρ25+32μρ25+μρ23+1ρ13ρ2+1ρ1ρ23-3μ-22ρ13-1-μ2ρ23. The triangular points are the solutions of (9) with η≠0.
Since 1/c2≪1 and in the case 1/c2→0 and in the absence of oblateness and radiation (i.e., A2=δ=0), one can obtain ρ1=ρ2=1, and we assume in the relativistic R3BP with radiation pressure and oblateness that ρ1=1+x and ρ2=1+y where x,y≪1 may be depending upon the radiation, oblateness, and relativistic terms. Substituting these values in (4), solving them for ξ, η, and ignoring terms of second and higher powers of x and y, we get(11)ξ=x-y+1-2μ2,η=±32+x+y3.Now substituting the values of ρ1, ρ2, ξ, η from the above equations (9) with η≠0 and neglecting second and higher terms in, y, 1/c2, A2, δ, we have(12)321-μ+32A21-μ+δ2μ-1x-3μ2+321+32μA2+δ1-μy+34A21-μ+δ21-μ+1c2-9μ16+27μ216-9μ38+-74+51μ8-109μ216+35μ316δ+34-137μ32+195μ232-63μ316A2=0,31-μ+3μ-1δx+3μ+152μA2y+32·A21-μ+1-μδ+1c2218μ1-μ+32+117μ16-121μ216+3μ32A2+-72+3μ4+39μ28-17μ38δ=0.
Solving these equations for x and y, we get(13)x=-μ2+3μ8c2-12+39-16μ+51μ248c2A2+-13+124c223-22μδ,y=-1-μ5-3μ8c2--15+97μ-68μ2+30μ348μc2A2+124c233-64μ+5μ+26μ2δ.Thus, the coordinates of the triangular points ξ,±η denoted by L4 and L5, respectively, are(14)ξ=1-2μ21+54c2+-12-15-58μ+52μ2+21μ348μc2A2-13+18c2103-14μ+53μ+26μ23δ,η=±321+112c2-5+6μ-6μ2+-36-3-15+136μ-84μ2+81μ3144μc2A2+39-1+14c228-43μ+52μ+13μ2δ.
4. Stability of L4
Let (a,b) be the coordinates of the triangular points L4.
We set ξ=a+α, η=b+β, (α,β≪1), in (5) of motion.
First, we compute the terms of their R.H.S, neglecting second and higher order terms of A2, δ, and their products. We have(15)∂W∂ξξ=a+α,η=b+β=Aα+Bβ+Cα˙+Dβ˙,where(16)A=341+12c22-19μ+19μ2+3+24μ8+-30+367μ-915μ2+717μ3+87μ432μc2A2+124μc2-30-61μ+444μ2-615μ3+234μ4δ,B=3341-2μ1-23c2+-3826μ-7-330-458μ+1498μ2-1533μ3+417μ496μc2·A2+36-1+μ+18μc2-10-26μ+101μ2-9μ3+40μ4δ,C=32c21-2μ+-346μ-1124c2A2+233c2μ-2δ,D=6-5μ+5μ22c2+34-57μ+45μ28c2A2+13c2-5+9μ-6μ2δ.Similarly, we obtain(17)∂W∂ηξ=a+α,η=b+β=A1α+B1β+C1α˙+D1β˙,where(18)A1=3341-2μ1-23c2-326μ-78+330-458μ+1498μ2-1533μ3+417μ496μc2·A2+-31+μ6+12c2-13312+101324μ-338μ2-5312μ+533μ3δ,B1=941+73c2-2+3μ-3μ2+338+30-347μ+915μ2-777μ3+111μ432μc2A2+121-3μ+14c2-454-11μ+52μ+1934μ2-592μ3δ,C1=12c2-4+μ-μ2-16-5μ+9μ28c2A2+13c25-7μδ,D1=-31-2μ2c2+3-11+46μ24c2A2+39c24-20μ+18μ2δ,(19)ddt∂W∂ξ˙ξ=a+α,η=b+β=Fα˙+B2β˙+C2α¨+D2β¨,where(20)F=32c21-2μ-346μ-1124c2A2+239c2μ-2δ,B2=-4-μ+μ22c2-16-5μ+9μ28c2A2+13c25-7μδ,C2=17-2μ+2μ24c2+19-2μ+6μ28c2A2+13c2-8+7μδ,D2=-34c21-2μ+3-1+14μ24c2A2+39c22-μδ,(21)ddt∂W∂η˙ξ=a+α,η=b+β=A3α˙+B3β˙+C3α¨+D3β¨,where(22)A3=6-5μ+5μ22c2+34-57μ+45μ28c2A2+13c2-5+9μ-6μ2δ,B3=-31-2μ2c2+3-11+46μ24c2A2+239c22-10μ+9μ2δ,C3=-31-2μ4c2+3-1+14μ24c2A2+39c22-μδ,D3=35-2μ+2μ24c2+13-6μ+18μ28c2A2+13c29μ-8δ.Thus, the variational equations of motion corresponding to (5), on making use of (7), can be obtained as (23)P1α¨+P2β¨+P3α˙+P4β˙+P5α+P6β=0,q1α¨+q2β¨+q3α˙+q4β˙+q5α+q6β=0,where(24)P1=1+C2,P2=D2,P3=F-C,P4=B2-21+34A2-32c21-13μ1-μ-D,P5=-A,P6=-B,q1=C3,q2=1+D3,q3=21+34A2-32c21-13μ1+μ2-C1+A3,q4=B3-D1,q5=-A1,q6=-B1.Then, the characteristic equation is (25)P1q2-P2q1λ4+P1q6+P5q2+P3q4-P6q1-P2q5-P4q3λ2+P5q6-P6q5=0.Substituting the values of Pi, qi, i=1,2,…,6 in (25), the characteristic equation (25) after normalizing becomes(26)λ4+bλ2+d=0,where(27)b=1-9c2+32-3μ+-101-48μ+51μ2+18μ38c2A2+61-6μ-57μ2+12μ312c2δ,d=27μ1-μ4+9μ-65+77μ-24μ2+12μ38c2+1174μ1-μ+3-61-5289μ+4602μ2-18μ3+846μ464c2·A2+32μ1-μ+-61+397μ-832μ2+240μ3+336μ432c2δ.When 1/c2→0 and in the absence of the oblateness and radiation (i.e., A1=δ=0), (26) reduces to its well-known classical restricted problem form (see, e.g., [10]): (28)λ4+λ2+27μ1-μ4=0.The discriminant of (26) is(29)Δ=-54c2-12698c2A2-42c2δμ4+108c2+638c2A2-28c2δμ3+27+117A2+6δ-6932c2-68018c2A2+1892c2δμ2+-27-123A2-6δ+5852c2+1653916c2A2-4058c2δμ+1+3A2-18c2-653A216c2+42724c2δ.Its roots are(30)λ2=-b±Δ2,where(31)b=1-9c2+32-3μ+-101-48μ+51μ2+18μ38c2A2+61-6μ-57μ2+12μ312c2δ.From (29), we have(32)dΔdμ=4-54c2-12698c2A2-42c2δμ3+3108c2+638c2A2-28c2δμ2+227+117A2+6δ-6932c2-68018c2A2+1892c2δ·μ+-27-123A2-6δ+5852c2+165316c2A2-4058c2δ<0∀μ∈0,12.From (32), it follows that Δ is decreasing in 0,1/2.
But(33)Δμ=0=1+3A2-18c2-653A216c2+42724c2δ>0,Δμ=1/2=-234-117A24-32δ+2074c2+32585A2128c2+47948c2δ<0.Since Δμ=0 and Δμ=1/2 are of opposite signs, and Δ is monotone continuous, there is one value of μ, for example, μc in the interval 0,1/2 for which Δ vanishes. Solving the equation Δ=0, using (29), we obtain critical value of the mass parameter as(34)μc=12-11869+191-1369A2-22769δ-1769486c2+197133+1549369536544c2A2+327543+2026769804816c2δ,μc=μ0-1769486c2+191-1369A2-22769δ+197133+1549369536544c2A2+327543+2026769804816c2δ,where μ0=0.03852… is Routh’s value.
We consider the following three regions of the values of μ separately:
When 0≤μ<μc, the values of λ2 given by (30) are negative and therefore all the four characteristic roots are distinct pure imaginary numbers. Hence, the triangular points are stable.
When μc<μ≤1/2, Δ<0, the real parts of the characteristic roots are positive. Therefore, the triangular points are unstable.
When μ=μc, Δ=0, the values of λ2 given by (30) are the same. This induces instability of the triangular points.
Hence, the stability region is(35)0<μ<μ0-1769486c2+191-1369A2-22769δ+197133+1549369536544c2A2+327543+2026769804816c2δ.
5. Discussion
In this section we discuss the triangular libration points in the framework of relativistic restricted three-body problem when the more massive primary body is a source of radiation and the less massive (secondary) body is oblate.
In analogy to classical problem, the positions and critical mass parameter of analogous triangular libration points (14) and (34), respectively, are obtained. It is important to note that these triangular libration points (14) cease to be classical ones. That is, they no longer form equilateral triangles with the primaries since ρ1≠ρ2≠1. Rather they form scalene triangles but turn into classical ones in the absence of relativistic, radiation, and oblateness factors. Degenerating the present relativistic R3BP to corresponding models, the results of this study are in accordance with those of Douskos and Perdios [32] but in disagreement with those of Bhatnagar and Hallan [31]. In the absence of relativistic factor, it is observed from (25) and (34) that the results of this study correspond to AbdulRaheem and Singh [20] when the bigger primary is luminous only and the smaller one is oblate only in the absence of small perturbations in the Coriolis and centrifugal forces.
In the absence of relativistic terms, our results are also in agreement with those of Abouelmagd and Sharaf [21] when the linear terms, in coefficient B1, are retained only and B2=0, where Bi=J-2iR22i(i=1,2), R2 is the mean radius, and J-2i(i=1,2) are dimensionless coefficients that characterize the size of nonspherical components of the potential of the smaller primary. It is also noticed that, in the absence of oblateness terms, the present results are in agreement with those of Singh and Bello [35]. However, for the locations of triangular points, there is disagreement with those of Katour et al. [34] when the bigger primary is radiating only and the smaller one is oblate only. It is seen from (34) that the effects of relativistic terms, radiation, and oblateness on the linear stability of triangular points are observable. They reduce the size of the stability region independently whereas the joint effect of oblateness and relativistic or radiation and relativistic expands it. This expansion can be explained by the presence of positive coefficients of coupling terms A2/c2 and δ/c2.
6. Numerical Applications
We now use the above mentioned analysis to locate triangular equilibrium points and compute the values of the critical mass parameters for some of Sun-Planet pairs of our solar system. The necessary data have been borrowed from Ragos et al. [37] and Subbarao and Sharma [13]. These values are shown in Tables 1–10.
Locations of triangular points for the Sun-Mars system: c=12424.24, μ=0.000000322700, and A2=0.0000000001×10-8.
δ
ξ
±η
Classical
Relativistic
Equation (14)
Classical
Relativistic
Equation (14)
0.0
0.4999996773
0.4999996813
0.4999996813
0.8660254040
0.8660254016
0.8660254016
0.0001
0.4999996773
0.4999996813
0.4999659297
0.8660254040
0.8660254016
0.8660063980
0.001
0.4999996773
0.4999996813
0.4996621656
0.8660254040
0.8660254016
0.8658353662
0.01
0.4999996773
0.4999996813
0.4966245245
0.8660254040
0.8660254016
0.8641250476
0.02
0.4999996773
0.4999996813
0.4932493676
0.8660254040
0.8660254016
0.8622246936
Critical Mass for the Sun-Mars system: c=12424.24, μ=0.000000322700, and A2=0.0000000001×10-8.
δ
μc
Classical
Relativistic
Equation (34)
0.0
0.0385208965
0.0385208946
0.0385208946
0.0001
0.0385208965
0.0385208946
0.0385200028
0.001
0.0385208965
0.0385208946
0.0385119771
0.01
0.0385208965
0.0385208946
0.0384317199
0.02
0.0385208965
0.0385208946
0.0383425452
Locations of triangular points. Sun-Jupiter system: c=22247.35, μ=0.000953692200, and A2=0.0000192887×10-8.
δ
ξ
±η
Classical
Relativistic
Equation (14)
Classical
Relativistic
Equation (14)
0.0
0.4990463078
0.4990463093
0.4990463093
0.8660254040
0.8660254033
0.8660254033
0.0001
0.4990463078
0.4990463093
0.4990129759
0.8660254040
0.8660254033
0.8660061583
0.001
0.4990463078
0.4990463093
0.4987129755
0.8660254040
0.8660254033
0.8658329534
0.01
0.4990463078
0.4990463093
0.4957129715
0.8660254040
0.8660254033
0.8641009050
0.02
0.4990463078
0.4990463093
0.4923796338
0.8660254040
0.8660254033
0.8621764067
Critical Mass. Sun-Mars system: c=22247.35, μ=0.000000322700, and A2=0.0000000001×10-8.
δ
μc
Classical
Relativistic
Equation (34)
0.0
0.0385208965
0.0385208959
0.0385208959
0.0001
0.0385208965
0.0385208959
0.0385200042
0.001
0.0385208965
0.0385208959
0.0385119784
0.01
0.0385208965
0.0385208959
0.0384317212
0.02
0.0385208965
0.0385208959
0.0383425464
Locations of triangular points. Sun-Saturn system: c=31050.90, μ=0.000285726000, and A2=0.0000018690×10-8.
δ
ξ
±η
Classical
Relativistic
Equation (14)
Classical
Relativistic
Equation (14)
0.0
0.4997142740
0.4997142745
0.4997142745
0.8660254040
0.8660254037
0.8660254037
0.0001
0.4997142740
0.4997142745
0.4996809411
0.8660254040
0.8660254037
0.8660061587
0.001
0.4997142740
0.4997142745
0.4993809404
0.8660254040
0.8660254037
0.8658329539
0.01
0.4997142740
0.4997142745
0.4963809336
0.8660254040
0.8660254037
0.8641009071
0.02
0.4997142740
0.4997142745
0.4930475927
0.8660254040
0.8660254037
0.8621764106
Critical Mass. Sun-Saturn system: c=31050.90, μ=0.00028572600, and A2=0.0000018690×10-8.
δ
μc
Classical
Relativistic
Equation (34)
0.0
0.0385208965
0.0385208961
0.0385208961
0.0001
0.0385208965
0.0385208961
0.0385200044
0.001
0.0385208965
0.0385208961
0.0385119786
0.01
0.0385208965
0.0385208961
0.0384317215
0.02
0.0385208965
0.0385208961
0.0383425467
Locations of triangular points. Sun-Uranus system: c=44056.13, μ=0.000043548000, and A2=0.0000000070×10-8.
δ
ξ
±η
Classical
Relativistic
Equation (14)
Classical
Relativistic
Equation (14)
0.0
0.4999564520
0.4999564525
0.4999564525
0.8660254040
0.8660254038
0.8660254038
0.0001
0.4999564520
0.4999564525
0.4999231189
0.8660254040
0.8660254038
0.8660061590
0.001
0.4999564520
0.4999564525
0.4996231167
0.8660254040
0.8660254038
0.8658329551
0.01
0.4999564520
0.4999564525
0.4966230945
0.8660254040
0.8660254038
0.8641009172
0.02
0.4999564520
0.4999564525
0.4932897365
0.8660254040
0.8660254038
0.8621764305
Critical Mass. Sun-Uranus system: c=44056.13, μ=0.000043548000, and A2=0.0000000070×10-8.
δ
μc
Classical
Relativistic
Equation (34)
0.0
0.0385208965
0.0385208963
0.0385208963
0.0001
0.0385208965
0.0385208963
0.0385200046
0.001
0.0385208965
0.0385208963
0.0385119788
0.01
0.0385208966
0.0385208963
0.0384317216
0.02
0.0385208965
0.0385208963
0.0383425469
Locations of triangular points. Sun-Neptune system: c=55148.85, μ=0.000051668900, and A2=0.0000000010×10-8.
δ
ξ
±η
Classical
Relativistic
Equation (14)
Classical
Relativistic
Equation (14)
0.0
0.4999483311
0.4999483311
0.4999483311
0.8660254040
0.8660254040
0.8660254040
0.0001
0.4999483311
0.4999483311
0.4999149976
0.8660254040
0.8660254040
0.8660061590
0.001
0.4999483311
0.4999483311
0.4996149964
0.8660254040
0.8660254040
0.8658329546
0.01
0.4999483311
0.4999483311
0.4966149845
0.8660254040
0.8660254040
0.8641009107
0.02
0.4999483311
0.4999483311
0.4932816379
0.8660254040
0.8660254040
0.8641009107
Critical Mass. Sun-Neptune system: c=55148.85, μ=0.000051668900, and A2=0.0000000010×10-8.
δ
μc
Classical
Relativistic
Equation (34)
0.0
0.0385208965
0.0385208964
0.0385208964
0.0001
0.0385208965
0.0385208964
0.0385200047
0.001
0.0385208965
0.0385208964
0.0385119789
0.01
0.0385208966
0.0385208964
0.0384317217
0.02
0.0385208965
0.0385208964
0.0383425469
7. Conclusion
Linear stability of triangular points has been studied when the more massive primary is radiating and the smaller one is an oblate spheroid well. It is discovered that their locations and critical mass are affected by radiation and oblateness factors. It is also observed that both factors have destabilizing behavior; that is, an increase in their values reduces the size of stability region independently. However the joint effect of radiation, oblateness, and relativistic terms expands the size of the stability region. It is also noticed that the expressions for A, D, A2, and C2 in Bhatnagar and Hallan [31] differ from the present study when the radiation pressure and oblateness are absent. Consequently, the expression of P1, P3, P4, and P5 and the characteristic equation are also different. This led them [31] to infer that triangular points are unstable, contrary to Douskos and Perdios [32] and our results.
Our results are also in disagreement with those of Katour et al. [34]. One major distinction is that the expression of the mean motion n which they used in their study differs from our own. It seems that in their expression for the mean motion. In addition to that they have not included the coupling term A2/c2 in their study, while we do. Also, they have not analyzed stability, while we do.
It is observed from Tables 1–10 that the numerical exploration conducted on some of the Sun-Planet pairs reveals that the effect of radiation is prominent in the size of stability region and the positions, while their oblateness effects are comparatively negligible.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
NewtonS. I.1729Translated into English by A. MottePoincareH.1892Paris, FranceGauthier-VillarsChencinerA.Three body problem2007210, article 211110.4249/scholarpedia.2111ValtonenM.KarttunenH.2006Cambridge, UKCambridge University Press10.1017/cbo9780511616006MR2223553RenzettiG.Exact geodetic precession of the orbit of a two-body gyroscope in geodesic motion about a third mass20121091–4555910.1007/s11038-012-9402-22-s2.0-84871290768NordtvedtK.Jr.Testing relativity with laser ranging to the moon196817051186118710.1103/PhysRev.170.11862-s2.0-36049057836IorioL.Orbital motions as gradiometers for post-Newtonian tidal effects20141, article 310.3389/fspas.2014.00003EulerL.De motu rectilineo trium corporum se mutuo attrahention176511144LagrangeJ. L.Essai sur le probléme des trois corps17726272292SzebehelyV.1967New York, NY, USAAcademic PressSetoN.MutoT.Relativistic astrophysics with resonant multiple inspirals20108110610300410.1103/PhysRevD.81.103004SchnittmanJ. D.The Lagrange equilibrium points L4 and L5 in black hole binary system201072413910.1088/0004-637X/724/1/39SubbaraoP. V.SharmaR. K.A note on the stability of the triangular points of equilibrium in the restricted three-body problem1975433381383MarkellosV. V.PapadakisK. E.PerdiosE. A.Non-linear stability zones around triangular equilibria in the plane circular restricted three-body problem with oblateness1996245115716410.1007/bf006378112-s2.0-0007113745RadzievskiiV. V.The restricted problem of three-bodies taking account of light pressure195027250256MR0038750RadzievskiiV. V.The spatial case of the restricted problem of three radiating and gravitating bodies195330265273MR0058330SharmaR. K.FrickleW.TelekiG.On linear stability of triangular libration points of the photogravitational restricted three-body problem when the more massive primary is an oblate spheroid1982435Reidel Publishing Company345436SimmonsJ. F. L.McDonaldA. J. C.BrownJ. C.The restricted 3-body with radiation pressure198535214518710.1007/bf01227667MR7869252-s2.0-0000181347SharmaR. K.The linear stability of libration points of the photogravitational restricted three-body problem when the smaller primary is an oblate spheroid1987135227128110.1007/BF00641562ZBL0645.700062-s2.0-0003069521AbdulRaheemA.SinghJ.Combined effects of perturbations, radiation, and oblateness on the stability of equilibrium points in the restricted three-body problem200613131880188510.1086/4993002-s2.0-33645238621AbouelmagdE. I.SharafM. A.The motion around the libration points in the restricted three-body problem with the effect of radiation and oblateness2013344232133210.1007/s10509-012-1335-8ZBL1276.850032-s2.0-84874571700EinsteinA. S.GrommerJ.19271-2Berlin , GermanyWalter de GruyterSitzungsberichte der Preussischen Akademie der Wissenschaftende SitterW.On Einstein's theory of garviation and its astronomical consequences191676699728KottlerF.Gravitation und relativitts theorie192262, article 22aChazyJ.19281Paris, FranceGauthier-VillarsChazyJ.19301Paris, FranceGauthier-VillarsEddingtonA. S.ClarkG. L.The problem of n bodies in general relativity theory1938166465KrefetzE.Restricted three-body problem in the post-Newtonian approximation19677247147310.1086/110252BrumbergV. A.1972Moscow, RussiaNaukaBrumbergV. A.1991New York, NY, USAHilgerMR1200453BhatnagarK. B.HallanP. P.Existence and stability of L4,5 in the relativistic restricted three-body problem199869327128110.1023/a:1008271021060MR1645113DouskosC. N.PerdiosE. A.On the stability of equilibrium points in the relativistic restricted three-body problem200282431732110.1023/a:1015296327786MR19057702-s2.0-0036099820Abd El-SalamF. A.Abd El-BarS. E.On the triangular equilibrium points in the photogravitational relativistic restricted three-body problem2014349112513510.1007/s10509-013-1629-52-s2.0-84890343517KatourD. A.Abd El-SalamF. A.ShakerM. O.Relativistic restricted three body problem with oblatness and photo-gravitational corrections to triangular equilibrium points2014351114314910.1007/s10509-014-1826-x2-s2.0-84897539656SinghJ.BelloN.Effect of radiation pressure on the stability of L4,5 in the relativistic R3BP2014351248349010.1007/s10509-014-1858-22-s2.0-84899947844McCuskeyS. W.1963New York, NY, USAAddison-WesleyRagosO.PerdiosE. A.KalantonisV. S.VrahatisM. N.On the equilibrium points of the relativistic restricted three-body problem20014753413341810.1016/s0362-546x(01)00456-4MR19792352-s2.0-0035418866