Made-to-measure N-body systems
Abstract
We describe an algorithm for constructing N-body realisations of equilibrium stellar systems. The algorithm complements existing orbit-based modelling techniques using linear programming or other optimization algorithms. The equilibria are constructed by integrating an N-body system while slowly adjusting the masses of the particles until the time-averaged density field and other observables converge to a prescribed value. The procedure can be arranged to maximise a linear combination of the entropy of the system and the statistic for the observables. The equilibria so produced may be useful as initial conditions for N-body simulations or for modelling observations of individual galaxies.
Cite
@article{arxiv.astro-ph/9605061,
title = {Made-to-measure N-body systems},
author = {D. Syer and S. Tremaine},
journal= {arXiv preprint arXiv:astro-ph/9605061},
year = {2015}
}
Comments
8 pages, tex, figures included, accepted by MNRAS, also available at http://www.mpa-garching.mpg.de/~syer/papers/