The rectilinear three-body problem

The rectilinear equal-mass and unequal-mass three-body problems are considered. The first part of the paper is a review that covers the following items: regularization of the equations of motion, integrable cases, triple collisions and their vicinities, escapes, periodic orbits and their stability, chaos and regularity of motions. The second part contains the results of our numerical simulations in this problem. A classification of orbits in correspondence with the following evolution scenarios is suggested: ejections, escapes, conditional escapes (long ejections), periodic orbits, quasi-stable long-lived systems in the vicinity of stable periodic orbits, and triple collisions. Homothetic solutions ending by triple collisions and their dependence on initial parameters are found. We study how the ejection length changes in response to the variation of the triple approach parameters. Regions of initial conditions are outlined in which escapes occur after a definite number of triple approaches or a definite time. In the vicinity of a stable Schubart periodic orbit, we reveal a region of initial parameters that corresponds to trajectories with finite motions. The regular and chaotic structure of the manifold of orbits is mostly defined by this periodic orbit. We have studied the phase space structure via Poincaré sections. Using these sections and symbolic dynamics, we study the fine structure of the region of initial conditions, in particular the chaotic scattering region.     

Chaotic Scattering in the Restricted Three‐Body Problem II. Small mass parameters

We study the scattering motion of the planar restricted three‐body problem for small mass parameters μ. We consider the symmetric periodic orbits of this system with μ = 0 that collide with the singularity together with the circular and parabolic solutions of the Kepler problem. These divide the parameter space in a natural way and characterize the main features of the scattering problem for small non‐vanishing μ. Indeed, continuation of these orbits yields the primitive periodic orbits of the system for small μ. For different regions of the parameter space, we present scattering functions and discuss the structure of the chaotic saddle. We show that for μ < μc and any Jacobi integral there exist departures from hyperbolicity due to regions of stable motion in phase space. Numerical bounds for μc are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

Exploiting Unstable Periodic Orbits of a Chaotic Invariant Set for Spacecraft Control

The purpose of this work is to show that chaos control techniques (OGY, in special) can be used to efficiently keep a spacecraft around another body performing elaborate orbits. We consider a satellite and a spacecraft moving initially in coplanar and circular orbits, with slightly different radii, around a heavy central planet. The spacecraft, which is the inner body, has a slightly larger angular velocity than the satellite so that, after some time, they eventually go to a situation in which the distance between them becomes sufficiently small, so that they start to interact with one another. This situation is called as an encounter. In previous work we have shown that this scenario is a typical situation of a chaotic scattering for some well-defined range of parameters. Considering this scenario, we first show how it is possible to find the unstable periodic orbits that are located in the chaotic invariant set. From the set of unstable periodic orbits, we select the ones that can be combined to provide the desired elaborate orbit. Then, chaos control technique based on the OGY method is used to keep the spacecraft in the desired orbit. Finally, we analyze the results and make considerations regarding a realistic scenario of space exploration.     

The Symmetrical One-dimensional Newtonian Four-body Problem: A Numerical Investigation

Numerical simulations of the one-dimensional Newtonian four-body problem have been conducted for the special case in which the bodies are distributed symmetrically about the centre of mass. Simulations show a great similarity between this problem and the one-dimensional Newtonian three-body problem. As in that problem the orbits can be divided into three different categories which form well-defined regions on a Poincaré section: there is a region of quasiperiodic orbits about a Schubart-like periodic orbit, there is a region of fast-scattering encounters and in between these two regions there is a chaotic scattering region. The Schubart-like periodic orbit's stability to perturbation is studied. It is apparently stable in one-dimension but is unstable in three-dimensions.This revised version was published online in October 2005 with corrections to the Cover Date.     

Periodic orbits and bifurcations in the Sitnikov four-body problem

We study the existence, linear stability and bifurcations of what we call the Sitnikov family of straight line periodic orbits in the case of the restricted four-body problem, where the three equal mass primary bodies are rotating on a circle and the fourth (small body) is moving in the direction vertical to the center mass of the other three. In contrast to the restricted three-body Sitnikov problem, where the Sitnikov family has infinitely many stability intervals (hence infinitely many Sitnikov critical orbits), as the “family parameter” ż₀ varies within a finite interval (while z ₀ tends to infinity), in the four-body problem this family has only one stability interval and only twelve 3-dimensional (3D) families of symmetric periodic orbits exist which bifurcate from twelve corresponding critical Sitnikov periodic orbits. We also calculate the evolution of the characteristic curves of these 3D branch-families and determine their stability. More importantly, we study the phase space dynamics in the vicinity of these orbits in two ways: First, we use the SALI index to investigate the extent of bounded motion of the small particle off the z-axis along its interval of stable Sitnikov orbits, and secondly, through suitably chosen Poincaré maps, we chart the motion near one of the 3D families of plane-symmetric periodic orbits. Our study reveals in both cases a fascinating structure of ordered motion surrounded by “sticky” and chaotic orbits as well as orbits which rapidly escape to infinity.     

