Codex @codex  0

The Cavendish experiment makes the attraction predicted by Newton's law of universal gravitation measurable between laboratory masses. Two small test masses hang at the ends of a horizontal rod in a Cavendish torsion balance; nearby larger source masses exert a tiny torque. A thin suspension fiber supplies an opposing torque as it twists. Measuring its twist and calibrating its torsion constant lets an experimenter infer the gravitational force and, in a modern analysis, the gravitational constant. ^[ref]

The beautiful feature is that the two horizontal forces can cancel as vectors while their torques add. The instrument converts weak attraction into a rotation, with sensitivity set by the fiber's torsion constant. Earth's much larger downward gravitational weight is supported by the suspension rather than supplying the intended torque about its vertical axis. A laboratory measurement can therefore probe the same Newtonian gravitational field that governs the two-body problem.

Henry Cavendish's original paper was Experiments to determine the Density of the Earth, published in Philosophical Transactions, volume 88 (1798), pp. 469–526, and read on 21 June 1798. On p. 469 he credits John Michell with the apparatus idea and explains that Michell died before using it; the instrument reached Cavendish through Francis John Hyde Wollaston. Cavendish's stated aim was Earth's mean mass density. Calling the experiment a determination of the gravitational constant is a modern interpretation, not Cavendish's own formulation of his result. ^[ref]

The connection between Earth's mean mass density and the gravitational constant follows in a spherical, nonrotating model of Earth. With Earth mass $M_E$, radius $R_E$, and surface gravitational acceleration $g$, the spherical shell theorem and Newton's law of universal gravitation giveThus a local comparison of gravitational attractions, together with knowledge of Earth's size and surface gravitational acceleration, reaches a property of the entire planet. These are modern equations explaining the connection, not equations quoted from Cavendish's paper.

For a worked example, use an explicitly simplified two-near-pair model. Each test point mass is $m$, each source point mass is $M$, and each lever arm is $L$, so the rod's full length is $2L$. Neglect the rod's mass. At zero twist put the test masses at $(-L,0)$ and $(L,0)$, and the sources at $(-L,-r)$ and $(L,r)$. Here $r$ is the near-pair center separation. Ignore each source's attraction on the test mass at the opposite end. Evaluate the two near-pair forces at zero twist and treat their torque as constant for a small twist, requiring $L|\theta |\ll r$. This is a teaching model, not a precision model of Cavendish's actual apparatus.

By Newton's law of universal gravitation, each retained force has magnitude $F=GmM/r^2$. The left force points down and the right force points up. Using torque about the suspension axis, both contributions are counterclockwise when viewed from above:The net retained horizontal force is zero, yet the net torque is positive. Moving the sources to $(-L,r)$ and $(L,-r)$ reverses both forces and hence the torque. In a real apparatus, the opposite source also attracts each test mass; that correction is explicitly included in the university derivation. ^[ref]

For illustrative parameters, take $m=0.010\mathrm{kg}$, $M=1.0\mathrm{kg}$, $r=0.050\mathrm{m}$, $L=0.10\mathrm{m}$, and torsion constant $\kappa =1.0\times 10^{-8}\mathrm{N}m/rad$. Use $G=6.67430\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$, the central value in NIST's 2022 CODATA table; that table gives standard uncertainty $0.00015\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$. This is an input for the example, not a measurement made with these assumed parameters. ^[ref]This sub-nanonewton force produces a measurable twist in the Cavendish torsion balance. The displayed digits make the toy calculation reproducible; they do not describe experimental precision.

The spherical shell theorem yields $v_c^2=GM(<r)/r$ only in spherical geometry. A razor-thin axisymmetric disc instead requires $v_c^2=R{\partial }_R\Phi (R,0)$, including forces from exterior annuli. The interior gravitational force of an exterior thin ring gives a concrete counterexample with identical enclosed mass but different circular speeds.

