Two-body problem

Last updated September 14, 2024

Left: Two bodies of similar mass orbiting a common barycenter external to both bodies, with elliptic orbits. This model is typical of binary stars.
Right: Two bodies with a "slight" difference in mass orbiting a common barycenter. Their sizes and this type of orbit are similar to the Pluto–Charon system (in which the barycenter is external to both bodies), as well as the Earth–Moon system (in which the barycenter is internal to the larger body).

In classical mechanics, the two-body problem is to calculate and predict the motion of two massive bodies that are orbiting each other in space. The problem assumes that the two bodies are point particles that interact only with one another; the only force affecting each object arises from the other one, and all other objects are ignored.

Results for prominent cases
Gravitation and other inverse-square examples
Inapplicability to atoms and subatomic particles
Reduction to two independent, one-body problems
Center of mass motion (1st one-body problem)
Displacement vector motion (2nd one-body problem)
Two-body motion is planar
Energy of the two-body system
Central forces
See also
References
Bibliography
External links

The most prominent example of the classical two-body problem is the gravitational case (see also Kepler problem), arising in astronomy for predicting the orbits (or escapes from orbit) of objects such as satellites, planets, and stars. A two-point-particle model of such a system nearly always describes its behavior well enough to provide useful insights and predictions.

A simpler "one body" model, the "central-force problem", treats one object as the immobile source of a force acting on the other. One then seeks to predict the motion of the single remaining mobile object. Such an approximation can give useful results when one object is much more massive than the other (as with a light planet orbiting a heavy star, where the star can be treated as essentially stationary).

However, the one-body approximation is usually unnecessary except as a stepping stone. For many forces, including gravitational ones, the general version of the two-body problem can be reduced to a pair of one-body problems, allowing it to be solved completely, and giving a solution simple enough to be used effectively.

By contrast, the three-body problem (and, more generally, the n-body problem for n ≥ 3) cannot be solved in terms of first integrals, except in special cases.

Results for prominent cases

Gravitation and other inverse-square examples

The two-body problem is interesting in astronomy because pairs of astronomical objects are often moving rapidly in arbitrary directions (so their motions become interesting), widely separated from one another (so they will not collide) and even more widely separated from other objects (so outside influences will be small enough to be ignored safely).

Under the force of gravity, each member of a pair of such objects will orbit their mutual center of mass in an elliptical pattern, unless they are moving fast enough to escape one another entirely, in which case their paths will diverge along other planar conic sections. If one object is very much heavier than the other, it will move far less than the other with reference to the shared center of mass. The mutual center of mass may even be inside the larger object.

For the derivation of the solutions to the problem, see Classical central-force problem or Kepler problem.

In principle, the same solutions apply to macroscopic problems involving objects interacting not only through gravity, but through any other attractive scalar force field obeying an inverse-square law, with electrostatic attraction being the obvious physical example. In practice, such problems rarely arise. Except perhaps in experimental apparatus or other specialized equipment, we rarely encounter electrostatically interacting objects which are moving fast enough, and in such a direction, as to avoid colliding, and/or which are isolated enough from their surroundings.

The dynamical system of a two-body system under the influence of torque turns out to be a Sturm-Liouville equation.^[1]

Inapplicability to atoms and subatomic particles

Although the two-body model treats the objects as point particles, classical mechanics only apply to systems of macroscopic scale. Most behavior of subatomic particles cannot be predicted under the classical assumptions underlying this article or using the mathematics here.

Electrons in an atom are sometimes described as "orbiting" its nucleus, following an early conjecture of Niels Bohr (this is the source of the term "orbital"). However, electrons don't actually orbit nuclei in any meaningful sense, and quantum mechanics are necessary for any useful understanding of the electron's real behavior. Solving the classical two-body problem for an electron orbiting an atomic nucleus is misleading and does not produce many useful insights.

Reduction to two independent, one-body problems

The complete two-body problem can be solved by re-formulating it as two one-body problems: a trivial one and one that involves solving for the motion of one particle in an external potential. Since many one-body problems can be solved exactly, the corresponding two-body problem can also be solved.

Let $x 1$ and $x 2$ be the vector positions of the two bodies, and m₁ and m₂ be their masses. The goal is to determine the trajectories $x 1 (t)$ and $x 2 (t)$ for all times t, given the initial positions $x 1 (t = 0)$ and $x 2 (t = 0)$ and the initial velocities $v 1 (t = 0)$ and $v 2 (t = 0)$ .

