2021ApJ...916...25R

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⁷ M_☉ <= M_BH <= 3 x 10⁷ M_☉. 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.