# Stability Analysis of the Rhomboidal Restricted Six-Body Problem

Authors: M. A. R. Siddique; A. R. Kashif; M. Shoaib; S. Hussain
Journal: Advances in Astronomy
 (2021)

---

## Abstract

We discuss therestricted rhomboidal six-body problem (RR6BP), which has four positive masses at the vertices of the rhombus, and the fifth mass is at the intersection of the two diagonals. These masses always move in rhomboidal CC with diagonals 2a and 2b. The sixth body, having a very small mass, does not influence the motion of the five masses, also called primaries. The masses of the primaries are m1=m2=m0=m and m3=m4=m˜. The masses m and m˜ are written as functions of parameters a and b such that they always form a rhomboidal central configuration. The evolution of zero velocity curves is discussed for fixed values of positive masses. Using the first integral of motion, we derive the region of possible motion of test particle m5 and identify the value of Jacobian constant C for different energy intervals at which these regions become disconnected. Using semianalytical techniques, we show the existence and uniqueness of equilibrium solutions on the axes and off the axes. We show that, for b∈1/3,1.1394282249562009, there always exist 12 equilibrium points. We also show that all 12 equilibrium points are unstable.


---

## Body

## 1. Introduction

In therestricted n-body problem, a body of infinitesimal mass moves under the gravitational influence of n−1 massive bodies called primaries. The force exerted by the infinitesimal mass on the primaries can be ignored as it has a negligible mass, whereas primaries revolve perpetually in concentric circles around their common center of mass by retaining a particular central configurations (CCs). The CC play a vital role in the understanding of the n-body problem in celestial mechanics. It can be used to find simple or special solutions of the n-body problem since the geometry formed by the arrangement of the primaries remains constant for all time.The CC forn>3 has been investigated extensively after being identified as the problem of the 21st century by Arnol’d [1]. Xia [2] used the method of analytical continuation to find exact number of central configurations for n positive masses. Hampton and Moeckel [3] studied the finiteness of relative equilibria of the four-body problem and showed that the number of relative equilibria is always between 32 and 8472 for the Newtonian four-body problem. The homographic solutions to the rhombus four-body problem are the variational minimizers of the Lagrangian action confined on a holonomically constrained rhombus loop space, according to Mansur and Offin [4]. The approach used by Mansur et al. [5] to demonstrate spectral instability for the entire parameter ranges of mass ratio and eccentricity e≥0 for the homographic family of rhombus solutions within the four degree of freedom parallelogram four-body problem is based on a topological invariant called the Maslov index. Shi and Xie [6], using analytical approach, have shown that there is exactly one family of concave and one family of convex central configurations in addition to the family of equilateral triangle configurations. Llibre and Mello [7] classified the central configurations of the four-body problem. Liu and Zhou [8] investigated the four-body problem with three masses forming a Lagrangian triangle, used the bifurcation diagram of linearly stable and unstable regions, and found two linearly stable subregions with respect to α, β, and e. Deng et al. [9] investigated the CC of the four-body problem with equal masses and showed that, for the planar Newtonian four-body problem having adjacent equal masses i.e., m1=m2≠m3=m4 and equal lengths for the two diagonals, any convex noncollinear CC must have a symmetry and must be an isosceles trapezoid. They also showed that when the length between m1 and m4 equals the length between m2 and m3, the CC is also an isosceles trapezoid. Marchesin [10] investigated the rhomboidal configuration stability with a central mass and two pairs of equal masses. The mass m0 is at the center of the configuration, and the equilibria obtained in this case were all shown to be unstable. Shoaib et al. [11] established the central configuration for the rhomboidal 5-body problem and highlighted the regions in the phase plane where it is possible to have central configuration. On the axis of symmetry, Shoaib et al. [12] considered a symmetric five-body problem with three unequal collinear masses. The remaining two masses were symmetrically arranged on both sides of the axis of symmetry, and areas of feasible central configurations were identified analytically and numerically for the rhomboidal and triangular four- and five-bodies. Using Levi-Civita type transformations, the equations of motion were regularized, and the phase space for chaotic and periodic orbits was explored using the Poincaré surface of sections. Zotos [13] numerically explored the restricted four-body problem with three equal masses with a dynamically stable triangular configuration and found that the linearly stable Lagrange points only exist when one of the three masses has a considerably larger mass. Dewangan et al. [14] investigated the elliptic restricted four-body problem by considering radiation and oblateness effects; they considered a bigger primary as a radiation source and the other primaries to be of equal masses as oblate spheroid. They found that the equilibrium points to be linearly stable. Liu et al. [15] studied the four-body problem and found that the boundaries of possible motions obey the change in parameter c2E, that is, if the value of c2E is less than or equal to a critical value c2Ecr, then the system is stable. Ismail et al. [16] studied the four-body problem by considering the effects of radiation pressure and oblateness and used the Lyapunov function to show the stability of equilibrium points. Wang and Gao [17] did a numerical study of the restricted five-body problem regarding the zero velocity surface and transfer trajectory by considering four equal masses (primaries) forming a regular tetrahedron configuration and the fifth (infinitesimal) mass moving under the gravitational influence of the four primaries. They numerically simulated the zero velocity surface of the infinitesimal mass in the three-dimensional space and designed the transfer trajectory of the infinitesimal mass. Suraj et al. [18] studied the five-body problem to investigate the effects of perturbation parameter on the positions, motion, and stability of the libration points due to the variable mass of the fifth mass. Li and Liao [19] obtained 695 families of Newtonian periodic planar collisionless orbits of three-body systems with equal mass and zero angular momentum numerically. In the planar restricted three-body problem, Sosnitskii [20] investigated Lagrange stability and proved a theorem on the Lagrange stability of the infinitesimal particle particularly for the circular restricted three-body problem. Ding et al. [21] proved that there exists another Eulerian collinear central configuration at any instant in the three-body problem. Using the port-Hamiltonian approach, Liu and Dong [22] reformulated the Circular Restricted Three-Body Problem (CRTBP) and obtained the closed loop Hamiltonian by designing a control strategy (based on energy shaping and dissipation injection) as a candidate of the Lyapunov function that assures asymptotic stability. This control technique also demonstrated global stability within the CRTBP model’s application region. Corbera et al. [23] established that every four-body convex center configuration with perpendicular diagonals will have a kite configuration. Lara and Bengochea [24] investigated the symmetric periodic orbits of the four-body system both theoretically and numerically. Alvarez-Ramírez and Medina [25] studied the planar restricted five-body problem in which the four primaries form a axisymmetric four-body central configuration and described the equilibrium points that depend on the mass parameters of the primaries.In the six-body problem, Mello et al. [26] demonstrated the existence of three new families of stacked spatial central configurations. Alsaedi et al. [27] used variational methods and computational algorithm to investigate the six-body problem and its new style periodic solutions. Idrisia and Ullah [28] studied the CC of the restricted six-body problem with the central body and showed that all the libration (equilibrium) points exist on the concentric circles C1, C2, and C3 having center at the origin. The libration points that lie on circles C1 and C3 are unstable, while there are some stable libration points on circle C2.In the present paper, we consider a restricted six-body problem, where four of the primaries are at the vertices of a rhombus and the fifth mass is at the intersection of the two diagonals. Three of the primaries have equal masses ofm1=m2=m0=m and are located on the horizontal axis, and two other equal mass primaries are located on the vertical axis with masses m3=m4=m˜. Before dividing the equations of motion for the restricted 6-body in Section 3, we find continuous families of central configurations for the rhomboidal 5-body problem. The rest of the paper is organized as follows: we investigate the Hill region and possible region of motion of m5 according to the Jacobi constant in Section 4. In Sections 5 and 6, we show the existence, uniqueness, and stability of equilibrium points, respectively. Conclusions are given in Section 7.


