Past exam of the physics B course of the University of Cambridge / 2013 / Paper 1 / Section C / C11 / Solution

The shell theorem gives enclosed halo mass $4\pi \rho r^3/3$. Define $\alpha =4\pi \rho /3$. Then $F_r=-GMm/r^2-\alpha Gmr$. Interior to the uniform halo its potential energy may be writtenOuter shells alter the additive constant, not this interior force. An actually infinite uniform mass distribution has no finite potential referenced to infinity; only local potential differences are needed here.

The central force has zero torque, so $\mathbf{J}=\mathbf{r}\times m\dot{\mathbf{r}}$ is conserved and the motion is planar. With $J=m{r}^2\dot{\varphi }$, the effective radial potential and equation areFor circular radius $r_0$, $J=m{r}_0^2\Omega$ givesSet $r=r_0+\xi$ at fixed $J$. The derivative of the right-hand side at equilibrium is $-3{\Omega }^2+2GM/r_0^3-\alpha G$, soBoth frequencies are positive, but radial and azimuthal motion no longer close with the Kepler ratio. At first order in eccentricity $e$, one may write $r=r_0(1-e\cos \omega t)$ and $\varphi =\Omega t+(2\Omega e/\omega )\sin \omega t$. For a weak halo this is a nearly Keplerian ellipse whose apsis rotates at angular rate $\Omega -\omega$.

Put ${\Omega }_K=\sqrt{GM/r_0^3}$ and $\epsilon =\alpha r_0^3/M\ll 1$. Expansion yieldsThe precession is retrograde. In cycles per unit time the precession frequency is $\dot{\varpi }/(2\pi )$; per radial period the apsidal shift is $2\pi (\Omega /\omega -1)\simeq -3\pi \epsilon$. This weak-halo interpretation avoids confusing a focus-centred Kepler ellipse with the centred ellipse of a purely harmonic potential.