physics2u
Depth
PU-201 / L11 · lecture

The Two-Body Problem and Radial Reduction

~6 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)

BEFORE YOU START

You should be comfortable with Lagrangian mechanics, the identification of cyclic coordinates, and the conservation of angular momentum in central force fields.

Why this matters

Complex orbital motion becomes tractable when reduced to a one-dimensional problem.

The two-body problem, such as a planet orbiting a star, appears to involve two interacting masses moving in three dimensions. However, by separating the motion of the centre of mass from the relative motion, we reduce the system to a single particle moving in a central potential. This reduction allows us to apply the powerful tools of one-dimensional mechanics, specifically the effective potential, to analyse stability, turning points, and the shape of orbits without solving the full coupled equations of motion.

From Two Bodies to One Particle

Imagine two masses, \( m_1 \) and \( m_2 \), interacting via a force that depends only on the distance between them. Newton’s third law tells us the forces are equal and opposite. If we look at the system from a frame moving with the centre of mass, the total momentum is zero. The motion of the two bodies is perfectly correlated: as one moves, the other moves in a specific way to keep the centre of mass fixed. Instead of tracking two separate positions, we can define a single vector \( \mathbf{r} = \mathbf{r}_2 - \mathbf{r}_1 \) that points from mass 1 to mass 2. The dynamics of this relative vector are mathematically identical to the dynamics of a single particle with a specific "reduced mass" \( \mu \) moving in a potential field created by the other mass. This is the essence of the reduction: we trade a complex two-body interaction for a simpler one-body problem in a central field.

The Reduced Mass

To perform this reduction rigorously, we introduce the reduced mass \( \mu \), defined by the relation \( \frac{1}{\mu} = \frac{1}{m_1} + \frac{1}{m_2} \). This quantity has units of mass and captures the inertial resistance of the relative coordinate. If one mass is much larger than the other, \( \mu \) approaches the smaller mass, recovering the familiar case of a test particle in a fixed potential. The kinetic energy of the two-body system can be decomposed into the kinetic energy of the centre of mass plus the kinetic energy of the relative motion, \( T = \frac{1}{2} M \dot{\mathbf{R}}^2 + \frac{1}{2} \mu \dot{\mathbf{r}}^2 \), where \( M = m_1 + m_2 \). Since the potential energy \( V \) depends only on the separation \( r = |\mathbf{r}| \), the Lagrangian for the relative motion is \( L = \frac{1}{2} \mu \dot{\mathbf{r}}^2 - V(r) \). This is formally identical to the Lagrangian of a single particle of mass \( \mu \) in a central potential.

Radial Reduction and the Effective Potential

Because the potential is central, the angular momentum \( \ell \) is conserved. In polar coordinates \( (r, \phi) \), the kinetic energy is \( T = \frac{1}{2} \mu (\dot{r}^2 + r^2 \dot{\phi}^2) \). The angular momentum is \( \ell = \mu r^2 \dot{\phi} \). We can eliminate \( \dot{\phi} \) from the energy equation by substituting \( \dot{\phi} = \ell / (\mu r^2) \). The total energy \( E \) becomes:

\[ E = \frac{1}{2} \mu \dot{r}^2 + \frac{\ell^2}{2 \mu r^2} + V(r) \]

This equation looks exactly like the energy of a one-dimensional particle with kinetic energy \( \frac{1}{2} \mu \dot{r}^2 \) moving in a potential \( V_{\mathrm{eff}}(r) \). The effective potential is defined as:

\[ V_{\mathrm{eff}}(r) = V(r) + \frac{\ell^2}{2 \mu r^2} \]

The term \( \frac{\ell^2}{2 \mu r^2} \) is the centrifugal barrier. It arises from the conservation of angular momentum and prevents the particle from reaching \( r=0 \) if \( \ell \neq 0 \). The radial motion is governed by the equation \( \mu \ddot{r} = -\frac{dV_{\mathrm{eff}}}{dr} \). This reduction is crucial because it allows us to analyse the radial dynamics using one-dimensional techniques, such as finding turning points where \( E = V_{\mathrm{eff}}(r) \).

Interpreting the Effective Potential

The shape of \( V_{\mathrm{eff}}(r) \) determines the allowed regions of motion. For an attractive inverse-square potential, \( V(r) = -\frac{k}{r} \), the effective potential is \( V_{\mathrm{eff}}(r) = -\frac{k}{r} + \frac{\ell^2}{2 \mu r^2} \). This function has a minimum at \( r_0 = \frac{\ell^2}{\mu k} \). If the total energy \( E \) is less than zero, the motion is bounded, corresponding to an elliptical orbit. The particle oscillates between a minimum radius \( r_{\min} \) and a maximum radius \( r_{\max} \), which are the roots of \( E = V_{\mathrm{eff}}(r) \). If \( E \) is greater than zero, the motion is unbounded, corresponding to a hyperbolic trajectory. The centrifugal term dominates at small \( r \), creating a repulsive barrier, while the attractive term dominates at large \( r \). This interplay dictates the stability of circular orbits and the existence of bound states.

