Question 1

What is the restricted three-body problem?

Accepted Answer

The restricted three-body problem considers a massless test particle moving under the gravitational influence of two massive primary bodies (m1 and m2) that orbit their common centre of mass (barycentre) in circular orbits. Because the test particle has negligible mass, it does not affect the primaries' motion, greatly simplifying the equations while retaining rich dynamical complexity including chaos.

Question 2

What are Lagrange points L1 through L5?

Accepted Answer

Lagrange points are the five equilibrium positions in the corotating reference frame where the net effective force on a test particle is zero. L1, L2, and L3 lie on the line joining the two primaries and are unstable saddle points. L4 and L5 form equilateral triangles with the primaries and are stable when the mass ratio q = m2/(m1+m2) is less than the Routh critical value (~0.0385), explaining the Trojan asteroids at Jupiter's L4 and L5.

Question 3

What is the Jacobi integral and why is it conserved?

Accepted Answer

The Jacobi integral CJ = 2U*(x,y) - v² is a conserved quantity in the restricted three-body problem, where U* is the effective potential in the rotating frame and v is the particle speed. It is conserved because the rotating frame problem, while not energy-conserving in the usual sense, has a time-independent Hamiltonian in the rotating frame. The Jacobi constant determines which regions of space are accessible to the particle (zero-velocity curves bound the allowed regions).

Question 4

What are zero-velocity curves (Roche equipotentials)?

Accepted Answer

Zero-velocity curves are contours of the effective potential U* = -(x²+y²)/2 - (1-μ)/r1 - μ/r2 at which a particle with given Jacobi constant CJ would have zero kinetic energy (v=0). They define the boundaries of regions forbidden to the particle. As CJ decreases, the forbidden regions shrink; when CJ equals the potential at L1, the inner Roche lobe opens and mass transfer between the two primaries becomes possible — a key mechanism in binary star evolution.

Question 5

What is a Roche lobe in binary star systems?

Accepted Answer

A Roche lobe is the region around each star in a binary system within which material is gravitationally bound to that star in the corotating frame. The two Roche lobes meet at the L1 Lagrange point. When an evolving star expands to fill its Roche lobe, mass spills through L1 onto its companion — a process called Roche lobe overflow. This drives mass transfer in cataclysmic variables, X-ray binaries, and Type Ia supernova progenitors.

Question 6

Why is the three-body problem chaotic?

Accepted Answer

The three-body problem exhibits sensitive dependence on initial conditions (chaos) because the equations of motion are nonlinear and have no general closed-form solution. Even in the restricted case, trajectories near the Lagrange saddle points L1, L2, L3 are exponentially unstable. Small differences in initial position or velocity lead to exponentially diverging trajectories over time, characteristic of positive Lyapunov exponents. This is why long-term prediction of asteroid and comet orbits near resonances is fundamentally limited.

Question 7

What integration method is used and why RK4?

Accepted Answer

This simulation uses the classical fourth-order Runge-Kutta (RK4) method. RK4 achieves O(h⁵) local truncation error per step, making it significantly more accurate than simpler Euler or RK2 methods for the same step size. The Jacobi integral serves as a built-in accuracy check: a well-integrated trajectory should conserve CJ to within a small tolerance. For long-term orbital studies, symplectic integrators (like Verlet or Yoshida) would be preferred because they conserve a modified Hamiltonian exactly.

Question 8

What is the mass ratio parameter μ (mu)?

Accepted Answer

The mass ratio parameter μ = m2/(m1+m2) is the only free parameter in the restricted three-body problem (in units where total mass and orbital separation equal 1). It ranges from 0 (negligible secondary) to 0.5 (equal masses). Real examples: Earth-Moon system μ ≈ 0.0121, Sun-Jupiter μ ≈ 0.000953, Sun-Earth μ ≈ 3×10⁻⁶. Equal-mass binaries (μ = 0.5) have symmetric Roche lobes.

Question 9

How does the corotating (synodic) reference frame work?

Accepted Answer

The corotating (synodic) frame rotates at the same angular velocity ω as the two primaries, so the primaries appear stationary in this frame. In exchange, two fictitious forces appear: the centrifugal force (−ω×(ω×r)) and the Coriolis force (−2ω×v). These are absorbed into the effective potential U* and the Coriolis coupling terms (ẍ−2ẏ and ÿ+2ẋ) in the equations of motion. The Lagrange equilibria are precisely the points where gravitational and centrifugal forces exactly cancel.

Question 10

What real-world systems does this model describe?

Accepted Answer

The restricted three-body problem describes: (1) asteroid motion in the Sun-Jupiter system, explaining Kirkwood gaps and Trojan asteroids; (2) spacecraft trajectory design using libration point orbits (JWST at Sun-Earth L2, SOHO at L1); (3) mass transfer in binary stars including cataclysmic variables and X-ray binaries; (4) the stability of moons around planets in binary asteroid systems; and (5) tidal disruption of stellar clusters passing near galactic centres.

🪐 Restricted Three-Body Problem — Roche Lobes

About this simulation

🔬 What it shows

🎮 How to use

💡 Did you know?

Frequently asked questions

What did you find?