The origin of chaotic behaviour in the Miranda-Umbriel 3 : 1 resonances

We have investigated the pericentric resonances through which Miranda and Umbriel are believed to have passed when, due to tidal evolution, their orbital mean motions reached a 3 : 1 commensurability. Our investigation is based upon a perturbative treatment. The predictions of this theory are in good agreement with the results of numerical integrations concerning both the extend of the chaotic layers generated by the separatrices of the primary resonances and the location of the secondary resonances. The effect of tidal evolution is discussed on the bases of the adiatatic invariant theory and its extension to separatrix crossing. We recover qualitatively the mean features of the numerical experiments of Tittermore and Wisdom (1988–1989), Dermott et al (1988) and Malhotra and Dermott (1989).     

Surfaces of section in the Miranda-Umbriel 3:1 inclination problem

The recent numerical simulations of Tittemore and Wisdom (1988, 1989, 1990) and Dermottet al. (1988), Malhotra and Dermott (1990) concerning the tidal evolution through resonances of some pairs of Uranian satellites have revealed interesting dynamical phenomena related to the interactions between close-by resonances. These interactions produce chaotic layers and strong secondary resonances. The slow evolution of the satellite orbits in this dynamical landscape is responsible for temporary capture into resonance, enhancement of eccentricity or inclination and subsequent escape from resonance. The present contribution aims at developing analytical tools for predicting the location and size of chaotic layers and secondary resonances. The problem of the 3:1 inclination resonance between Miranda and Umbriel is analysed.     

Lunar shape does not record a past eccentric orbit

The Moon has long been known to have an overall shape not consistent with expected past tidal forces. It has recently been suggested (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) that the present lunar moments of inertia indicate a past high-eccentricity orbit and, possibly, a past non-synchronous spin-orbit resonance. Here I show that the match between the lunar shape and the proposed orbital and spin states is much less conclusive than initially proposed. Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) spin and shape evolution scenarios also completely ignore the physics of the capture into such resonances, which require prior permanent deformation, as well as tidal despinning to the relevant resonance. If the early lunar orbit was eccentric, the Moon would have been rotating at an equilibrium non-synchronous rate determined by it eccentricity. This equilibrium supersynchronous rotation would be much too fast to allow a synchronous spin-orbit lock at e = 0.49, while the capture into the 3:2 resonance is possible only for a very constrained lunar eccentricity history and assuming some early permanent lunar tri-axiality. Here I show that large impacts in the early history of the Moon would have frequently disrupted this putative resonant rotation, making the rotation and eccentricity solutions of Garrick-Bethell et al. (Garrick-Bethell, I., Wisdom, J., Zuber, M.T. [2006]. Science 313, 652-655) unstable. I conclude that the present lunar shape cannot be used to support the hypothesis of an early eccentric lunar orbit.     

Adiabatic chaos in the Prometheus–Pandora system

The chaotic orbital motion of Prometheus and Pandora, the 16th and 17th satellites of Saturn, is studied. Chaos in their orbital motion, as found by Goldreich & Rappaport and Renner & Sicardy, is due to interaction of resonances in the resonance multiplet corresponding to the 121:118 commensurability of the mean motions of the satellites. It is shown rigorously that the system moves in adiabatic regime. The Lyapunov time (the 'time horizon of predictability' of the motion) is calculated analytically and compared to the available numerical–experimental estimates. For this purpose, a method of analytical estimation of the maximum Lyapunov exponent in the perturbed pendulum model of non-linear resonance is applied. The method is based on the separatrix map theory. An analytical estimate of the width of the chaotic layer is made as well, based on the same theory. The ranges of chaotic diffusion in the mean motion are shown to be almost twice as big compared to previous estimates for both satellites.     

Local and global diffusion in the Arnold web of a priori unstable systems

Using the numerical techniques developed by Froeschlé et al. (Science 289 (5487): 2108–2110, 2000) and by Lega et al. (Physica D 182: 179–187, 2003) we have studied diffusion and stochastic properties of an a priori unstable 4D symplectic map. We have found two different kinds of diffusion that coexist for values of the perturbation below the critical value for the Chirikov overlapping of resonances. A fast diffusion along some resonant lines that exist already in the unperturbed case and a slow diffusion occurring in regions of the phase space far from such resonances. The latter one, although the system does not satisfy the Nekhoroshev hypothesis, decreases faster than a power law and possibly exponentially. We compare the diffusion coefficient to the indicators of stochasticity like the Lyapunov exponent and filling factor showing their behavior for chaotic orbits in regions of the Arnold web where the secondary resonances appear, or where they overlap.     

