The Effective Potential and the Orbit Equation

The Centrifugal Barrier

Imagine a particle moving in a plane under a central force. Because the force always points towards the origin, it exerts no torque about the origin. Consequently, the angular momentum \(L\) is conserved. This conservation implies that the particle cannot simply fall straight into the center unless it has zero angular momentum. Instead, the conservation of \(L\) creates a "centrifugal barrier" that prevents the particle from reaching \(r=0\) if \(L \neq 0\). We can interpret this barrier as an additional potential energy term that grows rapidly as \(r\) decreases. This allows us to treat the radial motion as if the particle were moving in one dimension under the influence of a modified potential, which we call the effective potential. The total energy is then the sum of the radial kinetic energy and this effective potential energy.

Effective Potential and Radial Motion

The total kinetic energy in polar coordinates is \(T = \frac{1}{2}\mu \dot{r}^2 + \frac{1}{2}\mu r^2 \dot{\theta}^2\). Using the conservation of angular momentum \(L = \mu r^2 \dot{\theta}\), we can substitute \(\dot{\theta} = \frac{L}{\mu r^2}\) into the angular part of the kinetic energy. This yields \(\frac{1}{2}\mu r^2 \left(\frac{L}{\mu r^2}\right)^2 = \frac{L^2}{2\mu r^2}\). This term acts exactly like a potential energy. Thus, the total energy \(E\) can be written as:

where \(U_{\text{eff}}(r) = U(r) + \frac{L^2}{2\mu r^2}\). The radial motion is now governed by the single equation \(\dot{r}^2 = \frac{2}{\mu}(E - U_{\text{eff}}(r))\). The particle can only exist where \(E \geq U_{\text{eff}}(r)\). The points where \(E = U_{\text{eff}}(r)\) are the turning points, or apsides, where the radial velocity \(\dot{r}\) vanishes. If \(E\) is less than the minimum of \(U_{\text{eff}}\), the motion is bounded (an orbit). If \(E\) is greater than the maximum of \(U_{\text{eff}}\) (for attractive forces), the motion is unbounded (scattering).

The Binet Orbit Equation

While the effective potential tells us about the radial extent of the motion, it does not directly give the shape of the path. To find the orbit \(r(\theta)\), we use a change of variable. Let \(u = 1/r\). We need to express the acceleration in terms of \(u\) and \(\theta\). The radial acceleration in polar coordinates is \(\ddot{r} - r\dot{\theta}^2\). Using the chain rule and the conservation of \(L\), one can derive that the radial force equation \(\mu(\ddot{r} - r\dot{\theta}^2) = F(r)\) transforms into a second-order differential equation for \(u(\theta)\). This is the Binet equation:

This equation is powerful because it relates the geometry of the orbit directly to the force law. For a specific force law, such as the inverse-square law \(F(r) = -k/r^2\), the right-hand side becomes a constant, leading to a simple harmonic oscillator equation in \(u\). The solution to this linear differential equation yields the conic sections (ellipses, parabolas, hyperbolas) that characterize Keplerian orbits.

Stability of Circular Orbits

A circular orbit occurs at a radius \(r_0\) where the effective potential has a local minimum, i.e., \(\frac{dU_{\text{eff}}}{dr}\big|_{r_0} = 0\). This condition implies that the centripetal force required for circular motion is exactly provided by the central force. The stability of this orbit depends on the curvature of the effective potential at \(r_0\). If \(\frac{d^2U_{\text{eff}}}{dr^2}\big|_{r_0} > 0\), the minimum is stable, and small perturbations produce small radial oscillations about \(r_0\) (for an inverse-square force the perturbed orbit closes into an ellipse; for other force laws it generally precesses). If the second derivative is negative, the orbit is unstable. For the inverse-square law, all circular orbits are stable. However, for other force laws, such as \(F(r) \propto -1/r^3\), the effective potential may not have a stable minimum, leading to unstable orbits.