Abstract

We explore the central configuration of the rhomboidal restricted six-body problem in Newtonian gravity, which has four primaries (where ) at the vertices of the rhombus , , , and , respectively, and a fifth mass is at the point of intersection of the diagonals of the rhombus, which is placed at the center of the coordinate system (i.e., at the origin ). The primaries at the rhombus’s opposite vertices are assumed to be equal, that is, and . After writing equations of motion, we express , and in terms of mass parameters and . Finally, we find the bounds on and for positive masses. In the second part of this article, we investigate the motion and different features of a test particle (sixth body ) with infinitesimal mass that moves under the gravitational effect of the five primaries in the rhomboidal configuration. All four cases have 16, 12, 20, and 12 equilibrium points, with case-I, case-II, and case-III having stable equilibrium points. A significant shift in the position and the number of equilibrium points was found in four cases with the variations of mass parameters and . The regions for the possible motion of test particles have been discovered. It has also been observed that as the Jacobian constant increases, the permissible region of motion expands. We also have numerically verified the linear stability analysis for different cases, which shows the presence of stable equilibrium points.

1. Introduction

The system of point masses is considered to be under motion due to the mutual gravitational force as explained by Newton’s gravitational law. Dynamical systems with point masses have been extensively studied, and various models have been proposed for research aiming at approximating the behaviour of real celestial systems. There are many reasons for studying the -body problem since it is known that approximately two-thirds of the stars in our Galaxy exist as part of multistellar systems.

In celestial mechanics, the central configuration (CC) plays an important role in the investigation of -body problems. It is the configuration that can be used to find simple or special solutions of the -body problem since the shape of the figure formed by the arrangement of the bodies remains constant for all time. That is why the CCs lead to homographic solutions. The first three collinear homographic solutions for were found by Euler in 1767, and two equilateral triangular homographic solutions were found by Lagrange [1] in 1772 for . The central configuration for has been investigated a lot. Xia [2] used the analytical continuation approach to find the exact number of CCs for an open set of positive masses in 1991. Llibre and Mello [3] investigated 7-body problem for the existence of families of central configurations. In 2017, Deng et al. [4] investigated the central configuration of four-body problem with equal masses and showed that for the planar Newtonian four-body problem having adjacent equal masses, that is, , and equal lengths for the two diagonals, any convex noncollinear central configuration 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 central configuration is also an isosceles trapezoid. In the same year, Marchesin [5] explored the stability of rhomboidal configurations with a mass in the center. He considered a system with five bodies at , , , , and with masses , and that are moving in a plane. They investigated the effect of altering the mass of the central body on the stability of such a configuration for several values of mass parameter . Saari [6] studied the role and properties of central configuration for -body problem, and the other work on central configuration is given in Refs. [715] (for latest publications on CCs, cf. Suraj et al. [16], Zotos, and Papadakis [17], and Cornelio et al. [18]).The three-body problem has been a great challenge for the scientists as it needed some special assumptions for the third body. The ever first problem solved using restriction for the infinitesimal body is “restricted three-body problem” (RTBP). The RTBP originally arose from the work of Newton, but unfortunately, Newton was unable to solve the problem throughout his life. Alexis Clairaut [19, 20] in 1747 solved the problem with some approximations, and in 1757, he calculated the return of “Halley’s Comet.” Ollongren [21] in 1988 was the first to introduce the restricted five-body problem. Ollongren investigated the motion of the fifth infinitesimal mass that moved under the gravitational attraction of the other four massive bodies. He discovered nine equilibrium points, three of which were stable, and the rest were linearly unstable. In 2007, Papadakis and Kavanos [22] extended the work of Ollongren and investigated the restricted five-body problem numerically, and they showed that the counts of equilibrium points depend on the radiation factors as well as on the mass parameter. Marchesin and Vidal [23] in 2013 studied the stability of five-body problem by assuming that the central body makes a generalized force on the other four bodies. In 2018, Aggarwal et al. [24] investigated the effect of small perturbation in Coriolis and centrifugal forces on the presence of equilibrium points in the restricted four-body problem. The three primaries were placed at the vertices of equilateral triangle and the fourth infinitesimal mass moves under the gravitational effect of the primaries. They showed the existence of 8 to 10 equilibrium points and found that all the liberation points are unstable. Mello et al. [25], in 2009, introduced the stacked CC for the spatial six-body problem, and they illustrated three novel families of stacked spatial central configurations for the six-body problem with a regular tetrahedron model. Recently, in 2020, Idrisi and Ullah [26] investigated the existence and linear stability of equilibrium points in a six-body problem with a square CC and a central mass, and discovered that there are twelve equilibrium points, four of which are stable for a given value of the mass parameter . Furthermore, in 2020 [27], Ansari et al. studied six body CC by placing the primary at the vertices of a square and a fifth mass at the center of a circle, which is thought to be the orbit of motion for the rest of the primaries. Using the Jeans law and the Mesheherskii space time transformation, they discovered twelve equilibrium points for the sixth small body of variable mass, all of which are unstable. Kulesza et al. [28] in the restricted rhomboidal five-body problem found that the number of libration points depends on the semidiagonal ratio and they showed that there can be eleven, thirteen, and even fifteen libration points. In 2013, Marchesin and Vidal [29] studied the restricted rhomboidal five-body problem and found that no chaos exists and the behaviour of the fifth mass is quite predictable. They also showed the existence of the periodic solutions of infinitesimal mass and also studied numerically their linear horizontal stability. The rhomboidal six-body CC was presented by Siddique et al. in 2021 [30] by putting the four primaries at the corners of the rhombus and the fifth mass at the intersection of the two diagonals of rhombus. They computed the region of probable motion of the sixth mass using the first integral of motion and the values of the Jacobian constant for various energy intervals and found the limitation on the region of motion for infinitesimal mass. They demonstrated the existence and uniqueness of equilibrium solutions on and off the axes using semi-analytic approaches, and established that there are always twelve unstable equilibrium points when and . Recently, Ref. [31] examined the basins of convergence by deploying the well-known Newton–Raphson iterative scheme, associated with the libration points (indeed, act as attractors), in the restricted rhomboidal six-body problem.

