Central Black Holes in Late-type Galaxies

Chandra Survey of Nearby Galaxies: A Significant Population of Candidate Central Black Holes in Late-type Galaxies

Rui She11affiliation: Department of Engineering Physics and Center for Astrophysics, Tsinghua University, Beijing 100084, China , Luis C. Ho22affiliation: Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100087, China 33affiliation: Department of Astronomy, School of Physics, Peking University, Beijing 100087, China , and Hua Feng11affiliation: Department of Engineering Physics and Center for Astrophysics, Tsinghua University, Beijing 100084, China

Based on the Chandra data archive as of March 2016, we have identified 314 candidate active galactic nuclei in 719 galaxies located closer than 50 Mpc, among them late-type (Hubble types Sc and later) galaxies that previously had been classified from optical observations as containing star-forming (H ii) nuclei. These late-type galaxies comprise a valuable subsample to search for low-mass ( ) central black holes. For the sample as a whole, the overall dependence of the fraction of active nuclei on galaxy type and nuclear spectral classification is consistent with previous results based on optical surveys. We detect 51 X-ray cores among the 163 H ii nuclei and estimate that, very conservatively, 74% of them with luminosities above  erg s are not contaminated by X-ray binaries; the fraction increases to 92% for X-ray cores with a luminosity of  erg s or higher. This allows us to estimate a black hole occupation fraction of % in these late-type, many bulgeless, galaxies.

Subject headings:
galaxies: active — galaxies: nuclei — galaxies: Seyfert — X-rays: galaxies

1. Introduction

Nearby galaxies offer us the unique laboratories to search for and study weak active galactic nuclei (AGNs), including those emanating from low-mass central black holes ( ). Although their formation mechanism is still poorly constrained, low-mass black holes may have played a key role in forming supermassive black holes (SMBHs) in the early universe. The early manifestation of SMBHs at high redshifts (Mortlock et al., 2011; Wu et al., 2015) requires the growth of these giants to start with “seed” black holes that are significantly more massive than stellar-mass black holes of typical mass 10 (Volonteri, 2010). Some of the seed black holes may have survived until today and appear as low-mass black holes lurking at the centers of low-mass (late-type spiral and dwarf) galaxies (Volonteri et al., 2008; Greene, 2012; Reines & Comastri, 2016). Present-day low-mass galaxies, to the extent that they have remained relatively unevolved since their formation, are the optimal sites to search for such relic seed black holes. Finding these low-mass black holes, knowing their occupation fraction in low-mass galaxies, and measuring their mass function may help constrain their formation mechanism and distinguish between the two competing models: whether they are Population III remnants (Heger et al., 2003) or due to direct collapse (Haehnelt & Rees, 1993).

According to the scaling relations between black hole mass and host galaxy properties (Kormendy & Ho, 2013), the best place to unveil low-mass black holes is in late-type bulgeless galaxies. Optical and X-ray searches have been conducted in the past. Systematic searches for low-mass black holes using the Sloan Digital Sky Survey uncovered several hundred new objects residing in low-mass galaxies (Greene & Ho, 2004, 2007; Dong et al., 2012; Reines et al., 2013), but the statistical completeness of these optical searches are difficult to quantify (Greene & Ho, 2009), as AGN samples selected by optical emission lines tend to have black holes radiating at moderate to high fractions of their Eddington limits. On the contrary, X-ray observations, especially of high resolution and high sensitivity afforded by the Chandra Advanced CCD Imaging Spectrometer (ACIS; Weisskopf et al., 2002), can reveal very weak nuclear activity in galaxies, even with brief exposures (ks; e.g., Ho et al., 2001; Gallo et al., 2008; Miller et al., 2012). Desroches & Ho (2009) performed the first dedicated search for low-mass black holes with Chandra and detected 17 candidates from 64 late-type spirals. Later on, cross-matching nearby dwarf galaxies with Chandra and XMM-Newton catalogs of X-ray point-like sources increased the sample size to 93 (Lemons et al., 2015; Pardo et al., 2016; Nucita et al., 2017). Using stacked Chandra images, Mezcua et al. (2016) revealed a population of low-mass black holes that are too faint to be detected individually in dwarf galaxies. The hard X-ray NuSTAR observations of Chen et al. (2017), albeit of a small sample, also suggested a non-negligible fraction of accreting low-mass central black holes that may have been missed in optical surveys.

Here we report the search of low-mass black holes from an even larger sample of nearby late-type galaxies. We have recently conducted an X-ray survey of a large sample of nearby galaxies. The main goals of this survey, described in She et al. (2017, hereafter Paper I), are to systematically search for and study AGNs in all galaxies within 50 Mpc contained in the Chandra archives observed with ACIS. The sample contains 719 galaxies, of which 314 are identified as AGN candidates by the cross-correlation of X-ray point-like sources and the near-infrared or optical stellar nuclei of these galaxies. The technical details of the survey and the X-ray properties of the AGN candidates are presented in Paper I. With 210 keV luminosities less than 10 erg s for most objects, these sources occupy the regime of low-luminosity AGNs (for a review, see Ho, 2008)

This paper focuses on the identification of a population of previously hidden AGNs in the centers of late-type, mostly bulgeless galaxies previously classified as H ii nuclei in optical surveys. Most of them have very low luminosities and are mostly likely associated with low-mass black holes. Unlikely their optically identified counterparts (e.g., Greene & Ho, 2004, 2007; Dong et al., 2012; Reines et al., 2013), this X-ray–selected population, having significantly lower accretion rates, is much more numerous, allowing us to establish a new, firmer estimate of the occupation fraction of low-mass black holes in the nearby universe. For completeness, the paper also addresses the overall demographics of AGNs in the full sample of nearby galaxies, which covers a broad range of Hubble types.