## 2. Rhomboidal Central Configurations

In this section, we prove the existence and uniqueness of central configuration of a rhomboidal 5-body problem for positive masses. The mass ratio is written as a function of “a” and “b” which can be used to find regions of central configuration for the rhomboidal 5-body problem. The classical equation of motion for then-body problem has the following form:(1)mir..i=∑j=0,j≠inmimjrj−rirj−ri3,where the units are chosen so that the gravitational constant is equal to one. A central configuration is a particular configuration of the n bodies where the acceleration vector of each body is proportional to its position vector, and the constant of proportionality is the same for the n bodies. Therefore, a CC is a configuration that satisfies the following equation:(2)−ω2ri−c=∑j=0,j≠inmimjrj−rirj−ri3,where ω is angular speed and (3)c=∑i=1nmiri∑i=1nmi,which is the center of mass of the n bodies. We take the position vector (see Figure 1) of the five primaries mj, where j=0,1,…,4 as(4)r0=0,r1=aeiωt,r2=−aeiωt,r3=bieiωt,r4=−bieiωt.Figure 1
The restricted six-body problem.Fori=1and3 in equation (2), the CC equation for masses m1 and m3 is(5)−ω2a,0=m−a,0a3+m−2a,08a3+m˜−a,ba2+b23/2+m˜−a,−ba2+b23/2,(6)−ω20,b=m0,−bb3+ma,−ba2+b23/2+m−a,−ba2+b23/2+m0,−2b8b3,where m1=m2=m0=m and m3=m4=m˜. The positions of the primaries are taken from equation (4). Writing equations (5) and (6) in components’ form, we get(7)ω2=5m4a3+2m˜a2+b23/2,(8)ω2=m1b3+2a2+b23/2+m˜4b3.Because of the symmetry of the problem, the equations form2 and m4 are identical to equations (7) and (8). The equation for m0 is identically zero. We have extended the astrophysical model given in [29] by keeping the mass m0 at the center. The results of the model given in [29] can be reproduced by taking m0=0 in our model. Without loss of generality, we take ω=1 and solve equations (7) and (8) simultaneously for ma,b and m˜a,b:(9)ma,b=Nma,bDma,b,m˜a,b=Nm˜a,bDma,b,where(10)Nma,b=4a3a2+b22a2+b23/2−8b3,(11)Nm˜a,b=4a2+b23/2−4a7−8a5b2+5a4b3−4a3b4+10a2b5−8a3b3a2+b2+5b7,(12)Dma,b=−32a7−64a5b2−32a3b4+5a6+15a4b2−64a3b3+15a2b4+5b6a2+b2.Definingτ=m˜/m, then equations (9)–(11) give(13)τa,b=4−5b3/a31+b2/a23/2+8b3/a38b3/a3−1+b2/a23/2.Again, letμ=b/a in equation (13); we get an alternate form of equation (13):(14)τμ=8μ3+4−5μ31+μ23/28μ3−1+μ23/2.Lemma 1.
τμ given by (13) is a continuous and strictly decreasing positive function in the interval 1/3,1.1394282249562009,limμ⟶1/3τμ=∞, and τ1.1394282249562009=0.Proof.
One can easily check the continuity andlimμ⟶1/3τμ=∞ and τ1.1394282249562009=0 directly from equation (14). It means the value of τμ which is the mass ratio will remain positive i.e., τ∈0,∞ as the ratio of the distance parameters varies between 1/3,1.1394282249562009. Now we need to prove that the function is decreasing and positive for the given interval. Taking the derivative of τμ, (15)τ′μ=5κ′μ,κ′μ=15μ2μ2+1μ2+15/2−8μ5+1μ2+13/2−8μ32,where τ′μ is a constant multiple of κ′λ given in [29] by Kulesza et al. The proof will therefore be similar and can be seen in [29].
Here, we present some particular cases for different values ofτ and μ. We show these different shapes of rhombus for different values of τ and μ in Figure 2.(i)
Ifτ=1, then μ=1; from equation (13), we get the shape of true square with m=m˜=1 and a=b=1(ii)
Ifτ=0.67, then μ=1.04232; from equation (13), we get m˜=0.67m and b=1.04232a(iii)
Ifτ=0.97, then μ=1.00367; from equation (13), we get m˜=0.97m and b=1.00367a
There are two cases in which the RR6BP central configuration degenerates.(i)
Ifμ=1.1394282249562009, then τ=0; this implies m˜=0 and b=1.1394282249562009a(ii)
Ifμ=1/3, then τ=∞; this implies m=0 and b=1/3aFigure 2
Some possible shapes for the rhombus for different values ofτ and μ.


