Abstract

The three-body problem is a fundamental concept in celestial mechanics that study the motion of three large entities or heavenly bodies under the influence of gravity. In this case, each body exerts gravitational forces on the other two, creating a complicated interaction of forces that makes calculating their orbits extremely difficult. Unlike the simpler two-body problem, in which only two bodies interact gravitationally, the addition of a third body causes chaotic behaviour, making exact solutions impossible in many circumstances. For understanding of the nature of the three-body problem, it’s essential to grasp some key concepts in orbital mechanics. Newton’s law of universal gravitation states that every particle in the universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely proportional to the square of the distance between their centers. This gravitational force regulates how celestial bodies orbit each other. In a three-body system, such as a planet-moon-star configuration or a binary star system with a distant companion, each body’s motion is influenced not only by the gravitational pull of the other two bodies but also by their changing positions relative to one another over time. This leads to a dynamic and chaotic interaction where the orbits of the bodies can become highly irregular and unpredictable. One of the key challenges in solving the three-body problem lies in its nonlinearity. The gravitational force between the bodies depends on their instantaneous positions and velocities, resulting in a set of coupled differential equations that cannot be easily solved analytically. Even small changes in the initial conditions can lead to vastly different outcomes, a phenomenon known as sensitive dependence on initial conditions, famously illustrated by chaos theory. Because of this sensitivity, three-body systems can exhibit very unpredictable behaviour over time. Orbits may exhibit irregularities, such as precession, oscillations, and even chaotic wandering. While some special cases of the three-body problem can be solved analytically or approximated using simplified models, the general case remains an active area of research in theoretical physics and computational astrodynamics. To study three-body systems, scientists employ a variety of techniques, including numerical simulations, perturbation theory, and statistical analysis. These methods help researchers understand the underlying dynamics of complex orbital systems and make predictions about their behavior over time. On spite of the challenges posed by the three-body problem, we have to explore more deeply in this field which is essential for understanding the mechanics of celestial bodies in the universe.

Introduction

The three-body problem is famous for being impossible to solve. But actually it’s been solved many times, and in ingenious ways. Some of those solutions are incredibly useful, and some are incredibly bizarre. Newton's Principia Physics – and arguably all of science changed forever in 1687 when Isaac Newton published his Principia. Within it were equations of motion and gravity that transformed our erratic-seeming cosmos into a perfectly tuned machine of clockwork predictability. Given the current positions and velocities of the bodies of the solar system, Newton’s equations could be used in principle be used to calculate their locations at any distant time, future or past. I say “in principle” because the reality is not so simple. Despite the beauty of Newton’s equations, they lead to a simple solution for planetary motion in only one case – when two and only two bodies orbit each other sans any other gravitational influence in the universe. Add just one more body and in most cases all motion becomes fundamentally chaotic – there exists no simple solution. This is the three-body problem, and we’ve been trying to solve it for 300 years. What does it mean to find a solution to the three-body problem? Newton’s laws of motion and his law of universal gravitation give us a set of differential equations. In some cases these can be solved with Newton’s other great invention – calculus – to give a simple equation. Plug numbers into that equation and its solved. Those numbers are the starting positions and velocities of your gravitating bodies, plus a value for time. The equations will then give you the state of the system at that time, no matter how far in the past or future. We call such a simple, exactly-solvable equation an analytic expression. That just means it can be written out with a finite number of mathematical operations and functions. In the case of two gravitating bodies, the solutions to Newton’s laws are just the equations for the path traveled by the bodies – be it the parabola of a thrown ball, the circle or ellipse of a planetary orbit, or the hyperbola of an interstellar comet – in general, conic sections – the shapes you get when you slice up a cone.

The Three-Body Problem

These solutions were so simple that Johanne Kepler figured out much about the elliptical solution for planetary motion 70 years before Newton’s laws were even known. And after the Principia was published many sought simple, analytic solutions for more complex systems, with systems of three gravitating bodies being the natural next step. But the additional influence of even a single extra body appeared to make an exact solution impossible. The three body problem became the obsession for many great mathematicians – but over the following three centuries, solutions have been found for very few specialized cases. Why? Well, in the late 1800s, mathematicians Ernst Bruns and Henri Poincaré convincingly asserted that no general analytic solutions exists. The reality of the three-body problem is that the evolution of almost all starting configurations is dominated by chaotic dynamics. Future states are highly dependent on small changes in the initial conditions. Orbits tend towards wild and unpredictable patterns, and almost inevitably one of the bodies is eventually ejected from the system. But despite the apparent hopelessness, there was much profit in learning to predict the gravitational motion of many bodies. For most of the three centuries since Newton, predicting the motion of the planets and the moon was critical for nautical navigation. Now it’s essential to space travel. How do we do it? Well, just because the three body problem for the most part has no useful analytic solution, approximate solutions can be found.

Approximate Solutions

