On the Structure of Symplectic Mappings. The Fast Lyapunov Indicator: a Very Sensitive Tool

It is already known (Froeschlé, Lega and Gonczi, 1997) that the Fast Lyapunov Indicator (FLI), that is the computation on a relatively short time of the largest Lyapunov indicator, allows to discriminate between ordered and weak chaotic motion. We have found that, under certain conditions, the FLI also discriminates between resonant and non-resonant orbits, not only for two-dimensional symplectic mappings but also for higher dimensional ones. Using this indicator, we present an example of the Arnold web detection for four and six-dimensional symplectic maps. We show that this method allows to detect the global transition of the system from an exponentially stable Nekhoroshevs like regime to the diffusive Chirikovs one.     

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.     

Discrete Virial Theorem

We reexamine the classical virial theorem for bounded orbits of arbitrary autonomous Hamiltonian systems possessing both regular and chaotic orbits. New and useful forms of the virial theorem are obtained for natural Hamiltonian flows of arbitrary dimension. A discrete virial theorem is derived for invariant circles and periodic orbits of natural symplectic maps. A weak and a strong form of the virial theorem are proven for both flows and maps. While the Birkhoff Ergodic Theorem guarantees the existence of the relevant time averages for both regular and chaotic orbits, the convergence is very rapid for the former and extremely slow for the latter. This circumstance leads to a simple and efficient measure of chaoticity. The results are applied to several problems of current physical interest, including the Hénon–Heiles system, weak chaos in the standard map, and a 4D Froeschlé map.     

The Lyapunov characteristic exponents as indicators of stochasticity in the restricted three-body problem

The Lyapunov characteristic numbers are used to continue the study of a larger sample of orbits in the three-body three-dimensional restricted problem than that previously studied by Froeschle [6]. These numbers are found to be sensitive indicators of stochasticity and the results confirm those found by Froeschlé.Equipe de Recherche Associee au CNRS.     

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.     

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.     

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.     

Symmetric and asymmetric 3:1 resonant periodic orbits with an application to the 55Cnc extra-solar system

We study the dynamics of 3:1 resonant motion for planetary systems with two planets, based on the model of the general planar three body problem. The exact mean motion resonance corresponds to periodic motion (in a rotating frame) and the basic families of symmetric and asymmetric periodic orbits are computed. Four symmetric families bifurcate from the family of circular orbits of the two planets. Asymmetric families bifurcate from the symmetric families, at the critical points, where the stability character changes. There exist also asymmetric families that are independent of the above mentioned families. Bounded librations exist close to the stable periodic orbits. Therefore, such periodic orbits (symmetric or asymmetric) determine the possible stable configurations of a 3:1 resonant planetary system, even if the orbits of the two planets intersect. For the masses of the system 55Cnc most of the periodic orbits are unstable and they are associated with chaotic motion. There exist however stable symmetric and asymmetric orbits, corresponding to regular trajectories along which the critical angles librate. The 55Cnc extra-solar system is located in a stable domain of the phase space, centered at an asymmetric periodic orbit.     

Stable Chaos versus Kirkwood Gaps in the Asteroid Belt: A Comparative Study of Mean Motion Resonances

We have shown, in previous publications, that stable chaos is associated with medium/high-order mean motion resonances with Jupiter, for which there exist no resonant periodic orbits in the framework of the elliptic restricted three-body problem. This topological “defect” results in the absence of the most efficient mechanism of eccentricity transport (i.e., large-amplitude modulation on a short time scale) in three-body models. Thus, chaotic diffusion of the orbital elements can be quite slow, while there can also exist a nonnegligible set of chaotic orbits which are semiconfined (stable chaos) by “quasi-barriers” in the phase space. In the present paper we extend our study to all mean motion resonances of order q≤9 in the inner main belt (1.9-3.3 AU) and q≤7 in the outer belt (3.3-3.9 AU). We find that, out of the 34 resonances studied, only 8 possess resonant periodic orbits that are continued from the circular to the elliptic three-body problem (regular families), namely, the 2/1, 3/1, 4/1, and 5/2 in the inner belt and the 7/4, 5/3, 11/7, and 3/2 in the outer belt. Numerical results indicate that the 7/3 resonance also carries periodic orbits but, unlike the aforementioned resonances, 7/3-periodic orbits belong to an irregular family. Note that the five inner-belt resonances that carry periodic orbits correspond to the location of the main Kirkwood gaps, while the three outer-belt resonances correspond to gaps in the distribution of outer-belt asteroids noted by Holman and Murray (1996, Astron. J.112, 1278-1293), except for the 3/2 case where the Hildas reside. Fast, intermittent eccentricity increase is found in resonances possessing periodic orbits. In the remaining resonances the time-averaged elements of chaotic orbits are, in general, quite stable, at least for times t∼250 Myr. This slow diffusion picture does not change qualitatively, even if more perturbing planets are included in the model.     

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.     

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.     

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 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.     

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.     

Using FLI maps for preliminary spacecraft trajectory design in multi-body environments

