Skip to content

Central Force Problems | Physics

6.1 Reduction to One Dimension

For a particle of mass mm in a central potential V(r)V(r)Using polar coordinates (r,ϕ)(r, \phi) in the Plane of motion:

L=12m(r˙2+r2ϕ˙2)V(r)L = \frac{1}{2}m(\dot{r}^2 + r^2\dot{\phi}^2) - V(r)

The angular momentum l=mr2ϕ˙l = mr^2\dot{\phi} is conserved. Substituting ϕ˙=l/(mr2)\dot{\phi} = l/(mr^2):

L=12mr˙2+l22mr2V(r)L = \frac{1}{2}m\dot{r}^2 + \frac{l^2}{2mr^2} - V(r)

The effective potential is Veff(r)=V(r)+l22mr2V_{\mathrm{eff}(r) = V(r) + \frac{l^2}{2mr^2}}Where the second term is The centrifugal barrier.

6.2 Effective Potential Analysis

Definition. The effective one-dimensional energy is:

E=12mr˙2+Veff(r)E = \frac{1}{2}m\dot{r}^2 + V_{\mathrm{eff}(r)}

Since EE and ll are conserved, the radial motion is completely determined by Veff(r)V_{\mathrm{eff}(r)}.

Circular orbits occur at radii r0r_0 where Veff(r0)=0V_{\mathrm{eff}'(r_0) = 0}:

V(r0)l2mr03=0V'(r_0) - \frac{l^2}{mr_0^3} = 0

The orbit is stable if Veff(r0)>0V_{\mathrm{eff}''(r_0) \gt 0} and unstable if Veff(r0)<0V_{\mathrm{eff}''(r_0) \lt 0}.

For the Kepler problem V(r)=k/rV(r) = -k/r:

Veff(r)=kr+l22mr2V_{\mathrm{eff}(r) = -\frac{k}{r} + \frac{l^2}{2mr^2}}

Veff(r)=kr2l2mr3=0    r0=l2mkV_{\mathrm{eff}'(r) = \frac{k}{r^2} - \frac{l^2}{mr^3} = 0 \implies r_0 = \frac{l^2}{mk}}

Veff(r0)=2kr03+3l2mr04=m3k2l4>0V_{\mathrm{eff}''(r_0) = -\frac{2k}{r_0^3} + \frac{3l^2}{mr_0^4} = \frac{m^3k^2}{l^4} \gt 0}

So the circular orbit is always stable for the Kepler problem.

{/* prettier-ignore */}

The effective potential Veff(r)=k/r+l2/(2r2)V_{\mathrm{eff}}(r) = -k/r + l^2/(2r^2) (blue) combines the attractive 1/r-1/r well with the centrifugal barrier 1/r2\propto 1/r^2. Adjust sliders k (force constant), l (angular momentum), and E (total energy) to explore bound and unbound orbits.

6.3 The Orbit Equation

Starting from conservation of energy and angular momentum, we derive the orbit equation. Let u=1/ru = 1/r and use d/dt=(l/mr2)d/dϕ=(lu2/m)d/dϕd/dt = (l/mr^2)\, d/d\phi = (lu^2/m)\, d/d\phi:

r˙=drdt=lmr2drdϕ=lmdudϕ\dot{r} = \frac{dr}{dt} = \frac{l}{mr^2}\frac{dr}{d\phi} = -\frac{l}{m}\frac{du}{d\phi}

Substituting into the energy equation:

E=l22m(dudϕ)2+l2u22m+V(1u)E = \frac{l^2}{2m}\left(\frac{du}{d\phi}\right)^2 + \frac{l^2 u^2}{2m} + V\left(\frac{1}{u}\right)

Differentiating with respect to ϕ\phi:

l2mdudϕd2udϕ2+l2umdudϕ1u2V(1u)dudϕ=0\frac{l^2}{m}\frac{du}{d\phi}\frac{d^2u}{d\phi^2} + \frac{l^2 u}{m}\frac{du}{d\phi} - \frac{1}{u^2}V'\left(\frac{1}{u}\right)\frac{du}{d\phi} = 0

Dividing by l2u/(m)l^2 u'/(m) (assuming u0u' \neq 0):

d2udϕ2+u=ml2u2V(1u)\frac{d^2u}{d\phi^2} + u = -\frac{m}{l^2 u^2}V'\left(\frac{1}{u}\right)

This is the Binet equation.

6.4 The Kepler Problem

For V(r)=k/rV(r) = -k/r (gravitational or Coulomb potential):

V(1u)=dVdr=kr2=ku2V'\left(\frac{1}{u}\right) = \frac{dV}{dr} = \frac{k}{r^2} = ku^2

The Binet equation becomes:

d2udϕ2+u=mkl2\frac{d^2u}{d\phi^2} + u = \frac{mk}{l^2}

This is an inhomogeneous ODE with solution:

u(ϕ)=mkl2(1+ecos(ϕϕ0))u(\phi) = \frac{mk}{l^2}(1 + e\cos(\phi - \phi_0))

Where ee is the eccentricity and ϕ0\phi_0 is a constant. Setting ϕ0=0\phi_0 = 0:

r(ϕ)=l2/(mk)1+ecosϕ=p1+ecosϕr(\phi) = \frac{l^2/(mk)}{1 + e\cos\phi} = \frac{p}{1 + e\cos\phi}

Where p=l2/(mk)p = l^2/(mk) is the semi-latus rectum.

The eccentricity is determined by the energy:

e=1+2El2mk2e = \sqrt{1 + \frac{2El^2}{mk^2}}

Classification of orbits:

EnergyEccentricityOrbit Type
E<0E \lt 0e<1e \lt 1Ellipse (bound)
E=0E = 0e=1e = 1Parabola (marginally bound)
E>0E \gt 0e>1e \gt 1Hyperbola (unbound)

Bertrand’s Theorem: The only central potentials that give closed orbits for all bound states are V(r)1/rV(r) \propto 1/r (Kepler/Coulomb) and V(r)r2V(r) \propto r^2 (harmonic oscillator).

Key results for Kepler orbits:

  • Orbits are conic sections (ellipses, parabolas, or hyperbolas).
  • The semi-major axis aa satisfies E=k/(2a)E = -k/(2a) for bound orbits.
  • The period is T=2πma3/kT = 2\pi\sqrt{ma^3/k} (Kepler’s third law).

6.5 Worked Example: Satellite Orbit

Problem. A satellite of mass mm orbits Earth (M=5.97×1024kgM = 5.97 \times 10^{24}\,\mathrm{kg}) in an elliptical orbit with perigee (closest approach) rp=7000kmr_p = 7000\,\mathrm{km} and apogee (farthest point) ra=42000kmr_a = 42000\,\mathrm{km}. Find the eccentricity, semi-major axis, and orbital period.

Solution

The semi-major axis is a=(rp+ra)/2=(7000+42000)/2=24500kma = (r_p + r_a)/2 = (7000 + 42000)/2 = 24500\,\mathrm{km}.

From the orbit equation, at perigee (ϕ=0\phi = 0): rp=p/(1+e)r_p = p/(1 + e) And at apogee (ϕ=π\phi = \pi): ra=p/(1e)r_a = p/(1 - e).

Therefore:

e=rarpra+rp=42000700042000+7000=35000490000.714e = \frac{r_a - r_p}{r_a + r_p} = \frac{42000 - 7000}{42000 + 7000} = \frac{35000}{49000} \approx 0.714

The energy is E=k/(2a)E = -k/(2a) where k=GMm=6.674×1011×5.97×1024×m=3.986×1014mm3/s2k = GMm = 6.674 \times 10^{-11} \times 5.97 \times 10^{24} \times m = 3.986 \times 10^{14} m\,\mathrm m^3/\mathrm s^2.

The period (independent of mass mm):

T=2πma3k=2πa3GM=2π(2.45×107)33.986×1014T = 2\pi\sqrt{\frac{ma^3}{k}} = 2\pi\sqrt{\frac{a^3}{GM}} = 2\pi\sqrt{\frac{(2.45 \times 10^7)^3}{3.986 \times 10^{14}}}

=2π1.471×10223.986×1014=2π3.691×1072π×607538170s10.6hours= 2\pi\sqrt{\frac{1.471 \times 10^{22}}{3.986 \times 10^{14}}} = 2\pi\sqrt{3.691 \times 10^7} \approx 2\pi \times 6075 \approx 38170\,\mathrm s \approx 10.6\,\mathrm{hours}

\blacksquare

6.6 Worked Example: Rutherford Scattering

Problem. A particle of mass mm and energy E>0E \gt 0 is scattered by a repulsive Coulomb potential V(r)=k/rV(r) = k/r (k>0k \gt 0). Find the scattering angle Θ\Theta as a function of the impact parameter bb.

Solution

The angular momentum is l=mvbl = mv_\infty b where v=2E/mv_\infty = \sqrt{2E/m}. The eccentricity is:

e=1+2El2mk2=1+(2Ebk/m)2e = \sqrt{1 + \frac{2El^2}{mk^2}} = \sqrt{1 + \left(\frac{2Eb}{k/m}\right)^2}

The orbit is a hyperbola. The asymptotic angles satisfy rr \to \inftyI.e., 1+ecosϕ=01 + e\cos\phi = 0Giving ϕ±=π±arccos(1/e)\phi_{\pm} = \pi \pm \arccos(1/e). The scattering angle is:

Θ=π2arccos(1/e)=2arcsin(1/e)\Theta = \pi - 2\arccos(1/e) = 2\arcsin(1/e)

Using sin(Θ/2)=1/e\sin(\Theta/2) = 1/e:

cotΘ2=2Ebk/m=mv2bk/m\cot\frac{\Theta}{2} = \frac{2Eb}{k/m} = \frac{mv_\infty^2 b}{k/m}

This is the Rutherford scattering formula relating the scattering angle to the impact parameter. \blacksquare

flowchart TD
    A[6_Central Force Problemsx] --> B[Key Concepts]
    A --> C[Core Principles]
    A --> D[Practical Applications]
    B --> E[Fundamental definitions]
    C --> F[Design patterns]
    D --> G[Real-world usage]

Intuition

Central force problems describe motion under forces that point toward or away from a fixed center. The key insight is that angular momentum is conserved, confining motion to a plane. The effective potential combines the real potential with the centrifugal barrier, creating an energy landscape that determines the orbit shape. For attractive inverse-square forces, the orbits are conic sections: circles, ellipses, parabolas, or hyperbolas. Rutherford scattering shows how alpha particles deflected by the Coulomb force reveal the atomic nucleus. The orbit equation transforms the problem from tracking position versus time to tracking the trajectory shape directly, eliminating time as a variable.

Common Mistakes

Mistake 1: Confusing the effective potential with the actual potential The effective potential Veff(r)=V(r)+l2/(2mr2)V_{\mathrm{eff}}(r) = V(r) + l^2/(2mr^2) includes the centrifugal barrier term l2/(2mr2)l^2/(2mr^2), which is not a real force but a consequence of angular momentum conservation. Students often forget this term and attempt to find circular orbits from V(r0)=0V'(r_0) = 0 instead of the correct Veff(r0)=0V_{\mathrm{eff}}'(r_0) = 0.

Mistake 2: Assuming all central force orbits are closed Bertrand’s theorem states that only the Kepler potential (V1/rV \propto 1/r) and the harmonic oscillator potential (Vr2V \propto r^2) produce closed orbits for all bound states. For other central potentials, orbits generally precess and are not closed. The orbit equation still applies, but the trajectory does not repeat.

Mistake 3: Mixing up the semi-latus rectum with the semi-major axis The semi-latus rectum p=l2/(mk)p = l^2/(mk) determines the shape of the orbit, while the semi-major axis aa determines the energy via E=k/(2a)E = -k/(2a). For a circular orbit p=ap = a, but for eccentric orbits p<ap < a. Using pp in place of aa in energy formulas gives incorrect results.

Cross-References

  • Lagrangian Mechanics: The Lagrangian formulation provides the effective potential analysis for central force problems.
  • Hamiltonian Mechanics: The Hamiltonian formalism gives the total energy and phase space representation of central force orbits.

Advanced Content

This section provides detailed coverage of advanced concepts, including full derivations, proofs, and extended examples.

Derivations and Proofs

Complete mathematical derivations and proofs are provided where appropriate. Each step is explained to ensure understanding of the underlying reasoning.

Extended Examples

Advanced examples demonstrate the application of concepts to complex problems. These examples go beyond standard exam questions to develop deeper understanding.

Research Connections

This material connects to current research and advanced applications in the field. Understanding these connections provides context for the study material.

Prerequisites

Ensure you have mastered the prerequisite material before attempting this advanced content.