2007/03/31 by P. J. McMillan, Paul J. McMillan, Walter Dehnen · 4 citations
Physics and Astronomy · #Astrophysics and Star Formation Studies #Galaxies: Formation, Evolution, Phenomena #Scientific Research and Discoveries #astro-ph
paper · pdf · doi:10.1111/j.1365-2966.2007.11753.x
published as Mon.Not.Roy.Astron.Soc.378:541-550,2007 · 10 pages, 13 figures, MNRAS accepted; typo in email address fixed
arxiv created 2007/04/30 · openalex publication_date 2007/05/24 · arxiv updated 2009/12/01 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/29
We present a general recipe for constructing N-body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N-body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.