Section 2 gives a brief introduction of the sample and measurements. Section 3 presents the detection rates of X-ray AGN candidates in different kinds of galaxies and their possible dependence on the presence of a large-scale bar. Black hole masses and Eddington ratios are estimated in Section 4. Section 5 highlights the properties of X-ray AGN candidates associated with H ii nuclei. We summarize our results in Section 7.

2. Sample and Data Analysis

Paper I gives a detailed description of the sample and data reduction. Here we just give a brief summary of the most important aspects. The sample is assembled based on the full list of Chandra imaging observations, with both ACIS-S and ACIS-I, as of March 2016. All Chandra-observed galaxies with a distance of less than 50 Mpc were included in our sample, resulting in a total of 719 galaxies.

The Hubble types of the galaxies come from the NASA/IPAC Extragalactic Database (NED),111http://ned.ipac.caltech.edu which are taken largely from the Third Reference Catalogue of Bright Galaxies (RC3; de Vaucouleurs et al. 1991). A total of 285 galaxies have a large-scale bar structure (with Hubble type SAB or SB). Optical nuclear spectral classifications are available for 249 galaxies from the Palomar spectroscopic survey of bright, northern galaxies (Ho et al., 1995, 1997c). Nuclear spectral classifications for another 44 galaxies are found in the catalog of Véron-Cetty & Véron (2010). For the 125 galaxies whose spectral classifications are not included in these two catalogs, we performed our own classification using optical emission-line ratios. Methods of spectral fitting and classification can be found in Paper I. To summarize: 418 out of 719 galaxies in our sample have an optical nuclear spectral classification (for a review see Ho, 2008), including 59 Seyferts, 66 low-ionization nuclear emission-line regions (LINERs; Heckman, 1980), 41 transition objects (emission-line nuclei with [O i] strengths intermediate between those of H ii nuclei and LINERs; see Ho et al. 1993), 163 H ii nuclei, and 89 absorption-line nuclei (those without emission lines).

Point-like X-ray sources are detected using CIAO task wavdetect, and then AGN candidates are identified by detection of an X-ray core astrometrically coincident with the near-infrared/optical position of the galaxy nucleus, taking various errors into account. This results in 314 X-ray AGN candidates, of which 228 are located in the sky plane less than 1 from the near-infrared/optical nucleus. X-ray spectral fitting is conducted for AGN candidates with enough photons; otherwise, hardness ratios are obtained to characterize the shape of the spectrum. The X-ray luminosity in the 210 keV band of each AGN candidate is derived from spectral fitting if available, or from count rates otherwise. For 250 objects, Eddington ratios are calculated based on black holes masses derived from the  relation222The  relation we used in Paper I is . (Kormendy & Ho, 2013, as supplemented in Paper I).

3. Demographics of AGNs

The numbers of galaxies and AGN candidates as a function of Hubble type are shown in Table 1 and in Figure 1 (left), which also displays the X-ray core detection rate with a 90% error range for each type. For early-type galaxies (ESbc), the X-ray core detection rate is as high as 60% (250/434), whereas for late-type galaxies (ScIm), the X-ray core detection rate drops to 25% (54/212). For galaxies with a Hubble type of SmIm, at the extreme end of the late-type sequence, the fraction is less than 10%. The X-ray core detection rate in all galaxies is 43.7%. These AGN detection rates are broadly consistent with those derived from optical spectroscopic observations (Ho et al., 1997d; Ho, 2008), which reported AGN fractions of 50%70% in early-type galaxies and 15% in late-type systems. There appears to be a slight overabundance of X-ray AGN candidates among ScSdm galaxies (see Table 1) compared to optical AGN detection rate found in the Palomar survey. X-ray observations are more sensitive to weak AGNs than optical observations. The results in this work largely confirm the findings of former X-ray investigations (Ho, 2008, 2009; Desroches & Ho, 2009; Zhang et al., 2009; Grier et al., 2011). Compared to previous similar X-ray studies, this work targets a much larger sample of galaxies and utilizes high-quality, high-resolution Chandra data, which were uniformly analyzed. We also perform a more comprehensive analysis on position errors, which is critical to evaluate whether any given compact X-ray source is likely to be associated with the nucleus of the galaxy.

Type All Galaxies Unbarred Galaxies Barred Galaxies
Rate (%) Rate (%) Rate (%)
E 142 80 56
S0 126 70 55 37 26 69 51 27 52
S0/a-Sab 71 47 65 22 17 75 38 24 62
Sb-Sbc 95 53 55 29 16 54 60 35 58
Sc-Scd 100 39 39 34 10 30 58 29 50
Sd-Sdm 52 10 20 11 4 38 37 6 17
Sm-Im 60 5 9 4 2 50 40 2 7
I0 6 4 62
pec 6 1 25
Unknown 61 5 9
All 719 314 43 137 75 54 284 123 43

Note. – The number of galaxies and AGN candidates in different bins of Hubble type, for all galaxies, galaxies without a bar (SA), and galaxies with a bar (SAB and SB). is the number of galaxies, and is the number of AGN candidates. For AGN detection rates, the errors are given at a 90% confidence level.

Table 1AGN numbers and detection rates as a function of Hubble type
Figure 1.— Number of galaxies (open histograms), number of X-ray AGNs (hatched histograms), and the fraction of AGNs (blue dots with 90% errors) as a function of the galaxy Hubble type (left) and optical nuclear spectral classification (right).

