Advances in Astronomy
Volume 2015 (2015), Article ID 473483, 21 pages
http://dx.doi.org/10.1155/2015/473483

Research Article

Equilibrium Points and Related Periodic Motions in the Restricted Three-Body Problem with Angular Velocity and Radiation Effects

E. A. Perdios,¹ V. S. Kalantonis,¹ A. E. Perdiou,² and A. A. Nikaki³

¹Department of Electrical & Computer Engineering, University of Patras, 26500 Patras, Greece
²Department of Civil Engineering, University of Patras, 26500 Patras, Greece
³Department of Computer Engineering & Informatics, University of Patras, 26500 Patras, Greece

Received 24 April 2015; Accepted 2 June 2015

Academic Editor: Zdzislaw E. Musielak

Copyright © 2015 E. A. Perdios et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The paper deals with a modification of the restricted three-body problem in which the angular velocity variation is considered in the case where the primaries are sources of radiation. In particular, the existence and stability of its equilibrium points in the plane of motion of the primaries are studied. We find that this problem admits the well-known five planar equilibria of the classical problem with the difference that the corresponding collinear points may be stable depending on the parameters of the problem. For all planar equilibria, sufficient parametric conditions for their stability have been established which are used for the numerical determination of the stability regions in various parametric planes. Also, for certain values of the parameters of the problem for which the equilibrium points are stable, the short and long period families have been computed. To do so, semianalytical expressions have been found for the determination of appropriate initial conditions. Special attention has been given to the continuation of the long period family, in the case of the classical restricted three-body problem, where we show numerically that periodic orbits of the short period family, which are bifurcation points with the long period family, are connected through the characteristic curve of the long period family.

1. Introduction

The restricted three-body problem is the most celebrated problem of Celestial Mechanics (see, among others, Szebehely [1], Marsden and Ross [2], and references therein) and due to its interdisciplinary applications it has also attracted the interest of many researchers from different fields of science (see, e.g., Prosmiti et al. [3], Vrahatis et al. [4], Sano [5]). During the last century, several modifications of this classical problem have been introduced in order to make it more relevant and applicable to certain systems of Dynamical Astronomy. In our present study, we will consider the combination of two such modifications in which the variation of the angular velocity (the Chermnykh’s problem) and the radiation effects of the primaries (the photogravitational problem) are considered.

The Chermnykh’s problem is a generalization of the Euler’s problem of two fixed gravitational centers and the restricted three-body problem, in which the third body of negligible mass moves in the orbital plane of a dumbbell which is rotating with constant angular velocity around the center of mass (Chermnykh [6]). This problem admits applications in both Dynamical Astronomy and Celestial Mechanics as well as in Chemistry (see Goździewski and Maciejewski [7] and references therein) and has been studied recently by several authors. In particular, Goździewski and Maciejewski [7] investigated the nonlinear Lyapunov stability of the triangular equilibrium points. Perdios and Ragos [8] studied asymptotic motions around the collinear equilibrium points while Papadakis [9] considered planar and three-dimensional periodic motions around them. Perdiou et al. [10] studied three-dimensional motions through bifurcations of the Sitnikov family. In recent years, certain variants of this problem have been also proposed by other authors. Specifically, Jiang and Yeh [11] investigated the equilibrium points of the Chermnykh’s problem in the case where the value of the angular velocity is determined by the gravitational potential of the belt which was considered to be a part of the potential function. For the same variant Abouelmagd et al. [12] studied the linear stability of the triangular equilibrium points as well as periodic solutions around them when the primaries are radiating oblate spheroids. Recently, Singh and Leke [13, 14] considered also the oblateness of the test particle as well as the effect of small perturbations in the Coriolis and centrifugal forces.

On the other hand, an interesting modification of the classical problem is the photogravitational restricted three-body problem in which the repulsive force of the radiation is also considered in the potential function and it was introduced for studying the specific three-body problem of Sun, planet, and a dust particle (Radzievskii [15]). This is a simplified model, since it does not take into account the perturbation force of the Poynting-Robertson effect but only the radiation pressure as a reduction factor of the magnitude of the inverse square law forces of the primaries upon the massless body. However, it may be useful for the study of motion of extremely small particles such as dust grains and interplanetary drifters (see, e.g., Kalvouridis [16], Bewick et al. [17], and references therein). In the framework of this simplified model many works have been done in the last years. For example, Simmons et al. [18] studied the existence and stability of the equilibrium points for all values of the radiation factors and the mass parameter, Goździewski et al. [19] considered the nonlinear stability of the triangular equilibrium points, Elipe and Lara [20] studied three-dimensional periodic orbits, and Perdios et al. [21] investigated the bifurcations of plane to three-dimensional periodic orbits, while Perdios [22] determined series of horizontally critical symmetric periodic orbits as well as their vertical stability. Also, Zimovshikov and Tkhai [23] studied the stability of the collinear equilibrium points in the elliptic case of the corresponding problem, Abouelmagd [24] considers the additional effects of the oblateness as well as the Coriolis and centrifugal forces on the location and stability of the equilibrium points, while Singh and Taura [25] investigated the linear stability of the triangular equilibrium points by taking into account the oblateness of the smaller primary and the gravitational potential from a circumbinary belt. Recently, Yárnoz et al. [26] examined the feasibility of manipulating asteroid material by means of solar radiation pressure, while Verrier et al. [27] studied the evolution of the halo family under the effects of radiation pressure in the Sun-Earth system.

