in the present - day view , mergers are a common and important way of building up massive elliptical galaxies . with the presence of supermassive black holes in the centers of most galaxies ( @xcite ;
@xcite ; @xcite ) , merger processes are thought to have left signatures on the central structure of the remnant galaxies .
@xcite suggest that a binary black hole that is formed by the merger of two galaxies scours the stars from the center of the newly created system as the binary shrinks .
the energy liberated by the hardening of the binary evacuates the central part through the ejection of surrounding stars , causing less light in the center
.    _ hubble space telescope _
( _ hst _ ) observations of galaxy nuclei see two kinds of behavior in the central light profiles of ellipticals , traditionally classified as `` power - law '' and `` core '' galaxies @xcite .
the power - law profiles took their name from their approximation by a single power law at small radii ( @xmath6 or 20 ) .
more modern interpretations have emphasized that these profiles can be better understood as the inward continuation of the galaxy s overall srsic profile , usually modified by an additional , nuclear - scale stellar component @xcite .
the core galaxy , in contrast , displays a surface brightness profile with a distinct break from a steep outer slope to a shallower inner cusp .
core profiles mainly occur in very luminous elliptical galaxies and are considered the result of dissipationless mergers of two galaxies that have central black holes .
much numerical work has been done to explore the binary black hole idea .
the work of @xcite suggests that the stellar mass ejected from the core is comparable to the mass of the central black hole @xmath0 .
@xcite find that core radius roughly doubles after each major merger .
@xcite , using isothermal king models , suggest that the core size should scale with the final @xmath0 .
simulations performed by @xcite show for the first time that merging two galaxies with steep cusps can result in a merger remnant with a shallow power - law cusp in the inner part ( core ) .
@xcite is the first to quantify the magnitude of the mass deficit ( or the mass ejected by the binary ) with respect to the galaxy s merger history . by following the binary evolution up to the stalling radius , he formulates the mass deficit to be @xmath7@xmath0  with @xmath8 is the number of major merger events .
an investigation into the post - coalescence stage of the binary suggests that the mass deficit could still be larger .
this is because the coalescing binary emits gravitational waves which impart a kick on the merged black hole and eject it from the center .
due to dynamical friction , the resulting oscillations of the black hole would gradually damp , returning it to the center of the galaxy and ejecting more stars from the core .
this process could increase the mass deficit existing before the binary coalescence by up to 5@xmath0  for a single merger @xcite .
several publications in the past have performed function fitting to galaxy profiles in order to quantify the core parameters such as the core size , the deficit in light or mass , and to analyse the correlation between these core properties and @xmath0  or the global properties of the galaxies ( @xcite ; @xcite ; @xcite , kormendy et al .
( 2009 , hereafter kfcb09 ) , kormendy & bender ( 2009 , hereafter kb09 ) , @xcite ; @xcite ; see also @xcite ) .
these studies are based on different samples , datasets , approaches , parametrizations , and assumptions , which apparently lead to different , partly systematic , core quantities .
the most important difference is the amount of central mass deficit with respect to the black hole mass that can be expected in core galaxies . @xcite and @xcite find an average mass deficit of around twice the black hole mass , @xcite derive central mass deficits of about 0.54 @xmath0  while find their mean mass deficit to be 10.5 @xmath0 . aside from that ,
each of these studies has a drawback in at least one of three aspects crucial to achieve high accuracy .
the first is the radial range of the surface brightness profiles .
@xcite , @xcite , and @xcite fit galaxy profiles that extend typically to only 1020 .
moreover , it appears that some galaxy profiles in @xcite , which are taken from lauer et  al .
( 2005 , hereafter l05 ) , suffer from oversubtraction of the sky that biases the core - srsic fits ( see section [ previousfits ] ) .
a full galaxy profile is necessary to properly measure the outer profile , to optimally remove the sky background , and thus to provide better constraints on the light and mass deficit .
the second is the determination of the stellar mass - to - light ratio ( @xmath9 ) to convert the light deficit to the mass deficit . in the past studies ,
either different values are assigned to the galaxies according to their luminosities ( @xcite , ) or a constant value is assumed for all the galaxies in the sample ( @xcite , @xcite ) .
this of course introduces additional uncertainties and biases in deriving the mass deficits .
the last important aspect lies in the black hole masses , which are mostly derived from outdated black hole - bulge relations , especially the ones between @xmath0  and the velocity dispersion of the host bulge ( @xmath0@xmath3relation ) .
the majority of core galaxies are bright massive elliptical galaxies . for this class of galaxies
@xcite find that the @xmath10-@xmath11  relation predicts a higher density of massive black holes ( @xmath12 @xmath13 ) compared to the @xmath0@xmath3  relation .
furthermore , the @xmath0@xmath3  relation used in these papers also fails to predict the existence of the largest black holes ( @xmath14@xmath13 ) that are found in distant quasars @xcite in the sample of local galaxies .
the recent finding of black holes of 10 billion solar masses @xcite , which might be the missing link between the local black holes and the massive distant quasars , further supports the indications that the previous @xmath0@xmath3  relations for massive bright early type galaxies need to be reevaluated ( e.g. , * ? ? ?
in addition , most previous studies have used @xmath0  measurements ( or @xmath0@xmath3relations based on @xmath0  measurements ) which did not include dark - matter halos in the modeling .
values typically increase by @xmath15% @xcite or 30% ( * ? ? ?
* hereafter paper i ) when dark matter is included in the modeling , but the increase can be as large as a factor of 6 when the black hole s sphere of influence is not well resolved @xcite . in the particular sample that we use in this paper , which includes bcgs , we find that the measured black hole mass is on average more than 3 times higher than that predicted by the @xmath0@xmath3  relation ( e.g. @xcite , @xcite , @xcite .    in this paper , we combine improved core - property measurements with new measurements and present the analysis of possible relations between the core and black hole mass or the galaxy properties .
we use only directly measured black hole masses and consequently avoiding scatter and possible biases inherent in using proxy measurements .
we use dynamically - determined , individual stellar mass - to - light ratios ( @xmath2 ) to estimate the mass deficit , instead of relying on a scaling function or a constant value .
we also use full galaxy profiles in the fitting and the quantification of the core parameters .
our analysis is based on a sample of 20 galaxies ; seven of these galaxies are taken from , in which the black hole masses and @xmath2  are derived taking into account the presence of dark matter halos . for completeness , we add a further 3 galaxies to the sample ; they are not included in the main sample due to the lack of reliable literature sources for the black hole masses . in order to derive the mass deficit ,
we assume that the original ( pre - scouring ) profile is the extrapolated srsic  profile that fits the light distribution at large radii . in section  [ themethod ] ,
we describe the method that we use to identify the core galaxies and to fit the observed surface brightness profiles . in section  [ corefitting ]
, we present the best - fitting models for each galaxy ; cases where galaxies are best fit with the addition of an outer envelope are discussed in section  [ multi - component - fits ] .
we compare the fitting results with previous studies in section  [ previousfits ] .
section  [ coredeficit ] presents the details of how the luminosity and mass deficit in the core are computed based on the best - fitting models .
the core - related quantities derived in section  [ corefitting ] and [ coredeficit ] are compared with the black hole mass , velocity dispersion and galaxy luminosity in section [ bhcore ] .
the last section discusses and summarizes the results .
all of the photometric profiles examined in this work are quite extended , reaching out to at least 9 kpc .
therefore , we base our analysis on a function that is designed to fit not just the inner or the outer part of galaxies , but rather the galaxy as a whole . to identify core galaxies , we follow the criteria in @xcite and @xcite , who define the core as `` a downward deviation from the inward extrapolation of the outer ( srsic ) profile '' .
this involves fitting the galaxy in question with both the srsic  and the core - srsicfunctions .
the srsic  profile @xcite is written as : @xmath16 \right\}. \label{sersiceq}\ ] ] @xmath17 is the intensity at @xmath18 , the projected half - light radius .
@xmath19 is called the srsic  index which describes the shape or curvature of the light profile .
the quantity @xmath20 is a function of @xmath19 , defined in such a way that @xmath18 encloses half of the total luminosity .
we approximate @xmath21 by using the asymptotic expansion of @xcite , eq . 18 for @xmath22 .
many other approximation formulas are available and are summarized in @xcite . for @xmath23 ,
the srsic  function reduces to an exponential function and for @xmath24 , it becomes the de vaucouleurs profile @xcite .
the core - srsic  function introduced by @xcite and @xcite is expressed as : @xmath25^{\gamma/\alpha}\,{\rm exp}\left[-b_n\left(\frac{r^\alpha+r_b^\alpha}{r_e^\alpha}\right)^{1/n\alpha}\right ] \label{coresersiceq1}\ ] ] with @xmath26 .
\label{coresersiceq2}\ ] ] this profile consists of a srsic  profile in the outer part , specified by the projected half - light radius @xmath18 and the srsic  index @xmath19 , and a power - law profile in the inner part with a slope of @xmath27 . the change from one to another regime occurs at the break radius @xmath28 and the sharpness of the transition is specified by the parameter @xmath29 .
higher values of @xmath29 mean sharper transitions .
@xmath20 is approximated in the same way as for the srsic  function .
we summarize four criteria in @xcite that we use to identify core galaxies below . in section [ corefitting ]
we examine each criterion .    1
.   a characteristic pattern should be visible in the residuals when fitting a srsic  profile to an idealized core galaxy .
this serves as a qualitative evidence .
2 .   a core - srsic  function gives a significantly better fit than a srsicmodel .
quantitatively , it is expected that the reduced @xmath30 value for the srsic  fit @xmath31 is larger than twice that of the core - srsicfit @xmath32 .
3 .   to avoid an ambiguous core detection , the potential core must be well - resolved by the data .
quantitatively , @xmath28 should be @xmath33 where @xmath34 is the second innermost datapoint in the observed profile .
4 .   the power - law slope has to be consistently less than the slope of the srsic  fit inside @xmath28 .
we decided to adopt an additional criterion based on resolution limits .
trujillo et al .  did not include such a limit , partly because many of their profiles did not probe radii smaller than the _ hst _ resolution limit due to dust , etc . ; however , they did note that one of their `` possible core '' galaxies ( ngc 3613 ) had a break radius smaller than the resolution limit of 0.16  originally suggested by @xcite .
the latter limit is @xmath35 times the fwhm of most of the pc data used by @xcite ; generalizing this , we suggest the following additional criterion :    1
.   the potential core should have @xmath36 the fwhm of the imaging data used for the innermost part of the profile . for most of the galaxies considered in this paper , which use optical wfpc2 or acs / wfc data
, this corresponds to @xmath37 ; for two galaxies using nicmos2 data ( ngc 4261 and ngc 4374 ) , the limit would be @xmath38 , and for ngc 1550 , where we use sinfoni data with a fwhm of 0.17 , the limit would be 0.34 .
we construct a one - dimensional light profile along the circularized radius ( mean profile ) , instead of along the semi major radius as is commonly done .
the circularized radius is defined as @xmath39 where @xmath40 and @xmath41 are the semi major and semi minor axes radius of the isophote .
this way , we take into account the variation of ellipticity with radius , giving a fair representation of the isophotes .
the use of this mean profile is consistent with the implementation of the psf ( point spread function ) correction ( section [ psfconv ] ) .
the fitting is done using a one - dimensional profile fitting software written by p. erwin .
an initial model is generated based on user - given starting values .
it is compared to the observed profile and the @xmath30 is calculated .
the best - fitting solution is achieved by minimizing the @xmath30 of the fit , using the levenberg - marquardt algorithm @xcite .
the user provides a starting value and its range of allowed values for each free parameter of the model .
there are three free parameters in the srsic  model , i.e. @xmath19 , @xmath18 , @xmath17 and there are six free parameters in the core - srsic  model , i.e. @xmath19 , @xmath28 , @xmath18 , @xmath42 , @xmath29 , @xmath27 .
since we are mainly interested in the inner region of the galaxy , the blurring of the central light profile by the psf becomes important and has to be taken into account . for several galaxies ,
the photometric profiles are derived from deconvolved images ( see table  [ tab : coresample ] ) .
we fit these profiles as described in the previous paragraph . for the other galaxies we incorporate a psf correction routine , described in the next section , into the fitting program such that the initial model is convolved first before being compared to the observed light profile . in this case
, the minimised @xmath30 represents the difference between the convolved model and the data .
the psf convolution is done at every iteration .
equal to 1/40 of the half - light radius of the srsic model .
dashed lines and diamonds show the case of a circular srsic  image ; the solid lines and rectangles are for an elliptical srsic image ( ellipticity = @xmath43 ) .
the red lines are the unconvolved srsic  models .
the black symbols are the results of one - dimensional convolution ( equation [ equ5]-[equ7 ] ) and the blue lines are the radial profiles of the convolved images , derived using ellipse - fitting . *
bottom panel : * the flux ratio or the difference in surface brightness ( @xmath44 ) between the one - dimensionally convolved srsic  profile ( @xmath45 or @xmath46 ) and the radial profile of the convolved srsic  image ( @xmath47 or @xmath48 ) for both the circular and elliptical cases .
+ ]    for galaxies whose light profiles are not seeing - corrected , we take the psf into account in the fitting by means of convolution .
we use the psf corresponding to the galaxy image used to derive the innermost galaxy isophotes .
for the galaxies from , the psf image is the same one which was used in the deprojection .
for the other galaxies , the psf image is generated using tinytim by supplying the appropriate filter and camera . from the image , we generate a circularized profile of the psf by averaging its radial profile at different azimuthal angles ( done using midas ) .
we further normalize the psf profile such that the total flux is one .
the implementation of the convolution follows @xcite ( their appendix a ) .
below , we review the method .
the psf convolution is carried out by multiplying the galaxy surface brightness profile and the psf in the fourier space . on the sky ,
the galaxy surface brightness is a function of a two - dimensional position * r * and so is the psf .
if @xmath49 and @xmath50 denote the surface brightness and the psf distribution , then the convolved surface brightness profile becomes : @xmath51 where @xmath52    since we use a one - dimensional , circularized surface brightness profile , i or p are a function of only the radius @xmath53 . in this case , @xmath54 where @xmath55 is the bessel function of order zero .
then , equation [ eqq ] becomes : @xmath56 with @xmath57 @xmath58    in this way , although the surface - brightness and the psf profiles are one dimensional , the convolution is implemented such that the light coming from different angles is also taken into account .
this is important because the light is scattered , due to the seeing , not only in the radial but also in the azimuthal direction .
equations [ equ5]-[equ7 ] are also called the hankel transform .
we compute the convolved surface brightness profile numerically , making use of integration routines in the gnu scientific library .    fig .
[ testconv ] shows how well this convolution works in the case of circular ( upper profiles ) and non - circular ( lower profiles ) images .
for the former , we generate a galaxy image which follows a srsic  law with @xmath24 and zero ellipticity everywhere using the ` makeimage ` program from the two - dimensional image - fitting package imfit ( erwin , in preparation ) . for the psf
, we create a circular gaussian image with @xmath3 = 1/40 of the half - light radius of the galaxy model .
we extract the radial profile of both images and use these to perform the one - dimensional convolution as described above . as an independent check
, we separately convolve the galaxy image with the gaussian image in imfit . the radial profile of this convolved image ( @xmath47 )
is compared to the result of the 1d convolution ( @xmath45 ) .
the ratio is shown in the lower panel of the figure by the dashed line .
this process is then repeated for a srsic image with an ellipticity of 0.4 .
the circularized profile of the convolved elliptical image is shown by the solid blue line and this is compared to the 1d convolution result ( the rectangles ) ; the ratio of both is displayed as the solid line in the lower panel .
the method works well in the case of circular galaxy and shows a deviation of less than 0.02 mag arcsec@xmath59 for an elliptical image ; we have repeated this exercise for different srsic indices , resulting in similar conclusions . since the ellipticity of 0.4 is larger than the ones measured in the majority of galaxies in our sample , we expect smaller errors in the actual galaxy case .
the list of galaxies examined in this section is given in table  [ tab : coresample ] , along with the data sources .
we compiled this list taking into account all galaxies which had both dynamically measured central black hole mass and evidence for cores , but in the end retained only those with clear evidence for a core _ and _ good enough data for measuring the core properties .
the recent compilation of @xcite lists 9 more galaxies which we do not include : a1836-bcg , a3565-bcg ( ic 4296 ) , ngc 524 ngc3607 , ngc4473 , ngc5077 , ngc5128 , ngc5576 , ngc7052 .
strong nuclear dust in a1836-bcg , ic 4296 , ngc3607 , ngc 5128 and ngc7052 makes it very difficult to robustly compute the core properties , so these galaxies are not included .
ngc  4473 has been shown not to be a core galaxy by @xcite , , and @xcite , while the same is true for ngc  5576 @xcite ; see also @xcite .
in addition , @xcite were only able to classify ngc  5077 as a `` possible core '' galaxy , so we exclude it as well . ngc  524 is a particularly problematic case , since there is evidence for very complex morphology
e.g. , @xcite fit a ground - based image using a srsic bulge , an exponential disk , and two `` lens '' components .
this means that a simple core - srsic fit  or even core - srsic + outer envelope , as we use for some galaxies  is probably not appropriate .
in addition , the best available high - resolution data is _ hst_-nicmos2 , and the fit to this data by @xcite suggests a core break radius @xmath60 , which is smaller than our adopted resolution limit for nicmos2 data ( see section [ identifycore ] ) .
( wfpc2 images are not usable due to the extensive circumnuclear dust lanes . )
we reluctantly exclude this galaxy as well .
ngc1374 represents a case where the core detection is ambiguous , so it is discussed independently in a separate section .
this galaxy was classified as a non - core galaxy in @xcite using a limited dataset from , and also by @xcite using a more extended acs - wfc profile from @xcite ; in section  [ n1374 ] we show the fit using our surface brightness profile which is slightly more extended in radius and which does not suffer the possible sky - subtraction uncertainties of the acs - only turner et al .
table  [ tab : csparams ] summarizes the list of galaxies for which the presence of core is confirmed .
the best - fitting core - srsic  parameters are also given in the same table .
the best - fit srsic  and core - srsic  functions of these galaxies are shown in fig .
[ fittingplots ] and the criteria which lead to the core identification are discussed below .      when possible , we use all the available datapoints in the fitting . for all galaxies ,
the core - srsic  function gives a much better fit to the whole profile than the srsic .
in several galaxies , however , we find that the fit residuals become rather large at several outermost datapoints .
since these datapoints do not seem to be well - described by a srsicfunction , they are excluded from the fit .
we keep only the datapoints where the fit residual is less than 0.15 mag .
the deviation from the srsic model at large radii by more than 0.15 mag may be caused by imperfect sky subtraction or it may represent a real departure from the srsic  profile , for example when where an outer halo is present ; see section  [ multi - component - fits ] for details on fits using extra components .
the fit is then repeated using only the remaining datapoints , after excluding the outliers .
the fit residuals , their rms and the @xmath30 values are presented in fig .
[ fittingplots ] . in this figure ,
the surface brightness profiles shown are not corrected for galactic extinction .    in all the galaxies in fig .
[ fittingplots ] we see a pattern in the residuals of the srsic  fit that qualitatively suggests the presence of the core
compare to fig .  3 in @xcite .
this residual pattern disappears when the galaxy is fitted with the core - srsic  profile .
so the first criterion in section [ identifycore ] is secured . in each case , the rms of the residuals for the best srsic  fit is much larger than that for the core - srsic .
the @xmath30 values for each galaxy are calculated by adopting the residual rms of the best core - srsic  fit as the error on the surface brightness measurements . as discussed in ,
this estimate for the real error is reasonable as long as the rms is of order of the profile measurement error , i.e. , a few hundredths of a mag arcsec@xmath59 , which is indeed the case for all the galaxies that we examine . using this error approximation , the @xmath30 of the core - srsic  fit
is artificially set to 1.0 .
note that the rms and @xmath30 values are dependent on the radial extent and the radial sampling ( logarithmic sampling is used here ) .
we find that the ( reduced ) @xmath30 value for each srsic  fit is much more than twice that of the core - srsic .
the average ratio of @xmath61 is around 49 , with a minimum of 6.8 , easily fulfilling the second criterion . to assess the third criterion
, one can compare the position of the second innermost datapoint to the position of the dotted vertical line which marks the break radius .
it is thus clear that we have cases of resolved cores .
the fourth point is also satisfied , since inside the break radius the core - srsic  function shows a consistent deficit in light compared to the extrapolated srsicprofile ( shown by the blue line ) . for the cases where we exclude some datapoints ,
the core is already well identified when fitting all the datapoints , i.e. fulfilling the four criteria simultaneously .
limiting the fitting range only improves the distinction .
lastly , the break radius for each galaxy is larger than twice the fwhm of the imaging data used to construct the innermost part of the surface brightness profile , described more quantitatively in section [ identifycore ] .    among the 23 galaxies ,
several exhibit more complex structures . inside @xmath62 ,
the surface - brightness profile of ic1459 seems to show a small excess of light .
we model this nuclear excess as an extra gaussian component in the center on top of the core - srsic  function , which after psf convolution provides a good fit to the data .
( in the cases of ngc1407 and ngc4552 , it is not clear that their inner excesses can be modeled with a gaussian , so we exclude the data at @xmath63 for ngc1407 and @xmath64 for ngc4552 . ) for six galaxies ( including ngc1407 ) , we find that significantly better fits are achieved by including one or two extra srsic  or exponential components , representing outer halos .
this approach is discussed in more detail in section  [ multi - component - fits ] .
the parameter values of the best - fit core - srsic  model for each galaxy are provided in table  [ tab : csparams ] .
the uncertainties of the core - srsicparameters are determined through a monte carlo simulation ( see next section ) .
the srsic  indices for these galaxies are almost all larger than 4 , something which is common for giant ellipticals ( e.g. , @xcite , ) .
the exceptions are four of the galaxies which we model with multiple components ( see section  [ multi - component - fits ] ) .
the best - fit @xmath29 parameter for ngc7768 is found to be 16.6 .
this large value indicates a very sharp transition between the srsic  and the power - law profile .
when the fit is repeated with @xmath65 , which is a good approximation for @xmath66 @xcite , the best - fit values for the other five parameters hardly change .
since for large @xmath29 ( @xmath67 10 ) the profile becomes insensitive to the exact value of @xmath29 , we prefer to use the latter fit ( fixing @xmath29 to 100 ) to the light profile of ngc7768 .
this choice prevents wildly varying values of @xmath29 that are not physically meaningful in the monte carlo simulations .    for six galaxies whose
surface brightness profiles are taken from ( i.e. , ngc4261 , ngc4374 , ngc4472 , ngc4486 , ngc4552 , and ngc4649 ) , the best - fitting srsic  model of is also overplotted on the right - hand panels of fig .
[ fittingplots ] .
our inward - extrapolated srsic  components and their srsic  fits do not coincide , reflecting the differences in the way the core is defined .
note that use major - axis profiles while we use circularized profiles , resulting in a seemingly poorer fit of kfcb s srsic  models in fig .
[ fittingplots ] . for a more thorough comparison between our and their fitting procedure , see appendix [ comparisonkfcb09 ] .
ngc3091 and ngc7619 show small values in @xmath29 , i.e. @xmath68 .
@xcite produce a best - fitting core - srsic  model of ngc7619 by keeping @xmath29 constant at 5 , which enforces a rather sharp transition from srsic  to power - law profiles . doing this exercise using our circularized profile of ngc7619 results in a smaller @xmath19 ( 7.71 ) , larger @xmath28 ( 0.73 ) and @xmath27 ( 0.47 ) and a larger rms ( 0.033 ) .
it is worth emphasizing that @xcite s profile for ngc7619 is not as extended as ours and stops short of @xmath69 .
ngc1374 is the one galaxy in the sample of with conflicting arguments about the existence of a core .
classify ngc1374 as a core galaxy , based on a nuker profile fit to their _ hst _ imaging data .
@xcite , however , do not confirm the core in this galaxy using the same data , and @xcite come to the same no - core conclusion using a more extended acs - wfc  based profile from @xcite .
in contrast , @xcite argues for the presence of the core .    from fitting core - srsic  models to our data , we find no solid evidence that a core is present in this galaxy .
we tested different fitting ranges by varying the maximum radii of the fitted datapoint , but there is no strong visibility of a core in any case .
in the top row of fig .
[ fitting1374 ] , we show the srsic  and core - srsic  fits to the whole profile and to a limited profile . for the latter , we limit the fitting range to @xmath70 ( for a fair comparison to the fitting using data , see below ) . judging by eye , there is no significant difference in fit quality between the srsic  and the core - srsicprofiles , although the overall residual rms and @xmath30 of the core - srsic  are smaller and @xmath71 .
when fitting the whole range , the srsic  profile actually provides a better fit in the center and instead of seeing a core , the core - srsic  fit indicates an extra light on top of the srsic  component ( this might be consistent with the extra nuclear component used by @xcite in their fits to this galaxy ) . there is a slight improvement and a vague indication of a core when we fit only up to 7 .
however , if we base our core classification on @xcite s criteria , this fit fails at the second criterion .
our next attempt is to fit the data of to confirm or reject their finding .
their profile has been corrected for the psf effect , so we directly fit it with the models without convolution . as shown in fig .
[ fitting1374 ] , we detect a core in this galaxy , albeit weak , by discarding the two outermost datapoints ( @xmath72 ) . extending the fitting range of data
does not give the same results , i.e the existence of the core becomes inconclusive and the core disappears altogether when all the datapoints are included in the fit .
fit their light profile with a nuker model and it appears that they also made the compromise of neglecting the 34 last datapoints in the fit , in order to detect the core .
moreover , our best - fitting break radius is smaller than 0.1 , less than the inner resolution limit that we adopt as an additional criterion of core detection in section [ identifycore ] .
one might conclude that our data lack spatial resolution and that this situation is aggravated by the psf convolution in the fitting process .
it is , however , important to note that the core disappears when extended datapoints are included in either our profile or the profile .
graham et al .
( 2003 ) pointed out that fitting a nuker model to a limited light profile , such as the case of ngc1374 in , could lead to a false conclusion that a core exists . having no definite evidence of a depleted core in ngc1374 , we exclude it from the core sample .
[ cols= " < , < , < , < , < , < , < , < " , ]     taking the break radius ( @xmath28 ) as a measure of the core size , we examine its relationships with galaxy velocity dispersion , luminosity and the black hole mass .
these are displayed in fig .
[ rbplots ] .
the fit to the datapoints ( also for all the other fits in this section ) is based on a bayesian method with gaussian errors and intrinsic scatter @xmath73 , described in @xcite , and the results are shown by the black lines . for each diagram
we perform non - symmetrical regression using either of the parameter to predict the other .
we further compute the bisector of the two fitting results ( with flipped @xmath74 and @xmath75-axes ) to give a symmetrical representation of the correlation .
the coefficients of the fits involving @xmath28 are shown in table  [ tab : core - relations ] .
all the fits are linear fits to the logarithm of the relevant quantities expressed in the following form : @xmath76 where @xmath75 is @xmath77 ( in kpc ) and @xmath74 is @xmath78 , @xmath79 , or @xmath80 . for
the `` direct '' fits @xmath81 , @xmath79 , and @xmath82 are the independent variables in the fit , while in the `` reverse '' fits @xmath28 is the independent variable .
the errors on the coefficient values are taken as the rms of the a - posteriori parameter probability distribution .
we also list the spearman correlation coefficients and the corresponding probabilities for each relation in table  [ tab : core - relations ] .
@xcite analysed the same relationships using bisector fits . in the @xmath28@xmath81 diagram ,
our fit predicts a slightly higher @xmath28 for a given velocity dispersion , though the fits are compatible within the errors . in the @xmath28@xmath83 diagram ,
our fit coincides with their fit for high - luminosity galaxies ( @xmath84 ) . in the case of the @xmath28@xmath0 relation ,
graham & dullo performed fits using both measured black hole masses ( their equation  12 ) and proxy black hole masses estimated from the @xmath85@xmath3 and @xmath0@xmath86 relations ( their equation  13 and 14 ) ; only the former is actually plotted ( their fig .
20 ) . since our fits
use only measured black hole masses , we compare our results with their equation  12 ( lower right panel of figure  6 ) . for this relation ,
our fit suggests larger @xmath0 for the same value of break radius ( or , conversely , smaller break radius for a given @xmath0 ) .
we note that the fit in dullo & graham s equation  12 is based on only 7 galaxies , all but one of which have black hole masses @xmath87 , versus the 20 galaxies used in our fit , with 14 having @xmath88 .
the measured rms scatter is almost identical for the @xmath28@xmath81 and @xmath28@xmath0 fits ( 0.24 dex versus 0.23 dex ) .
the tightest relationship turns out to be between @xmath28 and @xmath83 , where the rms is 0.20 dex .
the correlation coefficients indicate that all three relations are strong and statistically significant , with the strongest being that between @xmath89 and @xmath0 ( @xmath90 , @xmath91 ) . both @xcite and @xcite
use a larger number of datapoints than ours in examining the @xmath89@xmath92 and @xmath89@xmath93 relationships , and they find opposing results as to which parameter ( @xmath92 or @xmath93 ) is better correlated with @xmath89 ; our results suggest that @xmath93 is a better predictor of @xmath89 than @xmath92 .
( note that the @xmath89@xmath93 correlation uses the _ total
_ galaxy luminosity ; if for the multi - component galaxies we use the luminosity of just the core - srsic  component instead , the correlation is weaker and only marginally significant , with @xmath94 and @xmath95 . ) despite the moderately large scatter , the @xmath89@xmath0 relation is apparently the strongest one among the three , suggesting that @xmath85 is the most relevant parameter for determining @xmath89 .
[ deficitplots ] explores relationships between the light or mass deficit in the core and the black hole mass or the stellar velocity dispersion .
the upper left panel shows the plot of the light deficit in the @xmath96-band versus the stellar velocity dispersion calculated within the effective radius @xmath81 .
the correlation between stellar velocity dispersion of host galaxies and @xmath97 determined in is overplotted in red .
note that the velocity dispersion used to construct the latter correlation is defined differently : it is based on the velocity dispersion averaged within a slit aperture with a length of two half - light radii .
although there is probably no significant difference between these two methods of determining @xmath81 , this means that the individual datapoints can not be strictly compared to the relation . for seven galaxies that overlap with s sample ,
we indicate their values of light deficit using gray diamonds . on average ,
our values are more negative ( i.e. more luminous ) than those of by 0.30 mag , and the rms of the difference is 1.15 mag . in the upper right panel of fig .
[ deficitplots ] , the plot between @xmath97 and the black hole masses is shown .
similarly , we overplot the correlation between these two quantities , as derived in , in red .
coefficients of the fits plotted in fig .
[ deficitplots ] are listed in table  [ tab : mdef - relations ] .
the meaning of the coefficients is the same as in equation  [ eqn : basic ] , except now @xmath75 corresponds to either @xmath98  or @xmath99 . in the direct fits , the independent variables are either @xmath100 or @xmath4 .
compared to , our @xmath97@xmath81 correlation is steeper and our @xmath97@xmath0  exhibits a shift towards a larger @xmath0  at a given value of light deficit .
this highlights the difference between estimates based on the @xmath0@xmath3  relation and the more reliable dynamical estimates of @xmath0  that we use .
the scatter is notably higher in @xmath97@xmath0  diagram , opposite to what is found in .
the most interesting relation from a theoretical point of view is that between the _ mass _ deficit and the black hole mass , since this is the quantity which simulations predict ( see references in the introduction ) .
the lower panels of fig .  [ deficitplots ] show this relation , and also the relation between mass deficit and velocity dispersion ; the coefficients of the fits are in the bottom half of table  [ tab : mdef - relations ] . the large intrinsic scatter ( and the large rms values ) that we derive indicate that the linear fit is not representing the data well .    of the seven correlations we have examined ( figs . [ rbplots ] and [ deficitplots ] ) ,
those between @xmath28 and @xmath0 and between @xmath97 and @xmath81 appear to be the strongest ( @xmath101 and @xmath102 , respectively ) .
moreover , the rms scatter in @xmath103 for the fit between @xmath0 and @xmath28 is only 0.30 dex ( it is 0.35 dex for the correlation with @xmath104 and 0.33 dex for that with the light deficit ) .
this is comparable to or even lower than the typical rms scatter for the @xmath0@xmath3 relation ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
this suggests that measuring @xmath89 can allow one to estimate @xmath0 ( in core galaxies ) with a precision similar to that obtained using the @xmath85@xmath3 relation , without the need for spectroscopy .
we give here below the corresponding equation ( derived from the fits in table  [ tab : core - relations ] ) : @xmath105    the ratio between @xmath104  and @xmath0  seems to vary strongly , but is biased toward the low end of the range .
most galaxies have a mass deficit between 1 and 5 @xmath0 , and the median value is 2.2 @xmath0 . in three galaxies ( ic1459 , ngc3842 and ngc5813 ) ,
the mass deficit is just below the black hole mass .
ngc3608 has a remarkably low mass deficit : only 0.21 @xmath0 .
the largest mass deficit that we find is just below 40 @xmath0 , occurring in ngc1399 .
this unusually large value should perhaps be taken with a grain of salt , however , given the large uncertainty in the black hole mass of ngc1399 .
the two publications which have presented dynamical measurements of the black hole mass @xcite provide values which differ by almost a factor of 3 , and neither measurement included the dark matter halo in the modelling , which could likely result in underestimating @xmath0  and thereby overestimating the ratio between the mass deficit and @xmath0 .
the next - largest mass deficit is 18 times @xmath0 , in ngc4261 , with the one other galaxy having a mass deficit larger than 10@xmath0  being ngc4486 .
the broad range of @xmath106 ( between 0.2 and 18 if we exclude ngc1399 ) is in better agreement with @xcite than with .
the former also shows a number of galaxies having mass deficits of 14 @xmath0  while the latter has no galaxies in their sample with mass deficits below 5 @xmath0 . considering that and @xcite use mostly outdated @xmath0values and an old @xmath0@xmath3  relationship for their analysis ,
it is likely that their ratios of @xmath106 are simply too high .
( for the seven of our galaxies in common with those in @xcite , we find that their mass deficits are somewhat higher than ours
median ratio @xmath107  while their black hole masses are only @xmath108% of ours . )
the simulations of @xcite indicate that mass deficits up to @xmath109 can be generated from single dry mergers .
therefore , our results can probably be explained with a reasonably low number of major dry mergers .
low @xmath106 ratios could also be obtained if gas to refill the core via new star formation is available ( i.e. , the merger is `` wet '' ) .
why does @xmath28 correlate so well with @xmath0 , while the light and the mass deficit do not ? we note that our estimate of @xmath110 , from which the mass deficit is in turn derived , is based on the assumption that the inward extrapolation of the outer srsic fit is the correct representation of the inner projected density profile of the galaxy before the binary black hole scouring takes place .
this might not be the case for real galaxies , and the large scatter about the relations of fig .
[ deficitplots ] could simply reflect this ( see , e.g. * ? ? ?
the size of the core may be a more robust quantity in this respect .
it will be interesting to investigate whether the internal kinematics of the galaxies change at this radius .
we have presented measurements for cores in 23 elliptical galaxies , where the cores are clearly detected based on the criteria suggested in @xcite .
these detections and measurements were made using core - srsic fits to extended surface - brightness profiles , the majority of which go out significantly beyond the half - light radius ; for a quarter of the sample , we found that the best fits were with the sum of core - srsic  and an outer envelope , the latter represented by an exponential or srsic  component . for 20 of these galaxies ,
reliable @xmath0 measurements are available from the literature or from our own work ; the majority of these are derived from recent stellar dynamical modelling .
we computed the light deficit in the core ( assumed to be the result of supermassive black hole mergers which follow galaxy mergers ) assuming that the original profile , before the core was formed , followed a srsic model matching the outer part of the galaxy , based on the core - srsic fit ( or the outer part of the core - srsic  component in the case of multi - component fits ) .
we then derived the corresponding stellar mass deficit using stellar mass - to - light ratios from dynamical models , available for all of the galaxies in the sample .
we find that the size of the core ( break radius @xmath89 ) is most strongly correlated with the black hole mass , while the second strongest correlation ( and the linear fit with the lowest scatter ) is with the luminosity of the host galaxy ; the weakest correlation is that between the break radius and the galaxy s stellar velocity dispersion .
the low scatter of the @xmath28@xmath0 relation ( 0.30 dex when @xmath89 is the independent variable ) opens up the possibility of accurately estimating black hole masses in core galaxies without the need for spectroscopy .
the stellar mass deficits we derive range between 0.2 and 39 times the black hole mass ; for about 75% of the sample , the range is 15 @xmath0 , and the median value for the whole sample is 2.2 @xmath0 . given the results of recent simulations ,
these values can be explained by a reasonably small number of ( dry ) mergers .
we thank the anonymous referee for a number useful comments that helped us significantly improve the paper .
spr was supported by the dfg cluster of excellence `` origin and structure of the universe '' .
pe was supported by the deutsche forschungsgemeinschaft through the priority programme 1177 `` galaxy evolution '' .
some of the data presented in this paper were obtained using sinfoni at the very large telescope vlt and from the mikulski archive for space telescopes ( mast ) .
the vlt is operated by the european southern observatory on cerro paranal in the atacama desert of northern chile .
stsci is operated by the association of universities for research in astronomy , inc .
, under nasa contract nas5 - 26555 .
this research is based in part on data from the isaac newton group archive , which is maintained as part of the casu astronomical data centre at the institute of astronomy , cambridge .
funding for the creation and distribution of the sdss archive has been provided by the alfred p. sloan foundation , the participating institutions , the national aeronautics and space administration , the national science foundation , the u.s .
department of energy , the japanese monbukagakusho , and the max planck society .
the sdss web site is http://www.sdss.org/.    the sdss is managed by the astrophysical research consortium ( arc ) for the participating institutions .
the participating institutions are the university of chicago , fermilab , the institute for advanced study , the japan participation group , the johns hopkins university , the korean scientist group , los alamos national laboratory , the max - planck - institute for astronomy ( mpia ) , the max - planck - institute for astrophysics ( mpa ) , new mexico state university , university of pittsburgh , university of portsmouth , princeton university , the united states naval observatory , and the university of washington .
finally , this work 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 .
we constructed the circularized @xmath96-band photometry of ic1459 by matching _
hst _ and ground based data .
inside 12 we used an _ hst _
wfpc2 f814w profile ( proposal i d : 5454 , pi : marijn franx ) derived using the software of @xcite while masking the central dusty features . at larger radii
we used ccd data from @xcite ( see also @xcite ) , calibrated with photoelectric photometry available in hyperleda .
the circularized @xmath111-band photometry of ngc1399 was constructed by matching the extended profile considered in @xcite to the profile derived using the iraf software from an _ hst _ acs f606w image ( proposal i d : 10129 , pi : thomas puzia ) of the galaxy in the inner 12 .
the light profile for this galaxy is based on that used in for the bh mass measurement . to facilitate the exploration of possible outer - envelope models ( section  [ multi - component - fits ] )
, we replaced the inner 7.5 of the previous profile with data from a deconvolved version of the _
hst_-acs f435w image ( 20 iterations of lucy - richardson deconvolution ) .
this allowed us to fit the profile without psf convolution .
the light profile was constructed using wfpc2 f814w images from the _ hst _ archive ( proposal i d : 5512 , pi : sandra faber ) and sdss @xmath112-band images .
the pc chip image of _ hst _ was deconvolved using 20 rounds of lucy - richardson deconvolution and then the profile was measured through ellipse fitting .
that profile was used for @xmath113 . at larger radii
, we used an ellipse - fit profile measured using a mosaic of sdss @xmath112-band images , constructed using swarp @xcite .
the sdss @xmath112-band magnitudes were converted to cousins @xmath114 using an @xmath115 color of 0.44 ( measured from ellipse fits to single sdss @xmath53 and @xmath112 images ) and the transformation of @xcite .
the surface - brightness profile for ngc3608 is a combination of ellipse fits to an f555w wfpc2 image ( proposal i d : 5454 , pi : marijn franx ) for @xmath116 and ellipse fits to an sdss @xmath117-band image for @xmath118 .
the sdss image was calibrated to @xmath96-band using the sdss @xmath117 and @xmath53-band images and the color transformations of lupton .
we matched the two datasets using fixed - ellipse - fit profiles to the sdss image and the wfpc2 mosaic , using the overlap between 685  to determine the best scaling and background subtraction to apply to the wfpc2 data .
the surface brightness profile for ngc3842 is calibrated to the @xmath96-band and is assembled from three pieces .
the innermost part ( @xmath119 ) comes from the published f555 wfpc2 pc1 profile of .
the profile in the intermediate radii is derived from an ellipse fit to a f555w wfpc2 mosaic image ( galaxy center on wf3 chip ; proposal i d : 265 , pi : douglas geiser ) .
this is matched to an ellipse fit to an sdss g - band image , including determining the level of sky background of the wfpc2 image .
the combination then consists of wfpc2 data at @xmath120 and sdss data at @xmath121 .
the final surface brightness profile is constructed by matching the combined profile and the profile .
for ngc4291 , the profile is a combination of ellipse fits to an wfpc2 f702w image [ proposal i d : 6357 , pi : walter jaffe ; see @xcite ] and ellipse fits to a sloan @xmath53-band isaac newton telescope wide field camera ( int - wfc ) image retrieved from the isaac newton group ( ing ) archive .
the latter was a 120s exposure originally obtained on 2004 january 24 .
the int - wfc image was calibrated to @xmath96-band using aperture photometry from @xcite .
we matched the two datasets using fixed - ellipse - fit profiles to the int - wfc image and the wfpc2 mosaic , using the overlap between 1080 to determine the best scaling and background subtraction to apply to the wfpc2 data .
the final profile used data from the wfpc2 pc chip for @xmath70 and data from the int - wfc image for @xmath72 .
the profile for this galaxy was constructed by combining data from wfpc2 f814w images from the _ hst _ archive ( proposal i d : 6579 , pi : john tonry ) and an @xmath122-band image from the wide field imager of the eso - mpi 2.2 m telescope .
the sky background for the wfpc2 data was determined by matching ellipse - fit profiles of the full wfpc2 mosaic with a profile from the wfi image .
the pc chip image of _ hst _ was deconvolved using 20 rounds of lucy - richardson deconvolution and then the profile was measured through ellipse fitting ; the resulting profile was used for @xmath123 . at larger radii
, we used an ellipse - fit profile from the wfi image .
the profile for this galaxy is a combination of ellipse fits to an f814w wfpc2 image ( proposal i d : 5454 , pi : marijn franx ) and an sdss @xmath112-band image . since patchy dust was present in the nucleus , we combined the f814w image with the corresponding f555w image to perform a simple dust extinction correction on the f814w image .
sky backgrounds for the two filters were determined by matching ellipse - fit profiles of the full wfcp2 mosaic image with ellipse - fit profiles from the sky - subtracted sdss @xmath117 and @xmath112 images .
the f555w and f814w images were also used , along with the prescriptions in @xcite , to generate @xmath96-band and @xmath114-band calibrations for each , iterating from an intial guess for the central ( @xmath113 ) @xmath124 color until the measured color converged .
we then formed a @xmath124 pc - chip image and used this to color - correct the f814w pc image , assuming a simple screen model for the dust .
the final profile used the ( extinction - corrected ) wfpc2 pc - chip data for @xmath125 and sdss @xmath112-band data at larger radii .
since ngc5846 had available the same kind of data as ngc5813 _ and _ similar patchy dust extinction in the circumnuclear region , we followed an essentially identical process for this galaxy .
the final profile used the ( extinction - corrected ) wfpc2 pc - chip data for @xmath126 ( proposal i d : 5920 , pi : jean brodie ) and sdss @xmath112-band data at larger radii ( the difference in transition radius compared with ngc5813 is due to the larger extent of circumnuclear dust in ngc5846 ) .
we used a wfpc2 f814w image ( proposal i d : 7281 , pi : roberto fanti ) and an sdss @xmath112-band image to construct the surface brightness profile based on combined ellipse fits .
the part of the profile inside 11 comes from the pc chip of wfpc2 ; and the innermost 1.5 is from the pc chip of wfpc2 after performing 20 rounds of lucy - richardson deconvolution .
the light profile is calibrated to the @xmath122-band .
this was done by first converting the sdss @xmath112-band profile to sdss @xmath53 by matching with a profile from the sdss @xmath53-band image .
the sdss @xmath53-band magnitudes were then converted to cousins @xmath122 using a @xmath127 color of 0.95 ( measured from the sdss g and r images ) and the conversion of @xcite .
the light profile of ngc7768 is mostly based on wfpc2 f814w images ( proposal i d : 8184 , pi : p. ct ) , which were combined to make both a pc image and a mosaic .
the pc image was corrected for dust extinction using the f814w image and f450w wfpc2 images , but there was still considerable residual extinction inside 1 , so an _ hst_-nicmos2 f160w image was used for the inner region .
the final profile is obtained through ellipse fits to the following : ( 1)nicmos f160w image for @xmath128 , ( 2)wfpc2 f814w pc image for @xmath40 between 0.93 and 10 and ( 3)wfpc2 f814w mosaic for @xmath129 .
the f814w image was corrected for sky background by matching a fixed - ellipse - fit profile with a similar fixed - ellipse - fit profile from an @xmath53-band image with the isaac newton telescope wide field camera ( 100s , 3 dec 2002 ) , taken from the isaac newton group archive .
the profiles were calibrated to @xmath96-band using aperture photometry from @xcite ( specifically , the int - wfc @xmath53-band image was calibrated to @xmath96 , and then the wfpc2 f814w profile was matched to the int - wfc @xmath53-band profile . )
of all the galaxies presented in fig .  [ rbplots ] , there are five galaxies which are also part of sample , i.e. ngc4472 , ngc4486 ( m87 ) and ngc4649 , ngc4261 and ngc4374 . using the photometric data available in , we derive the missing light using the core - srsic  function as described in sections [ themethod ] and [ coredeficit ] . from these five measurements , we find that some of our @xmath110 fall outside the uncertainties provided by .
the difference is largest for ngc4486 and ngc4261 ( by @xmath1300.7 magnitude or a factor of two in luminosity ) and it is not systematic .
the disagreement probably reflects the different details that were involved in the fitting and the derivation of the missing light :    \1 .
we use the circularized radius ; uses semi - major axis radius .
galaxies are not perfectly round or circular in projection .
non - zero ellipticity anywhere in the galaxy , not just in the center , will induce a deviation in @xmath110 since the light deficit in both and this work depends also on the srsic  fit to the outer part of the galaxy .
the magnitude of the deviation is difficult to estimate , since this depends very much on the characteristics of each galaxy .
a different fitting range introduces more discrepancies .
because we use circularized radius , we can use only the datapoints for which the ellipticity is given . in the photometric profiles of ,
the ellipticity is often not specified for a few innermost or outermost radii ; therefore we can not use those points in our fits .
in addition , in the particular case of ngc4486 , used a srsic fit ( from , using semi - major axis ) which excluded data outside 400 , because it apparently gave rise to a srsic index which was too large ( almost 12 ) due to the cd halo .
however , we found no difficulties in fitting a core - srsic function to the entire profile ( out to @xmath131 ) , using the circularized radius ; we obtain a srsic index of 8.9 ( almost identical to s value ) , albeit with different values of @xmath18 and @xmath132 .
we fit the galaxy profile with a core - srsic  profile , fit the outer part of the galaxy with a srsic  profile .
the critical point is the transition region between the srsic  and the core . in the way
defines the core , there is no transition region : the core immediately starts when the srsic  model no longer fits the datapoints in the inner part .
we use a core - srsic  parametrization and since we fit the @xmath29 parameter , instead of fixing it to infinity , we allow some part of the galaxy to be partly srsic , partly power - law .
this leads to a different inner limit of the radial range that is fitted only by the srsic  model , which in the end results in a different best - fitting srsic  profile .
the light deficit in is obtained by calculating the difference of the integrated intensity between the extrapolated srsic  function and the actual ( observed ) profile after psf deconvolution .
the ellipticity of the srsic  model is fixed to the ellipticity at the inner limit of the fitted radial range . since some of our light profiles are not psf - deconvolved , we instead compute the difference between the best - fitting core - srsic  function and the extrapolated srsic  component .
we also implicitly assume an ellipticity of zero when integrating the extrapolated srsic  function , since we have used a circularized profile in the first place .