On quasi-satellite periodic motion in asteroid and planetary dynamics

Applying the method of analytical continuation of periodic orbits, we study quasi-satellite motion in the framework of the three-body problem. In the simplest, yet not trivial model, namely the planar circular restricted problem, it is known that quasi-satellite motion is associated with a family of periodic solutions, called family f, which consists of 1:1 resonant retrograde orbits. In our study, we determine the critical orbits of family f that are continued both in the elliptic and in the spatial models and compute the corresponding families that are generated and consist the backbone of the quasi-satellite regime in the restricted model. Then, we show the continuation of these families in the general three-body problem, we verify and explain previous computations and show the existence of a new family of spatial orbits. The linear stability of periodic orbits is also studied. Stable periodic orbits unravel regimes of regular motion in phase space where 1:1 resonant angles librate. Such regimes, which exist even for high eccentricities and inclinations, may consist dynamical regions where long-lived asteroids or co-orbital exoplanets can be found.     

On the 2/1 resonant planetary dynamics – periodic orbits and dynamical stability

The 2/1 resonant dynamics of a two-planet planar system is studied within the framework of the three-body problem by computing families of periodic orbits and their linear stability. The continuation of resonant periodic orbits from the restricted to the general problem is studied in a systematic way. Starting from the Keplerian unperturbed system, we obtain the resonant families of the circular restricted problem. Then, we find all the families of the resonant elliptic restricted three-body problem, which bifurcate from the circular model. All these families are continued to the general three-body problem, and in this way we can obtain a global picture of all the families of periodic orbits of a two-planet resonant system. The parametric continuation, within the framework of the general problem, takes place by varying the planetary mass ratio ρ. We obtain bifurcations which are caused either due to collisions of the families in the space of initial conditions or due to the vanishing of bifurcation points. Our study refers to the whole range of planetary mass ratio values  [ρ∈ (0, ∞)]  and, therefore we include the passage from external to internal resonances. Thus, we can obtain all possible stable configurations in a systematic way. As an application, we consider the dynamics of four known planetary systems at the 2/1 resonance and we examine if they are associated with a stable periodic orbit.     

The rectilinear three-body problem as a basis for studying highly eccentric systems

The rectilinear elliptic restricted three-body problem (TBP) is the limiting case of the elliptic restricted TBP when the motion of the primaries is described by a Keplerian ellipse with eccentricity \(e'=1\), but the collision of the primaries is assumed to be a non-singular point. The rectilinear model has been proposed as a starting model for studying the dynamics of motion around highly eccentric binary systems. Broucke (AIAA J 7:1003–1009, 1969) explored the rectilinear problem and obtained isolated periodic orbits for mass parameter \(\mu =0.5\) (equal masses of the primaries). We found that all orbits obtained by Broucke are linearly unstable. We extend Broucke’s computations by using a finer search for symmetric periodic orbits and computing their linear stability. We found a large number of periodic orbits, but only eight of them were found to be linearly stable and are associated with particular mean motion resonances. These stable orbits are used as generating orbits for continuation with respect to \(\mu \) and \(e'<1\). Also, continuation of periodic solutions with respect to the mass of the small body can be applied by using the general TBP. FLI maps of dynamical stability show that stable periodic orbits are surrounded in phase space with regions of regular orbits indicating that systems of very highly eccentric orbits can be found in stable resonant configurations. As an application we present a stability study for the planetary system HD7449.     

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.     

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.     

Numerical continuation of periodic orbits from the restricted to the general 3-dimensional 3-body problem

In this paper several monoparametric families of periodic orbits of the 3-dimensional general 3-body problem are presented. These families are found by numerical continuation with respect to the small massm ₃, of some periodic orbits which belong to a family of 3-dimensional periodic orbits of the restricted elliptic problem.     

Classical periodic orbits and quantum mechanical eigenvalues and eigenfunctions

We have calculated several families of classical periodic orbits in simple Hamiltonian systems of two degrees of freedom and the corresponding quantum mechanical eigenvalues and eigenfuctions. We have found that in most cases the eigenfunctions have their maxima and minima on some simple periodic orbits. These periodic orbits are of several resonant types and can be either stable or unstable. In the latter case the quantum Poincaré surfaces of section are very different from the classical Poincaré surfaces of section.     

Periodic Orbits in the Restricted Three Body Problem with Radiation Pressure

