The Restricted Three-Body Problem: Chaos in Celestial Mechanics

The setup: a co-rotating frame and one mass parameter

The restricted three-body problem asks how a massless test particle moves under the gravity of two much larger bodies, the primaries, which themselves orbit their common centre of mass on fixed circular paths. Because the third body is too light to perturb the primaries, we can work in a co-rotating reference frame that spins at the same rate as their orbit — the two primaries sit still, at the cost of introducing centrifugal and Coriolis pseudo-forces. In normalised units the whole system reduces to a single mass parameter μ = m₂/(m₁+m₂).

ẍ − 2ẏ = ∂Ω/∂x        ẏ̈ + 2ẋ = ∂Ω/∂y
Ω(x,y) = ½(x²+y²) + (1−μ)/r₁ + μ/r₂   (effective potential)

These coupled, non-linear equations have no general solution in elementary functions. Henri Poincaré's study of exactly this system in the late nineteenth century revealed that trajectories can be enormously sensitive to their starting conditions — the historical birthplace of chaos theory, and still a cornerstone of modern celestial mechanics.

The Jacobi constant and zero-velocity curves

Although energy in the usual sense is not conserved in the rotating frame, the system does have one conserved quantity: the Jacobi constant C = 2Ω − v². Because C never changes along a trajectory, it confines the small body to regions where v² ≥ 0. The boundaries of these regions are the zero-velocity curves — invisible walls the body can approach but never cross at a given energy, and mission designers exploit them to plot fuel-efficient low-energy transfer routes.

The five Lagrange points

The effective potential Ω has five stationary points where its gradient vanishes — the Lagrange points L1 to L5. The three collinear points, lying on the line joining the primaries, are dynamically unstable, like a ball balanced on a saddle. The two triangular points L4 and L5 form equilateral triangles with the primaries and are stable whenever the primaries' mass ratio exceeds roughly 24.96 to 1 — the reason thousands of Trojan asteroids cluster around Jupiter's L4 and L5. Real missions live in this geometry: the James Webb Space Telescope orbits the Sun–Earth L2 in a halo orbit, and SOHO sits near Sun–Earth L1 for an uninterrupted view of the Sun.

Chaos is deterministic, not random

Two trajectories that begin a hair's breadth apart in phase space diverge exponentially, a property measured by a positive Lyapunov exponent. The motion is never random — the same starting state always yields the same path — but because any real initial condition is known only approximately, long-term prediction degrades rapidly. Regular, quasi-periodic orbits and wildly chaotic ones can coexist in the same system, separated by delicate boundaries mapped with Poincaré sections. The three-body problem does have infinitely many trajectories, and even special exact solutions such as the figure-eight orbit — what it lacks is a single general formula.

Frequently asked questions

Try it live

Everything above runs in your browser — open Restricted Three-Body Problem — Roche Lobes, drag the mass ratio μ and Jacobi constant C_J sliders, and click the canvas to launch new particles and watch chaos emerge. Nothing is installed, nothing is uploaded.