Introduction and results
The well-known Newtonian n-body problem concerns with the motion of n mass points with positive mass mi moving under their mutual attraction in Rd in accordance with Newton’s law of gravitation.
The equations of the motion of the n-body problem are
where we have taken the unit of time in such a way that the Newtonian gravitational constant be one, and ri€Rd (i=1…,n) denotes the position vector of the i-body,
is the Euclidean distance between the i-body and the j-body.
The solutions of the 2-body problem (also called the Kepler problem) has been completely solved, but the solutions for the n-body for n>2, is still an open problem.
For the Newtonian n-body problem the simplest possible motions are such that the configuration formed by the n-bodies is constant up to rotations and scaling, such motions are called the homographic solutions of the n-body problem, and are the unique known explicit solutions of the n-body problem when n>2. Only some special configurations of particles are allowed in the homographic solutions of the n-body problem, called by Wintner [1] central configurations. Also, central configurations are of utmost importance when studying bifurcations of the hypersurfaces of constant energy and angular momentum, for more details see Meyer [2] and Smale [3]. These last years some central configurations have been used for different missions of the spacecrafts in the solar system, see for instance [4,5].
More precisely, let
be the total mass and the center of masses of the n bodies, respectively.
A configuration r=(r1,…,rn) is called a central configuration if the acceleration vectors of the n bodies are proportional to their positions with respect to the center of masses with the same constant λ of proportionality, i.e.
where λ is the constant of proportionality.
Equations (1) are strongly nonlinear and to find the explicit central configurations (r1,…,rn) in function of the masses m1,….mn when n>3 is an unsolved problem.
There is an extensive literature on the study of central configurations, see for instance Euler [6], Lagrange [7], Hagihara [8], Llibre [9,10], Meyer [2], Moeckel [11], Moulton [11], Saari [12], Smale [3], ..., and the papers quoted in these references.
In this paper we are interested in the planar central configurations of a circular restricted 4-body problem. Of course, for the central configurations of the 4-body problem there are many partial results, see for instance the papers [13-66].
We note that the set of central configurations is invariant under translations, rotations, and homothecies with respect their center of mass. It is said that two central configurations are equivalent if after having the same center of mass (doing a translation if necessary) we can pass from one to the other through a rotation around its common center of mass and a homothecy. This defines a relation of equivalence in the set of central configurations. From now on when we talk about a central configuration, we are talking on a class of central configurations under this relation of equivalence.
The objective of the present article is to study the central configurations of the circular restricted 4-body problem with three equal primaries in the collinear central configuration of the 3-body problem. We recall that for the 3-body problem when the three masses are equal there is a unique collinear central configuration, where the mass in the middle equidistant from the other two, of course the equal masses can be permuted in the positions of this configuration.
As in any circular restricted problem the objective is to describe the motion of the infinitesimal mass with respect to the primaries. Usually this problem is studied in a rotating system of coordinates where the positions of the primaries remain fixed, see for more details on the restricted problems the book of Szebehely [67].
More precisely, taking the unit of mass equal to the masses of the three primaries and since a central configuration is invarinat under rotations and homothecies through its center of mass without loss of generality we can assume that the position vector rj of the three primaries with masses m1= m2= m3= 1 are
We denote the position of the infinitesimal mass m4 = 0 by r4 = (x4,y4) = (x,y). Then our main result is the following one.
Theorem 1 The circular restricted 4-body problem with three primaries of equal masses m1= m2= m3= 1 with position vectors given in (2), and one infinitessimal mass m4 = 0 with position vector r4 =(x4,y4) =(x,y) have the following six central configurations with r4 = Pj for j=1 being:
(i)
, where the value of the coordinate is a root of the polynomial
;
(ii)
;
(iii)
, where the value of the coordinate is a root of the polynomial
;
(iV)
, where the value of the coordinate is a root of the polynomial
;
(v)
;
(vi)
.
Figure 1.
Figure 1: The six central configurations of the circular restricted 4-body probl qem with three equal primaries in the collinear central configuration of the 3-body problem. The three primaries are indicated with the big circles, and the position of the infinitesimal mass in the corresponding six central configurations is indicated with a small circle.
The proof of Theorem 1 in given in the next section.
Proof of theorem 1
From (1) we obtain the following eight equations for the central configurations of the 4-body problem in the plane
(3)
Where c= (c1, c2). Substituting in (3) the expressions (2), m1= m2= m3= 1, m4 = 0 and r4 =(x4,y4) =(x,y), corresponding to our circular restricted 4-body problem these eight equations reduce to
Therefore
, and the position vector of r4 =(x4,y4) in order to have a central configuration of the circular restricted 4-body problem must be a real solution of the system
In Figure 2 we have shwon the curves e4 (x,y) =0 and e8 (x,y) =0, and in Figure 3 the intersection of these two curves. We see that these two curves intersect in six points inside the rectangle
. Computing the coordenates of these six points numerically using the Newton method (see for instance [68]), we get the six points Pj which appear in the statement of Theorem 1. Of course we have omitted the three points where are located the three primaries in the intersections of the two curves e4 (x,y) =0 and e8 (x,y) =0, because there really these two curves are not defined. Now we shall prove that these six points obtained numerically really are solutions of the system e4 (x,y) =0 and e8 (x,y) =0.
Figure 2: The curves e4(x,y) = 0 and e8(x,y) = 0 in the rectangle (x,y) ∈ [−2.2,2.2] × [−2.2,2.2].
We note that equations e4 (x,y) =0 and e8 (x,y) =0 are invariant if we change x by -x, and y by -y, so if (x,y) is a solution of the system e4 (x,y) =0 and e8 (x,y) =0, also (-x,y), (x, -y) and (-x,-y) are solutions. So in order to prove Theorem 1 we only need to study the solutions of system e4 (x,y) =0 and e8 (x,y) =0 satisfying
and
. Moreover, from Figure 3 we see that all the solutions are of the form (x,0) or (0,y), and since in the origin (0,0) there is one primary, we must look only for the solutions (x,0) or (0,y) with x > 0 and y> 0.
Figure 3: In this picture we can see the six intersection points between the two curves e4(x,y) = 0 and e8(x,y) = 0 different from the positions of the primaries, which provide the six central configurations of the circular restricted 4-body problem with three primaries of equal masses at the collinear configuration of the 3body problem and an infinitisimal mass.
First we look for the solutions (0,y) with y>0, then system e4 (x,y) =0 and e8 (x,y) =0 reduce to
or equivalently to
Squaring the both sides of the this equation we get the equation
This polynomial equation has only two real roots
but only the second root satisfies equation (4). This provides the solution P1 of Theorem 1, and consequently also the solution P2.
Now we look for the solutions (x,0) with (x>0) of the system e4 (x,y) =0 and e8 (x,y) =0. For these solutions the system reduce to
squaring the both sides of the previous equality we obtain
Writting this equation with a common denominator, which only vanishes at the positions of the primeries, its numerator equal zero can be written as
Squaring again the both sides of the this equality we get
The real zero x = 1 is not good because it correspond to the position of a primary. The unique real root of the polynomial -4+5x3-12x4-10x5+5x7 is 1.7576799791694022.. which also is a zero of equation (5), and consequently provides the central configuration P3, and by the symmetries of the equations of the central configurations also provides the central configuration P6.
The unique real root of the polynomial -4+8x2+21x3-4x4-10x5+5x7 is 0.4946664910171345… which also is a zero of equation (5), and consequently provides the central configuration P4, and due to the symmetries of the equations of the central configurations also provides the central configuration P5.
The real roots of the polynomials -4+8x2-11x3-4x4-10x5+5x7 and 4+16x2+5x3+4x4-10x5+5x7 are not zeros of the equation (5). This completes the proof of Theorem 1.