Out of the 418 galaxies with reliable optical spectral classifications, 243 (58%) have an X-ray AGN candidate. The X-ray AGN detection rate corresponding to each optical spectral classification is shown in Table 2 and in Figure 1 (right). The frequency of X-ray AGNs is high in optical-selected Seyferts and LINERs, 90% for type 1 (broad-line) and 80% for type 2 (narrow-line) sources. About 70% of the transition objects contain an X-ray nucleus. Pure H ii nuclei have the lowest detection rate of X-ray counterparts, but still %. Approximately 60% of absorption-line nuclei harbor an X-ray AGN candidate. These results are broadly consistent with those in Ho (2008) and Zhang et al. (2009). For all X-ray AGN candidates with an optical spectral classification, 21% are Seyferts (7% type 1 and 14% type 2), 24% are LINERs (7% type 1 and 17% type 2), 12% are transition objects, 21% are H ii nuclei, and 22% are absorption-line nuclei. As a result, more than 40% of X-ray AGN candidates (H ii and absorption-line nuclei) in this work were previously completely missed in optical surveys.

It is worth noting that, most of absorption-line nuclei are early-type galaxies (77/89 are ellipticals and lenticulars), all of which should harbor a central black hole (Kormendy & Ho, 2013). The dearth of cool gas in these galaxies presumably leads to very low mass accretion rates for their central black holes, resulting in very low accretion-powered luminosities and consequently undetected optical line emission, especially in the presence of strong contamination by host galaxy starlight. For the 16 Seyferts and LINERs without a detection of an X-ray core, the most likely culprit is sensitivity. The exposure time is less than 6 ks for most of them (10/16), and for 13 of them the 3 upper limit of the X-ray luminosity in 210 keV is higher than  erg s. For those having H measurements, the 1090 percentile range of their X-ray (upper limit) to H luminosity ratios is 0.210, consistent with 0.3100 for X-ray detected Seyferts and LINERs in the same percentile range. One object, NGC 4625, has a long exposure (55 ks) and a stringent X-ray upper limit of  erg s. We note that for this object, the spectroscopic classification from Véron-Cetty & Véron (2010) may be questionable. In Moustakas et al. (2010), it was classified as a star-forming nucleus based on the [N ii]/H ratio alone. The situation for H ii nuclei or transition objects is difficult to say. As many or most of these objects—especially the H ii nuclei—reside in late-type or even bulgeless galaxies, they may very well lack a central black hole, and hence no X-ray core is expected in them.

Spectral Class Galaxies AGNs Rate (%)
Seyfert 1 19 18 90
Seyfert 2 40 34 83
Seyfert 1 and 2 59 52 86
LINER 1 16 16 94
LINER 2 50 41 80
LINER 1 and 2 66 57 85
Transition 41 28 67
H ii nucleus 163 51 31
Absorption-line 89 55 61
All 418 243 58

Note. – Errors are quoted at 90% confidence level.

Table 2AGN numbers and detection rates as a function of optical spectral classification

3.1. The Influence of Bars

Among disk galaxies, those without a bar have a marginally higher detection rate of AGN candidates compared to barred (SAB or SB) systems: % versus %. This is the exact opposite of what is expected, if bars, as often supposed (e.g., Heller & Shlosman 1994; see a review in Kormendy & Kennicutt 2004), help deliver gas to the center of a galaxy. However, when studying the possible correlations between nuclear activity and a large-scale bar, the dependence of bar frequency on Hubble type should be taken into account, as bars are more prevalent in late-type spirals. To account for this effect, we compare the X-ray AGN detection rates and their uncertainties (given at a 90% confidence level) of barred and unbarred galaxies separated by Hubble type (Figure 2; Table 1). For galaxies with Hubble types S0Sab or SdIm, the X-ray AGN detection rates are slightly higher for unbarred galaxies than barred galaxies, whereas for Hubble types SbScd, unbarred galaxies have marginally lower X-ray AGN rates than barred galaxies. We note that the error bars in Figure 2 are relative large, and any apparent differences in AGN fractions for barred and unbarred are not statistically significant. We find no compelling evidence that bars influence AGN activity, at least not within the context of the current sample of mostly low-luminosity sources in nearby galaxies. Similarly, Ho et al. (1997a) found that bars have a negligible effect on optically-selected AGN activity in the Palomar sample.

Figure 2.— AGN detection rates for unbarred (SA; black) and barred (SAB and SB; blue) galaxies. Errors are quoted at 90% confidence level.

4. Black Hole Mass and Eddington Ratio

Paper I describes our methodology for deriving black hole masses (), X-ray luminosities (; 210 keV), and Eddington ratios (  /; adopted). The statistics of these properties are summarized in Table 3 for the sample as a whole, as well as for various subgroups divided by Hubble type and optical spectral classification. Figure 3 plots the distributions of  and .

Sample  ()  (erg s)
Median Min Max Median Min Max Median Min Max
All 250 7.83 3.45 9.95 314 39.26 35.94 43.11 250 5.18 7.97 0.16
E 75 8.55 4.84 9.61 80 39.08 35.94 41.60 75 6.17 7.96 3.24
S0 64 8.07 5.58 9.95 70 39.27 36.56 42.71 64 5.65 7.57 1.79
S0/a-Sab 35 7.64 6.35 8.80 47 39.50 38.45 43.11 35 4.35 6.67 1.08
Sb-Sbc 40 7.12 5.59 8.55 53 39.98 36.88 42.36 40 4.35 7.97 0.16
Sc-later 31 5.64 3.45 8.10 54 38.92 36.79 42.04 31 3.78 5.32 1.02
Seyfert 44 7.54 4.03 9.16 52 40.61 37.87 43.11 44 3.22 7.39 0.16
LINER 51 8.26 6.39 9.46 57 39.56 38.02 41.60 51 5.37 7.03 3.55
Transition 28 7.84 3.45 8.95 28 39.17 37.56 41.62 28 5.32 6.74 1.69
H ii 35 6.32 3.48 8.69 51 38.87 36.56 41.05 35 4.14 6.03 1.23
Absorption-line 49 8.27 4.84 9.95 55 38.98 35.94 40.91 49 6.34 7.97 3.24

Note. – Black hole masses are estimated from the  relation (see details in Paper I).  is the X-ray luminosity in the 210 keV band. Eddington ratios () are based on the X-ray luminosity and a bolometric correction factor of =16.

