The web of three-planet resonances in the outer Solar System: II. A source of orbital instability for Uranus and Neptune

The motion of the giant planets from Jupiter to Neptune is chaotic with Lyapunov time of approximately 10 Myr. A recent theory explains the presence of this chaos with three-planet mean-motion resonances, i.e. resonances among the orbital periods of at least three planets. We find that the distribution of these resonances with respect to the semi-major axes of all the planets is compatible with orbital instability. In particular, they overlap in a region of 10⁻³ AU with respect to the variation of the semi-major axes of Uranus and Neptune. Fictitious planetary systems with initial conditions in this region can undergo systematic variations of semi-major axes. The true Solar System is marginally in this region, and Uranus and Neptune undergo very slow systematic variations of semi-major axes with speed of order 10⁻⁴ AU/Gyr.     

A study of the propagation of electromagnetic waves in Titan’s atmosphere with the TLM numerical method

A numerical modeling of the electromagnetic characteristics of Titan’s atmosphere is carried out by means of the TLM numerical method, with the aim of calculating the Schumann resonant frequencies of Saturn’s satellite. The detection and measurement of these resonances by the Huygens probe, which will enter Titan’s atmosphere at the beginning of 2005, is expected to show the existence of electric activity with lightning discharges in the atmosphere of this satellite. As happens with the Schumann frequencies on Earth, losses associated with electric conductivity will make these frequencies lower than theoretically expected, the fundamental frequency being located between 11 and 15 Hz. This numerical study also shows that the strong losses associated to the high conductivity make it impossible for an electromagnetic wave with a frequency of 10 MHz or lower, generated near the surface, to reach the outer part of Titan’s atmosphere.     

Saturn’s wayward shepherds: the peregrinations of Prometheus and Pandora

Saturn’s narrow F ring is flanked by two nearby small satellites, Prometheus and Pandora, discovered in Voyager images taken in 1980 and 1981 (Synnott et al., 1983, Icarus 53, 156-158). Observations with the Hubble Space Telescope (HST) during the ring plane crossings (RPX) of 1995 led to the unexpected finding that Prometheus was ∼19° behind its predicted orbital longitude, based on the Synnott et al. (1983) Voyager ephemeris (Bosh and Rivkin, 1996 Science 272, 518-521; Nicholson et al., 1996, Science 272, 509-515). Whereas Pandora was at its predicted location in August 1995, McGhee (2000, Ph.D. thesis, Cornell University) found from the May and November 1995 RPX data that Pandora also deviates from the Synnott et al. (1983) Voyager ephemeris. Using archival HST data from 1994, previously unexamined RPX images, and a large series of targeted WFPC2 observations between 1996 and 2002, we have determined highly accurate sky-plane positions for Prometheus, Pandora, and nine other satellites found in our images. We compare the Prometheus and Pandora measurements to the predictions of substantially revised and improved ephemerides for the two satellites based on an extensive analysis of a large set of Voyager images (Murray et al., 2000, Bull. Am. Astron. Soc. 32, 1090; Evans, 2001 Ph.D. thesis, Queen Mary College). From December 1994 to December 2000, Prometheus’ orbital longitude lag was changing by −0.71° year⁻¹ relative to the new Voyager ephemeris. In contrast, Pandora is ahead of the revised Voyager prediction. From 1994 to 2000, its longitude offset changed by +0.44° year⁻¹, showing in addition an ∼585 day oscillatory component with amplitude Δλ_CR0 = 0.65 ± 0.07° whose phase matches the expected perturbation due to the nearby 3:2 corotation resonance with Mimas, modulated by the 71-year libration in the longitude of Mimas due to its 4:2 resonance with Tethys. We determine orbital elements for freely precessing equatorial orbits from fits to the 1994-2000 HST observations, from which we conclude that Prometheus’ semimajor axis was 0.31 km larger, and Pandora’s was 0.20 km smaller, than during the Voyager epoch. Subsequent observations in 2001-2002 reveal a new twist in the meanderings of these satellites: Prometheus’ mean motion changed suddenly by an additional −0.77° year⁻¹, equivalent to a further increase in semimajor axis of 0.33 km, at the same time that Pandora’s mean motion changed by +0.92° year⁻¹, corresponding to a change of −0.42 km in its semimajor axis. There is an apparent anticorrelation of the motions of these two moons seen in the 2001-2002 observations, as well as over the 20-year interval since the Voyager epoch. This suggests a common origin for their wanderings, perhaps through direct exchange of energy between the satellites as the result of resonances, possibly involving the F ring.     

