Application of the stroboscopic method to the stellar three-body problem

The stellar three-body problem has been approached by directly integrating the equations of motion in the orbital elements. The problem is set up in a barycentric chain with the orbits as perturbed ellipses.The integration was performed using the semi-analytical stroboscopic method, which is particularly useful for solving differential systems that depend on several slow variables and one fast, angular variable. Perturbation theory is applied and the solution for a particular order is obtained by way of successive approximations on the fundamental period of the fast variable.     

Numerical models for the study of motion of lunar satellites

The present study deals with numerical modeling of the elliptic restricted three-body problem as well as of the perturbed elliptic restricted three-body (Earth-Moon-Satellite) problem by a fourth body (Sun). Two numerical algorithms are established and investigated. The first is based on the method of the series solution of the differential equations and the second is based on a 5th-order Runge-Kutta method. The applications concern the solution of the equations and integrals of motion of the circular and elliptical restricted three-body problem as well as the search for periodic orbits of the natural satellites of the Moon in the Earth-Moon system in both cases in which the Moon describes circular or elliptical orbit around the Earth before the perturbations induced by the Sun. After the introduction of the perturbations in the Earth-Moon-Satellite system the motions of the Moon and the Satellite are studied with the same initial conditions which give periodic orbits for the unperturbed elliptic problem.     

Periodic solutions of the singly averaged hill problem

Based on the ideas of Lyapunov’s method, we construct a family of symmetric periodic solutions of the Hill problem averaged over the motion of a zero-mass point (a satellite). The low eccentricity of the satellite orbit and the sine of its inclination to the plane of motion of the perturbing body are parameters of the family. We compare the analytical solution with numerical solutions of the averaged evolutionary system and the rigorous (nonaveraged) equations of the restricted circular three-body problem.     

A note on independent variables for restricted three-body problems

In studies of the elliptic restricted three-body problem, the true anomaly of the motion of the primaries is often used as the independent variable. The equations of motion then show invariancy in form from the circular case. It is of interest whether other independent variables exist, such that the invariant form of the equations is maintained. It is found that true anomaly is the only such variable.     

On the motion around the centre of mass of a viscously elastic sphere in the restricted elliptic three-body problem

In the restricted three-body problem we consider the motion of a viscously elastic sphere (planet) with its centre of mass moving in a conditionally-periodic orbit. The approximate equations describing the rotational motion of the sphere in terms of the Andoyer variables are obtained by the method of the separation of motion and averaging; the evolution of the motion is also analysed.     

The circular restricted three-body problem in curvilinear coordinates

We derive the exact equations of motion for the circular restricted three-body problem in cylindrical curvilinear coordinates together with a number of useful analytical relations linking curvilinear coordinates and classical orbital elements. The equations of motion can be seen as a generalization of Hill’s problem after including all neglected nonlinear terms. As an application of the method, we obtain a new expression for the averaged third-body disturbing function including eccentricity and inclination terms. We employ the latter to study the dynamics of the guiding center for the problem of circular coorbital motion providing an extension of some results in the literature.     

A global regularisation of the gravitationalN-body problem

The work of Aarseth and Zare (1974) is extended to provide aglobal regularisation of the classical gravitational three-body problem: by transformation of the variables in a way that does not depend on the particular configuration, we obtain equations of motion which are regular with respect to collisions between any pair of particles. The only cases excepted are those in which collisions between more than one pair occur simultaneously and those in which at least one of the masses vanishes. However, by means of the same principles the restricted problem is regularised globally if collisions between the two primaries are excluded. Results of numerical tests are summarised, and the theory is generalised to provide global regularisations, first, for perturbed three-body motion and, second, for theN-body problem. A way of increasing the number of degrees of freedom of a dynamical system is central to the method, and is the subject of an Appendix.     

Studies in the application of recurrence relations to special perturbation methods

