Abstract

We discuss the restricted 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 and . 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 and . The masses and are written as functions of parameters and 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 and identify the value of Jacobian constant 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 , there always exist 12 equilibrium points. We also show that all 12 equilibrium points are unstable.

1. Introduction

In the restricted n-body problem, a body of infinitesimal mass moves under the gravitational influence of 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 -body problem in celestial mechanics. It can be used to find simple or special solutions of the -body problem since the geometry formed by the arrangement of the primaries remains constant for all time.

The CC for has been investigated extensively after being identified as the problem of the century by Arnol’d [1]. Xia [2] used the method of analytical continuation to find exact number of central configurations for 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 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 . 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., 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 and equals the length between and , 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 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 , that is, if the value of is less than or equal to a critical value , 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 , , and having center at the origin. The libration points that lie on circles and are unstable, while there are some stable libration points on circle .

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 of and are located on the horizontal axis, and two other equal mass primaries are located on the vertical axis with masses . 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 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 the n-body problem has the following form:where the units are chosen so that the gravitational constant is equal to one. A central configuration is a particular configuration of the bodies where the acceleration vector of each body is proportional to its position vector, and the constant of proportionality is the same for the bodies. Therefore, a CC is a configuration that satisfies the following equation:where is angular speed and which is the center of mass of the bodies. We take the position vector (see Figure 1) of the five primaries , where as

For in equation (2), the CC equation for masses and iswhere and . The positions of the primaries are taken from equation (4). Writing equations (5) and (6) in components’ form, we get

Because of the symmetry of the problem, the equations for and are identical to equations (7) and (8). The equation for is identically zero. We have extended the astrophysical model given in [29] by keeping the mass at the center. The results of the model given in [29] can be reproduced by taking in our model. Without loss of generality, we take and solve equations (7) and (8) simultaneously for and :where

Defining , then equations (9)–(11) give

Again, let in equation (13); we get an alternate form of equation (13):

Lemma 1. given by (13) is a continuous and strictly decreasing positive function in the interval ,, and .

Proof. One can easily check the continuity and and directly from equation (14). It means the value of which is the mass ratio will remain positive i.e., as the ratio of the distance parameters varies between . Now we need to prove that the function is decreasing and positive for the given interval. Taking the derivative of , 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 , then ; from equation (13), we get the shape of true square with and (ii)If , then ; from equation (13), we get and (iii)If , then ; from equation (13), we get and There are two cases in which the RR6BP central configuration degenerates.(i)If , then ; this implies and (ii)If , then ; this implies and

3. Equations of Motion

In this section, we describe the motion of the infinitesimal body, , 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 . 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 arewhere dot represents the derivative with respect to time. From here onward, without loss of generality, we take the value of . The equations of motion of in the synodical coordinates and arewherewhich is the effective potential. The mutual distances of from the primaries in the corotating frame are

The Jacobian constant is given by Curtis [30]:

For a given value of the Jacobi constant, is only a function of position in the rotating frame. Since cannot be negative, it must be true that

The boundaries between forbidden and allowed regions of motion are found by setting , i.e.,

It is now trivial to show that is the first integral of motion of system (11) by proving that .

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 are known as the zero velocity curves.

The zero velocity curves when and when are given in Figure 4, and the corresponding Hill regions are given in Figure 3. We also give the region of possible motion of for six different values of Jacobi constants in Figures 57 for mass parameters . The shaded regions represent the forbidden regions of motion for the infinitesimal mass . It is numerically confirmed that the permitted regions are completely disconnected for , as shown in Figures 57, for the above values of . For the increasing values of , the allowed region of motion (white region) becomes partially connecting at and completely connected at . It can be seen from Figures 57 that the transition of motion from totally disconnected to completely connected occurs in six stages for . For these values of , can freely move in the gravitational field of the CC region for and cannot reach any of the primaries for .

5. Equilibrium Solutions

Equilibrium solutions of the RR6BP are the solutions of and . The derivative and of the effective potential, given in equation (19), are found as follows:

5.1. Equilibrium Solutions on the Coordinates’ Axes

Since the potential given in equation (19) is invariant under the symmetry , , and , we will restrict our computation to the first quadrant: and . 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 ; then, equations (24) and (25) are given as

To solve , divide into subintervals and .

5.2. ;

Rewrite the right hand side of equation (25) by taking into account that :

At , ; and, at , ; therefore, by the mean value theorem, there is at least one zero of when . The derivative of is given by

The only term that can make negative for is dominated by the rest of the term in equation (26); therefore, . This proves the existence of the unique equilibrium solution inside rhombus on the -axis.

5.3.

Now, consider the case . Use equation (26) and rewrite as

At , ; and, at , ; therefore, by mean value theorem, there is at least one zero of when . The derivative of is given by

Following the same procedure as given for , one can easily prove the uniqueness of the equilibrium solution for .

We discuss here two special cases of CC for . When , then the masses on the horizontal axis are zero, i.e., . In this case, we only get the two equilibrium points along y-axis. When , then . 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.

5.4. Equilibrium Solutions off the Coordinates Axes

It is numerically confirmed that, for 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 , therefore the equilibrium points change their positions around the primaries for changing values of . It is numerically confirmed that majority of equilibrium points are around the primaries along the horizontal axis if is around as the masses on the horizontal axis are dominant (see Figures 9(a)). For , the equilibrium points concentrated around the primaries on the vertical axis (see Figure 9(b)).

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 by , and consider a small displacement from the point such that and 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:

The matrix form of the linearized equations is

These equations can also be written in the following matrix form:where

The characteristic polynomial for iswhere and . Let ; then, equation (35) reduces to

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).must be real and negative. For , we must have

We will numerically identify regions when either condition (i) or condition (ii) is satisfied. Figures 10(a) and 10(b) (left to right) give the stability region in case (i) and case (ii), and their projections for , , and are shown in Figures 11 and 12, respectively.

We have tested a large number of equilibrium points for many values of and found all of them unstable. In other words, the intersection of within and is an empty set. Representative examples are given in Tables 13.

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 that and . The primaries always move in the rhomboidal configuration. The sixth mass , where , is moving in the gravitational field of the five primaries. To get rid of the time dependency, the equation of motion of is written in the rotating coordinate system. It is shown that, for the mass parameter , 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 . We have identified the values of at which the permissible regions of motion are partially disconnected and fully disconnected.

Data Availability

No data were used in the study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.