Effects of resonances on the stability of retrograde satellites

The stability evolution of family f of the planar circular restricted three-body problem in the Earth–Moon case is explored numerically using the Poincaré surface of section. It is shown that third order resonances are the main cause of the reduction of the stability region of retrograde satellites. Several branches of family f are also computed and these are seen by the configuration of their family characteristics to roughly determine the stability region. Previous results on smaller mass ratios of primaries are thus extended to the Earth–Moon system.     

Exoplanetary systems: The role of an equilibrium at high mutual inclination in shaping the global behavior of the 3-D secular planetary three-body problem

On the basis of a high-order (order 12) expansion of the perturbative potential in powers of the eccentricities and the inclinations, we analyze the secular interactions of two non-coplanar planets which are not in mean-motion resonance. The model is based on the planetary three-body problem which can be reduced to two degrees of freedom by the well-known elimination of the nodes [Jacobi, C.G.J., 1842. Astron. Nachr. XX, 81-102]. We introduce non-singular canonical variables which bring forward the symmetries of the problem. The main dynamical features depend on the location and stability of the equilibria which are easily found with our analytical model. We find that there exists an equilibrium when both eccentricities are zero. When the mutual inclination is small, this equilibrium is stable, but for larger mutual inclination it becomes unstable, generating a large chaotic zone and, by bifurcation, two regular regions, the so-called Kozai resonances. This analytical study which depends on only two parameters (the ratio of the semi-major axes and the mass ratio of the planets) makes possible a large survey of the problem and enables us to identify and quantify its main dynamical features, periodic orbits, regular and chaotic zones, etc. The results of our analytical model are illustrated and confirmed by numerical integrations.     

Normal modes of synchronous rotation

The dynamics of synchronous rotation and physical librations are revisited in order to establish a conceptually simple and general theoretical framework applicable to a variety of problems. Our motivation comes from disagreements between the results of numerical simulations and those of previous theoretical studies, and also because different theoretical studies disagree on basic features of the dynamics. We approach the problem by decomposing the orientation matrix of the body into perfectly synchronous rotation and deviation from the equilibrium state. The normal modes of the linearized equations are computed in the case of a circular satellite orbit, yielding both the periods and the eigenspaces of three librations. Libration in longitude decouples from the other two, vertical modes. There is a fast vertical mode with a period very close to the average rotational period. It corresponds to tilting the body around a horizontal axis while retaining nearly principal-axis rotation. In the inertial frame, this mode appears as nutation and free precession. The other vertical mode, a slow one, is the free wobble. The effects of the nodal precession of the orbit are investigated from the point of view of Cassini states. We test our theory using numerical simulations of the full equations of the dynamics and discuss the disagreements among our study and previous ones. The numerical simulations also reveal that in the case of eccentric orbits large departures from principal-axis rotation are possible due to a resonance between free precession and wobble. We also revisit the history of the Moon's rotational state and show that it switched from one Cassini state to another when it was at 46.2 Earth radii. This number disagrees with the value 34.2 derived in a previous study.     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

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.     

On the oscillations in Mercury's obliquity

Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. Which is the obliquity's dynamical behavior deriving from a complete spin–orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model (acronym of Spin–Orbit N-bodY Relativistic model) integrating the spin–orbit N-body problem applied to the Solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin–orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin–orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin–orbit motion of Mercury in the 2-body problem case (Sun–Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S_2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (S_n case). We find that the remaining amplitude of the oscillations in the S_n case is one order of magnitude larger than in the S_2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state.