Robe’s Restricted Three-Body Problem with Variable Masses and Perturbing Forces

The linear stability of equilibrium points of a test particle of infinitesimal mass in the framework of Robe’s circular restricted threebody problem, as in Hallan and Rana, together with effect of variation in masses of the primaries with time according to the combined Meshcherskii law, is investigated. It is seen that, due to a small perturbation in the centrifugal force and an arbitrary constant κ of a particular integral of the Gylden-Meshcherskii problem, every point on the line joining the centers of the primaries is an equilibrium point provided they lie within the shell. Further, a number of pairs of equilibrium points lying on the ξζ-plane and forming triangles with the centers of the shell and the second primary exist, for some values of κ > 1. The points collinear with the center of the shell are found to be stable under some conditions and the range of stability depends on the small perturbations and κ, while the triangular points are unstable. Illustrative numerical exploration is given to indicate significant improvement of the problem in Hallan and Rana.


Introduction
In general, the classical restricted three-body problem (RTBP) assumes that the masses of the participating bodies are constant and do not change with time during mechanical motion.The study of two bodies with variable mass has received attention over the years.One of such equation that describes a particular mass variation law is the Gylden-Meshcherskii problem [1,2].This problem extends the classical two-body problem of Newton and Kepler by considering time varying gravitating masses, an important extension of the two body problem for modeling cometary motion and cosmological phenomena and also for describing the evolution of binary stars during secular mass loss owning to photon and corpuscular activity.The slow loss of isotropic mass acts like a perturbation and the problem could be dealt with using a typical method of the theory of perturbations.
Meshcherskii [3] studied the mechanics of the bodies of variable mass.Gelf 'gat [4] examined the restricted problem of three-body of variable masses in which the primary bodies move within the framework of the Gylden-Meshcherskii problem (GMP) and established the existence of five libration points analogous to the classical libration points.Singh and Leke [5] studied the stability of the equilibrium points when the luminous primaries move within the framework of the GMP and vary their masses in accordance with the unified Meshcherskii law.
The classical RTBP assumes that the infinitesimal mass moves under only the mutual gravitational force of the primaries, but in practice Coriolis and centrifugal forces are effective there, and small perturbations affect these forces.Recent interesting studies of this kind of the problem that have been carried out include Hallan and Mangang [6], Singh et al. [7], and Singh and Leke [5].
Robe [8] considered a new kind of restricted three-body problem, in which one of the primaries of mass  1 , is a rigid spherical shell, filled with homogenous, incompressible fluid of density  1 and the other one is a point mass  2 outside the shell and moving around the first primary in Keplerian orbit.A third body of infinitesimal mass  3 is a small solid sphere of density  3 that moves inside the shell.He discussed the linear stability of an equilibrium point of the problem obtained in two cases.In the first case,  2 moves around  1 in circular orbits; in the second case the orbit is elliptic, but the shell is empty or densities of  1 and  3 are equal.Robe's problem can be used to study the small oscillation of the Earth inner core taking into account the Moon's attraction, the stability of the Earth's center, and the motion of artificial satellites in the Earth atmosphere under the action of floating force of the atmosphere.Given its significance and applications, further studies under different assumptions have been carried out since then.
Shrivastava and Garain [9] investigated the effect of small perturbation in the Coriolis and centrifugal forces on the location of libration point in the Robe [8] circular restricted problem of three bodies when the shell is empty.Hallan and Rana [10] studied the linear stability of this same problem and showed that the equilibrium point is linearly stable under some conditions.Later, its nonlinear stability was investigated by Hallan and Mangang [6].
In this paper, we extend the work of Hallan and Rana [10] assuming that the masses of the primaries vary with time in accordance with the unified Meshcherskii [3] law and their motion determined by the GMP.The emergence of more equilibrium points near the center of the shell and infinite number of triangular equilibrium points is observed.

