Resonance in Planetary Orbits

Stability: We define stability within the solar system as a state wherein the orbits of bodies remain well separated, and the planets will remain bound to the Sun for an infinite amount of time.

Chaos: Two trajectories that are arbitrarily close in phase space will diverge exponentially with time. In a chaotic system, the timescale for divergence is independent of the precise values of the initial conditions. We can characterize the deviation between the two trajectories as d, where

The scaling variable c is the Lyapunov coefficient, and we define its inverse to be the Lyapunov timescale, tc. Strictly speaking, the system is so sensitive to initial conditions that detailed, long-term behavior is lost within several iterations of this characteristic timescale. We expect a 100% discrepancy within 20 tc for a change of 10-10 in initial conditions. For Terrestrial planets, tc is roughly 5 Myr, a relatively short amount of time that implies significant disorder within 100 Myr (contrast this with tc = 20 Myr for Jupiter). The timescale for large changes in principal orbital elements, however, is often many orders of magnitude larger than tc. Furthermore, chaos can be strongly constrained by patterns of resonance.

Consider the case of a simple one-dimensional harmonic oscillator - our friend the spring. We can express the equations of motion as

where a force Fd drives a mass m which oscillates at frequency o, at a driving frequency d.

This can produce a large-amplitude, long period response if o is roughly equal to d, even if the amplitude of the driving force is small. For the case o = d, the motion is quite different.

This resonant driving force can produce secular (steadily increasing, not periodic) growth. The simplest resonances to visualize are mean motion resonances, for which the orbital periods of two bodies form a ratio of the form

where N is an integer. Well-known examples of resonances found in the solar system include the following.

The effect of one body on a neighboring object can be calculated by evaluating the potential as the sum of (1) the Keplerian motion about the Sun, (2) a disturbing function Rdist made up of direct terms for the pairwise interaction, and (3) indirect terms associated with the back-reaction of the bodies on the Sun. For the case of two bodies of mass much less than their primary,

The rough limit to the range of influence for a planet, or a moon, of mass m2 in orbit around a primary of mass m1 is defined as the extent of the Hill sphere:

A test particle at the edge of the Hill sphere will experience a gravitational force from the planet equal in magnitude to the difference in the tidal force from the Sun on the planet and on the test particle. It extends to the radius of the L1 Lagrange Point, and in the limit where m2 << m1 it is equivalent to the Roche lobe (known to us in the context of binary star mass transfer).