Dynamical Modeling of Velocity Profiles: The Dark Halo around the Elliptical Galaxy NGC 2434
Abstract
We describe a powerful technique to model and interpret the stellar lineofsight velocity profiles of galaxies. It is based on Schwarzschild's approach to build fully general dynamical models. A representative library of orbits is calculated in a given potential, and the nonnegative superposition of these orbits is determined that best fits a given set of observational constraints. The most significant new feature of our implementation is that we calculate and fit the full velocity profile shapes, represented by a GaussHermite series. This allows us to constrain the orbital anisotropy in the fit. We also use an objective χ^{2} measure for the quality of fit, taking into account the error on each observational constraint. Given χ^{2} from the observational constraints, the technique assesses the relative likelihood of different orbit combinations in a given potential, and of models with different potentials. In our implementation only projected, observable quantities are included in the fit, aperture binning and seeing convolution of the data are properly taken into account, and smoothness of the models in phase space can be enforced through regularization. This scheme is valid for any geometry.
In a first application of this method, we focus here on spherical geometry; axisymmetric modeling is described in companion papers by Cretton et al. and van der Marel et al. We test the scheme on pseudodata drawn from an isotropic Hernquist model and then apply it to the issue of dark halos around elliptical galaxies. We model radially extended stellar kinematical data for the E0 galaxy NGC 2434 obtained by Carollo et al. This galaxy was chosen because it may be nearly round, in which case the present spherical modeling is applicable. Models with constant masstolight ratio are clearly ruled out, regardless of the orbital anisotropy. To study the amount of dark matter needed to match the data, we considered a sequence of cosmologically motivated ``star + halo'' potentials. These potentials are based on the CDM simulations by Navarro et al., but also account for the accumulation of baryonic matter; they are specified by the stellar masstolight ratio Υ_{*,B} and the characteristic halo velocity, V_{200}. The star + halo models provide an excellent fit to the data, with Υ_{*,B} = 4.35 +/ 0.35 (in Bband solar units) and V_{200} = 450 +/ 100 km s^{1}. The bestfitting potential has a circular velocity V_{c} that is constant to within ~10% between 0.2 and 3 effective radii and is very similar to the bestfitting logarithmic potential, which has V_{c} = 300 +/ 15 km s^{1}. In NGC 2434 roughly half of the mass within an effective radius is dark. In comparison, our models without a dark halo estimate a masstolight ratio for the stellar population that is twice as large. If NGC 2434 is a significantly flattened system seen nearly faceon, it would be considerably more difficult to limit the gravitational potential without further observational constraints.
 Publication:

The Astrophysical Journal
 Pub Date:
 October 1997
 DOI:
 10.1086/304733
 arXiv:
 arXiv:astroph/9702126
 Bibcode:
 1997ApJ...488..702R
 Keywords:

 Galaxies: Elliptical and Lenticular;
 cD;
 Galaxies: Halos;
 Galaxies: Individual: NGC Number: NGC 2434;
 Galaxies: Kinematics and Dynamics;
 Galaxies: Structure;
 Astrophysics
 EPrint:
 41 pages, Latex file with 11 PostScripts figures. Submitted to the Astrophysical Journal