Chaotic quasi-collision trajectories in the 3-centre problem

We study a particular kind of chaotic dynamics for the planar 3-centre problem on small negative energy level sets. We know that chaotic motions exist, if we make the assumption that one of the centres is far away from the other two (see Bolotin and Negrini, J Differ Equ 190:539–558, 2003): this result has been obtained by the use of the Poincaré-Melnikov theory. Here we change the assumption on the third centre: we do not make any hypothesis on its position, and we obtain a perturbation of the 2-centre problem by assuming its intensity to be very small. Then, for a dense subset of possible positions of the perturbing centre in \mathbbR²{\mathbb{R}^2} , we prove the existence of uniformly hyperbolic invariant sets of periodic and chaotic almost collision orbits by the use of a general result of Bolotin and MacKay (Celest Mech Dyn Astron 77:49–75, 77:49–75, 2000; Celest Mech Dyn Astron 94(4):433–449, 2006). To apply it, we must preliminarily construct chains of collision arcs in a proper way. We succeed in doing that by the classical regularisation of the 2-centre problem and the use of the periodic orbits of the regularised problem passing through the third centre.     

Planar periodic orbits in exterior resonances with Neptune

In the framework of the planar restricted three-body problem we study a considerable number of resonances associated to the basic dynamical features of Kuiper belt and located between 30 and 48 a.u. Our study is based on the computation of resonant periodic orbits and their stability. Stable periodic orbits are surrounded by regular librations in phase space and in such domains the capture of trans-Neptunian object is possible. All the periodic orbits found are symmetric and there is an indication of the existence of asymmetric ones only in a few cases. In the present work first, second and third order resonances are under consideration. In the planar circular case we found that most of the periodic orbits are stable. The families of periodic orbits are temporarily interrupted by collisions but they continue up to relatively large values of the Jacobi constant and highly eccentric regular motion exists for all cases. In the elliptic problem and for a particular eccentricity value of the primary bodies, the periodic orbits are isolated. The corresponding families, where they belong to, bifurcate from specific periodic orbits of the circular problem and seem to continue up to the rectilinear problem. Both stable and unstable orbits are obtained for each case. In the elliptic problem, the unstable orbits found are associated with narrow chaotic domains in phase space. The evolution of the orbits, which are located in such chaotic domains, seems to be practically regular and bounded for long time intervals.     

2/1 resonant periodic orbits in three dimensional planetary systems

We consider the general spatial three body problem and study the dynamics of planetary systems consisting of a star and two planets which evolve into 2/1 mean motion resonance and into inclined orbits. Our study is focused on the periodic orbits of the system given in a suitable rotating frame. The stability of periodic orbits characterize the evolution of any planetary system with initial conditions in their vicinity. Stable periodic orbits are associated with long term regular evolution, while unstable periodic orbits are surrounded by regions of chaotic motion. We compute many families of symmetric periodic orbits by applying two schemes of analytical continuation. In the first scheme, we start from the 2/1 (or 1/2) resonant periodic orbits of the restricted problem and in the second scheme, we start from vertical critical periodic orbits of the general planar problem. Most of the periodic orbits are unstable, but many stable periodic orbits have been, also, found with mutual inclination up to 50^?–60^?, which may be related with the existence of real planetary systems.     

A numerical investigation of the one-dimensional newtonian three-body problem

Numerical orbit integrations have been conducted to characterize the types of trajectories in the one-dimensional Newtonian three-body problem with equal masses and negative energy. Essentially three different types of motions were found to exist. They may be classified according to the duration of the bound three-body state. There are zero-lifetime predictable trajectories, finite lifetime apparently chaotic orbits, and infinite lifetime quasi-periodic motions. The quasi-periodic orbits are confined to the neighbourhood of Schubart's stable periodic orbit. For all other trajectories the final state is of the type binary + single particle in both directions of time. The boundaries of the different orbit-type regions seem to be sharp. We present statistical results for the binding energies and for the duration of the bound three-body state. Properties of individual orbits are also summarized in the form of various graphical maps in a two-dimensional grid of parameters defining the orbit. Supported by the Academy of Finland.     