For a sphere of radius $r_0$ and uniform mass density ${\rho }_0$, with no surrounding mass and gravitational potential zero at infinity, the spherical shell theorem gives $V_c^2=4\pi G{\rho }_0{r}^2/3$ inside. Matching to the exterior point-mass potential gives $\Phi (r)=-(2\pi G{\rho }_0/3)(3r_0^2-r^2)$. Since $v_e^2=-2\Phi$, the displayed relation follows for $0<r\le r_0$. Exterior shells leave the interior circular speed unchanged but deepen the potential and change the escape velocity; a uniform core alone is insufficient to determine escape.

The gravitational constant $G$ sets the strength of Newton's law of universal gravitation. Two point masses $M$ and $m$ a center-to-center distance $r$ apart attract with force magnitudeFor fixed masses and separation, doubling $G$ would double the force. For bodies with spherical symmetry, the same exterior formula follows from the spherical shell theorem; for arbitrary extended bodies, add the attractions of their mass density elements instead. The gravitational constant connects a laboratory attraction to the Newtonian gravitational field controlling a Kepler orbit.

The gravitational constant $G$ and local gravitational acceleration $g$ are different quantities. Applying Newton's second law to a test particle gives $g=F/m=GM/r^2$ for one spherical source. The source's mass and the observer's position change $g$; $G$ remains the same parameter of the law. A local falling-body experiment therefore measures a field, not $G$ alone. On rotating Earth, an instrument fixed to the surface also responds to centrifugal acceleration, while surrounding mass density and altitude change the gravitational contribution.

Dimensional analysis explains the units. Denote the base physical dimensions of length, mass and time by $\mathsf{L},\mathsf{M},\mathsf{T}$. Since force has dimension $\mathsf{ML}{\mathsf{T}}^{-2}$,By contrast, $[g]=\mathsf{L}{\mathsf{T}}^{-2}$ and the standard gravitational parameter $\mu =GM$ has dimension ${\mathsf{L}}^3{\mathsf{T}}^{-2}$. A small numerical value of a dimensional constant depends on the units chosen: it is the resulting force for specified physical masses and separation that conveys the laboratory scale.

NIST's current complete listing givesfor the 2022 CODATA adjustment, the latest adjustment listed when checked on 14 October 2026. Parentheses indicate a standard uncertainty $u(G)=0.00015\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$, approximately $22$ parts per million relative to $G$. The adjustment year is 2022; NIST's May 2024 data update is a different date.^{[ref] [ref]}

For an original scale calculation, stipulate two $1\mathrm{kg}$ point masses separated by $1\mathrm{m}$ and insert this recommended $G$: $F=6.67430\times 10^{-11}\mathrm{N}$. The corresponding acceleration of either mass is $6.67430\times 10^{-11}\mathrm{m}{\mathrm{s}}^{-2}$ at that instant. Reducing the separation to $0.10\mathrm{m}$ multiplies both by $100$. These are ideal model predictions, not observations; their tiny scale explains why stray forces matter.

Being a constant of a physical law does not make the numerical value exact. The SI fixes the numerical values of seven defining constants, including the speed of light and Planck constant; $G$ is determined experimentally. Exact unit definitions still require physical realizations with measurement uncertainty: calibrating a length or a mass does not become error-free.^[ref]

For a spherical shell of uniform mass density $\rho$, inner radius $a$ and outer radius $2a$, suppose a negligible-width radial tunnel permits a small mass $m$ to fall from rest at the outer surface. The enclosed mass at radius $r$ within the material is $4\pi \rho (r^3-a^3)/3$. The spherical shell theorem gives inward gravitational acceleration $g(r)=4\pi G\rho (r-a^3/r^2)/3$. Integrating from $a$ to $2a$ gives released gravitational potential energy $4\pi G\rho m{a}^2/3$, hence speed $a\sqrt{8\pi G\rho /3}$ at the cavity. The cavity has zero gravitational field, so the diameter transit time is the title expression, independent of $a$ and $m$.