Equations of Motion
The relative motion of the two bodies  1 and  2 whose masses vary with time according to Gylden [1] and Meshcherskii [2] can be described in polar coordinates as where = 0 is a constant of the area integral; θ = () is the angular velocity of revolution of the bodies;  is the gravitational constant while  is the distance between the primaries.Equation (1) has the particular solutions of the type where  1 and  2 are constants given by Also, (1) has a particular integral [4]  =  2 : 0 <  < ∞.
Now, let  * 1 (),  2 (), and  3 be the masses of the first primary (which is a shell and a fluid of density  1 ); second primary and the infinitesimal body, respectively.The masses of the primaries are assumed to vary with time isotropically and the first one contains the infinitesimal mass of density  3 , while the second one is a point mass located outside the first and describes a circular orbit around it.Let us suppose that  is the radius of the shell;  1 ,  2 , and  3 are the centers of the shell,  2 () and  3 , respectively.Let the positions vector between the center of the mass of the first and second primary be ⃗  and those between the infinitesimal body and the primaries ⃗  13 and ⃗  23 , respectively.Following the terminologies of Robe [8], taking into account the gravitational and buoyancy forces exerted by the fluid, and the attraction of  2 (), the equations of motion of the third body of density  3 in a rotating coordinate system under the effects of small perturbations in the Coriolis and centrifugal forces have the forms where Now, (6) does not fully autonomize to a system of equations with constant coefficients when the masses vary according to the combined Meshcherskii law and their motion determined by the GMP, except when the density parameter is zero (i.e.,  = 0).In this regard, the equations of motion of the autonomized system after choosing units of measurements have the forms where where  is the mass parameter and is such that 0 <  < 1 while  physically coincides with the sum of the masses of the primaries.
As  → −∞, () → −∞ and as  → (1 − ), () → +∞; consequently, () vanishes only once in the interval (−∞, 1 − ) and indicates that (13) has only one root in this interval.Now, when there is no perturbation in the centrifugal force, the only solution of ( 13) is  = −.When a small perturbation is given, the solution of (13) can be assumed as Substituting ( 14) in (10), solving and ignoring second-and higher-order term in , we get where These points lie to the right or left of the center of the shell depending on whether the perturbation is positive or negative.Further, due to the centrifugal force and the arbitrary constant  (0 <  < ∞), every point on the line joining the centers of the primaries is an equilibrium point, provided these points lie within the shell.

Triangular Points.
The triangular points of the autonomized systems are the solutions of ( 11) when  = 0,  ̸ = 0, and  ̸ = 1.Solving, we get This gives the position of a pair of equilibrium points which lie in the -plane, forming triangles with the center of the shell and the second primary.These solutions depend on the mass ratio, small perturbation in the centrifugal force, and  and exist only when  > 1. Numerically, it is observed that these equilibrium points do not exist in the Earth-Moon system when  < 1.0036 and  > 1.01 for any value of the perturbation in the centrifugal force chosen very small.We note that these points do not exist in any study of the Robe [8] when the density parameter is zero.Equations ( 12) have no solution, since 0 <  < 1.Hence, we conclude that all the equilibrium points perturbed by small perturbations in the Coriolis and centrifugal forces lie on the -plane.

Stability of Equilibrium Points
The test for stability of the equilibrium points of the nonautonomous system using the definition of a Lyapunov stable solution [11] on the points near the center of the shell yields lim This at once proves that these points are unstable according to Lyapunov's theorem.The same applies to the triangular equilibrium points.
In order to examine the linear stability of equilibrium points of the autonomized system, we displace the infinitesimal mass from its position ( 0 ,  0 ,  0 ) by applying small displacements , V,  to the position (, , ) and then linearize (7) to obtain the variational equations where the partial derivatives are computed at the equilibrium point under consideration.