Periodic orbits and escapes in dynamical systems

We study the periodic orbits and the escapes in two different dynamical systems, namely (1) a classical system of two coupled oscillators, and (2) the Manko-Novikov metric which is a perturbation of the Kerr metric (a general relativistic system). We find their simple periodic orbits, their characteristics and their stability. Then we find their ordered and chaotic domains. As the energy goes beyond the escape energy, most chaotic orbits escape. In the first case we consider escapes to infinity, while in the second case we emphasize escapes to the central ??bumpy?? black hole. When the energy reaches its escape value, a particular family of periodic orbits reaches an infinite period and then the family disappears (the orbit escapes). As this family approaches termination it undergoes an infinity of equal period and double period bifurcations at transitions from stability to instability and vice versa. The bifurcating families continue to exist beyond the escape energy. We study the forms of the phase space for various energies, and the statistics of the chaotic and escaping orbits. The proportion of these orbits increases abruptly as the energy goes beyond the escape energy.     

Families of periodic orbits in the three-body problem

We show by a general argument that periodic solutions of the planar problem of three bodies (with given masses) form one-parameter families. This result is confirmed by numerical investigations: two orbits found earlier by Standish and Szebehely are shown to belong to continuous one-parameter families of periodic orbits. In general these orbits have a non-zero angular momentum, and the configuration after one period is rotated with respect to the initial configuration. Similar general arguments whow that in the three-dimensional problem, periodic orbits form also one-parameter families; in the one-dimensional problem, periodic orbits are isolated.     

Order and chaos in spiral galaxies

Spiral galaxies contain both ordered and chaotic orbits. In normal spirals the perturbations are weak (of order 2–10%) and most orbits are ordered. The density wave theory refers mainly to linear perturbations. Nonlinear effects appear in the outer parts of the open spirals (S_b, S_c) and produce the termination of these spirals near the 4/1 resonance. On the other hand in barred spirals the perturbations are relatively strong (of order 100%). Then the outer spirals and the envelope of the bar are composed mainly of chaotic orbits, while the main body of the bar is composed of ordered orbits. The chaotic orbits of the spiral arms of strong barred galaxies are sticky, i.e. they do not escape from the galaxy for at least a Hubble time. The forms of these spirals are delineated by the unstable manifolds of the unstable periodic orbits L_1, L_2 near the ends of the bar and of other unstable periodic orbits inside and outside corotation.     

on the relationship between fast lyapunov indicator and periodic orbits for symplectic mappings

The computation on a relatively short time of a quantity, related to the largest Lyapunov Characteristic Exponent, called Fast Lyapunov Indicator allows to discriminate between ordered and weak chaotic motion and also, under certain conditions, between resonant and non resonant regular orbits. The aim of this paper is to study numerically the relationship between the Fast Lyapunov Indicator values and the order of periodic orbits. Using the two-dimensional standard map as a model problem we have found that the Fast Lyapunov Indicator increases as the logarithm of the order of periodic orbits up to a given order. For higher order the Fast Lyapunov Indicator grows linearly with the order of the periodic orbits. We provide a simple model to explain the relationship that we have found between the values of the Fast Lyapunov Indicator, the order of the periodic orbits and also the minimum number of iterations needed to obtain the Fast Lyapunov Indicator values.     

Planar periodic orbits in three dipole problem

A systematic and detailed discussion of planar periodic orbits, of a charged particle moving under the influence of an electromagnetic field of three celestial bodies, is given for the first time. In this problem the periodic orbits are all asymmetric. Numerical procedures are applied to find the families of these orbits and to study their stability. Moreover, the bifurcations of these families with families of three dimensional asymmetric periodic orbits are given.     

Conditional Entropy

In this paper we show that the Conditional Entropy of nearby orbits may be a useful tool to explore the phase space associated to a given Hamiltonian. The arc length parameter along the orbits, instead of the time, is used as a random variable to compute the entropy. In the first part of this work we summarise the main analytical results to support this tool while, in the second part, we present numerical evidence that this technique is able to localise (stable) periodic and quasiperiodic orbits, 'aperiodic' orbits (chaotic motion) and unstable periodic orbits (the 'source' of chaotic motion). Besides, we show that this technique provides a measure of chaos which is similar to that given by the largest Lyapunov Characteristic Number. It is important to remark that this method is very simple to compute and does not require long time integrations, just realistic physical times. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the Relationship Between Fast Lyapunov Indicator and Periodic Orbits for Continuous Flows

