{'fraction_non_alphanumeric': 0.04363708164675684, 'fraction_numerical': 0.03406061802744595, 'mean_word_length': 4.240300051733057, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Effect of Spatial Gradients in Stellar Mass-to-Light Ratio on Black Hole Mass Measurements

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease).


INTRODUCTION
Dynamical measurements of black hole masses (M • ) at the centers of nearby galaxies have opened new avenues for studying galaxy evolution, by exposing correlations between M • and various host galaxy properties. Ongoing studies of these black hole scaling relations explore whether the growth and quenching of black holes and galaxies are causally linked, or only loosely connected by broad trends in cosmic evolution. One way to improve on the determination of the black hole scaling relations and the cosmic scatter in M • is to reduce uncertainties and systematic biases in measurements of M • in local galaxies.
Among the approximately 70 nearby galaxies with dynamically determined M • , about 50 M • have been measured using stellar kinematic data and orbit superposition models (McConnell & Ma 2013 and references therein). A number of assumptions made in the models can contribute to potential biases in the inferred black hole masses. Recent advances include modeling stellar mass and dark matter as separate components (e.g., Gebhardt & Thomas 2009), and accounting for triaxiality (e.g., van den Bosch & de Zeeuw 2010).
Another assumption that has been used in most previous orbit models is that the stellar mass-to-light ratio, Υ, is constant throughout a galaxy. This assumption is clearly a simplification since spatial gradients have been observed for many galaxy properties: e.g., color, metallicity, age, and α-element enhancement, which can all be associated with spatial variations in Υ. Early-type galaxies are often bluer towards their outer regions. Decreasing metallicity outward is likely the dominant cause for the color gradients in these galaxies (e.g., Strom et al. 1976;Tamura et al. 2000), with the age gradient playing a minor or insignificant role (e.g., Saglia et al. 2000;Tamura & Ohta 2004). This picture is supported by recent photometric studies of large statistical samples of early-type galaxies in SDSS (Tortora et al. 2010a) as well as spectroscopic studies of individual galaxies (Spolaor et al. 2009;Rawle et al. 2010;Kuntschner et al. 2010). The spatial coverage of these studies is typically limited to the central regions of the galaxies within one effective radius ( r eff ). Measurements of ten earlytype galaxies have been extended to larger radii using deep observations with integral field spectrographs (Weijmans et al. 2009;Greene et al. 2012). The central gradients in metallicity and absorption line strengths are seen to continue to 2-4 r eff . Gradients in color and inferred metallicity in galaxies' globular cluster systems persist to ∼ 8 r eff (Arnold et al. 2011;Forbes et al. 2011).
Tight positive correlations between Υ and color have been found in both late-type and early-type galaxies (e.g., Bell & de Jong 2001;Bell et al. 2003;Tortora et al. 2011), with redder colors corresponding to larger Υ. By fitting synthetic spectral models to the observed optical color gradients, Tortora et al. (2011) have found mild negative gradients in Υ for local early-type galaxies with stellar ages > 6 Gyr, over a large range of stellar masses (M ⋆ ∼ 10 9.0 -10 11.3 M ⊙ ). Color and mass profiles have also been recently investigated in galaxies at 0.5 < z < 2.5 (Szomoru et al. 2013).
In this letter, we assess the amount of systematic error in M • resulting from the standard assumption of constant Υ in the orbit models. We modify the orbit superposition code by Gebhardt et al. (2000) and introduce a spatial gradient in Υ to the stellar component of the potential. We investigate the effects on M • for three giant elliptical galaxies that we have previously analyzed (assuming constant Υ): NGC 3842, NGC 6086, andNGC 7768 (McConnell et al. 2011a,b, 2012). All three are massive ellipticals and brightest cluster galaxies (BCGs) with high stellar mass M ⋆ and velocity dispersion σ: M ⋆ = 1.55 × 10 12 , 1.43 × 10 12 , 1.16 × 10 12 M ⊙ and σ = 270, 318, 257 km s −1 for NGC 3842, NGC 6086, NGC 7768, respectively. The black hole masses were predominantly constrained by stellar kinematics from the integral-field spectrographs OSIRIS (Larkin et al. 2006) and GMOS (Allington-Smith et al. 2002;Hook et al. 2004) Stellar kinematics at large radii were recorded with the Mitchell Spectrograph (Hill et al. 2008) or adopted from the literature (Carter et al. 1999;Loubser et al. 2008).
A few prior studies of M • have explored variations in Υ. Recent investigations of S0 galaxies NGC 1332, NGC 3368, and NGC 3489 modeled bulge and disk components with separate Υ values (Nowak et al. 2010;Rusli et al. 2011). Nowak et al. (2007 assumed a separate Υ value for the nuclear disk in NGC 4486A. Gebhardt et al. (2000) assessed the V − I color gradients in NGC 3379 and applied the corresponding Υ gradient to one set of models; Cappellari et al. (2002) applied a similar treatment to IC 1459. For both galaxies, models with a Υ gradient and models with uniform Υ yielded statistically consistent measurements of M • . Here we extend from prior investigations by examining systematic trends in measured M • over a range of Υ gradients.

STELLAR ORBIT MODELS AND Υ GRADIENTS
We use the axisymmetric orbit modeling algorithm of Gebhardt et al. (2000Gebhardt et al. ( , 2003, Thomas et al. (2005), and Siopis et al. (2009). In this model, a galaxy is described by the density profile where ν(r, θ) is the luminosity density derived from the galaxy's deprojected surface brightness profile, and ρ dm (r) represents a spherical NFW or cored logarithmic (LOG) dark matter profile. Equation (1) is used to construct a gravitational potential, and kinematic models of the galaxy are constructed by propagating test particles through the potential and computing their timeaveraged velocities ("orbits") throughout a polar grid. Orbital weights are varied to fit the observed kinematics of the galaxy but are constrained such that the sum of weighted orbits exactly reproduces ν(r, θ), which is decoupled from variations in Υ(r). The model velocity distributions are then compared to the observed kinematics, yielding a goodness-of-fit statistic χ 2 . Each instance of the model adopts a single density profile and outputs a single χ 2 value. The best-fit values and confidence limits for M • and other parameters in Equation (1) are determined by analyzing the distribution of χ 2 from many models. Previous investigations using orbit superposition models typically have set Υ(r) to a uniform value, Υ 0 . To -Color vs. radius in three giant elliptical galaxies based on photometric data from the Hubble Space Telescope and Kitt Peak National Observatory. In each panel, the black dashed line represents the best log-linear fit to all data points, and the solid grey line represents the best log-linear fit to data with r ≤ r eff . For NGC 6086, our V −I data extend to only 0.5 r eff (17 ′′ ). B −R data from Peletier et al. (1990) extend to 1.2 r eff and have a logarithmic slope of −0.06. assess the impact of a gradient in Υ, we set For each value of α, we sample a two-dimensional grid of values for M • and the normalization of Υ. The galaxies' color gradients indicate likely values of α. Figure 1 shows the color profile for each galaxy in our sample. The color gradients δ(B-R)/δ(log 10 r) and δ(V -I)/δ(log 10 r) range from −0.02 to −0.06. In Figure  2 of Tortora et al. (2011), the matching gradients in Υ are ≈ −0.1 to 0 for SDSS galaxies. Their Figure 3 shows the Υ gradients as a function of stellar mass and velocity dispersion and indicates different mean values of α for galaxies with different stellar ages. Early-type galaxies with central age older than 6 Gyr have α ∼ −0.2 to −0.1 over a large range of σ, but show a slight upturn toward α = 0 at M ⋆ > 10 11 M ⊙ .
For NGC 3842, we have an additional two hours of on-source data from May 2012 observations with the Mitchell Spectrograph, allowing for a direct measurement of the Υ gradient. As described in detail in Greene et al. (2012), we construct elliptically averaged radial bins to maximize the S/N in the spectra. We have a median S/N of 40 per pixel in the largest bin, at r ≈ 1.1 r eff . First, weak emission features are iteratively corrected (Graves & Faber 2010) and velocity dispersions are measured at each bin using pPXF (Cappellari & Emsellem 2004). Then the stellar population properties are measured from the Lick EW+EZ Ages algorithm (Graves & Schiavon 2008). We use the α−enhanced models of Schiavon (2007), assuming a Salpeter IMF, to determine the luminosityweighted Υ at each radius. We find the resulting Υ gradient to be α = −0.13 ± 0.30, broadly consistent with the color gradient in NGC 3842. Uncertainties in the inner bins are dominated by systematic uncertainties in emission line and velocity dispersion measurements, while the outer bins include additional uncertainties from sky subtraction.
Gradients in the stellar initial mass function (IMF) are potentially an additional source of Υ gradients within our galaxies. The IMF varies with the [Mg/Fe] line ratio and the corresponding inferred timescale for star formation (Conroy & van Dokkum 2012). Our line index measurements from Mitchell Spectrograph data of NGC 3842 are consistent with no radial gradient in [Mg/Fe]. Therefore we do not expect IMF gradients to alter the likely range of α.
Given the uncertainties discussed above, we have chosen to test α ranging from −0.2 to 0.1. We employ α as a spherical, three-dimensional gradient in the models. Integrating our model stellar mass profiles along the line of sight, we find the projected gradients to be ≈ 80% of the three-dimensional α. The corresponding range of projected gradients, −0.16 to 0.08, is appropriate with regards to the inferred gradients from photometry and spectroscopy. Table 1 and Figures 2 and 3 illustrate our best-fit values of M • and Υ from modeling NGC 3842, NGC 6086, and NGC 7768 with different index α for the Υ gradient. The stellar orbit models optimize the normalization of Υ, while the radial dependence is set by Equation (2) for fixed α. Each galaxy's distance D, effective radius r eff , and black hole radius of influence, r inf ≡ GM • /σ 2 , are listed in Table 1. Since r inf is much smaller than r eff (with the ratio ranging from 0.006 to 0.03), we list Υ at three radii ( r inf , 0.1 r eff , and 1 r eff ) for comparsion. We also list the average luminosity-weighted stellar mass-tolight ratio Υ , equal to the total stellar mass divided by the total luminosity. Table 1 and Figure 2 show that within r inf , there is a mild degeneracy between stellar mass and black hole mass. As α increases from negative toward positive values, Υ( r inf ) decreases, causing the best-fit M • to in-  crease. At large radii the trend with α is reversed, such that the best-fit Υ( r eff ) and Υ both increase as α increases. The intermediate radius of ∼ 0.1 r eff is roughly the pivot point at which the best-fit Υ has similar values for different α. Our kinematic data for each galaxy fall mainly within 1 r eff and are most thoroughly represented at small radii (r < ∼ 3 ′′ , or ∼ 0.1 r eff ), where the OSIRIS and GMOS integral-field spectrographs yield measurements at multiple position angles. This kinematic coverage drives the normalization of Υ, such that best-fitting models for each value of α have similar enclosed masses near 0.1 r eff for a given galaxy. This is also the scale on which stars dominate the total enclosed mass: much within ∼ 0.1 r eff , the black hole begins to contribute, whereas the dark matter begins to dominate at ∼ 0.5 − 2 r eff (Figure 4). Figure 3 shows a clear positive correlation between α and the best-fit M • value. As α decreases from 0 to −0.2, the best-fit M • decreases by 28% in NGC 3842 and 27% in NGC 6086. For these galaxies, the bias in M • over this range in α is comparable to the statistical uncertainty in M • for each suite of models with fixed α. NGC 7768 exhibits a slightly smaller bias of 17% in M • , versus statistical uncertainties ∼ 40% from individual trials.

RESULTS
For each galaxy, the value of χ 2 for the best-fit model  -Contours of χ 2 as a function of Υ( r inf ) and M•, for orbit superposition models with different slopes for the logarithmic Υ gradient: α = −0.1 (blue dashed), 0 (red solid), and +0.1 (green dotted). For each α, the thick and thin contours represent 68% and 90% confidence levels (corresponding to ∆χ 2 = 1.0 and 2.71), respectively.
(for a given α) changes by less than 1.0 when α increases from −0.2 to 0. By contrast, χ 2 for the best-fit model increases by 2.5 as α increases from 0 to +0.2 for NGC 3842, and χ 2 increases by 1.0 and 2.8 for NGC 6086 and NGC 7768, respectively, as α increases from 0 to +0.1. Models with a positive gradient in Υ are therefore mildly disfavored by our kinematic data. This trend is consistent with the range of α inferred from the three The y-axis is the ratio of the best-fit M• value relative to the best-fit value for models with uniform Υ (α = 0). Error bars represent the statistical errors from individual trials (Table 1), also normalized relative to M•(α = 0). galaxies' negative color gradients shown in Figure 1. Figure 4 shows the relative contributions of stars, dark matter, and the black hole to the total enclosed mass as a function of deprojected radius in NGC 3842. Similar trends are seen for NGC 6086 and NGC 7768. For simplicity, we have modeled a single dark matter halo for each galaxy, for all values of α. As α decreases from 0 to −0.2, our kinematic data prefer a model with more (less) stellar mass at r ≪ 0.1 r inf (r ≫ 0.1 r inf ). Since we have assumed a fixed dark matter halo, the resulting enclosed stellar mass fraction is decreased at large radii.
If a range of dark matter profiles were sampled, we would expect changes in Υ(r) to alter the mass and/or shape of the best-fitting dark matter halo. Introducing a negative Υ gradient would redistribute stellar mass toward smaller radii, and the models could reconstruct the outer mass profile by fitting a more massive halo, amplifying the trend for dark matter in Figure 4. Our models are are not sensitive to the inner slope of the dark matter profile, as the enclosed mass at small radii is dominated by stars and the black hole.
Since the dark matter halo and Υ gradient both cause a non-uniform total mass-to-light ratio, one might suspect that optimizing the dark matter profile for each value of α would decrease the variations in Υ normalization and M • . While our kinematic data for NGC 3842, NGC 6086, and NGC 7768 do not extend far enough to tightly constrain the dark matter profile, we have performed a simple test by varying α in models of NGC 3842 with no dark matter. Even without dark matter, the overall trend of Υ and M • versus α is similar to those depicted in Figures 2a and 3.

DISCUSSION
We have run stellar orbit superposition models of three giant elliptical galaxies, assuming a range of gradients defined by α = d log(Υ)/d log(r) in the stellar mass-tolight ratio Υ. We find that the best-fit black hole mass M • varies systematically with the strength and sign of α. As α decreases from 0 to −0.2, the best-fit values of M • decrease by 28%, 27%, and 17% for NGC 3842, NGC 6086, and NGC 7768, respectively. As α increases from 0 to +0.1, the best-fit values of M • increase by 14%, 22%, and 17% for the three galaxies. For comparison, individual trials yield statistical errors ∼ 30%-40% in M • . Color and line strength gradients in these three galaxies suggest α ∼ −0.2 to −0.1.
Our results suggest that gradients in Υ may be a nonnegligible source of systematic error in measurements of M • and M ⋆ . The overall effect of this error on the correlations between M • and host galaxy properties will depend on the distribution of Υ gradients across the sample of galaxies with existing M • measurements.
All three galaxies in our study are BCGs with M ⋆ > 10 12 M ⊙ and velocity dispersions σ > 250 km s −1 . Extrapolating the trend illustrated in Figure 3 to other values of α will indicate the potential level of bias in dynamical measurements of M • that assume uniform Υ. A systematic variation in α with galaxy stellar mass or velocity dispersion can have an impact on the slopes of the M • -M bulge and M • -σ scaling relations, which are currently determined assuming α = 0 (e.g., McConnell & Ma 2013 and references therein). For early-type galaxies, Tortora et al. (2011) find a non-monotonic trend in the dependence of α on M ⋆ and σ, where α has a minimum around M ⋆ ∼ 1-3 × 10 10 M ⊙ (σ ∼ 100-160 km s −1 ) and increases mildly toward the low-and high-M ⋆ (σ) ends. The scatter, however, is large. If α decreases (increases) with decreasing M bulge or σ, our results here indicate that the inferred M • would be systematically lower (higher) for smaller galaxies, thereby steepening (flattening) the slopes of the scaling relations. In some galaxies where kinematic measurements superbly resolve r inf , lower degeneracies between M • and Υ will reduce systematic errors in M • .
In addition to biasing measurements of M • , gradients in Υ may impact other attempts to decompose mass profiles into multiple components. In particular, strong and weak lensing provide complementary data to stellar kinematics for probing the total mass profiles of galaxies and galaxy clusters. In contrast to our small-scale kinematic data for NGC 3842, NGC 6086, and NGC 7768, lensing studies typically constrain the total enclosed mass at larger radii of ∼ 1 − 100 r eff . In this case, the inner slope of the dark matter profile can be degenerate with the stellar mass profile. Current methods estimate the stellar mass component using stellar kinematics and/or population synthesis, under the assumption of spatially uniform Υ (e.g., Jiang & Kochanek 2007;Sand et al. 2008;Auger et al. 2010;Tortora et al. 2010b;Newman et al. 2011Newman et al. , 2012. Allowing for a gradient in Υ in these studies could affect the decomposition of the stellar and dark matter components. For instance, in the presence of a negative gradient in Υ, models sensitive to the enclosed mass slope near 1 r eff might prefer a less cuspy dark matter profile, to counteract the inward steepening in stellar mass. Gravitational lensing is better suited for probing lensing galaxies at z ∼ 0.2 and beyond, whereas direct dynamical measurements of M • are currently feasible only in local galaxies within ∼ 150 Mpc. Other dynamical tracers such as globular clusters have been used to infer galaxy mass profiles out to several effective radii (e.g. Romanowsky et al. 2009;Murphy et al. 2011). A joint analysis combining the kinematics from our small-scale studies with those from the large-scale tracers can provide more comprehensive measurements of the mass profiles of stars and dark matter in early-type galaxies, from the black hole sphere of influence out to several effective radii.