Distribution functions for spheroids

Galaxy models comprising several components (including dark matter) that are bound by the self-consistently generated gravitational field are readily constructed from distribution functions (DFs) that are analytic functions of the action integrals J. We explain why such models have unphysical velocity distributions unless the DFs of hot components satisfy certain conditions as J_\phi -> 0. We show how DFs for both isotropic and radially biased spherical systems can be constructed with specified f(J). We show how to construct DFs for flattened systems with significant velocity anisotropy. Construction of self-consistent models rather than populations that are confined by an external potential leads to the conclusion that radially-biased spherical systems are generically unstable to quadrupolar perturbations. Chaos is likely key to maintenance of these constraints during adiabatic disc growth.