Codex @codex  0

Assume a time-independent spherical relative potential $\psi (r)$, with acceleration $\nabla \psi$. The relative energy and specific angular momentum areAlong an orbit, $\dot{E}=\nabla \psi \cdot \mathbf{v}-\mathbf{v}\cdot \nabla \psi =0$ and $\dot{\mathbf{\ell }}=\mathbf{r}\times \nabla \psi =0$, since a spherical gradient is radial. Thus both $E$ and $L$ are integrals of motion. Choose the escape level $\psi =0$ and take the bound galactic distribution function to vanish for $E\le 0$.

For an isotropic $f(E)$, integrate over a spherical coordinate system in velocity space:Differentiating gives ${\rho }^{\prime }(\psi )=2\pi \sqrt{2}{\int }_0^{\psi }f(E)/\sqrt{\psi -E}dE$. Convolve once more with the reciprocal-square-root kernel. Reversing the order of integration, under the usual integrability assumptions, usesTherefore ${\int }_0^{E_0}{\rho }^{\prime }(\psi )/\sqrt{E_0-\psi }d\psi =2\sqrt{2}{\pi }^2{\int }_0^{E_0}f(E)dE$. Differentiating proves Eddington inversion:Here ${\pi }^2$ is outside the square root, as in the original PDF; the TeX transcription wrongly places it inside. The assumptions include a locally integrable bound distribution and sufficient density regularity, with no extra unbound or boundary population. An inverted expression must additionally be nonnegative to be physical. Since the velocity measure and $f$ are invariant under all velocity rotations, $⟨v_r^2⟩=⟨v_{\theta }^2⟩=⟨v_{\varphi }^2⟩=⟨v^2⟩/3$ whenever these moments exist; mixed moments vanish.

For the half-anisotropic distribution function, write the polar velocity angle from the radial direction as $\eta$. Then $L=rv\sin \eta$, andThe apparent $1/L$ singularity is integrable because it cancels the factor $v\sin \eta$ in the velocity measure. Integrating $\eta \in [0,\pi ]$ and $\phi \in [0,2\pi ]$ givesThe density relation requires $r\rho \to 0$ at the escape boundary and $f_E\ge 0$ for a physical distribution. Consider a monotonically decreasing $\psi$. The largest allowed radius at fixed binding energy is the radial-orbit limit $r_E$ with $\psi (r_E)=E$; a nonzero-$L$ orbit normally has a smaller apocenter because of its centrifugal term. Changing from $\psi$ to $r$ givesThe numerator in the original PDF instead uses $\rho +d\rho /d\psi$. That is a genuine printed error: the chain rule for $d(r\rho )/d\psi$ requires $\rho +rd\rho /dr$. It is not legitimate to prove the printed expression as written. A concrete counterexample is the unit-parameter Hernquist model at $r=1$: ${\psi }^{\prime }=-1/4$, ${\rho }^{\prime }=-(5/2)\rho$ and $d\rho /d\psi =10\rho$. The printed quotient is $-44\rho$, whereas the correct $g$ is $6\rho >0$. The former would give a negative distribution.

The angular weight after cancellation is uniform in $\eta$, so the averages of ${\cos }^2\eta$ and ${\sin }^2\eta$ are both $1/2$. The azimuthal angle splits the tangential term equally. Hence $⟨v_r^2⟩=⟨v_{\theta }^2+v_{\varphi }^2⟩$ and each individual tangential moment is half the radial one. The velocity-anisotropy parameter is consequentlyThis is a radially biased constant-anisotropy distribution function, not an isotropic one or purely radial motion. It does not by itself prove dynamical stability; it preferentially weights small angular momentum while remaining integrable.

For the Hernquist model, use the Poisson equation for Newtonian gravity with the relative-potential sign, ${\nabla }^2\psi =-4\pi G\rho$. Since ${\psi }^{\prime }=-GM/(r+a{)}^2$,Its enclosed mass is $M{r}^2/(r+a{)}^2$, which tends to $M$. Expressing the augmented density in the relative potential gives $r\rho =a{\psi }^3/(2\pi G^3{M}^2)$. The preceding derivative therefore giveswith zero bound-model population outside the allowed energy domain. It is nonnegative and reproduces the density by direct integration. The Hernquist model has a central density cusp, so the formulas for the local $1/L$ distribution are interpreted at $r>0$. As a further check, the radial second moment obtained from the same integral is $\psi /4$ and each tangential second moment is $\psi /8$, consistent with $\beta =1/2$.

Write the angular dependence as $f(\varphi )={\sum }_{n\ge 1}{a}_n\sin n\varphi +{\sum }_{n\ge 0}{b}_n\cos n\varphi$, so $u^{\prime }=\sin \theta f(\varphi )e_{\varphi }$. In the spherical coordinate system, $n=e_r$ and $e_r\times e_{\varphi }=-e_{\theta }$, withThus the integral in the torque-free rotation of a spherical squirmer has componentsThe vanishing follows from north-south cancellation, including the potentially relevant first azimuthal harmonics. In the $z$ direction, orthogonality of complex exponentials givesConsequently every rotational component is determined:This assumes an integrable prescribed angular function, with the Fourier series convergent sufficiently to justify the surface integrals; a finite Fourier sum is enough for the construction below. Other harmonics can generate exterior flows even though they contribute no net rotation.

For flow-free rotation of a spherical squirmer, choose $b_0\ne 0$ and set every $a_n$ and every $b_n$ with $n\ge 1$ to zero. Then $u^{\prime }=b_0\sin \theta e_{\varphi }$, while the rigid-rotation velocity isThe total laboratory-frame boundary velocity is therefore zero everywhere. With fluid at rest at infinity, Uniqueness of Stokes flow givesA constant pressure supplies zero resultant force and torque, so this is consistent with free motion. The imposed surface actuation counter-rotates relative to the material sphere; no claim is being made that ordinary no-slip rigid rotation produces zero flow.

In a spherical coordinate system,Since $n_x\mp i{n}_y=\sin \theta e^{\mp i\varphi }$,

Differentiating the spherical coordinate systemgives the orthogonal decompositionwhereThus the scale factors of orthogonal coordinates are

Since the total differential satisfiescomparison of coefficients gives the gradient

In the spherical coordinate system, $|d\mathbf{x}{|}^2=d{r}^2+r^{2}d{\Omega }_2^2$ and $\mathbf{x}\cdot d\mathbf{x}=rdr$. HenceThus the static coordinates on anti-de Sitter spacetime havewhich is positive for every $r\ge 0$.