## 3. Equations of Motion

In this section, we describe the motion of the infinitesimal body,m5, under the gravitational attraction of the five primaries. We assume that the sixth body has a significantly smaller mass compared to the masses of the primary m5≪m0,m1,m2,m3,m4. On this basis, the sixth body acts as an infinitesimal test particle, and therefore, it does not influence the motion of the five primaries. In the RR6BP, the equations of motion of m5 are(16)r..5=−∑j=04mjr5−rjr5−rj3,where dot represents the derivative with respect to time. From here onward, without loss of generality, we take the value of a=1. The equations of motion of m5 in the synodical coordinates x and y are(17)x¨−2y˙=Ux,(18)y¨+2x˙=Uy,where(19)Ux,y=x2+y22+m1r50+1r51+1r52+m˜1r53+1r54,which is the effective potential. The mutual distances of m5 from the primaries in the corotating frame are(20)r50=x2+y2,r51=x−12+y2,r52=x+12+y2,r53=x2+y−b2,r54=x2+y+b2.The Jacobian constant is given by Curtis [30]:(21)C+U=12x2.+y2.=v2.For a given value of the Jacobi constant,v2 is only a function of position in the rotating frame. Since v2 cannot be negative, it must be true that(22)C+U≥0.The boundaries between forbidden and allowed regions of motion are found by settingv2=0, i.e.,(23)C+U=0.It is now trivial to show thatCx,y is the first integral of motion of system (11) by proving that C˙x,y=0.


