Logarithmic-potential two-integral distribution ($F=A{L}_z^2{e}^{-4E/V^2}+B{e}^{-2E/V^2}$)

The singular axisymmetric logarithmic gravitational potential admits this stationary galactic distribution function withThe zero of specific orbital energy is the printed logarithmic convention. Products of Gaussian integrals normalize its two spatial density terms. For $1/2\le q^2\le 1$ the two terms are separately nonnegative; mildly prolate models require a bound on their sum.

For the logarithmic-potential two-integral distribution, $q^2>1$ makes $A<0$. Accessible phase space satisfies $L_z^2{e}^{-2E/V^2}\le V^2/e$, attained by equatorial circular motion. Hence the galactic distribution function is nonnegative precisely when $B+A{V}^2/e\ge 0$, giving the displayed upper bound. A larger prolate flattening parameter produces negative phase-space density even though the spatial mass density remains positive.