Astrophys. J., 602, 66-92 (2004/February-2)
Difficulties with recovering the masses of supermassive black holes from stellar kinematical data.
VALLURI M., MERRITT D. and EMSELLEM E.
Abstract (from CDS):
We investigate the ability of three-integral, axisymmetric, orbit-based modeling algorithms to recover the parameters defining the gravitational potential (mass-to-light ratio Υ and black hole mass Mø) in spheroidal stellar systems using stellar kinematical data. We show that the potential estimation problem is generically underdetermined when applied to long-slit kinematical data of the kind used for most black hole mass determinations to date. A range of parameters (Υ, Mø) can provide equally good fits to the data, making it impossible to assign best-fit values. The indeterminacy arises from the large variety of orbital solutions that are consistent with a given mass model. We demonstrate the indeterminacy using a variety of data sets derived from realistic models, as well as published observations of the galaxy M32. The indeterminacy becomes apparent only when a sufficiently large number of distinct orbits are supplied to the modeling algorithm; if too few orbits are used, spurious minima appear in the χ2(Υ,Mø) contours, and these minima do not necessarily coincide with the parameters defining the gravitational potential. We show that the range of degeneracy in Mødepends on the degree to which the data resolve the radius of influence rhof the black hole. For FWHM/2rh≳0.5, where FWHM refers to the instrumental resolution, we find that only very weak constraints can be placed on Mø. In the case of M32, our reanalysis demonstrates that when a large orbit library is used, data published prior to 2000 (FWHM/2rh~0.25) are equally consistent with black hole masses in the range 1.5x106M☉<Mø<5x106 M☉, with no preferred value in that range. Exactly the same data can reproduce previous published results with smaller orbit libraries. While the Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) data for this galaxy (FWHM/2rh~0.06) may overcome the degeneracy in Mø, HST data for most galaxies do not resolve the black hole's sphere of influence and in these galaxies the degree of degeneracy allowed by the data may be greater than previously believed. We investigate the effect of regularization, or smoothness constraints, on the degree of degeneracy of the solutions. Enforcing smoothness reduces the range of acceptable models, but we find no indication that the true potential can be recovered simply by enforcing smoothing. For a given smoothing level, all solutions in the minimum χ2 valley exhibit similar levels of noise; as the smoothing is increased, there is a systematic shift in the midpoint of the χ2 valley, until at a high level of smoothing the solution is biased with respect to the true solution. These experiments suggest both that the indeterminacy is real (i.e., that it is not an artifact associated with nonsmooth solutions) and that there is no known empirical way to choose the smoothing parameter to ensure that the correct solution is selected.
Galaxies: Elliptical and Lenticular, cD - Galaxies: Nuclei - Galaxies: Structure - Stellar Dynamics
View the reference in ADS
To bookmark this query, right click on this link: simbad:2004ApJ...602...66V and select 'bookmark this link' or equivalent in the popup menu