With the lack of analytical solutions of Equations (
6), (
7), (
10), and (
11), we used
MATLAB to determine the points on the
,
parameter plane that satisfy these equations. The regions where the kite central configurations of four bodies can exist are already known from [
6] and are also shown here as the colored areas in
Figure 3 (for the convex) and in
Figure 4 (for the two concave cases). Therefore, we only considered solutions in the given domains. The resulting curves appear in the
,
parameter plane as the graphs, i.e., the
pairs, of such implicit functions
that satisfy the above equations. In the following, we describe our results separately for the convex and concave cases.
3.1. Convex Case
The solutions in the convex case are displayed in
Figure 3. All convex central configurations, with the assumptions of [
6], are restricted to the interior of a triangle (colored area in
Figure 3), bordered by the critical lines
,
, and
, along which the masses take extremal values as indicated in the figure.
The square configuration [
13] with four equal masses at
(point
G in
Figure 3) can be considered as a generating configuration of the three-equal-mass cases, since both the
and
solutions start from here. (The solution points
are obtained either from Equation (
6) or from Equation (
10) by substituting
.)
In the
case (blue curve in
Figure 3), the masses of the bodies
A,
E, and
are equal, and as
increases from
to
and
from
to
, the initial square configuration is replaced by a convex deltoid, which is more and more elongated along its line of symmetry. The endpoint
in
Figure 3 is on the critical line
, and the value of
for this point is obtained from Equation (
6) by substituting
. One observes that
is zero along this critical line, thus, at the endpoint
of the blue curve, the four-body problem is reduced to a three-body problem with the bodies
A,
E, and
forming an equilateral triangle and having the common mass
. We note here an important property, namely that in the convex
case, for each
there exists a unique
satisfying the condition of three equal masses. In other words, it means that the function
, associated with the blue curve, is injective, for which it is enough to be strictly monotone. The proof of this is given in
Appendix B.1.
In the
case, the solution is represented by a red curve in
Figure 3. In contrast to the previous case, the corresponding function
has a local minimum at
,
(see
Appendix B.2), implying it to be noninjective. It follows from this finding that between this minimum point and the endpoint of the red curve on the critical line
(point
in
Figure 3), for each value of
there exist two values of
for which
. This interval of
is rather narrow, but the corresponding interval of
is wide:
. For example, for
there is a
with corresponding masses
,
; and another
with corresponding masses
,
. The values of
and
at the endpoint
can be computed from Equation (
10) by substituting
as obtained from the equation of the respective critical line. At this endpoint, the masses are
,
; consequently,
B,
E, and
become zero-mass bodies, and the configuration is reduced to a one-body problem with
A of mass 1.
In [
11], the central configurations of four bodies with three equal masses were studied by a different method and from a different point of view. As a main result, the authors determined the number of all possible types of configurations depending on the ratio of the nonequal mass with respect to the equal mass of the bodies. Here, we study the problem of three equal masses with respect to the general case of the axisymmetric central configurations of four bodies, with an emphasis on the physical aspects of the configurations. Though the methods and the scopes are different in the two studies—in [
11], the problem was studied in itself—some of the results can be compared. The abovementioned points
and
appear in our approach as intersections of the solution curves and two of the critical lines in
Figure 3. They are also among the 7 special points of [
11]:
and
. By the transformation
,
between the angles
,
and the coordinates
k,
l of [
11], the identity of the two solutions can be seen.
3.2. Concave Cases
The results in the concave cases are shown in
Figure 4. According to [
6], the concave kite central configurations are confined to the interiors of two triangles (colored areas in
Figure 4,
and
denoting the first and second concave cases, respectively). The bordering critical lines are
,
,
, and
, where the masses take extremal values as indicated in the figure.
The common point of the two triangles, labeled by
S, refers to a singular point, where
and
are individually undetermined but for the sum of them holds the condition
. From the definition of
, it also follows that
at the singular point, thus the masses of
A,
E, and
are equal for any values of
, satisfying the above condition. Therefore, the blue curve, corresponding to the
case, necessarily goes through
S. However, the red curve that corresponds to the
case also crosses this point, since
is also possible here—that is, when all four masses are equal:
. In both cases,
A,
E, and
form an equilateral triangle (since
) and
B occupies its center (since
). We note that apart from the singular point
S, there exists one more point allowing four equal masses (see also [
13]), namely, at
,
, where the blue and red curves cross each other. (The above values were obtained by equating the left-hand sides of Equations (
7) and (
11).) The two curves stay very close after point
S; therefore, we show a magnification of the two intersections in
Figure 5.
The conditions and have solutions both in the first and second concave cases.
The blue curve in
Figure 4 represents the solutions of
. Its starting point is labeled by
, and it is on the critical line
in the first concave region. The other coordinate
was computed from Equation (
7), using the substitution
. At this point,
is zero (as is all along the critical line), and thus the common mass of the bodies
A,
E, and
is necessarily
. The configuration is a Lagrangian equilateral triangle. After leaving
, the blue curve (associated with function
) first decreases to a minimum at
,
, then it increases again. Since the initial value of
is reached again at the singular point
S, the property that for each value of
there exist two values of
for which
, holds along the blue curve in the whole first concave region. For example, for
, there is a
with corresponding masses
,
; and another
with corresponding masses
,
. (This noninjectivity of
in the first concave region is shown in
Appendix B.3) In the second concave case, the blue curve continues to increase (the proof is given in
Appendix B.3 by showing the positive sign of the derivative of
in this region), and it ends on the critical line
. The coordinates
,
of this endpoint (
in
Figure 4) can be computed from Equation (
7) with the substitution
. In this point,
and
, in accordance with the property of the respective critical line, thus the masses of
A,
E, and
are 0 and the system is reduced to a one-body problem of
B with unit mass.
The red curve in
Figure 4 represents the solutions of
. Its starting point is labeled by
and it is on the critical line
. The other coordinate
was computed from Equation (
11) using the substitution
. At this point,
is zero (as is all along the critical line), and thus the common mass of the bodies
B,
E, and
is
. The configuration is a Eulerian collinear three-body problem. In contrast to the previous
case, the function
—associated with the red curve—strictly monotonically increases in all points of its domain of definition, implying to be injective both in the first and second concave regions. (See
Appendix B.4 for the proof.) The endpoint (
in
Figure 4) falls on the critical line
, thus its other coordinate
is obtained from Equation (
11) with the above substitution. In this point,
, in accordance with the property of the respective critical line, thus the masses of
B,
E, and
are
again, but now they form a Lagrangian equilateral triangle.
The endpoints
,
,
, and
, as the intersections of the solution curves and critical lines in
Figure 4, are related to the special points of [
11]:
,
,
, and
. In the case of
(concerning the points
and
), the transformation between the angle coordinates
,
and the coordinates
k,
l of [
11] is
,
. In the case of
(regarding the points
and
), the transformation is
,
. We note that the point
of [
11] corresponds to the singular point
S in
Figure 4.