Hamiltonian Stability of Spin–Orbit Resonances in Celestial Mechanics

The behaviour of ‘resonances’ in the spin-orbit coupling in celestial mechanics is investigated in a conservative setting. We consider a Hamiltonian nearly-integrable model describing an approximation of the spin-orbit interaction. The continuous system is reduced to a mapping by integrating the equations of motion through a symplectic algorithm. We study numerically the stability of periodic orbits associated to the above mapping by looking at the eigenvalues of the matrix of the linearized map over the full cycle of the periodic orbit. In particular, the value of the trace of the matrix is related to the stability character of the periodic orbit. We denote by ε_* (p/q) the value of the perturbing parameter at which a given elliptic periodic orbit with frequency p/q becomes unstable. A plot of the critical function ε_* (p/q) versus the frequency at different orbital eccentricities shows significant peaks at the synchronous resonance (for low eccentricities) and at the synchronous and 3:2 resonances (at higher eccentricities) in good agreement with astronomical observations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Chaotic Dynamics of Satellite Systems

We consider the problem of calculating the Lyapunov time (the characteristic time of predictable dynamics) of chaotic motion in the vicinity of separatrices of orbital resonances in satellite systems. The primary objects of study are the chaotic regimes that have occurred in the history of the orbital dynamics of the second and fifth Uranian satellites (Umbriel and Miranda) and the first and third Saturnian satellites (Mimas and Tethys). We study the dynamics in the vicinity of separatrices of the resonance multiplets corresponding to the 3 : 1 commensurability of mean motions of Miranda and Umbriel and the multiplets corresponding to the 2 : 1 commensurability of mean motions of Mimas and Tethys. These chaotic regimes have most probably contributed much to the long-term orbital evolution of the two satellite systems. The equations of motion have been numerically integrated to estimate the Lyapunov time in models corresponding to various epochs of the system evolution. Analytical estimates of the Lyapunov time have been obtained by a method (Shevchenko, 2002) based on the separatrix map theory. The analytical estimates have been compared to estimates obtained by direct numerical integration.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 364–374.Original Russian Text Copyright © 2005 by Mel’nikov, Shevchenko.     

On the chaotic orbital dynamics of the planet in the system 16 Cyg

The chaotic orbital dynamics of the planet in the wide visual binary star system 16 Cyg is considered. The only planet in this system has a significant orbital eccentricity, e = 0.69. Previously, Holman et al. suggested the possibility of chaos in the orbital dynamics of the planet due to the proximity of 16 Cyg to the separatrix of the Lidov–Kozai resonance. We have calculated the Lyapunov characteristic exponents on the set of possible orbital parameters for the planet. In all cases, the dynamics of 16 Cyg is regular with a Lyapunov time of more than 30 000 yr. The dynamics is considered in detail for several possible models of the planetary orbit; the dependences of Lyapunov exponents on the time of their calculation and the time dependences of osculating orbital elements have been constructed. Phase space sections for the system dynamics near the Lidov–Kozai resonance have been constructed for all models. A chaotic behavior in the orbital motion of the planet in 16 Cyg is shown to be unlikely, because 16 Cyg in phase space is far from the separatrix of the Lidov–Kozai resonance at admissible orbital parameters, with the chaotic layer near the separatrix being very narrow.     

Unusual rotation modes of minor planetary satellites

An analysis of the character of the possible dynamics of all hitherto known planetary satellites shows two satellites—Amalthea (J5) and Prometheus (S16)—to have the most unusual structure of the phase space of possible rotational motion. These are the only satellites whose phase space of planar rotation may host synchronous resonances of three different kinds: the α resonance, the β resonance, and a mode corresponding to the period doubling bifurcation of the α resonance. We analyze the stability of these states against the tilt of the rotational axis.     

A note concerning the 2:1 and the 3:2 resonances in the asteroid belt

The dynamical behavior of asteroids inside the 2:1 and 3:2 commensurabilities with Jupiter presents a challenge. Indeed most of the studies, either analytical or numerical, point out that the two resonances have a very similar dynamical behavior. In spite of that, the 3:2 resonance, a little outside the main belt, hosts a family of asteroids, called the Hildas, while the 2:1, inside the main belt, is associated to a gap (the Hecuba gap) in the distribution of asteroids.In his search for a dynamical explanation for the Hecuba gap, Wisdom (1987) pointed out the existence of orbits starting with low eccentricity and inclination inside the 2:1 commensurability and going to high eccentricity, and thus to possible encounters with Mars. It has been shown later (Henrard et al.), that these orbits were following a path from the low eccentric belt of secondary resonances to the high eccentric domain of secular resonances. This path crosses a bridge, at moderate inclination and large amplitude of libration, between the two chaotic domains associated with these resonances.The 3:2 resonance being similar in many respects to the 2:1 resonance, one may wonder whether it contains also such a path. Indeed we have found that it exists and is very similar to the 2:1 one. This is the object of the present paper.     