Fast Lyapunov Indicator (FLI) maps are presented as a tool for solving spacecraft preliminary trajectory design problems in multi-body environments with long-term stability requirements. In particular, the FLI maps are shown to provide a global overview of the dynamics in the restricted three-body problem that can guide mission designers in selecting long-term stable regions of phase space which are inherently more robust to model parameter perturbations. The FLI is also shown to numerically detect the normally hyperbolic manifolds associated with unstable periodic orbits. These, in turn, provide a global map of the principal heteroclinic connections between the various resonance regions which form the basic backbone of dynamical transfers design. Examples of maps and transfers are provided in the restricted three-body problem modeling the Jupiter–Europa system.     

Estimate of the Transition Value of Librational Invariant Curves

We investigate the break-down threshold of librational invariant curves. As a model problem, we consider a variant of a mapping introduced by M. Hénon, which well describes the dynamics of librational motions surrounding a stable invariant point. We verify in concrete examples the applicability of Greene's method, by computing the instability transition values of a sequence of periodic orbits approaching an invariant curve with fixed noble frequency. However, this method requires the knowledge of the location of the periodic orbits within a very good approximation. This task appears to be difficult to realize for a libration regime, due to the different topology of the phase space. To compute the break-down threshold, we tried an alternative method very easy to implement, based on the computation of the fast Lyapunov indicators and frequency analysis. Such technique does not require the knowledge of the periodic orbits, but again, it appears very difficult to have a precision better than Greene's method for the computation of the critical parameter.     

Ordered and chaotic Bohmian trajectories

We discuss the issue of ordered and chaotic trajectories in the Bohmian approach of Quantum Mechanics from points of view relevant to the methods of Celestial Mechanics. The Bohmian approach gives the same results as the orthodox (Copenhagen) approach, but it considers also underlying trajectories guided by the wave. The Bohmian trajectories are rather different from the corresponding classical trajectories. We give examples of a classical chaotic system that is ordered quantum-mechanically and of a classically ordered system that is mostly chaotic quantum mechanically. Then we consider quantum periodic orbits and ordered orbits, that can be represented by formal series of the “third integral” type, and we study their asymptotic properties leading to estimates of exponential stability. Such orbits do not approach the “nodal points” where the wavefunction ψ vanishes. On the other hand, when an orbit comes close to a nodal point, chaos is generated in the neighborhood of a hyperbolic point (called X-point). The generation of chaos is maximum when the X-point is close to the nodal point. Finally we remark that high order periodic orbits may behave as “effectively ordered” or “effectively chaotic” for long times before reaching the period.     

Complex statistics in Hamiltonian barred galaxy models

We use probability density functions (pdfs) of sums of orbit coordinates, over time intervals of the order of one Hubble time, to distinguish weakly from strongly chaotic orbits in a barred galaxy model. We find that, in the weakly chaotic case, quasi-stationary states arise, whose pdfs are well approximated by q-Gaussian functions (with 1 <?q < 3), while strong chaos is identified by pdfs which quickly tend to Gaussians (q =?1). Typical examples of weakly chaotic orbits are those that ??stick?? to islands of ordered motion. Their presence in rotating galaxy models has been investigated thoroughly in recent years due to their ability to support galaxy structures for relatively long time scales. In this paper, we demonstrate, on specific orbits of 2 and 3 degree of freedom barred galaxy models, that the proposed statistical approach can distinguish weakly from strongly chaotic motion accurately and efficiently, especially in cases where Lyapunov exponents and other local dynamic indicators appear to be inconclusive.     

On the Calculation of Lyapunov Indicators with Post-stabilization in a Weyl Field

We present details of a work aiming at the overestimation of Lyapunov exponents defined by the geodesic deviation equations in the previous work. The geodesic deviation vector with post-stabilization is used to compute the fast Lyapunov indicator, considered to be a very sensitive tool for discrimination between ordered or weakly chaotic motions. We make a detailed study of the dynamics in the superposed Weyl field between a black hole and shell of octopoles by using the fast Lyapunov indicator with the Poincare surface of section. In particular, we examine the effects on the dynamics around the fixed points, of varying one of the three parameters (specific energy E, specific angular momentum L and octopolar moment O), while keeping the other two fixed, and identify the intervals of the varying parameter where the motion is regular or chaotic.     

The 1/1 resonance in extrasolar systems

We present families of symmetric and asymmetric periodic orbits at the 1/1 resonance, for a planetary system consisting of a star and two small bodies, in comparison to the star, moving in the same plane under their mutual gravitational attraction. The stable 1/1 resonant periodic orbits belong to a family which has a planetary branch, with the two planets moving in nearly Keplerian orbits with non zero eccentricities and a satellite branch, where the gravitational interaction between the two planets dominates the attraction from the star and the two planets form a close binary which revolves around the star. The stability regions around periodic orbits along the family are studied. Next, we study the dynamical evolution in time of a planetary system with two planets which is initially trapped in a stable 1/1 resonant periodic motion, when a drag force is included in the system. We prove that if we start with a 1/1 resonant planetary system with large eccentricities, the system migrates, due to the drag force, along the family of periodic orbits and is finally trapped in a satellite orbit. This, in principle, provides a mechanism for the generation of a satellite system: we start with a planetary system and the final stage is a system where the two small bodies form a close binary whose center of mass revolves around the star.