Codex @codex  0

At fixed radius, use spherical coordinates in velocity space with polar angle $\alpha$ measured from the radial direction:The galactic distribution function and volume element give angular weightAll dependence on $v$, $r$, and $g(E)$ cancels from ratios of second moments. Symmetry in the tangential plane givesWriting the angular integrals as beta function integrals and using the Gamma function recurrence,Therefore, for every admissible energy factor $g$,Equivalently, the velocity-anisotropy parameter is the constant $\beta =1-p/2$.

At fixed energy, the circular orbit has $L_c^2=GMa$, whileAn isotropic galactic distribution function has conditional angular-momentum density $n(L)dL\propto LdL$. ThereforeNormalization on $0\le e<1$ gives the thermal eccentricity distribution

With the convention in the question, the constant velocity-anisotropy parameter is $\beta =1-{\sigma }_t^2/{\sigma }_r^2$. Since $V_c^2=rd\Phi /dr$, the Spherical Jeans equation becomesFor $\nu ={\nu }_0{r}^{-\alpha }$ and $V_c^2=V_0^2{r}^{2\gamma }$, its general integrating factor solution isThe homogeneous term represents a boundary pressure. For an extended scale-free system the physical boundary condition normally removes it, leavingPositivity requires $\alpha -2\beta -2\gamma >0$. A complete physical model must also have a nonnegative galactic distribution function and sensible inner and outer boundary behaviour; for example, strong radial anisotropy is restricted by density-slope--anisotropy inequalities.

Observationally, the tracer density can be estimated from star counts only after correcting distances, extinction, survey selection, and incompleteness. Spectroscopy supplies mainly line-of-sight velocities; proper motions add transverse information but become less precise for distant halo stars. The equation shows the mass--anisotropy--density degeneracy directly: the same measured ${\sigma }_r$ can result from a larger $V_c$, a steeper tracer slope $\alpha$, or a different $\beta$. Even globally constant power laws therefore do not determine the galactic mass profile unless some of these quantities are independently constrained.

If $\alpha$, $\beta$, or $\gamma$ changes near a break radius, the solution at one radius also depends on the outer boundary integral. A break in observed dispersion may be attributed to a mass-profile feature, a tracer-density break, or a change in orbital anisotropy. Separate tracer populations, full three-dimensional velocities, higher velocity moments, and measurements over a wide radial range help break this degeneracy.

More flexible alternatives model a nonnegative solution of the Collisionless Boltzmann equation itself. An action-based galactic distribution function gives an analytic or parametrized $f(\mathbf{J})$; a Schwarzschild orbit-superposition model assigns nonnegative weights to an orbit library; and a made-to-measure stellar-dynamical model adjusts particle weights to reproduce observations. These methods retain more phase-space information than Jeans moments, although their flexibility introduces model choices and regularization.

A galactic distribution function $f(\mathbf{x},\mathbf{v},t)$ is the stellar mass or number per six-dimensional phase space volume,Its velocity moments give the spatial density, mean velocity, and velocity dispersion; integrating those quantities along the line of sight and weighting by luminosity produces surface-brightness and line-of-sight-velocity observables. A model is compared with data only after the same projection, selection function, and instrumental convolution have been applied.

Jeans theorem states that every steady solution of the Collisionless Boltzmann equation depends on phase-space coordinates only through integrals of motion. Conversely, every nonnegative function of isolating integrals is a steady collisionless distribution function on the region where those integrals are defined.

For the stated power law in relative energy, isotropy givesWith the requested change of variables $v=\sqrt{2\Psi }\cos \theta$, the density becomeswhere the last equality uses the Beta function and Gamma function. Thusfor $n>1/2$.

The normalized second velocity moment isThe same substitution and the Beta-function recurrence giveConsequently the one-dimensional isotropic velocity dispersion is ${\sigma }^2=\overline{v^2}/3=\Psi /(n+1)$, proving the required linear dependence on the relative potential.