Using the rectangular equations of motion for the restricted three-body problem a comparison is made of the integration of these equations by the Encke method and by a set of perturbational equations. Each set of differential equations is integrated using Taylor series expansions where the coefficients of the powers of time are determined by recurrence relations. It is shown that for very small perturbations the use of the perturbational equations is more efficient than the use of the Encke method. A discussion is also given of when Cowell's method is more efficient than either of these techniques.     

Out-of-plane motion about libration points: Nonlinearity and eccentricity effects

Out-of-plane motion about libration points is studied within the framework of the elliptic restricted three-body problem. Nonlinear motion in the circular restricted problem is given to third order in the out-of-plane amplitudeA _z by Jacobi elliptic functions. Linear motion in the elliptic problem is studied using Mathieu's and Hill's equations. Additional terms needed for a complete third-order theory are found using Lindsted's method. This theory is constructed for the case of collinear libration points; for the case of triangular points, a third-order nonlinear solution is given separately in terms of Jacobi elliptic functions.     

On the periodically evolving orbits in the singly averaged Hill problem

We continue to analyze the periodic solutions of the singly averaged Hill problem. We have numerically constructed the families of solutions that correspond to periodically evolving satellite orbits for arbitrary initial values of their eccentricities and inclinations to the plane of motion of the perturbing body. The solutions obtained are compared with the numerical solutions of the rigorous (nonaveraged) equations of the restricted circular three-body problem. In particular, we have constructed a periodically evolving orbit for which the well-known Lidov-Kozai mechanism manifests itself, just as in the doubly averaged problem.     

Conditionally periodic solutions near the libration points of the restricted circular three-body problem

Families of conditionally periodic solutions have been found by a slightly modified Lyapunov method of determining periodic solutions near the libration points of the restricted three-body problem. When the frequencies of free oscillations are commensurable, the solutions found are transformed into planar or spatial periodic solutions. The results are confirmed by numerically integrating the starting nonlinear differential equations of motion.     

Comment on the paper “On the triangular libration points in photogravitational restricted three-body problem with variable mass” by Zhang, M.J., et al.

In a recent paper, published in Astrophys. Space Sci. (337:107, 2012) (hereafter paper ZZX) and entitled “On the triangular libration points in photogravitational restricted three-body problem with variable mass”, the authors study the location and stability of the generalized Lagrange libration points L ₄ and L ₅. However their study is flawed in two aspects. First they fail to write correctly the equations of motion of the variable mass problem. Second they attribute a variable mass to the third body of the restricted three-body model, a fact that is not compatible with the assumptions used in deriving the mathematical formulation of this model.     

Explicit evolution relations with orbital elements for eccentric,inclined, elliptic and hyperbolic restricted few-body problems

Planetary, stellar and galactic physics often rely on the general restricted gravitational $N$ -body problem to model the motion of a small-mass object under the influence of much more massive objects. Here, I formulate the general restricted problem entirely and specifically in terms of the commonly used orbital elements of semimajor axis, eccentricity, inclination, longitude of ascending node, argument of pericentre, and true anomaly, without any assumptions about their magnitudes. I derive the equations of motion in the general, unaveraged case, as well as specific cases, with respect to both a bodycentric and barycentric origin. I then reduce the equations to three-body systems, and present compact singly- and doubly-averaged expressions which can be readily applied to systems of interest. This method recovers classic Lidov–Kozai and Laplace–Lagrange theory in the test particle limit to any order, but with fewer assumptions, and reveals a complete analytic solution for the averaged planetary pericentre precession in coplanar circular circumbinary systems to at least the first three nonzero orders in semimajor axis ratio. Finally, I show how the unaveraged equations may be used to express resonant angle evolution in an explicit manner that is not subject to expansions of eccentricity and inclination about small nor any other values.     

Hill stability of the planar three-body problem: General and restricted cases

The equations of motion of the planar three-body problem split into two parts, called an external part and an internal part. When the third mass approaches zero, the first part tends to the equations of the Kepler motion of the primaries and the second part to the equations of motion of the restricted problem.We discuss the Hill stability from these equations of motion and the energy integral. In particular, the Jacobi integral for the circular restricted problem is seen as an infinitesimal-mass-order term of the Sundman function in this context.     

