long before the development of general relativity , john michell wondered about the gravitational influence of objects on the light they emit , and how one might go about finding objects so dense that light could not escape their surfaces . `` [ i]f any other luminous bodies should happen to revolve about them we might still perhaps from the motions of these revolving bodies infer the existence of the central ones with some degree of probability , as this might afford a clue to some of the apparent irregularities of the revolving bodies , which would not be easily explicable on any other hypothesis '' @xcite .
more than 200 years later , precisely this method has been employed to determine the mass of the black hole ( bh ) at the center of our milky way .
sgr  a@xmath10 is one of the most tightly constrained bhs in the universe , with a total mass uncertainty ( statistical and systematic ) of less than 10% ( see * ? ? ?
* and references therein ) .
using stellar motions to measure the bh masses in more distant galaxies is complicated by our current inability to spatially resolve the orbits of individual stars .
instead , one must rely on the luminosity - weighted line - of - sight velocity distribution ( losvd ) within different spatial resolution elements , and model the most likely contribution to the gravitational potential from the unseen bh ( beyond that provided by the stars themselves and any dark matter ) .
the most common numerical approach to constraining the bh mass ( @xmath0 ) and mass - to - light ratio of the stars ( @xmath11 ) is the orbit superposition method of @xcite , which simultaneously optimizes the fit to the surface brightness distribution and the line - of - sight kinematics of stars within the nuclear region .
application of that technique has yielded more than 50 bh mass estimates ( e.g. , * ? ? ?
* ; * ? ? ?
a complimentary mass scale has been developed for the bhs in active galactic nuclei ( agns ) , based on the technique of reverberation mapping ( rm ; * ? ? ?
* ; * ? ? ?
measuring the size - scale and velocity of gas near the bh has provided mass estimates for @xmath1250 bhs @xcite , but systematic uncertainties remain due to the unknown geometry and dynamics of the gas .
this uncertainty is traditionally encapsulated in a scaling ( or projection ) factor , @xmath13 an empirical calibration for the rm mass scale @xmath14 has been established by assuming that those agns also lie on the tight correlation between bh mass and bulge stellar velocity dispersion ( the @xmath15 relation ) found in quiescent galaxies .
this was first done by @xcite , who obtained an average rm mass scale of @xmath16 .
subsequent determinations have ranged from @xmath17 @xcite to @xmath18 @xcite , with the most recent value being @xmath19 @xcite . for ngc  4151 , this gives a current rm - derived bh mass of @xmath20  ( based on the revised estimate of * ? ? ?
it is important to note that the error generally quoted on a rm - derived bh mass in an individual galaxy is the formal uncertainty in the virial product and does not include the uncertainty in the nominal mean scale factor @xmath21 .
what has been missing is an independent bh mass measurement for a sample of galaxies with rm estimates .
early results have suggested that the empirically calibrated rm masses are consistent with the values produced by other techniques @xcite , but the few agns near enough to be probed by stellar dynamics have not yielded particularly tight constraints ( in part , owing to the complicating factor of the bright agn overwhelming the stellar absorption features at the smallest galactic radii ) . while only a few galaxies lend themselves to multiple methods for measuring the masses of their central bhs , such independent bh mass determinations are crucial to isolating the systematic errors inherent in each method .
it is only via independent measurements that it is possible to derive robust bh masses .
ngc  4151 is a particularly interesting galaxy since it is one of the rare objects for which there already exist two independent dynamical measurements for the mass of the bh .
@xcite use optical spectroscopy of stellar absorption lines and orbit superposition modeling to derive an upper limit for the bh mass of @xmath22 when the galactic bulge is assumed to be edge - on .
when the bulge is assumed to have the same inclination as the large - scale disk , they derive a bh mass of 4 - 5@xmath23 , which they caution is likely to be a biased solution since it is associated with a very noisy @xmath1 surface that also gives a poor fit to the data ( in terms of absolute @xmath1 value ) .
@xcite use adaptive optics ( ao ) to obtain high spatial resolution , near - infrared spectroscopy of the molecular , ionized , and highly ionized gas in ngc  4151 . from the kinematics of gas in the vicinity of the bh
, they derive a dynamical mass of @xmath24 for the bh .
being somewhat smaller than the value obtained by @xcite , this suggests a need to revisit the stellar dynamical bh mass measurement .    two serious concerns with the result of @xcite
motivate the need for the high spatial resolution observations that are presented in this paper .
first , based on the rm mass and the bh mass derived by @xcite , the sphere of influence of the bh is estimated to be @xmath25 , significantly smaller than the spatial resolution of the spectra used by @xcite .
second , low spatial resolution integral field spectroscopy of a larger field of view @xcite shows that the true kinematic major axis of the bulge ( derived from the 2d velocity field ) is oriented @xmath26 from the axis assumed by @xcite .
the observations and modeling presented in this paper are designed to overcome both of these drawbacks and to provide a more robust comparison with the gas dynamical measurement of @xcite .    to enhance the comparison between the rm and stellar dynamical mass scales for bhs
, we have obtained near - ir integral field spectroscopy with ao in order to study the stellar dynamics in the inner regions of the reverberation - mapped agn , ngc  4151 . in conjunction with our previous optical spectroscopy and multi - wavelength imaging data ,
we are able to place the tightest dynamical constraints to date on the bh mass in ngc  4151 . in section  2
, we describe the observational data . in section  3
, we discuss our stellar dynamical modeling .
section  4 describes our results .
section  5 contains further discussion and our conclusions .
ngc  4151 is a nearby ( @xmath27 ) galaxy with a prominent bulge , faint spiral arms and a large scale bar ( see sections 3.4 and 3.5 for more details on the morphology and distance of the galaxy ) .
the galaxy s nuclear activity has been studied extensively since the original paper by seyfert ( 1943 ) ; in particular , the analysis of uv and optical emission lines have been used for determining the bh mass via rm ( e.g. , * ? ? ?
the analysis presented here uses new near - ir spectroscopy performed with ao , as well as previously published optical long - slit spectroscopy , and space- and ground - based imaging ( primarily from * ? ? ?
* hereafter paper  i ) to provide a mass estimate independent of the rm value
.      observations of ngc  4151 and three giant stars ( used as velocity templates ) were obtained with the near - infrared integral field spectrograph ( nifs ; * ? ? ? * ) on gemini north between 16 feb 2008 and 24 feb 2008 in queue mode ( program i d gn-2008a - q-41 ) .
the nifs image - slicing integral field unit ( ifu ) provides @xmath28 near - ir spectra across a @xmath29 field , with a spatial sampling of @xmath30 , or @xmath31 `` spaxels '' ( spatial pixels ) , where the smaller pixel scale is along each of the 29 image slices .
because previous observations of the k - band stellar absorption lines in ngc  4151 show the features to be weak in the central regions @xcite , we focus on the co  ( 6 - 3 ) bandhead at a rest - frame wavelength of 1.62  @xmath32 m , which has been found to be quite prominent in previous nifs spectra @xcite .
the @xmath2-band spectra cover a wavelength range of @xmath33m@xmath34 m ( centered at 1.6490  @xmath32 m ) , with a dispersion of 1.60    pixel@xmath35 across the 2040 spectral pixels , and a velocity resolution of 56.8  .
to take advantage of the spatial sampling afforded by nifs , the observations were taken with the altair ao system @xcite , using the agn in ngc  4151 as the ao guide star .
the observations employed the altair field lens @xcite to provide a more uniform point spread function ( psf ) across the nifs field .
the level of ao correction depends upon the prevailing seeing conditions at the time .
the median uncorrected seeing was @xmath36 , and more than 95% of the data were obtained in seeing better than @xmath37 , which provided an ao - corrected nuclear fwhm between @xmath38 and @xmath39 . to estimate the psf
, we model the contribution from the agn to the spectrum in each spaxel across the field of view , taking advantage of the agn s strong [ ]  @xmath40 m emission line , at an observed wavelength of @xmath41 m .
although the [ ] emission is spatially extended in ngc  4151 , the off - nuclear emission line is narrower than the nuclear line profile , and tight constraints can be placed on the amount of agn spectrum that ends up in each spaxel ( see   4.2 of * ? ? ?
the analysis was restricted to a wavelength region around the [ ] line ( @xmath42m@xmath43 m ) , and to control the effects of noise , the spectra were box - car smoothed over 5 pixels .
a slice was taken through the center of the resulting 2d map , cutting along the smaller pixel scale ( @xmath44 ) .
the psf profile shows the standard core+halo structure delivered by ao systems , which we model as a 2-part gaussian , where the core has @xmath45 , the halo has @xmath46 , and each component has equal amplitude .
the data for ngc  4151 were acquired at a position angle of @xmath47 , which aligns the nifs slitlets perpendicular to the radio jet ( matching the orientation used by * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the typical observing sequence included a 9-point dither pattern with separations of 0206 , and offsets to a sky field ( 210  away ) approximately every third frame .
individual exposures were 120  s long , and the observations were split into sequences of 80 minutes in duration ( in some instances , shortened by changing weather conditions ) .
arc lamps were taken in the middle of each observing sequence , and telluric standards were obtained before and after each sequence .
the telluric stars were chosen to be 2mass sources to simultaneously provide photometric calibration : hd  98152 was observed before ngc  4151 , and hd  116405 afterwards .
darks , flats , and ronchi mask calibrations were obtained during the daytime after each night s science observations . in total , 253 on - source frames
were obtained , for a total exposure time of 30360  s. a total of 112 sky frames were also acquired .
the giant stars we observed were included so as to provide a measure of the spectral broadening of the instrument , allowing the intrinsically narrow - featured stars to serve as templates that could be broadened to match the actual stellar velocity distribution of ngc  4151 .
these velocity template stars were observed without ao , and were selected to span a range of spectral type : hd  35833 ( g0 ) , hd  40280 ( ) , hip  60145 ( m0 ) .
the observations of each star consist of four on - source and two sky exposures , each with an integration time of 5.3  s.      the nifs data were reduced with a combination of version 1.9.1 of the gemini iraf package and our own iraf scripts .
calibrations were processed independently for each night .
the observations of ngc  4151 were grouped by proximity to a telluric standard , typically dividing the observing sequence in half .
the data were flat - fielded , dark - subtracted , wavelength - calibrated , and spatially rectified , then combined , corrected for telluric features , and flux - calibrated .
the flux calibration provided by the two telluric / flux standards differs systematically by @xmath48% , implying a 0.3  mag discrepancy in the 2mass photometry for the two stars .
because the absolute flux levels in the nifs data play no role in our subsequent analysis , we do not attempt to correct for this offset .
since the image quality of some frames is significantly worse than others , we select a subsample of individual frames with fwhm values below a certain upper limit . with a threshold of 016 in the fwhm ( measured at a uniform wavelength of 1.6598  @xmath32 m to avoid emission features in the spectra ) ,
236 frames are included in the final datacube .
we spatially shift the data to account for the dither pattern and then median - combine the images , yielding a datacube containing 29@xmath4969  spatial pixels and 2040 spectral pixels .
the near - ir spectra at each spatial position in the final nifs datacube are combined onto a bin of size @xmath1202-square ( 2@xmath495 raw spaxels ) , resulting in a cube of 15@xmath4915 spaxels .
this choice of binning the spaxels arose from two main considerations .
first , prior to carrying out the axisymmetric dynamical modeling , the kinematic data need to be symmetrized about the kinematic major axis and reflected about the rotation axis , to avoid modeling inconsistencies .
this process is best accomplished with square bins , especially when the symmetry axes of the galaxy are not aligned with the axes of the spectrograph ( as is generally the case ) .
second , since the halo of the psf has @xmath46 and this component comprises 50% of the overall psf amplitude , the spatial resolution of the data is not significantly altered by using larger bins .
we reduce the observations of the telluric / flux standard stars and the velocity template stars with procedures similar to those we apply to the ngc  4151 data .
we median - combine the four frames for each star , extract spectra with a fixed aperture of 1@xmath505 , and then we normalize each extracted spectrum for the velocity template stars by a low - order fit to the continuum .
in addition to the new data from nifs , we utilize the results from the optical spectra of ngc  4151 that were analyzed in paper  i. the data consist of long - slit spectra of the calcium triplet region ( @xmath51 ) taken with the mayall 4  m telescope at kitt peak national observatory in 2001 , and with the 6.5  m mmt observatory at mt .
hopkins in 2004 . for the data from kitt peak , the spectroscopic slit was 2wide and was placed at a position angle of 135@xmath52 . to enhance the signal - to - noise
, spectra were extracted with an aperture that increased in size away from the galaxy center , extending to radii of @xmath5314 .
the spectral resolution was @xmath54  .
for the observations with the mmt , the slit was 1  wide , at a position angle of 69@xmath52 .
spectra were extracted in @xmath55 bins to a distance of 12on either side of the agn .
the spectral resolution was @xmath56  .
the seeing conditions for the two datasets were @xmath121.8  and @xmath123 , respectively .
we refer the reader to paper  i for full details .
to model the surface brightness profile of ngc  4151 in the near - ir , we make use of the @xmath2-band imaging from the early data release of the ohio state university bright spiral galaxy survey ( osubsgs ; * ? ? ?
the data were obtained with the 1.8  m perkins telescope at lowell observatory .
the image has a plate scale of 15  pixel@xmath35 , and we estimate the seeing to be 2 .
we photometrically calibrate the osubsgs imaging using 2mass point - sources in the field , deriving a photometric zero - point of @xmath57  mag .
to bridge the gap in the surface brightness information between the small field - of - view of the ao - assisted nifs data and the large - scale low resolution imaging data of the osubsgs , we use the _ hubble space telescope _
images from paper  i. these data were taken with the f550 m filter on the high resolution channel of the advanced camera for surveys ( acs ) , and cover a field - of - view of @xmath1225@xmath4929 .
we refer the reader to paper  i for additional details .
we also make use of @xmath58- and @xmath59-band optical imaging from the sloan digital sky survey ( sdss ) to assist in constraining the mass - to - light ratio ( @xmath11 ) .
the images have an exposure time of 53.9  s in each band , with 0396 pixels and @xmath1210-seeing in both filters .
we use the photometrically calibrated and sky - subtracted images that sdss makes available .
as in paper  i , our modeling relies on simultaneously fitting both the luminosity distribution of the bulge and the observed integrated line - of - sight velocity distributions ( losvd ) of stars obtained with various spectrographs over a range of spatial scales .
we detail below the methods of deriving those quantities from the data described above .      as discussed in [ sec : nifs_details ] , the nifs instrument produces a grid of spatial pixels ( `` spaxels '' ) of size @xmath30 which were combined onto a bin of size 02-square , resulting in a cube of 15@xmath4915 spaxels . the spectrum from each binned spaxel
is then analyzed using the `` penalized pixel fitting '' ( ppxf ) package for idl by @xcite .
the ppxf routine computes the velocity shift , the velocity dispersion , and the higher - order moments of the stellar losvds ( parameterized as gauss - hermite ( gh ) polynomials ) that , when convolved with the velocity template star spectrum , provide the best match to the observed data .
the software allows emission features to be masked in the fit , and also provides for multiplicative and additive polynomials to be included ( which we use to account for the contribution of the agn continuum ) .
our ppxf analysis fits the losvd up to 4th - order gh terms , including 2nd - order multiplicative and additive polynomials , and fits the spectra over the wavelength range @xmath60m@xmath61 m , with four small wavelength gaps around agn emission lines and residual sky features .
figure  [ fig.spec ] shows four example spectra and ppxf fits , from different positions within the datacube .
we obtain the smallest uncertainties on the kinematic parameters when we fit the agn spectra with the m0 velocity template , and allowing simultaneous use of all three templates does not measurably improve the results .
we also compared the velocity measurements to those obtained when using nifs velocity templates ( spectral types , , , and a ) acquired by watson et al .
( 2008 ; gemini program gn-2006b - sv-110 ) .
the m - type stars provided better matches to the relative strengths of the absorption lines in each ngc  4151 spectrum , and the kinematic uncertainties were smaller than for the k - type stars .
compared to the kinematics derived from our m0 velocity template , the differences were typically less than 1@xmath4 , suggesting our fits are robust against the effects of template - mismatch .
our results are consistent with those of @xcite , who found better matches to @xmath2-band galaxy spectra with m - type stars than with k - type stars .
figure  [ fig : vel_fields_orig ]  ( top row ) shows maps of the line - of - sight velocity ( @xmath62 ) , velocity dispersion ( @xmath63 ) , and the 3rd and 4th gh moments ( @xmath64 ) in 02 spaxels .
the @xmath65 coordinates are aligned with the major and minor kinematic axes , respectively .
the color scale is in km / s for @xmath62 and @xmath63 , and in dimensionless units for @xmath66 and @xmath67 .
we then bi - symmetrize the 2d losvd map ( fig .
[ fig : vel_fields_orig ]  bottom row ) around the kinematic minor axis at a position angle of @xmath68 ( @xmath69 relative to the nifs y - axis ) , using an idl routine developed by @xcite .
this procedure enhances the s / n by combining the measurements on either side of the kinematic minor axis in a symmetric ( @xmath63 and @xmath67 ) , or anti - symmetric ( @xmath62 and @xmath66 ) way .
data point is multiplied by -1 and averaged with the @xmath66 value at a position mirrored on the opposite side of the bh .
] such a procedure is necessary in order to obtain sensible results from an axisymmetric dynamical modeling code , which assumes such an underlying symmetry .
the kinematic data were bi - symmetrized in a weighted mean , using the individual uncertainties reported by ppxf to derive the weights .
the final kinematic uncertainties are computed as the errors in the weighted mean values .
figure  [ fig : nifs_sau ] shows the line - of - sight stellar velocity fields from nifs ( left ) and sauron ( right ) .
the latter data are from @xcite and show a similar large scale rotational velocity to that seen in the inner 3@xmath493  region observed with nifs .
although the large scale line - of - sight velocity fields in both panels show the same global rotational velocity , the central region of the sauron velocity field within  5@xmath49 5  ( right ) is dominated by agn continuum .
therefore , the innermost parts allow no velocity measurement , and appear as a black oval region in the image .
thus , the sauron data do not provide constraints on the bh mass , although they could provide constraints on the mass - to - light ratio ( @xmath11 ) .
however , we do not include the sauron data in our modeling , relying instead on the long - slit kinematics obtained from mmt and kpno to provide the constraints on @xmath11 at large radii .
this is partly a pragmatic choice since the use of ifu data on both small and large scales would significantly increase the total number of kinematic constraints that would need to be fitted and , therefore , would require a significantly larger orbit library than we use here , making the problem significantly more computationally expensive than we are able to handle .
we calculate the @xmath2-band surface brightness profile using our nested series of images to probe different spatial scales . for large radii ( @xmath70 )
, we use the photometrically calibrated osubsgs images . at intermediate radii ( @xmath71 )
, we use the acs optical image , allowing a significant region of overlap with the osubsgs data .
the acs flux level is then rescaled to provide the best match to the @xmath2-band data in the overlap region .
our use of the optical data could potentially bias the surface brightness profile we derive if there are significant variations in extinction or stellar population at radii of @xmath72 .
however , there are no independent indications of stellar population changes in that part of ngc  4151 , and the radial profile of the acs image is relatively smooth at those distances .
those radii are also well isolated from the psf of the agn in the acs image .    in the central regions ,
we derive a stellar flux profile directly from the nifs datacube .
one of the outputs of the ppxf modeling is the derived amplitude of the stellar template that best matches the galaxy data . as our template spectra are continuum - normalized , these amplitudes constitute relative measurements of the stellar flux in each spaxel , calculated in a way that naturally omits the agn contribution to the light profile .
the unbinned data are used to probe radii of 0@xmath73 ( excluding the agn - dominated central spaxel ) . the overall nifs stellar profile
is then scaled in amplitude to match the overlapping acs data from larger radii .
we then fit the full @xmath2-band surface brightness profile using the multi - gaussian expansion ( mge ) package @xcite .
the mge routine fits the smallest necessary number of nested 2d gaussians with a common center and radially varies their amplitudes and axis ratios to obtain the best possible smooth fit to the observed 2d surface brightness distribution , taking account of the psf .
the deprojection of a 2d surface brightness distribution to derive the 3d luminosity distribution suffers from the well known degeneracy arising from `` konus densities '' that are invisible in projection @xcite .
a number of different methods have been used to derive the 3d luminosity density , none of which guarantee a unique deprojection .
the mge method has the advantage that , once an angle of inclination is assumed , each 2d gaussian has a unique 3d luminosity density whose potential is easily computed and which is positive everywhere @xcite .
the code outputs the amplitude ( total counts ) , width in pixels ( @xmath74 ) , and observed axis ratio ( @xmath75 ) for each gaussian .
we convert those values to a central surface brightness , @xmath76 , for each gaussian as : @xmath77 where @xmath78 is the @xmath2-band surface brightness of the sun at a distance where 1  pc@xmath79 subtends 1  arcsec@xmath79 . for an absolute @xmath2-band magnitude of 3.32@xmath530.03  mag @xcite ,
this gives @xmath80  mag  arcsec@xmath81 .
the surface brightness is computed from the mge fits as : @xmath82 where @xmath83 is the osubsgs @xmath2-band image zero - point of 22.2  mag , @xmath84 is the total counts , @xmath85 is the plate scale of @xmath86  pixel@xmath35 , and the @xmath2-band extinction is @xmath87=0.016  mag ( * ? ? ? * retrieved from the nasa / ipac extragalactic database ) . for the mge modelling ,
a 1-component psf was used , having an effective width of @xmath88 .    the output of our mge fit to the surface brightness distribution
is given in table  [ tab : mge_sb ] .
for each gaussian in the mge fit , the first column gives an identifying number , the second column gives the @xmath2-band luminosity density , the third column gives the standard deviation @xmath4 ( in arcsec ) , and the fourth column gives the flattening @xmath75 ( ellipticity ) .
it is clear from the values of @xmath75 in the third column that the bulge surface brightness distribution is close to circular since most of the gaussians have @xmath89 ( except for the innermost and outermost gaussians ) suggesting that the bulge is probably close to spherical . in figure
[ fig : sb_profile ] , we show the composite radial surface brightness profile obtained by patching together the three imaging datasets .
lcccgaussian & surface density , @xmath90 & gaussian @xmath4 & axis ratio @xmath75 + number & ( @xmath91/pc@xmath92 ) & ( ) & + ( 1 ) & ( 2 ) & ( 3 ) & ( 4 ) +   + 1 & 9.3963e+04 & 0.126 & 0.7906 + 2 & 1.3200e+05 & 0.232 & 1.0000 + 3 & 7.8100e+04 & 0.549 & 1.0000 + 4 & 2.0951e+04 & 1.421 & 1.0000 + 5 & 6.7191e+03 & 2.799 & 1.0000 + 6 & 3.7009e+03 & 6.703 & 1.0000 + 7 & 6.8485e+02 & 14.513 & 0.9550 + 8 & 3.1906e+02 & 64.987 & 0.7052 +    [ tab : mge_sb ]      to provide a rough constraint on the relationship between the luminosity distribution and the stellar mass distribution , we combine optical imaging from sdss with the osubsgs @xmath2-band image .
we construct ( @xmath93 ) and ( @xmath94 ) color maps of ngc  4151 , finding the bulge colors of ( @xmath93)=1.1@xmath530.1  mag , and ( @xmath94)=2.4@xmath530.2  mag to be quite uniform across the bulge ( aside from the central few arcseconds , where the agn dominates the emission ) . @xcite calculated @xmath2-band mass - to - light ratios ( @xmath95 ) in ( @xmath93)-(@xmath94 ) color - space , based on a monte carlo analysis of @xcite stellar population synthesis models . examining our measured colors in the @xcite color - space
indicates that @xmath8 .
our initial exploration with dynamical models spanning the entire range @xmath96 showed that the best fit values were always in the range @xmath97 , hence this range of values is used in the models below .
ngc  4151 was classified by @xcite as ( r)sab(rs)ab : a disk galaxy with a pseudo - outer ring ( r ) , a `` mixed '' bar ( ab ) ( i.e. , possibly barred but is not strongly so ) , with mixed spiral / inner ring classification ( rs ) , and an intermediate to large bulge ( ab ) .
figure  [ fig : slit_orient ]  ( left ) shows the palomar observatory sky survey ( poss ) image of ngc  4151 in the red optical filter . in part because of the very circular bulge , the classification of this galaxy as barred has been questioned by authors as far back as @xcite , who speculated that the elongation was not a bar , but was in fact made up of material that was ejected from the nucleus in a former explosive event . however , we believe that ngc  4151 may be an example of a galaxy containing a `` barlens '' , an oval - shaped structure that is likely to be the vertically thick aspect of a bar seen from a particular viewing angle @xcite .    in paper
i , we assumed that the kinematic major axis is aligned with the bar , since spectral information was only available along two long - slits .
the kpno slit ( represented by the red line in fig .
[ fig : slit_orient ]  left ) was oriented along the bar at a position angle of 135@xmath52 @xcite .
the mmt slit ( green line in fig .
[ fig : slit_orient ]  left ) was placed at a position angle of 69@xmath52 .
the middle panel of figure  [ fig : slit_orient ] shows the region within the white box from the left panel .
the grey box in the middle panel represents the orientation and shape of the nifs field of view .
it is quite clear from a comparison of the nifs velocity field and the sauron velocity field in fig .
[ fig : nifs_sau ] , that the rotation axis of the disk and bulge are not perpendicular to the bar as was assumed in paper  i. this is an important difference between the models presented here and the models presented in paper  i.      to first order , the bh mass derived by the stellar dynamical modeling depends linearly on the assumed distance to the galaxy where r is the linear radius within which the velocity ( or velocity dispersion ) is measured , and the conversion between angular and linear radii is inversely proportional to distance in the nearby universe . ] . for consistency with the analysis in paper  i
, we adopt a distance for ngc  4151 of 13.9  mpc ( based on the recession velocity of 998  measured by * ? ? ?
however , as noted in paper  i , the actual distance to the galaxy is rather uncertain , with published estimates spanning a range of 10  mpc to 30  mpc .
recently , @xcite compiled a group - averaged tully - fisher distance of 11.2@xmath531.1  mpc in the extragalactic distance database . however , that group distance is based on just four galaxies , which have a standard deviation in their distance moduli of more than four times the individual distance modulus errors .
furthermore , it is unclear whether the contribution to the group average distance from ngc  4151 ( of @xmath124  mpc ) accounts for the agn contribution to the galaxy brightness .
thus , we can not regard the @xcite value as being any more reliable than other distances , and until a more accurate distance can be obtained , we continue to use the value of 13.9  mpc .
we apply what is now the standard stellar dynamical modeling approach of orbit superposition ( e.g * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) to fit the observed losvds and luminosity distribution of ngc  4151 .
the orbit superposition ( or `` schwarzschild '' ) method of dynamical modeling uses line - of - sight kinematical and imaging data to constrain the mass distributions of external galaxies .
the implementation of the code used here differs from the method detailed in ( * ? ? ?
* hereafter vme04 ) in two small aspects .
it is ( a ) designed to fit generalized axisymmetric models that are obtained from direct deprojection of the surface brightness profiles of the individual galaxies using a mge code , following the method outlined in @xcite ; ( b ) modified to fit integral field kinematics in addition to long - slit kinematics .
we restrict our models to the surface brightness distribution within 50 i.e. , the bulge region  using the mge decomposition from table  [ tab : mge_sb ] . for each dynamical model , we assume a constant @xmath2-band mass - to - light ratio ( @xmath95 ) for the stars value for the stars and the best - fit model values of @xmath95 described below are completely consistent with the expectations derived from the stellar photometry in sec .  [
sec : mtol ] .
even if we were to include such a halo , its contribution to the inner 15 region , where the nifs kinematics exist , would likely be small .
as we will show , the nifs data entirely determine our final bh mass , justifying our decision to not include a dark matter halo .
] , with values varying from @xmath98 .
the resulting density profile is then converted to a gravitational potential using an axisymmetric multipole expansion scheme @xcite . the gravitational potential of a point mass ( @xmath0 ) , representing the central bh , is then added .
models are constructed for every pair of parameters ( @xmath0 , @xmath95 ) , building an orbit library for each combination .
the large - scale disk of ngc  4151 is inclined to the line of sight at @xmath99 @xcite . however , since the bulge appears very close to circular in projection ( see fig .  [
fig : slit_orient ] and tab .
[ tab : mge_sb ] ) , it is probably nearly spherical or being viewed close to face on .
thus , we run models on a 2d grid of ( @xmath100 ) values for three different inclination angles @xmath101 ( edge - on ) , @xmath102 and @xmath103 ( close to face - on ) .
we also run a limited set of models for which the value of @xmath0 is held fixed while the inclination angle of the model is varied with values @xmath104 , and @xmath105 , in an effort to set constraints on the inclination of the bulge to the line - of - sight .    for inclination angles of @xmath101 and @xmath106 , we run orbit libraries for 18 bh masses between @xmath107 , using a baseline @xmath108 . these models
are then scaled to create libraries for 18 different @xmath95 values , producing a total of 324 models for each inclination angle .
we also run orbit libraries for @xmath102 , using 10 bh masses across the same mass range and 18 @xmath95 values for each @xmath0 ( 180 models ) .
the range of bh mass values was set by initial experimentation with a coarse grid on a wide range of values .
the lowest two values of @xmath0 are set to 10@xmath109 and 10@xmath110 , both of which consistently give poor fits to the data .
the upper limit of @xmath0 was raised until it was clear that models once again gave very poor fits to the data .
the combined gravitational potential of the bulge and the bh is used to integrate a large library of @xmath111 stellar orbits selected on a grid with @xmath112 energy values , @xmath113 angular momentum values ( @xmath114 ) at each energy , and @xmath115 pseudo third integral values ( @xmath116 ) at each angular momentum value ( see vme04 for details ) .
for instance , models with @xmath117 used @xmath118 , @xmath119 , @xmath120 .
each orbit is integrated for 100 orbital periods and its average `` observed '' contribution to each of the `` apertures '' in which kinematical data are available is stored .
the orbital kinematics are subjected to the same type of instrumental effects as the real data ( i.e. psf convolution , pixel binning ) .
self - consistent orbit - superposition models are constructed by linearly co - adding orbits in each library to simultaneously fit several elements within the observed apertures ( imaging and spectroscopic ) : the 3d mass distribution ( for each value of @xmath95 ) , the projected surface density distribution , and the losvds .
the fits use a non - negative least squares optimization algorithm ( nnls ; * ? ? ?
* ) , which gives the weighted superposition of the orbits that best reproduce both the self - consistency constraints ( i.e. , the 3d mass distribution and surface brightness distribution ) and the observed kinematical constraints ( @xmath121 ) .
vme04 showed that as the ratio of orbits to constraints decreases , the error bars decrease artificially . in the majority of the models presented here we have 1240 constraints that are fitted by a library of 9936 orbits .
this size of orbit library gives an orbits - to - constraints ratio of 8 which is only slightly above the ratio ( 5 ) at which vme04 find that solutions begin to be biased . in section  [ sec : robust_size ] we carry out a limited exploration with a 50% larger orbit library and find that our results are unchanged .    smoothing ( or `` regularization '' ) of the orbital solutions
is typically performed in this type of analysis , and a variety of methods have been employed to that end .
regularization mainly helps to reduce the sensitivity of the solutions to noise in the data by requiring a smoother sampling of the orbit libraries used in the solution .
following @xcite , we run a trial series of models with a regularization scheme that requires local smoothing in phase space .
we find that for small smoothing parameters ( @xmath122 ) the resultant best - fit @xmath0 values do not differ from the models without regularization .
however , the @xmath1 value for the best - fit model is larger and the error bars on @xmath0 are smaller .
vme04 showed that it is impossible ( without computationally expensive approaches such as generalized cross - validation ) to determine the ideal regularization parameter . not regularizing the solutions can yield large error bars , but choosing too large
a value of the smoothing parameter can introduce a bias in the best - fit value of @xmath0 . in this paper , we present models without regularization , since imposing regularization constraints significantly increases the total number of constraints and therefore decreases even further the orbits - to - constraints ratio .
for the libraries used here , models without regularization provide the most conservative error bars on the estimated value of @xmath0 .
since the mass and surface mass density constraints are _ self - consistency constraints _ rather than observational constraints , they do not have observational errors associated with them .
there are two ways in which the absence of observational errors are handled in the literature . @xcite and vme04 include self - consistency constraints in the nnls minimization problem with errors corresponding to predefined level of relative accuracy in the fit ( e.g. 5% -10% ) , while others @xcite solve the orbit superposition problem by requiring the self - consistency constraints ( mass and surface brightness ) to be fitted as separate as linear ( in)equality constraints with predetermined absolute accuracy .
@xcite state that they require the self - consistency constraints to be fitted to an accuracy of 0.02 .
we also attempted to fit self - consistency constraints as linear equalities but find that that since the values of the mass constraints vary by 5 orders of magnitude ( from @xmath123 to 10 ) , using a fixed numerical accuracy fails to fit the mass constraints at small radii ( where the logarithmic radial bins have the smallest mass values per bin and where the fit to mass constraints most strongly affects the estimated @xmath0 ) .    for the models presented here
we include the self - consistency constraints in the optimization problem and require the optimization problem to fit each of the 472 self - consistency constraints ( 280 mass constraints and 192 aperture surface brightness constraints ) with an error corresponding to a relative accuracy - per - constraint of 5% . for our preferred model with @xmath101 ,
the fractional accuracy per self - consistency constraint varies from 0.27% ( for the best fit model ) to 2.2% for the worst fit model , i.e. , always better than our required accuracy .
figure  [ fig : mass_err ] shows contours of constant fractional error in the fit to the self - consistency constraints ( mass and surface brightness in apertures ) over the 2d parameter space ( @xmath0 , @xmath95 ) for models with @xmath101 .
the star is the location of the model with 0.27% accuracy per self - consistency constraint , and the contours are spaced at intervals of 0.1% relative accuracy per constraint .    in previous applications of our schwarzschild modeling code ( vme04 , * ? ? ?
* paper  i ) , the error in the fit to the self - consistency constraints was essentially random over the entire 2d parameter space ( @xmath124 ) .
although the self - consistency constraints were always included in the @xmath1 , they contributed so little to the total @xmath1 that the best fit solution was always driven by the fit to the kinematical constraints .
we will show that this is not true for our current ngc  4151 dataset . despite the fact that the fits to the 3d mass distribution and surface brightness distribution are very good , because the fit vary systematically rather than randomly across the parameter space , the best fit solution is altered by the decision of whether or not to include the fit to the self - consistency constraints in the @xmath1 .
based on the arguments above , we compute the @xmath1 in two different ways : ( 1 ) the total @xmath125 includes all the constraints in the nnls optimization problem : 280 mass constraints , 192 surface brightness constraints , and 768 kinematical constraints ( giving a total of @xmath126 constraints ) ; ( 2 ) @xmath127 which only considers the quality of the fit to @xmath128 kinematical constraints ( @xmath129 and @xmath67 in 192 apertures : 136 nifs spaxels , 15 kpno apertures , and 41 mmt apertures ) .
for each of the three assumed angles of inclination for the bulge , we construct models by fitting these constraints using an orbit library consisting of @xmath117  orbits . for one inclination angle ( @xmath101 ) we do a limited exploration with a library consisting of @xmath130 orbits . after generating orbit libraries for a grid of parameters ( @xmath131 )
the best fit parameters are determined by the minimum in the 2d @xmath1-contour plot .
in this section , we present the results of our modeling . in section  [ sec : define_chi2 ]
we describe results for inclination angles @xmath101 , @xmath102 and @xmath103 , showing 2d contour plots of @xmath125 and @xmath127 using a full 2d grid of models for each of these inclination angles . in section  [ sec : inclination ]
we examine the dependence on inclination angle for 4 additional models using a single @xmath0 . in section  [ sec : robust ] we explore the robustness of our solutions by examining their sensitivity to the size of the orbit library and to the effect of restricting the set of kinematic constraints to those from the nuclear region only .
we also present some of the best - fit solutions and show how well different pairs of @xmath100 fit both the 1d and 2d kinematical data . in section 5.4
, we discuss our results in the context of the relation .
figure  [ fig:2dchi2 ] shows 2d contour plots of @xmath125 in the plane of model parameters @xmath0 and @xmath95 for @xmath101 ( top ) , @xmath102 ( middle ) , and @xmath103 ( bottom ) .
the grid of points shows the model parameters ( @xmath100 ) for which the orbit libraries is constructed and the optimization problem is solved .
the first 6 contour levels in all the @xmath1 contour plots that follow correspond to @xmath132 , 68.3% confidence ) , @xmath133 , 90% confidence ) , @xmath134 , 95.4% confidence ) , @xmath135 , 99% confidence ) , @xmath136 , 99.73% ) , and @xmath137 , 99.99% ) , with subsequent contours equally spaced in @xmath138 .    for @xmath101 ,
the minimum value of @xmath125 is obtained for @xmath139 and @xmath140 , with @xmath141 . for inclination angle of @xmath102 ,
the best fit model has @xmath142 and @xmath143 , with @xmath144 . for inclination angle of @xmath103 ,
the best fit model has @xmath145 and @xmath140 , with @xmath146 .
table  [ tab : bestfit ] summarizes the results and includes minimum @xmath1 values and the 1@xmath4 and 3@xmath4 error bars .
the errors on a given parameter are obtained by marginalizing the 2d @xmath147 values over the other parameter .
figure  [ fig:1dchi2 ]  ( top ) shows the marginalized 1d @xmath138 versus @xmath0 for the three different inclination angles , as indicated by the legends .
the horizontal dotted lines are at @xmath148 ( 1@xmath4 , 68.5% confidence interval for 1 degree - of - freedom ) , @xmath149 ( 2@xmath4 , 90% confidence interval ) , and @xmath150 ( 3@xmath4 , 99% confidence interval ) .
the 1@xmath4 and 3@xmath4 errors on @xmath0 and @xmath95 are given in table  [ tab : bestfit ] .
likewise , figure  [ fig:1dchi2 ]  ( bottom ) shows the marginalized 1d @xmath138 versus @xmath95 for each of the three inclinations of models to the line - of - sight .
it is important to note that with the standard methods used to derive the velocity profiles ( vps ) from the spectra , the errors in the gauss - hermite ( gh ) coefficients are not independent .
the errors associated with the even moments ( @xmath4 and @xmath67 ) are correlated , and errors associated with the odd moments ( @xmath151 and @xmath66 ) are also correlated ( see , e.g. , * ? ? ?
* ; * ? ? ?
@xcite point out that an important consequence of correlated errors in the gh coefficients is that the errors obtained by comparing data to dynamical models using the @xmath1 estimator are significantly underestimated . in this paper
we quote both 1-@xmath4 and 3-@xmath4 errors , and often consider the 3-@xmath4 errors as better representations of the true modeling errors .
lcccccccc & & & & + & & & & & & & & + ( 1 ) & ( 2 ) & ( 3 ) & ( 4 ) & ( 5 ) & ( 6 ) & ( 7 ) & ( 8) & ( 9 ) +   +   + 90@xmath52 & 9936 & 709.1 & 4.68 & @xmath152 & @xmath153 & 0.304 & @xmath154 & @xmath155 + & & & & & & & + 60@xmath52 & 9936 & 767.5 & 5.42 & @xmath156 & @xmath157 & 0.305&@xmath158 & @xmath159 + & & & & & & & + 23@xmath52 & 9936 & 1050.8 & 3.69 & @xmath160 & @xmath161 & 0.304 & @xmath158 & @xmath155 + & & & & & & & +   +   + 90@xmath52 & 9936 & 186.7 & 7.76 & @xmath162 & @xmath163 & 0.332 & @xmath164 & @xmath165 + & & & & & & & + 60@xmath52 & 9936 & 193.1 & 8.62 & @xmath166 & @xmath167 & 0.326&@xmath159 & @xmath168 + & & & & & & & + 23@xmath52 & 9936 & 348.2 & 8.50 & @xmath169 & @xmath170 & 0.329 & @xmath171 & @xmath172 + & & & & & & & +    [ tab : bestfit ]    the best fit solutions for both @xmath0 and @xmath95 for the three different inclination angles are consistent with each other within 3@xmath4 . however @xmath173 is the smallest for the edge - on model ( @xmath101 ) , increases for @xmath102 and is significantly larger for @xmath103 implying that it may be possible to constrain the inclination of the bulge ( see section  [ sec : inclination ] ) .
marginalizing @xmath125 over the three values of inclination gives a best fit @xmath174 and @xmath175 .
figure  [ fig:2dchi2_kin ] shows 2d contour plots of @xmath176 for the same sets of models as in figure  [ fig:2dchi2 ] , with inclination angles @xmath177 ( from top to bottom ) .
we emphasize that the model solutions that are used to generate figure  [ fig:2dchi2_kin ] are identical to those used to generate figure  [ fig:2dchi2 ]  only the quantities used in computing the @xmath1 of the fit to the data differ .
figure  [ fig:1dchi2_kin ] shows the same data marginalized over @xmath95 ( top ) and marginalized over @xmath0 ( bottom ) .
it is clear from these two figures that the best - fit solutions obtained using @xmath178 give both a larger @xmath179 and a larger @xmath180 .
we will discuss the possible causes for this in section  [ sec : robust ] and section  [ sec : discuss ] .
in contrast to the results of paper  i , where @xmath101 models gave an upper limit of @xmath181 and @xmath103 models gave a best fit @xmath0 between @xmath182 , figures  [ fig:1dchi2 ] and [ fig:1dchi2_kin ] show that best - fit @xmath0 is relatively insensitive to inclination ( although the actual best - fit value depends on how @xmath1 was computed ) .
table  [ tab : bestfit ] shows that the minimum values of @xmath183 and @xmath184 for @xmath101 are lower than the corresponding @xmath1 values for the other two inclination angles .
we now try to constrain the inclination angle of the bulge of ngc 4151 by assuming a constant bh mass of @xmath185 . with this value of @xmath0
we compute orbit libraries for four additional inclination angles : @xmath186 , @xmath187 , @xmath188 , and @xmath105 .
figure  [ fig : chi_inc ] shows the values of @xmath183 for models with this value of @xmath0 and @xmath189 as a function of the inclination of the bulge to the line - of - sight .
it is clear that @xmath183 decreases steadily with increasing inclination angle with a minimum at @xmath101 .
( a similar dependence on inclination angle is obtained for @xmath190 , although all the models give much worse fits as assessed by the @xmath125 values . )    the inclination angle of @xmath103 that is inferred from the velocity field is strongly disfavored . since the disk of this galaxy is clearly seen to be close to face - on , the edge - on orientation of the bulge preferred by our models is rather puzzling : the surface brightness isophotes of the bulge are nearly circular , hence the best - fit edge - on orientation implies that the bulge is nearly spherical in shape ; but in this case the rotation axis of the bulge
must be misaligned with the rotation axis of the disk by @xmath191 .
previous determinations of the orientations of the various components of ngc  4151 have shown a remarkable diversity . while the circular isophotes of the large - scale disk suggest the nearly face - on orientation with @xmath103 @xcite  which motivated our testing of that value of inclination  the components associated with the agn
have been inferred to lie further from our line of sight .
models of the bi - conical narrow - line region ( nlr ) in the uv and optical @xcite , and in the near - ir @xcite have typically found inclinations of @xmath192 .
closer to the agn , @xcite fit the absorption using a pc - scale nuclear disk inclined at @xmath193 . on even smaller scales ,
an inclination of @xmath194 for the agn accretion disk was obtained from modeling the x - ray spectrum @xcite , while @xcite suggest that the radio jet in ngc  4151 lies very close to the plane of the sky and is interacting directly with circumnuclear gas in the plane of the galaxy ( i.e. , @xmath195 ) .
this lack of alignment between the different spatial scales could be the consequence of a merger in the galaxy s recent past .
it has been claimed that the periodic variability in the nuclear activity of ngc  4151 is indicative of a bh binary , with a separation on the order of @xmath196  au @xcite .
while distinguishing between a single bh and a close binary is far beyond the capability of our dataset , the remnants of a merger could also contribute to the complicated dynamics we observe on the size scales of the bulge . in section  [ sec : robust ] , we will show that the kinematics in the inner 5  show evidence for a bar in this region , implying that the internal dynamics of the bulge are likely not that of a spherical or axisymmetric system .
a full understanding of why the kinematics of our model prefers an edge - on configuration to the more inclined configuration requires the system to be modeled by a bar dynamical modeling code , which does not exist at the present time .    since the result of
this section show that models with edge - on inclination ( @xmath101 ) are greatly preferred over more face - on models , in the rest of this paper we confine our analysis to models with @xmath101 .
lllllcccccc & & & & & & + & & & & & & & & & & + ( 1 ) & ( 2 ) & ( 3 ) & ( 4 ) & ( 5 ) & ( 6 ) & ( 7 ) & ( 8) & ( 9 ) & ( 10 ) & ( 11 ) +   + @xmath183 & 15092 & 1240 & 641.9 & 0.52 & 4.68 & @xmath197 & @xmath198 & 0.31 & @xmath158 & @xmath199 + & & & & & & & & & & + @xmath178&15092 & 768 & 170.2 & 0.22 & 7.32 & @xmath200 & @xmath201 & 0.31 & @xmath199 & @xmath202 + & & & & & & & & & & +   +   + @xmath183 & 9936 & 960 & 583.1 & 0.61 & 3.82 & @xmath203 & @xmath204 & 0.32 & @xmath199 & @xmath205 + & & & & & & & & & & + @xmath178 & 9936 & 680 & 101.7 & 0.14 & 3.70 & @xmath206 & @xmath207 & 0.35 & @xmath199 & @xmath202 +    [ tab : robustbestfit ]    vme04 showed that the solution derived from the schwarzschild method was sensitive to two important factors : ( a ) the size of the orbit libraries used in obtaining the best - fit bh mass ; and ( b ) the spatial resolution of the kinematic data and whether or not they resolved the sphere - of - influence of the bh .
they showed that biased solutions could result when the size of the orbit libraries was too small for the available data .
they also showed that when the kinematic data did not resolve the sphere - of - influence of the bh , spurious results could arise . in the next two sections we perform two tests to assess the robustness of our solutions .
in particular , our goal is to assess which of the two values of @xmath0 obtained previously ( 5@xmath208 or 8@xmath209 )  is more robust .      in this section ,
we construct a set of 288 models with @xmath101 but with orbit libraries containing @xmath210 orbits ( as compared to @xmath211 used thus far ) to fit the same set of 1240 constraints that were fitted in section  [ sec : define_chi2 ] .
large orbit libraries are constructed for 16 bh masses between 10@xmath109  and 15@xmath212  and for 18 values of @xmath213 .
figure  [ fig:1dchi_i9015k ] shows the 1d @xmath183  ( solid curves ) and @xmath178  ( dashed curves ) obtained by marginalizing over the mass - to - light ratio @xmath95 ( top panel ) and obtained by marginalizing over @xmath0 ( bottom panel ) .
the resulting best - fit values for @xmath0 and @xmath95 , and their 1@xmath4 and @xmath214 errors , are shown in the top two rows of table  [ tab : robustbestfit ] .
the reduced @xmath1 values are given in column five .
it is clear from a comparison with results in table  [ tab : bestfit ] that the best - fit @xmath0 values are relatively insensitive to an increase in the size of the orbit library by @xmath1250% .
the best - fit value of @xmath95 obtained by using @xmath178 ( dashed curves in the bottom panel of fig .
[ fig:1dchi_i9015k ] ) is somewhat lower than the value of @xmath215 obtained with the smaller library , but is still within 2@xmath4 of this value .
the test with the larger orbit libraries confirms the result we obtained with @xmath117 orbits per library , implying that our best - fit parameters are not biased by having too small an orbit library .
furthermore , the discrepancy between the best - fit values obtained using @xmath147 and @xmath178 persists even with larger orbit libraries .
figure  [ fig : kinenifs_90 ] shows examples of model fits to the 2 dimensional nifs velocity fields ( from left to right , , and the gh moments @xmath66 and @xmath67 ) for the edge - on ( @xmath101 ) case with @xmath210 .
note that while the best - fit solutions listed in tables  [ tab : bestfit ] & [ tab : robustbestfit ] are determined by marginalizing over one parameter in the 2d contour surface ( which itself is smoothed with a kernel ) , the 2d grid of models is discrete .
the derived best - fit value ( shown as stars on the 2d @xmath1 maps of figs .
[ fig:2dchi2 ] and [ fig:2dchi2_kin ] ) does not overlap with any model on the grid , but lies between 4 grid points .
we therefore show velocity fields and 1d kinematics for the model whose parameters are closest to the best - fit solution .    in figure
[ fig : kinenifs_90 ] , the top row shows velocity fields for @xmath216 ( the model closest to the best - fit obtained with for @xmath210 ) and the middle row shows velocity fields for @xmath217 , @xmath218 ( the model closest to the best - fit obtained with for @xmath210 ) .
( the bottom row of this figure is for a model with @xmath219 , @xmath218 , obtained when only nifs kinematical data are fitted , and will be discussed in section  [ sec : robust_nuclear ] . )
the fitted velocity in figure  [ fig : kinenifs_90 ] should be compared with the bi - symmetrized velocity fields in the bottom row of figure  [ fig : vel_fields_orig ] .
the four white pixels in each panel of the figure were not fitted because it was determined that their extremely high velocity dispersion values ( @xmath220  ) were spurious ( see fig .  [ fig : vel_fields_orig ] ) since their values were larger than the central velocity dispersion despite being far from the center .    in figure
[ fig : kinenifs_90 ] the odd - moments of the losvd ( , @xmath66 ) are fairly well fitted by both the low @xmath0 model ( top row ) and the high @xmath0 model ( middle row ) . however , the nifs velocity dispersion  for the model with @xmath221 , @xmath108 ( top row ) is much lower than  in the data ( bottom of fig .
[ fig : vel_fields_orig ] ) . the model with the larger value of @xmath217
does a slightly better job of reproducing the higher  values within @xmath222 , but fails to fit the  at the edges of the nifs field , where model values are significantly lower ( i.e. bluer pixel colors ) than the observed values seen in figure  [ fig : vel_fields_orig ] ( bottom ) .
this is because  from nifs in the inner 15 region is overall larger than that obtained with kpno ( @xmath223  ) and mmt ( @xmath224  ) .
inconsistent velocity dispersion values in the _ same physical region
_ drives the solution to fit the data with the smallest error bars ( i.e. , mmt data ) thereby under - estimating  from nifs . also note that neither model is able to fit both the low central and high outer @xmath67 values seen in the bi - symmetrized nifs velocity fields ( fig .
[ fig : vel_fields_orig ] bottom row , right - most panel ) .
the difference between the model fits and the data is more clearly seen in figure  [ fig : kinematics_90_1d ] which shows the one - dimensional line - of - sight velocity ( ) , velocity dispersion ( ) and the gh moments of the losvd ( @xmath66 and @xmath67 ) that we obtain for the edge - on ( @xmath101 ) model .
the red and blue curves show the best fits for same models as in the top and middle panels of figure  [ fig : kinenifs_90 ] with values and @xmath95 indicated by the line legends in the top left panel of the first row ( these two models used libraries with @xmath210 orbits ) .
the green curves in the top four nifs panels are obtained when only the nifs data are fitted ( with @xmath211 ) , this model is discussed in the next section .    the open circles in the panels show the three different kinematic datasets used in the modeling , along with their 1@xmath4 uncertainties .
( the top four panels show kinematics in the nifs spaxels lying closest to the major axis of the model ; the middle four panels show data from the kpno long - slit ; the bottom four panels show data from the mmt long - slit . )
the most striking feature of the fits is that the red curve for @xmath225 ( red curve )  the model closest to the best - fit solution derived from the minimum in @xmath183 ( from figure .
[ fig:1dchi2 ] )  clearly _ underestimates _ the nifs velocity dispersion ( top right panel ) for the central - most points .
the uniformly large central @xmath226 values in the inner @xmath227 region are much better fitted by larger values of @xmath228 and @xmath218 ( blue curve ) .
both the models do an equally good job of fitting , , and @xmath66 from the mmt data , but overestimate  from the kpno data .
another important point to note is that the values of the gh parameters @xmath67 are _ negative _ in the innermost regions of each of the three datasets .
neither model is able to fit the negative @xmath67 values measured with the mmt spectrograph nor the central dips to negative values in @xmath67 seen in the nifs data and kpno data .    in spherical and axisymmetric galaxies ,
@xmath67 is negative when the velocity distribution is tangentially biased ; it is positive when orbits are predominantly radial ; and it is zero when the velocity distribution in isotropic @xcite .
however , negative @xmath67 values have also been associated specifically with the presence of a bar . in recent @xmath229-body simulations of barred galaxies by @xcite ,
the losvds have negative @xmath67 values in the region where the bar dominates .
the negative @xmath67 parameters arise due to the kinematic properties of the special orbits that constitute the bar , and due to the fact that from an external observer s point of view , the bar has a large tangential velocity component due to its pattern speed . even after a bar buckles and forms a boxy bulge ,
negative @xmath67 values can be seen in face - on systems @xcite .
another signature of the kinematics associated with a nuclear bar is undulations in the rising part of the inner rotation curve @xcite , which we also see in the mmt data for ngc  4151 . finally ,
while @xmath66 is always anti - correlated with the line - of - sight velocity  in axisymmetric systems , @xcite show that @xmath66 is _ correlated _ with  along the major axis of the bar . for the data obtained with kpno ( where the slit lies along the major axis of the large scale bar ) @xmath66
is correlated with  over the radial range @xmath230 ( figure  [ fig : kinematics_90_1d ] ) , providing additional kinematical evidence for a bar .
therefore it appears that despite the very circular photometric contours of the bulge in ngc  4151 , and the visual appearance of a rather weak bar at radii outside the bulge , there are several clear kinematical signatures of a small scale bar in the vicinity of the nuclear bh .
this is consistent with the properties of `` barlens '' features that have recently been identified in real and simulated galaxies @xcite . as pointed out by @xcite , barlenses frequently masquerade as classical bulges because of their nearly circular central isophotes , especially when viewed nearly face - on , but have bar - like kinematics
.    the combination of bars and bhs can have significant implications for the central dynamics of galaxies . @xcite
analyze a suite of @xmath229-body simulations of barred galaxies in which point masses representing bhs were grown adiabatically .
they find that the growth of a bh of a given mass causes a 5 - 10% larger increase in @xmath63 in a barred galaxy than in an otherwise identical axisymmetric galaxy .
@xcite find that if a disk galaxy with a bulge has a pre - existing black hole , the formation and evolution of a bar can result in a 15 - 40% increase in the central velocity dispersion .
both studies find that the increase in the observed line - of - sight velocity dispersion in the barred systems is a consequence of three separate factors : ( a ) mass inflow due to angular momentum transport by the bar , ( b ) velocity anisotropy due to the presence of bar orbits , and ( c ) weak dependence on orientation of position angle of the bar .
if ngc  4151 does have a bar and is thus non - axisymmetric , that fact also has implications for our modeling results . in stellar dynamical modeling of axisymmetric galaxies ,
fitting a negative @xmath67 requires a larger fraction of tangential orbits , which contribute little to the line - of - sight velocity dispersion . fitting a given velocity dispersion with such an orbit population
will imply a larger enclosed mass than would be required if @xmath67 was zero or positive ( which is typically the case in axisymmetric systems with bhs ) .
@xcite argue that if an axisymmetric stellar dynamical modeling code is used to derive the value of @xmath0 in a barred galaxy with a high central @xmath231 and a negative @xmath67 , the central bh mass will be systematically overestimated .
if the large values of @xmath231 in the nuclear regions and the negative @xmath67 values in the mmt data are a consequence a nuclear bar , then the discrepancy between @xmath127 and @xmath125 could reflect this predicted bias in @xmath0 .      in this section
we examine the consequences of fitting only the nifs kinematical data ( ignoring the kinematical constraints at larger radii obtained from the lower spatial resolution long - slit spectra from kpno and mmt ) .
there are three reasons why the nifs data alone may provide better constraints on the mass of the bh .
first , the nifs data have higher spatial resolution ( 02 ) than the long slit data ( 1 - 2 ) . for
a @xmath232 , the sphere - of - influence of the bh ( assuming @xmath233  ; see section  [ sec : msigma ] below ) is @xmath1216  pc , which corresponds to @xmath12 024 for our assumed distance of 13.9  mpc for this galaxy .
therefore the nifs data used in our modeling are barely able to resolve the sphere - of - influence of the bh , but are quite close to doing so .
if @xmath234 , the sphere - of - influence of the bh would be 26  pc ( 038 ) and the bh sphere - of - influence would be resolved by the nifs data .
the long slit data do not provide constraints on @xmath0 and are included mainly to provide constraints on the mass - to - light ratio of the stars , which is assumed to be independent of radius .
second , as discussed previously , the nifs velocity dispersion within @xmath235 is significantly larger than the velocity dispersion on same physical scale obtained from the long - slit data .
however , since the error bars on the long - slit mmt data are about half those from the nifs data , the optimization code tries to fit them better than the nifs data .    in figure
[ fig : kinenifs_90 ] ( top two panels ) , the need to fit the lower velocity dispersion values from mmt manifests as low velocity dispersion values ( blue ) at the edges of the nifs field .
these low dispersion values outside the sphere - of - influence of the bh also force the model to adopt a lower mass - to - light ratio , which is then compensated for in the inner region by requiring a much larger to fit the nifs velocity dispersion . neglecting the long - slit data
would allow the optimization code to raise the mass - to - light ratio in the inner @xmath532   region ( where the velocity dispersion is quite flat ) .
finally , in figure  [ fig : kinematics_90_1d ] , we see kinematic evidence for the bar in mmt and kpno long - slit data ( negative @xmath67 values and undulations in the rising part of the rotation curve seen in mmt data ; correlation between @xmath66 and @xmath236  in kpno data ) occurring on scales of 3 - 5 , well outside the nuclear region . except for the central value of @xmath67 in the nifs dataset ( which , as we see from fig .
[ fig.spec ] , is contaminated by the agn ) , all the other @xmath67 values are positive , which is what is expected for orbits around a central bh . neglecting the long - slit kinematical data could help overcome biases introduced by the bar kinematics . assuming that the bh mass should primarily be constrained by kinematic constraints at the smallest radii we now test the effect of not including the long - slit data in the nnls optimization problem .
the expectation is that neglecting the long - slit data should bring the best - fit solutions obtained from @xmath183 and @xmath178 closer to each other .
we construct a full sets of models ( 18 values of @xmath0 and 18 values of @xmath95 ) in which we only fitted the nuclear kinematical data obtained with the nifs spectrograph , with an orbit library of @xmath211 .
figure  [ fig : nifs_only ] shows @xmath138 values obtained by marginalizing over @xmath95 ( top ) and by marginalizing over @xmath0 ( middle ) .
the solid curves show @xmath237 while the dashed curves show @xmath238 .
the dotted horizontal lines show the 1@xmath4 , 2@xmath4 and 3@xmath4 confidence levels , respectively ( for 1 degree - of - freedom ) .
the top plot shows that the minimum @xmath178 ( dashed curve ) has decreased quite significantly and is now in agreement with the minimum @xmath183 ( solid curve ) .
the best - fit values of @xmath95 ( middle panel ) are still inconsistent within 1@xmath4 but are consistent within 3@xmath4 .
the bottom - most panel shows the error on the fit to the mass and surface brightness distribution .
contours are spaced at intervals of 0.1% error per constraint .
it is clear that the fit to the self - consistency constraints is uniform over a much larger range of @xmath0 and @xmath95 values implying that the fit to the mass is no longer a significant factor in determining the location of the minimum .
the best fit black hole mass and m / l ratio obtained using @xmath237 and @xmath238 , along with the reduced @xmath138 and 1@xmath4 and 3@xmath4 error bars , are shown in the bottom two rows of table  [ tab : robustbestfit ] . since the best fit @xmath0 values and their 1@xmath4 error bars from both methods of computing @xmath1 are nearly the same we use the average of the two best - fit values as the solution and sum their 1@xmath4 errors in quadrature to obtain @xmath239 .
figure  [ fig : nifs_only ] shows that the @xmath95 values obtained from @xmath178 are consistent with the best fit value obtained with @xmath183 at the 3@xmath4 level , but both are well within the expected range @xmath240 .
therefore we use the average of the two best - fit m / l values and sum their 3@xmath4 errors in quadrature to get @xmath241 .
the bottom panel of figure  [ fig : kinenifs_90 ] shows a 2d kinematic map for a model with @xmath219 and @xmath242 , with parameters closest to the best - fit solution obtained with @xmath125 ( see 3rd row of table  [ tab : robustbestfit ] ) .
this model achieves adequately large values of @xmath231 ( see red / orange regions ) over the central nifs field by raising the m / l ratio very slightly from @xmath243 to @xmath244 .
the 1d kinematic fit of this model to the nifs data along the kinematic major axis of the model is shown by the green curves in figure  [ fig : kinematics_90_1d ] , from which it is also clear that the lower value of @xmath0 is able to fit the high velocity dispersion values even better than the high bh mass when the model attempted to fit both the nifs and long - slit kinematic data .
this test largely confirms our hypothesis that the two main causes of the large value of @xmath0 obtained in the previous section were ( a ) inconsistencies between @xmath231 values from nifs and the long - slit data in the same spatial region , ( b ) the need to fit the low @xmath67 values and high central @xmath231  values with a low m / l ratio .
thus , we see that by not including the kinematics data with poorer spatial resolution ( which are clearly unaffected by the bh as can be seen in figure  [ fig : kinenifs_90 ] ) we are able to simultaneously fit the nifs kinematics and the mass and surface brightness distributions with a lower value of @xmath0 .      in this section ,
we determine how the new nifs kinematic measurements in the nuclear region of ngc  4151 affect the galaxy s location in the @xmath15 relation .
the aperture used to define the appropriate @xmath4 value for the relation is a historically contentious issue ( e.g. * ? ? ?
* ; * ? ? ?
* ) , but the nifs data naturally lend themselves to one of those standards : @xmath245 , which is the velocity dispersion within an aperture of radius @xmath246/8 , where @xmath246 is the bulge effective radius @xcite .
( despite having to adopt one particular standard , we note that @xcite found no systematic differences between estimates of @xmath245 ( inside of @xmath246/8 ) and estimates of @xmath4 that extended out to @xmath246 , likely because of the steep radial surface brightness gradient in most bulges . )
the @xmath246 value for ngc  4151 has been measured to be @xmath247 @xcite , while the nifs field - of - view extends to a radius of between @xmath86 and @xmath248 .
thus , we sum the nifs datacube across all of the spaxels into a single spectrum , and fit for the dispersion using ppxf .    from the nifs data , we find @xmath249  .
this is somewhat larger than either the @xmath250  measured by @xcite , or the @xmath251  measured by @xcite . both of those earlier measurements used the calcium triplet absorption lines , but whereas the former covered a similar region of the galaxy as our nifs spectra ( @xmath252 ) , the latter used an aperture of @xmath253 .
previous studies of bulge kinematics have found near - ir and optical data to give consistent results ( see * ? ? ? * and references therein ) , while the discrepancies seen in some luminous ir galaxies are in the sense of _ smaller _ @xmath4 values from the co - bandheads @xcite .
our new measurement of @xmath254  is @xmath1225% above the previous velocity dispersion measurements , while our best - fit stellar dynamical estimate of @xmath255 is 17% below our earlier stellar dynamical estimate @xcite .
previously , ngc  4151 had been an outlier from the agn @xmath15 relation in the direction of low-@xmath4 or high-@xmath0 .
the new @xmath245 and @xmath0 measurements helps to bring ngc  4151 closer to the best - fit relation of @xcite , though it does remain on the low-@xmath4 side .
finally , we note that if the bh masses in barred galaxies have been overestimated due to the use of axisymmetric dynamical modeling codes , such galaxies would lie even further below the standard relation than has been found previously ( * ? ? ? * and references therein ) .
this demonstrates the importance of developing robust bar dynamical modeling approaches in the future .
we have conducted ao - assisted near - ir integral field spectroscopy of the local seyfert galaxy ngc  4151 , using the nifs instrument on gemini north .
we used an axisymmetric orbit - superposition code to fit the observed surface brightness distribution within 50 , nifs kinematics within @xmath5315 , and long slit kinematics along two different position angles in order to estimate the best - fit values of the mass of the bh ( @xmath0 ) and the mass - to - light ratio of the stars ( @xmath95 ) .
models were constructed for 10 - 18 values of @xmath0 , 18 values of @xmath95 and 3 different inclination angles .
the main results of this paper are summarized below .
we use the @xmath1 estimator to determine the model which gives the best - fit to the data but find that the solution depends on whether the @xmath1 includes both self - consistency constraints ( which are not strictly speaking `` data '' ) and kinematic constraints or only kinematic constraints .
when we use both types of constraints and all the available kinematical data ( from low - resolution long - slit data and high resolution nifs data ) the best - fit solution is @xmath256 and @xmath175 , for a model with @xmath101 .
however , this best fit model gives a central @xmath231 that is too low to fit the data obtained with the nifs instrument .
when we use @xmath127 ( determined by only considering the fit to the kinematic data ) , the best fit model has @xmath257 and @xmath258 , and provides a better fit to the nuclear kinematical data from nifs spectrograph , although it gives a slightly worse fit to the mass and surface brightness distributions .
an interesting point worth noting is that @xcite found a best - fit value of @xmath259 in ngc 4151 when they fitted the kinematics of the h@xmath260 line - emitting gas within 1  of the center ( see figure 43 in their paper ) .
however , if they included the kinematics of gas within 2  they obtained a larger value of @xmath261 , although with a somewhat worse reduced-@xmath1 .
it is intriguing that both our stellar dynamical modeling and their gas dynamical modeling , while preferring lower values of @xmath262 , also suggest that _ higher _ values of @xmath0 may be obtained when lower spatial resolution data are included in the fit .
models are generated for three different inclination angles of the bulge @xmath263 and show that the best - fit @xmath0 and @xmath95 values are relatively insensitive to inclination angle .
we examine an additional 4 inclination angles for fixed bh mass values ( @xmath264 ) and fixed mass - to - light ratio ( @xmath265 ) .
we find that edge - on models ( @xmath101 ) give the smallest @xmath1 values , while models with the inclination of the large scale disk ( @xmath103 ) are strongly disfavored .
this suggests that the bulge must be nearly spherical but its rotation axis is significantly misaligned from the rotation axis of the disk .
however , it is well known from previous work @xcite that inclination of a spheroid is extremely difficult to determine via dynamical modelling , especially in the presence of a bar . this determination of the inclination of the bulge
should therefore be treated with caution .
a detailed examination of the kinematics in the inner 5  of ngc  4151 shows evidence for bar kinematics , which manifest as high values of @xmath231 , low or negative values of @xmath67 and @xmath66 values which are correlated with @xmath236 .
we hypothesize that the discrepancy between the @xmath1 obtained from including all the constraints and that obtained by only including the kinematic constraints is likely a consequence of two factors : ( a ) a stellar bar , which is not possible to model with our axisymmetric code , and ( b ) inconsistencies between the velocity dispersion values obtained from the low - spatial resolution long - slit data and the high spatial resolution nifs data over the same spatial range .
this hypothesis was tested by fitting only the nifs data .
when we fit only the kinematical data obtained with the nifs spectrograph ( neglecting all kinematical constraints beyond @xmath266 - i.e those data that show kinematical evidence for a bar and data that is discrepant with high - resolution nifs data ) , the non - negative optimization problem was able to find a much better fit to the nifs data and the same best - fit @xmath0 was obtained by both methods of computing @xmath1 .
the stellar - dynamical modeling carried out in this paper gives a best - fit value of @xmath267 ( obtained by averaging over the @xmath125 and @xmath127 values in the lower two rows of table  [ tab : robustbestfit ] ) which is consistent at the 1@xmath4 level with the reverberation mapping - based mass of @xmath268 ( 1@xmath4 errors ) obtained by using data from @xcite , but relying on a recently updated empirical calibration of the rm mass scale by @xcite .
the best fit m / l ratio @xmath269 ( 3@xmath4 error ) is consistent with the photometrically derived mass - to - light ratio of @xmath270 .    the new nifs data yields @xmath254 , the velocity dispersion within @xmath271 , which is @xmath1225% larger than previous values , while the new best estimate of @xmath272 17% smaller than but fully consistent with both our previous stellar dynamical upper limit @xcite and the dynamical estimate based on h@xmath260 line - emitting gas @xcite .
the larger value of @xmath245 and smaller value of @xmath0 help to bring ngc  4151 closer to the best - fit relation of @xcite .
the analysis in this paper demonstrates that biased estimates of bh masses can arise when an axisymmetric orbit superposition code is used to model a galaxy with a weak , but kinematically identifiable barred galaxy , possibly resulting in an over - estimate in @xmath0 if the m / l ratio is constrained primarily by data beyond the sphere - of - influence of the bh .
this confirms the prediction made from the recent analysis of @xmath229-body simulations of barred galaxies with bhs by @xcite , and points to the need for new dynamical modeling tools capable of modeling a stellar bar .
when such codes are applied to the existing sample of barred galaxies , our results suggest that it will likely enhance the discrepancy between barred and unbarred galaxies in the @xmath15 relation @xcite .
our modeling of this complex system suggests that the standard practice of fitting kinematical constraints over a large range of radii with a constant mass - to - light ratio ( for exceptions to this practice , see * ? ? ?
* ; * ? ? ?
* ) can bias the mass of the central bh and the derived mass - to - light ratio . to some extent
, the biases can be overcome by only considering the kinematics very close to the central bh , although this gives larger errors on the estimated solutions of @xmath0 and @xmath95 .
our tests of robustness demonstrate that the use of such high spatial resolution nuclear kinematical data in such modeling is very valuable .
however , the complications of the stellar dynamical bar imply that additional work is required in order to perform the crucial test of the rm mass scale calibration in this galaxy . only with improved non - axisymmetric modeling methods and/or bh mass measurements in other reverberation - mapped agns
will we ultimately be able to assess whether the growth history of the bhs that we see accreting at present is systematically different from those currently in quiescence .
we thank remco van den bosch for valuable input and for providing his symmetrization routine .
we thank gaelle dumas and eric emsellem for providing us with the sauron integral field kinematics .
we thank linda watson for providing copies of her nifs stellar velocity template spectra .
we thank the anonymous referee for suggestions that helped to improve our manuscript .
we are grateful to michele cappellari for making his analysis routines publicly available and for his continuing efforts in adding software features .
m.  valluri and jb were supported by national science foundation grant ast-0908346 , and m.  valluri acknowledges support from the university of michigan s elizabeth crosby award .
bmp and rwp are grateful to the national science foundation for support of this work through grant ast-1008882 to the ohio state university .
mb gratefully acknowledges support by the national science foundation under grant ast-1253702 .
m.  vestergaard acknowledges support from a freja fellowship granted by the dean of the faculty of natural sciences at the university of copenhagen and a marie curie international incoming fellowship .
the research leading to these results has received funding from the people programme ( marie curie actions ) of the european union s seventh framework programme fp7/2007 - 2013/ under rea grant agreement no .
300553 ( mv ) .
the dark cosmology centre is funded by the danish national research foundation .    based on observations obtained at the gemini observatory , which is operated by the association of universities for research in astronomy , inc .
, under a cooperative agreement with the nsf on behalf of the gemini partnership : the national science foundation ( united states ) , the science and technology facilities council ( united kingdom ) , the national research council ( canada ) , conicyt ( chile ) , the australian research council ( australia ) , ministrio da cincia e tecnologia ( brazil ) and ministerio de ciencia , tecnologa e innovacin productiva ( argentina ) .
this research was supported in part through computational resources and services provided by advanced research computing and the flux cluster at the university of michigan , ann arbor .
this work made use of data from the ohio state university bright spiral galaxy survey , which was funded by grants ast-9217716 and ast-9617006 from the united states national science foundation , with additional support from the ohio state university .
based on observations made with the nasa / esa hubble space telescope , obtained from the data archive at the space telescope institute .
stsci is operated by the association of universities for research in astronomy , inc . under the nasa contract
nas 5 - 26555 .
funding for the sdss and sdss - ii has been provided by the alfred p. sloan foundation , the participating institutions , the national science foundation , the u.s .
department of energy , the national aeronautics and space administration , the japanese monbukagakusho , the max planck society , and the higher education funding council for england .
the sdss is managed by the astrophysical research consortium for the participating institutions .
the participating institutions are the american museum of natural history , astrophysical institute potsdam , university of basel , university of cambridge , case western reserve university , university of chicago , drexel university , fermilab , the institute for advanced study , the japan participation group , johns hopkins university , the joint institute for nuclear astrophysics , the kavli institute for particle astrophysics and cosmology , the korean scientist group , the chinese academy of sciences ( lamost ) , los alamos national laboratory , the max - planck - institute for astronomy ( mpia ) , the max - planck - institute for astrophysics ( mpa ) , new mexico state university , ohio state university , university of pittsburgh , university of portsmouth , princeton university , the united states naval observatory , and the university of washington .
this research has made use of the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration .
this research has made use of the simbad database , operated at cds , strasbourg , france .
, m.  c. , denney , k.  d. , cackett , e.  m. , dietrich , m. , fogel , j.  k.  j. , ghosh , h. , horne , k. , kuehn , c. , minezaki , t. , onken , c.  a. , peterson , b.  m. , pogge , r.  w. , pronik , v.  i. , richstone , d.  o. , sergeev , s.  g. , vestergaard , m. , walker , m.  g. , & yoshii , y. 2006 , , 651 , 775    , m.  c. , walsh , j.  l. , barth , a.  j. , baliber , n. , bennert , v.  n. , canalizo , g. , filippenko , a.  v. , ganeshalingam , m. , gates , e.  l. , greene , j.  e. , hidas , m.  g. , hiner , k.  d. , lee , n. , li , w. , malkan , m.  a. , minezaki , t. , sakata , y. , serduke , f.  j.  d. , silverman , j.  m. , steele , t.  n. , stern , d. , street , r.  a. , thornton , c.  e. , treu , t. , wang , x. , woo , j .- h . , & yoshii , y. 2009 , , 705 , 199    , m.  c. , walsh , j.  l. , barth , a.  j. , yoshii , y. , woo , j .- h . ,
wang , x. , treu , t. , thornton , c.  e. , street , r.  a. , steele , t.  n. , silverman , j.  m. , serduke , f.  j.  d. , sakata , y. , minezaki , t. , malkan , m.  a. , li , w. , lee , n. , hiner , k.  d. , hidas , m.  g. , greene , j.  e. , gates , e.  l. , ganeshalingam , m. , filippenko , a.  v. , canalizo , g. , bennert , v.  n. , & baliber , n. 2010 , , 716 , 993                                  , g. , de vaucouleurs , a. , corwin , jr .
, h.  g. , buta , r.  j. , paturel , g. , & fouqu , p. 1991
, third reference catalogue of bright galaxies .
volume i : explanations and references . volume ii : data for galaxies between 0@xmath273 and 12@xmath273 .
volume iii : data for galaxies between 12@xmath273 and 24@xmath273 .
( springer - verlag : new york )            , p.  b. , frogel , j.  a. , pogge , r.  w. , quillen , a.  c. , berlind , a.  a. , davies , r.  l. , depoy , d.  l. , gilbert , k.  m. , houdashelt , m.  l. , kuchinski , l.  e. , ramrez , s.  v. , sellgren , k. , stutz , a. , terndrup , d.  m. , & tiede , g.  p. 2002
, , 143 , 73          , k. , richstone , d. , tremaine , s. , lauer , t.  r. , bender , r. , bower , g. , dressler , a. , faber , s.  m. , filippenko , a.  v. , green , r. , grillmair , c. , ho , l.  c. , kormendy , j. , magorrian , j. , & pinkney , j. 2003 , , 583 , 92                          , c.  l. , merritt , d. , olling , r. , valluri , m. , bender , r. , bower , g. , danks , a. , gull , t. , hutchings , j. , kaiser , m.  e. , maran , s. , weistrop , d. , woodgate , b. , malumuth , e. , nelson , c. , plait , p. , & lindler , d. 2001 , , 550 , 668          , p .- y . , cappellari , m. , emsellem , e. , bournaud , f. , michel - dansac , l. , alatalo , k. , blitz , l. , bois , m. , bureau , m. , davies , r.  l. , davis , t.  a. , de zeeuw , p.  t. , duc , p .- a . , khochfar , s. , krajnovi , d. , kuntschner , h. , morganti , r. , mcdermid , r.  m. , naab , t. , oosterloo , t. , sarzi , m. , scott , n. , serra , p. , weijmans , a .- m . , & young , l.  m. 2012 , , 424 , 1495          , p.
j. , hart , j. , conroy , p.  g. , pfitzner , m.  l. , bloxham , g.  j. , jones , d.  j. , downing , m.  d. , dawson , m. , young , p. , jarnyk , m. , & van harmelen , j. 2003 , in society of photo - optical instrumentation engineers ( spie ) conference series , vol .
4841 , society of photo - optical instrumentation engineers ( spie ) conference series , ed .
m.  iye & a.  f.  m. moorwood , 15811591              , c.  a. , valluri , m. , peterson , b.  m. , pogge , r.  w. , bentz , m.  c. , ferrarese , l. , vestergaard , m. , crenshaw , d.  m. , sergeev , s.  g. , mchardy , i.  m. , merritt , d. , bower , g.  a. , heckman , t.  m. , & wandel , a. 2007 , , 670 , 105 ( paper  i )                              , s. , gebhardt , k. , bender , r. , bower , g. , dressler , a. , faber , s.  m. , filippenko , a.  v. , green , r. , grillmair , c. , ho , l.  c. , kormendy , j. , lauer , t.  r. , magorrian , j. , pinkney , j. , & richstone , d. 2002 , , 574 , 740                      , j .- h . ,
treu , t. , barth , a.  j. , wright , s.  a. , walsh , j.  l. , bentz , m.  c. , martini , p. , bennert , v.  n. , canalizo , g. , filippenko , a.  v. , gates , e. , greene , j. , li , w. , malkan , m.  a. , stern , d. , & minezaki , t. 2010 , , 716 , 269