This paper deals with the Restricted Three Body Problem (RTBP) in which we assume that the primaries are radiation sources and the influence of the radiation pressure on the gravitational forces is considered; in particular, we are interested in finding families of periodic orbits under theses forces. By means of some modifications to the method of numerical continuation of natural families of periodic orbits, we find several families of periodic orbits, both in two and three dimensions. As starters for our method we use some known periodic orbits in the classical RTBP. This revised version was published online in August 2006 with corrections to the Cover Date.     

Resonant motion in the restricted three body problem

The resonant structure of the restricted three body problem for the Sun- Jupiter asteroid system in the plane is studied, both for a circular and an elliptic orbit of Jupiter. Three typical resonances are studied, the 2 : 1, 3 : 1 and 4 : 1 mean motion resonance of the asteroid with Jupiter. The structure of the phase space is topologically different in these cases. These are typical for all other resonances in the asteroid problem. In each case we start with the unperturbed two-body system Sun-asteroid and we study the continuation of the periodic orbits when the perturbation due to a circular orbit of Jupiter is introduced. Families of periodic orbits of the first and of the second kind are presented. The structure of the phase space on a surface of section is also given. Next, we study the families of periodic orbits of the asteroid in the elliptic restricted problem with the eccentricity of Jupiter as a parameter. These orbits bifurcate from the families of the circular problem. Finally, we compare the above families of periodic orbits with the corresponding families of fixed points of the averaged problem. Different averaged Hamiltonians are considered in each resonance and the range of validity of each model is discussed.     

Error analyses of resonant orbits for geodesy

An error analysis of resonant orbits for geodesy indicates that attempts to use resonance to recover high order geopotential coefficients may be seriously hampered by errors in the geopotential. This effect, plus the very high correlations (up to .99) of the resonant coefficients with each other and the orbital period in single satellite solutions, makesindividual resonant orbits of limited value for geodesy. Multiple-satellite, single-plane solutions are only a slight improvement over the single satellite case. Accurate determination of high order coefficients from low altitude resonant satellites requires multiple orbit planes and small drift-periods to reduce correlations and effects of errors of non-resonant geopotential terms. Also, the effects of gravity model errors on low-altitude resonant satellites make the use of tracking arcs exceeding two to three weeks of doubtful validity. Because high-altitude resonant orbits are less affected by non-resonant terms in the geopotential, much longer tracking arcs can be used for them.     

A Map for a Group of Resonant Cases in a quartic Galactic Hamiltonian

We present a map for the study of resonant motion in a potential made up of two harmonic oscillators with quartic perturbing terms. This potential can be considered to describe motion in the central parts of non-rotating elliptical galaxies. The map is based on the averaged Hamiltonian. Adding on a semi-empirical basis suitable terms in the unperturbed averaged Hamiltonian, corresponding to the 1:1 resonant case, we are able to construct a map describing motion in several resonant cases. The map is used in order to find thex − p _x Poincare phase plane for each resonance. Comparing the results of the map, with those obtained by numerical integration of the equation of motion, we observe, that the map describes satisfactorily the broad features of orbits in all studied cases for regular motion. There are cases where the map describes satisfactorily the properties of the chaotic orbits as well.     

Doubly-symmetric motions in the elliptic problem

Poincaré's continuation method is applied to the elliptic restricted problem for the computation of four families of doubly symmetric three-dimensional periodic orbits emanating from similar orbits corresponding to zero eccentricity. The results are given in four tables and the orbits are characterized in regard to their stability.     

Periodic orbit families near the 4:1 jovian resonance

An enlarged averaged Hamiltonian is introduced to compute several families of periodic orbits of the planar elliptic 3-body problem, in the Sun–Jupiter–Asteroid system, near the 4:1 resonance. Four resonant critical point families are found and their stability is studied. The families of symmetric periodic orbits of the elliptic problem appear near the corresponding fixed points computed in this model. There is a good agreement for moderate eccentricity of the asteroid for three of these families, whereas the remaining family cannot be considered as a family of periodic orbits of the real model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Resonant orbits of grains and the formation of satellites

Grains, an abundant constituent of the former solar system, will have had a high probability of being driven into orbits resonant with major bodies already formed. This arises because of the presence of gas drag and Poynting-Robertson drag on small grains, providing the dissipation necessary to concentrate matter into special orbits. Since the mean density in resonant orbits can be built up by such a process without limit, these may become the favored orbits for gravitational contraction to gather material into major bodies. Satellite formation processes may therefore depend upon the buildup of resonant lanes of dust grains around the parent body. Saturn's rings are possibly one example of such lanes, though an unsuitable one for the final step of satellite formation on account of their being too close to Saturn.     

