Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 62 / 4

For a spherical system, use the spherical shell theorem with shell mass $dm=4\pi \rho (s)s^{2}ds$. Inner shells contribute $-Gdm/r$ to the gravitational potential, while an outer shell contributes its constant interior potential $-Gdm/s$. With $\Phi (\infty )=0$, assuming the required integrals converge,On differentiation, the two terms containing $\rho (r)r$ cancel. Hence ${\Phi }^{\prime }=GM(r)/r^2$, where $M(r)=4\pi {\int }_0^r\rho (s)s^{2}ds$. Radial balance for a circular orbit gives $v_c^2=r{\Phi }^{\prime }=GM(r)/r$.

For a razor-thin axisymmetric astrophysical disk, the element of mass is $dm=\Sigma (s)sdsd{\varphi }^{\prime }$, where $\Sigma$ is the surface density. Superposing Newtonian gravitational potentials therefore givesAxisymmetry permits $\varphi =0$. Unlike the spherical case, an exterior annulus generally exerts a radial force: enclosed mass does not determine a disc rotation curve.

For the Legendre expansion of thin-disk gravity, split the radial integral at $s=R$. In the inner part the kernel expands in $s/R$, while in the outer part it expands in $R/s$. The angular integrals of odd Legendre polynomials vanish, and those of even degree $n$ equal $2{\alpha }_n$, whereIntroduce $I_n(R)={\int }_0^R\Sigma (s)s^{n+1}ds$ and $J_n(R)={\int }_R^{\infty }\Sigma (s)s^{-n}ds$. The gravitational potential becomesWhen differentiating each bracket, the moving-limit terms cancel: $R^{-(n+1)}{I}_n^{\prime }=\Sigma (R)$ and $R^n{J}_n^{\prime }=-\Sigma (R)$. Since $v_c^2=R{\Phi }^{\prime }$, this givesThe zeroth term has ${\alpha }_0=\pi$ and $M(R)=2\pi I_0$. Separating it provesThe inner correction is inward, while the exterior correction is outward. For a smooth surface density, the original potential singularity at coincident points is integrable. Its radial force is understood through a symmetric Cauchy principal value or a vanishing-thickness regularization; the paired interior and exterior terms above retain the cancellation at $s=R$. This avoids treating the two singular local force contributions separately.

For an exponential galactic disk, write $\Sigma (s)={\Sigma }_0{e}^{-s/R_d}$, so $M_{\infty }=2\pi {\Sigma }_0{R}_d^2$ and $M(R)=M_{\infty }[1-(1+R/R_d)e^{-R/R_d}]$. At large $R$, the missing mass and exterior-ring terms are exponentially small. The leading nonspherical interior term is $k=1$, with ${\alpha }_2=\pi /4$ and $I_2(\infty )=6{\Sigma }_0{R}_d^4$. Consequently the exponential-disk Keplerian asymptotic isThe positive leading correction shows that the rotation curve approaches the Keplerian limit from above. The finite-order large-radius expansion is sufficient here; an infinite moment expansion need not converge for a disk extending to infinity.

A useful special example is the Mestel disk, with surface density $\Sigma (s)=C/s$ for $s>0$, $C>0$. It has $M(R)=2\pi CR$. For every positive even $n$,All the nonspherical corrections cancel, givingThus the Mestel disk has a perfectly flat galaxy rotation curve.

This example has infinite total mass and a singular central surface density. The absolute gravitational potential cannot be set to zero at infinity, but the radial force exists as a limit of disks with increasing outer cutoff. Potential differences are $\Phi (R)-\Phi (R_0)=v_M^2\log (R/R_0)$. The earlier potential integral is therefore interpreted up to a radius-independent divergent constant for this example; the force calculation remains valid. Truncating the Mestel disk gives a more physical finite system but changes the exact flat curve near its edges.