Two-body problem

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 also through any conservative force obeying the inverse-square law, with electrostatic attraction being an 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. Inapplicability to subatomic systems Although we can describe the repulsion and attraction of particles in the atom using the same laws that apply to the leaves of an electroscope, we need quantum mechanics to predict how the particles will behave under those forces. Most behavior of subatomic systems 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. for two-body problem; Jacobi coordinates are \boldsymbol{R}=\frac {m_1}{M} \boldsymbol{x}_1 + \frac {m_2}{M} \boldsymbol{x}_2 and \boldsymbol{r} = \boldsymbol{x}_1 - \boldsymbol{x}_2 with M = m_1+m_2 \ . Let and be the vector positions of the two bodies, and m1 and m2 be their masses. The goal is to determine the trajectories and for all times t, given the initial positions and and the initial velocities and . When applied to the two masses, Newton's second law states that {{NumBlk||\mathbf{F}_{12}(\mathbf{x}_{1},\mathbf{x}_{2}) = m_{1} \ddot{\mathbf{x}}_{1} |Equation }} {{NumBlk||\mathbf{F}_{21}(\mathbf{x}_{1},\mathbf{x}_{2}) = m_{2} \ddot{\mathbf{x}}_{2} |Equation }} where F12 is the force on mass 1 due to its interactions with mass 2, and F21 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 () 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 between the masses changes with time. The solutions of these independent one-body problems can be combined to obtain the solutions for the trajectories and . 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 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 is also constant (conservation of momentum). Hence, the position 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 and where 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 and not of their absolute positions and ; 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 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 and 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 and the angular momentum 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 is the relative position (with these written taking the center of mass as the origin, and thus both parallel to ) the rate of change of the angular momentum equals the net torque \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 for any vectors and pointing in the same direction, \mathbf{N} \ = \ \frac{d\mathbf{L}}{dt} = \mathbf{r} \times \mathbf{F} \ , with . 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 and the angular momentum vector is constant (conserved). Therefore, the displacement vector and its velocity are always in the plane perpendicular to the constant vector . == Energy of the two-body system ==

If the force is conservative then the system has a potential energy , 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 and 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 and that separately contain the kinetic energy of each body: \begin{align} 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{align} == Central forces ==

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