## 4. The Hill Regions

The region of permitted motion is also known as the Hill region, and the curves found by equation (16) for various values of C are known as the zero velocity curves.The zero velocity curves whenb=0.67,m=0.07717,andm˜=0.7879 and when b=0.97,m=0.4588,andm˜=0.5766 are given in Figure 4, and the corresponding Hill regions are given in Figure 3. We also give the region of possible motion of m5 for six different values of Jacobi constants C in Figures 5–7 for mass parameters b=0.67,0.97,and1.13. The shaded regions represent the forbidden regions of motion for the infinitesimal mass m5. It is numerically confirmed that the permitted regions are completely disconnected for C≥−2.20,−2.85,and−2.56, as shown in Figures 5–7, for the above values of b. For the increasing values of C, the allowed region of motion (white region) becomes partially connecting at C=−2.574,−3.328,and−3.342 and completely connected at C=−3.2,−3.9,and−3.5. It can be seen from Figures 5–7 that the transition of motion from totally disconnected to completely connected occurs in six stages for b=0.67,0.97,and1.13. For these values of b, m5 can freely move in the gravitational field of the CC region for C≥−2.2,−2.85,and−2.56 and cannot reach any of the primaries for C≤−3.2,−3.9,and−3.5.Figure 3
The effective potential. (a)b=0.67;m=0.0778;m˜=0.7879. (b) b=0.97;m=0.4588;m˜=0.5766.
(a)(b)Figure 4
The evolution of zero velocity curves. (a)b=0.67;m=0.07717;m˜=0.7879. (b) b=0.97;m=0.4588;m˜=0.5766.
(a)(b)Figure 5
The regions of motion ofm5 (white region) for energy interval C∈−3.2,−2.20 when b=0.67 from the top-left corner to the bottom-right corner. (a) C=−3.2. (b) C=−3.157. (c) C=−2.574. (d) C=−2.2832. (e) C=−2.269. (f) C=−2.20.
(a)(b)(c)(d)(e)(f)Figure 6
The regions of motion ofm5 (white region) for energy interval C∈−3.5,−2.85 when b=0.97 from the top-left corner to the bottom-right corner. (a) C=−3.5. (b) C=−3.47. (c) C=−3.328. (d) C=−3.1868. (e) C=−3.121. (f) C=−2.85.
(a)(b)(c)(d)(e)(f)Figure 7
The regions of motion ofm5 (white region) for energy interval C∈−3.9,−2.56 when b=1.13 from the top-left corner to the bottom-right corner. (a) C=−3.9. (b) C=−3.3842. (c) C=−3.342. (d) C=−2.655. (e) C=−2.6360. (f) C=−2.56.
(a)(b)(c)(d)(e)(f)