Table 3AGN Candidate Statistics
Figure 3.— Distributions of  and  for galaxies binned by Hubble type and spectral classification.

The 210 keV X-ray luminosities of the AGN candidates in this work lie in the range   erg s, with a median value of  erg s (Table 3), which indicates most of these AGN candidates are in the regime of low-luminosity AGNs. The most luminous of the nuclei are found in early-type spirals, whereas galaxies with Hubble types Sc and later have the least luminous nuclei. Seyferts and LINERs are more luminous than H ii nuclei, and  decreases along the sequence Seyferts LINERs transition objects H ii nuclei. This result, drawn entirely from our new X-ray measurements, mirrors the activity sequence proposed by Ho (2009), which was based on H luminosities and a set of largely independent X-ray data.

As previously noted (Ho, 2008, 2009; Gallo et al., 2010; Miller et al., 2012), the nuclei of nearby galaxies generally exhibit an inverse correlation between the Eddington ratio and black hole mass (Figure 3). This is due, in part, to the fact that black holes with higher mass can be detected more easily. The black hole masses in this work cover a very large range, from to   with a median value of . As expected, late-type (Sc and later) galaxies have the smallest black hole masses, while elliptical galaxies have the largest. As the majority of late-type galaxies have a spectral classification (if available) of H ii nuclei, it is not surprising that H ii nuclei have the smallest black hole masses, while LINERs and absorption-line nuclei occupy the upper end of the mass distribution. We note that the uncertainties of the black hole mass estimates, derived from the  relation, can be considerable. Kormendy & Ho (2013) show that classical bulges and ellipticals obey a tight  relation with an intrinsic scatter of only 0.28 dex. By contrast, the  relation of pseudobulges—those with a protracted evolutionary history driven by secular, internal processes—exhibits much larger scatter and a significantly lower zero point compared to classical bulges (see also Greene et al., 2016). A practical difficulty we face is that the vast majority of the galaxies in our sample lack reliable bulge classifications, and hence we cannot use a bulge-type–specific  relation to estimate black hole masses. Under these circumstances, we have little choice but to derive a single  relation for all bulge types (pseudobulges combined with classical bulges and ellipticals), which has a larger intrinsic scatter of 0.44 dex (Paper I). As the typical uncertainty of the bulge stellar velocity dispersions is only 5%–15% (Ho et al., 2009), the final uncertainties of the black hole masses are dominated by the intrinsic scatter of the  relation. Moreover, because the zero point of the  relation of pseudobulges is lower than that of classical bulges and ellipticals, using a single, combined  relation will tend to systematically overestimate the black hole masses for pseudobulges and bulgeless galaxies, while ellipticals and classical bulges will have their black hole masses systematically underestimated. Conversely, Eddington ratios will be underestimated in pseudobulges but overestimated in classical bulges and ellipticals.

The Eddington ratios spread over an even larger dynamic range, from to , with a median value of . Our sample contains predominantly black holes accreting at extremely sub-Eddington rates. The AGN candidates in ellipticals, having the most massive black holes, exhibit the lowest values of , whereas late-type galaxies (Sc and later) have relatively higher  as a result of their lower black hole masses. Among objects traditionally recognized to host black holes, at a given black hole mass, Seyferts, as a class, have the highest Eddington ratios (median ), followed by LINERs and transition objects (median ), and then absorption-line nuclei (median ). The AGN candidates newly identified in H ii nuclei, with median , lie in between Seyferts and LINERs/transition objects.

Apart from the new population of AGN candidates in H ii nuclei, the results above are broadly consistent with those previously obtained from optical observations (Ho et al., 1997d, 2003). Ho (2009), using a compilation of more heterogeneous and generally less sensitive X-ray measurements from various literature sources, also derived very similar statistical results for the X-ray luminosities and Eddington ratios of the galaxies in the Palomar survey. This investigation, however, represents a significant improvement in sample size and data quality, as all of the X-ray data come from Chandra/ACIS observations. It is clear that—with sufficiently high sensitivity and angular resolution—X-ray observations allow us to detect nuclear accretion at much lower levels of activity than possible by optical techniques. In this work, for AGN candidates in H ii and absorption-line nuclei, the median value of X-ray luminosity is  erg s, which corresponds to  erg s,333Low-luminosity AGNs have a bolometric correction of 220 16 (Ho, 2008) about an order of magnitude lower than the median  for the Palomar objects (Ho et al., 1997d,  erg s). In addition to sheer sensitivity, Chandra/ACIS observations, which have the unique capability of “filtering out” most of the contamination from host galaxy starlight, also afford more reliable detections of faint AGNs, which can be easily masked or confused by dust and star-forming regions at other wavelengths. This is well illustrated in Figure 4 (see also Ho et al., 2001). The next section highlights a population of AGN candidates that was previously hidden among optically classified H ii nuclei.

5. AGNs Candidates in H ii Nuclei

H ii nuclei, those whose dominant source of excitation derives from photoionization by young, massive stars, are extremely prevalent. The Palomar optical spectroscopic survey of nearby galaxies found that 40% of all nearby galaxies host an H ii nucleus, and for late-type galaxies the fraction is even higher, reaching 80% (Ho et al., 1997c, d). Compared to H ii regions in galactic disks, however, H ii nuclei exhibit systematically stronger low-ionization optical forbidden lines, suggesting that they contain an extra source of excitation not present in normal star-forming regions (Kennicutt et al., 1989; Ho et al., 1997b). While Kennicutt et al. (1989) suggested that the additional excitation may come from a hidden AGN, Ho et al. (1997b) argued that shocks are more likely, although they noted that X-ray observations would be effective in uncovering optically hidden AGNs, should they be present.