In this study, we focus on the rhomboidal restricted six-body problem, which has four primaries at the corners of the rhombus and a fifth primary where the diagonals of the rhombus cross, which is in the center of the coordinate system, that is, . There is a sixth test mass that moves under the gravitational effect of the five primaries. The primaries at the opposite vertices of rhombus are assumed to be equal, that is, and . In section 2, we investigate the continuous families of five primaries central configurations for different values of the mass parameters and . In section 3, we derived the equation of motion for infinitesimal mass and discovered the Jacobian constant . In section 4, we explore the Hill’s and permissible regions of motion of the . Sections 5 and 6 analyze the positions and stability of the equilibrium points for the sixth body, . Concluding remarks are given in section 7.

2. Rhomboidal Central Configurations

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 central configuration is a configuration that satisfies the following equation:where is an angular speed andis the center of mass on -bodies. We choose the coordinates (for positions, see Figure 1) of five primaries , where as

Expanding equation (2) for and using equation (3) and letting and , we get the following equation for central configurations for the masses and :and for the masses and , the equation of motion is

The mass is stationery so its equation of motion ends up zero. We take sum of all primaries is equal to unity that is

Taking and solving equations (5)–(7) givewhere

Lemma 1. For any and , is always negative.

Proof. Let us define the following:To prove , we need to prove . For this, differentiate with respect to and find the critical points of in . There is only one critical point of (i.e., ) in and . One can easily see that is monotonically increasing function because in and is monotonically decreasing function because in (see Figure 2). When , . So for , and hence, for and .
To prove and positive, we need to prove and must be negative for and . Because and are nonlinear algebraic functions of and , so it is difficult to solve these inequalities. For this, we draw the region (shaded region of Figure 3) where both and are negative. From the graph, we can easily find the approximate bounds for and where and are negative.
Using equation (7) to equation (9), one can easily see that and , and for and . In Figure4 we show the region of central configuration for which are positive. Because Figure 4 is symmetric about the line , we divide here Figure 4 into two parts: upper region and lower region . Here, we discuss only the central configuration for the region of . The upper CC region is approximately surrounded by the following three interpolating polynomials:Figure 5 shows the changes of masses , , and with variations of parameter . We can clearly see that as varies in , the mass is increasing, and other masses are becoming zeros.