Equilibrium Points Near the Center of the Shell.
The values of the partial derivatives computed at the equilibrium point (  , 0, 0) are Substituting ( 21) in (20), knowing that  = 1 + ,  = 1 +   at once yields The third equation of ( 22) does not depend on the other two and it shows that the motion parallel to the -axis is stable when  < 1/(1 − (1 + 3  )) and unstable when the reverse holds.Now, the characteristic equation corresponding to the first two equations of ( 22) is given by where Here,  > 0,  <=> 0 when  <=> (4 The roots of (22) are where Now, in the case  =   = 0, (26) is zero only when   = (4/9) + (4 √ 9 − 8/9).Therefore, when the small perturbations are present,  = 0 when where  = 9 − 8. Equation ( 27) gives the various values of , which exist for different values of  and also depends on the small perturbations given.Thus, we denote them by   (see Table 1).These are the critical mass parameters and describe the joint effect of the involved parameters on the stability of motion around the equilibrium point (  , 0, 0).When  =   = 0,  = 1, it becomes the same as that obtained by Robe [8].
Table 1 gives the critical mass values for 0 <  < ∞, and || = 0.001, |  | = 0.002 when the perturbations are considered as positive or negative.We observe that (see Table 1) when  ≤ 0 the critical mass does not exist whether  < 0,  > 0,   > 0,   < 0 but exists for  > 0 and tends to zero as  approaches infinity.Therefore, the values of   increases or decreases with increase in  and whether the perturbations are positive or negative.Also, for 8/9 <  ≤ 1.33, the values increase with increase in , become unity at the peak, and begin to decrease to zero.Hence, stability holds in all the regions of  except for  = 1.33 when  < 0,   > 0 and  = 1.33 when  > 0,   < 0 because in this case  < 1 <   .
Figures 1(a)-1(f) give the graphs of   for  = 0.001,   = 0.002 as  approaches infinity.These graphs indicate that the region of stability is increasing with increase in .
When  = 1,  =   = 0,  increases from 4/9 to 8/9 and to unity; the discriminant increases from −(16/9) to 0 and to 1, respectively.Now, the nature of the roots in (25) depends on the discriminant, small perturbations in the Coriolis and centrifugal forces, mass ratio, and the constant  of a particular integral of the GMP.In view of this, we consider the three regions of  together with the changes in .When  <   , the discriminant is negative; therefore, for  > 0 or  < 0, the real parts of the two of the roots are positive and equal in both cases and the equilibrium points are unstable.When  =   and  < 0, two of the roots are real and equal, while the other two are negative and equal as well, and so the equilibrium points are unstable.In the case when  = 0 all the roots are zero and the points are unstable.Also, when  > 0, all four roots are imaginary, with two positive and equal and the other two negative and equal.In this case we have positive stable resonance.Lastly, when   < , the discriminant is positive and if  > 0, then the roots are distinct and imaginary and the equilibrium points are stable.In the situation  < 0, the roots are real and distinct and the equilibrium points are unstable due to a positive root.Thus, we conclude that the equilibrium points near the center of the shell perturbed by the small perturbation in the centrifugal force and  are unstable for  <   and   ≤  < 1 and stable for   ≤  < 1, depending on the arbitrary constant  of the GMP, the mass ratio , and small perturbations given in the Coriolis and centrifugal forces.Now, we consider the straight line    (Figure 2) whose equation in (,   ) plane is and divides the plane into two parts giving four regions I, II, III, and IV.When  = 1, in (29), we get (43/25)  −(10/3) = 0, which is the same as found by Hallan and Rana [10].Now, firstly for the points in the region I, (29) is positive only when 8/9 <  < 1.34 and negative for 1.34 ≤  < ∞.Therefore, for the point (−,  ) in the region I and 8/9 <  < 1.34,   > (4/9)(1 + √ 9 − 8) and so the range of stability is decreasing as  is increasing.Conversely, in the same region I, when 1.34 ≤  < ∞ then   < (4/9)(1+ √ 9 − 8) and so the range of stability is increasing as  is increasing.Secondly, for the point (,   ) in the region II, (29) is always positive for 8/9 <  < ∞ and the range of stability is decreasing in this region.In the third region III, (29) is negative when 8/9 <  < 1.34 and positive when 1.34 ≤  < ∞.Consequently, for the point (, −  ) the range of stability is increasing and decreasing as  is increasing, respectively.Fourthly, for the region IV containing the point (−,−  ), here, (29) is negative for 8/9 <  < 1, positive for 1 <  < 1.34, and again negative when 1.34 ≤  < ∞.Therefore, in this region, the range of stability increases, decreases, and finally increases as  is increasing.Lastly, for the points on the line   , the range of stability remains unaltered.Now, for the point lying on the -axis,   = 0; therefore, when  > 0, (29) is negative when 8/9 <  < 1.34, and positive for  ≥ 1.34; in these cases, the range of stability is increasing and decreasing, respectively, as  is increasing.Thus, in our case, the Coriolis force is both a stabilizing and destabilizing force when the centrifugal force is kept constant, because of the presence of the constant  of a particular of the GMP.Conversely, for the point on the -axis,   = 0, we see that, when  < 0, (29) is positive when 8/9 <  < 1.34, and negative for  ≥ 1.34; in these cases, the range of stability is decreasing and increasing respectively with increase in , and so the Coriolis force is again both a stabilizing and destabilizing force when the centrifugal force is kept constant.
For the point lying on the   -axis,  = 0; hence, when   > 0, (29) is positive throughout the physically possible interval 8/9 <  < ∞.Here, the range of stability decreases as  is increases, establishing that the centrifugal force is a destabilizing force when the Coriolis force is absent.Similarly, for   < 0, the range of stability increases for 8/9 <  < 1.34, decreases for 1.34 ≤  < 33.9, and afterwards increases for 33.9 ≤  < ∞.Hence, the centrifugal force has both stabilizing and destabilizing behavior when the Coriolis force is kept constant.

Stability of the Triangular Points.
To determine the stability of the triangular equilibrium points, we consider the following partial derivatives: Substituting (30) in the variational equations (20), we obtain the characteristic equation of the triangular points expressed as where The linear stability of the triangular points is determined by the roots of the characteristic equation (31).We explored the equilibrium points, the partial derivatives, and values of the coefficients of equation (31) numerically using the software package Mathematica, for  > 1,  = 1.001,  = 1.002, and 0 <  < 1 and observed that the signs of the quantities  1 ,  2 , and  3 could be negative or positive depending on the interval where  lies.The following cases arise: (i)  1 < 0,  2 > 0, and  3 < 0; there is one change in sign, which implies there is exactly one positive root according to the Descartes rule of sign, (ii)  1 > 0,  2 < 0, and  3 < 0; there is also one change in sign and at least one positive root, (iii)  1 > 0,  2 < 0, and  3 > 0; there are two changes in sign indicating there two positive, two negative, and two complex roots.
A positive root and positive real part of the complex roots induce instability at the triangular point.Hence, we conclude that the triangular points of Robe's problem of two variable mass bodies when the shell is empty are unstable.

Conclusion
We have established the equations of motion of an infinitesimal body in the frame of Robe's [8] circular restricted threebody problem under effects of small perturbations in the Coriolis and centrifugal forces, when the density parameter is zero and the masses of the primaries vary isotropically with time in accordance with the Meshcherskii [3] combined law.The loss or gain of mass to the masses of the primaries acts like a perturbation and the problem is dealt with using the typical perturbation method.The equations of motion (7) of the autonomized system is different from those in Shrivastava and Garain [9] and Hallan and Rana [10].
We also investigated the possible equilibrium points and found that the equilibrium point at the center of the shell found by Robe [8] is shifted due to a small perturbation in the centrifugal force and the parameter  to the left or right of the center of the shell accordingly as whether the perturbation is positive or negative.Also, because of the small perturbation in the centrifugal force, which allows the appearance of the parameter , every point inside the shell on the line joining the centers of the primaries is an equilibrium point and contrary to those of Shrivastava and Garain [9] and Hallan and Rana [10], where there is only one.Further, for some values of the parameter  > 1, a number of pair of equilibrium points forming triangles with the center of the shell and the second primary exists.When  = 1, these points do not exist and the only one equilibrium point is that found by Shrivastava and Garain [9].When there is no perturbation in the centrifugal force, this point fully coincides with that of Robe [8].
The linear stability of the equilibrium point near the center of the shell is unstable for  <   and   ≤  < 1 and stable for   ≤  < 1, depending on the arbitrary constant  of the GMP, the mass ratio , and small perturbations given in the Coriolis and centrifugal forces.When  = 1, the range of stable motion is accordant with Hallan and Rana [10].The increase, decrease, or unchanged range of stability depends on the constant  of the GMP and whether these perturbations are positive or negative.When there are no perturbations the range of stable motion reduces to that of Robe [8].The triangular points are unstable due to a positive root and positive real part of the complex roots of the characteristic equation of sixth degree.
We have assumed that the density parameter is zero, not because previous authors [8][9][10] have done so, but because the equations of motion with variable coefficients does not fully autonomize to that with constant coefficients when the masses vary isotropically and their motion is governed by the GMP.
Modern concept about the formation of celestial bodies and their evolution lead to the necessity of investigating dynamics problems in celestial mechanics in which the motion of the bodies is determined by nonsteady gravitational and nonreactive forces.This study may be used to investigate the dynamic problem in the Earth-Moon system since the distance between them is changing and the masses are also changing due to meteoric activities, space dust, outgassing, and so forth.

Table 1 :
Numerical computation of the critical mass values   .