This work confirms that a significant fraction—31% (51/163)—of optically classified H ii nuclei have an X-ray core. The general and X-ray properties of this special subsample of AGN candidates are listed in Table 4. To highlight the morphology of the host galaxies, we assembled the highest available resolution optical ( band if available, other bands otherwise) images from NED. These are displayed in Figure 4, superposed with the Chandra/ACIS view of the central 1515 area to illustrate the clear association of the compact X-ray core with the stellar nucleus. The very late-type morphologies of the host galaxies are immediately obvious, as is the fact that many of them clearly lack a bulge altogether. Figure 5 examines the Veilleux & Osterbrock (1987) optical line-intensity ratio diagnostic diagrams for the sample of H ii nuclei included in this study. It is interesting that the subset with detected X-ray cores do not show preferentially enhanced low-ionization lines of [N ii], [S ii], or [O i]. Thus, while some H ii nuclei do harbor optically hidden AGNs, they are sufficiently weak that they do not contribute significantly to the ionization budget of the optical line emission, and, in particular, they do not appear to be the primary cause for the excitation enhancement reported by Kennicutt et al. (1989) and Ho et al. (1997b).

Name Hubble Type Model
() (erg s) ( cm) (erg s)
IC 342 6.07 38.11 PL+VME 1.92 38.38
M 61 6.35 39.93 PL+ME 1.53 38.91
M 101 3.48 39.14 PL 1.88 38.07
NGC 45 1.60 37.22
NGC 891 6.04 38.85L 9.04 38.38
NGC 922 0.24 39.02
NGC 925 6.00 38.36 0.45 38.75
NGC 1073 3.60 38.36 0.57 39.27
NGC 1370 5.96 0.52 38.73
NGC 1493 0.32 38.41
NGC 1603 7.03 0.13 38.99
NGC 1637 PL 3.01 38.27
NGC 1808 7.64 PL+VME 0.74 39.50
NGC 2276 6.34 40.83 0.02 38.87
NGC 2500 5.04 38.00 0.83 38.17
NGC 2748 6.32 39.99 0.35 38.32
NGC 2782 8.11 41.51 PABS*PL+ME+G 0.55 54.27 41.05
NGC 2798 7.01 PL+ME 2.62 39.56
NGC 2993 PL+ME 1.37 0.14 40.03
NGC 3077 4.20 38.93 PL 1.69 1.42 38.22
NGC 3184 4.85 39.72 0.54 37.94
NGC 3198 5.00 39.62 0.95 38.76
NGC 3310 6.35 40.53 PL 1.41 0.44 40.25
NGC 3367 5.64 41.34 PL+ME 0.98 40.71
NGC 3521 7.35 39.10b PL 1.92 0.64 38.92
NGC 3665 8.69 39.75b 0.75 39.70
NGC 3877 6.41 40.29 0.11 37.99
NGC 4038 0.06 38.22
NGC 4102 8.00 41.50 PL+ME+G 1.52 40.47
NGC 4136 4.58 38.30 1.01 38.14
NGC 4194 6.83 PL+VME 1.44 0.66 40.32
NGC 4217 6.54 38.84 2.59 38.45
NGC 4490 4.95 38.83 PL 2.40 0.92 39.63
NGC 4526 8.45 39.61 PL 1.10 0.13 39.33
NGC 4559 5.14 38.58 0.59 38.23
NGC 4561 PL 1.45 39.56
NGC 4654 5.07 39.97 1.49 39.44
NGC 4666 1.09 39.04
NGC 4670 0.76 39.28
NGC 4900 5.54 39.48L 0.62 38.52
NGC 5102 0.68 36.56
NGC 5248 6.73 39.76 0.77 38.15
NGC 5253 0.01 36.90
NGC 5483 1.14 39.03
NGC 5775 6.75 39.54 1.00 38.62
NGC 6764 PL+ME 1.82 40.05
NGC 6946 5.43 40.14 PL 2.23 0.58 38.83
NGC 7320 0.86 38.40
NGC 7714 7.14 PL 1.62 0.16 40.55
UGC 5720 0.01 39.14
UGCA 166 PL 1.72 39.34

Note. – Column 1: Galaxy name. Column 2: Hubble type from NED. Column 3: Log of black hole masses estimated using the  relation. More details are in Paper I. Column 4: Extinction corrected H luminosity from the Palomar survey (Ho et al., 1997c, 2003). The typical uncertianties is 10-30%, and the letters “b”,“c”,“L” and “u” indicate an uncertainty of 30-50%, 100%, a 3 lower limit, and a 3 upper limit, respectively. Column 5: Spectral fitting models, see Paper I for details. Here “PL”,“ME”,“VME”,“PABS”,“BB”,“G” denote for a powerlaw, mekal, vmekal, pcfabs, bbody, and gaussian (for iron K line) component in XSPEC, respectively. Column 6: Photon index . Column 7: Intrinsic absorption from spectral fitting or hardness ratio otherwise. Column 8: X-ray luminosity in 2-10 keV band.

Table 4Properties of AGN Candidates in H ii Nuclei
Figure 4.— Optical and Chandra/ACIS images (upper right corner) of AGN candidates in H ii nuclei. The horizontal segment in the bottom left corner of each panel represents a scale of from RC3, and each X-ray image subtends , smoothed with a Gaussian of pixel in ds9. The cross marks the near-infrared/optical nuclear position, and the green ellipse in the X-ray image is the source region given by wavdetect. North is up and east is to the left.

Fig. 4 (Continued)

Fig. 4 (Continued)

Figure 5.— Optical emission-line flux ratio diagnostic diagrams for the sample of H ii nuclei in this study (upper panels) and histogram of strength of low-ionization lines (lower panels).