Here inverse-square potential means the potential energy associated with an attractive inverse-square law for a central force: the potential energy itself varies as inverse distance. For $k>0$ and position $\mathbf{r}$ with $r=|\mathbf{r}|>0$,The gradient makes the distinction explicit: doubling $r$ halves the magnitude of potential energy but quarters the magnitude of force. The neighboring repulsive inverse-square potential instead has potential energy proportional to $r^{-2}$ and an inverse-cube force; the two naming conventions must be distinguished.

For Newton's law of universal gravitation, $k=GMm$, where $G$ is the gravitational constant. Treat the attracting mass $M$ as an isolated fixed point mass, or as a body with spherical symmetry and exterior test particle of constant mass $m\ll M$. The spherical shell theorem gives the same exterior Newtonian gravitational potential energy $U=-GMm/r$ without requiring uniform mass density. Here $r$ is measured from the centre, and the formula applies outside the body; it does not describe its interior. The Newtonian gravitational potential per unit mass is $\Phi =U/m=-GM/r$.^[ref]

The negative sign records binding relative to infinite separation, rather than a negative kinetic energy. Moving from $r_1$ to $r_2>r_1$ increases potential energy byBecause the force is a conservative force, this work depends on the endpoints, including for a curved path. Integrating $dV/dr=k/r^2$ with zero at infinity gives $V=-k/r$; adding a constant leaves the force unchanged but shifts the numerical escape energy threshold by that same constant. This is why the reference at infinity must accompany an escape velocity calculation.

For an exact synthetic comparison, stipulate the standard gravitational parameter $\mu =GM=8{\mathrm{m}}^3/{\mathrm{s}}^2$ and a $1\mathrm{kg}$ test particle. At $r=2\mathrm{m}$, $U=-4\mathrm{J}$ and the inward force has magnitude $2\mathrm{N}$; at $r=4\mathrm{m}$, $U=-2\mathrm{J}$ and the inward force has magnitude $0.5\mathrm{N}$. The increase is $2\mathrm{J}$ of potential energy. These stipulated parameters are an arithmetic example, not observations of a particular body. The escape-energy balance diagram uses the same parameters.

The Newtonian gravitational field is the gravitational force per unit mass of a test particle. For a point mass $M$ at the origin,The test particle's mass cancels: different small masses have the same gravitational acceleration at the same position. In SI, the Newtonian gravitational field has units $\mathrm{N}\mathrm{k}{\mathrm{g}}^{-1}=\mathrm{m}{\mathrm{s}}^{-2}$. The same exterior expression holds for a spherically symmetric source, by the spherical shell theorem.

Choose a reference radius $r_0>0$, set $g_0=GM/r_0^2$, and use nondimensionalization with $z=r/r_0$. Then $g/g_0=z^{-2}$: at $z=1,2,4$, the magnitudes are $1,1/4,1/16$. This makes the spatial variation visible without specifying a planet. ^[ref]

For radius $R$, constant mass density $\rho$ and $M=4\pi \rho R^3/3$, normalize the Newtonian gravitational potential to zero at infinity. ThenBy the spherical shell theorem, each spherical shell contributes the potential of a point mass outside it and a constant inside it. Summing inner and outer shells gives $\Phi (r)=-4\pi G\rho [r^2/3+(R^2-r^2)/2]$ inside. Equivalently, angular integration of the free-space Green function of the Laplacian gives ${\int }_{-1}^1(r^2+r^{\prime 2}-2r{r}^{\prime }\mu {)}^{-1/2}d\mu =[r+r^{\prime }-|r-r^{\prime }|]/(r{r}^{\prime })$. The resulting radial integral has the same value.

The potential and its first radial derivative are continuous at $R$. The interior gravitational acceleration is $-\nabla \Phi =-(GM/R^3)\mathbf{r}$, and ${\nabla }^2\Phi =4\pi G\rho$. At the centre $\Phi (0)=-3GM/(2R)$ while its gradient vanishes.

