Journal Article

Multimass schemes for collisionless <i>N</i>-body simulations

Mimi Zhang and John Magorrian

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 387, issue 4, pages 1719-1726
Published in print July 2008 | ISSN: 0035-8711
Published online July 2008 | e-ISSN: 1365-2966 | DOI: http://dx.doi.org/10.1111/j.1365-2966.2008.13362.x
Multimass schemes for collisionless N-body simulations

Show Summary Details

Preview

We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f0. Our method uses importance sampling to find the sampling DF fs that minimizes the mean-square formal errors in a given set of projections of the DF f0. The result is a multimass N-body realization of the galaxy model in which ‘interesting’ regions of phase space are densely populated by lots of low-mass particles, increasing the effective N there, and less interesting regions by fewer, higher mass particles.

As a simple application, we consider the case of minimizing the shot noise in estimates of the acceleration field for an N-body model of a spherical Hernquist model. Models constructed using our scheme easily yield a factor of ∼100 reduction in the variance at the central acceleration field when compared to a traditional equal-mass model with the same number of particles. When evolving both models with a real N-body code, the diffusion coefficients in our model are reduced by a similar factor. Therefore, for certain types of problems, our scheme is a practical method for reducing the two-body relaxation effects, thereby bringing the N-body simulations closer to the collisionless ideal.

Keywords: methods: N-body simulations; methods: numerical; galaxies: kinematics and dynamics

Journal Article.  5257 words.  Illustrated.

Subjects: Astronomy and Astrophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.