Although we are cautious to call these “candidate” AGNs, we have good reason to believe that this is the correct interpretation. With a median 210 keV luminosity of  erg s, 90% of the AGN candidates in H ii nuclei exceeds  erg s. X-ray binaries associated with a potential nuclear star cluster may contaminate our AGN candidates. On the one hand, the luminosity function of low-mass X-ray binaries (LMXBs) shows a sharp cutoff at a few times  erg s and will have few sources above  erg s (Gilfanov, 2004), suggesting that most of the sources above  erg s cannot be LMXBs. On the other hand, the typical age of  yr for nuclear star clusters (Rossa et al., 2006) is too old to form high mass X-ray binaries (HMXBs), which could reach  erg s (Swartz et al., 2011), unless there is ongoing star formation. For the latter reason, Desroches & Ho (2009) argued that the HMXB contamination to the nuclear X-ray sources is negligible, and only the contribution from LMXBs is considered. Here, to be conservative, we estimate both contributions to the detected X-ray cores.

For HMXBs, their luminosity function scales with the star formation rate (Mineo et al., 2012), which itself depends on the H luminosity,   (erg s) M yr(Kennicutt & Evans, 2012). Out of the 163 H ii nuclei, 90 of them have  measurements. For galaxies in Mineo et al. (2012), H ii regions contribute 55% of the total  emission (Greenawalt et al., 1998; Hoopes et al., 1999; Thilker et al., 2002), and only 10% of HMXBs (in 10 erg s) are coincident with H ii regions (with position offsets less than ). As the H ii nuclei have H luminosities not too different from that emitted by other H ii regions in the optical disk of these galaxies, we use /0.55 to correct for the total H luminosity in the whole galaxy, and renormalize the HMXB luminosity function by a factor of 0.1 to account for their association with the nuclear H ii region.

For LMXBs, their luminosity function is scaled with the total stellar mass. We treat nuclear clusters as globular clusters, as they share remarkably similar properties (Walcher et al., 2005). It is known that the LMXB production rate per unit stellar mass in globular clusters is roughly 100 times more efficient than in the field (Clark, 1975; Peacock et al., 2017). Moreover, metal-rich globular clusters are known to produce more LMXBs than metal-poor ones, roughly by a factor of 3 (Kundu et al., 2007; Kim et al., 2013). The metallicities of nuclear clusters in late-type galaxies are found to be similar to those of metal-rich clusters (Brodie & Strader, 2006; Rossa et al., 2006). Thus, we take the LMXB luminosity function in globular clusters from Kim et al. (2009) and the normalization from Peacock et al. (2017) as 0.11 LMXB with  erg s per   for metal-rich globular clusters. The typical mass for the nuclear clusters in our sample ( ) is estimated from the 39 objects that are also shown in Georgiev et al. (2016). The expected luminosity functions for both HMXBs and LMXBs are shown in Figure 6.

Figure 6.—  distributions of AGN candidates in H ii nuclei and expected X-ray luminosity function of HMXBs and LMXBs.

Those constitute strong evidence that, in general, 40% of the X-ray cores may be contaminated by X-ray binaries in the luminosity range of  erg s, and at luminosities above  erg s the contamination drops to 8%. The total contamination at luminosities of  erg s or higher is 26%. Thus, the majority of the X-ray cores should be previously unrecognized low-luminosity AGNs, especially at luminosities above  erg s.

The ratios of / for the AGN candidates in H ii nuclei (with a median value of 0.3; see Figure 7) are considerably lower than those of other AGN candidates in Seyferts, LINERs or transition objects, which have a median / . A simple explanation for this order-of-magnitude difference in / is that the majority (%) of the H emission from these H ii nuclei—as expected, given their optical classification—does not come from AGN photoionization but from star-forming activity. The integrated optical emission from star-forming regions vastly overwhelm the weak AGN signal from these relatively low-mass central black holes.444Most of the optical classifications derive from observations from the Palomar (Ho et al., 1995) or 6dF (Jones et al., 2009) surveys, which employed an aperture size much larger than that of the X-ray observations. The Palomar spectra were mostly extracted using a aperture, and the 6dF survey used a fiber of diameter 6.7; typical radius of the Chandra/ACIS point-spread function is for Chandra/ACIS.

An important advantage of our study is that a sizable portion of our sample has enough counts to enable X-ray spectral analysis, and hence we can evaluate whether the X-ray sources have spectral properties consistent with other know low-luminosity AGNs (e.g., Ho, 2008). Of the 21 (41%) of AGN candidates associated with H ii nuclei for which we could analyze using XSPEC, 11 can be well fitted with a single, absorbed power-law, nine requires a thermal component (mekal), and one (NGC 2782) requires a partially covered absorption component (pcfabs). Except for NGC 2782, 60% of the photon indices are in the range 1.42.0 with a typical value of of 1.7, which is only marginally flatter than the typical value for low-luminosity AGNs (Ho, 2008). Two sources (NGC 2782 and NGC 4102) show an iron K emission line.

Lastly, timing analysis can, in principle, yield additional constraints on the physical origin of the X-ray emission. Unfortunately, due to the low count rates and short durations of most of the observations, the evidence based on timing is inconclusive. We searched for variability on short timescales for the 11 sources in the sample with more than 500 photons in the 0.38 keV band. We analyzed their light curves with a time bin of 2000 s and calculated the for each light curve by assuming a constant count rate. Based on these light curves, the normalized excess variance (; Vaughan et al., 2003) are also calculated. Only one object, NGC 3310, shows significant short timescale variability (=76/24, ); for other objects, only upper limits of can be given. We also checked whether our Chandra data can detect variability in some of the brighter,well-known AGNs in our full sample (e.g., NGC 4261, NGC 5506, NGC 4945, NGC 4051), whose short timescale variability has been previously studied in. Some vary, but others do not. We note that the count rate and the total duration of the observations are crucial for searching for short timescale variabilities. The vast majority of the archival Chandra data we used simply cannot yield meaningful constraints on X-ray variability.

