Abstract

Studies have suggested that there is a strong correlation between the masses of nuclear star clusters (NSCs) and their host galaxies, a correlation which is said to be an extension of the well-known correlations between supermassive black holes (SMBHs) and their host galaxies. But careful analysis of disk galaxies—including 2D bulge/disk/bar decompositions—shows that while SMBHs correlate with the stellar mass of the bulge component of galaxies, the masses of NSCs correlate much better with the total galaxy stellar mass. In addition, the mass ratio for NSCs in spirals (at least those with Hubble types Sc and later) is typically an order of magnitude smaller than the mass ratio of SMBHs. The absence of a universal “central massive object” correlation argues against common formation and growth mechanisms for both SMBHs and NSCs. We also discuss evidence for a break in the NSC-host galaxy correlation, galaxies with Hubble types earlier than Sbc appear to host systematically more massive NSCs than do types Sc and later.

1. Introduction

As far as we can tell, all massive galaxies in the local universe harbor supermassive black holes (SMBHs, with masses ). The masses of these SMBHs correlate strongly with several global properties of the host galaxies, particularly with the central velocity dispersion [1, 2] and with the bulge luminosity or mass (e.g., [3, 4]). These correlations imply that the processes which drive galaxy growth and the processes which drive black hole growth are intimately linked—perhaps even the same processes.

It is now also clear that many galaxies, particularly later-type spirals, host luminous nuclear star clusters (NSCs; e.g., [5, 6]), with masses in the range ; see the review by Böker [7] for more details. Recently, several authors have argued that NSCs and central SMBHs have the same host-galaxy correlations, in particular, that SMBHs and NSCs have the same correlation with bulge luminosity and mass [811] (but see Balcells et al. [12]). The suggestion, then, is that NSCs and SMBHs are in a sense members of the same family of “Central Massive Objects” (CMOs), and thus that they may have grown via the same mechanisms (e.g., [1316]).

We argue, however, that one should be cautious about assuming that NSCs and SMBHs are really part of the same family, with the same host-galaxy relationships. To begin with, the samples of Wehner and Harris [8] and Ferrarese et al. [9], which were used to make the CMO argument, were almost entirely early-type galaxies—mostly ellipticals and dwarf ellipticals. These are galaxies which are, in essence, “pure bulge” systems, so one could just as easily argue for a correlation with total galaxy mass. But we know that SMBHs in disk galaxies correlate better with just the bulge, and not with the total galaxy mass or light (e.g., [17, 18]). Given that there have been previous claims that NSCs in spiral galaxies correlate with the total galaxy light (e.g., [19]), we are prompted ask the question: do nuclear clusters in disk galaxies correlate with the bulge (like SMBHs), or with the whole galaxy?

2. Samples, Methodology, and Data Sources

Although current studies suggest that the relation is tighter and has less intrinsic scatter than the relation (e.g., [20]), velocity dispersion is not the ideal host-galaxy measure to use here, for three reasons. First, most of the best-determined NSC masses are based directly on the measured velocity dispersion of the NSC (e.g., [21]), which is often indistinguishable from that of the surrounding bulge; this means a (spurious) correlation between NSC mass, and central velocity dispersion is only to be expected. Second, some NSCs are found in galaxies with no detectable bulge at all (see discussion in Section 3). Finally, it is difficult to see how one should discriminate between a velocity dispersion due to the bulge versus one due to the whole galaxy. But discriminating between bulge and whole-galaxy luminosities and masses is much simpler. So we choose instead to compare NSCs and their host galaxies with the relation, which means comparing NSC masses with the stellar masses of host galaxies and their bulges.

For NSCs, we emphasize galaxies where the NSC masses have been dynamically measured, since this is the most direct analog to well-determined SMBH masses (i.e., those with direct dynamical mass measurements from stellar, gas, or maser kinematics, where the SMBH sphere of influence is resolved). In addition, dynamical measurements avoid possible problems with multiple stellar populations; the latter can potentially bias stellar masses estimated from broadband colors. Spectroscopic studies [11, 22] have shown that NSCs often contain multiple stellar populations; this renders mass estimates based on single stellar population (SSP) models (e.g., those used by [9]) somewhat uncertain. The NSCs we focus on are taken primarily from the sample of Walcher et al. [21], with additional data from Ho and Filippenko [23], Böker et al. [24], Kormendy and Bender [25], Matthews et al. [26] and Gebhardt et al. [27], Barth et al. [28], Seth et al. [29], and Kormendy et al. [30]; we use the estimate of Launhardt et al. [31] for the Milky Way’s NSC. This gives us a total of 18 galaxies with dynamically determined NSC masses. These cover Hubble types S0–Sm, but the sample is in fact heavily biased towards later types; over three-quarters are Hubble types Second or later. As an additional, secondary sample, we include 15 galaxies from Rossa et al. [11], where the masses are estimated by fits of multiple SSP models to high-resolution spectroscopy. Most of these galaxies are Sc and later, but a few earlier-type spirals (Sa–Sb) are also included.

