Physical Sciences and Mathematics Commons^™

An Adaptive Algorithm For N-Body Field Expansions, Md Weinberg

Astronomy Department Faculty Publication Series

An expansion of a density field or particle distribution in basis functions that solve the Poisson equation both provides an easily parallelized N-body force algorithm and simplifies perturbation theories. The expansion converges quickly and provides the highest computational advantage if the lowest order potential-density 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 Sturm-Liouville …