
SOME HISTORICAL NOTES
From: Stephen Wolfram, A New Kind of Science Notes for Chapter 7: Mechanisms in Programs and Nature
Section: Chaos Theory and Randomness from Initial Conditions
Page 972
Threebody problem. The twobody problem was analyzed by Johannes Kepler in 1609 and solved by Isaac Newton in 1687. The threebody problem was a central topic in mathematical physics from the mid1700s until the early 1900s. Various exact results were obtained  notably the existence of stable equilateral triangle configurations corresponding to socalled Lagrange points. Many approximate practical calculations, particularly on the EarthMoonSun system, were done using series expansions involving thousands of algebraic terms. (It is now possible to get most results just by direct numerical computation using for example NDSolve.) From its basic setup the threebody system conserves standard mechanical quantities like energy and angular momentum. But it was thought it might also conserve other quantities (or socalled integrals of the motion). In 1887, however, Heinrich Bruns showed that there could be no such quantities expressible as algebraic functions of the positions and velocities of the bodies (in standard Cartesian coordinates). In the mid1890s Henri Poincaré then showed that there could also be no such quantities analytic in positions, velocities and mass ratios. And from these results the conclusion was drawn that the threebody problem could not be solved in terms of algebraic formulas and integrals. In 1912 Karl Sundman did however find an infinite series that could in principle be summed to give the solution  but which converges exceptionally slowly. And even now it remains conceivable that the threebody problem could be solved in terms of more sophisticated standard mathematical functions. But I strongly suspect that in fact nothing like this will ever be possible and that instead the threebody problem will turn out to show the phenomenon of computational irreducibility discussed in Chapter 12 (and that for example threebody systems are universal and in effect able to perform any computation). (See also page 1138.)
In Henri Poincaré’s study of the collection of possible trajectories for threebody systems he identified sensitive dependence on initial conditions (see above), noted the general complexity of what could happen (particularly in connection with socalled homoclinic tangles), and developed topology to provide a simpler overall description. With appropriate initial conditions one can get various forms of simple behavior. The pictures below show some of the possible repetitive orbits of an idealized planet moving in the plane of a pair of stars that are in a perfect elliptical orbit.
The pictures below show results for a fairly typical sequence of initial conditions where all three bodies interact. (The two bodies at the bottom are initially at rest; the body at the top is given progressively larger rightward velocities.) What generically happens is that one of the bodies escapes from the other two (like t or sometimes t^(2/3)). Often this happens quickly, but sometimes all three bodies show complex and apparently random behavior for quite a while. (The delay before escaping is reminiscent of resonant scattering.
Stephen Wolfram, A New Kind of Science (Wolfram Media, 2002), page 972.
© 2002, Stephen Wolfram, LLC

