Central Force Problems and Kepler Orbits

Understanding the motion of planets, satellites, and charged particles under an inverse-square law force is a cornerstone of classical mechanics. This framework reduces seemingly complex two-body dynamics to a single-particle problem, providing profound insight into orbital shapes, stability, and scattering phenomena that govern systems from the solar system to atomic nuclei.

Reducing the Two-Body Problem

Many fundamental physical problems, such as a planet orbiting a star or an electron interacting with a nucleus, involve two bodies interacting via a central force—a force whose magnitude depends only on the distance between the bodies and points along the line connecting them. Solving the equations of motion for two coupled bodies is challenging, but a powerful simplification exists. We separate the motion into two parts: the motion of the system's center of mass (which moves at constant velocity if no external forces act) and the relative motion of one body with respect to the other.

The key step is introducing the reduced mass, defined as $\mu =\frac{m_1{m}_2}{m_1+m_2}$. The dynamics of the relative separation vector $\vec{r}={\vec{r}}_1-{\vec{r}}_2$ are then governed by a single equation: $\mu \ddot{\vec{r}}=\vec{F}(r)$, where $\vec{F}(r)$ is the central force law (e.g., gravity). This transformation reduces the complex two-body central force problem to an equivalent one-body problem of a particle of mass $\mu$ moving in a fixed central force field centered at the origin.

The Effective Potential and Energy Analysis

For any central force, angular momentum $\vec{L}$ is conserved. This confines the motion to a plane. Using polar coordinates $(r,\varphi )$ in this plane, we can express the total energy $E$ of the equivalent one-body system as:

$$E=\frac{1}{2}\mu {\dot{r}}^2+\underset{\text{Effective Potential }{U}_{\text{eff}}(r)}{\underbrace{\frac{L^2}{2\mu r^2}+V(r)}}.$$

The term $\frac{L^2}{2\mu r^2}$ is the centrifugal barrier arising from angular momentum conservation. The effective potential $U_{\text{eff}}(r)$ is the sum of the actual potential $V(r)$ and this centrifugal term. Analyzing motion is now simplified: the radial kinetic energy is $\frac{1}{2}\mu {\dot{r}}^2=E-U_{\text{eff}}(r)$. The particle can only move in regions where $E\ge U_{\text{eff}}(r)$. Turning points occur where $E=U_{\text{eff}}(r)$. This formalism allows us to classify orbits without solving the full equations of motion.

Classifying Orbit Types

The shape and stability of an orbit are determined by the interplay between total energy $E$ and the shape of $U_{\text{eff}}(r)$. For an attractive inverse-square force like gravity, $V(r)=-\frac{k}{r}$, where $k=G{m}_1{m}_2>0$.

• Bound Orbits (E < 0): The particle's energy is less than the maximum of the centrifugal barrier. It oscillates between a minimum ($r_{\text{min}}$) and maximum ($r_{\text{max}}$) distance—this is a closed elliptical orbit with the force center at one focus. A circular orbit is the special case where $E$ equals the minimum of $U_{\text{eff}}(r)$ and $r$ is constant.
• Unbound Orbits (E ≥ 0): If $E\ge 0$, the particle comes in from infinity, reaches a single turning point (periapsis), and escapes back to infinity. This is a hyperbolic orbit (or parabolic for $E=0$), typical of scattering or fly-by events.

Deriving Kepler's Laws from Newtonian Gravity

Applying the one-body reduction and effective potential to the specific case of Newton's law of universal gravitation yields Kepler's Laws as mathematical consequences.

1. Law of Ellipses: For bound systems ($E<0$) under an inverse-square attractive force, the direct solution of the orbit equation shows the path is a conic section. Closed, periodic orbits are ellipses with the central body (e.g., the Sun) at one focus.
2. Law of Equal Areas: This is a direct statement of conservation of angular momentum. The area swept out per unit time, $\frac{dA}{dt}=\frac{L}{2\mu }$, is constant.
3. Law of Periods: The square of the orbital period $T$ is proportional to the cube of the semi-major axis $a$ of the ellipse: $T^2=\frac{4{\pi }^2}{G(m_1+m_2)}{a}^3$. This derives from integrating the area of the ellipse and using the constant areal velocity. Note the dependence on the sum of the masses, a correction Kepler could not make.

Scattering for Inverse-Square Forces

The concept of scattering cross-section quantifies the likelihood of deflection when a particle on an unbound ($E>0$) hyperbolic orbit interacts with a central force. The key parameter is the impact parameter $b$, the initial perpendicular distance between the particle's path and a parallel line through the force center. For a repulsive inverse-square force (like Coulomb repulsion), the scattering angle $\Theta$ is given by the Rutherford formula:

$$\Theta =2\arctan \left(\frac{k^{\prime }}{2Eb}\right),$$

where $k^{\prime }$ is the force constant (e.g., $k^{\prime }=\frac{q_1{q}_2}{4\pi {ϵ}_0}$ for electrostatics). The differential scattering cross-section $\frac{d\sigma }{d\Omega }$ describes the effective area for scattering into a given solid angle. The famous Rutherford scattering result is:

$$\frac{d\sigma }{d\Omega }={\left(\frac{k^{\prime }}{4E}\right)}^2{\csc }^4\left(\frac{\Theta }{2}\right).$$

This formula's ${\csc }^4(\Theta /2)$ dependence confirmed the nuclear model of the atom, as it arises only from an inverse-square repulsive force.

Common Pitfalls

• Ignoring the reduced mass: Treating the central body as infinitely massive ($\mu \approx m_{\text{small}}$) is an excellent approximation for a planet around the Sun, but it fails for systems of comparable mass like binary stars. Always reduce the problem first.
• Misinterpreting the effective potential minimum: The minimum of $U_{\text{eff}}(r)$ corresponds to a stable circular orbit. Many students mistakenly believe all bound orbits (ellipses) occur at this minimum. In elliptical orbits, $E$ is slightly greater than the minimum, causing radial oscillation.
• Confusing energy and force: A particle can be in a bound orbit even if the actual force $F(r)$ is attractive everywhere (like gravity). Binding is determined by total energy $E$ relative to zero (the value at infinite separation), not the sign of the force.
• Applying Kepler's Third Law incorrectly: The constant of proportionality in $T^2\propto a^3$ depends on the total mass of the two-body system. Using a constant from our solar system (e.g., Earth-Sun) for a different system (e.g., a moon-planet) will yield an incorrect period.

Summary

• The two-body central force problem is reducible to an equivalent one-body problem using center-of-mass motion and the reduced mass $\mu$.
• Conservation of angular momentum allows the definition of an effective potential $U_{\text{eff}}(r)$, which combines the real potential with a centrifugal barrier. Plotting $U_{\text{eff}}$ against $r$ classifies orbits as bound (elliptical) or unbound (parabolic/hyperbolic) based on total energy.
• For Newtonian gravity, this analysis rigorously derives Kepler's Laws: elliptical orbits, constant areal velocity (angular momentum conservation), and the period-law relation $T^2\propto a^3$.
• Unbound orbits describe scattering. The relationship between impact parameter $b$ and scattering angle $\Theta$ leads to the concept of a cross-section, with the famous Rutherford formula confirming the inverse-square law nature of electrostatic forces.