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.     

Resonant Periodic Orbits of Trans-Neptunian Objects

We study two and three-dimensional resonant periodic orbits, usingthe model of the restricted three-body problem with the Sun andNeptune as primaries. The position and the stability character ofthe periodic orbits determine the structure of the phase space andthis will provide useful information on the stability and longterm evolution of trans-Neptunian objects. The circular planarmodel is used as the starting point. Families of periodic orbitsare computed at the exterior resonances 1/2, 2/3 and 3/4 withNeptune and these are used as a guide to select the energy levelsfor the computation of the Poincaré maps, so that all basicresonances are included in the study. Using the circular planarmodel as the basic model, we extend our study to more realisticmodels by considering an elliptic orbit of Neptune and introducingthe inclination of the orbit. Families of symmetric periodicorbits of the planar elliptic restricted three-body problem andthe three-dimensional problem are found. All these orbitsbifurcate from the families of periodic orbits of the planarcircular problem. The stability of all orbits is studied. Althoughthe resonant structure in the circular problem is similar for allresonances, the situation changes if the eccentricity of Neptuneor the inclination of the orbit is taken into account. All theseresults are combined to explain why in some resonances there aremany bodies and other resonances are empty.     

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.     

Some results on the global dynamics of the regularized restricted three-body problem with dissipation

We perform an analysis of the dynamics of the circular, restricted, planar three-body problem under the effect of different kinds of dissipation (linear, Stokes and Poynting–Robertson drags). Since the problem is singular, we implement a regularization technique in the style of Levi–Civita. The effect of the dissipation is often to decrease the semi-major axis; as a consequence the minor body collides with one of the primaries. In general, it is quite difficult to find non-collision orbits using random initial conditions. However, by means of the computation of the Fast Lyapunov Indicators (FLI), we obtain a global view of the dynamics. Precisely, we detect the regions of the phase space potentially belonging to basins of attraction. This investigation provides information on the different regions of the phase space, showing both collision and non-collision trajectories. Moreover, we find periodic orbit attractors for the case of linear and Stokes drags, while in the case of the Poynting–Robertson effect no other attractors are found beside the primaries, unless a fourth body is added to counterbalance the dissipative effect.     

Periodic orbits of the planar collision restricted 3-body problem

Using the continuation method we prove that the circular and the elliptic symmetric periodic orbits of the planar rotating Kepler problem can be continued into periodic orbits of the planar collision restricted 3-body problem. Additionally, we also continue to this restricted problem the so called “comet orbits”. An erratum to this article can be found at     

On the families of periodic orbits which bifurcate from the circular Sitnikov motions

In this paper we deal with the circular Sitnikov problem as a subsystem of the three-dimensional circular restricted three-body problem. It has a first analytical part where by using elliptic functions we give the analytical expressions for the solutions of the circular Sitnikov problem and for the period function of its family of periodic orbits. We also analyze the qualitative and quantitative behavior of the period function. In the second numerical part, we study the linear stability of the family of periodic orbits of the Sitnikov problem, and of the families of periodic orbits of the three-dimensional circular restricted three-body problem which bifurcate from them; and we follow these bifurcated families until they end in families of periodic orbits of the planar circular restricted three-body problem. We compare our results with the previous ones of other authors on this problem. Finally, the characteristic curves of some bifurcated families obtained for the mass parameter close to 1/2 are also described.     

The elliptic restricted problem at the 3 : 1 resonance

Four 3 : 1 resonant families of periodic orbits of the planar elliptic restricted three-body problem, in the Sun-Jupiter-asteroid system, have been computed. These families bifurcate from known families of the circular problem, which are also presented. Two of them, I _c , II _c bifurcate from the unstable region of the family of periodic orbits of the first kind (circular orbits of the asteroid) and are unstable and the other two, I _e , II _e , from the stable resonant 3 : 1 family of periodic orbits of the second kind (elliptic orbits of the asteroid). One of them is stable and the other is unstable. All the families of periodic orbits of the circular and the elliptic problem are compared with the corresponding fixed points of the averaged model used by several authors. The coincidence is good for the fixed points of the circular averaged model and the two families of the fixed points of the elliptic model corresponding to the families I _c , II _c , but is poor for the families I _e , II _e . A simple correction term to the averaged Hamiltonian of the elliptic model is proposed in this latter case, which makes the coincidence good. This, in fact, is equivalent to the construction of a new dynamical system, very close to the original one, which is simple and whose phase space has all the basic features of the elliptic restricted three-body problem.     

New results on the possible chaotic motion of Enceladus

