On Robe ’ s Circular Restricted Problem of Three Variable Mass Bodies

This paper investigates the motion of a test particle around the equilibrium points under the setup of the Robe’s circular restricted three-body problem in which themasses of the three bodies vary arbitrarily with time at the same rate.The first primary is assumed to be a fluid in the shape of a spherewhose density also varies with time.Thenonautonomous equations are derived and transformed to the autonomized form. Two collinear equilibrium points exist, with one positioned at the center of the fluid while the other exists for the mass ratio and density parameter provided the density parameter assumes value greater than one. Further, circular equilibrium points exist and pairs of out-of-plane equilibrium points forming triangles with the centers of the primaries are found. The out-of-plane points depend on the arbitrary constant κ, of the motion of the primaries, density ratio, and mass parameter. The linear stability of the equilibrium points is studied and it is seen that the circular and out-of-plane equilibrium points are unstable while the collinear equilibrium points are stable under some conditions. A numerical example regarding out-of-plane points is given in the case of the Earth, Moon, and submarine system.This study may be useful in the investigations of dynamic problem of the “ocean planets” Kepler-62e and Kepler-62f orbiting the star Kepler-62.


Introduction
The classical restricted three-body problem (RTBP) constitutes one of the most important problems in dynamical astronomy.The study of this problem is of great theoretical, practical, historical, and educational relevance.The investigation of this problem in its several versions has been the focus of continuous and intense research activity for more than two hundred years.The study of this problem in its many variants has had important implications in several scientific fields including, among others, celestial mechanics, galactic dynamics, chaos theory, and molecular physics.The RTBP is still a stimulating and active research field that has been receiving considerable attention of scientists and astronomers because of its applications in dynamics of the solar and stellar systems, lunar theory, and artificial satellites.
A different kind of restricted three-body problem was formulated by Robe [1], a set up in which the first primary is a rigid spherical shell filled with homogenous, incompressible fluid of density  1 , and the second primary is a mass point outside the shell and moving around the first primary in a Keplerian orbit, while the infinitesimal mass is a small sphere of density  3 moving inside the shell and is subject to the attraction of the second primary and the buoyancy force due to the fluid.
In estimating buoyancy force, Robe [1] assumed that the pressure field of the fluid  1 has spherical symmetry around the center of the shell, and he considered only one out of the three components of the pressure field, which is due to the own gravitational field of the fluid  1 .
A. R. Plastino and A. Plastino [2] took into account all these components of pressure field.But in their study, they assumed the hydrostatic equilibrium figure of the first primary as Roche's ellipsoid.They found that when the density parameter  is zero, every point inside the fluid is an equilibrium point; otherwise, the ellipsoid's center is the only equilibrium point.They also examined the linear stability of equilibrium points.Hallan and Rana [3] investigated the existence of all equilibrium point and their stability in the Robe's [1] problem.It was seen that the Robe's elliptic restricted three-body problem has only one equilibrium point for all values of the density parameter and the mass parameter, while the Robe's circular restricted three-body problem can have two, three, or infinite numbers of equilibrium points.As regards to the stability of these equilibria, they found that the collinear equilibrium points are stable while triangular and circular points are always unstable.Recently, Kaur and Aggarwal [4] investigated the Robe's problem of 2 + 2 bodies and applied it to the study of the motion of two submarines in the Earth-Moon system.Singh and Hafsah [5] examined the Robe's circular restricted three-body problem when the first primary is a fluid in the shape of an oblate spheroid and the second primary is a triaxial rigid body.
The classical restricted three-body problem assumes that the masses of celestial bodies are constant.However, the phenomenon of isotropic radiation or absorption in stars led scientists to formulate the restricted problem of three bodies with variable mass.As an example, we could mention the motion of rockets, black holes formation, motion of a satellite around a radiating star surrounded by a cloud and varying its mass due to particles of the cloud, and comets loosing part or all of their mass as a result of roaming around the Sun (or other stars) due to their interaction with the solar wind which blows off particles from their surfaces.The problem of the motion of astronomical objects with variable mass has many interesting applications in stellar, galactic, and planetary dynamics.
The study of two bodies with variable masses seems to have been first investigated by Dufour [6] where he examined the astronomical phenomena of variable mass relating the secular variation of lunar acceleration with the increase of the Earth's mass due to the impact of meteorites.Later, Gylden [7] wrote the differential equations of motion for the problem when the masses are subject to variation.The integrable case to this differential equation was then given by Meshcherskii [8,9].The problem was later known as the Gylden-Meshcherskii problem.A characterization of this problem was studied by Singh and Leke [10].The effect of the isotropic variation of the mass of the star in a planetary system and the possible ejection of a planet from the system were studied by Veras et al. [11].Recently, Singh and Leke [12] investigated the existence and stability of equilibrium points in the Robe's restricted three-body problem with variable masses.
Besides the Gylden-Meshcherskii problem, there are other different cases of two bodies with variable masses, which are classified according to the presence or absence of reactive forces, to whether the bodies move in an inertial frame or not, and so on (see [13]).For instance, when the particles are at rest in an inertial coordinate system, this case may be used to study the orbits of a celestial body moving through a static atmosphere, whose particles attach to it or detach from it as it moves.The restricted three-body problem with nonisotropic variation of the masses has been studied by Bekov [14], Bekov et al. [15], and Letelier and Da Silva [16].A simple example of this kind of problem is the system of two variable primaries and a rocket.In this case, it is the thrust from the rocket that defines the force that acts on the test particle, in addition to the gravitational attraction from the two primaries, while the rocket does not affect the orbits of the primaries.
In this paper, the existence and stability of equilibrium points under the frame of the Robe problem [1], when the participating bodies vary their masses at the same rate, is studied.Here, we assume that the primaries move in a stationary medium, from which they absorb or lose mass; the first primary being a fluid in the shape of a sphere and the test particle which is a small sphere located inside the fluid also gain or lose mass to the fluid.Hence, there is no need to assume a rigid spherical shell.This study may be useful in the investigations of dynamic the problem of water-planetary system discovered by Kepler spacecraft.These "ocean planets" are orbiting the star Kepler-62 and are designated Kepler-62e and Kepler-62f.The existence of these Earth-size planets covered completely by a water envelope (water planets) has long fascinated scientists and the general public.The model of this problem can also be used to study the small oscillation of the Earth's inner core taking into account the Moon's attraction during the course of evolution.
This paper is orginzed as follows: Section 2 contains the equations of motion; the equilibrium points are investigated in Section 3; Section 4 investigates the linear stability of the equilibrium points; Section 5 discusses the obtained results and the conclusions.