About regions of convergence of expansions of differential equations of the three-dimensional restricted three-body problem in the vicinity of collinear libration points

The analysis of regions of convergence of expansions of the right-hand sides of the differential equations of motion in the vicinity of the collinear libration points in the circular restricted three-body problem in powers of coordinates of the infinitesimal body due to F. R. Moulton (Moulton, 1920) is shown to be erroneous and his results are corrected. The generalisation of Moulton's results to analogous expansions of the equations in the elliptic problem of three bodies made by R. W. Farquhar (Farquhar, 1968) is shown to be groundless.     

The stellar three-body problem

Two applications of von Zeipel's method to the stellar three-body problem eliminate the short period terms and establish two new integrals of the motion beyond the classical integrals. The remaining time averaged problem with only the second order Hamiltonian has one additional integral and can be solved. The motion with the third order averaged Hamiltonian included is more complex, in that there may be additional resonances, and the additional integral does not exist in all cases.     

A keplerian method to estimate perturbations in the restricted three-body problem

We develop a new and fast method to estimate perturbations by a planet on cometary orbits. This method allows us to identify accurately the cases of large perturbations in a set of fictitious orbits. Hence, it can be used in constructing perturbation samples for Monte Carlo simulations in order to maximize the amount of information. Furthermore, the estimated perturbations are found to yield a good approximation to the real perturbation sample. This is shown by a comparison of the perturbations obtained by the new estimator with the results of numerical integration of regularized equations of motion for the same orbits in the same dynamical model: the three-dimensional elliptic restricted three-body problem (Sun-Jupiter-comet).     

Stable Manifolds and Homoclinic Points Near Resonances in the Restricted Three-Body Problem

The restricted three-body problem describes the motion of a massless particle under the influence of two primaries of masses 1− μ and μ that circle each other with period equal to 2π. For small μ, a resonant periodic motion of the massless particle in the rotating frame can be described by relatively prime integers p and q, if its period around the heavier primary is approximately 2π p/q, and by its approximate eccentricity e. We give a method for the formal development of the stable and unstable manifolds associated with these resonant motions. We prove the validity of this formal development and the existence of homoclinic points in the resonant region. In the study of the Kirkwood gaps in the asteroid belt, the separatrices of the averaged equations of the restricted three-body problem are commonly used to derive analytical approximations to the boundaries of the resonances. We use the unaveraged equations to find values of asteroid eccentricity below which these approximations will not hold for the Kirkwood gaps with q/p equal to 2/1, 7/3, 5/2, 3/1, and 4/1. Another application is to the existence of asymmetric librations in the exterior resonances. We give values of asteroid eccentricity below which asymmetric librations will not exist for the 1/7, 1/6, 1/5, 1/4, 1/3, and 1/2 resonances for any μ however small. But if the eccentricity exceeds these thresholds, asymmetric librations will exist for μ small enough in the unaveraged restricted three-body problem.     

Periodic orbits of the planetary type and their stability

A review is presented of periodic orbits of the planetary type in the general three-body problem and fourbody problem and the restricted circular and elliptic tnreebody problem. These correspond to planetary systems with one Sun and two or three planets (or a planet and its satellites), the motion of asteoids and also planetary systems with two Suns. The factors which affect the stability of the above configurations are studied in connection with resonance or additional perturbations. Finally, the correspondence of the periodic orbits in the restricted three-body problem with the fixed points obtained by the method of averaging or the method of surface of section is indicated.     

The perturbations of the orbital elements of Trojan asteriods

An asymptotic solution for the cylindrical coordinates of Trojan asteroids is derived by using a three-variable expansion method in the elliptic restricted three-body problem. The perturbations of the orbital elements are obtained from this solution by applying the formulas of the two-body problem. The main perturbations of the mean motion are studied in detail.