Planetary systems are fundamentally chaotic—minute perturbations can lead to radically different outcomes over billion-year timescales. Stability depends on orbital spacing, mass ratios, and proximity to resonances. Some systems remain stable indefinitely; others destabilize on Gyr timescales, triggering planet ejections or collisions. Numerical simulations and analytical stability criteria predict long-term evolution.
From your work on N-body planetary dynamics, you know that even three gravitationally interacting bodies have no general closed-form solution—orbits must be computed numerically, and small differences in initial conditions can produce wildly divergent outcomes. This sensitivity to initial conditions is the hallmark of chaos, and it pervades real planetary systems. Our own solar system is chaotic: Mercury has a small but nonzero probability of colliding with Venus over the next several billion years, and the inner planets' orbits are predictable only to about 50–100 million years into the future. The question of planetary system stability is therefore not "is this system stable forever?" but rather "how long before instability manifests, and what triggers it?"
The primary analytical tool for assessing stability is the concept of Hill stability and its extensions. Two adjacent planets are Hill-stable if their orbits cannot cross—meaning no collision or close encounter is possible regardless of how the eccentricities evolve. The criterion depends on the planets' masses relative to the star and their orbital separation measured in units of their mutual Hill radius, R_H = a[(m₁ + m₂)/(3M*)]^(1/3), where a is the semi-major axis and M* is the stellar mass. Empirically, systems with spacing greater than about 3.5 mutual Hill radii tend to be long-term stable, while tighter configurations are vulnerable. From your study of multi-planet system architecture, you know that observed exoplanet systems cluster near this stability boundary—they are packed as tightly as dynamical stability allows.
Mean-motion resonances play a dual role in stability. When orbital periods form simple integer ratios (2:1, 3:2, 5:3), the resulting periodic gravitational kicks can either stabilize or destabilize orbits depending on the geometry. Resonance capture during planetary migration can lock planets into stable configurations—as seen in the TRAPPIST-1 system, where seven planets maintain a resonant chain. But resonances can also pump eccentricities to orbit-crossing values, as Jupiter's resonances do to asteroids in the Kirkwood gaps. Whether a resonance stabilizes or destabilizes depends on the relative phases of the planets (the resonant arguments) and whether energy dissipation mechanisms maintain the lock.
In practice, long-term stability is assessed through large suites of numerical integrations. Researchers run thousands of simulations with slightly varied initial conditions and track which fraction survive for billions of orbits. Instability typically manifests suddenly after a long quiescent period—eccentricities grow slowly through secular interactions until orbits finally cross, triggering a rapid cascade of close encounters, ejections, or collisions within a few thousand years. This pattern explains why the solar system appears orderly today despite being formally chaotic: we may simply be living during the quiescent phase of a system that will eventually destabilize. Understanding these timescales is essential for interpreting the architectures of observed exoplanet systems and reconstructing the dynamical histories that produced them.