The derivation

◆ DERIVATION EMBED · D-130 · not re-derived here

Effective Potential and Radial Reduction

For a single particle of mass m moving in a central potential V(r), the planar two-degree-of-freedom problem in (r, φ) reduces to a single one-dimensional equation for the radial coordinate, m r̈ = −dVeff/dr, in which the conserved angular momentum ℓ is folded into an effective potential Veff(r) = V(r) + ℓ2/(2m r2). The angular motion is eliminated, not ignored: it survives entirely inside the centrifugal term.

m r̈ = −dVeff/dr ,   Veff(r) = V(r) + ℓ2/(2m r2)
Open full derivation →
Tier 3 · watch out

A common misconception is that the centrifugal term \( \frac{\ell^2}{2 \mu r^2} \) is a real force acting on the particle. It is not. It is a mathematical artifact of using polar coordinates and eliminating the angular degree of freedom. In an inertial frame, the only real force is the central force \( \mathbf{F} = -\nabla V \). The centrifugal term appears in the radial equation of motion because the basis vectors \( \hat{\mathbf{r}} \) and \( \hat{\boldsymbol{\phi}} \) are rotating. Confusing this effective potential term with a physical force can lead to errors when analysing non-inertial frames or generalised coordinates. C

Worked example

Problem. Two masses, \( m_1 = 2.0 \times 10^{24} \ \mathrm{kg} \) and \( m_2 = 6.0 \times 10^{24} \ \mathrm{kg} \), interact gravitationally. Calculate the reduced mass \( \mu \). If the angular momentum of the relative motion is \( \ell = 1.0 \times 10^{40} \ \mathrm{kg \ m^2/s} \), find the radius \( r_0 \) of the circular orbit, assuming the gravitational constant \( G = 6.67 \times 10^{-11} \ \mathrm{N \ m^2/kg^2} \).

  1. Calculate the reduced mass: \( \mu = \frac{m_1 m_2}{m_1 + m_2} = \frac{(2.0 \times 10^{24})(6.0 \times 10^{24})}{2.0 \times 10^{24} + 6.0 \times 10^{24}} = \frac{12.0 \times 10^{48}}{8.0 \times 10^{24}} = 1.5 \times 10^{24} \ \mathrm{kg} \).
  2. Gravity gives \( V(r) = -k/r \) with \( k = G m_1 m_2 \) (the full product of the masses, not \( G\mu M \)). A circular orbit sits at the minimum of \( V_{\mathrm{eff}} = -\frac{k}{r} + \frac{\ell^2}{2\mu r^2} \): \( \frac{\mathrm{d}V_{\mathrm{eff}}}{\mathrm{d}r} = \frac{k}{r^2} - \frac{\ell^2}{\mu r^3} = 0 \), so \( r_0 = \frac{\ell^2}{\mu k} \).
  3. Substitute values: \( k = (6.67 \times 10^{-11})(2.0 \times 10^{24})(6.0 \times 10^{24}) = 8.004 \times 10^{38} \ \mathrm{J \ m} \).
  4. Calculate \( r_0 \): \( r_0 = \frac{(1.0 \times 10^{40})^2}{(1.5 \times 10^{24})(8.004 \times 10^{38})} = \frac{1.0 \times 10^{80}}{1.2006 \times 10^{63}} \approx 8.33 \times 10^{16} \ \mathrm{m} \).
\[ r_0 \approx 8.33 \times 10^{16} \ \mathrm{m} \]

Units check. The numerator is \( \mathrm{kg^2 m^4/s^2} \). The denominator is \( \mathrm{kg} \cdot \mathrm{kg m^3/s^2} = \mathrm{kg^2 m^3/s^2} \). The ratio is \( \mathrm{m} \), which is correct for a radius.

CHECK YOURSELF

Three quick tests

  1. What is the reduced mass of two equal masses \( m \)?
    Answer\( \mu = m/2 \).
  2. Why does the effective potential have a centrifugal term?
    AnswerIt arises from the conservation of angular momentum, which constrains the radial motion.
  3. For an attractive potential, what does a minimum in \( V_{\mathrm{eff}} \) represent?
    AnswerA stable circular orbit.
PROBLEMS

Problem set for PU-201 →

Each derivation in this unit ends with five graded problems with full solutions.

NOTES

Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.

Derivations cited: D-130 · Threads: force, energy, symmetry