2.2 Two-body and many-body problems

Newtonian Mechanics in Celestial Systems

Celestial mechanics is built on Newtonian gravity: a set of rules that lets you predict how objects move under gravitational attraction. The two-body problem is the cleanest version of this, and it's one of the few gravitational problems you can solve exactly. Once you add a third body or more, exact solutions vanish and you need new strategies.

Two-body problem setup and solution

The two-body problem asks: given two masses interacting only through gravity, how do they move? The key insight is that you can reduce this to an equivalent one-body problem using the concept of reduced mass.

Starting from Newton's laws:

• First law: Objects stay at rest or in uniform motion unless a force acts on them.
• Second law: $F = ma$ connects net force to acceleration.
• Third law: Forces between two bodies are equal in magnitude and opposite in direction.

Gravitational force between two masses is given by Newton's law of universal gravitation:

$F = G\frac{m_1 m_2}{r^2}$

where $G$ is the gravitational constant ($6.674 \times 10^{-11} \, \text{N m}^2 \text{kg}^{-2}$) and $r$ is the distance between the two masses.

Reducing to one body:

To solve the two-body problem, you transform it into an equivalent single-body problem. Here's how:

1. Define the center of mass of the system. In the center-of-mass frame, the total momentum is zero, which simplifies the math considerably.

2. Introduce the relative position vector $\vec{r} = \vec{r}_2 - \vec{r}_1$, which tracks the separation between the two bodies.

3. Define the reduced mass: $\mu = \frac{m_1 m_2}{m_1 + m_2}$. This lets you treat the problem as a single particle of mass $\mu$ orbiting a fixed center.

4. The equation of motion for the relative coordinate becomes:

$\frac{d^2\vec{r}}{dt^2} = -G(m_1 + m_2)\frac{\vec{r}}{r^3}$

This is mathematically identical to a single particle moving in a central force field, which is a problem with a known solution.

Conservation laws constrain the solution:

• Energy conservation fixes the total mechanical energy (kinetic + gravitational potential). The sign of the total energy determines the orbit type.
• Angular momentum conservation forces the motion into a single plane. This is why orbits are flat, not three-dimensional spirals.

Orbit types depend on total energy $E$:

• $E < 0$: Elliptical orbit (bound). Planets around the Sun follow these.
• $E = 0$: Parabolic trajectory (marginally unbound). The object has exactly escape velocity.
• $E > 0$: Hyperbolic path (unbound). Spacecraft flybys and interstellar objects like 'Oumuamua follow these.

All three are conic sections, unified by the same gravitational equation. The eccentricity $e$ distinguishes them: $e < 1$ for ellipses, $e = 1$ for parabolas, $e > 1$ for hyperbolas.

Challenges of the many-body problem

Once you go beyond two bodies, the problem changes fundamentally. There is no general closed-form solution for three or more gravitationally interacting bodies. This was proven in the context of the three-body problem and remains one of the central facts of celestial mechanics.

Why it's so hard:

The equations of motion become a coupled set of nonlinear differential equations. Each body's acceleration depends on the positions of all other bodies, and those positions are themselves changing due to the same mutual forces. You can't decouple the equations the way you can with two bodies.

Chaos enters the picture. Many-body gravitational systems are generically chaotic, meaning tiny differences in initial conditions grow exponentially over time. For the solar system, this means orbital predictions become unreliable beyond roughly 50-100 million years, even though the system appears stable on human timescales.

Approximate approaches:

• Perturbation theory treats the extra bodies as small disturbances to a known two-body solution. This works well when one mass dominates (like the Sun in our solar system) and the perturbations are small, but it breaks down for strong interactions.
• The restricted three-body problem is a useful special case: two massive bodies orbit each other, and a third body with negligible mass moves in their combined gravitational field. This setup gives rise to the five Lagrange points and is used to study satellite orbits and Trojan asteroids (like Jupiter's Trojans at the L4 and L5 points).

Stability in many-body systems depends on the system's architecture:

• Hierarchical systems, where orbits are nested at very different scales (e.g., a moon orbiting a planet orbiting a star), tend to be stable over long periods.
• Resonant configurations can either stabilize orbits by preventing close encounters or destabilize them by pumping up eccentricities. More on this below.

Numerical methods for gravitational systems

Since analytic solutions don't exist for most many-body problems, astrophysicists rely on numerical simulation. The goal is to integrate the equations of motion forward in time, step by step, to approximate the true trajectories.

Main simulation approaches:

• Direct N-body integration computes the gravitational force between every pair of particles at each timestep. This is the most accurate method but scales as $O(N^2)$, making it expensive for large $N$.
• Tree codes (like the Barnes-Hut algorithm) group distant particles together and approximate their combined gravitational effect. This reduces the scaling to roughly $O(N \log N)$.
• Particle-mesh methods map mass onto a grid and solve for the gravitational potential using fast Fourier transforms. These are efficient for large-scale cosmological simulations but sacrifice resolution at small scales.

Keeping simulations accurate over long times:

• Symplectic integrators are designed specifically for Hamiltonian systems (which gravitational systems are). They don't conserve energy exactly at each step, but they prevent systematic energy drift over millions of orbits. This makes them the standard choice for long-term solar system integrations.
• Adaptive time-stepping shrinks the timestep during close encounters (when forces change rapidly) and expands it during calm phases. Without this, you'd either waste computation on quiet periods or miss the physics during close passes.
• Error monitoring tracks conserved quantities like total energy and angular momentum. If these drift beyond acceptable thresholds, something has gone wrong with the integration.

Applications across scales:

Resonance in orbital dynamics

Orbital resonance occurs when two orbiting bodies exert periodic gravitational influence on each other because their orbital periods form a ratio of small integers. These repeated nudges, always occurring at the same points in the orbits, can accumulate over time and dramatically reshape orbital configurations.

Types of resonance:

• Mean motion resonances involve the orbital periods themselves. If body A orbits exactly twice for every one orbit of body B, they're in a 2:1 mean motion resonance.
• Secular resonances involve the slow precession rates of orbits rather than the orbital periods. The $\nu_6$ secular resonance at the inner edge of the asteroid belt, for example, is set by the precession rate of Saturn's orbit and is responsible for clearing asteroids from that region.

Common resonance ratios observed in nature include 2:1, 3:2, and 4:3. These appear across planetary systems, moon systems, and even exoplanetary architectures.

What resonances do to orbits:

Resonances can pump up eccentricity, making orbits more elongated, or alter inclination, tilting orbital planes. Whether this stabilizes or destabilizes a system depends on the specific resonance.

Resonance capture happens when orbits evolve slowly (through migration or tidal effects) and bodies drift into resonant configurations where they become locked. This is thought to have occurred during the early solar system when the giant planets migrated through the protoplanetary disk.

Notable examples:

• Jupiter's Galilean moons (Io, Europa, Ganymede) are locked in a 4:2:1 Laplace resonance. For every four orbits of Io, Europa completes two and Ganymede completes one. This resonance drives tidal heating in Io, making it the most volcanically active body in the solar system.
• TRAPPIST-1 hosts seven Earth-sized planets, several of which form a near-resonant chain. This architecture likely formed through convergent migration in the protoplanetary disk.
• Kepler-223 contains four planets in a resonant chain, providing evidence that resonance capture during migration is a common process in planet formation.