New results are shown for the possible long term evolution of the orbit of Enceladus perturbed by Dione in the 2:1 resonance in the planar elliptic restricted problem. The numerical integrations of the secular variation equations (Ferraz-Mello and Dvorak, 1986) yield different types of chaotic behaviour. We also studied the onset of chaos for this dynamical system. Additionally, some illustrative pictures of the whole phase trajectory show the complex structure of phase space.     

Stable Chaos in High-Order Jovian Resonances

In a previous publication (Tsiganis et al. 2000, Icarus146, 240-252), we argued that the occurrence of stable chaos in the 12/7 mean motion resonance with Jupiter is related to the fact that there do not exist families of periodic orbits in the planar elliptic restricted problem and in the 3-D circular problem corresponding to this resonance. In the present paper we show that nonexistence of resonant periodic orbits, both for the planar and for the 3-D problem, also occurs in other jovian resonances—namely the 11/4, 22/9, 13/6, and 18/7—where cases of real asteroids on stable-chaotic orbits have been identified. This property may provide a “protection mechanism”, leading to semiconfinement of chaotic orbits and extremely slow migration in the space of proper elements, so that diffusion is practically unrelated to the value of the Lyapunov time, T_L, of chaotic orbits. However, we show that, in more complicated dynamical models, the long-term evolution of chaotic orbits initiated in the vicinity of these resonances may also be governed by secular resonances. Finally, we find that stable-chaotic orbits have a characteristic spectrum of autocorrelation times: for the action conjugate to the critical argument the autocorrelation time is of the order of the Lyapunov time, while for the eccentricity- and inclination-related actions the autocorrelation time may be longer than 10³T_L. This behavior is consistent with the trajectory being sticky around a manifold of lower-than-full dimensionality in phase space (e.g., a 4-D submanifold of the 5-D energy manifold in a three-degrees-of-freedom autonomus Hamiltonian system) and reflects the inability of these “flawed” resonances to modify secular motion significantly, at least for times of the order of 200 Myr.     

Origin and continuation of 3/2, 5/2, 3/1, 4/1 and 5/1 resonant periodic orbits in the circular and elliptic restricted three-body problem

We consider a planetary system consisting of two primaries, namely a star and a giant planet, and a massless secondary, say a terrestrial planet or an asteroid, which moves under their gravitational attraction. We study the dynamics of this system in the framework of the circular and elliptic restricted three-body problem, when the motion of the giant planet describes circular and elliptic orbits, respectively. Originating from the circular family, families of symmetric periodic orbits in the 3/2, 5/2, 3/1, 4/1 and 5/1 mean-motion resonances are continued in the circular and the elliptic problems. New bifurcation points from the circular to the elliptic problem are found for each of the above resonances, and thus, new families continued from these points are herein presented. Stable segments of periodic orbits were found at high eccentricity values of the already known families considered as whole unstable previously. Moreover, new isolated (not continued from bifurcation points) families are computed in the elliptic restricted problem. The majority of the new families mainly consists of stable periodic orbits at high eccentricities. The families of the 5/1 resonance are investigated for the first time in the restricted three-body problems. We highlight the effect of stable periodic orbits on the formation of stable regions in their vicinity and unveil the boundaries of such domains in phase space by computing maps of dynamical stability. The long-term stable evolution of the terrestrial planets or asteroids is dependent on the existence of regular domains in their dynamical neighbourhood in phase space, which could host them for long-time spans. This study, besides other celestial architectures that can be efficiently modelled by the circular and elliptic restricted problems, is particularly appropriate for the discovery of terrestrial companions among the single-giant planet systems discovered so far.     

The Phase Space Structure Around L4 in the Restricted Three-Body Problem

The phase space structure around L ₄ in the restricted three-body problem is investigated. The connection between the long period family emanating from L ₄ and the very complex structure of the stability region is shown by using the method of Poincarés surface of section. The zero initial velocity stability region around L ₄ is determined by using a method based on the calculation of finite-time Lyapunov characteristic numbers. It is shown that the boundary of the stability region in the configuration space is formed by orbits suffering slow chaotic diffusion.     

Symmetric and Nonsymmetric Periodic Orbits in the Exterior Mean Motion Resonances with Neptune

We present families of periodic orbits and their stability for the exterior mean motion resonances 1:2, 1:3 and 1:4 with Neptune in the framework of the planar circular restricted three-body problem. We found that in each resonance there exist two branches of symmetric elliptic periodic orbits with stable and unstable segments. Asymmetric periodic orbits bifurcate from the corresponding symmetric ones. Asymmetric periodic orbits are stable and the motion in their neighbourhood is a libration with respect to the resonant angle variable. In all the families of asymmetric periodic orbits the eccentricity extends to high values. Poincaré sections reveal the changes of the topology in phase space.     

