Numerical stability of a family of Osipkov-Merrit models

Jul, 1997
12 pages
Published in:
  • Astrophys.J. 490 (1997) 136-142
e-Print:
Report number:
  • DF-9701

Citations per year

199820052012201920241026
Abstract: (arXiv)
We have investigated the stability of a set of non-rotating anisotropic spherical models with a phase-space distribution function of the Osipkov-Merritt type. The velocity distribution in these models is isotropic near the center and becomes radially anisotropic at large radii. They are special members of the family studied by Dehnen and Tremaine et al. where the mass density has a power-law cusp ρrγ\rho\propto r^{-\gamma} at small radii and decays as ρr4\rho\propto r^{-4} at large radii. The radial-orbit instability of models with γ\gamma = 0, 1/2, 1, 3/2, and 2, was studied using an N-body code written by one of us and based on the `self-consistent field' method developed by Hernquist and Ostriker. These simulations have allowed us to delineate a boundary in the (γ,ra)(\gamma,r_{a})-plane that separates the stable from the unstable models. This boundary is given by 2Tr/Tt=2.31±0.272T_{r}/T_{t} = 2.31 \pm 0.27, for the ratio of the total radial to tangential kinetic energy. We also found that the stability criterion df/dQ0df/dQ\le 0, recently raised by Hjorth, gives lower values compared with our numerical results.
Note:
  • AASTEX, 22 pages, 11 figures, Figs. 5 available from author. Accepted for publication in Astrophysical Journal Report-no: DF-9701 Journal-ref: Astrophys.J. 490 (1997) 136-142