## 5. Equilibrium Solutions

Equilibrium solutions of the RR6BP are the solutions ofUxx,y=0 and Uyx,y=0. The derivative Ux and Uy of the effective potential, given in equation (19), are found as follows:(24)Ux=x−mxx2+y23/2−mx−ax−a2+y23/2+x+1x+12+y23/2−m˜x1y−b2+x23/2+1y+b2+x23/2,(25)Uy=y−myx2+y23/2−my1x−12+y23/2+1x+12+y23/2−m˜y−by−b2+x23/2+b+yy+b2+x23/2.


### 5.1. Equilibrium Solutions on the Coordinates’ Axes

Since the potential given in equation (19) is invariant under the symmetry x,−y, −x,y, and −x,−y, we will restrict our computation to the first quadrant: x≥0 and y≥0. Initially, we study the existence and number of equilibrium solutions on the axes and then off the coordinate axes. To study the equilibrium solutions on the y-axis, let x=0; then, equations (24) and (25) are given as(26)Ux=0,Uy=y−myy23/2−2my1+y23/2−m˜y−by−b23/2+b+yy+b23/2.To solveUy=0, divide y into subintervals 0<y<b and y>b.


### 5.2.0<y<b; 1/3<b<1.1394282249562009

Rewrite the right hand side of equation (25) by taking into account that y∈0,b:(27)f1y=y−my2−2my1+y23/2+m˜1y−b2−1y+b2.Aty≈0, f1y<0; and, at y≈b, f1y>0; therefore, by the mean value theorem, there is at least one zero of f1y when y∈0,b. The derivative of f1y is given by(28)df1ydy=1+2m1y3+1y2+13/2−3y2y2+15/2+2m˜1y+b3−1y−b3.The only term−3y2/y2+15/2 that can make df1y/dy negative for y∈0,b is dominated by the rest of the term in equation (26); therefore, df1y/dy>0. This proves the existence of the unique equilibrium solution inside rhombus on the y-axis.


### 5.3.y>b

Now, consider the casey>b. Use equation (26) and rewrite as(29)f2y=y−my2−2my1+y23/2−m˜1y−b2+1y+b2.Aty≈b, f2y<0; and, at y≈∞, f2y>0; therefore, by mean value theorem, there is at least one zero of f2y when y∈b,∞. The derivative of f2y is given by(30)df2ydy=1+2m1y3+1y2+13/2−3y2y2+15/2+2m˜1y+b3+1y−b3.Following the same procedure as given for0<y<b, one can easily prove the uniqueness of the equilibrium solution for y>b.We discuss here two special cases of CC forb∈1/3,1.1394282249562009. When b=1/3, then the masses on the horizontal axis are zero, i.e., m=0. In this case, we only get the two equilibrium points along y-axis. When b=1.1394282249562009, then m˜=0. In this case, we get four equilibrium points along x-axis. The positions of the masses and the corresponding equilibrium points for these two cases are shown in Figure 8. The stability of these cases will be discussed in Section 6.Figure 8
(a) Equilibrium points (red color) alongy-axis when m=0. (b) Equilibrium points (red color) along x-axis when m˜=0.
(a)(b)


### 5.4. Equilibrium Solutions off the Coordinates Axes

It is numerically confirmed that, forb=0.67,0.97,and1.13, there are always a total of 12 equilibrium points. As shown in Figure 9, four of the equilibrium points are on the x-axis, four on the y-axis, and remaining four of the equilibrium points are off the axes. Since the gravitational field is a function of mass parameters mb, therefore the equilibrium points change their positions around the primaries for changing values of b. It is numerically confirmed that majority of equilibrium points are around the primaries along the horizontal axis if b is around 1/3 as the masses on the horizontal axis are dominant (see Figures 9(a)). For b≥1, the equilibrium points concentrated around the primaries on the vertical axis (see Figure 9(b)).Figure 9
(a–c) 12 equilibrium points (red colors) for different values ofb=0.67,0.97,and1.13, respectively.
(a)(b)(c)


## 5.1. Equilibrium Solutions on the Coordinates’ Axes