For example, if the bodies are far enough apart then we can approximate a many-body system as a series of two-body systems. For example, each planet of our solar system can be thought of as a separate two-body system with the Sun. That gives you a series of simple elliptical orbits, like those predicted by Kepler. But those orbits eventually shift due to the interactions between the planets. Another useful approximation is when one of the three bodies has a very low mass compared to the other two. We can ignore the minuscule gravitational influence of the smaller body and assume that it moves within the completely solvable two-body orbits of its larger companions. We call this the reduced three-body problem. It works very well for tiny things like artificial satellites around the Earth. It can also be used to approximate the orbits of the moon relative to the Earth and Sun, or the Earth relative to the Sun and Jupiter. These approximate solutions are useful, but ultimately fail to predict perfectly. Even the smallest planetary bodies have some mass, and the solar system as a whole has many massive constituents. The Sun, Jupiter and Saturn alone are automatically a three-body system with no analytic solution, before we even add in the Earth. But the absence of an analytic solution doesn’t mean the absence of any solution. To get an accurate prediction for most three-body systems you need to break the motion of the system into many pieces, and solve them one at a time.

Numerical Integration

A sufficiently small section of any gravitational trajectory can be approximated with an exact, analytical solution – perhaps a straight line or a segment of a two-body path around the center of mass of the entire system, assuming everything else stays fixed. If you break up the problem into tiny enough path segments or time steps, then the small motions of all bodies in the system can be updated step by step. This method of solving differential equations one step at a time is called numerical integration, and when applied to the motion of many bodies it’s an N-body simulation. With modern computers, N-body simulations can accurately predict the motion of the planets into the distant future or solve for millions of objects to simulate the formation and evolution of entire galaxies. But these numerical solutions didn’t begin with the invention of artificial computers. Before that, these calculations had to be done by hand – in fact by many hands. The limitations of approximate solutions, the laboriousness of pre-computer numerical integration, and also the legendary status of the three-body problem inspired generations of physicists and mathematicians to continue to seek exact, analytic solutions. And some succeeded – albeit in very specialized cases.

Euler and Lagrange

The first was Leonhard Euler, who found a family of solutions for three bodies orbiting around a mutual center of mass, where all bodies remain in a straight line – essentially in a permanent eclipse. Joseph-Louis Lagrange found solutions in which the three bodies form an equilateral triangle. In fact, for any two bodies orbiting each other, the Euler and Lagrange’s solutions define 5 additional orbits for a third body that can be described with simple equations. These are the only perfectly analytical solutions to the three body problem that exist. Place a low-mass object on any of these 5 orbits and it will stay there indefinitely, tracking the Earth’s orbit around the Sun. We now call these the Lagrange points, and they’re useful places to park our spacecraft. There was a bit of a gap after Euler and Lagrange because to discover new specialized three-body solutions, we had to search the vast space of possible orbits using computers. The key was to find three-body systems that had periodic motion – they evolve – sometimes in complex ways – back to their starting configuration. In the 70s, Michel Henon and Roger Broucke found a family of solutions involving two masses bouncing back and forth in the center of a third body’s orbit.

The Shape Sphere

In the 90s Cris Moore discovered a stable figure-8 orbit of three equal masses. The numerical discovery of the figure-8 solution was proved mathematically by Alain Chenciner and Richard Montgomery, and insights gained from that proof led to a boom in the discovery of new periodic three-body orbits. Some of these periodic solutions are incredibly complex, but Montgomery came up with a fascinating way to depict them in the absence of simple equations. It’s called the shape-sphere, and it works like this. Imagine the bodies in the 3-body system are the vertices of a triangle, whose center is the center of mass of the system. The evolution of the system can be expressed through the changing shape of that triangle. We throw away certain information — The size of the triangle and its orientation, keeping only information about the relative lengths of the edges, or equivalently the angles between the edges. Now we map that information on the surface of a sphere. We only need the 2-D surface because if we know 2 internal angles of the triangle we also know the 3rd. So, the equator of the sphere represents both angles being zero- that’s a fully collapsed triangle — The 3-bodies are in a straight line, as in Euler’s solutions. The poles are equilateral triangles – so, Lagrange’s solutions. All other orbits move on this sphere as the triangle defined by the orbits evolves. It turns out that the periodic motion on the shape-sphere appears much simpler and easier to analyze than the motion of the bodies themselves. Now hundreds of stable 3-body orbits are known – although it should be noted that besides the Euler and Lagrange solutions, none of these are likely to occur in nature. So their practical use may be limited. Very recently, a new approach to solving the three-body problem has appeared, which transforms the chaotic nature of three-body interactions into a useful tool, rather than a liability. Nicholas Stone and Nathan Leigh published this in Nature in December 2019. The thing about chaotic motion is that the state of the system seems to get randomly shuffled over time. The motion is actually perfectly deterministic – defined between one instant and the next – but can be thought of as approximately random over long intervals. Such a pseudo-random system will, over time, explore all possible configurations consistent with some basic properties like the energy and angular momentum of the system. The system explores what we Call a phase space – a space of possible arrangements of position and velocity. Well, for a pseudo-random system, statistical mechanics lets us calculate the probability of the system being in any part of that phase space at any one time. How is this useful? Well, actually, almost all three-body systems eject one of the bodies, leaving a nice, stable two-body system – a binary pair. Stone and Leigh found that they could identify the regions of phase space where these ejections were likely – and by doing so they could map the range of likely orbital properties for the two objects left behind after the ejection. This looks to be incredibly useful for understanding the evolution of dense regions of the universe, where three-body systems of stars or black holes may form and then disintegrate very frequently. One last thing about the three-body problem. Henri Poincare thought the general case could not be solved. In fact he was wrong. In 1906, not so long after Poincare stern proclamation, finnish mathematician Karl Sundman found a solution to the general three-body problem. It was a converging infinite series that added together an endless chain of terms to solve the orbital calculation. Because the series converged, which successive terms diminished to effectively nothing, so in principle the equation could be written out on paper. However the convergence of Sundman’s series is so slow that it would 10^8 million terms to converge for a typical calculation in celestial mechanics. That is a lot of sheets of paper. So there you have it – the three-body problem is perfectly solved uselessly, or for seemingly useless and bizarre orbits. And it can be approximately solved for all useful and practical purposes – with enough precision to work just fine. Good to know next time you’re in a Chaotic orbit, trying to astronavigate around two other gravitating denizens of space-time.

