Stability of rotating spherical stellar systems

The stability of rotating isotropic spherical stellar systems is investigated by using N-body simulations. Four spherical models with realistic density profiles are studied: one of them fits the luminosity profile of globular clusters, while the remaining three models provide good approximations to the surface brightness of elliptical galaxies. The phase-space distribution function f(E) of each one of these non-rotating models satisfies the sufficient condition for stability df/dE < 0. Different amounts of rotation are introduced in these models by changing the sign of the z-component of the angular momentum for a given fraction of the particles. Numerical simulations show that all these rotating models are stable to both radial and non-radial perturbations, irrespective of their degree of rotation. These results suggest that rotating isotropic spherical models with realistic density profiles might generally be stable. Furthermore, they show that spherical stellar systems can rotate very rapidly without becoming oblate.