Orbits and manifolds near the equilibrium points around a rotating asteroid

We study the orbits and manifolds near the equilibrium points of a rotating asteroid. The linearised equations of motion relative to the equilibrium points in the gravitational field of a rotating asteroid, the characteristic equation and the stable conditions of the equilibrium points are derived and discussed. First, a new metric is presented to link the orbit and the geodesic of the smooth manifold. Then, using the eigenvalues of the characteristic equation, the equilibrium points are classified into 8 cases. A theorem is presented and proved to describe the structure of the submanifold as well as the stable and unstable behaviours of a massless test particle near the equilibrium points. The linearly stable, the non-resonant unstable, and the resonant equilibrium points are discussed. There are three families of periodic orbits and four families of quasi-periodic orbits near the linearly stable equilibrium point. For the non-resonant unstable equilibrium points, there are four relevant cases; for the periodic orbit and the quasi-periodic orbit, the structures of the submanifold and the subspace near the equilibrium points are studied for each case. For the resonant equilibrium points, the dimension of the resonant manifold is greater than 4, and we find at least one family of periodic orbits near the resonant equilibrium points. As an application of the theory developed here, we study relevant orbits for the asteroids 216 Kleopatra, 1620 Geographos, 4769 Castalia and 6489 Golevka.     

Families of first kind periodic solutions of the spatial and planar elliptic restricted three-body problem with collision

We consider a restricted three-body problem where the primaries are moving in an elliptic collision orbit and the infinitesimal mass moves in a three dimensional space. This paper is devoted to prove analytically the existence of several families of symmetric periodic solutions as continuation of Keplerian circular orbits. In our approach the perturbing parameter is related with the energy of the primaries.     

Periodic Orbits of a Collinear Restricted Three-Body Problem

In this paper we study symmetric periodic orbits of a collinear restricted three-body problem, when the middle mass is the largest one. These symmetric periodic orbits are obtained from analytic continuation of symmetric periodic orbits of two collinear two-body problems.     

Symmetric and asymmetric librations in extrasolar planetary systems: a global view

We present a global view of the resonant structure of the phase space of a planetary system with two planets, moving in the same plane, as obtained from the set of the families of periodic orbits. An important tool to understand the topology of the phase space is to determine the position and the stability character of the families of periodic orbits. The region of the phase space close to a stable periodic orbit corresponds to stable, quasi periodic librations. In these regions it is possible for an extrasolar planetary system to exist, or to be trapped following a migration process due to dissipative forces. The mean motion resonances are associated with periodic orbits in a rotating frame, which means that the relative configuration is repeated in space. We start the study with the family of symmetric periodic orbits with nearly circular orbits of the two planets. Along this family the ratio of the periods of the two planets varies, and passes through rational values, which correspond to resonances. At these resonant points we have bifurcations of families of resonant elliptic periodic orbits. There are three topologically different resonances: (1) the resonances (n + 1):n, (2:1, 3:2, ...), (2) the resonances (2n + 1):(2n-1), (3:1, 5:3, ...) and (3) all other resonances. The topology at each one of the above three types of resonances is studied, for different values of the sum and of the ratio of the planetary masses. Both symmetric and asymmetric resonant elliptic periodic orbits exist. In general, the symmetric elliptic families bifurcate from the circular family, and the asymmetric elliptic families bifurcate from the symmetric elliptic families. The results are compared with the position of some observed extrasolar planetary systems. In some cases (e.g., Gliese 876) the observed system lies, with a very good accuracy, on the stable part of a family of resonant periodic orbits.     

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange’s planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting–Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.     

The three-dimensional general three-body problem: Determination of periodic orbits

The three-dimensional general three-body problem is formulated suitably for the numerical determination of periodic orbits either directly or by continuation from the three-dimensional periodic orbits of the restricted problem. The symmetry properties of the equations of motion are established and the algorithms for the numerical determination of families of periodic orbits are outlined. A normalization scheme based on the concept of the invariable plane is introduced to simplify the process. All three types of symmetric orbit, as well as the general type of asymmetric orrbit, are considered. Many threedimmensional p periodic orbits are given.