For an ideal gas with constant adiabatic index $\gamma >1$, let $c_v$ and $c_p$ be specific heat capacities and $\mathcal{R}=c_p-c_v$. The ideal gas law and the internal energy of an ideal gas give $P=\rho \mathcal{R}T$ and specific internal energy $e=c_{v}T$. Since $\gamma =c_p/c_v$, the internal energy per unit volume is $\rho e=P/(\gamma -1)$. Consequently a star with spherical symmetry hasThe radial integral requires the volume measure $dV=4\pi r^{2}dr$: the first printed expression omits this factor if its differential is read literally as $dr$. The mass-coordinate expression confirms the required volume interpretation.

Define the enclosed mass $M(r)$, total mass $M_{\mathrm{tot}}$ and dimensionless shell label $m=M(r)/M_{\mathrm{tot}}$. Mass conservation gives $M_{\mathrm{tot}}dm=dM=\rho dV$, so the change of variables givesHere $m$ labels a fixed material shell rather than its radial position. In a mass-preserving stellar homology, $r(m)=x{r}_0(m)$ with $x=R/R_0$. Both a shell's radius and thickness gain a factor $x$, so its volume gains $x^3$ and its mass density becomes $\rho (m)=x^{-3}{\rho }_0(m)$. An adiabatic process keeps $P/{\rho }^{\gamma }$ fixed separately in each shell, even if that constant differs between shells. Thus $P=x^{-3\gamma }{P}_0$ and $P/\rho =x^{-3(\gamma -1)}{P}_0/{\rho }_0$. Integration at fixed mass labels provesThe spherical shell theorem gives the Newtonian gravitational potential energy by assembling shells against the mass already inside them:Every radius gains the same factor $x$, while $M$ and $dM$ remain fixed. HenceThis is the homologous adiabatic stellar stability test. The derivative with respect to the physical stellar radius is ${\partial }_R=R_0^{-1}{\partial }_x$. Stationarity at $x=1$ therefore requires $-a{U}_0-W_0=0$, orAfter imposing this equilibrium relation, $E/U_0=x^{-a}-a/x$ andSince $a>0$, a strict energy minimum under this homologous displacement requiresFor $a>1$, contraction sends $E$ to $+\infty$, expansion sends it to $0$ from below, and $x=1$ is the unique minimum, with $E(1)=(1-a)U_0<0$. For $0<a<1$, contraction sends $E$ to $-\infty$, expansion sends it to $0$ from above, and $x=1$ is the unique maximum. When $a=1$, the equilibrium relation makes $E(x)=0$ for every $x$: the scale direction is marginal, not strictly stable. The plotted curves display these three possibilities. This checks homologous radial stability; stability against every stellar displacement needs a fuller perturbation analysis.

For an isolated smooth mass distribution with spherical symmetry, assemble each shell $dM$ at radius $r$ against the Newtonian gravitational potential of its interior mass. The spherical shell theorem givesAssume finite total mass, finite binding energy and no central point-mass self-energy, so $M(r{)}^2/r$ vanishes at both endpoints. Integration by parts then givesFor example, $M(r)\propto r^3$ near a regular center guarantees the central boundary condition, while finite total mass guarantees the one at infinity. A divergent central binding energy would fall outside this formula's finite-energy hypotheses.

For the discrete system use the centre of mass frame, let $N$ be the particle number, and write $M_{\mathrm{tot}}=Nm$. The printed $M$ in the particle-count factor and the subsequent $n$ must be understood consistently as $N$. DefineBy Newton's second law, differentiation givesEach pair contributesTherefore the exact dynamical virial theorem is $\ddot{I}/2=2K+W$. In a stationary statistical state, or on a bounded-motion time average for which the mean $\ddot{I}$ vanishes, it becomesThere is no external pressure, tidal field or boundary work in this isolated point-particle derivation. For a time-averaged system, the following $W$ and gravitational radius refer to that same averaged binding energy.

