An Adaptive Algorithm for NBody Field Expansions
Abstract
An expansion of a density field or particle distribution in basis functions that solve the Poisson equation both provides an easily parallelized Nbody force algorithm and simplifies perturbation theories. The expansion converges quickly and provides the highest computational advantage if the lowest order potentialdensity pair in the basis looks like the unperturbed galaxy or stellar system. Unfortunately, there are only a handful of such bases in the literature that limit this advantage. This paper presents an algorithm for deriving these bases to match a wide variety of galaxy models. The method is based on efficient numerical solution of the SturmLiouville equation and can be used for any geometry with a separable Laplacian. Two cases are described in detail. First, for the spherical case, the lowest order basis function pair may be chosen to be exactly that of the underlying model. The profile may be (1) cuspy or have a core and (2) truncated or of infinite extent. Second, the method yields a threedimensional cylindrical basis appropriate for studying galactic disks. In this case, the vertical and radial bases are coupled; the lowest order radial part of the basis function can be chosen to match the underlying profile only in the disk plane. Practically, this basis is still a very good match to the overall disk profile and converges in a small number of terms. The ease of combining several bases makes this force solver ideally suited to multicomponent simulations, such as those of disks embedded in halos.
 Publication:

The Astronomical Journal
 Pub Date:
 January 1999
 DOI:
 10.1086/300669
 arXiv:
 arXiv:astroph/9805357
 Bibcode:
 1999AJ....117..629W
 Keywords:

 CELESTIAL MECHANICS;
 STELLAR DYNAMICS;
 GALAXIES: STRUCTURE;
 GALAXY: STRUCTURE;
 METHODS: NUMERICAL;
 Astrophysics
 EPrint:
 16 pages, 5 figures, submitted to AJ