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