Abstract
The mass of a supermassive black hole (M (BH)) is a fundamental property that can be obtained through observational methods. Constraining M (BH) through multiple methods for an individual galaxy is important for verifying the accuracy of different techniques and for investigating the assumptions inherent in each method. However, there exist only a few galaxies where multiple M (BH) measurement techniques can be applied. NGC 4151 is one of these rare galaxies for which multiple methods can be used: stellar and gas dynamical modeling because of its proximity (D = 15.8 +/- 0.4 Mpc from Cepheids), and reverberation mapping because of its active accretion. In this work, we reanalyzed H-band integral field spectroscopy of the nucleus of NGC 4151 from Gemini NIFS, improving the analysis at several key steps. We then constructed a wide range of axisymmetric dynamical models with the new orbit-superposition code Forstand. One of our primary goals is to quantify the systematic uncertainties in M (BH) arising from different combinations of the deprojected density profile, inclination, intrinsic flattening, and mass-to-light ratio. As a consequence of uncertainties on the stellar luminosity profile arising from the presence of the active galactic nucleus, our constraints on M (BH) are rather weak. Models with a steep central cusp are consistent with no black hole; however, in models with more moderate cusps, the black hole mass lies within the range of 0.25 x 10(7) M (circle dot) less than or similar to M (BH) less than or similar to 3 x 10(7) M (circle dot). This measurement is somewhat smaller than the earlier analysis presented by Onken et al. but agrees with previous M (BH) values from gas dynamical modeling and reverberation mapping. Future dynamical modeling of reverberation data, as well as IFU observations with JWST, will aid in further constraining the mass of the central supermassive black hole in NGC 4151.