Collective resonant phenomena on small bodies in the solar system

For both asteroids and meteor streams, and also for comets, resonances play a major role for their orbital evolutions but on different time scales. For asteroids both mean motion resonances and secular resonances not only structure the phase space of regular orbits but are mainly at the origin for the inherent chaos of planet crosser objects.For comets and their chaotic routes temporary trapping into orbital resonances is a well known phenomenon. In addition for slow diffusion through the Kuiper belt resonances are the only candidates for originating a slow chaos.Like for asteroids, resonances with Jupiter play a major role for the orbital evolution of meteor streams. Crossing of separatrix like zones appears to be crucial for the formation of arcs and for the dissolution of streams. In particular the orbital inclination of a meteor stream appears to be a critical parameter for arc formation. Numerical results obtained in an other context show that the competition between the Poynting-Robertson drag and the gravitational interaction of grains near the 2/1 resonance might be very important in the long run for the structure of meteor streams.     

Complete spin and orbital evolution of close-in bodies using a Maxwell viscoelastic rheology

In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide the equations of motion both in the body frame and in the inertial frame. With this study, we generalize preexisting models to the spatial case and to arbitrary multipole orders using a formalism taken from quantum theory. We also provide the vectorial expression of the secular tidal torque expanded in Fourier series. Applying this model to close-in exoplanets, we observe that if the relaxation time is longer than the revolution period, the phase space of the system is characterized by the presence of several spin-orbit resonances, even in the circular case. As the system evolves, the planet spin can visit different spin-orbit configurations. The obliquity is decreasing along most of these resonances, but we observe a case where the planet tilt is instead growing. These conclusions derived from the secular torque are successfully tested with numerical integrations of the instantaneous equations of motion on HD 80606 b. Our formalism is also well adapted to close-in super-Earths in multiplanet systems which are known to have non-zero mutual inclinations.     

A web of secondary resonances for large <Emphasis Type="Italic">A</Emphasis>/<Emphasis Type="Italic">m</Emphasis> geostationary debris

The dynamics of space debris with very high A/m near the geostationary orbit is dominated by the gravitational coefficient C ₂₂ and the solar radiation pressure. An analysis of the stability of the orbits by the chaos indicator MEGNO and frequency analysis map FAM shows chaotic layers around the separatrix and reveals a web of sub-structures associated to resonances with the annual period of the Sun. This succession of stable thin islands and chaotic layers can be reproduced and explained by a quite simple toy model, based on a pendulum approach, perturbed, through the eccentricity, by the external (Sun) frequency. The use of suitable action-angle variables in the circulation and libration regions of the pendulum allows to point out new resonances between the geostationary libration angle and the Sun’s longitude. They correspond very well (positions, shape, width) to the structures visible on the FAM representations.     

Long-term evolution of the spin of Mercury: I. Effect of the obliquity and core-mantle friction

The present obliquity of Mercury is very low (less than 0.1°), which led previous studies to always adopt a nearly zero obliquity during the planet’s past evolution. However, the initial orientation of Mercury’s rotation axis is unknown and probably much different than today. As a consequence, we believe that the obliquity could have been significant when the rotation rate of the planet first encountered spin-orbit resonances. In order to compute the capture probabilities in resonance for any evolutionary scenario, we present in full detail the dynamical equations governing the long-term evolution of the spin, including the obliquity contribution.The secular spin evolution of Mercury results from tidal interactions with the Sun, but also from viscous friction at the core-mantle boundary. Here, this effect is also regarded with particular attention. Previous studies show that a liquid core enhances drastically the chances of capture in spin-orbit resonances. We confirm these results for null obliquity, but we find that the capture probability generally decreases as the obliquity increases. We finally show that, when core-mantle friction is combined with obliquity evolution, the spin can evolve into some unexpected configurations as the synchronous or the 1/2 spin-orbit resonance.     

The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincaré variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator Mean Exponential Growth factor of Nearby Orbits (MEGNO), a measure of the maximal Lyapunov exponent, that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best-fitting solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best-fitting solutions are found in dynamically active region of the phase space. The long-term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myr of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behaviour.     