Figure 7.— Correlation of X-ray luminosity and extinction-corrected H luminosity. The black open circles denote X-ray AGN candidates, with upper limits marked by arrows. Among the X-ray–detected galaxies, those classified optically as H ii nuclei are marked with green squares. The dashed lines represent /=100 and 1. The typical uncertainties of  and  are plot at the bottom right corner.

6. Implications for Low-mass Central Black Holes

Most of the host galaxies of AGN candidates in H ii nuclei have late types (28 out of the 51 have Hubble types ScIm), whose bulge component, if it exists at all, is certainly not dominant. The AGN candidates in these H ii nuclei have the lowest median value of  ( ) compared to other spectral classes (Seyferts, LINERs, transition objects and absorption-line nuclei; see Table 3). As they occupy the low-mass end of the mass distribution of central black hole, these AGN candidates in H ii nuclei are very important for understanding the overall demography of “seed” supermassive black holes. Previous systematic searches for low-mass ( ) black holes based on optical selection (Greene & Ho, 2004, 2007; Barth et al., 2008; Dong et al., 2012; Reines et al., 2013; Moran et al., 2014) were only sensitive to sources radiating at a significant fraction of the Eddington limit. It is highly uncertain whether these optically selected samples give an accurate census of the incidence of low-mass black holes in the nearby universe. This study, following the initial suggestion by Desroches & Ho (2009) based on a Chandra study of a much more limited sample of late-type galaxies, confirms that the true frequency of low-mass black holes, dominated by objects accreting at sub-Eddington rates, is significantly higher. Among the 51 detected X-ray H ii nuclei, 46 of them have an X-ray luminosity above  erg s. Taking into account the contamination from X-ray binaries (26% above  erg s), we regard 21% as a reasonable lower limit on the true incidence of AGNs in late-type galaxies, as the contamination by X-ray binaries (both LMXBs and HMXBs) very likely has been overestimated, and some low-luminosity AGNs are not detected due to our sensitivity limit.

7. Conclusions

We identified 314 AGN candidates out of 719 galaxies within 50 Mpc (Paper I). In this work, we study the demographics of AGNs in nearby galaxies, and reveal 50 hidden AGN candidates in H ii nuclei. The main conclusions are summarized below.

  1. In early-type galaxies, about 60%–70% of them are found to have an AGN candidate, and this fraction drops to less than 20% for very late-type galaxies (Sdm and Sm). About 85% of Seyferts and LINERs show an X-ray core. This fraction is 70% for transition objects and 60% for absorption-line nuclei. 30% of H ii nuclei have an X-ray core. We do not find a significant correlation between the presence of a large-scale bar and AGN activity.

  2. The X-ray luminosity decreases in a sequence of Seyferts LINERs transition objects H ii nuclei. Late-type spirals have the lowest luminosities.

  3. In general, more massive black holes tend to display a smaller Eddington ratio. At a given black hole mass, the Seyferts have the highest Eddington level.

  4. In 163 H ii nuclei, 51 out of them are identified to have an X-ray core. Most of these X-ray sources are not X-ray binaries. The morphology of these galaxies is of late-type and most likely has no or very small bulges. These objects are candidate low-mass black holes according to the  relation, suggesting that the black hole occupation fraction in late-type galaxies is at least 21%.