Since the potential given in equation (19) is invariant under the symmetry x,−y, −x,y, and −x,−y, we will restrict our computation to the first quadrant: x≥0 and y≥0. Initially, we study the existence and number of equilibrium solutions on the axes and then off the coordinate axes. To study the equilibrium solutions on the y-axis, let x=0; then, equations (24) and (25) are given as(26)Ux=0,Uy=y−myy23/2−2my1+y23/2−m˜y−by−b23/2+b+yy+b23/2.To solveUy=0, divide y into subintervals 0<y<b and y>b.


## 5.2.0<y<b; 1/3<b<1.1394282249562009

Rewrite the right hand side of equation (25) by taking into account that y∈0,b:(27)f1y=y−my2−2my1+y23/2+m˜1y−b2−1y+b2.Aty≈0, f1y<0; and, at y≈b, f1y>0; therefore, by the mean value theorem, there is at least one zero of f1y when y∈0,b. The derivative of f1y is given by(28)df1ydy=1+2m1y3+1y2+13/2−3y2y2+15/2+2m˜1y+b3−1y−b3.The only term−3y2/y2+15/2 that can make df1y/dy negative for y∈0,b is dominated by the rest of the term in equation (26); therefore, df1y/dy>0. This proves the existence of the unique equilibrium solution inside rhombus on the y-axis.


## 5.3.y>b

Now, consider the casey>b. Use equation (26) and rewrite as(29)f2y=y−my2−2my1+y23/2−m˜1y−b2+1y+b2.Aty≈b, f2y<0; and, at y≈∞, f2y>0; therefore, by mean value theorem, there is at least one zero of f2y when y∈b,∞. The derivative of f2y is given by(30)df2ydy=1+2m1y3+1y2+13/2−3y2y2+15/2+2m˜1y+b3+1y−b3.Following the same procedure as given for0<y<b, one can easily prove the uniqueness of the equilibrium solution for y>b.We discuss here two special cases of CC forb∈1/3,1.1394282249562009. When b=1/3, then the masses on the horizontal axis are zero, i.e., m=0. In this case, we only get the two equilibrium points along y-axis. When b=1.1394282249562009, then m˜=0. In this case, we get four equilibrium points along x-axis. The positions of the masses and the corresponding equilibrium points for these two cases are shown in Figure 8. The stability of these cases will be discussed in Section 6.Figure 8
(a) Equilibrium points (red color) alongy-axis when m=0. (b) Equilibrium points (red color) along x-axis when m˜=0.
(a)(b)


## 5.4. Equilibrium Solutions off the Coordinates Axes

It is numerically confirmed that, forb=0.67,0.97,and1.13, there are always a total of 12 equilibrium points. As shown in Figure 9, four of the equilibrium points are on the x-axis, four on the y-axis, and remaining four of the equilibrium points are off the axes. Since the gravitational field is a function of mass parameters mb, therefore the equilibrium points change their positions around the primaries for changing values of b. It is numerically confirmed that majority of equilibrium points are around the primaries along the horizontal axis if b is around 1/3 as the masses on the horizontal axis are dominant (see Figures 9(a)). For b≥1, the equilibrium points concentrated around the primaries on the vertical axis (see Figure 9(b)).Figure 9
(a–c) 12 equilibrium points (red colors) for different values ofb=0.67,0.97,and1.13, respectively.
(a)(b)(c)


## 6. Stability Analysis