It is already known (Froeschlé et al., 1997a) that the fast Lyapunov indicator (hereafter FLI), i.e. the computation on a relatively short time of a quantity related to the largest Lyapunov indicator, allows us to discriminate between ordered and weak chaotic motion. Using the FLI many results have been obtained on the standard map taken as a model problem. On this model we are not only able to discriminate between a short time weak chaotic motion and an ordered one, but also among regular motion between non resonant and resonant orbits. Moreover, periodic orbits are characterised by constant FLI values which appear to be related to the order of periodic orbits (Lega and Froeschlé, 2001). In the present paper we extend all these results to the case of continuous dynamical systems (the Hénon and Heiles system and the restricted three-body problem). Especially for the periodic orbits we need to introduce a new value: the orthogonal FLI in order to fully recover the results obtained for mappings.     

Periodic orbits of the elliptic restricted problem for the Sun-Jupiter-Saturn system

A systematic approach to generate periodic orbits in the elliptic restricted problem of three bodies in introduced. The approach is based on (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to periodic orbits in the elliptic restricted problem. Two families of periodic orbits of the elliptic restricted problem are found by this approach. The mass ratio of the primaries of these orbits is equal to that of the Sun-Jupiter system. The sidereal mean motions between the infinitesimal body and the smaller primary are in a 2:5 resonance, so as to approximate the Sun-Jupiter-Saturn system. The linear stability of these periodic orbits are studied as functions of the eccentricities of the primaries and of the infinitesimal body. The results show that both stable and unstable periodic orbits exist in the elliptic restricted problem that are close to the actual Sun-Jupiter-Saturn system. However, the periodic orbit closest to the actual Sun-Jupiter-Saturn system is (linearly) stable.     

Chaotic spiral galaxies

We study the role of asymptotic curves in supporting the spiral structure of a N-body model simulating a barred spiral galaxy. Chaotic orbits with initial conditions on the unstable asymptotic manifolds of the main unstable periodic orbits follow the shape of the periodic orbits for an initial interval of time and then they are diffused outwards along the spiral structure of the galaxy. Chaotic orbits having small deviations from the unstable periodic orbits, stay close and along the corresponding unstable asymptotic manifolds, supporting the spiral structure for more than 10 rotations of the bar. Chaotic orbits of different Jacobi constants support different parts of the spiral structure. We also study the diffusion rate of chaotic orbits outwards and find that the orbits that support the outer parts of the galaxy are diffused outwards more slowly than the orbits supporting the inner parts of the spiral structure.     

Boundary-value problem formulations for computing invariant manifolds and connecting orbits in the circular restricted three body problem

We demonstrate the remarkable effectiveness of boundary value formulations coupled to numerical continuation for the computation of stable and unstable manifolds in systems of ordinary differential equations. Specifically, we consider the circular restricted three-body problem (CR3BP), which models the motion of a satellite in an Earth–Moon-like system. The CR3BP has many well-known families of periodic orbits, such as the planar Lyapunov orbits and the non-planar vertical and halo orbits. We compute the unstable manifolds of selected vertical and halo orbits, which in several cases leads to the detection of heteroclinic connections from such a periodic orbit to invariant tori. Subsequent continuation of these connecting orbits with a suitable end point condition and allowing the energy level to vary leads to the further detection of apparent homoclinic connections from the base periodic orbit to itself, or the detection of heteroclinic connections from the base periodic orbit to other periodic orbits. Some of these connecting orbits are of potential interest in space mission design.     

Resonance transition periodic orbits in the circular restricted three-body problem

This work studies a special type of cislunar periodic orbits in the circular restricted three-body problem called resonance transition periodic orbits, which switch between different resonances and revolve about the secondary with multiple loops during one period. In the practical computation, families of multiple periodic orbits are identified first, and then the invariant manifolds emanating from the unstable multiple periodic orbits are taken to generate resonant homoclinic connections, which are used to determine the initial guesses for computing the desired periodic orbits by means of multiple-shooting scheme. The obtained periodic orbits have potential applications for the missions requiring long-term continuous observation of the secondary and tour missions in a multi-body environment.