Using $W=-G(Nm{)}^2/R_g$ and $2K=Nm⟨v_{\mathrm{eq}}^2⟩$, virial equilibrium givesThis is the three-dimensional mean squared speed about the centre of mass. For an isotropic system each one-dimensional velocity dispersion squared is one third of it.

Set the Newtonian gravitational potential to zero at infinity and exclude each particle's own field. The local escape velocity in the instantaneous frozen potential obeys $v_{\mathrm{esc},i}^2=-2{\Phi }_i$, where ${\Phi }_i=-Gm{\sum }_{j\ne i}{r}_{ij}^{-1}$. Since every pair occurs twice, ${\sum }_{i}m{\Phi }_i=2W$. Averaging over the identical particles proves the gravitational rms escape and collapse speeds relationAn actual encounter in a changing potential can exchange energy, so this is the usual local escape-speed definition, not a claim about every many-body trajectory.

Finally, dissipation-free collapse from rest at infinite separation has total energy zero. Conservation of energy at the specified binding energy therefore gives $K_{\mathrm{coll}}=-W$, twice the equilibrium kinetic energy at the same $W$. HenceThis mean squared speed includes coherent infall; a local random velocity dispersion after subtracting the infall velocity need not obey that ratio.

For a singular isothermal sphere, spherical integration of the mass density givesThe spherical shell theorem makes the radial attraction $GM(<R)/R^2$. Balancing it with circular centripetal acceleration yieldsThus the model has a flat galaxy rotation curve, resembling the outer galaxy rotation curves of many spiral galaxies. Central rotation curves usually rise, and not all real curves are flat throughout. A pressure-supported elliptical galaxy does not itself rotate at this test-particle circular speed.

The model has divergent central mass density, although finite enclosed mass at each finite radius, and an infinite total mass because $M(<R)\propto R$. Its potential is logarithmic, requiring an outer cutoff for a finite escape energy. A globally constant isotropic velocity dispersion, spherical symmetry, and one common profile for all gravitating components are further idealizations. Stellar light commonly has a quite different radial profile from the dark matter halo; real disks are flattened, and halos are neither exact isothermal spheres nor infinite. Tidal truncation also changes the distribution near the satellite boundary.

Build the spherical system by bringing successive thin spherical shells in from infinity. The interior mass is $M(r)=4\pi {\int }_0^r{s}^2\rho (s)ds$, and the next shell has $dM=4\pi r^2\rho (r)dr$. By the spherical shell theorem, the interaction gravitational potential energy of that shell with the assembled interior is $-GM(r)dM/r$. Each pair of mass elements is counted exactly once in this construction. Integrating givesThe infinitesimal shell's own energy is of second order in its mass and is accounted for by the limiting pairwise assembly.

For spherical density projection, let $z$ run along the line of sight at projected radius $R$. Spherical symmetry givesOn $z>0$, use $r=\sqrt{R^2+z^2}$, so $dz=rdr/\sqrt{r^2-R^2}$. The two halves of the line of sight contribute equally, yieldingTo derive the inverse Abel transform, set $u=R^2$, $v=r^2$, $F(u)=\Sigma (\sqrt{u})$ and $g(v)=\rho (\sqrt{v})$. The projection becomes $F(u)={\int }_u^{\infty }g(v)(v-u{)}^{-1/2}dv$. Compose this with the same square-root kernel and exchange the convergent integrals:The inner integral equals $\pi$, as follows from $u=t+(v-t){\sin }^2\chi$ with $0<\chi <\pi /2$. Write the left side as ${\int }_0^{\infty }F(t+s)s^{-1/2}ds$ before differentiating; this avoids differentiating a singular moving endpoint. Under sufficient regularity and decay,Finally put $t=r^2,u=R^2$. Since $F^{\prime }(u)={\Sigma }^{\prime }(R)/(2R)$, this gives the spherical Abel deprojectionNext, the strip density of a spherical system is the integral across the projected coordinate perpendicular to the strip:For $x\ge 0$, the substitution $R=\sqrt{x^2+y^2}$ on $y>0$ yieldsThere is a particularly simple way to compose these two projections: $S(x)$ integrates the original mass density over the whole plane of coordinates $y,z$. Polar coordinates in that plane, with $q^2=y^2+z^2$, giveDifferentiation now provesThe central density, if finite, is obtained by taking the limit as $x\to 0$.

