Basis function expansions for galactic dynamics: Spherical versus cylindrical coordinates

Abstract

Aims. The orbital structure of galaxies is strongly influenced by the accuracy of the force calculation during orbit integration. We explore the accuracy of force calculations for two expansion methods and determine which one is preferable for orbit integration. Methods. We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations.We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations. Results. For the generalized Dehnen density, the CylSP method is more accurate than the HO method for nearly all inner power-law indices and shapes at all radii. For $N$-body realizations of the Dehnen model, or snapshots of an $N$-body simulation, the CylSP method is more accurate than the HO method in the central region for the oblate, prolate, and triaxial Hernquist profiles if the particle number is more than $5\times 10^5$. For snapshots of the Hernquist models with spherical shape, the HO method is preferred. (Abridged)

14 pages, 17 figures, Accepted for publication in A&A

References in corpus (16)

• A Classical Morphological Analysis of Galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S4G)
• Triaxial orbit based galaxy models with an application to the (apparent) decoupled core galaxy NGC 4365
• NMAGIC: Fast Parallel Implementation of a Chi-Squared-Made-to-Measure Algorithm for Modelling Observational Data
• Dearth of dark matter or massive dark halo? Mass-shape-anisotropy degeneracies revealed by NMAGIC dynamical models of the elliptical galaxy NGC 3379
• Orbital decomposition of CALIFA spiral galaxies
• Predicted Space Motions for Hypervelocity and Runaway Stars: Proper Motions and Radial Velocities for the GAIA Era
• The Next Generation Virgo cluster Survey. V. modelling the dynamics of M87 with the Made-to-Measure method
• Evaluating the ability of triaxial Schwarzschild modelling to estimate properties of galaxies from the Illustris simulation
• Applying Schwarzschild's orbit superposition method to barred or non-barred disc galaxies
• Baryonic impact on the dark matter orbital properties of Milky Way-sized haloes
• Expansion Techniques for Collisionless Stellar Dynamical Simulations
• N-body simulations for testing the stability of triaxial galaxies in MOND
• Self-consistent models of triaxial galaxies in MOND gravity
• New biorthogonal potential--density basis functions
• Lopsidedness of self-consistent galaxies by the external field effect of clusters
• Made-to-measure galaxy modelling utilising absorption line strength data