Total stellar masses are based on -band total magnitudes from 2MASS [32] or from Malhotra et al. [33] for M31 and M33 (which are too large for accurate sky subtraction of 2MASS images), combined with color-based mass-to-light () ratios from Bell et al. [34]. For the latter, we use optical colors from the literature (primarily from HyperLeda (http://leda.univ-lyon1.fr/)) or from direct measurements on Sloan Digital Sky Survey (SDSS, [35]) images. The bulge masses are derived using bulge-to-total () values determined individually for each galaxy by 2D image decomposition, using the BUDDA software package [36, 37], which incorporates bulge and disk components and optional bars and central point sources (the latter can be used for both nuclear star clusters and AGN). Note that we explicitly define “bulge” to be the “photometric bulge”, that is, the excess light (and stellar mass) which is not part of the disk, bar, or nuclear star cluster. We defer questions of how SMBH (or nuclear cluster) mass relates to so-called “pseudobulges” versus “classical bulges” (e.g., [38, 39]) to a later analysis.

Full 2D decompositions, as described above, were used for all S0 and spiral SMBH host galaxies. For the NSC host galaxies, we follow the same approach, with one simplification. Since we have found that ratios for unbarred galaxies do not change dramatically if we use 1D surface-brightness profile decompositions instead of 2D image decompositions, we use the former for genuinely unbarred galaxies; we are careful to exclude (or separately model) the NSC contribution to the surface-brightness profile in these cases. Galaxies which do possess bars are subjected to full 2D decompositions; see the following section for details.

2.1. Bulge-Disk Decompositions

As noted previously, we use 2D image decompositions via the BUDDA software package to determine the ratios, and thus the bulge stellar masses, for SMBH host galaxies and for barred NSC host galaxies. For the NSC galaxies, we use HST data wherever possible, to enable the NSC itself to be properly modeled as a separate source. However, we have found that when the NSC is sufficiently luminous, and when the bulge is sufficiently low-contrast, we can achieve reasonable decompositions with ground-based images; these are sometimes preferable if they are near-IR (to minimize the effects of dust extinction and recent star formation) and/or large enough to include the entire galaxy (to allow better recovery of the disk component).

We have completed decompositions for the galaxies with dynamically determined NSC masses (we use the published 2D decomposition of [28] for NGC 3621); in the special case of the Milky Way, we assume a bulge mass of and a total stellar mass of , based on arguments in Dehnen and Binney [40], Klypin et al. [41], and Flynn et al. [42]. Decompositions for the spectroscopic sample are still in progress but are mostly complete; since our primary analysis (e.g., computing the relation) is based on the dynamical masses, the incompleteness of the spectroscopic sample does not affect our results.

Full details of the individual decompositions will be published elsewhere (Erwin and Gadotti 2012a, in prep.). An example of one of the 2D decompositions is given in Figure 1 for the galaxy NGC 7418, where we fit an -band image from the Ohio State University Bright Spiral Galaxy Survey [43, OSU BSGS] using an exponential disk (95.6% of the light), a Sérsic bulge (1.6% of the light), a bar (2.4% of the light), and a point source for the NSC (0.5% of the light). (In this particular galaxy, the disk appears to be truncated, but this has little effect on the decomposition; including a broken-exponential profile for the disk changes the ratio from 0.016 to 0.017.) This illustrates the importance of including a separate bar component in the decomposition when the galaxy is barred; the bar has almost twice the luminosity of the bulge, and, without it, the bulge luminosity (and stellar mass) would certainly be overestimated. In fact, a 1D decomposition for this galaxy gives a value almost twice as large (0.030); similar results were found for four other barred galaxies in the sample, with mean values a factor of 2.1 times larger when the bar was omitted; see also [37].

3. Comparing Black Holes and Nuclear Star Clusters

Although black-hole–bulge correlations are sometimes described as correlations between the black hole mass and the host galaxy mass (or luminosity as a proxy for mass), this is really only true for elliptical galaxies, where the entire galaxy is the “bulge.” Kormendy and Gebhardt [17] explicitly compared -band total and bulge luminosities for SMBH hosts and showed that the latter provided a much better correlation. Most recently, Kormendy et al. [18] have show for a larger, updated sample that SMBH masses in disk galaxies correlated much better with (classical) bulge -band luminosity than with the luminosity of the disk component; this naturally suggests that total-galaxy luminosity is unlikely to correlate well with SMBH mass when the galaxy is disk dominated.

In Figure 2, we compare SMBH masses with total galaxy stellar mass (Figure 2(a)) and with bulge stellar mass (Figure 2(b)), based on our careful bulge/disk/bar decompositions (see Table 1). Error bars include the effects of uncertainties in the distance and in the and ratios. As expected, the correlation between SMBH mass and bulge mass is much stronger than any correlation with total galaxy mass; the Spearman correlation coefficients are for the relation versus 0.29 for the relation, with the latter correlation not being statistically significant.

We also plot linear fits of as a function of ; these fits are made using galaxies with well-determined distances (filled points) to minimize distance-based uncertainties, using the Bayesian-based approach of D’Agostini [44], which explicitly incorporates errors in both variables and intrinsic scatter in the black hole mass (see also [45, 46]). By “well-determined distances,” we mean those determined using direct methods such as surface-brightness fluctuations and Cepheid stars, or redshift-based distances, where (to avoid large relative uncertainties due to peculiar motions.) The best-fitting relation for the whole sample (black line; the fit to just the elliptical galaxies, shown by the red line, is almost identical) is with intrinsic scatter in SMBH mass of dex; the best-fitting relation for the bulges of disk galaxies only (blue line) is with intrinsic scatter . (The errors are based on bootstrap resampling.)

We apply exactly the same methodology to NSC-host galaxies in Figure 3, plotting NSC mass versus total galaxy stellar mass in the left panel and versus bulge stellar mass in the right panel. Since several of the NSC host galaxies are genuinely bulgeless systems (without even a distinct “pseudobulge”), we plot their bulge masses as upper limits (). As the figure shows, NSC mass clearly correlates better with total stellar mass than it does with bulge mass. (The respective correlation coefficients are versus 0.38; the bulge-mass correlation is not statistically significant.) Fitting NSC mass versus total stellar mass, using the same methodology as for the SMBH fits, gives the following relation: with intrinsic scatter dex. Note that the slope is formally indistinguishable from unity; that is, the ratio does not appear to depend on itself.

It is important to note that the difference in correlation coefficients actually underestimates the true difference between the two relations, because the correlation was computed assuming that bulgeless spirals still have nominal bulges (using ). In the combined sample of dynamical and spectroscopic NSC masses, we can identify at least three galaxies which have no detectable bulge. In two of these (NGC 1493 and NGC 2139), our 2D decomposition assigned stellar light to a bar in addition to a pure exponential disk; in 1D decompositions (or simple bulge + disk 2D decompositions), light from the bar might be (wrongly, we would argue) interpreted as “bulge” light. For the other galaxy (NGC 300), however, there is no ambiguity; this is an unbarred spiral galaxy with a surface brightness profile consisting of only an exponential disk and the NSC (see, e.g., Figure 8 of [47]).

The existence of nuclear star clusters in genuinely bulgeless spirals is simply an additional, direct confirmation of our basic conclusion: nuclear star cluster masses scale with the total stellar mass of their host galaxies, not with the bulge mass. This means that NSCs and SMBHs do not follow a common host-galaxy correlation.

We have also investigated whether other galaxy parameters might correlate with NSC mass, or even with residuals from the relation. In particular, we have compared NSC mass with rotation velocity and with total baryonic mass (stellar mass plus atomic gas from HI measurements). In both cases, correlations exist, but they are not as strong as the correlation with total stellar mass. No particular correlations with residuals of the relation are seen.

Closer inspection of the left-hand panel of Figure 3 suggests that the spectroscopic masses (black diamonds) tend to be offset from the NSC- relation, in the sense that they have larger NSC masses for the same total stellar mass. This could, in principle, be evidence of a systematic overestimation of NSC masses in the spectroscopic sample, but of the four galaxies in common between Walcher et al. [21] and Rossa et al. [11] only one has a (slightly) higher spectroscopic mass, while the other three have spectroscopic masses slightly lower than the dynamical masses. There is, however, another difference to consider; the spectroscopic sample tends to have earlier Hubble types.

This brings us to something which Seth et al. [48] pointed out several years ago, using a larger dataset of NSCs and host galaxies, with NSC masses based (mostly) on colors or assumed ratios. They noted that NSCs in late-type spirals tended to have lower relative masses () than early-type spirals and ellipticals. (Rossa et al. [11] pointed out a similar trend in absolute NSC mass for their smaller sample of NSCs in early- and late-type spirals.) Figure 4 makes this explicit by plotting versus Hubble type for the galaxies in Seth et al.’s compilation, plus seven galaxies from our updated dynamical-mass sample which were not in their sample. We have also added galaxy stellar-mass estimates for 16 galaxies that did not have masses in Seth et al., using total -band magnitudes from 2MASS and either colors from HyperLeda or measured colors from SDSS images to derive the -band via Bell et al. [34].

What is curious about Figure 4 is not just that the ratio depends on Hubble type, but that it actually appears to do so in a bimodal fashion; Hubble types Sb and earlier have relatively large NSC masses, while Sc and later-type galaxies have significantly smaller relative NSC masses. Plotted on top of the figure are simple fits of a function where the ratio can take two constant values, one for Hubble types and the other for , with a simple linear transition between and . Fits to just the dynamical + spectroscopic masses (red dashed line) and to the entire sample (gray dashed line) are similar, indicating that Sb and earlier Hubble types form one class, with , and Sc and later types form a different group, with almost an order of magnitude smaller (~0.0003). The corresponding best-fit values of (, ) are (3.51, 4.05) for the dynamical + spectroscopic masses and (3.10, 5.01) for the complete sample. As a crude check on whether this split is statistically significant, we performed Kolmogorov-Smirnov tests on the values of for galaxies with and galaxies with . The K-S test gives a probability for the two sets of ratios coming from the same parent population if we use only the dynamical + spectroscopic masses, or if we use the entire set of NSC masses.

Do NSCs in late-type spirals differ from those in early-type spirals, S0s, and ellipticals in any sense other than average mass? The available evidence is ambiguous. Böker [7] notes that NSC sizes appear to be independent of Hubble type. On the other hand, Rossa et al. [11] compared stellar populations of NSCs in early- and late-type spirals using fits to their spectroscopy and noted that the NSCs in late-type spirals did tend to have younger stellar populations and (slightly) lower metallicities. (They also argued against any observational effects that might produce systematic overestimates of NSC mass in early spirals.) This does at least suggest that different star-formation histories may lie behind the mass differences in NSCs.

We also plot the ratio for SMBH host galaxies (thick gray dotted line in Figure 4, based on (1)). What this indicates is that the NSC—host-galaxy relationship for Sb and earlier types is consistent with the SMBH relation, if all of the galaxy mass is in the bulge. Since most of the galaxies used for the original CMO studies [8, 9] were dwarf and giant ellipticals (or S0 galaxies with high ratios), it is easy to see why the “NSC = SMBH’’ connection could be made. But this is clearly true only for very bulge-dominated systems.

Appendix

Data Tables

In Tables 1 and 2, we list the basic data parameters for SMBH and NSC hosts. References for the NSC masses are in the captions for Table 2. For the SMBH masses, the numbers in column 5 of Table 1 translate into the following references: 1 = Gillessen et al. [49]; 2 = Bender et al. [50]; 3 = Verolme et al. [51]; 4 = Krajnović et al. [52]; 5 = Gebhardt et al. [53]; 6 = Bower et al. [54]; 7 = Lodato and Bertin [55]; 8 = Atkinson et al. [56]; 9 = Rusli et al. [57]; 10 = Nowak et al. [58]; 11 = Houghton et al. [59]; 12 = Sarzi et al. [60]; 13 = Devereux et al. [61]; 14 = Davies et al. [62]; 15 = Barth et al. [63]; 16 = Nowak et al. [39]; 17 = Kondratko et al. [64]; 18 = Güeltekin et al. [65]; 19 = de Francesco et al. [66]; 20 = Hicks and Malkan [67]; 21 = Miyoshi et al. [68]; 22 = Ferrarese et al. [69]; 23 = Cretton and van den Bosch [70]; 24 = Walsh et al. [71]; 25 =Macchetto et al. [72]; 26 = Nowak et al. [73]; 27 = Shen and Gebhardt [74]; 28 = de Francesco et al. [75]; 29 = Neumayer et al. [76]; 30 = Capetti et al. [77]; 31 = Ferrarese and Ford [78]; 32 = van der Marel and van der Bosch [79]; 33 = Cappellari et al. [80]; 34 = Bontà et al [81].

Acknowledgments

The authors would like to thank Eva Noyola and Anil Seth for useful and interesting conversations, along with comments from two referees which improved the manuscript. They would also like to thank the organizers of the 2008 “Nuclear Star Clusters across the Hubble Sequence” workshop in Heidelberg and the 2010 “ESO Workshop on Central Massive Objects” in Garching for helping motivate and inspire this research. This work was supported by Priority Programme 1177 (“Witnesses of Cosmic History: Formation and evolution of black holes, galaxies and their environment”) of the Deutsche Forschungsgemeinschaft. 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 Ohio State University. It is also based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archive at the Space Telescope Science Institute. STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 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 research made use of the Lyon-Meudon Extragalactic Database (LEDA; part of HyperLeda at http://leda.univ-lyon1.fr/), and 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.