An iterative method for the construction of Nbody galaxy models in collisionless equilibrium
Abstract
We describe a new iterative approach for the realization of equilibrium Nbody systems for given density distributions. Our method uses elements of Schwarzschild's technique and of the madetomeasure method, but is based on a different principle. Starting with some initial assignment of particle velocities, the difference of the timeaveraged density response produced by the particle orbits with respect to the initial density configuration is characterized through a merit function, and a stationary solution of the collisionless Boltzmann equation is found by minimizing this merit function directly by iteratively adjusting the initial velocities. Because the distribution function is in general not unique for a given density structure, we augment the merit function with additional constraints that single out a desired target solution. The velocity adjustment is carried out with a stochastic process in which new velocities are randomly drawn from an approximate solution of the distribution function, but are kept only when they improve the fit. Our method converges rapidly and is flexible enough to allow the construction of solutions with third integrals of motion, including disc galaxies in which radial and vertical dispersions are different. A parallel code for the calculation of compound galaxy models with this new method is made publicly available.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 October 2014
 DOI:
 10.1093/mnras/stu1421
 arXiv:
 arXiv:1402.1623
 Bibcode:
 2014MNRAS.444...62Y
 Keywords:

 methods: numerical;
 stars: kinematics and dynamics;
 galaxies: haloes;
 galaxies: kinematics and dynamics;
 galaxies: structure;
 Astrophysics  Cosmology and Nongalactic Astrophysics
 EPrint:
 18 pages, 15 figures, 1 table. Accepted by MNRAS. Minor changes from submitted version. The source code of GALIC is publicly available now and can be downloaded here: http://www.hits.org/tap/galic/