Equations of Motion
Let  1 be the mass of the first primary which is a fluid in the shape of a sphere of radius R with center at  1 having density  1 and volume  1 .Also, let  2 be the mass of the second primary with center at  2 which describes a circular orbit around the first one.Both masses are assumed to vary with time as they travel in a static medium which acts as a sink or source of mass.Now, let  3 be the mass of the test particle whose mass is very small compared with the masses of the primaries, with center at  3 , having density  3 and volume  3 .We suppose that its mass varies with time, also as it moves about in the fluid, it gains or loses mass to the medium.Let the positions vector between the center of the fluid and the centers of the second and the test particle be ⃗  12 and ⃗  13 , respectively, and let that between the test particle and the second primary be ⃗  23 .Following Robe [1] and knowing that the masses, distances, and densities vary with time, the forces acting on the third body are the force of attraction of  2 ; the gravitational force ⃗   exerted by the fluid, that is ⃗   = −(4R 3 /3   →  3  1 .We adopt a rotating coordinate system  with origin at the center of mass, , of the primaries,  pointing towards the second primary, and  being the orbital plane of  2 .The equations of motion of the test particle, taking into account the forces acting on it, have the following form [1,14]: where  2 13 = ( −  1 ) 2 +  2 +  2 ,  2 23 = ( −  2 ) 2 +  2 +  2 .The barycentric coordinates  1 and  2 are connected with the distance between the primaries by the following equations: where  1 () =  1 (),  2 () =  2 (),  1 () =  1  1 ;  1 = 4R 3 /3, while  is the gravitational constant and the over dot denotes differentiation with respect to time .Now, in order to obtain useful dynamical predictions, we transform (, , , ) to the autonomous form (, , , ).Following [15], the time dependence of the masses is described [15] by the function (): As in [14], the particular solutions for the case with variable parameters  1 ,  2 , and  3 in the form of the Eddington-Jeans laws with indices  = 3 and  = 6 is expressed with the help of the function  : where ,  1 , and  2 are constants.The exponent  = 1 falls in the stellar range while  = 2 and  = 3 result, respectively, in the first and second law of Meshcherskii [8,9] mass variations.Equations of (4) indicate that the laws of variation of the three masses are the same.Also, the dynamical system has the particular solution of the following type [15,17]: where  > 0 is a constant, and  ̸ = 0 is a constant of the area integral.
Finally, in addition, we assume that the densities of the fluid and the test particle vary such that where  01 and  03 are the densities of the medium and the test particle, respectively, at initial time  0 .Substituting (3) to ( 6) in ( 1) and reducing it throughout by  7 / 7 0 , we get where Here,  12 is constant and connects the parameter  by the relation Equations ( 5) and ( 9) indicate that the ratio of the product of the distances between the center of the primaries, mass of the test particle and the sum of the masses with the gravitational constant to the constant of the area integral, always remains a constant in both the autonomous and the non autonomous systems.Now, we choose units for the distance and time, such that at initial time  0 ,  12 = 1,  0 =  = 1, respectively.Putting these in (9), for the unit of sum of the masses, we get  0 =  = .
Next, without loss of generality, we introduce the mass parameter defined as  =  20 / 0 (0 <  < 1) and also assume that the pressure field of the fluid of density  01 maintains a spherical symmetry around the center of the fluid such that  13 = R. Also, we have  1 = − and  2 = 1 − .With the help of these units, the system of (7) takes the following form: where and the dash signifies differentiation with respect to the new time .Equations ( 10) are the autonomized equations of motion of the test particle of our problem.These equations are different from that in [18] and analogous to the equations in Hallan and Rana [3] only differing due to the second term that appears in the force function and the parameter .

The Equilibrium Points
The equilibrium points represent stationary solutions of the RTBP.These solutions are the singularities of the manifold of the components of the velocity and the coordinates and are found by setting   =   =   =   =   =   = 0 in the equations of motion (10).That is, they are the solutions of the equations Ω  = Ω  = Ω  = 0, which are The solutions are categorized as follows.
(1) The solutions of first equation of ( 12) with  =  = 0 yield the collinear equilibrium points.These points lie on the line joining the center of the first and second primary.
(2) When the first and the second equation of ( 12) are solved with  = 0, we get the circular points.These points lie in the spherical fluid and form circles.
(3) The solutions of first and the third equations of (12) with  = 0 results in the out-of-plane equilibrium points.These solutions are valid provided they lie inside the fluid.We shall consider them in Sections 3.1, 3.2, and 3.3, respectively.

Collinear Points.
From first equation of ( 12), with  =  = 0, we get Hence, (13) has three roots, the first being  = − and is always a solution whether  = 0 or is not, for 0 <  < ∞.The two remaining roots are found by considering the second equality in (13), which gives where Δ = (4 +  − 4).Now, the solutions exist only for  > 1, however, the second solution is greater than 1 −  and consequently will lie outside the shell, so we ignore it.The third solution is and is less than 1 − , the -coordinate of the second primary.Hence, there are two collinear equilibrium points which lie on the line joining the centers of the primaries.Therefore, the coordinate (−, 0, 0) is always an equilibrium point.For  > 1, there exist an equilibrium point ( 3 , 0, 0) which lies to the left or right of the first primary depending upon whether  < 1 + 2 or  > 1 + 2.When  = 1 + 2, the only equilibrium point is the center of the shell.

Circular
Points.These solutions are found by solving the first and second equations of ( 12) with  = 0. Solving the second equation of (12) gives Substituting (16) in the first equation of ( 12) yields  = 1 −  and consequently  3 23 = 1.Hence, we have the solution which gives the coordinate of any point on the circle (17) with center (1 − , 0, 0) which is the center of the second primary and radius one which is the distance between the centers of the fluid and the second primary.Thus, the solution gives us an infinite number of equilibrium points, provided they lie inside the fluid.

Out-of-Plane
Points.The out-of-plane equilibrium points are found by solving first and third equations of (12) with  = 0. Solving first for  23 in the third equation of ( 12), we get Substituting (18) in the first equation of ( 12) and simplifying results in Knowing that  3 23 = ( +  − 1) 2 +  2 , substituting ( 18) and (19) in it, and solving for , we get Equations ( 19) and (20) give the position (, 0, ) of real outof-plane equilibrium points provided  ≥ 1/(1 − ).Should  < 1/(1 − ), then no real out-of-plane points exist as (20) turns out to be imaginary or complex quantity.When  = 1, these points fully coincide with that of Hallan and Rana [3].When  = 0, (i.e.,  01 =  03 ), the coordinates of the out-ofplane points become the same with that in Singh and Leke [12].Hence, it is seen that the equilibrium points are fully analogous to those found by Hallan and Rana [3] except for the out-of-plane equilibrium points which are affected by the parameter , the density parameter, and the mass ratio.The positions of the out-of-plane points are given in Tables 1, 2 and 3 and their graphical representations in Figures 1, 2, and 3 for a test particle in the Earth-Moon system, when the density parameter  is negative, zero, and positive, respectively.We summarize our numerical effort as follows.In Table 1, when  = −0.001this implies that  03 <  01 .In this case for 0 <  ≤ 0.999 and 1.009 <  < ∞, out-of-plane points do not exist, but however exist in the interval 0.999001 ≤  ≤ 1.999.In Table 2,  = 0 and so  03 =  01 .In this case the out-of-plane points exist only when 1.001 ≤  ≤ 1.01 and do not exist in the remaining entire range of , while for Table 3,  = 0.001 and so  03 >  01 .Here, real out-of-plane solutions exist, for the values of the parameter , in the interval 1.002 ≤  ≤ 1.0111 and are nonexistent for any value of  outside this range.Hence, though  has a large range of values, however, the physically meaningful range is  ∈ [0.999001, 1.0111], at which the out-of-plane points exist (see Figure 4).However, these ranges may differ for different density parameters , which are determined by the densities of the fluid and the test particle.

Stability of Equilibrium Points
To examine the stability of an equilibrium configuration, that is, its ability to restrain the body motion in its vicinity, we apply small displacement , V,  to the coordinates ( 0 ,  0 ,  0 ) of the third body, to the positions,  =  0 + ,  =  0 + V and  =  0 + .If its motion rapidly departs from the vicinity of the point, we call such a position of equilibrium an unstable one.If however the body merely oscillates about the point, it is said to be a stable position.Now, we linearize (10) to obtain the variational equations: where the superscript 0 indicates that the partial derivatives are to be evaluated at the equilibrium points.

Collinear Points.
Robe [1] discussed the stability of the equilibrium point at the center of the shell, and Hallan and Rana [3] have also discussed that in the case of the noncollinear points when  < 0. Hence, we shall discuss here the stability of the equilibrium point ( 3 , 0, 0) near the center of the fluid.To do this, we let solutions of the first two equations of (21) be  =  exp(), V =  exp(), where , , and  are constants.Finding first and second derivatives of the solutions, substituting them in the first two equations of (21), and simplifying, we obtain the matrix which has a nonzero solution when Expanding the determinant, we get This is the characteristic equation corresponding to the variational equations (21) when motion is considered in the -plane.Now, the values of the second order partial derivatives computed at the point ( 3 , 0, 0), with the substitution  =  = 0 are as follow: where  1 = (4( − 1) 3 )/[ + √(4 +  − 4)] 3 .
When  = 1, the equations in system (24) fully coincide with those of Hallan and Rana [3].
The roots of (28) are where Δ is the discriminant of ( 27) and is always positive for any .
When  = 1, the value of Δ fully coincides with that of Hallan and Rana [3] which ought to be Δ = 4{4( −  1 ) + 9 2 1 }.Now, since Δ > 0, if in (28) the quantity in parenthesis is positive; that is,  > 0, which occurs when we see that  > √ Δ. Hence from (28), we have − ± √ Δ < 0, and so both values of Λ  are negative and consequently, the roots (28) are distinct pure imaginary.Therefore, the equilibrium point ( 3 , 0, 0) is stable provided equation (30) holds; otherwise, it is unstable.The partial derivatives at these points are When these are substituted in the variational equations ( 21) with  = 0, we have The last equation of (32) shows that motion is stable along the -axis.The characteristic equation of the first two equations of (32) is Its roots are  1,2 = 0,  3,4 = ± √ −4 + 3, and so the equilibrium points are unstable due to multiple zero roots.

Out-of-Plane Points.
For the stability of the out-of-plane equilibrium points, we consider the following partial derivatives: The characteristic equation in this case is gotten by substituting the trial solutions  = exp(), V = exp(),  =  exp() in the variational equations ( 21) to get where Substituting (34) in ( 35), we at once have Its roots are where The roots (37) are computed numerically for motion of a test particle (a submarine) under the gravitational attraction of the Earth-Moon system when  > 0 and  < 0 ( 01 <  03 and  01 >  03 , resp.) for 0 <  < ∞.Therefore, we take  = 0.01 and  ∈ [0.999001, 1.0111] with the following values in each case following Kaur and Aggarwal [4].In the case when  01 <  03 , we take Aside from these examples, we also consider the case when  = ±0.001(see Tables 1 and 3), so that a wider generalization can be reached regarding the characteristic roots (38) which consequently determines whether the equilibrium point is a stable one or not.
Using the software package Mathematica; the six characteristic roots are presented in Tables 4, 5, 6, and 7 numerically for  = 0.01, for different density parameters and a wide range of the parameter .We seek to find the case where all the six roots are pure imaginary quantities or complex figures with negative real parts.If this happens, then the solutions will be bounded and motion will be stable; otherwise, they will be unstable.
From these tables, we see that that for a specific set of values of these parameters at least one of the roots among all has a positive real part or a complex root with the existence of a positive real part.Therefore, this causes the solutions to be unbounded and consequently producing unstable equilibrium points.Hence, we conclude that the out-of-plane equilibrium points are unstable equilibrium points due to a positive root and positive real part in complex roots.This agrees with the result of Singh [5].
The equilibrium solutions of the nonautonomous system with variable coefficients are in general unstable points according to the Lyapunov's theorem of stable solutions [19].

Discussion and Conclusion
We have derived the equations of motion and established the possible equilibrium points of the third body of infinitesimal mass in a setup of Robe's [1] restricted three-body problem when the three participating bodies all vary their masses arbitrarily at the same rate and the density of the fluid and third body also vary as the masses.We find that the nonautonomous equations of motion are different from that of the restricted problem of three variable mass bodies derived by Bekov et al. [15], while the autonomized equations of our study are also different from those of Robe [1] and Hallan and Rana [3] due to the assumptions we have introduced here.The autonomized system of ( 6) is different from that in Hallan and Rana [3] only due to an additional term which appears due to the parameter .When  = 1, the equations become fully analogous to theirs.
The equilibrium points are sought, and it is seen that the point at the center of the fluid is always an equilibrium point of the Robe problem.An equilibrium point near the center of the fluid, points on the circle (circular points), and two out-of-plane points on the -plane also exist, with all similar to that found in Hallan and Rana [3]; though the later have several points in our case while there exists only a pair in that of Hallan and Rana [3].
The linear stability of the equilibrium points of the autonomized have been studied and the outcomes are analogous with the stability results in Hallan and Rana [3], in the sense that the equilibrium points collinear with the centers of the fluid and the second primary are stable under some given conditions which depend on the mass ratio, density parameter, and the parameter .The circular are unstable due to the presence of multiple zero roots, while the out-of-plane equilibrium points are also unstable due to a positive root and a positive real part of the complex roots.
In our recent paper, Singh and Leke [12], the motion of a test particle around the equilibrium points was generalized to include the effect of mass variations of the primaries which vary isotropically in accordance with the unified Meshcherskii law, when the motion of the primaries is determined by the Gylden-Meshcherskii problem.Here, ejection or attachment form or to the surfaces of the primaries do not create reactive forces.Also, we have taken the first primary as   a rigid spherical shell filled with a fluid of constant density and volume and containing the test particle, while in the present study, we have assumed that the first primary is a fluid in the shape of a sphere with nonisotropic mass and density variation, given that, the test particle is contained in the fluid.The second primary and the test particle both have their masses varying arbitrarily with time at the same rate as the first primary.
In the previous study, the autonomized dynamical system with constant coefficients is gotten, only when the shell is empty or when the densities of the medium and the test particle are equal, while in the present study, such limitation do not arise.In the present study, we found two collinear equilibrium points on the line joining the centers of the fluid and the second primary with one at the center of the fluid and the other away from it.Further, circular equilibrium points exist on the -plane and pairs of out-of-plane points which depend on the arbitrary constant  ∈ (0, ∞), and density and mass parameters are found on -plane.However, in the previous study, there is only one collinear equilibrium point located at the center of the rigid shell and a pair of the out-of plane equilibrium points which exist only for  > 1.
The linear stability analysis however turns out to be same as the equilibrium points on the line collinear with the centers of the primary of the autonomized system which are conditionally stable; while the equilibrium points on the plane and the circular points are unstable.
The result of our research work can be summarized as follows.The restricted problem under the framework of the Robe's [1] problem with three variable mass bodies, which vary arbitrarily with time at the same rate, has the equilibrium points which; are, the points near the center of the fluid, points on the circle (circular points), and pairs of out-ofplane points.These equilibrium points are analogous to the problem studied by Hallan and Rana [3] and their stability results are similar, though we have assumed that the first primary is a spherical fluid and that the masses of the test particle and primaries vary with time.The idea of the submarine's mass changing with time may perhaps seem impossible because it is a rigid body.However, from Archimedes' principle, a submarine floating or submerging depends on buoyancy which is controlled by the ballast tanks found between the submarine's inner and outer hulls.Submarine resting on the water surface has positive buoyancy which means it is less dense than water.When this happens, the ballast tanks are empty.For it to submerge, vents on top of the tanks are opened and water floods in thereby making the submarine denser than the sea water.In this case, the submarine has negative buoyancy.For it to float again, the water in the tanks will be forced out and the submarine becomes less dense than the water and eventually floats.In this study, the buoyancy force depends on the mass of the submarine  3 , the densities of the submarine, and the medium, as well as the volume of the submarine.And so by the simple relation, connecting mass, density, and volume, we understand that, when water is allowed into the ballast, as the mass of the submarine is changing so is the density until it becomes denser than the medium and sinks.When water is let out of the tanks, the mass and density of the submarine reduce and as a matter of fact, it will float when the medium is denser.Neutral buoyancy is attained when the weight of the submarine equals the amount of water it displaces.The submarine will neither rise nor sink in this state.Hence, with changing mass comes a changing density while the volume remains a constant.
This study may be useful in the investigations of the dynamic problem of Earth-size planets covered completely by a water envelope (water planets), and also the study of the small oscillation of the Earth's inner core taking into account the Moon's attraction during the course of evolution.The problem discussed in this paper is highly idealized and therefore calls for more research.