Classical MechanicsBSc · Notes

1. Separating radial and angular motion

A central-force orbit lies in a fixed plane, but the particle still changes both radius and angle. Angular-momentum conservation lets us express the angular speed in terms of radius. The remaining problem can be studied as motion along a radial line in a modified potential.

Throughout this chapter, μ\mu is the particle or reduced mass, V(r)V(r) is a time-independent potential, and ℓ=μr2θ˙>0\ell=\mu r^2\dot\theta>0 is fixed angular momentum magnitude. The radial case ℓ=0\ell=0 is discussed separately.

2. Polar acceleration from the moving basis

Prerequisite. In the orbital plane, define r^=(cos⁡θ,sin⁡θ)\hat{\mathbf r}=(\cos\theta,\sin\theta) and θ^=(−sin⁡θ,cos⁡θ)\hat{\boldsymbol\theta}=(-\sin\theta,\cos\theta). Differentiation gives

dr^dt=θ˙θ^,dθ^dt=−θ˙r^.\frac{d\hat{\mathbf r}}{dt}=\dot\theta\hat{\boldsymbol\theta},\qquad \frac{d\hat{\boldsymbol\theta}}{dt}=-\dot\theta\hat{\mathbf r}.(17.1)

Since r=rr^\mathbf r=r\hat{\mathbf r}, velocity is r˙=r˙r^+rθ˙θ^\dot{\mathbf r}=\dot r\hat{\mathbf r}+r\dot\theta\hat{\boldsymbol\theta}. Differentiate each term:

r¨=r¨r^+r˙θ˙θ^+(r˙θ˙+rθ¨)θ^−rθ˙2r^.\ddot{\mathbf r}=\ddot r\hat{\mathbf r}+\dot r\dot\theta\hat{\boldsymbol\theta} +(\dot r\dot\theta+r\ddot\theta)\hat{\boldsymbol\theta}-r\dot\theta^2\hat{\mathbf r}.(17.2)

Collecting components gives

a=(r¨−rθ˙2)r^+(rθ¨+2r˙θ˙)θ^.\boxed{\mathbf a=(\ddot r-r\dot\theta^2)\hat{\mathbf r}+(r\ddot\theta+2\dot r\dot\theta)\hat{\boldsymbol\theta}.}(17.3)

Newton's law supplies μ(r¨−rθ˙2)=F(r)=−V′(r)\mu(\ddot r-r\dot\theta^2)=F(r)=-V'(r) and μ(rθ¨+2r˙θ˙)=0\mu(r\ddot\theta+2\dot r\dot\theta)=0. Multiplying the angular equation by rr identifies d(μr2θ˙)/dt=0d(\mu r^2\dot\theta)/dt=0, in agreement with the torque proof.

3. Radial equation and effective potential

Substitute θ˙=ℓ/(μr2)\dot\theta=\ell/(\mu r^2) in the radial equation:

μr¨=ℓ2μr3−V′(r).\mu\ddot r=\frac{\ell^2}{\mu r^3}-V'(r).(17.4)

Define the effective potential

Veff(r)=V(r)+ℓ22μr2.\boxed{V_{\rm eff}(r)=V(r)+\frac{\ell^2}{2\mu r^2}.}(17.5)

Its derivative is Veff′=V′−ℓ2/(μr3)V'_{\rm eff}=V'-\ell^2/(\mu r^3), so μr¨=−Veff′\mu\ddot r=-V'_{\rm eff}. The additional term is the angular kinetic energy after fixing angular momentum. Bringing the particle inward at fixed ℓ\ell requires faster angular motion and therefore more angular kinetic energy. The resulting radial barrier is often called the centrifugal barrier. It is not an additional physical outward interaction in the inertial-frame vector equation.

The conserved relative energy is

E=12μr˙2+12μr2θ˙2+V(r)=12μr˙2+Veff(r).E=\frac12\mu\dot r^2+\frac12\mu r^2\dot\theta^2+V(r) =\boxed{\frac12\mu\dot r^2+V_{\rm eff}(r).}(17.6)

Allowed radii satisfy E≥Veff(r)E\geq V_{\rm eff}(r) because r˙2≥0\dot r^2\geq0. The corresponding radial phase curves are pr=±2μ[E−Veff(r)]p_r=\pm\sqrt{2\mu[E-V_{\rm eff}(r)]}.