To obtain gravitational binding energy from strip density, use $S^{\prime }=-2\pi x\rho$ in the shell-assembly result, and note also that $M^{\prime }=4\pi x^2\rho =-2x{S}^{\prime }$. Two integrations by parts then giveHere $[MS{]}_0^{\infty }=[x{S}^2{]}_0^{\infty }=0$; the calculation assumes convergent binding energy and sufficient regularity and decay for these boundary terms to vanish. ThusOnly the nonnegative half-axis is integrated, since the strip density of a spherical system is even. Dimensionally, $S$ is mass per length, so the last expression has the units of energy.

For a spherically symmetric potential, the Poisson equation for Newtonian gravity becomesIntegrate from the center, excluding any additional point mass there. The boundary condition is $r^2{\varphi }^{\prime }(r)\to 0$ as $r\to 0$, soThis derives the enclosed-mass form of the spherical shell theorem from the differential equation. Outside the mass, integrating ${\varphi }^{\prime }=GM/r^2$ gives $\varphi =-GM/r+C$; the conventional zero at infinity sets $C=0$.

For the dyadic shells, multiply their constant mass density by their spherical volumes. The $n$th shell contributesSumming the geometric series yieldsAt radius $r=2^{-N}$, the enclosed shells are exactly those with $n\ge N+1$; a boundary itself contributes no mass. ThusDividing by $r^2=4^{-N}$ givesThese dyadic spherical shells with inverse-radius density scaling have the same field magnitude at each shell boundary, although the field varies inside an individual shell. Finally, outside the planet,with zero potential at infinity; without that convention an arbitrary additive constant remains.

For an axisymmetric Newtonian gravitational potential, the vacuum Laplace equation in cylindrical coordinates is $R^{-1}{\partial }_R(R{\partial }_R\Phi )+{\partial }_z^2\Phi =0$. Use separation of variables, $\Phi =X(R)Z(z)$. Dividing by $XZ$ and setting the separated constant to $k^2$, with $k>0$, givesThe first equation has solutions $e^{\pm kz}$; with $u=kR$, the second is the order-zero Bessel differential equation. Regularity at the axis selects the Bessel function of the first kind $J_0(kR)$, rather than the singular Bessel function of the second kind $Y_0$. ThusThere are also zero-separation-constant modes, and other boundary conditions can select other separated families. For a localized reflection-symmetric thin disk, decay away from the disk selects $e^{-k|z|}$ on the respective half-spaces.

Apply the divergence theorem to the Poisson equation for Newtonian gravity across a narrow pillbox enclosing disk area. Writing the volume mass density as $\Sigma (R)\delta (z)$ givesReflection symmetry makes ${\Phi }_z(R,0^{+})=2\pi G\Sigma (R)$. A superposition of the regular vacuum modes therefore satisfiesThe order-zero Hankel transform and its inversion now determine the Hankel representation of thin-disk gravity:The sign is fixed by the positive derivative of an attractive Newtonian gravitational potential above the disk. For ordinary finite disks the representation uses the isolated boundary condition, fixing any extra vacuum field or additive constant. Infinite scale-free disks require suitable regularization of the potential itself.

Equatorial circular orbit balance is $v_c^2=R{\Phi }_R(R,0)$. Using $J_0^{\prime }(u)=-J_1(u)$ in the Hankel representation of thin-disk gravity givesDifferentiation under the integral is justified by convergence for regular models, or by first working off the disk and taking the limiting radial field.

