We know that the four most popular methods of proving the existence of periodic orbits are:
(i) the method of analytic continuation,
(ii) the process of equating Fourier coefficients of equal frequencies,
(iii) the application of fixed point theorem given by Poincare,
(iv) the method of power series.
Giacaglia  used the method of analytic continuation to examine the existence of periodic orbits of collision in the Restricted Three-body Problem (R3BP). Bhatnagar  generalized the problem in elliptic case. The problem of Giacaglia  was further extended by Bhatnagar  in the R4BP by taking the primaries at the vertices of an equilateral triangle. With different perturbations like oblateness, triaxiality, photogravitation, Pointing-Robertson drag effects of the primaries, the existence of periodic orbits of collision in the R3BP and in the R4BP, have been studied by different authors in two and three-dimensional co-ordinate system during the period of last three decades of the 20th century but nobody established the proper mathematical model of the R4BP. Recently Ceccaroni and Biggs  has studied the autonomous coplanar CR4BP by taking the third primary of comparatively inferior mass at the triangular equilibrium point of R3BP and with an extension to low-thrust propulsion for application to the future science mission.
In present paper, we have proposed to study the existence of periodic orbits of first kind in the Autonomous Four-body Problem by the method of analytic continuation. By using Poincare surfaces of section (PSS), the conditions for the existence of periodic orbits given by Duboshin  have been confirmed. For collision case, we have applied the criterion given by Levi-Civitas   and it is satisfied by our model.
2. Equations of Motion
Let be the three massive bodies of masses respectively, where and the fourth body of mass be at. These bodies are moving in the same plane under some restrictions as follows:
The fourth body at of mass is assumed to be of infinitesimal mass not influencing the motion of but motions of is being influenced by the motions of. Further, we have assumed that the mass at is taken small enough, so that it can’t influence the motion of the dominating primaries and and it is placed at any one of the triangular libration points (Lagrangian Points) of the classical restricted three body problem. Since the third primary can’t influence the motions of and, so the centre of rotation of the system remains at the barycentre of two main primaries and. Also, it is supposed, all the primaries are moving in the same plane in circular orbits around the bary-centre of massive primaries and with the same angular velocity and the fourth body is moving under the gravitational field and plane of motion of three primaries then to check the nature of motion of infinitesimal mass.
Let the line joining and be taken as the x-axis and their mass centre (bary-centre) O, as the origin. Let the line through O and perpendicular to lying in the plane of motion of the primaries be taken as the y-axis. Let the positions of masses be and respectively. Let be the position vector of and be the displace-
Figure 1. Configuration of four-body problem.
ments of and relative to as shown in Figure 1, then
Let be the gravitational forces exerted on by the primaries respectively, then
where is the gravitational constant.
The total gravitational force acting on by the three primaries is given by
Let be the magnitude of angular velocity and be the unit vector normal to the plane of motion of the primaries, then.
The Equation of motion of the infinitesimal mass in synodic frame is
Since the synodic frame are revolving with constant angular velocity about the bary-centre, hence and thus Equation (4) reduces to
In cartesian form, the equations of motion of the infinitesimal mass in the gravitational field of three primaries, are given by
Also the linear velocity of the infinitesimal mass on its orbit; is given by
If are two components of, then from Equation (7),
If mass of the infinitesimal body is supposed to be unity, then the kinetic energy of the infinitesimal mass is given by
Let be the momenta corresponding to the co-ordinates respectively, then
Combination of Equations ((9) and (10)) yields
The gravitational potential of the body of mass at any point of outside it, is given by
then, total gravitational potential at due to three primaries is given by
The Hamiltonian of the infinitesimal body of unit mass is given by
Let be the reduced mass of the second primary and be the reduced mass of the third primary, then from the definition of reduced mass, we have
The coordinates of are given by
Clearly, which implies that forms an equilateral triangle of sides of unit length. We know that is very small in comparison of masses of the other two primaries, so we can choose as the order of i.e.,. Now choosing unit of time in such a manner that and and taking, then the Hamilton canonical equations of motion of the infinitesimal body are given by
is the reduced Hamiltonian corresponding to canonically conjugate variables and.
3. Regularization at the Singularity
In our Hamiltonian given in Equation (15), there are three singularities. To examine the existence of periodic orbits of collision with the first primary, we have to eliminate the singularity. For this, let us define an extended generating function S given by
where is the momenta associated with new co-ordinate.
Thus the Hamiltonian given in Equation (15), can be written in terms of new variables, as
Let us introduce pseudo time by the differential equation
Thus the regularized Hamilton-canonical equations of motion of the infinitesimal body corresponding to the Hamiltonian, are given by
where the regularized Hamiltonian is given by
Let us write, then
4. Generating Solution (i.e., Solutions When)
For generating solutions, we shall choose for our Hamiltonian function, so in order to solve the Hamilton-Jacobi equation associated with, let us write
where is an arbitrary constant.
Since is not involved explicitly in: hence by using Equation (27) in Equation (25), the Hamilton-Jacobi equation may be written as
Putting, then the Equation (27) becomes
It may be noted that this differential equation is exactly the same as in Giacaglia  and Bhatnagar   and therefore the solution of Equation (29) can be written by the method of separation of variables, as
where is an arbitrary constant.
Let us introduce a new quantity by then from Equation (30), we get
Combination of Equations (29) and (30) yields
where is the smaller root of the roots of the equation.
From Equation (33), we conclude that for general solution; we need only two arbitrary constants as and. Therefore the solution of Equation (30) may be regarded as a general solution.
Let us introduce the parameters by the relations
where is the semi-major axis, is the eccentricity and is the latus-rec- tum of the elliptic orbit of the infinitesimal body.
It may be noted that for and is the other root of the equation.
We introduce a parameter by the relation
From Equations (33), (35) and (36), we get
Again from Equation (25)
Thus the equations of motion associated with are given as
where denotes the differentiation with respect to.
Now from we get.
and [Using Equation (38)]
Thus from the above relations, we have
From Equation (32), we get
From Equation (30),
[Using Equation (34)]
If we take and as arbitrary constants, the solutions may be written as
From the second equation of system (41), we get the argument as
hence for the problem generated by Hamiltonian (regularized two-body problem in rotating co-ordinate system), we have
The variables can now be expressed in terms of the canonical elements for, as
where is given by the first equation of system (42).
When and, then
where is given by the second equation of system (42).
The original synodic cartesian co-ordinates in a non-uniformly rotating system are obtained from Equations (18) and (20), when, as
The sidereal cartesian co-ordinates are obtained by considering the transformations
where is given by
where is a constant.
In terms of canonical variables introduced, the complete Hamiltonian may be written as, where can be obtained from Equation (26) after changing into canonical variables.
The equations of motion for the complete Hamiltonian are
Equation (49) forms the basis of a general perturbation theory for the present problem. The solution described by Equations ((44) and (45)) and is periodic if and g have commensurable frequencies, i.e., if
where and are integers.
The periods of are and respectively, so that in case of commensurability, the period of the solution is or.
5. Existence of Periodic Orbits When
Here we shall follow the method given by Chaudhary  to prove the existence of periodic orbits. Let then from Equation (44), when, we have
Integrating these equations with respect to, we get
These are the generating solutions of two-body problems. The generating solution will be periodic with the period, if
when are integers, so that are commensurable.
Following Poincare  , the general solution in the neighbourhood of the generating solution, may be given as
where is the new independent variable given by
The necessary and sufficient conditions for the existence of periodic solution are
Restricting our solution only up to the first order infinitesimals, the equations of motion may be written as
Expanding in ascending powers of, Equation (54) may be written as
Rejecting the second order term, integrating and putting the value of in Equation (51), we get
The Equation (55) gives
Equation (45) gives
By solving the Equations (54)-(57), we can find the values of, as analytic function of, reducing to zero with, if the conditions for periodic orbits given by Duboshin  are satisfied i.e.,
(i), and (58)
(ii), together (59)
where is the zero degree terms of given in Equation (26).
From Equation (43),
From the Equation (26)
Taking only zero order terms i.e., for
Now from equations of system (52)
and from Equation (63)
Here if either or and if either or.
But and don’t imply each other, so is only the case for which and will be simultaneously zero.
Now choosing suitably, then (say)
As so from Equation (63), we have
Using Equation (65), we get
Thus the conditions for the existence of periodic orbits given by Duboshin  are satisfied i.e., in the region of motion of the infinitesimal body, periodic orbits exist.
6. Poincare Surfaces of Section (PSS)
In this previous section, we have shown that Duboshin’s condition  for the existence of periodic orbits when, are satisfied. So to justify the mathematical model given in Equations (58)-(60), we have applied the method of Poincare surfaces of section (PSS) to the reduced equations of motion
together with the Jacobi Integral
To study the motion of the infinitesimal body by PSS, it is necessary to know its position and velocity which correspond to a point in four- dimensional phase space. By defining a plane, in the resulting three- dimensional space, the values of and can be plotted. Every time the particle has, whenever the trajectory intersects the plane in a particular direction say.
The techniques of PSS suggest to determine the regular or chaotic nature of the trajectories. If there are smooth, well-defined island then the trajectory is likely to be regular and the islands correspond to oscillation around a periodic orbit. As the curves shrink down to a point, the points represent a periodic orbit as per Kolmogorov-Arnold-Moser (KAM) theory. Any fuzzy distribution of points in surfaces of section, implies that trajectory is chaotic. In Figure 2, for and Poincare surfaces of section have been plotted in which atleast seven points are visible towards which the regular trajectories shrink, hence by KAM theory, periodic orbits exists. Again Figure 3 represents a Poincare surfaces of section for and in which atleast nine points are visible towards which the regular trajectories shrink, so we can say that the periodic orbits exist in the region of motion of infinitesimal mass. Other than the neighbourhood of these points, the quasi-periodic and chaotic regions are seen in the PSS. In Figure 4, in PSS for and, atleast ten shrinking regions of regular curves to a point are visible, i.e., that the degree of existence of periodic orbits increases in the region of motion of the infinitesimal mass. Thus by increasing the values of the Jacobi’s constant, the chances of existence of periodic orbits increase. Thus the Duboshin conditions and PSS both confirms the existence of periodic orbits when. In Figure 5, regions plot of ZVC (Zero Velocity Curves) for is shown, in which central white circle represents regions of no motion and coloured annulus represents the regions of periodic orbits. Figure 6 depicts the contour plot of ZVC for.
7. Periodic Orbits of Collision When
Levi-Civita   proved that the invariant relation for collision orbits can be analytically continued from the one that corresponds to the problem of two bo-
Figure 2. Poincare Surface of Section for.
Figure 3. Poincare Surface of Section for.
Figure 4. Poincare Surface of Section for.
Figure 5. Region Plot of ZVCs for.
Figure 6. Contour Plot of ZVCs for.
dies. Bhatnagar   has developed this as. For the present paper when, the condition must be
for sufficiently small and.
Further, he has proved that, in particular, such relation is uniform integral of the differential equation of motion along any collision orbit. He has also proved this integral is a power series in terms of the distance from the origin and the series is convergent through the radius of convergence is generally small. In section (5), we have shown that periodicity is conserved by analytic continuation. Let us show that the condition of collision is also conserved by analytic continuation.
Figure 7 shows the geometrical configuration of collision orbits. In order to show the validity of that continuation, we shall consider orbits corresponding to the case when. When, the orbits starts as an ejection
from the origin and return to it after. Bhatnagar   and Levi-Civita   finds the condition for collision as
Figure 7. Geometrical configuration of collision orbits.
Therefore, the condition of Equation (69) became,
Thus from Equations ((70) and (71))
Here the Equation (71) corresponds to the Equation (68), so it is easy to say that the collision orbits exist.
8. Discussions and Conclusion
In section 1 of this paper, historical background has been sketched with original and previous contributions. In section 2, the equations of motion of the infinitesimal mass moving under the gravitational field of the three primaries situated at the vertices of an equilateral triangle taken by Ceccaroni and Biggs  . In this the reduced Hamiltonian has been derived for regularization in the next section 3. In this section, the regularized Hamiltonian has been established. In section 4, generating solutions have been found by taking as the corresponding Hamiltonian. In this section generating solution forms a basis for general solution by the process of analytic continuation. In section 5, using Duboshin’s criterion  for the existence of periodic orbits has been satisfied following the method of Choudhary  . For confirmation of the existence of periodic orbits in section 4, we have analyzed PSS in section 6 and justified that the region of motion regular trajectory shrinking towards a point represents the periodic orbits and other region of the PSS represents quasi-periodic and chaotic belt in the region of motion. In section 7, the periodic orbits of collision for have been shown. In our discussion, we have shown that our condition of collision orbit has a resemblance with the condition given by Bhatnagar   .
Bary-Centre: It is the center of mass of two or more bodies that are orbiting each other, or the point around which they both orbit.
Synodic Co-ordinate System: The co-ordinate system, in which the xy-plane rotates in the positive direction with an angular velocity equal to that of the common velocity of one primary with respect to the other keeping the origin fixed, is called synodic co-ordinate system.
Reduced Mass: Mass ratio of the smaller primary to the total mass of the primaries or the non-dimensional mass of the smaller primary is known as reduced mass of the smaller primary.
Regularization: The process of elimination of the singularity from the force function is known as regularization.
Submit or recommend next manuscript to SCIRP and we will provide best service for you:
Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.
A wide selection of journals (inclusive of 9 subjects, more than 200 journals)
Providing 24-hour high-quality service
User-friendly online submission system
Fair and swift peer-review system
Efficient typesetting and proofreading procedure
Display of the result of downloads and visits, as well as the number of cited articles
Maximum dissemination of your research work
Submit your manuscript at: http://papersubmission.scirp.org/
Or contact email@example.com
 Ceccaroni, M. and Biggs, J. (2012) Low-Thrust Propulsion in a Coplanar Circular Restricted Four body Problem. Celestial Mechanics and Dynamical Astronomy, 112, 191-219.
 Choudhry, R.K. (1966) Existence of Periodic Orbits of the Third Kind in the Elliptic Restricted Problem of the Three Bodies and the Stability of the Generating Solution. Proceedings of the National Academy of Sciences India Section A, 36, 249-264.