To study the stability of the equilibrium points obtained in the previous section, we will follow the standard linearization procedure by linearizing the equation of motion of the infinitesimal mass. Let the location of an equilibrium point in the RR6BP be denoted byx,y, and consider a small displacement X,Y from the point such that x+X and y+Y are the new position of the infinitesimal. Using Taylor’s series expansion in equations (17) and (18), we obtain a new set of second-order linear differential equations:(31)X¨−2Y˙=XUxx+YUxy,Y¨+2X˙=XUxy+YUyy.The matrix form of the linearized equations is(32)X˙Y˙X¨Y¨=00100001UxxUxy02UxyUyy0−2XYX˙Y˙.These equations can also be written in the following matrix form:(33)Ψ.=AΨ,where(34)Ψ=X˙Y˙X¨Y¨,A=00100001UxxUxy02UxyUyy0−2.The characteristic polynomial forA is(35)Λ4+αΛ2+β=0,where α=4−Uxx−Uyy and β=UxxUyy−Uxy2. Let Λ2=λ; then, equation (35) reduces to(36)λ2+αλ+β=0.Now, in order for a Lagrange point to be linearly stable to a small perturbation, all four roots,Λ, of equation (35) must be purely imaginary. Thus, in turn, it implies that the two roots of equation (36).(37)λ±=−α±α2−4β2,must be real and negative. For λ±<0, we must have(38)iα>0and0<β≤α24,or,iiα>0andα2−4β=0.We will numerically identify regions when either condition (i) or condition (ii) is satisfied. Figures10(a) and 10(b) (left to right) give the stability region in case (i) and case (ii), and their projections for b=0.67, b=1, and b=1.13 are shown in Figures 11 and 12, respectively.Figure 10
(a-b) Shaded regions represent stability regions for casesi and ii, respectively.
(a)(b)Figure 11
Projection of stability regions for fixed values of the parameterb=0.67,1,and1.13.
(a)(b)(c)Figure 12
Case (ii): projection of stability regions for fixed values of the parameterb=0.67,1,and1.13.
(a)(b)(c)We have tested a large number of equilibrium points for many values ofb and found all of them unstable. In other words, the intersection of Ux=0andUy=0 within α>0and0<β≤α2/4 and α>0andα2−4β=0 is an empty set. Representative examples are given in Tables 1–3.Table 1
Equilibrium points and stability analysis forb=0.67.
Equilibrium pointEigenvaluesStabilityL1,2=±1.342241,0±1.660648,±1.525769iUnstableL3,4=±0.678932,0±1.943106,±1.168505iUnstableL5,6=0,±1.472145±1.218669,±1.351617iUnstableL7,8=0,±0.180864±6.347782,±4.5925990iUnstableL9,10,11,12=±1.023953,±0.479505∓0.913372±0.998959iUnstableTable 2
Equilibrium points and stability analysis forb=0.97.
Equilibrium pointEigenvaluesStabilityL1,2=±1.637980,0±1.480313,±1.458732iUnstableL3,4=±0.509436,0±3.583445,±2.544469iUnstableL5,6=0,±1.660744±1.454216,±1.451670iUnstableL7,8=0,±0.464669±4.047077,±2.908229iUnstableL9,10,11,12=±0.950136,±0.874620∓0.909938±0.998607iUnstableTable 3
Equilibrium points and stability analysis forb=1.13.
Equilibrium pointEigenvaluesStabilityL1,2=±1.751564,0±1.5080962,±1.499923iUnstableL3,4=±0.495101,0±4.863781,±3.594754iUnstableL5,6=0,±1.413901±1.797668,±1.600502iUnstableL7,8=0,±0.874617±2.535781,±1.929722iUnstableL9,10,11,12=±0.423616,±1.103657∓0.903618±0.9990127iUnstable


## 7. Conclusions

We have modeled and studied the Rhomboidal Restricted Six-Body Problem which has four masses at the vertices of the rhombus, and the fifth mass is at the intersection of the two diagonals placed at the origin. It is assumed thatm1=m2=m0=m and m3=m4=m˜. The primaries always move in the rhomboidal configuration. The sixth mass m5≪mi, where i=0,…,4, is moving in the gravitational field of the five primaries. To get rid of the time dependency, the equation of motion of m5 is written in the rotating coordinate system. It is shown that, for the mass parameter b∈1/3,1.1394282249562009, the total number of equilibrium points is always 12. Out of the 12 equilibrium points, 4 are on the x-axis and 4 are on the y-axis. As the mass parameters vary, the equilibrium points also change their positions around the five primaries, and the number of equilibrium points remains same. The linear stability analysis revealed that none of the equilibrium points are stable. The permissible region of motion is also discussed according to the variation of Jacobian constant of C. We have identified the values of C at which the permissible regions of motion are partially disconnected and fully disconnected.



---

*Source: 5575826-2021-07-01.xml*