Structure of the 3:1 Jovian resonance: A comparison of numerical methods

The motion of asteroids near the 3:1 Jovian resonance in the restricted planar case is studied using three numerical methods: (a) integrating the full equations of motion, (b) integrating the averaged equations of motion, and (c) using an algebraic mapping recently developed by Wisdom (1982, Astron. J.87, 577–593). The relative merits of each method are investigated. It is concluded that in the regular regions of the phase space, methods b and c give excellent agreement with each other and that provided the maximum eccentricity e_max < 0.4 differences with the exact solution (method a) are <6% in e_max and <27% in the period of the oscillations. The additions of higher order terms in the expansion of the averaged Hamiltonian provides marginally better agreement with the full integration. This is probably due to the slow convergence of the expansion of the disturbing function at large eccentricities (e > 0.3). In chaotic regions of the phase there is little agreement between the orbital elements at any given time calculated by each method. However, all methods reflect the qualitative behavior of the chaotic trajectories and give good agreement on the bounds of the motion. Since the map is at least 200 times faster than solving the full equations of motion it is an efficient method of rapidly exploring accessible regions of the chaotic phase space.     

Normal Form Invariants Around Spin-orbit Periodic Orbits

We consider a model of spin-orbit interaction, describing the motion of an oblate satellite rotating about an internal spin-axis and orbiting about a central planet. The resulting second order differential equation depends upon the parameters provided by the equatorial oblateness of the satellite and its orbital eccentricity. Normal form transformations around the main spin-orbit resonances are carried out explicitly. As an outcome, one can compute some invariants; the fact that these quantities are not identically zero is a necessary condition to prove the existence of nearby periodic orbits (Birkhoff fixed point theorem). Moreover, the nonvanishing of the invariants provides also the stability of the spin-orbit resonances, since it guarantees the existence of invariant curves surrounding the periodic orbit.     

The chaotic rotation of Hyperion

A plot of spin rate versus orientation when Hyperion is at the pericenter of its orbit (surface of section) reveals a large chaotic zone surrounding the synchronous spin-orbit state of Hyperion, if the satellite is assumed to be rotating about a principal axis which is normal to its orbit plane. This means that Hyperion's rotation in this zone exhibits large, essentially random variations on a short time scale. The chaotic zone is so large that it surrounds the ½ and 2 states, and libration in the 3/2 state is not possible. Stability analysis shows that for libration in the synchronous and ½ states, the orientation of the spin axis normal to the orbit plane is unstable, whereas rotation in the 2 state is attitude stable. Rotation in the chaotic zone is also attitude unstable. A small deviation of the principal axis from the orbit normal leads to motion through all angles in both the chaotic zone and the attitude unstable libration regions. Measures of the exponential rate of separation of nearby trajectories in phase space (Lyapunov characteristic exponents) for these three-dimensional motions indicate the the tumbling is chaotic and not just a regular motion through large angles. As tidal dissipation drives Hyperion's spin toward a nearly synchronous value, Hyperion necessarily enters the large chaotic zone. At this point Hyperion becomes attitude unstable and begins to tumble. Capture from the chaotic state into the synchronous or ½ state is impossible since they are also attitude unstable. The 3/2 state does not exist. Capture into the stable 2 state is possible, but improbable. It is expected that Hyperion will be found tumbling chaotically.     

Planetary perturbations on Mercury’s libration in longitude

Two space missions dedicated to Mercury (MESSENGER and BepiColombo) aim at understanding its rotation and confirming the existence of a liquid core. This double challenge requires much more accurate models for the spin-orbit resonant rotation of Mercury. The purpose of this paper is to introduce planetary perturbations on Mercury’s rotation using an analytical method and to analyse the influence of the perturbations on the libration in longitude. Applying a perturbation theory based on the Lie triangle, we were able to re-introduce short periodic terms into the averaged Hamiltonian and to compute the evolution of the rotational variables. The perturbations on Mercury’s forced libration in longitude mainly come from the orbital motion of Mercury (with an amplitude around 41 arcsec that depends on the momenta of inertia). It is completed by various effects from Jupiter (11.86 and 5.93 year-periods), Venus (with a 5.66 year-period), Saturn (14.73 year-period), and the Earth (6.58 year-period). The amplitudes of the oscillations due to Jupiter and Venus are approximately 33% and 10% of those from the orbital motion of Mercury and the amplitudes of the oscillations due to Saturn and the Earth are approximately 3% and 2%. We compare the analytical results with the solution obtained from the spin-orbit numerical model SONYR.