Gaussian angular-momentum suppression in a stellar distribution ($f\propto e^{E/{\sigma }^2-L^2/(2r_a^2{\sigma }^2)}$)

In a spherical Newtonian gravitational potential, use binding energy $E=\psi -v^2/2$ and specific angular momentum $L$. The displayed stationary galactic distribution function has locally Gaussian radial and transverse velocities: ${\sigma }_r^2={\sigma }^2$, ${\sigma }_{\theta }^2={\sigma }_{\varphi }^2={\sigma }^2/(1+r^2/r_a^2)$. Its mass density is $\rho ={\rho }_c{e}^{[\psi (r)-\psi (0)]/{\sigma }^2}/(1+r^2/r_a^2)$. The velocity-anisotropy parameter is ${\beta }_{\mathrm{an}}=r^2/(r^2+r_a^2)$, so the core is nearly isotropic and the halo nearly radial. These exact Gaussian integrals use the untruncated distribution over all velocities; imposing an escape-energy cutoff changes them.

For a self-consistent Gaussian angular-momentum suppression in a stellar distribution, put $y=[\psi (r)-\psi (0)]/{\sigma }^2$ and $s=\log (r/r_a)$. At large radius the Poisson equation for Newtonian gravity reduces to $y_{ss}+y_s\sim -\lambda e^y$, where $\lambda =4\pi G{\rho }_c{r}_a^2/{\sigma }^2$. The asymptotic balance $y_s\sim -1/s$ gives $e^y\sim 1/(\lambda s)$, the displayed density, and circular speed $v_c^2\sim {\sigma }^2/\log (r/r_{\ast })$. A sufficiently small regular core can have an intermediate isotropic singular isothermal sphere envelope with $\rho \simeq {\sigma }^2/(2\pi G{r}^2)$ and $v_c^2\simeq 2{\sigma }^2$. At the actual regular centre, however, $\rho \to {\rho }_c$ and $v_c\propto r$; the intermediate power law must not be mistaken for a central limit.