Publ. Astron. Soc. Japan, 35, 547-558 (1983)

Integration of the Collisionless Boltzmann Equation for Spherical Stellar Systems

Department of Astronomy, University of Kyoto, Sakyo-ku, Kyoto 606
(Received 1983 January 21; accepted 1983 July 8)


The splitting scheme for numerical integration of the collisionless Boltzmann equation, appropriate to spherical stellar systems is described. Both an equilibrium model (the Plummer model) and a nonequilibrium model (a uniform-density sphere) are tested. The results presented show that the scheme is stable and accurate. The evolution of the distribution function for the collapse of a uniform-density sphere has been reproduced in detail. It is found that in the core the velocity distribution after violent relaxation is not Maxwellian but partially degenerate.

Key words: Numerical experiments; Stellar dynamics; Violent relaxation.