3. Equation of Motion of Infinitesimal Body

In this section, we derive the equation of motion of the infinitesimal body, , that moves under the gravitational attraction of the five primaries. The sixth body has a significantly smaller mass as compared to the masses of the primaries, that is, where . Due to this fact, the sixth body acts as an infinitesimal test particle, and therefore, it does not have any influence on the motion of the five primaries. In the above scenario, the equation of motion of the test particle is

The dot represents the derivative with reference to time. We discuss here the dynamics of the infinitesimal mass with respect to the five primaries. The equations of motion of the in synodical coordinates and [32] arewhereis the effective potential, where the mutual distances are

The effective potential is shown in Figure 6 (right) for the four different cases of rhomboidal restricted six-body problem. We define the first Jacobian-type integral by

By demonstrating that , it is now straightforward to establish that is the first integral of motion of system (18).

Equation (18) shows that . Then, will establish a boundary between both the allowed and forbidden regions and presents the zero velocity curves for various values of .

4. The Spheres of Influence

Spheres of influence or gravitational spheres of influence are areas surrounding celestial objects where other celestial objects experience the greatest pull and can become satellites of the huge celestial object relative to its mass; these regions are also known as Hill’s regions, after George William Hill. The zero velocity curves (ZVC) are the contours of the Jacobian constant as mentioned above, and they are available in Figure 6 on the left, with their corresponding effective potentials on the right. The Hill’s regions are tightly packed circular regions surrounding primaries; in Figure 6 (i), for case-I when and , the Hill’s region for stretches between and , whereas for , it extends between and . Figure 6 (i) also shows that the spheres of influence shrink as the particle mass decreases, since has a relatively tiny Hill’s region as compared to the other primaries. The pink and green lines in Figures 6 (i), (iii), (v), and (vii) on the left depict the contours of equations (19) and (20). The black dots indicate masses and are represented by the symbol ; , the orange dots represent equilibrium points and are denoted by ; . Figures 6 (ii), (iv), (vi), and (viii) on right illustrate the effective potentials for each of the four cases, where the closely packed circular regions around the potentials of each mass, which resemble chimneys, can be seen clearly. The permissible regions of motion for the infinitesimal mass are shown in Figures 710 for each of the four cases for different values of mass parameters and . Six figures for each of the four cases show the forbidden (shaded) and permissible (white) regions of motion of for various values of . The regions of motion that are permitted have the same behaviour in all four cases; that is, the permissible region increases as the value of decreases. It has also been numerically confirmed that the permissible regions of motion are connected and are around the primaries for higher values of , that is, at , and . The permissible motion regions are mainly around the four primary and are isolated from one another at , and for the above values of , the motion of will be bounded around any of the four primaries and will not be able to escape form there normally. The permitted regions become totally detached when , and the infinitesimal mass may freely move in the gravitational field of the primaries.

5. Equilibrium Solutions

All rates of change should be zero for equilibrium solutions; hence, the right-hand side of the system in (15) can be set to zero, that is, and , and the solution of resulting equations will give the problem’s equilibrium solutions. and are

We considered four cases with different values of and that show a significant change in the number and location of equilibrium points, namely, and for case-I; and for case-II, and for case-III; and and for case-IV. The corresponding equilibrium points can be seen in subfigures (i), (iii), (v), and (vii) of Figure 6.

5.1. Equilibrium Solutions: On the Coordinate Axes

