N-body problem

Knowing three orbital positions of a planet's orbit – positions obtained by Sir Isaac Newton from astronomer John Flamsteed – Newton was able to produce an equation by straightforward analytical geometry, to predict a planet's motion; i.e., to give its orbital properties: position, orbital diameter, period and orbital velocity. Having done so, he and others soon discovered over the course of a few years, those equations of motion did not predict some orbits correctly or even very well. Newton realized that this was because gravitational interactive forces amongst all the planets were affecting all their orbits. The aforementioned revelation strikes directly at the core of what the n-body issue physically is: as Newton understood, it is not enough to just provide the beginning location and velocity, or even three orbital positions, in order to establish a planet's actual orbit; one must also be aware of the gravitational interaction forces. Thus came the awareness and rise of the -body "problem" in the early 17th century. These gravitational attractive forces do conform to Newton's laws of motion and to his law of universal gravitation, but the many multiple (-body) interactions have historically made any exact solution intractable. Ironically, this conformity led to the wrong approach. After Newton's time the -body problem historically was not stated correctly because it did not include a reference to those gravitational interactive forces. Newton does not say it directly but implies in his Principia the -body problem is unsolvable because of those gravitational interactive forces. Newton said) The version finally printed contained many important ideas which led to the development of chaos theory. The problem as stated originally was finally solved by Karl Fritiof Sundman for and generalized to by L. K. Babadzanjanz and Qiudong Wang. == General formulation ==