When applied to the two masses, Newton's second law states that

\mathbf {F} _{12}(\mathbf {x} _{1},\mathbf {x} _{2})=m_{1}{\ddot {\mathbf {x} }}_{1}

(Equation 1)

\mathbf {F} _{21}(\mathbf {x} _{1},\mathbf {x} _{2})=m_{2}{\ddot {\mathbf {x} }}_{2}

(Equation 2)

where F₁₂ is the force on mass 1 due to its interactions with mass 2, and F₂₁ is the force on mass 2 due to its interactions with mass 1. The two dots on top of the x position vectors denote their second derivative with respect to time, or their acceleration vectors.

Adding and subtracting these two equations decouples them into two one-body problems, which can be solved independently. Adding equations (1) and ( 2 ) results in an equation describing the center of mass (barycenter) motion. By contrast, subtracting equation (2) from equation (1) results in an equation that describes how the vector $r = x 1 - x 2$ between the masses changes with time. The solutions of these independent one-body problems can be combined to obtain the solutions for the trajectories $x 1 (t)$ and $x 2 (t)$ .

Center of mass motion (1st one-body problem)

Let $\mathbf {R}$ be the position of the center of mass (barycenter) of the system. Addition of the force equations (1) and (2) yields $m_{1}{\ddot {\mathbf {x} }}_{1}+m_{2}{\ddot {\mathbf {x} }}_{2}=(m_{1}+m_{2}){\ddot {\mathbf {R} }}=\mathbf {F} _{12}+\mathbf {F} _{21}=0$ where we have used Newton's third law $F 12 = - F 21$ and where ${\ddot {\mathbf {R} }}\equiv {\frac {m_{1}{\ddot {\mathbf {x} }}_{1}+m_{2}{\ddot {\mathbf {x} }}_{2}}{m_{1}+m_{2}}}.$

The resulting equation: ${\ddot {\mathbf {R} }}=0$ shows that the velocity $\mathbf {v} ={\frac {dR}{dt}}$ of the center of mass is constant, from which follows that the total momentum $m 1 v 1 + m 2 v 2$ is also constant (conservation of momentum). Hence, the position $R (t)$ of the center of mass can be determined at all times from the initial positions and velocities.

Displacement vector motion (2nd one-body problem)

Dividing both force equations by the respective masses, subtracting the second equation from the first, and rearranging gives the equation ${\ddot {\mathbf {r} }}={\ddot {\mathbf {x} }}_{1}-{\ddot {\mathbf {x} }}_{2}=\left({\frac {\mathbf {F} _{12}}{m_{1}}}-{\frac {\mathbf {F} _{21}}{m_{2}}}\right)=\left({\frac {1}{m_{1}}}+{\frac {1}{m_{2}}}\right)\mathbf {F} _{12}$ where we have again used Newton's third law $F 12 = - F 21$ and where $r$ is the displacement vector from mass 2 to mass 1, as defined above.

The force between the two objects, which originates in the two objects, should only be a function of their separation $r$ and not of their absolute positions $x 1$ and $x 2$ ; otherwise, there would not be translational symmetry, and the laws of physics would have to change from place to place. The subtracted equation can therefore be written: $\mu {\ddot {\mathbf {r} }}=\mathbf {F} _{12}(\mathbf {x} _{1},\mathbf {x} _{2})=\mathbf {F} (\mathbf {r} )$ where $\mu$ is the reduced mass $\mu ={\frac {1}{{\frac {1}{m_{1}}}+{\frac {1}{m_{2}}}}}={\frac {m_{1}m_{2}}{m_{1}+m_{2}}}.$

Solving the equation for $r (t)$ is the key to the two-body problem. The solution depends on the specific force between the bodies, which is defined by $\mathbf {F} (\mathbf {r} )$ . For the case where $\mathbf {F} (\mathbf {r} )$ follows an inverse-square law, see the Kepler problem.

Once $R (t)$ and $r (t)$ have been determined, the original trajectories may be obtained $\mathbf {x} _{1}(t)=\mathbf {R} (t)+{\frac {m_{2}}{m_{1}+m_{2}}}\mathbf {r} (t)$ $\mathbf {x} _{2}(t)=\mathbf {R} (t)-{\frac {m_{1}}{m_{1}+m_{2}}}\mathbf {r} (t)$ as may be verified by substituting the definitions of R and r into the right-hand sides of these two equations.