In this paper, the equilibrium points as well as the periodic orbits around them are studied in the framework of the photogravitational Chermnykh’s restricted three-body problem with the aim to investigate the changes which may result in these basic dynamical features due to the three parameters of this model-problem (the angular velocity and the radiation factors and of the primaries). In particular, we study the existence and location of the planar equilibrium points as well as their linear stability. An interesting result is that there are combinations of the parameters of the problem for which the collinear equilibrium points may be stable. This is established by certain conditions which are determined when the triangular equilibrium points fall on the Ox-axis and transfer their stability to the collinear points. Several stability regions are given to show this result. For all cases of stable collinear and triangular equilibrium points we compute the corresponding short and long period families. We also give emphasis on the computation of the long period family emanating from each equilateral equilibrium point, in the framework of the classical restricted three-body problem for the Sun-Jupiter system, showing numerically that the long period family constitutes a bridge connecting orbits of the short period family which are bifurcation points of these two families complementing thusly the results of Henrard [28].

Specifically, the paper is organized as follows. In Section 2, we introduce the equations of motion of the considered model-problem. In Section 3, the existence, location, and stability of planar equilibrium points are investigated while, in Section 4, periodic motions around stable collinear and triangular equilibrium points are studied. Finally, in Section 5, we summarize our results and conclude.

2. Equations of Motion

We consider a barycentric, rotating, and dimensionless coordinate system Oxy whose origin is at the center of mass of two bodies and , known as primaries, where is the mass parameter. The two bodies perform circular orbits around the center of mass with angular velocity and are always on the Ox-axis (Figure 1). They are also considered to be radiating sources with mass reduction factors due to the radiation , , respectively (Chernikov [29]). Thus, the motion of the third body of negligible mass is influenced by the gravitational force as well as the repulsive force of light pressure. This additional force of radiation has applications in stellar dynamics, where it is a main component of the influence of a star on dust particles (Kunitsyn and Tureshbaev [30] and Das et al. [31]). Note that radiation factors may also have negative values meaning that gravity is strengthened by radiation, a case which will not be considered in our study. The equations of motion of the massless body are given by the following system (Schuerman [32], Chermnykh [6], and Perdios and Ragos [8]):

where we have abbreviated

while is the potential function which is given by the following formula:

are the distances of the third body of infinitesimal mass from the primary and secondary body, respectively. The system admits the Jacobi integral:

where is the Jacobi constant.

Figure 1: The rotating coordinate system (Oxy) with angular velocity relative to the inertial frame (OXY). The two primaries and the massless body are denoted by and , respectively.

3. Equilibrium Points and Their Stability

3.1. Collinear Equilibria

The Eulerian equilibrium points lie on the Ox-axis and may be determined from system (1) by setting the components of velocity and acceleration equal to zero as well as and by solving numerically for the nonlinear equation:

for certain values of the four parameters of the problem, , , , and . The case where , that is, the photogravitational version of the restricted three-body problem, has been studied extensively by Simmons et al. [18], where, for , , they show the existence of three collinear equilibrium points. The addition of the parameter of the angular velocity does not affect the total number of collinear equilibria and following the usual terminology existing in the literature, we name them as , . Their positions satisfy

where and are the locations of the larger and smaller primary, respectively. It is worth mentioning here the positions of these points, in the limiting case where the angular velocity parameter tends to infinity. So, supposing that , , and by elimination of the denominators in (6), we obtain that

which means that the only acceptable solution is ; that is, only the inner equilibrium point still exists and it is located at the origin, while the outer collinear points approach the primaries as the value of the angular velocity parameter increases and disappears when .

3.2. Triangular Equilibria

For zero velocity and acceleration components in the equations of motion (1) we obtain

and by solving this system for and , when , we are able to determine the positions of the Lagrangian equilibrium points. System (9) yields the following expressions:

and by substituting them into (4) we obtain, after some simple calculations, the exact positions of the triangular equilibrium points:

It ensues from (10) that since , these points exist for as well as for

or

where we set and The (limit) equality in (12), that is, , shows that the triangular equilibria have coalesced with the inner collinear equilibrium point while the equalities in (13), that is, or , indicate that the triangular points have coalesced with the outer collinear equilibrium points or , respectively. Summarizing, the necessary conditions for the existence of planar equilibrium points, with , which simultaneously have to be satisfied are

We mention here that in the limiting case where the angular velocity parameter tends to infinity these points do not exist as this may be directly seen from (10). Note also that for inequality (12) is in agreement with Schuerman [32] while (13) cannot be established in the photogravitational version since the triangular equilibria do not fall on any outer collinear point.

3.3. Stability of Equilibria

In order to study their stability we transfer the origin at an equilibrium point , :

and we expand the equations of motion (1) into first-order terms with respect to and obtaining the linearized system:

where the matrix of variations is

with

where we have abbreviated

while , are given by (2) and , by (4) evaluated at the equilibrium point . The characteristic equation is quadratic with respect to and is given by

For the collinear equilibrium points , and stability occurs when all roots of the characteristic equation for are imaginary. Therefore, for stability the following three necessary and sufficient conditions must be fulfilled simultaneously:

which ensure that the characteristic equation (20) has two negative real roots ; that is, the characteristic equation for has four pure imaginary roots. The characteristic equation (20) may be written in the following form:

and inequalities (21) lead to the double inequality for :

Also, by solving (6) for we get for any collinear equilibrium point that

Let us consider now the case of the outer collinear equilibrium point located in the open interval ; that is, From (24) and the right inequality of (23) we easily obtain that

and using this inequality we have from (6) that

Following the same procedure, we obtain from (24) and the left inequality of (23) and (6) that

Therefore, in order to have stability for the collinear equilibrium point the following inequalities have to be satisfied:

where we recall that we have set and Using a similar analysis for the outer collinear equilibrium point located in the interval we obtain that for stability the corresponding inequalities take the following form:

while for the inner collinear equilibrium point located in the interval we get them in the following form:

To study the stability of the triangular equilibrium points it suffices to consider one of them, for example, the equilibrium point , the position of which is given by (11) for the positive sign of the coordinate. So, by substituting its exact position into the characteristic equation (20) we obtain

Note here that, in the case of the triangular equilibrium points, the quantities in (20) are not zero. Therefore, in order for to be stable, inequalities (21) take now the following form:

which depend only on the four parameters of the problem. The latter of the above two inequalities arises from the discriminant of the characteristic equation (31) and may be written as a second degree polynomial with respect to :

where

When the discriminant is equal to zero, the characteristic equation has a double negative root which means that the triangular points are unstable. Thus, the roots of (33) provide the necessary condition for the mass parameter when stability alternates to instability and vice versa.

All the above analysis has been used to determine the stability regions of all equilibrium points of the model under consideration. So, in Figure 2 we present the stability regions in the plane for and for six different values of the mass parameter. Stability of the collinear equilibrium points is shown with light grey color, while dark grey color indicates stability of the triangular equilibrium points and white color corresponds to instability. Note here that the straight continuous lines correspond to the equalities of (12) and (13); that is, , indicating that the triangular equilibrium points have coalesced to a collinear equilibrium point. So, we can see from this figure that the unstable, in the case of the classical restricted three-body problem, collinear equilibrium points become stable when the stable triangular equilibrium points fall on the -axis. This means that the stability of the triangular equilibrium points is transferred to the collinear equilibrium points when and is lost from the dashed curves of the figure representing the rest of the inequalities of (28), (29), and (30). In particular, in the first frame of this figure which corresponds to equal masses of the primaries, we see in the top left, bottom left, and bottom right corners the stability regions (light grey) of the collinear equilibrium points , , and , respectively. For lower values of the mass parameter and for the same value of the angular velocity we observe, in the next frames, that the stability regions of all equilibria become larger and larger and especially in the last frame we see that the triangular equilibrium points are stable for all values of the parameters and , for which they exist.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: Stability regions of all equilibrium points in the plane for and , , , (from (a) to (c)) ,  , and   (from (d) to (f)).