We shall limit our investigation to the first quadrant, and , because the potential given in equation (15) is unchanged under the symmetry , , and . To determine the presence and number of equilibrium solutions on the y-axis, we set and then write equations (19) and (20) as

The axis is subdivided into and to figure out for equilibrium solutions. means collisions of with or from inside. We will not discuss collisions cases here.

5.2.

Equation (21) is rewritten as

For and , when , and when , ; thus, the mean value theorem implies that has at least one zero when . To verify the existence and uniqueness of equilibrium solutions inside the interval , we take derivative of equation (22):

in equation (23) is the only term that can make the derivative negative for , and as a result, , which indicates that the equilibrium points in the interval exist and are unique.

5.3.

When , equation (21) is rewritten as

For and , when , , and when , ; thus, the mean value theorem implies that has at least one zero when . The derivative of equation (24) with respect to is

The equilibrium solution’s uniqueness for can be easily proved using the same technique as for .

Then, there are total 16 equilibrium points for case-I and out of which are along axis, while are on axis; similarly for case-II, III, and IV, there are four equilibrium points on each coordinate axis for distinct values of and , but the number of equilibrium points is different for other cases, which can be seen from the subfigures (i), (iii), (v), and (vii) of Figure 6.

5.4. Equilibrium Solutions: Off the Coordinate Axes

The mass parameters and , which we have numerically confirmed, determine the number and location of the equilibrium points, which is illustrated in subfigures (i), (iii), (v), and (vii) of Figure 6. For case-I, there are eight off-coordinate axis equilibrium points; for case-II, there are 12 off-coordinate axis equilibrium points, the majority of which are clustered around the central mass . Cases III and IV contain 4 equilibrium points that cluster around and off the coordinate axis.

6. Stability Analysis

We represent the location of an equilibrium point in our problem by the coordinates and consider a small displacement from the equilibrium point to be in order to linearize around the equilibrium point, and then, the new location of the equilibrium points will be . The system in equations (14) is subjected to Taylor series expansion, yielding the following set of linearized equations:

The linearized equations in the matrix form is

Let

The equation in (27) can be written as

From , we havewhere and .

Let , then equation (30) reduces to

For an equilibrium point to be linearly stable under a slight disturbance, all four roots of equation (30) must be completely imaginary. As a result, the two roots of equation (31)must be real and negative. For , we must have

The regions of stability when either condition (i) or (ii) holds true have been determined and presented in Figures 1114. The stability areas corresponding to the four cases for different values of and are shown in Figures 15 and 16, where is fixed at , and the regions are projected across the whole domain of the mass parameter , that is, . We investigated the stability of the equilibrium points found in the preceding section for each of the four different cases, and many off-coordinate stable equilibrium points were discovered; the results are reported in Tables 14.

7. Conclusion

Mass parameters and are used to examine the restricted rhomboidal six-body problem in depth. A point mass, , is put at the intersection of two diagonals of the rhombus, and the remaining four point masses are arranged in such a way that they preserve a rhomboidal central configuration during their two-by-two motion. The point masses at the opposite vertices of a rhombus are supposed to be the same, that is, and . In a synodical coordinate system, the motion of the sixth test mass, , which has an infinitesimal mass relative to the other five masses, that is, , is investigated for both the equilibrium points and their linear stability. A significant shift in the position and number of equilibrium points was found in four cases with the variations of mass parameters and . Cases-I to IV have 16, 20, 12, and 12 equilibrium points respectively, with case-I and case-II having eight off-the-coordinate-axis stable equilibrium points, case-III having four off-the-coordinate-axis stable equilibrium points, and case-IV having no on or off-the-coordinate-axis stable equilibrium points. The regions for the possible motion of have been discovered, and it has been observed that as the Jacobian constant increases, the permissible region of motion expands, and the values of for which the regions of the possible motion become disconnected or partially disconnected have also been discovered. We also have numerically verified the stability regions for different cases, which shows the presence of stable equilibrium points.

Data Availability

No data were used for the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.