Two-body motion is planar

The motion of two bodies with respect to each other always lies in a plane (in the center of mass frame).

Proof: Defining the linear momentum $p$ and the angular momentum $L$ of the system, with respect to the center of mass, by the equations $\mathbf {L} =\mathbf {r} \times \mathbf {p} =\mathbf {r} \times \mu {\frac {d\mathbf {r} }{dt}},$

where $μ$ is the reduced mass and $r$ is the relative position $r 2 - r 1$ (with these written taking the center of mass as the origin, and thus both parallel to $r$ ) the rate of change of the angular momentum $L$ equals the net torque $N$ $\mathbf {N} ={\frac {d\mathbf {L} }{dt}}={\dot {\mathbf {r} }}\times \mu {\dot {\mathbf {r} }}+\mathbf {r} \times \mu {\ddot {\mathbf {r} }}\ ,$ and using the property of the vector cross product that $v \times w = 0$ for any vectors $v$ and $w$ pointing in the same direction,

$\mathbf {N} \ =\ {\frac {d\mathbf {L} }{dt}}=\mathbf {r} \times \mathbf {F} \ ,$ with $F = μ d 2 r / dt 2$ .

Introducing the assumption (true of most physical forces, as they obey Newton's strong third law of motion) that the force between two particles acts along the line between their positions, it follows that $r \times F = 0$ and the angular momentum vector $L$ is constant (conserved). Therefore, the displacement vector $r$ and its velocity $v$ are always in the plane perpendicular to the constant vector $L$ .

Energy of the two-body system

If the force $F (r)$ is conservative then the system has a potential energy $U (r)$ , so the total energy can be written as $E_{\text{tot}}={\frac {1}{2}}m_{1}{\dot {\mathbf {x} }}_{1}^{2}+{\frac {1}{2}}m_{2}{\dot {\mathbf {x} }}_{2}^{2}+U(\mathbf {r} )={\frac {1}{2}}(m_{1}+m_{2}){\dot {\mathbf {R} }}^{2}+{1 \over 2}\mu {\dot {\mathbf {r} }}^{2}+U(\mathbf {r} )$

In the center of mass frame the kinetic energy is the lowest and the total energy becomes $E={\frac {1}{2}}\mu {\dot {\mathbf {r} }}^{2}+U(\mathbf {r} )$ The coordinates $x 1$ and $x 2$ can be expressed as $\mathbf {x} _{1}={\frac {\mu }{m_{1}}}\mathbf {r}$ $\mathbf {x} _{2}=-{\frac {\mu }{m_{2}}}\mathbf {r}$ and in a similar way the energy E is related to the energies $E 1$ and $E 2$ that separately contain the kinetic energy of each body: ${\begin{aligned}E_{1}&={\frac {\mu }{m_{1}}}E={\frac {1}{2}}m_{1}{\dot {\mathbf {x} }}_{1}^{2}+{\frac {\mu }{m_{1}}}U(\mathbf {r} )\\[4pt]E_{2}&={\frac {\mu }{m_{2}}}E={\frac {1}{2}}m_{2}{\dot {\mathbf {x} }}_{2}^{2}+{\frac {\mu }{m_{2}}}U(\mathbf {r} )\\[4pt]E_{\text{tot}}&=E_{1}+E_{2}\end{aligned}}$

Central forces

For many physical problems, the force $F (r)$ is a central force, i.e., it is of the form $\mathbf {F} (\mathbf {r} )=F(r){\hat {\mathbf {r} }}$ where $r = | r |$ and $r̂ = r / r$ is the corresponding unit vector. We now have: $\mu {\ddot {\mathbf {r} }}={F}(r){\hat {\mathbf {r} }}\ ,$ where F(r) is negative in the case of an attractive force.

Related Research Articles

In astronomy, Kepler's laws of planetary motion, published by Johannes Kepler absent the third law in 1609 and fully in 1619, describe the orbits of planets around the Sun. The laws replaced the circular orbits and epicycles in the heliocentric theory of Nicolaus Copernicus with elliptical orbits and explained how planetary velocities vary. The three laws state that:

The orbit of a planet is an ellipse with the Sun at one of the two foci.
A line segment joining a planet and the Sun sweeps out equal areas during equal intervals of time.
The square of a planet's orbital period is proportional to the cube of the length of the semi-major axis of its orbit.

In celestial mechanics, an orbit is the curved trajectory of an object such as the trajectory of a planet around a star, or of a natural satellite around a planet, or of an artificial satellite around an object or position in space such as a planet, moon, asteroid, or Lagrange point. Normally, orbit refers to a regularly repeating trajectory, although it may also refer to a non-repeating trajectory. To a close approximation, planets and satellites follow elliptic orbits, with the center of mass being orbited at a focal point of the ellipse, as described by Kepler's laws of planetary motion.

In physics, reduced mass is a measure of the effective inertial mass of a system with two or more particles when the particles are interacting with each other. Reduced mass allows the two-body problem to be solved as if it were a one-body problem. Note, however, that the mass determining the gravitational force is not reduced. In the computation, one mass can be replaced with the reduced mass, if this is compensated by replacing the other mass with the sum of both masses. The reduced mass is frequently denoted by $(mu), although the standard gravitational parameter is also denoted by . It has the dimensions of mass, and SI unit kg.$

In science, work is the energy transferred to or from an object via the application of force along a displacement. In its simplest form, for a constant force aligned with the direction of motion, the work equals the product of the force strength and the distance traveled. A force is said to do positive work if it has a component in the direction of the displacement of the point of application. A force does negative work if it has a component opposite to the direction of the displacement at the point of application of the force.

In special relativity, four-momentum (also called momentum–energy or momenergy) is the generalization of the classical three-dimensional momentum to four-dimensional spacetime. Momentum is a vector in three dimensions; similarly four-momentum is a four-vector in spacetime. The contravariant four-momentum of a particle with relativistic energy $E$ and three-momentum $p = (p x, p y, p z) = γm v$ , where $v$ is the particle's three-velocity and $γ$ the Lorentz factor, is

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed by prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

In physics, circular motion is a movement of an object along the circumference of a circle or rotation along a circular arc. It can be uniform, with a constant rate of rotation and constant tangential speed, or non-uniform with a changing rate of rotation. The rotation around a fixed axis of a three-dimensional body involves the circular motion of its parts. The equations of motion describe the movement of the center of mass of a body, which remains at a constant distance from the axis of rotation. In circular motion, the distance between the body and a fixed point on its surface remains the same, i.e., the body is assumed rigid.

In physics and astronomy, the Reissner–Nordström metric is a static solution to the Einstein–Maxwell field equations, which corresponds to the gravitational field of a charged, non-rotating, spherically symmetric body of mass M. The analogous solution for a charged, rotating body is given by the Kerr–Newman metric.

In astrodynamics or celestial mechanics, an elliptic orbit or elliptical orbit is a Kepler orbit with an eccentricity of less than 1; this includes the special case of a circular orbit, with eccentricity equal to 0. In a stricter sense, it is a Kepler orbit with the eccentricity greater than 0 and less than 1. In a wider sense, it is a Kepler orbit with negative energy. This includes the radial elliptic orbit, with eccentricity equal to 1.

A circular orbit is an orbit with a fixed distance around the barycenter; that is, in the shape of a circle. In this case, not only the distance, but also the speed, angular speed, potential and kinetic energy are constant. There is no periapsis or apoapsis. This orbit has no radial version.

In celestial mechanics, the specific relative angular momentum of a body is the angular momentum of that body divided by its mass. In the case of two orbiting bodies it is the vector product of their relative position and relative linear momentum, divided by the mass of the body in question.

In celestial mechanics, a Kepler orbit is the motion of one body relative to another, as an ellipse, parabola, or hyperbola, which forms a two-dimensional orbital plane in three-dimensional space. A Kepler orbit can also form a straight line. It considers only the point-like gravitational attraction of two bodies, neglecting perturbations due to gravitational interactions with other objects, atmospheric drag, solar radiation pressure, a non-spherical central body, and so on. It is thus said to be a solution of a special case of the two-body problem, known as the Kepler problem. As a theory in classical mechanics, it also does not take into account the effects of general relativity. Keplerian orbits can be parametrized into six orbital elements in various ways.

In orbital mechanics, the universal variable formulation is a method used to solve the two-body Kepler problem. It is a generalized form of Kepler's Equation, extending it to apply not only to elliptic orbits, but also parabolic and hyperbolic orbits common for spacecraft departing from a planetary orbit. It is also applicable to ejection of small bodies in Solar System from the vicinity of massive planets, during which processes the approximating two-body orbits can have widely varying eccentricities, almost always $e \geq 1$ .

In classical mechanics, the Udwadia–Kalaba formulation is a method for deriving the equations of motion of a constrained mechanical system. The method was first described by Anatolii Fedorovich Vereshchagin for the particular case of robotic arms, and later generalized to all mechanical systems by Firdaus E. Udwadia and Robert E. Kalaba in 1992. The approach is based on Gauss's principle of least constraint. The Udwadia–Kalaba method applies to both holonomic constraints and nonholonomic constraints, as long as they are linear with respect to the accelerations. The method generalizes to constraint forces that do not obey D'Alembert's principle.

<span class="mw-page-title-main">Lagrangian mechanics</span> Formulation of classical mechanics

In physics, Lagrangian mechanics is a formulation of classical mechanics founded on the stationary-action principle. It was introduced by the Italian-French mathematician and astronomer Joseph-Louis Lagrange in his presentation to the Turin Academy of Science in 1760 culminating in his 1788 grand opus, Mécanique analytique.

In classical mechanics, the central-force problem is to determine the motion of a particle in a single central potential field. A central force is a force that points from the particle directly towards a fixed point in space, the center, and whose magnitude only depends on the distance of the object to the center. In a few important cases, the problem can be solved analytically, i.e., in terms of well-studied functions such as trigonometric functions.

Orbit modeling is the process of creating mathematical models to simulate motion of a massive body as it moves in orbit around another massive body due to gravity. Other forces such as gravitational attraction from tertiary bodies, air resistance, solar pressure, or thrust from a propulsion system are typically modeled as secondary effects. Directly modeling an orbit can push the limits of machine precision due to the need to model small perturbations to very large orbits. Because of this, perturbation methods are often used to model the orbit in order to achieve better accuracy.

In theoretical physics, relativistic Lagrangian mechanics is Lagrangian mechanics applied in the context of special relativity and general relativity.

The Clohessy–Wiltshire equations describe a simplified model of orbital relative motion, in which the target is in a circular orbit, and the chaser spacecraft is in an elliptical or circular orbit. This model gives a first-order approximation of the chaser's motion in a target-centered coordinate system. It is used to plan the rendezvous of the chaser with the target.

Accelerations in special relativity (SR) follow, as in Newtonian Mechanics, by differentiation of velocity with respect to time. Because of the Lorentz transformation and time dilation, the concepts of time and distance become more complex, which also leads to more complex definitions of "acceleration". SR as the theory of flat Minkowski spacetime remains valid in the presence of accelerations, because general relativity (GR) is only required when there is curvature of spacetime caused by the energy–momentum tensor. However, since the amount of spacetime curvature is not particularly high on Earth or its vicinity, SR remains valid for most practical purposes, such as experiments in particle accelerators.

References

↑ Luo, Siwei (22 June 2020). "The Sturm-Liouville problem of two-body system". Journal of Physics Communications. 4 (6): 061001. Bibcode:2020JPhCo...4f1001L. doi: 10.1088/2399-6528/ab9c30 .
↑ David Betounes (2001). Differential Equations . Springer. p. 58; Figure 2.15. ISBN 0-387-95140-7.

Bibliography

Landau LD; Lifshitz EM (1976). Mechanics (3rd. ed.). New York: Pergamon Press. ISBN 0-08-029141-4.
Goldstein H (1980). Classical Mechanics (2nd. ed.). New York: Addison-Wesley. ISBN 0-201-02918-9.

External links

Two-body problem at Eric Weisstein's World of Physics

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Luo, Siwei (22 June 2020). "The Sturm-Liouville problem of two-body system". Journal of Physics Communications. 4 (6): 061001. Bibcode:2020JPhCo...4f1001L. doi: 10.1088/2399-6528/ab9c30 .

[Betounes-2] David Betounes (2001). Differential Equations . Springer. p. 58; Figure 2.15. ISBN 0-387-95140-7.

[1]