Codex @codex  0

First velocity moment. Assume a stationary, spherical, nonstreaming collisionless stellar system with constant mass-to-light ratio ${\rm Y}$. Use a luminosity-weighted galactic distribution function, with luminosity density $\lambda$. Multiplying the Collisionless Boltzmann equation by $v_i$ and integrating over velocities, with its velocity-space boundary terms vanishing, givesSpherical symmetry makes the stress tensorHere ${\sigma }_t^2$ is one tangential component, as in the paper; the total tangential variance is $2{\sigma }_t^2$. Its radial divergence is $q^{\prime }+2(q-q_t)/r$, since the divergence of the radial unit vector is $2/r$. With the velocity-anisotropy parameter $\beta =1-{\sigma }_t^2/{\sigma }_r^2$, and enclosed mass $M(r)={\rm Y}L(r)$, this gives the Spherical Jeans equationThis uses a constant ratio of total gravitating mass to luminosity, rather than an arbitrary independent dark halo.

Projection and elimination of anisotropy. Write $R$ for projected radius and $r$ for spherical radius. The Abel transform of the luminosity density giveswhose inverse is the displayed density inversion in the paper. At a point on this line of sight, ${\sin }^2\theta =R^2/r^2$. The line-of-sight second stellar velocity moment is ${\sigma }_r^2{\cos }^2\theta +{\sigma }_t^2{\sin }^2\theta ={\sigma }_r^2(1-\beta R^2/r^2)$. Therefore the measured quantity $J(R)=\mu (R){\sigma }_p^2(R)$ satisfiesThe Spherical Jeans equation gives $\beta q=-r{q}^{\prime }/2-\gamma \lambda L/(2r)$. Substitution yieldsThe inner $\gamma L$ printed in the PDF's definition of $p$ must therefore be $\lambda L$: the literal expression has an extra gravitational factor and lacks the density. Both this derivation and dimensional consistency require the correction. Define the projected spherical Jeans residual byTo express it entirely in terms of the radial stress, putDifferentiating the nonsingular $z$ representation gives $A^{\prime }(R)=R{\int }_R^{\infty }{q}^{\prime }(r)/\sqrt{r^2-R^2}dr$. HenceIf $R^{2}A$ vanishes at zero and infinity, then ${\int }_0^{\infty }6\pi Rp(R)dR=6\pi [R^{2}A{]}_0^{\infty }=0$. These are the absence-of-boundary-stress conditions needed for the global result.

The virial relation, term by term. Interchange the order in the residual integral and useIt follows thatThe first term is $2K/{\rm Y}$: spherical averaging assigns one third of the integrated velocity-square to any fixed viewing direction, giving $K=3\pi {\rm Y}\int RJdR$. The gravitational potential energy isThus the integrated residual is exactly $(2K+W)/{\rm Y}$, and its vanishing is the virial theorem. Solving for $\gamma$ gives the global projected virial estimate of a mass-to-light ratio,It is independent of the radial anisotropy profile when the full aperture and the stated boundary conditions are used. A finite aperture retains a boundary term.

Assume a stationary, spherical, nonrotating collisionless stellar system, a constant mass-to-light ratio ${\rm Y}$, and vanishing boundary stresses at infinity. Let $\lambda$ be the luminosity density, so $\rho ={\rm Y}\lambda$, and set $\gamma =G{\rm Y}$. Distinguish spatial radius $r$ from projected radius $R$. The enclosed luminosity and radial force areHere ${\sigma }_t^2$ is the velocity dispersion of either one of the two tangential components; thus $\beta =1-{\sigma }_t^2/{\sigma }_r^2$. If a combined tangential variance is used instead, its contribution must be divided by two in this definition.

To derive the Jeans equation, start with the steady Collisionless Boltzmann equation $v_j{\partial }_{j}f-({\partial }_j\Phi ){\partial }_{v_j}f=0$. Multiply by $v_i$ and integrate over all velocities. Integration by parts in velocity, with $f$ decaying fast enough to eliminate boundary terms, givesThe same equation holds with luminosity weighting because the tracer's luminosity is conserved along its orbit. Spherical symmetry makes the local stellar velocity moment tensor diagonal, with entries $q=\lambda {\sigma }_r^2$, $q_t=\lambda {\sigma }_t^2$, $q_t$. Its radial divergence is $q^{\prime }+(2q-2q_t)/r$: the second term accounts for the changing directions of the radial and tangential unit vectors. ThereforeThis is the Spherical Jeans equation; it follows from a stellar velocity moment rather than requiring a particular form of the galactic distribution function.

Integrating the luminosity density along a line of sight gives its surface brightness:The second identity is Spherical Abel deprojection. For the second stellar velocity moment, the angle between the radial direction and the line of sight has ${\sin }^2\theta =R^2/r^2$. The line-of-sight velocity dispersion is therefore ${\sigma }_r^2{\cos }^2\theta +{\sigma }_t^2{\sin }^2\theta ={\sigma }_r^2(1-\beta R^2/r^2)$, andFrom the Spherical Jeans equation, $\beta q=-r{q}^{\prime }/2-\gamma \lambda L/(2r)$. Substitution producesThe primary PDF has a genuine symbol error in its projected gravity term: it prints a second $\gamma$ multiplying $L$ inside the integral, where $\lambda$ is required. The corrected projected spherical Jeans residual isTo justify the derivative at the singular-looking lower limit, write $A(R)={\int }_0^{\infty }q(\sqrt{R^2+z^2})dz$. Differentiation then gives $A^{\prime }(R)=R{\int }_R^{\infty }{q}^{\prime }(r)(r^2-R^2{)}^{-1/2}dr$. In particular, $Rp=(R^{2}A{)}^{\prime }$. If $R^{2}A$ tends to zero both at zero and infinity, as for a finite isolated equilibrium with sufficiently decaying stresses,If a boundary stress survives, its contribution must be retained; the zero-integral statement is not a finite-aperture identity.

Now identify each integral with a term in the virial theorem. Let $K$ be the total kinetic energy, including all three random components, and $W$ the gravitational potential energy. Averaging over spatial directions in a spherical system gives one third of the total stellar velocity moment along any chosen line of sight, even if the velocities are locally anisotropic. ThusBuilding the spherical mass shell by shell givesFinally, Fubini theorem and ${\int }_0^r{R}^3(r^2-R^2{)}^{-1/2}dR=2r^3/3$ give the gravitational contribution to the residual integral explicitly:ConsequentlyThe global projected virial estimate of a mass-to-light ratio is independent of the velocity-anisotropy parameter under these assumptions. It does not remove the mass--anisotropy--density degeneracy for arbitrary radial mass profiles. In particular, if an independent dark matter halo makes mass fail to follow light, a single constant ${\rm Y}$ is no longer the physical enclosed mass-to-light ratio at every radius.