In Figures 3, 4, and 5, we show, for comparison reasons, how the corresponding stability regions of the five equilibrium points become, for the specific values of the angular velocity parameter, , , and , by keeping the same values for the mass parameter. Note here that the case corresponds to the classical photogravitational restricted three-body problem.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: Stability regions of all equilibrium points in the plane for and , , , (from (a) to (c)) ,  , and   (from (d) to (f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: Stability regions of all equilibrium points in the plane for (i.e., the photogravitational restricted three-body problem) and , , , (from (a) to (c)) ,  , and   (from (d) to (f)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: Stability regions of all equilibrium points in the plane for and , , , (from (a) to (c)) ,  , and   (from (d) to (f)).

In Figure 3, we observe that for a larger value of the angular velocity parameter the stability regions of the collinear equilibrium points and are more restricted by the line (this line cuts both axes to larger values), whereas the stability region for the equilibrium point has been increased. In Figure 4,  , we see that cuts the axes , at the points and , respectively, which implies that the collinear points , are no longer stable (for this value of the angular velocity parameter the triangular equilibria coalesce only with the inner equilibrium point Finally, in Figure 5  , we also see that only the triangular and the inner collinear equilibrium points are stable but the stability region for the triangular points has been reduced substantially due to the limitation by the line

In conclusion, as the value of the parameter increases and up to values , the stability regions of the collinear equilibrium points and are restricted more and more by the curve , until these points lose completely their stability for . For the stability of the Lagrangian equilibrium points is transferred only to the inner Eulerian equilibrium point. In all the considered cases, the smaller the value of the mass parameter is, the larger the regions of the stability are.

In Figure 6, we show the stability regions of the triangular equilibrium points and the inner equilibrium point in the plane in the case where , for various values of the parameter . In this figure we do not obtain the stability regions of the outer Eulerian equilibrium points, since it illustrates the entry and exit points of the stability regions of Figures 2, 3, 4, and 5 when we move along the straight line for all values of the mass parameter. As in the previous figures, we depict the stability region of with light grey color and with dark grey we depict the corresponding region of the triangular points. As we see from Figure 6, when we move from the classical photogravitational problem, that is, , to lower values of the parameter , the stability region of the triangular equilibrium points increases, while the stability region of decreases due to (12); as we show (12) plays a fundamental role in the area of the corresponding regions for all equilibrium points and in this case takes the special form (see the horizontal continuous line, in all frames of this figure). Note here that when the second of inequalities (32), which arises from the discriminant of the characteristic equation (31), becomes

Solving this as a quadratic equation for we obtain its discriminant:

which is equal to zero when , that is, for the unique acceptable value of the mass parameter:

This is a remarkable result for this special case since the curve, which defines the stability region of the triangular equilibrium point, has a minimum value with respect to the mass parameter (see Figure 6) which is independent of the other parameters of the problem; that is, linear stability of triangular equilibria occurs at least for values of the mass parameter which are less than this unique value This value for the mass parameter has also been found by Schuerman [32] as a terminal value for the stability of the triangular equilibrium points of the photogravitational restricted three-body problem.

Figure 6: Stability regions of the triangular equilibrium points and the inner equilibrium point in the plane for , , and , respectively.

4. Periodic Solutions around Stable Equilibrium Points

The general solution of (16) is written in the following form:

where the coefficients , ,  , may be determined from the initial conditions. In case where all roots of the characteristic equation of the variational matrix (17) are purely imaginary, the motion is bounded and consists of two harmonic motions; the general solution depends on the eigenfrequencies , and can be written in the following form:

where , are the frequencies corresponding to the long and short periodic orbits, respectively. Therefore, the long periodic terms of (39) are those with the coefficients , ,  , while the short periodic terms are those with the coefficients , ,   Substituting (39) in (16) we obtain for the coefficients and ,  , that

Recall here that for the collinear equilibrium points and all the above formulae are simplified giving the following simple relations for the coefficients:

Now, for , we obtain from (39) the following:

In order to determine the long period solution we set ,  , in (39) obtaining

while for the short periodic solution we set , , leading to

To determine appropriate initial conditions for the numerical integration of the equations of motion (1), for the computation of the family of long periodic orbits, we consider a small orbital parameter so as and, without loss of generality, ; that is, we start our integration from the transferred Ox-axis to any stable equilibrium point. Hence, the first coefficient of (41) gives

and substituting it in the second coefficient of (41) we also get that

Accordingly, to determine suitable initial conditions for the numerical computation of the family of short periodic orbits we set and (we also choose to start our integration from the transferred Ox-axis), obtaining, thereupon, from the third coefficient of (40) the following:

and substituting it into the fourth coefficient of (40) one gets the following:

Therefore, by (43) and (46)–(49) we obtain the corresponding initial conditions:

where for we get the initial approximation of a long periodic orbit, while for the corresponding initial approximation of a short periodic orbit is obtained. For small values of the orbital parameter and fixed values of the parameters of the problem, expressions (50) may be used for the numerical determination of the first part of the families of short and long periodic orbits emanating from the stable Eulerian and Lagrangian equilibrium points, which are then continued numerically for the accurate computation of its members with a predetermined accuracy.

4.1. Numerical Results for (Classical Problem)

We will first focus on the determination of the long period family emanating from the Lagrangian equilibrium point (due to the symmetry of the problem with respect to the Ox-axis the same results may be obtained for the equilibrium point for the Sun-Jupiter case of the classical restricted three-body problem, that is, for the mass ratio , while the remaining parameters of the problem are equal to unity . Henrard [28] has given a schematic representation of the family evolution and stated that, in the Sun-Jupiter system, the family of long periodic orbits goes to a bifurcation with an orbit of the short period family travelled thirteen times and then it continues and bifurcates with another orbit of the short period family travelled fourteen times. This pattern continues to appear up to, as Henrard conjectured, the long period family’s termination on a family of homoclinic orbits. Almost recently, Henrard [33] has presented numerical results for the evolution of this family for and . Our aim here is to reproduce that schematic representation of Henrard [28] by precise numerical computations, so as to present the mechanism of the long period family’s evolution in more detail. For the vertical stability of the long and short period families, we may mention here the works by Perdios and Zagouras [34] and Hou and Liu [35].

In Figure 7, we present the evolution of this family by showing some of its members in the Oxy-plane. The long period family emanating from the triangular equilibrium point consists of nonsymmetric periodic orbits with respect to the Ox-axis, and its first members are small ovals around the corresponding equilibrium point (Figure 7(a)). Then, its periodic orbits create inner loops (Figure 7(b)) which gradually become tangled (Figures 7(c) and 7(d)). Afterwards, these loops start unraveling (Figures 7(e)–7(g)) and finally all loops coincide to a simple nonsymmetric periodic orbit, which is the bifurcation orbit of the short period family travelled thirteen times (Figure 7(h)). We denote this part of the long period family by

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 7: The evolution of the long period family by plotting sample periodic orbits. Note that the coordinate system has been transferred to the equilibrium point

In Figure 8, we show the remarkable next evolution of the family of long periodic orbits. We observe that as the family continues, its orbits unfold gradually (Figures 8(a) and 8(b)), while its next orbits are tangled (Figures 8(c) and 8(d)) which gradually unravel (Figure 8(e)) and in the sequel they form small ovals around the corresponding triangular equilibrium point (Figure 8(f)). Thereafter, the long periodic orbits start again to curl (Figures 8(g)–8(i)); then the tangled orbits progressively unravel (Figure 8(j)) and their loops tend to coincide to a simple nonsymmetric periodic orbit around (Figure 8(k)) until the long period family bifurcates with the short period family to an orbit travelled fourteen times (Figure 8(l)). We denote this part of the long period family by This pattern of orbits of branch also appears for the next part of the long period family, with the difference that now the new branch goes on an orbit of the short period family travelling it fifteen times. Accordingly, we denote this branch by (we do not present the corresponding orbits due to their similarity with the orbits of Figure 8). For the next branch of the long period family the same behavior has been observed and following the same notation we denote it by

(a)  

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

Figure 8: The evolution of the branch of the long period family by plotting sample periodic orbits. Note that the coordinate system has been transferred to the equilibrium point .

An alternative illustration of the evolution of the long period family is also shown in Figure 9, where we present the projection of its characteristic curve in the parametric plane of the space of initial conditions , where is the period. Note here that the ordinate is zero since we have transferred the coordinate system to and for the numerical computation of the corresponding family we always seek a crossing with the new Ox-axis. Also, in this figure, we denote by ,  , the bifurcation orbits of the short period family, where the subscript denotes the times that the orbit of the long family travels the corresponding orbit of the short family at the bifurcation point. Figures 9(b)–9(d) are appropriate magnifications of parts of Figure 9(a). As we see, in frame (a), the long period family emanates from the triangular equilibrium point and bifurcates with the short period family at the periodic orbit , which means that this orbit is travelled thirteen times by the corresponding periodic orbit of the long family (see Figure 9(b) for this junction). After this bifurcation point the long period family continues as and bifurcates the short period family at the orbit ; that is, this orbit is travelled fourteen times (see frame (c) for this junction). It then continues as and goes to the bifurcation orbit travelling it fifteen times (see frame (d) for this junction) and continues as This pattern of evolution of the long period family seems to continue but we have stopped our numerical computations at this point. Figure 9 shows that a connection of the bifurcation orbits ,  , of the short period family exists through the characteristic curve of the long period family. In Figure 10, we present this connection in the parametric plane For a detailed study on the evolution of the short period family for the Sun-Jupiter system we refer to Goodrich [36].

(a)

(b)

(c)

(d)

Figure 9: (a) Projection of the characteristic curve of the long period family in the plane. (b), (c), and (d) appropriate magnifications of certain parts of (a).

(a)

(b)

Figure 10: (a) Projection of the characteristic curve of the long period family in the plane. (b) Magnification of part of the left frame. In both frames, denotes the short period family.

In Table 1, we present initial conditions of representative periodic orbits of the long period family. In particular, for each nonsymmetric periodic orbit, we give its period , its abscissa on the Ox-axis of the transferred to the triangular equilibrium point coordinate system, the velocity components and , the value of the Jacobi constant , and the stability parameters and Note here that, according to Hénon [37, 38], a periodic orbit is stable if both and hold where

with , and , being the horizontal and vertical stability parameters, correspondingly. If only one of these inequalities hold, the orbit is considered horizontally or vertically stable, respectively (for the computation of the horizontal stability, see also the formulae given by Markellos [39]). In the case of symmetric periodic orbits and , obtaining thus and .

Table 1: Sample members of the long period family emanating from the triangular equilibrium point for the Sun-Jupiter system of the classical restricted three-body problem ( = 0.00095).

81.965981
0.2851750
0.1855400
−0.2593142
2.9747462
1.023708
0.951895

88.312060
0.1696324
0.1210896
−0.1151430
2.9964843
1.021820
0.929480

94.548562
0.3197985
0.1679201
−0.3395816
2.9546957
5.590865
0.932739

100.870333
0.2765249
0.1515870
−0.2739704
2.9733912
7.641533
0.918313

4.2. Numerical Results for , ,

We will now proceed to incorporate the other parameters of the problem and compute the corresponding long and short period families. The stability analysis of the equilibrium points, in Section 3, shows that additionally to the triangular equilibria these families may also be found at the collinear equilibrium points. Therefore, for our next results, the values of the parameters have been arbitrarily chosen but such that the corresponding equilibrium points are stable. We may mention here the work by Papadakis [40], in which periodic motions around the triangular equilibrium points of the Chermnykh’s problem for the Sun-Jupiter mass distribution were studied, as well as the corresponding short and long period families emanating from the inner equilibrium point for the critical value of the angular velocity parameter at which the triangular equilibrium points disappear by coalescing with transferring to it their stability. However, the addition in the Chermnykh’s problem of the new parameters and , corresponding to the new force of the radiation pressure, allows the existence of short and long period families at all collinear equilibrium points.

With regard to the triangular equilibrium points we have determined the corresponding families for the mass ratio , as in Section 4.1, while the values of the remaining parameters are , , and , respectively. The corresponding diagrams are shown in Figure 11. As we see from this figure, the evolution of the long period family follows the same mechanism as in the classical problem, but now the corresponding family goes to a bifurcation with an orbit of the short period family travelled twelve times and then it continues and bifurcates with another orbit of the short period family travelled thirteen times and so on. The asymmetric periodic orbits of the family fold and unfold as in the classical case (see Figures 7 and 8); therefore, we do not show their figures here. The characteristic curve of the short period family is shown in Figure 11(d) (black curve marked by ) and we note here that, as in the classical problem (see, e.g., Deprit et al. [41]), it terminates at a symmetric periodic orbit belonging to the Lyapunov family emanating from the equilibrium point (see Figure 12(a) for this bifurcation), at which it encounters with the corresponding short period family emanating from the other triangular equilibrium point Sample members of the short period family are presented in Figure 12(b) in which the bifurcation orbit is shown blue.

(a)  

(b)

(c)

(d)

Figure 11: Projection of the characteristic curve of the long period family for the parameter values , ,  , and (a) in the plane; (b), (c) magnifications of certain parts of the first frame showing the evolution in detail and (d) in the plane.

(a)

(b)

Figure 12: (a) The bifurcation of the short period family with the Lyapunov family emanating from the collinear equilibrium point . (b) Typical members of the short period family. The bifurcation orbit is shown blue.

In Tables 2 and 3, we give initial conditions of sample members of the long and short period families, respectively. Specifically, for each nonsymmetric periodic orbit of these two families, we give the full period , the abscissa on the new Ox-axis located at the triangular equilibrium point, the velocity components , , and the value of the Jacobi constant as well as the horizontal and vertical stability parameters and

Table 2: Sample members of the long period family emanating from the triangular equilibrium point for = 0.5, = 0.8, = 1.5, and = 0.00095.

50.433206
0.1294818
0.1565572
−0.2214219
2.4608218
0.817029
0.958225

54.583313
0.1250745
0.1515822
−0.2084973
2.4639318
0.634769
0.966236

58.728922
0.1157852
0.1421300
−0.1878419
2.4669604
1.401536
0.975416

Table 3: Sample members of the short period family emanating from the triangular equilibrium point for = 0.5, = 0.8, = 1.5, and = 0.00095.

4.206499
0.0220000
0.0239493
−0.0423981
2.4744822
1.894299
0.999646

4.205279
0.1520000
0.1352379
−0.3200771
2.4361058
1.848432
0.999648

4.197186
0.4156000
0.2395543
−1.0010119
2.0169264
1.607336
0.999336

4.187397
0.2010095
0.0259248
−1.2868579
0.9614449
−0.581246
0.999994

For the short and long period families emanating from the collinear equilibrium points when these are stable we present the corresponding diagrams in Figure 13, as examples of their evolution, for mass parameter and angular velocity and for pair values of the radiation factors and for , and , respectively. As we see, in Figure 13(a), the long period family emanating from the inner collinear equilibrium point terminates at this point while, in Figure 13(b), the long period family emanating from the outer collinear equilibrium point to the right of the smaller primary terminates at a bifurcation orbit which belongs to the short period family travelling it three times. The latter is also observed for the long period family emanating from the equilibrium point to the left of the larger primary (Figure 13(c)).

(a)

(b)

(c)

Figure 13: Projections, in the plane, of the characteristic curves of long (red color) and short (black color) period families emanating from the collinear equilibria for , and radiation factors: (a) and for , (b) and for , and (c) and for .

With regard to the short period families, at all collinear equilibrium points, they emanate from and terminate at their corresponding equilibria. In particular, the short period family emanating from starts with simple periodic orbits around the equilibrium point (see the ellipse in Figure 14(a)) but in its evolution a new vertical intersection with the Ox-axis appears (Figures 14(b) and 14(c)) continuing thus with double periodic orbits. Then, this new vertical intersection approaches the initial vertical intersection (Figure 14(d)) and at approximately the point these two vertical intersections coincide and then the family returns to the equilibrium point with the double period of the initial one; that is, it travels two times around The same pattern of evolution also appears for the short period family emanating from the collinear equilibrium point while for the corresponding family emanating from this starts with simple periodic orbits which become double (Figures 14(e)–14(g)). However, in its evolution the new vertical intersection does not approach the other vertical intersection but it disappears (Figure 14(h)); therefore, the family continues again with simple periodic orbits and it returns to the equilibrium point with the same period. Note here that all the corresponding short and long period families consist of periodic orbits which are symmetrical with respect to the Ox-axis.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 14: Sample members of the short period families emanating from the stable equilibrium points [(a)–(d)] and [(e)–(h)].

Initial conditions for representative periodic orbits of these families are presented in Table 4, where we give the half period, the abscissa on the Ox-axis, the velocity component ( at the half period of a symmetric periodic orbit), the value of the Jacobi constant , and the horizontal and vertical stability parameters and

Table 4: Sample members of the short and long period families emanating from the collinear equilibrium points for = 0.2 and = 0.5 and the following radiation factors: = 0.015 and = 0.05 for , = 0.9 and = 0.03 for , and = 0.03 and = 0.9 for .

6.721506
0.2232037
−0.0252754
0.1032007
−0.9880
0.9610

14.627908
0.1815907
0.0111466
0.1033552
0.7358
−0.1155

8.087888
1.7300000
−0.3327057
1.3965491
4.4041
−0.3814

19.595694
1.2147080
0.1069103
1.4042630
0.9273
0.8702

13.071411
−0.7991938
0.0721319
0.4596958
−0.8108
0.9720

15.920543
−0.6517328
−0.0378064
0.4589966
0.4090
−0.9978

5. Conclusions

Some basic dynamical features of the photogravitational Chermnykh’s restricted three-body problem were studied. In this problem, the primaries are considered to be radiating sources with mass reduction factors , , due to the radiation, and the massless body is moving in the orbital plane of the radiating primaries which is rotating with constant angular velocity around the center of mass.

With regard to the planar equilibrium points it was found that, as in the classical case, the problem admits five stationary solutions and, in particular, three collinear and two triangular. An interesting result is that the triangular equilibrium points may fall onto the Ox-axis, where they coalesce with the collinear equilibrium points. Specifically, they coalesce with the inner collinear point when and with the outer collinear points , when When these equations hold, the triangular equilibrium points transfer their stability to the corresponding collinear points. For all equilibria, stability criteria with respect to the parameters of the problem were established and several stability regions were plotted. An interesting result is that, in the special case of two primaries with equal mass reduction factors (due to the radiation), linear stability of the triangular equilibria occurs at least for values of the mass parameter which are less than the unique value , independently of the values of the other parameters of the problem. Comparing these results to previous results on the photogravitational restricted three-body problem and the corresponding Chermnykh’s problem we note that, for both of those problems, the triangular equilibrium points have been found to coalesce only with the inner collinear equilibrium point transferring to it their stability while for the classical case it is known that for each value of the mass parameter the triangular points are fixed at and

For all cases of stable collinear equilibrium points, the motion around them is bounded and is composed of two harmonic motions, giving rise to short and long period families, which consist of symmetric periodic orbits. Note that, for the corresponding Chermnykh’s or photogravitational problems and with respect to the specific equilibria, these families may be found only for the inner collinear equilibrium point, while for the classical problem, it is known that only the unstable Lyapunov families emanate from them. Certain examples of the short and long period families were given to show their existence and evolution. For the considered short period families of our examples we found that each one of them terminates at its corresponding equilibrium point.

Concerning the respective families emanating from the triangular equilibrium points, specific examples were drawn showing that their characteristic curves have a similar evolution with that of the corresponding families of the classical problem, which was also examined thoroughly. Specifically, the short period family terminates on the Lyapunov family emanating from the outer collinear equilibrium point to the left of the larger primary, while the long period family forms bridges connecting orbits of the short period family which are bifurcation points among these two families. The corresponding families consist of asymmetric periodic orbits. To the authors’ knowledge, related results for the Chermnykh or the photogravitational restricted three-body problem showing explicitly the evolution of the long period family do not exist in the literature.

A natural extension of the present study is to consider negative values for the radiation factors, a case for which at least new equilibrium points out of the orbital plane are allowed to exist. Another interesting case for study would be also to incorporate the Poynting-Robertson relativistic correction in the radiation force as well as to consider the belt gravitational potential as an additional term in the potential function.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

A. A. Nikaki acknowledges financial support under University of Patras “K. Karatheodory” research grant.

References

V. G. Szebehely, Theory of Orbits, the Restricted Problem of Three Bodies, Academic Press, 1967.

J. E. Marsden and S. D. Ross, “New methods in celestial mechanics and mission design,” Bulletin of the American Mathematical Society, vol. 43, pp. 43–73, 2005.

R. Prosmiti, S. C. Farantos, R. Guantes, F. Borondo, and R. M. Benito, “A periodic orbit analysis of the vibrationally highly excited LiNC/LiCN: a comparison with quantum mechanics,” The Journal of Chemical Physics, vol. 104, no. 8, pp. 2921–2931, 1996.

M. N. Vrahatis, A. E. Perdiou, V. S. Kalantonis et al., “Application of the characteristic bisection method for locating and computing periodic orbits in molecular systems,” Computer Physics Communications, vol. 138, no. 1, pp. 53–68, 2001.

M. M. Sano, “Dynamics starting from zero velocities in the classical Coulomb three-body problem,” Physical Review E, vol. 75, no. 2, Article ID 026203, 2007.

S. V. Chermnykh, “Stability of libration points in a gravitational field,” Leningradskii Universitet Vestnik Matematika Mekhanika Astronomiia, vol. 2, pp. 73–77, 1987.

K. Goździewski and A. J. Maciejewski, “Nonlinear stability of the Lagrangian libration points in the Chermnykh problem,” Celestial Mechanics & Dynamical Astronomy, vol. 70, no. 1, pp. 41–58, 1998.

E. A. Perdios and O. Ragos, “Asymptotic and periodic motion around collinear equilibria in Chermnykh's problem,” Astronomy and Astrophysics, vol. 414, no. 1, pp. 361–371, 2004.

K. E. Papadakis, “The 3D restricted three–body problem under angular velocity variation,” Astronomy & Astrophysics, vol. 425, no. 3, pp. 1133–1142, 2004.

A. E. Perdiou, A. A. Nikaki, and E. A. Perdios, “Periodic motions in the spatial Chermnykh restricted three-body problem,” Astrophysics and Space Science, vol. 345, no. 1, pp. 57–66, 2013.

I.-G. Jiang and L.-C. Yeh, “On the Chermnykh-like problems: I. The mass parameter μ = 0.5,” Astrophysics and Space Science, vol. 305, no. 4, pp. 341–348, 2006.

E. I. Abouelmagd, M. E. Awad, E. M. A. Elzayat, and I. A. Abbas, “Reduction the secular solution to periodic solution in the generalized restricted three-body problem,” Astrophysics and Space Science, vol. 350, no. 2, pp. 495–505, 2014.

J. Singh and O. Leke, “Periodic orbits in the Chermnykh–like restricted problem of oblate bodies with radiation,” Astrophysics and Space Science, vol. 350, no. 1, pp. 109–117, 2014.

J. Singh and O. Leke, “Motion in a modified Chermnykh's restricted three-body problem with oblateness,” Astrophysics and Space Science, vol. 350, no. 1, pp. 143–154, 2014.

V. V. Radzievskii, “The restricted three-body problem including radiation pressure,” Astronomicheskii Zhurnal, vol. 27, pp. 250–256, 1950.

T. J. Kalvouridis, “On some new aspects of the photo-gravitational Copenhagen problem,” Astrophysics and Space Science, vol. 317, no. 1-2, pp. 107–117, 2008.

R. Bewick, J. Sanchez, and C. R. McInnes, “The feasibility of using an L₁ positioned dust cloud as a method of space-based geoengineering,” Advances in Space Research, vol. 49, no. 7, pp. 1212–1228, 2012.

J. F. Simmons, A. J. McDonald, and J. C. Brown, “The restricted 3-body problem with radiation pressure,” Celestial Mechanics, vol. 35, no. 2, pp. 145–187, 1985.

K. Goździewski, A. J. Maciejewski, and Z. Niedzielska, “About stability of libration points in the restricted photogravitational three body problem,” Celestial Mechanics and Dynamical Astronomy, vol. 52, no. 2, pp. 195–201, 1991.

A. Elipe and M. Lara, “Periodic orbits in the restricted three body problem with radiation pressure,” Celestial Mechanics and Dynamical Astronomy, vol. 68, no. 1, pp. 1–11, 1997.

E. A. Perdios, S. S. Kanavos, and V. V. Markellos, “Bifurcations of plane to 3D periodic orbits in the photogravitational restricted three-body problem,” Astrophysics and Space Science, vol. 278, no. 4, pp. 407–413, 2001.

E. A. Perdios, “Critical symmetric periodic orbits in the photogravitational restricted three-body problem,” Astrophysics and Space Science, vol. 286, no. 3-4, pp. 501–513, 2003.

A. S. Zimovshikov and V. N. Tkhai, “Stability diagrams for a heterogeneous ensemble of particles at the collinear libration points of the photogravitational three-body problem,” Journal of Applied Mathematics and Mechanics, vol. 74, no. 2, pp. 158–163, 2010.

E. I. Abouelmagd, “The effect of photogravitational force and oblateness in the perturbed restricted three-body problem,” Astrophysics and Space Science, vol. 346, no. 1, pp. 51–69, 2013.

J. Singh and J. J. Taura, “Stability of triangular equilibrium points in the photogravitational restricted three-body problem with oblateness and potential from a belt,” Journal of Astrophysics and Astronomy, vol. 35, no. 2, pp. 107–119, 2014.

D. G. Yárnoz, J. P. S. Cuartielles, and C. R. McInnes, “Passive sorting of asteroid material using solar radiation pressure,” Journal of Guidance, Control, and Dynamics, vol. 37, pp. 1223–1235, 2014.

P. Verrier, T. Waters, and J. Sieber, “Evolution of the L₁ halo family in the radial solar sail circular restricted three-body problem,” Celestial Mechanics and Dynamical Astronomy, vol. 120, no. 4, pp. 373–400, 2014.

J. Henrard, “On Brown's conjecture,” Celestial Mechanics, vol. 31, no. 2, pp. 115–122, 1983.

Y. A. Chernikov, “The photogravitational restricted three-body problem,” Soviet Astronomy, vol. 14, pp. 176–181, 1970.

A. L. Kunitsyn and A. T. Tureshbaev, “On the collinear libration points in the photo-gravitational three-body problem,” Celestial Mechanics, vol. 35, no. 2, pp. 105–112, 1985.

M. K. Das, P. Narang, S. Mahajan, and M. Yuasa, “Effect of radiation on the stability of equilibrium points in the binary stellar systems: RW-Monocerotis, Krüger 60,” Astrophysics and Space Science, vol. 314, no. 4, pp. 261–274, 2008.

D. W. Schuerman, “The restricted three-body problem including radiation pressure,” The Astrophysical Journal, vol. 238, no. 1, pp. 337–342, 1980.

J. Henrard, “The web of periodic orbits at L ₄,” Celestial Mechanics and Dynamical Astronomy, vol. 83, no. 1–4, pp. 291–302, 2002.

E. A. Perdios and C. G. Zagouras, “Vertical stability of periodic solutions around the triangular equilibrium points,” Celestial Mechanics and Dynamical Astronomy, vol. 51, no. 1, pp. 75–81, 1991.

X. Y. Hou and L. Liu, “Vertical bifurcation families from the long and short period families around the equilateral equilibrium points,” Celestial Mechanics and Dynamical Astronomy, vol. 101, no. 3, pp. 309–320, 2008.

E. F. Goodrich, “Numerical determination of short-period Trojan orbits in the restricted three-body problem,” The Astronomical Journal, vol. 71, pp. 88–93, 1966.

M. Hénon, “Exploration numérique du problème restreint II—Masses égales, stabilité des orbites périodiques,” Annales d'Astrophysique, vol. 28, pp. 992–1007, 1965.

M. Hénon, “Vertical stability of periodic orbits in the restricted problem I. Equal masses,” Astronomy and Astrophysics, vol. 28, pp. 415–426, 1973.

V. V. Markellos, “On the stability parameters of periodic solutions,” Astrophysics and Space Science, vol. 43, no. 2, pp. 449–458, 1976.

K. E. Papadakis, “Motion around the triangular equilibrium points of the restricted three-body problem under angular velocity variation,” Astrophysics and Space Science, vol. 299, no. 2, pp. 129–148, 2005.

A. Deprit, H. Jacques, J. Palmore, J. F. Price, and D. H. Sadler, “The Trojan manifold in the system earth-moon,” Monthly Notices of the Royal Astronomical Society, vol. 137, no. 3, pp. 311–335, 1967.