r/pEnergy / (k/p)Veff: solid blueangular term: brownV = −1/r: dashed135E−½Minimum at r/p = 1; two turning points for −½ < E < 0.

E = −0.30: two radial turning points, at r/p ≈ 0.613 and 2.721.

Figure 17.1 · In units with μ = k = ℓ = 1, the effective potential is 1/(2r²) − 1/r. The horizontal line selects energy, and intersections are radial turning points.

4. Turning points and accessible regions

At a radial turning point, r˙=0\dot r=0 and E=VeffE=V_{\rm eff}. A simple root of that equation is an ordinary turning point: the slope of the potential supplies acceleration back into the allowed region. The particle's total speed is not zero there if ℓ≠0\ell\ne0; it still moves tangentially.

Several roots may exist for a complicated potential. The initial condition determines the connected allowed radial interval in which the motion lies. One cannot jump across a forbidden region where E<VeffE<V_{\rm eff}. A double root at an extremum can describe a circular orbit or asymptotic approach to an unstable circle, rather than an ordinary reversal.

If the radius stays in a finite interval, the orbit is radially bounded. If an accessible interval extends to arbitrarily large rr, escape is possible. The sign of energy alone decides neither question for a general potential: adding a constant to VV changes the sign without changing motion. The familiar E<0E<0 rule will apply to attraction −k/r-k/r with V(∞)=0V(\infty)=0.

Time can be reconstructed by separating the radial energy equation:

t−t0=±μ2∫r0rdsE−Veff(s).t-t_0=\pm\sqrt{\frac\mu2}\int_{r_0}^{r}\frac{ds}{\sqrt{E-V_{\rm eff}(s)}}.(17.7)

Choose the sign on each inward or outward branch. Near an ordinary turning point the square-root singularity is integrable. Near a double root it can produce infinite approach time.

5. Circular orbits and local radial stability

A circle has r=rcr=r_c constant, so both r˙\dot r and r¨\ddot r vanish. Therefore

Veff′(rc)=0,Ec=Veff(rc),V′(rc)=ℓ2μrc3.V'_{\rm eff}(r_c)=0,\qquad E_c=V_{\rm eff}(r_c),\qquad V'(r_c)=\frac{\ell^2}{\mu r_c^3}.(17.8)

To test radial stability at fixed ℓ\ell, set r=rc+ηr=r_c+\eta with small η\eta. Taylor expansion gives Veff′(rc+η)=Veff′(rc)+Veff′′(rc)η+O(η2)V'_{\rm eff}(r_c+\eta)=V'_{\rm eff}(r_c)+V''_{\rm eff}(r_c)\eta+O(\eta^2). The first term is zero, so

μη¨=−Veff′′(rc)η.\mu\ddot\eta=-V''_{\rm eff}(r_c)\eta.(17.9)

If Veff′′>0V''_{\rm eff}>0, radial disturbances oscillate with ωr2=Veff′′/μ\omega_r^2=V''_{\rm eff}/\mu. If Veff′′<0V''_{\rm eff}<0, they include exponentially growing solutions. If it is zero, the linear test is inconclusive and higher-order terms must be examined. Stability here concerns the circular orbit's radial behavior; a perturbed orbit can accumulate an angular phase difference without its radius diverging.

6. Applications

Example 1 · The inverse-square effective potential

For V=−k/rV=-k/r, k>0k>0,

Veff=ℓ22μr2−kr,Veff′=−ℓ2μr3+kr2.V_{\rm eff}=\frac{\ell^2}{2\mu r^2}-\frac kr,\qquad V'_{\rm eff}=-\frac{\ell^2}{\mu r^3}+\frac{k}{r^2}.(17.10)

The circular radius is rc=ℓ2/(μk)r_c=\ell^2/(\mu k). The second derivative is 3ℓ2/(μr4)−2k/r33\ell^2/(\mu r^4)-2k/r^3, which at rcr_c becomes k/rc3>0k/r_c^3>0. Thus the circle is radially stable. Its energy is Ec=−k/(2rc)=−μk2/(2ℓ2)E_c=-k/(2r_c)=-\mu k^2/(2\ell^2).

As r→0r\to0 the positive r−2r^{-2} term dominates, so nonzero-angular-momentum motion cannot hit the center at finite energy. As r→∞r\to\infty, Veff→0V_{\rm eff}\to0 from below. Energies Ec<E<0E_c<E<0 allow two finite turning points; E≥0E\geq0 allows arbitrarily large radius.