Motion of dust in mean motion resonances with planets

Effect of stellar electromagnetic radiation on the motion of spherical dust particle in mean motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting–Robertson (P–R) effect yields monotonic secular evolution of eccentricity when the particle is trapped in the resonance. Planar elliptic restricted three-body problem with the P–R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean motion orbital resonances. Evolutions of argument of perihelion in the planar circular and elliptical restricted three-body problems are shown. Numerical integrations show that an analytic expression for the secular time derivative of the particle’s argument of perihelion does not exist, if only dependence on semimajor axis, eccentricity and argument of perihelion is admitted. Connection between the shift of perihelion and oscillations in secular eccentricity is presented for the planar elliptic restricted three-body problem with the P–R effect. Period of the oscillations corresponds to the period of one revolution of perihelion. Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: (i) does not have any significant influence on the secular evolution of eccentricity, (ii) causes that the shift of perihelion is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.     

Bounded orbits in the elliptic restricted three-body problem

Special solutions of the planar rectilinear elliptic restricted 3-body problem are investigated for the limiting case e=1. Numerical integration is performed for primaries of equal masses. Starting values which define circular orbit solutions lead to bounded solutions if the initial radius a₀ is larger than 3.74 in units of the primaries' semimajor axis a. A comparison with the Eulerian two-fixedcentre problem is presented in order to understand qualitatively the characteristic features of bounded orbits and the transition to escape orbits.     

Applications of lagrangian coherent structures to expression of invariant manifolds in astrodynamics

This paper investigates the relationship between invariant manifold and Lagrangian coherent structure (LCS) in dynamical systems. LCS is defined as the ridge of finite-time Lyapunov exponent (FTLE) field, and is proving to be excellent platform for studies of stable and unstable manifold in flows with arbitrary time dependence. In this study, the LCS tool is applied to autonomous systems, simple pendulum and planar circular restricted three-body problem (PCR3BP), and also non-autonomous ones, double-gyre flow and bicircular problem (BCP). A comparison between LCS and invariant manifold is presented.     

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.     

From the circular to the spatial elliptic restricted three-body problem

We deal with the study of the spatial restricted three-body problem in the case where the small particle is far from the primaries, that is, the so-called comet case. We consider the circular problem, apply double averaging and compute the relative equilibria of the reduced system. It appears that, in the circular problem, we find not only part of the equilibria existing in the elliptic case, but also new ones. These critical points are in correspondence with periodic and quasiperiodic orbits and invariant tori of the non-averaged Hamiltonian. We explain carefully the transition between the circular and the elliptic problems. Moreover, from the relative equilibria of elliptic type, we obtain invariant 3-tori of the original system.     

Consecutive collision orbits in the limiting case μ=0 of the elliptic restricted problem

A timing condition for consecutive collision orbits in the planar, circular three-body problem has been extended to the elliptic restricted problem for =0. The expression developed relates eccentric anomalies at the time of collision. Some families of solutions are presented.     

Family <Emphasis Type="Italic">h</Emphasis> of periodic solutions of the restricted problem for big μ

The results of the computation of the family h of symmetric periodic solutions of the circular planar restricted three-body problem for μ = 0.3, 0.4, and 0.5 are presented. This family begins with retrograde circular orbits around a massive body. Associated with each value of μ are the table of critical orbits, the orbit pictures, the graphs of characteristics of the family in four coordinate systems, and the graphs of the period and traces (planar and vertical). Regularities on the family and its evolution as μ increased were observed.     

Numerical results to the Sitnikov-problem

We present numerical results of the so-called Sitnikov-problem, a special case of the three-dimensional elliptic restricted three-body problem. Here the two primaries have equal masses and the third body moves perpendicular to the plane of the primaries' orbit through their barycenter. The circular problem is integrable through elliptic integrals; the elliptic case offers a surprisingly great variety of motions which are until now not very well known. Very interesting work was done by J. Moser in connection with the original Sitnikov-paper itself, but the results are only valid for special types of orbits. As the perturbation approach needs to have small parameters in the system we took in our experiments as initial conditions for the work moderate eccentricities for the primaries' orbit (0.33e _primaries 0.66) and also a range of initial conditions for the distance of the 3 ^rd body (= the planet) from very close to the primaries orbital plane of motion up to distance 2 times the semi-major axes of their orbit. To visualize the complexity of motions we present some special orbits and show also the development of Poincaré surfaces of section with the eccentricity as a parameter. Finally a table shows the structure of phase space for these moderately chosen eccentricities.