A perturbation particle method for stability studies of stellar systems
Abstract
We discuss a numerical method for investigating the stability of stellar systems that uses an analytic distribution function to describe a stellar system in equilibrium and 'perturbation particles' to represent departures from that equilibrium state. The particles are used only to represent the perturbation, and statistical fluctuations due to the finite number of particles are therefore much less severe than in full Nbody codes. We provide a general description of the method, recipes for particular aspects of its implementation, and an example of its application to a simple model with known analytical solution.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 June 1993
 DOI:
 10.1093/mnras/262.4.1007
 Bibcode:
 1993MNRAS.262.1007W
 Keywords:

 Computational Astrophysics;
 Distribution Functions;
 Perturbation Theory;
 Stellar Motions;
 Stellar Systems;
 Many Body Problem;
 Particle Motion;
 Systems Stability;
 Astrophysics