Example 2 · Turning radii from energy

Use units in which μ=k=ℓ=1\mu=k=\ell=1, so Veff=1/(2r2)−1/rV_{\rm eff}=1/(2r^2)-1/r. Let E=−3/8E=-3/8. Multiplying E=VeffE=V_{\rm eff} by 8r28r^2 gives −3r2=4−8r-3r^2=4-8r, or 3r2−8r+4=03r^2-8r+4=0. Therefore

r=8±64−486=8±46,rmin⁡=23,rmax⁡=2.r=\frac{8\pm\sqrt{64-48}}6=\frac{8\pm4}{6},\qquad r_{\min}=\frac23,\quad r_{\max}=2.(17.11)

The allowed interval lies between these roots because the potential is below EE there. At either root, tangential speed is vθ=ℓ/(μr)=1/rv_\theta=\ell/(\mu r)=1/r, so the particle still moves. The inner-point speed is three times the outer-point speed.

Example 3 · Isotropic spring attraction

For V=κr2/2V=\kappa r^2/2, κ>0\kappa>0, the circular condition gives κrc=ℓ2/(μrc3)\kappa r_c=\ell^2/(\mu r_c^3), hence rc4=ℓ2/(μκ)r_c^4=\ell^2/(\mu\kappa). The second derivative is κ+3ℓ2/(μrc4)=4κ\kappa+3\ell^2/(\mu r_c^4)=4\kappa, so ωr=2κ/μ\omega_r=2\sqrt{\kappa/\mu}.

The circular angular frequency is Ω=ℓ/(μrc2)=κ/μ\Omega=\ell/(\mu r_c^2)=\sqrt{\kappa/\mu}. There are two radial oscillations per revolution near a circle. The large-rr potential grows without bound, so every finite-energy regular orbit is bounded; positive energy does not imply escape in this potential.

Example 4 · An unstable circular orbit

Take V=−a/r3V=-a/r^3 with a>0a>0. Then Veff′=3a/r4−ℓ2/(μr3)V'_{\rm eff}=3a/r^4-\ell^2/(\mu r^3), so rc=3μa/ℓ2r_c=3\mu a/\ell^2. At this radius,

Veff′′=−12arc5+3ℓ2μrc4=−3arc5<0.V''_{\rm eff}=-\frac{12a}{r_c^5}+\frac{3\ell^2}{\mu r_c^4}=-\frac{3a}{r_c^5}<0.(17.12)

The circle lies at a maximum. A small outward displacement is accelerated farther outward, and a small inward displacement inward. Here the attractive r−3r^{-3} potential dominates the centrifugal r−2r^{-2} term near the center, so a nonzero ℓ\ell does not guarantee a collision barrier for every potential.

7. Radial motion with zero angular momentum

If ℓ=0\ell=0, the effective potential is simply V(r)V(r). There is no angular barrier and no meaningful angular parameter for the orbit equation. In inverse-square attraction, inward radial motion can reach the singular center in finite time. The ideal point-particle force law supplies no regular continuation through that collision; a realistic collision model would be additional physics.

Exercises

Hints and solutions are collected separately.

  1. 17.1 Derive plane-polar acceleration by differentiating the basis vectors.

  2. 17.2 Derive the effective potential at fixed angular momentum.

  3. 17.3 Is the centrifugal term a new outward physical interaction?

  4. 17.4 For μ=k=ℓ=1\mu=k=\ell=1, find the inverse-square circular radius and energy.

  5. 17.5 For the same parameters, find turning radii at E=−2/9E=-2/9.

  6. 17.6 At a radial turning point with ℓ≠0\ell\ne0, is the total velocity zero?

  7. 17.7 Derive the local radial stability criterion for a circle.

  8. 17.8 For V=κr2/2V=\kappa r^2/2, find the circular radius at fixed ℓ\ell.

  9. 17.9 Why does energy sign alone not classify general central-force orbits?

  10. 17.10 For V=−a/r3V=-a/r^3, determine whether the centrifugal barrier always prevents collision.

  11. 17.11 Write the radial time integral and explain its sign.

  12. 17.12 Explain why a double root of E=VeffE=V_{\rm eff} requires special care.

Reading progress is saved on this device.

Search the notes

Search all four units and the course references.