Unraveling the Three-Body Problem: From Newton’s Laws to Numerical Solutions

In the 17th century Isaac Newton formulated the law of universal gravitation and the laws of motion, providing exact solutions for the two-body problem (e.g., Earth orbiting the Sun). Introducing a third body makes the problem extremely complex, and no general analytical solution exists.

Mathematical Description of the Three-Body Problem

Basic assumptions

Three point masses \(m_1\), \(m_2\), and \(m_3\).

Only mutual gravitational forces act between them.

No external forces or relativistic effects are considered.

Law of universal gravitation

The gravitational force between any two masses is

F_{ij} = G * m_i * m_j / r_{ij}^2 * \hat{r}_{ij}

where \(G\) is the gravitational constant, \(r_{ij}\) is the distance between the masses, and \(\hat{r}_{ij}\) is the unit vector pointing from one mass to the other.

Equations of motion

According to Newton’s second law, the acceleration of each mass is the sum of the gravitational forces from the other two masses, leading to three vector equations.

Derivation Process

Define position vectors : Let \(\mathbf{r}_i\) be the position vector of mass \(i\).

Compute interaction forces : For mass \(i\), the force from mass \(j\) is given by the formula above, and similarly for the third mass.

Form the system of differential equations : Each mass has three coordinates, giving nine second‑order equations. Converting each second‑order equation to two first‑order equations yields a total of eighteen first‑order ordinary differential equations.

These equations can be written in vector form as a state‑vector \(\mathbf{y}\) of 18 components, with \(\dot{\mathbf{y}} = \mathbf{f}(\mathbf{y}, t)\).

Numerical solution methods

Euler method : simple but low accuracy.

Runge‑Kutta method : commonly used high‑accuracy integrator.

Symplectic (e.g., Verlet) method : structure‑preserving algorithm suitable for Hamiltonian systems.

Example trajectories obtained by numerical integration illustrate the chaotic and unpredictable nature of the three-body problem.

The mathematical model of the three-body problem shows that even under simple physical laws, system behavior can become extremely complex. Although a general analytical solution is unavailable, numerical methods and studies of special cases continually deepen our understanding of this classic problem.