The -body problem considers point masses in an inertial reference frame in three dimensional space \mathbb{R}^3 moving under the influence of mutual gravitational attraction. Each mass has a position vector . Newton's second law says that mass times acceleration is equal to the sum of the forces on the mass. Newton's law of gravity says that the gravitational force felt on mass by a single mass is given by \mathbf{F}_{ij} = \frac{G m_i m_j}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^2} \cdot \frac{\left(\mathbf{q}_j - \mathbf{q}_i\right)}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|} = \frac{G m_i m_j \left(\mathbf{q}_j - \mathbf{q}_i\right)}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^3}, where is the gravitational constant and is the magnitude of the distance between and (metric induced by the norm). Summing over all masses yields the -body equations of motion:{{Equation box 1|cellpadding|border|indent=:|equation= m_i \frac{d^2\mathbf{q}_i}{dt^2} = \sum_{j=1 \atop j \ne i}^n \frac{G m_i m_j \left(\mathbf{q}_j - \mathbf{q}_i\right)}{\left\| \mathbf{q}_j - \mathbf{q}_i\right\|^3} = -\frac{\partial U}{\partial \mathbf{q}_i} |border colour=#0073CF|background colour=#F5FFFA}}where is the self-potential energy U = -\sum_{1 \le i Defining the momentum to be , Hamilton's equations of motion for the -body problem become \frac{d\mathbf{q}_i}{dt} = \frac{\partial H}{\partial \mathbf{p}_i} \qquad \frac{d\mathbf{p}_i}{dt} = -\frac{\partial H}{\partial \mathbf{q}_i}, where the Hamiltonian function is H = T + U and is the kinetic energy T = \sum_{i=1}^n \frac{\left\| \mathbf{p}_i \right\|^2}{2m_i}. Hamilton's equations show that the -body problem is a system of first-order differential equations, with initial conditions as initial position coordinates and initial momentum values. Symmetries in the -body problem yield global integrals of motion that simplify the problem. Translational symmetry of the problem results in the center of mass \mathbf{C} = \frac{\displaystyle\sum_{i=1}^n m_i \mathbf{q}_i}{\displaystyle\sum_{i=1}^n m_i} moving with constant velocity, so that , where is the linear velocity and is the initial position. The constants of motion and represent six integrals of the motion. Rotational symmetry results in the total angular momentum being constant \mathbf{A} = \sum_{i=1}^n \mathbf{q}_i \times \mathbf{p}_i, where × is the cross product. The three components of the total angular momentum yield three more constants of the motion. The last general constant of the motion is given by the conservation of energy . Hence, every -body problem has ten integrals of motion. Because and are homogeneous functions of degree 2 and −1, respectively, the equations of motion have a scaling invariance: if is a solution, then so is for any . The moment of inertia of an -body system is given by I = \sum_{i=1}^n m_i \mathbf{q}_i \cdot \mathbf{q}_i = \sum_{i=1}^n m_i \left\|\mathbf{q}_i\right\|^2 and the virial is given by . Then the Lagrange–Jacobi formula states that \frac{d^2I}{dt^2} = 2T - U. For systems in dynamic equilibrium, the longterm time average of is zero. Then on average the total kinetic energy is half the total potential energy, , which is an example of the virial theorem for gravitational systems. If is the total mass and a characteristic size of the system (for example, the radius containing half the mass of the system), then the critical time for a system to settle down to a dynamic equilibrium is t_\mathrm{cr} = \left(\frac{GM}{R^3}\right)^{-1/2}. == Special cases ==

Two-body problem Any discussion of planetary interactive forces has always started historically with the two-body problem. The purpose of this section is to relate the real complexity in calculating any planetary forces. Note in this Section also, several subjects, such as gravity, barycenter, Kepler's Laws, etc.; and in the following Section too (Three-body problem) are discussed on other Wikipedia pages. Here though, these subjects are discussed from the perspective of the -body problem. The two-body problem () was first solved by Isaac Newton in 1687 using geometric methods, but a complete solution was given in 1710 by Johann Bernoulli (1667–1748) by classical theory by assuming the main point-mass was fixed; this is outlined here. Consider then the motion of two bodies, say the Sun and the Earth, with the Sun fixed, then: \begin{align} m_1 \mathbf{a}_1 &= \frac{Gm_1m_2}{r_{12}^3}(\mathbf{r}_2-\mathbf{r}_1) &&\quad\text{Sun–Earth} \\ m_2 \mathbf{a}_2 &= \frac{Gm_1m_2}{r_{21}^3}(\mathbf{r}_1-\mathbf{r}_2) &&\quad\text{Earth–Sun} \end{align} The equation describing the motion of mass relative to mass is readily obtained from the differences between these two equations and after canceling common terms gives: \mathbf{\alpha} + \frac{\eta}{r^3} \mathbf{r} = \mathbf{0} Where • is the vector position of relative to ; • is the Eulerian acceleration ; • . The equation is the fundamental differential equation for the two-body problem Bernoulli solved in 1734. Notice for this approach forces have to be determined first, then the equation of motion resolved. This differential equation has elliptic, or parabolic or hyperbolic solutions. It is incorrect to think of (the Sun) as fixed in space when applying Newton's law of universal gravitation, and to do so leads to erroneous results. The fixed point for two isolated gravitationally interacting bodies is their mutual barycenter, and this two-body problem can be solved exactly, such as using Jacobi coordinates relative to the barycenter. Dr. Clarence Cleminshaw calculated the approximate position of the Solar System's barycenter, a result achieved mainly by combining only the masses of Jupiter and the Sun. Science Program stated in reference to his work: The Sun wobbles as it rotates around the Galactic Center, dragging the Solar System and Earth along with it. What mathematician Kepler did in arriving at his three famous equations was curve-fit the apparent motions of the planets using Tycho Brahe's data, and not curve-fitting their true circular motions about the Sun (see figure). Both Robert Hooke and Newton were well aware that Newton's Law of Universal Gravitation did not hold for the forces associated with elliptical orbits. In fact, Newton's Universal Law does not account for the orbit of Mercury, the asteroid belt's gravitational behavior, or Saturn's rings. Newton stated (in section 11 of the Principia) that the main reason, however, for failing to predict the forces for elliptical orbits was that his math model was for a body confined to a situation that hardly existed in the real world, namely, the motions of bodies attracted toward an unmoving center. Some present physics and astronomy textbooks do not emphasize the negative significance of Newton's assumption and end up teaching that his mathematical model is in effect reality. It is to be understood that the classical two-body problem solution above is a mathematical idealization. See also Kepler's first law of planetary motion. Three-body problem This section relates a historically important -body problem solution after simplifying assumptions were made. In the past not much was known about the -body problem for . The case has been the most studied. Many earlier attempts to understand the three-body problem were quantitative, aiming at finding explicit solutions for special situations. • In 1687, Isaac Newton published in the Principia the first steps in the study of the problem of the movements of three bodies subject to their mutual gravitational attractions, but his efforts resulted in verbal descriptions and geometrical sketches; see especially Book 1, Proposition 66 and its corollaries (Newton, 1687 and 1999 (transl.), see also Tisserand, 1894). • In 1767, Euler found collinear motions, in which three bodies of any masses move proportionately along a fixed straight line. The Euler's three-body problem is the special case in which two of the bodies are fixed in space (this should not be confused with the circular restricted three-body problem, in which the two massive bodies describe a circular orbit and are only fixed in a synodic reference frame). • In 1772, Lagrange discovered two classes of periodic solution, each for three bodies of any masses. In one class, the bodies lie on a rotating straight line. In the other class, the bodies lie at the vertices of a rotating equilateral triangle. In either case, the paths of the bodies will be conic sections. Those solutions led to the study of central configurations, for which for some constant . • A major study of the Earth–Moon–Sun system was undertaken by Charles-Eugène Delaunay, who published two volumes on the topic, each of 900 pages in length, in 1860 and 1867. Among many other accomplishments, the work already hints at chaos, and clearly demonstrates the problem of so-called "small denominators" in perturbation theory. • In 1917, Forest Ray Moulton published his now classic, An Introduction to Celestial Mechanics (see references) with its plot of the restricted three-body problem solution (see figure below). An aside, see Meirovitch's book, pages 413–414 for his restricted three-body problem solution. Moulton's solution may be easier to visualize (and definitely easier to solve) if one considers the more massive body (such as the Sun) to be stationary in space, and the less massive body (such as Jupiter) to orbit around it, with the equilibrium points (Lagrangian points) maintaining the 60° spacing ahead of, and behind, the less massive body almost in its orbit (although in reality neither of the bodies are truly stationary, as they both orbit the center of mass of the whole system—about the barycenter). For sufficiently small mass ratio of the primaries, these triangular equilibrium points are stable, such that (nearly) massless particles will orbit about these points as they orbit around the larger primary (Sun). The five equilibrium points of the circular problem are known as the Lagrangian points. See figure below: In the restricted three-body problem math model figure above (after Moulton), the Lagrangian points L4 and L5 are where the Trojan planetoids resided (see Lagrangian point); is the Sun and is Jupiter. L2 is a point within the asteroid belt. It has to be realized for this model, this whole Sun-Jupiter diagram is rotating about its barycenter. The restricted three-body problem solution predicted the Trojan planetoids before they were first seen. The -circles and closed loops echo the electromagnetic fluxes issued from the Sun and Jupiter. It is conjectured, contrary to Richard H. Batin's conjecture (see References), the two are gravity sinks, in and where gravitational forces are zero, and the reason the Trojan planetoids are trapped there. The total amount of mass of the planetoids is unknown. The restricted three-body problem assumes the mass of one of the bodies is negligible. For a discussion of the case where the negligible body is a satellite of the body of lesser mass, see Hill sphere; for binary systems, see Roche lobe. Specific solutions to the three-body problem result in chaotic motion with no obvious sign of a repetitious path. The restricted problem (both circular and elliptical) was worked on extensively by many famous mathematicians and physicists, most notably by Poincaré at the end of the 19th century. Poincaré's work on the restricted three-body problem was the foundation of deterministic chaos theory. In the restricted problem, there exist five equilibrium points. Three are collinear with the masses (in the rotating frame) and are unstable. The remaining two are located on the third vertex of both equilateral triangles of which the two bodies are the first and second vertices. Four-body problem Inspired by the circular restricted three-body problem, the four-body problem can be greatly simplified by considering a smaller body to have a small mass compared to the other three massive bodies, which in turn are approximated to describe circular orbits. This is known as the bicircular restricted four-body problem (also known as bicircular model) and it can be traced back to 1960 in a NASA report written by Su-Shu Huang. This formulation has been highly relevant in the astrodynamics, mainly to model spacecraft trajectories in the Earth-Moon system with the addition of the gravitational attraction of the Sun. The former formulation of the bicircular restricted four-body problem can be problematic when modelling other systems than the Earth-Moon-Sun, so the formulation was generalized by Negri and Prado to expand the application range and improve the accuracy without loss of simplicity. Planetary problem The planetary problem is the -body problem in the case that one of the masses is much larger than all the others. A prototypical example of a planetary problem is the Sun–Jupiter–Saturn system, where the mass of the Sun is about 1000 times larger than the masses of Jupiter or Saturn. In 1963, Vladimir Arnold proved using KAM theory a kind of stability of the planetary problem: there exists a set of positive measure of quasiperiodic orbits in the case of the planetary problem restricted to the plane. Central configurations A central configuration is an initial configuration such that if the particles were all released with zero velocity, they would all collapse toward the center of mass . Central configurations have played an important role in understanding the topology of invariant manifolds created by fixing the first integrals of a system. -body choreography Solutions in which all masses move on the same curve without collisions are called choreographies. A choreography for was discovered by Lagrange in 1772 in which three bodies are situated at the vertices of an equilateral triangle in the rotating frame. A figure eight choreography for was found numerically by C. Moore in 1993 and generalized and proven by A. Chenciner and R. Montgomery in 2000. Since then, many other choreographies have been found for . == Analytic approaches ==

For every solution of the problem, not only applying an isometry or a time shift but also a reversal of time (unlike in the case of friction) gives a solution as well. In the physical literature about the -body problem (), sometimes reference is made to "the impossibility of solving the -body problem" (via employing the above approach). However, care must be taken when discussing the 'impossibility' of a solution, as this refers only to the method of first integrals (compare the theorems by Abel and Galois about the impossibility of solving algebraic equations of degree five or higher by means of formulas only involving roots). Power series solution One way of solving the classical -body problem is "the -body problem by Taylor series". We start by defining the system of differential equations: \frac{d^2\mathbf{x}_i(t)}{dt^2}=G \sum_{k=1 \atop k\neq i}^n \frac{m_k \left(\mathbf{x}_k(t)-\mathbf{x}_i(t)\right)}{\left|\mathbf{x}_k(t)-\mathbf{x}_i(t)\right|^{3}}, As and are given as initial conditions, every is known. Differentiating results in which at which is also known, and the Taylor series is constructed iteratively. A generalized Sundman global solution In order to generalize Sundman's result for the case (or and ) one has to face two obstacles: • As has been shown by Siegel, collisions which involve more than two bodies cannot be regularized analytically, hence Sundman's regularization cannot be generalized. • The structure of singularities is more complicated in this case: other types of singularities may occur (see below). Lastly, Sundman's result was generalized to the case of bodies by Qiudong Wang in the 1990s. Since the structure of singularities is more complicated, Wang had to leave out completely the questions of singularities. The central point of his approach is to transform, in an appropriate manner, the equations to a new system, such that the interval of existence for the solutions of this new system is . Singularities of the -body problem There can be two types of singularities of the -body problem: • collisions of two or more bodies, but for which (the bodies' positions) remains finite. (In this mathematical sense, a "collision" means that two pointlike bodies have identical positions in space.) • singularities in which a collision does not occur, but does not remain finite. In this scenario, bodies diverge to infinity in a finite time, while at the same time tending towards zero separation (an imaginary collision occurs "at infinity"). The latter ones are called Painlevé's conjecture (no-collisions singularities). Their existence has been conjectured for by Painlevé (see Painlevé conjecture). Examples of this behavior for have been constructed by Xia and a heuristic model for by Gerver. Donald G. Saari has shown that for 4 or fewer bodies, the set of initial data giving rise to singularities has measure zero. == Simulation ==

While there are analytic solutions available for the classical (i.e. nonrelativistic) two-body problem and for selected configurations with , in general -body problems must be solved or simulated using numerical methods. which means that even small errors in integration may grow exponentially in time. Third, a simulation may be over large stretches of model time (e.g. millions of years) and numerical errors accumulate as integration time increases. There are a number of techniques to reduce errors in numerical integration. == Other -body problems ==

Most work done on the -body problem has been on the gravitational problem. But there exist other systems for which -body mathematics and simulation techniques have proven useful. In large scale electrostatics problems, such as the simulation of proteins and cellular assemblies in structural biology, the Coulomb potential has the same form as the gravitational potential, except that charges may be positive or negative, leading to repulsive as well as attractive forces. Fast Coulomb solvers are the electrostatic counterpart to fast multipole method simulators. These are often used with periodic boundary conditions on the region simulated and Ewald summation techniques are used to speed up computations. In statistics and machine learning, some models have loss functions of a form similar to that of the gravitational potential: a sum of kernel functions over all pairs of objects, where the kernel function depends on the distance between the objects in parameter space. Example problems that fit into this form include all-nearest-neighbors in manifold learning, kernel density estimation, and kernel machines. Alternative optimizations to reduce the time complexity to have been developed, such as dual tree algorithms, that have applicability to the gravitational -body problem as well. A technique in Computational fluid dynamics called Vortex Methods sees the vorticity in a fluid domain discretized onto particles which are then advected with the velocity at their centers. Because the fluid velocity and vorticity are related via a Poisson's equation, the velocity can be solved in the same manner as gravitation and electrostatics: as an -body summation over all vorticity-containing particles. The summation uses the Biot-Savart law, with vorticity taking the place of electrical current. In the context of particle-laden turbulent multiphase flows, determining an overall disturbance field generated by all particles is an -body problem. If the particles translating within the flow are much smaller than the flow's Kolmogorov scale, their linear Stokes disturbance fields can be superposed, yielding a system of 3 equations for 3 components of disturbance velocities at the location of particles. == See also ==