The research of LCH is supported by the National Key Program for Science and Technology Research and Development (2016YFA0400702) and the National Science Foundation of China (11473002, 11303008). HF acknowledges funding support from the National Natural Science Foundation of China under grant No. 11633003, and the National Program on Key Research and Development Project (Grant No. 2016YFA040080X).


  • Barth et al. (2008) Barth, A. J., Greene, J. E., & Ho, L. C. 2008, AJ, 136, 1179
  • Brodie & Strader (2006) Brodie, J. P., & Strader, J. 2006, ARA&A, 44, 193
  • Chen et al. (2017) Chen, C.-T. J., Brandt, W. N., Reines, A. E., et al. 2017, ApJ, 837, 48
  • Clark (1975) Clark, G. W. 1975, ApJ, 199, L143
  • de Vaucouleurs et al. (1991) de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., et al. 1991, Third Reference Catalogue of Bright Galaxies. (Springer)
  • Desroches & Ho (2009) Desroches, L.-B., & Ho, L. C. 2009, ApJ, 690, 267
  • Dong et al. (2012) Dong, X.-B., Ho, L. C., Yuan, W., et al. 2012, ApJ, 755, 167
  • Gallo et al. (2008) Gallo, E., Treu, T., Jacob, J., et al. 2008, ApJ, 680, 154
  • Gallo et al. (2010) Gallo, E., Treu, T., Marshall, P. J., et al. 2010, ApJ, 714, 25
  • Georgiev et al. (2016) Georgiev, I. Y., Böker, T., Leigh, N., Lützgendorf, N., & Neumayer, N. 2016, MNRAS, 457, 2122
  • Gilfanov (2004) Gilfanov, M. 2004, MNRAS, 349, 146
  • Greenawalt et al. (1998) Greenawalt, B., Walterbos, R. A. M., Thilker, D., & Hoopes, C. G. 1998, ApJ, 506, 135
  • Greene (2012) Greene, J. E. 2012, NatCo, 3, 1304
  • Greene & Ho (2004) Greene, J. E., & Ho, L. C. 2004, ApJ, 610, 722
  • Greene & Ho (2007) Greene, J. E., & Ho, L. C. 2007, ApJ, 670, 92
  • Greene & Ho (2009) Greene, J. E., & Ho, L. C. 2009, PASP, 121, 1167
  • Greene et al. (2016) Greene, J. E., Seth, A., Kim, M., et al. 2016, ApJ, 826, L32
  • Grier et al. (2011) Grier, C. J., Mathur, S., Ghosh, H., & Ferrarese, L. 2011, ApJ, 731, 60
  • Haehnelt & Rees (1993) Haehnelt, M. G., & Rees, M. J. 1993, MNRAS, 263, 168
  • Heckman (1980) Heckman, T. M. 1980, A&A, 87, 152
  • Heger et al. (2003) Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003, ApJ, 591, 288
  • Heller & Shlosman (1994) Heller, C. H., & Shlosman, I. 1994, ApJ, 424, 84
  • Ho (2008) Ho, L. C. 2008, ARA&A, 46, 475
  • Ho (2009) Ho, L. C. 2009, ApJ, 699, 626
  • Ho et al. (1993) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1993, ApJ, 417, 63
  • Ho et al. (1995) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1995, ApJS, 98, 477
  • Ho et al. (1997a) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997a, ApJ, 487, 591
  • Ho et al. (1997b) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997b, ApJ, 487, 579
  • Ho et al. (1997c) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997c, ApJS, 112, 315
  • Ho et al. (1997d) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997d, ApJ, 487, 568
  • Ho et al. (2003) Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 2003, ApJ, 583, 159
  • Ho et al. (2009) Ho, L. C., Greene, J. E., Filippenko, A. V., & Sargent, W. L. W. 2009, ApJS, 183, 1
  • Ho et al. (2001) Ho, L. C., Feigelson, E. D., Townsley, L. K., et al. 2001, ApJ, 549, L51
  • Hoopes et al. (1999) Hoopes, C. G., Walterbos, R. A. M., & Rand, R. J. 1999, ApJ, 522, 669
  • Jones et al. (2009) Jones, D. H., Read, M. A., Saunders, W., et al. 2009, MNRAS, 399, 683
  • Kennicutt & Evans (2012) Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531
  • Kennicutt et al. (1989) Kennicutt, Jr., R. C., Keel, W. C., & Blaha, C. A. 1989, AJ, 97, 1022
  • Kim et al. (2013) Kim, D.-W., Fabbiano, G., Ivanova, N., et al. 2013, ApJ, 764, 98
  • Kim et al. (2009) Kim, D.-W., Fabbiano, G., Brassington, N. J., et al. 2009, ApJ, 703, 829
  • Kormendy & Ho (2013) Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511
  • Kormendy & Kennicutt (2004) Kormendy, J., & Kennicutt, Jr., R. C. 2004, ARA&A, 42, 603
  • Kundu et al. (2007) Kundu, A., Maccarone, T. J., & Zepf, S. E. 2007, ApJ, 662, 525
  • Lemons et al. (2015) Lemons, S. M., Reines, A. E., Plotkin, R. M., Gallo, E., & Greene, J. E. 2015, ApJ, 805, 12
  • Mezcua et al. (2016) Mezcua, M., Civano, F., Fabbiano, G., Miyaji, T., & Marchesi, S. 2016, ApJ, 817, 20
  • Miller et al. (2012) Miller, B., Gallo, E., Treu, T., & Woo, J.-H. 2012, ApJ, 747, 57
  • Mineo et al. (2012) Mineo, S., Gilfanov, M., & Sunyaev, R. 2012, MNRAS, 419, 2095
  • Moran et al. (2014) Moran, E. C., Shahinyan, K., Sugarman, H. R., Vélez, D. O., & Eracleous, M. 2014, AJ, 148, 136
  • Mortlock et al. (2011) Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, Nature, 474, 616
  • Moustakas et al. (2010) Moustakas, J., Kennicutt, Jr., R. C., Tremonti, C. A., et al. 2010, ApJS, 190, 233
  • Nucita et al. (2017) Nucita, A. A., Manni, L., De Paolis, F., Giordano, M., & Ingrosso, G. 2017, ApJ, 837, 66
  • Pardo et al. (2016) Pardo, K., Goulding, A. D., Greene, J. E., et al. 2016, ApJ, 831, 203
  • Peacock et al. (2017) Peacock, M. B., Zepf, S. E., Kundu, A., et al. 2017, MNRAS, 466, 4021
  • Reines & Comastri (2016) Reines, A. E., & Comastri, A. 2016, PASA, 33, e054
  • Reines et al. (2013) Reines, A. E., Greene, J. E., & Geha, M. 2013, ApJ, 775, 116
  • Rossa et al. (2006) Rossa, J., van der Marel, R. P., Böker, T., et al. 2006, AJ, 132, 1074
  • She et al. (2017) She, R., Ho, L. C., & Feng, H. 2017, ApJ, 835, 223
  • Swartz et al. (2011) Swartz, D. A., Soria, R., Tennant, A. F., & Yukita, M. 2011, ApJ, 741, 49
  • Thilker et al. (2002) Thilker, D. A., Walterbos, R. A. M., Braun, R., & Hoopes, C. G. 2002, AJ, 124, 3118
  • Vaughan et al. (2003) Vaughan, S., Edelson, R., Warwick, R. S., & Uttley, P. 2003, MNRAS, 345, 1271
  • Veilleux & Osterbrock (1987) Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295
  • Véron-Cetty & Véron (2010) Véron-Cetty, M.-P., & Véron, P. 2010, A&A, 518, A10
  • Volonteri (2010) Volonteri, M. 2010, A&A Review, 18, 279
  • Volonteri et al. (2008) Volonteri, M., Lodato, G., & Natarajan, P. 2008, MNRAS, 383, 1079
  • Walcher et al. (2005) Walcher, C. J., van der Marel, R. P., McLaughlin, D., et al. 2005, ApJ, 618, 237
  • Weisskopf et al. (2002) Weisskopf, M. C., Brinkman, B., Canizares, C., et al. 2002, PASP, 114, 1
  • Wu et al. (2015) Wu, X.-B., Wang, F., Fan, X., et al. 2015, Nature, 518, 512
  • Zhang et al. (2009) Zhang, W. M., Soria, R., Zhang, S. N., Swartz, D. A., & Liu, J. F. 2009, ApJ, 699, 281
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description