Acknowledgement

We would like to express our sincere gratitude to all those who have contributed to the completion of this review paper on the 3 Body Problem: Explanation of the Phenomena of Three-Body Problem according to the Orbital Mechanical Behaviour System? Then, we would like to thank the researchers and scholars whose work has paved the way for our understanding of the complexities of the three-body problem in celestial mechanics. Their insights and discoveries have laid the foundation upon which this article is built. We are indebted to the pioneers of celestial mechanics, such as Isaac Newton, Johannes Kepler, Leonhard Euler, Joseph-Louis Lagrange, Ernst Bruns, and Henri Poincaré, whose groundbreaking work has shaped our understanding of gravitational dynamics and orbital mechanics. We would like to extend our appreciation to the authors and researchers whose publications and studies have provided valuable insights and information on the three-body problem, both in classical mechanics and quantum mechanics. Their contributions have enriched our understanding of this fascinating area of study. Then, we would like to acknowledge the contributions of modern computational astrodynamics researchers and numerical analysts, whose advanced techniques and simulations have allowed us to explore the behavior of three-body systems in unprecedented detail. Lastly, we would like to thank the scientific community as a whole for fostering an environment of collaboration, inquiry, and discovery. It is through the collective efforts of scientists, researchers, and educators that we continue to push the boundaries of knowledge and deepen our understanding of the universe. This article would not have been possible without the dedication, expertise, and passion of all those involved. We are grateful for their contributions and honored to have had the opportunity to share their insights with our readers.

.    .    .

References:

• Celletti, A., & Perozzi, E. (2010). Celestial Mechanics: The Waltz of the Planets. Springer Science & Business Media.
• Netchitailo, V. S. (2001). Introduction to the Three-Body Problem. Nova Science Pub Inc.
• Aarseth, S. J. (2003). Gravitational N-Body Simulations: Tools and Algorithms. Cambridge University Press.
• Lehto, O. (1998). The Mathematics of N-body Simulation in Celestial Mechanics. Cambridge University Press.
• Stone, N., & Leigh, N. (2019). Chaotic Dynamics of Three-Body Interactions. Nature, 576(7785), 211-215.
• Wisdom, J., & Holman, M. (1991). Symplectic Maps for the N-Body Problem. The Astronomical Journal, 102(4), 1528-1538.
• Belbruno, E. (2004). Capture Dynamics and Chaotic Motions in Celestial Mechanics: With Applications to the Construction of Low Energy Transfers. Princeton University Press.
• Chenciner, A., & Montgomery, R. (2000). A Remarkable Periodic Solution of the Three-Body Problem in the Case of Equal Masses. Annals of Mathematics, 152(3), 881-901.
• Heggie, D., & Hut, P. (2003). The Gravitational Million-Body Problem: A Multidisciplinary Approach to Star Cluster Dynamics. Cambridge University Press.
• Giordano, B., & Nakanishi, H. (2006). Computational Physics. Pearson/Addison-Wesley.
• Mikkola, S., Innanen, K., & Valtonen, M. (2004). The Three-Body Problem. Celestial Mechanics and Dynamical Astronomy, 90(3-4), 207-211.

I am Aayush Raj Dubey. I pursuing a bachelor’s degree from Dr A.P.J Abdul Kalam Technical University, Lucknow. I am interested in the field of Classical Mechanics, Quantum Mechanics, Celestial Mechanics, Orbital Mechanics, etc as a passion.