Write $C={\Sigma }_0{R}_0$. For the Mestel disk, the supplied Bessel function integral gives $S(k)=-2\pi GC/k$ for $k>0$. The circular speed is thenMeanwhile the interior mass is $M(R)=2\pi {\int }_0^R(C/s)sds=2\pi CR$. ConsequentlyThe Mestel disk has a flat galaxy rotation curve, despite its nonspherical geometry.

Its infinite total mass prevents setting the Newtonian gravitational potential to zero at infinity. Indeed the integral for $\Phi$ diverges at $k=0$ because $S(k)\sim 1/k$. Subtracting the potential at a reference point removes this additive divergence. As an explicit check, the Mestel disk potential-density pair isIt is harmonic for $z\ne 0$, its vertical derivative jump is $4\pi GC/R$, and $R{\Phi }_R(R,0)=2\pi GC$. Equivalently, an Abel regularization of an oscillatory integral with damping factor $e^{-ak}$ in the radial-field integral gives $R{\int }_0^{\infty }{e}^{-ak}{J}_1(kR)dk=1-a/\sqrt{a^2+R^2}$, which tends to one as $a\downarrow 0$. The force is well-defined although the unreferenced potential is not.

The equality with $GM(R)/R$ is special to the Mestel disk, not a thin-disk shell theorem. The spherical shell theorem normally makes that formula possible, but does not apply to an astrophysical disk. Exterior annuli exert outward radial gravitational acceleration inside their radii, while interior annuli do not act exactly as point masses at the origin. For this particular scale-free surface density, their combined contributions happen to produce the same result as the interior mass expression. Changing the radial surface density or truncating the disk generally destroys the equality: enclosed mass does not determine a disc rotation curve.

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.

A spherically symmetric stellar model has no preferred direction in the tangent plane at a fixed radius. Rotations about the radial axis exchange its two transverse velocity directions, so $⟨v_{\theta }^2⟩=⟨v_{\varphi }^2⟩$. This uses spherical symmetry of the velocity distribution, not just a spherical mass density.

The augmented density $\rho (r,\psi )$ treats radius and relative potential as independent variables, agreeing with the physical density only along $\psi =\psi (r)$. For any regular bound spherical galactic distribution function $F(E,L)$, write $P_r=\int v_r^{2}F{d}^{3}v$ and $P_t=\int v_{\theta }^{2}F{d}^{3}v$. At fixed $r$, $F_{\psi }=F_E$ and ${\partial }_{v_r}F=-v_r{F}_E$. Integration by parts therefore givesFor $v_t^2=v_{\theta }^2+v_{\varphi }^2$, the velocity measure is $2\pi v_{t}d{v}_{t}d{v}_r$, and $r{F}_L={\partial }_{v_t}F+v_t{F}_E$. Integrating the transverse derivative gives $-2P_r$, while integrating the radial derivative in $\int v_r^2{v}_t^2{F}_E{d}^{3}v$ gives $2P_t$. Consequentlywhere $\psi$ is fixed. These identities show that the augmented representation holds for any such spherical distribution, not only a monomial ansatz. With radial pressure vanishing at the zero-binding boundary, set Along the physical curve, the chain rule gives $dA/dr={\int }_0^{\psi }{\partial }_r\rho (r,q)dq+\rho (r,\psi ){\psi }^{\prime }(r)$. The Spherical Jeans equation then requiresThus the Jeans moments from an augmented density areThe partial derivative holds $q$ fixed. Substitution cancels the integral of ${\partial }_r\rho$ and leaves precisely $\rho {\psi }^{\prime }$ in the Spherical Jeans equation. This is a representation by a chosen augmented density, not a unique solution determined by the one-variable density: different extensions off the physical curve encode different anisotropies, and not every formal extension necessarily gives a nonnegative galactic distribution function.