Two-integral Schwarzschild models

Abstract

We describe a practical method for constructing axisymmetric two-integral galaxy models (with distribution functions of the form f(E,L_z), in which E is the orbital energy, and L_z is the vertical component of the angular momentum), based on Schwarzschild's orbit superposition method. Other f(E,L_z)-methods are mostly based on solving the Jeans equations or on finding the distribution function directly from the density, which often places restrictions on the shape of the galaxy. Here, no assumptions are made and any axisymmetric density distribution is possible. The observables are calculated (semi-)analytically, so that our method is faster than most previous, fully numerical implementations. Various aspects are tested extensively, the results of which apply directly to three-integral Schwarzschild methods. We show that a given distribution function can be reproduced with high accuracy and investigate the behaviour of the parameter that is used to measure the goodness-of-fit. Furthermore, we show that the method correctly identifies the range of cusp clopes for which axisymmetric two-integral models with a central black hole do not exist.

10 pages, 9 figures, Accepted for publication in MNRAS

Cited by in corpus (24)

• Axisymmetric Dynamical Models of the Central Regions of Galaxies
• TDCOSMO IV: Hierarchical time-delay cosmography -- joint inference of the Hubble constant and galaxy density profiles
• Triaxial orbit based galaxy models with an application to the (apparent) decoupled core galaxy NGC 4365
• The dynamical distance and intrinsic structure of the globular cluster omega Centauri
• Two-dimensional kinematics of SLACS lenses: III. Mass structure and dynamics of early-type lens galaxies beyond z ~ 0.1
• The counterrotating core and the black hole mass of IC1459
• Difficulties with Recovering The Masses of Supermassive Black Holes from Stellar Kinematical Data
• A SAURON study of M32: measuring the intrinsic flattening and the central black hole mass
• Mapping stationary axisymmetric phase-space distribution functions by orbit libraries
• Two-dimensional kinematics of SLACS lenses: II. Combined lensing and dynamics analysis of early-type galaxies at z = 0.08 - 0.33
• Dynamical modelling of stars and gas in NGC2974: determination of mass-to-light ratio, inclination and orbital structure by Schwarzschild's method
• Galaxy Kinematics with VIRUS-P: The Dark Matter Halo of M87
• Regularized orbit models unveiling the stellar structure and dark matter halo of the Coma elliptical NGC 4807
• Two-dimensional kinematics of SLACS lenses: I. Phase-space analysis of the early-type galaxy SDSS J2321-097 at z=0.1
• Recovery of the internal orbital structure of galaxies
• A Unifying Framework for Self-consistent Gravitational Lensing and Stellar Dynamics Analyses of Early-Type Galaxies
• The Fornax3D project: Assembly histories of lenticular galaxies from a combined dynamical and population orbital analysis
• Constraining the Mass Profiles of Stellar Systems: Schwarzschild Modeling of Discrete Velocity Datasets
• On the reliability of the black hole mass and mass to light ratio determinations with Schwarzschild models
• Spiral arm kinematics for Milky Way stellar populations
• Crash-testing the CAULDRON code for joint lensing and dynamics analysis of early-type galaxies
• The non-evolving internal structure of early-type galaxies: the case study SDSS J0728+3835 at z = 0.206
• Accuracy and precision of triaxial orbit models I: SMBH mass, stellar mass and dark-matter halo
• Order and chaos in a new 3D dynamical model describing motion in non axially symmetric galaxies