A Uniformly Selected Sample of Low-mass Black Holes in Seyfert 1 Galaxies. II. The SDSS DR7 Sample

, , , , and

Published 2018 April 19 © 2018. The American Astronomical Society. All rights reserved.
, , Citation He-Yang Liu et al 2018 ApJS 235 40 DOI 10.3847/1538-4365/aab88e

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0067-0049/235/2/40

Abstract

A new sample of 204 low-mass black holes (LMBHs) in active galactic nuclei (AGNs) is presented with black hole masses in the range of (1–20) × 105 M. The AGNs are selected through a systematic search among galaxies in the Seventh Data Release (DR7) of the Sloan Digital Sky Survey (SDSS), and careful analyses of their optical spectra and precise measurement of spectral parameters. Combining them with our previous sample selected from SDSS DR4 makes it the largest LMBH sample so far, totaling over 500 objects. Some of the statistical properties of the combined LMBH AGN sample are briefly discussed in the context of exploring the low-mass end of the AGN population. Their X-ray luminosities follow the extension of the previously known correlation with the [O iii] luminosity. The effective optical-to-X-ray spectral indices αOX, albeit with a large scatter, are broadly consistent with the extension of the relation with the near-UV luminosity L2500 Å. Interestingly, a correlation of αOX with black hole mass is also found, with αOX being statistically flatter (stronger X-ray relative to optical) for lower black hole masses. Only 26 objects, mostly radio loud, were detected in radio at 20 cm in the FIRST survey, giving a radio-loud fraction of 4%. The host galaxies of LMBHs have stellar masses in the range of 108.8–1012.4 M and optical colors typical of Sbc spirals. They are dominated by young stellar populations that seem to have undergone continuous star formation history.

Export citation and abstract BibTeX RIS

1. Introduction

Mounting evidence suggesting that many massive galaxies harbor supermassive black holes (SMBHs) with masses ranging from 106 M to 109 M at their centers in the present universe (e.g., Richstone et al. 1998; Kormendy & Ho 2013) has been accumulated. The mass of SMBHs correlates tightly with the parameters of the host bulges, such as mass, luminosity, and velocity dispersion, which hints at the co-evolution of SMBHs and their host galaxies (e.g., Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000a; Gültekin et al. 2009). How black holes (BHs) formed across the cosmic time remains far from being understood. It is commonly believed that SMBHs have grown from seed BHs mainly through merger-induced accretion (Volonteri 2010; Greene 2012). Secular processes also play a role in BH growth, especially for BHs located in low-redshift galaxies (Greene et al. 2008; Jiang et al. 2011; Conselice et al. 2014). However, even less is known regarding the properties of seed BHs, such as their origins, initial masses, and environments, simply because direct observations of seed BHs in the early universe are not feasible with current facilities. As an indirect approach, low-mass black holes (LMBHs) with masses of thousands to hundreds of thousands of solar masses, residing at the center of galaxies, may help give insight into the formation and evolution of the first seed BHs. Such BHs are also termed intermediate-mass black holes (IMBHs) in the literature. We refer to them as LMBHs to avoid confusion with ultraluminous X-ray sources (ULXs), which are off-nucleus point-like sources, some of which may be powered by IMBHs (Kaaret et al. 2017). LMBHs found in nearby dwarf stellar systems in the present universe seem to have formed early, but have not fully grown into SMBHs (Greene 2012). They may thus trace the seed BHs or the early phase of their growth. In addition, the study of the occupation fraction of present-day LMBHs can help discriminate among various seed BH models, such as light seeds as the end products of Population III stars or heavy seeds formed from the direct collapse of halos at high redshifts (e.g., Lodato & Natarajan 2006; Volonteri et al. 2008; Volonteri & Natarajan 2009). Interestingly, mergers of binary BHs with masses of the order of 105 M can produce gravitational-wave signatures (e.g., Hughes 2002) that will be the primary targets of LISA7 in the future. The detection of gravitational waves in these frequency regimes has become more promising after the prototypical event GW 150914, the first-ever gravitational-wave source, was detected by the Laser Interferometer Gravitational-wave Observatory (LIGO; e.g., Abbott et al. 2016).

The most direct and secure way to detect BHs at the center of galaxies is to seek the effect of their gravitational influence on the spatially resolved dynamics of closely surrounding gases and/or stars (e.g., Barth et al. 2001; Ghez et al. 2008). However, for LMBHs, this method is infeasible for distant galaxies since the BH gravitational sphere of influence can only be spatially resolved within the Local Group with existing facilities. A common practice is to search for accretion-powered radiative signatures of active galactic nuclei (AGNs) hosting such LMBHs. For AGNs exhibiting broad emission lines in their optical spectra, their BH masses can be estimated using MBH = fRυ2/G, assuming that the broad-line region (BLR) system is virialized and individual clouds are moving in Keplerian orbits. G is the gravitational constant. The velocity dispersion υ can be measured from the widths of the broad emission lines. The BLR radius R is estimated from the luminosity of the AGN continuum emission using the radius–luminosity relation derived from reverberation mapping studies of AGNs (e.g., Kaspi et al. 2005; Bentz et al. 2006, 2009, 2013). f is a scaling factor of the order of unity, depending on the distribution and inclination of the BLR cloud orbits to the line of sight, and can be calibrated from other independent measurements of the same BH systems. BH masses thus estimated are indirect and subject to systematics as large as ∼0.3 dex (e.g., Gebhardt et al. 2000b; Greene & Ho 2006; Grier et al. 2013). Nevertheless, they provide us with one of the fundamental parameters to study the vast majority of the AGN population. An archetypal LMBH AGN is NGC 4395 (Filippenko & Ho 2003), which is a bulgeless galaxy harboring a central BH of MBH ∼ 3 × 105 M measured through reverberation mapping observations of the C iv λ1549 line (Peterson et al. 2005). Multiwavelength observations in the X-ray (e.g., Iwasawa et al. 2000; Moran et al. 2005; Dewangan et al. 2008) and radio bands (e.g., Wrobel & Ho 2006) support the finding of a small BH mass in NGC 4395. Other convincing cases of LMBHs in AGNs include POX 52, a dwarf spheroidal galaxy with a BH of MBH ∼ 3 × 105 M (e.g., Barth et al. 2004; Thornton et al. 2008), and SDSS J160531.84+174826.1, a dwarf disk galaxy with a BH of MBH ∼ 7 × 104 M (Dong et al. 2007).

In general, a single-epoch optical spectrum is an effective tool for detecting the signature of nuclear BHs in active galaxies. AGNS with LMBHs as large samples were first searched by Greene & Ho (2004, 2007a) from the Sloan Digital Sky Survey (SDSS; York et al. 2000), yielding ∼200 broad-line sources with BH masses in the range of ∼105–106 M, estimated from the widths and luminosities of the broad Hα emission lines. Independently, we (Dong et al. 2012, hereafter Dong et al. 2012b) carried out a systematic and homogeneous search for LMBHs from the SDSS Fourth Data Release (DR4; Adelman-McCarthy et al. 2006), resulting in 309 AGNs with relatively lower BH masses and Eddington ratios. In recent surveys targeting low-mass type 2 AGNs, there were more than 100 sources in nearby dwarf galaxies identified as LMBH candidates using narrow-line diagnostics (e.g., Barth et al. 2008; Reines et al. 2013; Moran et al. 2014). According to Yuan et al. (2014), if the intrinsic Eddington ratio distribution of SMBHs obtained by Schulze & Wisotzki (2010) can be applied to the low-mass end, there should exist many more LMBHs in the universe. Clearly, an LMBH sample larger than the current ones homogeneously selected with well-understood completeness is essential to construct the intrinsic Eddington ratio and mass functions (W. J. Liu et al. 2018, in preparation). These will give more stringent constraints on seed BH models. Moreover, a statistical study of the interplay of LMBH AGNs and their host galaxies also requires a larger sample, especially at the lower-MBH and lower-accretion-rate end. In this work, we perform an extended search for more low-mass AGNs from the SDSS Seventh Data Release (DR7; Abazajian et al. 2009), based on our previous work using SDSS DR4 (Dong et al. 2012b).

It is difficult to search for AGNs with LMBHs since their optical spectra are mostly strongly dominated by starlight. Thus, careful subtraction of the starlight is essential for reliable measurement of the emission lines. In addition, the decomposition of the broad and the narrow components of the Balmer lines in LMBH AGNs is difficult as the broad lines in the spectra of LMBHs are relatively narrow and weak.8 In order to search for LMBHs, we have designed a set of elaborate codes and broad-line selection procedures as described in Dong et al. (2012b). In this work, we follow exactly the same method and apply it to SDSS DR7. Some 204 new LMBHs are found with MBH < 2 × 106 M (in accordance with Greene & Ho 2007a and Dong et al. 2012b), expanding the total SDSS sample of LMBHs to 513. The BH masses of the new sample range from 1 × 105 M to 2 × 106 M, with a median of 1.3 × 106 M, and the Eddington ratios range from ∼0.01 to ∼1. The data analysis and sample selection are outlined in Section 2, and the LMBH sample is described in Section 3. In Section 4, the sample properties are discussed, followed by a summary in Section 5. Throughout the paper, we assume a cosmology with H0 = 70 km s−1 Mpc−1, Ωm = 0.3, and ΩΛ = 0.7.

2. Data Analysis and Sample Selection

Our LMBH AGNs are selected following the data analysis procedures described in Dong et al. (2012b), which are only briefly summarized here (see Dong et al. 2012b for details). We start with the SDSS DR7 spectra classified as "galaxies" or "QSOs" by the SDSS pipeline, excluding objects that are also in SDSS DR4. The SDSS is a comprehensive imaging and spectroscopic survey using a dedicated 2.5 m telescope (Gunn et al. 2006) located at the Apache Point Observatory in Southern New Mexico. It utilizes a wide-field imager (Gunn et al. 1998), covering the sky in a drift-scan mode in five filters, ugriz (Fukugita et al. 1996), and a 640 fiber-fed pair of multi-object double spectrographs covering the wavelength 3800–9200 Å with a resolution of λλ ≈ 2000. The diameter of the optical fibers is 3'', and the instrumental dispersion is ∼69 km s−1 pixel−1. To ensure that the Hα is lying within the wavelength coverage range of the SDSS spectra, only sources with redshifts below 0.35 are considered. These result in a parent sample consisting of 337,988 "galaxies" and 4697 "QSOs." The spectra are corrected for Galactic extinction using the extinction map of Schlegel et al. (1998) and the reddening curve of Fitzpatrick (1999), and are then transformed to the rest frame using the redshifts provided by the SDSS pipeline.

The spectra of the parent sample are mostly dominated by starlight. As the first step, a pseudo-continuum is modeled and subtracted using the technique described in Zhou et al. (2006).9 The so-called pseudo-continuum is a linear combination of several components, including starlight, nuclear continuum, and the optical Fe ii multiplets. The Balmer continua and high-order Balmer emission lines are added if they can improve the reduced χ2 by at least 20%. The stellar component is modeled by six synthesized galaxy spectral templates built up from the library of simple stellar populations (Bruzual & Charlot 2003) using the Ensemble Learning Independent Component Analysis algorithm (Lu et al. 2006), which takes all the stellar features into account and can thus significantly mitigate the problem of overfitting. A power law is adopted to describe the AGN continuum. The optical Fe ii multiplets are modeled by two separate sets of templates constructed by Dong et al. (2008, 2011) based on the measurements of I Zw 1 by Véron-Cetty et al. (2004), one for broad lines and the other for narrow lines.

The next step is to fit the emission lines and to select broad-line candidates. We focus on the broad Hα line (HαB) since it is the strongest broad line in the optical spectra of AGNs. Our initial criteria for the addition of a broad component of Hα are as follows.10 (1) It would result in a significant decrease in χ2 with a chance probability of the F-test <0.05, (2) the width of HαB is relatively larger than those of narrow lines, particularly [O iii] λ5007, (3) HαB has a statistically significant flux, namely, Flux(HαB) > 3σ, where σ is the statistical noise, and (4) Flux(HαB) > 10−16 erg s−1 cm−2. After removing the continua, the residual spectra are first fitted in order to remove objects without any broad lines according to our selection criteria, resulting in a significant reduction of the number of objects needing refined fitting (∼303,000 objects are removed). Next, the Hα region is initially modeled using pure narrow-line profiles without a broad component. The narrow Hα + [N ii] λλ6548,6583 lines are fitted with a narrow-line model built up from the [S ii] λλ6716,6731 doublets or from the core of [O iii] λλ4959,5007 if [S ii] is weak. The profiles and redshifts of Hα and the [N ii] doublets are assumed to be the same as those of the narrow-line model obtained above. The centroid wavelengths of these narrow lines, as well as the flux ratios of the [N ii] doublets λ6583/λ6548 and the [O iii] doublets λ5007/λ4959, are fixed to their theoretical values, respectively. Then, a possible additional broad component of Hα is considered if it satisfies the broad-line criteria. This step results in ∼17,500 spectra leftover in which a candidate broad Hα component may be present.

The broad Hα lines of LMBH AGNs are generally relatively narrow and/or weak, and their fluxes are susceptible to the subtraction of the narrow lines. Thus, for the remaining ∼17,500 spectra, we employ a series of schemes to fit the narrow Hα line, including narrow-line profiles built up from narrow Hβ, [S ii], [O iii], or [N ii]. The result with the minimum reduced χ2 is adopted as the best fit. An example spectrum and the best-fit model are shown in Figure 1.

Figure 1.

Figure 1.

Illustration of the continuum and emission-line fitting for one of our LMBH AGNs as an example. Panel (a): the observed SDSS spectrum (black), the total model (blue), the decomposed components of the host galaxy (red), the AGN continuum (green), and the Fe ii multiplets (purple). Panel (b): emission-line profile fitting in the Hβ + [O iii] region. Panel (c): emission-line profile fitting in the Hα + [N ii] + [S ii] region. (The illustration of the continuum and emission-line fitting for all the LMBH AGNs are available in the figure set.) (The complete figure set (204 images) is available.)

Standard image High-resolution image

To ensure the reliability of the existence of a broad component, we further apply a stricter cut on the signal-to-noise ratio of the broad Hα flux, S/N(HαB) ≥ 5, where S/N(HαB) is the ratio of the broad Hα flux to the total uncertainties (σtotal), S/N(HαB) = Flux(HαB)/σtotal. σtotal is the quadrature sum of statistical noise (σstat), the uncertainties arising from subtraction of the continuum (σcont_sub), and the noise caused by the subtraction of the narrow lines (σNL_sub), namely, ${\sigma }_{\mathrm{total}}^{2}\,={\sigma }_{\mathrm{stat}}^{2}+{\sigma }_{\mathrm{NL}\_\mathrm{sub}}^{2}+{\sigma }_{\mathrm{cont}\_\mathrm{sub}}^{2}$. σNL_sub is estimated using the remaining n sets of fitting results that are worse than the best one and with the chance probability of the F-test greater than 0.1,

Equation (1)

As for σcont_sub, the term is negligible since the Hα absorption features are weak in most cases, and the continuum decomposition for each broad-line AGN candidate has been visually checked to ensure that the subtraction error is well below 1σ. With the cut of S/N(HαB), we obtain 6092 objects with reliable broad Hα detections, including 1850 "galaxies" and 4242 "QSOs."

3. LMBH AGN Sample from SDSS DR7

The optical spectra of low-mass AGNs are often substantially contaminated by starlight, and it is difficult to directly measure the continuum luminosities. Thus, the luminosities of the broad Balmer lines are adopted to estimate the BH masses. In this study, for ease of comparison, the BH masses are calculated using the formalism presented in Greene & Ho (2007a; hereafter GH07) following Dong et al. (2012b),

Equation (2)

The formalism was derived using the empirical correlations of $\lambda {L}_{\lambda }(5100\,\mathring{\rm A} )\mbox{--}{L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ and ${\mathrm{FWHM}}_{{\rm{H}}{\alpha }^{{\rm{B}}}}\mbox{--}{\mathrm{FWHM}}_{{\rm{H}}{\beta }^{{\rm{B}}}}$ from Greene & Ho (2005),11 along with the radius–luminosity relation reported by Bentz et al. (2006). Finally, a sample of 204 objects with BH masses less than 2 × 106 M12 are obtained. Table 1 summarizes the basic data of the sample. The emission-line parameters obtained from the best-fit models as described in Section 2, are tabulated in Tables 2 and 3, along with the BH mass and luminosity. The flux of Fe ii λ4570 is calculated by integrating the flux density of the corresponding Fe ii multiplets from 4434 to 4684 Å in the rest frame.

Table 1.  The SDSS DR7 LMBH Sample

ID SDSS Name z g g − r Ag
(1) (2) (3) (4) (5) (6)
1 J054248.74+004019.2 0.0518 17.49 0.98 1.14
2 J074345.47+480813.5 0.0181 16.24 0.74 0.22
3 J080807.13+563832.4 0.0990 18.93 0.55 0.18

Note. Column (1): identification number assigned in this paper. Column (2): official SDSS name in J2000.0. Column (3): redshift measured by the SDSS pipeline. Column (4): Petrosian g magnitude, uncorrected for Galactic extinction. Column (5): Petrosian g − r color. Column (6): Galactic extinction in the g band.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Table 2.  Emission-line Measurements

ID [O ii] λ3727 Fe ii λ4570 HβN HβB [O iii] λ5007 [O i] λ6300 HαN HαB [N ii] λ6583 [S ii] λ6716 [S ii] λ6731 ${\mathrm{FWHM}}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ FWHM[O iii] FWHM[S ii]
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15)
1 <−15.46 <−14.05 −15.57 <−15.54 −14.98 −15.73 −14.82 −14.70 −15.02 −15.59 −15.64 1986 227 133
2 <−14.54 −15.12 −14.92 −14.27 −15.16 −14.43 −13.95 −14.91 −15.00 −15.06 947 222 166
3 −14.80 −15.08 −14.99 −14.80 −14.38 −15.59 −14.39 −14.21 −14.62 −15.08 −15.20 1043 204 136

Note. Column (1): identification number assigned in this paper. Columns (2)–(12): emission-line fluxes (or 3σ upper limits) in log-scale, in units of erg s−1 cm−2. The measured emission-line fluxes are regarded to be reliable detections if they have significance greater than 3σ, or else the 3σ values will be adopted as the upper limits. Note that these are observed values only corrected for Galactic extinction, and no NLR or BLR extinction correction has been applied. The superscripts "N" and "B" in columns 4, 5, 8, 9, and 13 refer to the narrow and broad components of the line, respectively. Columns (13)–(15): line widths (FWHM) that are calculated from the best-fit models and have been corrected for instrumental broadening using the values measured from arc spectra and tabulated by the SDSS, in units of km s−1.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Table 3.  Luminosity, Mass, and Eddington Ratio Measurements

ID Mg(total) Mg(AGN) Mg(host) log ${L}_{{\rm{H}}{\alpha }^{B}}$ log L[O iii] λ5007 log MBH log Lbol/LEdd log M*
(1) (2) (3) (4) (5) (6) (7) (8) (9)
1 −20.54 −15.78 −20.53 40.10 39.82 6.24 −1.61 10.65a
2 −18.49 −15.38 −18.43 39.92 39.60 5.48 −1.01 9.80a
3 −19.68 −18.12 −19.38 41.18 41.02 6.14 −0.58 10.04c

Notes. Column (1): identification number assigned in this paper. Column (2): total g-band absolute magnitude. Column (3): AGN g-band absolute magnitude, estimated from the ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ given in column (5) and a conversion from ${L}_{{\rm{H}}\alpha }$ to Mg assuming fλ ∝ λ−1.56. Column (4): host galaxy g-band absolute magnitude, obtained by subtracting the AGN luminosity from the total luminosity. Column (5): Luminosity of broad Hα, in units of erg s−1. Column (6): luminosity of [O iii] λ5007, in units of erg s−1. Column (7): virial mass estimate of the BH calculated following GH07 and Dong et al. (2012b), in units of M. Column (8): Eddington ratio. Column (9): stellar mass of the host galaxy.

aEstimated using the color–M*/L relation in Into & Portinari (2013), where color is the SDSS r − i, and L is the luminosity of 2MASS Ks band. bEstimated using the method similar to that in note a, but using UKIDSS K-band magnitude. cDerived from the MPA–JHU catalog. dEstimated using the scaling relation between the WISE 3.4 μm luminosity and stellar mass provided by Wen et al. (2013), which is calibrated using the stellar masses from the MPA–JHU catalog.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Our sample has a median redshift of 0.1 (see Figure 2) and BH masses ranging from 1.1 × 105 M to 2.0 × 106 M. The FWHMs of HαB span a range of ∼500–2200 km s−1, with a median of 1000 km s−1. The broad Ha line luminosities ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ are in the range of ∼1039.3–1042.1 erg s−1, with a median lying at 1039.3 erg s−1. We caution that dust extinction may attenuate the observed broad Hα flux for some of the low-mass AGNs and may lower somewhat the measured ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$. Given that the measurement of the broad Hβ line in low-mass AGNs suffers large uncertainties, it is difficult to quantify the dust reddening effect. Following GH07 and Dong et al. (2012b), we do not correct for the dust reddening in this work. The Eddington ratio (Lbol/LEdd) is defined as the ratio of the bolometric luminosity (Lbol) to the Eddington luminosity (LEdd = 1.26 × 1038 MBH/M). Lbol is derived from the optical luminosities at 5100 Å using a bolometric correction factor of 9.8 (McLure & Dunlop 2004), Lbol = 9.8λLλ(5100 Å), where λLλ(5100 Å) is derived from the broad Hα luminosity using the scaling relation given in Greene & Ho (2005). The Eddington ratios thus estimated for the current LMBH sample range from 0.01 to 2. Figures 3 and 4 show the distributions of BH masses and Eddington ratios, as well as the luminosities and FWHMs of broad Hα, for our sample. Those of the Dong et al. (2012b) and GH07 samples are also plotted for comparison.

Figure 2.

Figure 2. Redshift distribution for our 204 low-mass AGNs. The dashed line denotes the median.

Standard image High-resolution image
Figure 3.

Figure 3. Distributions of the BH mass, Eddington ratio, luminosity, and the FWHM of broad Hα for the LMBH sample in this work (black solid histograms), GH07 (red dotted histograms), and Dong et al. (2012b) (blue dotted–dashed histograms). The vertical lines denote the corresponding medians. Numbers in brackets in the upper-left panel indicate the sample sizes.

Standard image High-resolution image
Figure 4.

Figure 4. Distributions of the LMBH sample in this work (filled circles) and Dong et al. (2012b) (open circles) on the Lbol/LEdd vs. MBH plane (a) and on the ${\mathrm{FWHM}}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ vs. ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ plane (b).

Standard image High-resolution image

As expected, the distributions of the DR7 LMBHs are similar to those in Dong et al. (2012b). The median MBH of the present sample is 1.3 × 106 M, comparable to the 1.2 × 106 M of Dong et al. (2012b). As for Lbol/LEdd, the medians are −0.59 and −0.64 in log-scale, respectively. The medians of ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ and ${\mathrm{FWHM}}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ for the current sample and Dong et al. (2012b) are very close (41.00 versus 40.99 and 2.99 versus 3.02, calculated in logarithm, respectively). In addition, the standard deviations of these quantities, including MBH, Lbol/LEdd, ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$, and ${\mathrm{FWHM}}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$, in the two samples are consistent within 0.07 dex.

4. Sample Properties

In combination with the 309 objects in Dong et al. (2012b), we expand the SDSS LMBH AGN sample to a total of 513 sources up to DR7. This is the largest sample of low-mass AGNs so far, which have uniform and homogeneous selection criteria and well-measured AGN parameters, thanks to the homogeneity and accurate spectrophotometry of the SDSS. It is worthwhile to investigate the statistical properties for the total sample. In this section, we present some of the ensemble properties of the LMBH AGNs and their host galaxies based on the SDSS data, as well as data from X-ray and radio surveys.

4.1. Narrow-line Diagnostic Diagrams

Compared with H ii galaxies, AGNs can emit a harder continuum which results in a distinct ionization condition in surrounding gas. Specific emission-line ratios can help distinguish the central radiating sources that ionize the circumnuclear medium. In practice, two-dimensional diagrams of certain narrow-line ratios have been widely applied to discriminate between H ii galaxies and type 2 AGNs (e.g., Baldwin et al. 1981; Veilleux & Osterbrock 1987; Ho et al. 1997a; Kewley et al. 2001, 2006; Kauffmann et al. 2003a). The so-called BPT diagrams involving the narrow-line ratios of Hα, Hβ, [O iii], [N ii], [S ii], and [O i] are shown in Figure 5. In general, the distributions of our LMBHs are consistent with those of Dong et al. (2012b). On the [O iii] λ5007/Hβ versus [N ii] λ6583/Hα diagram (panel (a)), the vast majority of these LMBHs are located in the region of either Seyfert galaxies or composite objects, and only a few sources (∼10) fall into the pure star-forming region. On the [O iiiλ5007/Hβ versus [S iiλλ6716,6731/Hα and [O iλ6300/Hα diagrams (panels (b) and (c)), about two-thirds of these objects are found in the AGN region, while the rest are located in the H ii region.

Figure 5.

Figure 5. Narrow-line diagnostic diagrams of [O iii] λ 5007/Hβ vs. [N ii] λ 6583/Hα (a), vs. [S ii] λλ6716,6731/Hα (b), and vs. [O i] λ 6300/Hα (c) for the LMBH sample in this study (black open circles) and the Dong et al. (2012b) sample (red filled circles). The extreme starburst classification line (blue dotted curve) from Kewley et al. (2001) and the Seyfert–LINER line (blue dotted–dashed line) obtained by Kewley et al. (2006) are adopted to separate H ii regions, AGNs, and low-ionization nuclear emission-line region sources (LINERs). In panel (a), the purple dashed line corresponds to the pure star formation line given by Kauffmann et al. (2003a), and the blue dotted horizontal line represents [O iii] λ5007/Hβ = 3, which is conventionally used to separate Seyfert galaxies and LINERs.

Standard image High-resolution image

These results confirm the AGN nature for the vast majority of our sample, though a few objects have narrow-line spectra similar to those of star-forming galaxies. Regarding those objects located in the pure H ii portion, their broad Hα fluxes, FWHMs, and luminosities are all similar to those of the whole parent sample, with chance probabilities >0.1 according to the Kolmogorov–Smirnov test (KS test). Visual checks of their SDSS spectra and the statistics (broad Hα flux with S/N > 5) both demonstrate that the broad Hα components of these objects are significant and reliable. A probable explanation is that the spectra of these objects are strongly contaminated by starlight from those star formation regions in the host galaxies, given the relatively large aperture of the fiber of 3''.

We note that about two dozen objects fall into the region of LINERs (Heckman 1980) according to the Seyfert–LINER demarcation lines of either [O iiiλ5007/Hβ = 3 in panel (a) or Kewley et al. (2006) in panels (b) and (c). LINERs are commonly found in early-type galaxies with classical bulges containing supermassive BHs especially with old stellar populations (Ho et al. 1997b). LMBHs tend to reside in late-type disk-dominated galaxies often harboring pseudobulges just on the opposite side. Therefore, such LMBH LINERs are of great interest in the study of LINERs as a population, though it is not a surprise that only a small number of low-mass AGNs show LINER-like spectra. Mostly manifesting themselves as low-luminosity AGNs, LINERs are thought to have relatively low accretion rates, and their accretion flows are dominated by ADAF (advection-dominated accretion flow; see, e.g., Narayan & Yi 1995; Kewley et al. 2006; Ho 2009). The combination of low BH mass and low accretion rate makes it even difficult to identify LINERs with LMBHs in the presence of host galaxy light. These may explain the low incidence of LINERs in the LMBH sample.

4.2. X-Ray Properties

X-ray is an important bandpass to study BH accretion. A number of LMBH AGNs in the GH07 sample have been observed in X-ray with Chandra (e.g., Greene & Ho 2007c; Desroches et al. 2009; Dong et al. 2012a) which gave a snapshot of the X-ray properties of LMBHs. Focusing on low-mass AGNs with very low Eddington ratios (Lbol/LEdd ∼ 0.01), i.e., the faintest AGN population known, Yuan et al. (2014) studied the X-ray properties of four objects observed by Chandra in the Dong et al. (2012b) sample, and suggested there should exist a large population of underluminous LMBHs in the nearby universe. Plotkin et al. (2016) performed a similar analysis on seven LMBHs with low Eddington ratios from GH07 using Chandra observations. The two well-studied LMBHs, NGC 4395 and POX 52, were found to show rapid and strong X-ray variabilities (e.g., Iwasawa et al. 2000; Moran et al. 2005; Dewangan et al. 2008; Thornton et al. 2008). Their spectral characteristics resemble those of classical Seyfert 1 galaxies. Miniutti et al. (2009) performed a detailed analysis of four LMBHs from Greene & Ho (2004) using XMM-Newton observations, finding that they are extremely variable, and three out of the four objects show soft excess in their spectra. Similarly, Ludlam et al. (2015) followed up a study with 14 LMBHs from GH07 with XMM-Newton, eight of which show soft excess emissions. In addition, LMBHs also present diverse timing properties. In general, LMBHs tend to show a larger X-ray variability amplitude than their more massive cousins. Pan et al. (2015) and Ludlam et al. (2015) showed the previously known inverse correlation between MBH and the normalized excess variance of X-ray variability, which flattens at MBH ∼ 106 M and thus vanishes for LMBH AGNs.

In this study, we briefly investigate the statistical X-ray properties of the total sample of 513 sources based on the ROSAT All-Sky Survey (RASS; Boller et al. 2016) and pointed observations (2RXP). Detailed X-ray spectral and timing analyses using XMM-Newton are deferred to later work.

4.2.1. X-Ray Detection

Of the 513 sources, 85 were detected by RASS and 32 in pointed observations, with 15 detected in both (for these sources, the data from the pointed observations are adopted). Thus, a total of 102 objects have X-ray detections with the ROSAT Position-Sensitive Proportional Counter (PSPC). Their X-ray fluxes and spectral indices in the 0.1–2.4 keV band are estimated following Schartel et al. (1996) and Yuan et al. (1998, 2008), assuming an absorbed power-law spectral shape. As the first step, the photon index Γ is calculated from the two hardness ratios,13 if available. The absorption column density NH is fixed at the Galactic value (Dickey & Lockman 1990) or set to a free parameter if no meaningful Γ is obtained using the Galactic NH. For those sources without meaningful hardness ratios, the mean value (Γ = 2.36) is adopted. Next, we estimate the X-ray fluxes in the 0.1–2.4 keV band from the count rates using the energy-to-counts conversion factor (ECF; ROSAT AO-2 technical appendix, 1991), which is calculated from the ROSAT PSPC effective area for each source, for the given power-law spectrum with Γ and NH obtained above. The X-ray fluxes in the band of 0.5–2.0 keV are tabulated in Table 4, which range from 1.8 × 10−14 to 2.9 × 10−12 erg s−1 cm−2. The corresponding luminosities span a range from 8.0 × 1039 to 1.1 × 1044 erg s−1, comparable to 9.3 × 1039–6.9 × 1043 erg s−1 for the GH07 sample. This indicates that the X-ray emission is mostly dominated by AGN radiation, since the X-ray luminosities of normal galaxies are generally below this level.

Table 4.  ROSAT Detections

ID Count Rate NH Γ f0.5–2.0 keV log L0.5–2.0 keV L2500 Å αOX
(1) (2) (3) (4) (5) (6) (7) (8)
Our new sample
6 0.074 ± 0.014 4.40 ${2.71}_{-0.23}^{+0.22}$ −12.24 ± 0.08 42.31 ± 0.08 27.57 −1.27
27 0.058 ± 0.014 3.10 ${2.14}_{-0.31}^{+0.28}$ −12.35 ± 0.11 43.94 ± 0.11 28.57 −0.92
34 0.043 ± 0.014 3.50 −12.49 ± 0.14 42.76 ± 0.14 27.96 −1.19
37 0.116 ± 0.019 3.60 ${2.49}_{-0.22}^{+0.20}$ −12.08 ± 0.07 42.88 ± 0.07 27.92 −1.15
39 0.174 ± 0.022 1.70 ${2.89}_{-1.33}^{+1.42}$ −12.32 ± 0.06 42.34 ± 0.06 27.63 −1.30
45 0.031 ± 0.012 1.50 ${1.56}_{-0.67}^{+0.49}$ −12.63 ± 0.16 43.29 ± 0.16 28.16 −0.96
47 0.061 ± 0.013 2.60 ${2.25}_{-0.34}^{+0.33}$ −12.40 ± 0.10 42.63 ± 0.10 27.64 −1.11
59 0.049 ± 0.014 1.60 −12.67 ± 0.12 42.77 ± 0.12 27.78 −1.12
75 0.005 ± 0.001 2.10 ${2.31}_{-0.33}^{+0.29}$ −13.59 ± 0.05 40.86 ± 0.05 26.43 −1.33
82 0.040 ± 0.012 1.80 −12.72 ± 0.12 42.92 ± 0.12 28.00 −1.14
87 0.072 ± 0.015 1.70 −12.49 ± 0.09 43.68 ± 0.09 28.60 −1.06
96 0.014 ± 0.002 1.40 1.73 ± 0.16 −13.05 ± 0.05 41.84 ± 0.05 26.87 −1.06
99 0.058 ± 0.000 1.60 −12.60 ± 0.00 43.29 ± 0.00 28.08 −1.02
100 0.176 ± 0.023 2.30 ${1.71}_{-0.17}^{+0.15}$ −11.84 ± 0.06 43.27 ± 0.06 27.88 −0.89
103 0.072 ± 0.014 1.40 ${2.94}_{-0.28}^{+0.23}$ −12.79 ± 0.09 43.31 ± 0.09 28.13 −1.11
104 0.051 ± 0.013 1.20 −12.73 ± 0.11 43.45 ± 0.11 28.41 −1.08
110 0.007 ± 0.002 2.60 2.39 ± 0.35 −13.39 ± 0.10 42.63 ± 0.10 28.35 −1.38
117 0.060 ± 0.013 1.30 2.04 ± 0.16 −12.52 ± 0.09 41.71 ± 0.09 27.06 −1.22
118 0.022 ± 0.009 1.10 −13.12 ± 0.18 42.66 ± 0.18 27.84 −1.18
121 0.027 ± 0.011 1.10 −13.04 ± 0.17 41.93 ± 0.17 27.51 −1.34
123 0.698 ± 0.038 1.20 ${2.83}_{-0.76}^{+0.79}$ −11.81 ± 0.02 43.55 ± 0.02 28.37 −1.11
124 0.098 ± 0.016 1.10 ${1.57}_{-0.22}^{+0.20}$ −12.19 ± 0.07 43.48 ± 0.07 28.10 −0.87
125 0.052 ± 0.012 1.10 −12.74 ± 0.10 42.59 ± 0.10 28.19 −1.35
126 0.062 ± 0.013 1.20 ${1.26}_{-0.21}^{+0.20}$ −12.30 ± 0.09 42.79 ± 0.09 28.15 −1.13
130 0.297 ± 0.023 1.20 ${2.43}_{-1.01}^{+1.02}$ −12.00 ± 0.03 42.74 ± 0.03 27.69 −1.11
132 0.053 ± 0.012 0.90 −12.79 ± 0.10 43.36 ± 0.10 28.26 −1.06
135 0.052 ± 0.000 1.10 ${2.88}_{-0.75}^{+0.77}$ −12.99 ± 0.00 42.03 ± 0.00 28.01 −1.56
138 0.083 ± 0.023 2.00 −12.38 ± 0.12 42.83 ± 0.12 28.00 −1.18
144 0.029 ± 0.010 1.40 ${1.93}_{-0.27}^{+0.26}$ −12.78 ± 0.15 43.06 ± 0.15 28.05 −1.05
150 0.028 ± 0.004 1.60 1.78 ± 0.18 −12.73 ± 0.06 42.97 ± 0.06 28.05 −1.07
155 0.098 ± 0.015 1.00 −12.49 ± 0.07 43.25 ± 0.07 28.14 −1.06
165 0.036 ± 0.014 1.90 −12.75 ± 0.17 43.30 ± 0.17 28.34 −1.11
166 0.095 ± 0.005 2.40 2.05 ± 0.11 −12.18 ± 0.02 42.74 ± 0.02 27.35 −0.93
176 0.023 ± 0.010 3.80 −12.73 ± 0.18 42.45 ± 0.18 27.49 −1.13
177 0.116 ± 0.022 3.90 ${2.37}_{-0.33}^{+0.29}$ −12.03 ± 0.08 42.55 ± 0.08 27.32 −1.04
186 0.047 ± 0.012 5.50 −12.32 ± 0.11 42.38 ± 0.11 27.62 −1.21
190 0.015 ± 0.001 3.50 ${1.27}_{-0.23}^{+0.18}$ −12.78 ± 0.04 41.56 ± 0.04 26.35 −0.92
192 0.111 ± 0.015 1.70 ${1.65}_{-0.18}^{+0.17}$ −12.09 ± 0.06 43.46 ± 0.06 28.15 −0.91
Dong et al. (2012b)
10 0.050 ± 0.013 3.9 ${2.53}_{-0.29}^{+0.27}$ −12.42 ± 0.11 42.76 ± 0.11 27.86 −1.18
15 0.030 ± 0.014 5.3 ${2.56}_{-0.44}^{+0.37}$ −12.55 ± 0.20 43.31 ± 0.20 27.94 −0.99
18 0.024 ± 0.004 3.6 ${1.55}_{-0.19}^{+0.16}$ −12.61 ± 0.13 42.48 ± 0.13 27.15 −0.90
21 0.067 ± 0.019 7.2 ${2.42}_{-0.16}^{+0.17}$ −12.11 ± 0.19 42.57 ± 0.19 27.82 −1.23
24 0.026 ± 0.011 8.1 −12.49 ± 0.16 43.16 ± 0.16 28.30 −1.16
27 0.026 ± 0.010 6.2 −12.54 ± 0.08 42.21 ± 0.08 27.55 −1.25
67 0.145 ± 0.025 2.8 ${2.91}_{-0.32}^{+0.31}$ −12.20 ± 0.13 42.83 ± 0.13 27.43 −1.04
68 0.046 ± 0.014 3.0 ${2.39}_{-0.29}^{+0.17}$ −12.52 ± 0.11 43.19 ± 0.11 28.03 −1.05
74 0.061 ± 0.015 4.9 ${3.83}_{-0.28}^{+0.27}$ −12.55 ± 0.08 42.75 ± 0.08 27.84 −1.36
78 0.105 ± 0.019 4.3 ${3.06}_{-0.45}^{+0.63}$ −12.18 ± 0.12 43.60 ± 0.12 28.11 −1.01
79 0.053 ± 0.014 3.7 ${3.96}_{-0.27}^{+0.34}$ −12.86 ± 0.08 43.98 ± 0.08 28.76 −1.23
87 0.113 ± 0.022 1.7 ${3.24}_{-0.87}^{+1.27}$ −12.66 ± 0.16 42.26 ± 0.16 27.69 −1.40
98 0.034 ± 0.013 1.6 −12.83 ± 0.07 42.23 ± 0.07 27.95 −1.39
99 0.107 ± 0.016 0.9 −12.48 ± 0.13 44.03 ± 0.13 28.40 −0.84
101 0.072 ± 0.000 0.7 −12.60 ± 0.20 42.17 ± 0.20 26.95 −1.01
102 0.040 ± 0.012 1.3 −12.82 ± 0.14 42.68 ± 0.14 27.89 −1.19
103 0.017 ± 0.008 2.5 −13.01 ± 0.14 43.09 ± 0.14 28.24 −1.16
105 0.033 ± 0.011 3.5 −12.60 ± 0.06 42.76 ± 0.06 27.53 −1.03
110 0.019 ± 0.002 3.0 ${2.61}_{-0.30}^{+0.31}$ −12.95 ± 0.12 42.47 ± 0.12 27.18 −1.03
112 0.029 ± 0.009 1.3 −12.96 ± 0.10 39.90 ± 0.10 25.45 −1.34
115 0.127 ± 0.017 0.7 ${2.53}_{-0.14}^{+0.15}$ −12.55 ± 0.05 42.54 ± 0.05 27.01 −0.93
117 0.002 ± 0.001 4.3 2.28 ± 0.05 −13.75 ± 0.06 41.79 ± 0.06 28.04 −1.58
123 0.046 ± 0.013 3.0 ${2.57}_{-0.22}^{+0.27}$ −12.56 ± 0.15 43.96 ± 0.15 28.85 −1.07
125 0.067 ± 0.016 1.4 ${2.37}_{-0.36}^{+0.42}$ −12.57 ± 0.12 43.23 ± 0.12 28.41 −1.18
135 0.005 ± 0.002 0.8 ${1.51}_{-0.06}^{+0.05}$ −13.50 ± 0.08 41.57 ± 0.08 27.09 −1.22
138 0.099 ± 0.000 1.9 ${1.46}_{-0.28}^{+0.33}$ −12.52 ± 0.15 43.62 ± 0.15 27.76 −0.84
149 0.032 ± 0.003 1.1 1.21 ± 0.13 −13.12 ± 0.12 42.43 ± 0.12 27.71 −1.26
153 0.115 ± 0.000 2.4 ${1.96}_{-0.41}^{+0.30}$ −12.19 ± 0.09 43.03 ± 0.09 27.89 −1.07
163 0.176 ± 0.021 1.1 ${2.52}_{-0.44}^{+0.52}$ −12.29 ± 0.17 42.51 ± 0.17 27.48 −1.13
164 0.226 ± 0.029 2.0 ${2.10}_{-0.12}^{+0.11}$ −11.86 ± 0.03 42.27 ± 0.03 27.77 −1.28
169 0.041 ± 0.014 1.9 −12.70 ± 0.13 42.40 ± 0.13 27.56 −1.18
174 0.044 ± 0.012 1.2 −12.72 ± 0.06 40.95 ± 0.06 26.84 −1.44
175 0.101 ± 0.019 1.6 1.56 ± 0.16 −12.12 ± 0.05 43.01 ± 0.05 27.60 −0.87
183 0.030 ± 0.011 1.0 −13.01 ± 0.09 42.46 ± 0.09 27.37 −1.08
184 0.045 ± 0.012 1.4 −12.75 ± 0.13 43.50 ± 0.13 28.54 −1.11
196 0.101 ± 0.022 2.0 ${1.44}_{-0.32}^{+0.30}$ −12.06 ± 0.14 42.23 ± 0.14 27.33 −1.06
203 0.021 ± 0.008 1.2 ${2.06}_{-0.96}^{+0.77}$ −13.01 ± 0.09 42.03 ± 0.09 27.67 −1.33
212 0.396 ± 0.027 0.9 −11.91 ± 0.07 43.67 ± 0.07 27.98 −0.85
214 0.032 ± 0.010 1.7 −12.81 ± 0.13 42.52 ± 0.13 27.98 −1.28
217 0.125 ± 0.016 1.8 ${2.36}_{-0.33}^{+0.31}$ −12.23 ± 0.16 42.77 ± 0.16 27.69 −1.08
221 0.121 ± 0.015 2.0 ${2.24}_{-0.16}^{+0.17}$ −12.17 ± 0.06 43.30 ± 0.06 28.12 −1.03
229 0.041 ± 0.008 1.6 −12.75 ± 0.07 42.97 ± 0.07 27.88 −1.07
230 0.027 ± 0.008 1.8 −12.89 ± 0.05 42.81 ± 0.05 27.98 −1.17
232 0.043 ± 0.014 2.8 ${2.28}_{-0.28}^{+0.37}$ −12.54 ± 0.08 42.95 ± 0.08 27.99 −1.12
235 0.055 ± 0.011 1.4 −11.90 ± 0.09 43.62 ± 0.09 27.64 −0.73
237 0.098 ± 0.016 1.5 ${2.89}_{-0.24}^{+0.29}$ −12.61 ± 0.07 43.04 ± 0.07 28.18 −1.23
239 0.047 ± 0.014 2.6 ${2.17}_{-0.41}^{+0.37}$ −12.55 ± 0.12 41.72 ± 0.12 26.70 −1.12
244 0.052 ± 0.005 1.3 ${1.69}_{-0.06}^{+0.05}$ −12.58 ± 0.02 42.76 ± 0.02 27.96 −1.15
246 0.024 ± 0.009 1.0 1.63 ± 0.15 −12.83 ± 0.06 42.35 ± 0.06 27.45 −1.07
249 0.095 ± 0.013 1.4 ${2.20}_{-0.46}^{+0.52}$ −12.36 ± 0.08 43.10 ± 0.08 28.14 −1.11
252 0.073 ± 0.012 1.4 ${2.49}_{-0.18}^{+0.16}$ −12.59 ± 0.12 42.72 ± 0.12 27.95 −1.22
253 0.107 ± 0.013 1.7 ${1.96}_{-0.42}^{+0.37}$ −12.18 ± 0.10 42.95 ± 0.10 28.04 −1.10
256 0.103 ± 0.018 1.5 ${2.87}_{-0.38}^{+0.39}$ −12.58 ± 0.14 42.67 ± 0.14 27.88 −1.26
259 0.038 ± 0.008 1.3 ${2.61}_{-0.45}^{+0.46}$ −12.95 ± 0.18 42.93 ± 0.18 27.73 −1.06
263 0.185 ± 0.028 2.4 ${3.37}_{-0.45}^{+0.40}$ −12.34 ± 0.14 42.76 ± 0.14 27.77 −1.26
269 0.042 ± 0.011 6.6 3.36 ± 0.10 −12.44 ± 0.06 42.99 ± 0.06 28.02 −1.26
271 0.895 ± 0.035 1.3 2.55 ± 1.21 −11.55 ± 0.00 42.89 ± 0.00 27.57 −1.03
276 0.088 ± 0.013 1.3 −12.40 ± 0.04 42.03 ± 0.04 27.10 −1.12
284 0.017 ± 0.003 2.5 ${2.65}_{-1.06}^{+0.78}$ −13.09 ± 0.15 42.24 ± 0.15 27.77 −1.36
286 0.017 ± 0.005 2.8 ${1.51}_{-0.61}^{+0.52}$ −12.81 ± 0.19 42.84 ± 0.19 27.68 −0.95
287 0.022 ± 0.005 3.1 ${2.20}_{-0.52}^{+0.53}$ −12.79 ± 0.00 42.62 ± 0.00 27.73 −1.14
289 0.034 ± 0.011 5.0 ${3.34}_{-1.18}^{+1.23}$ −12.66 ± 0.03 42.51 ± 0.03 27.98 −1.43
297 0.021 ± 0.009 4.9 ${3.20}_{-0.73}^{+0.75}$ −12.85 ± 0.00 43.05 ± 0.00 28.27 −1.30
301 0.045 ± 0.014 2.9 ${2.20}_{-0.88}^{+0.89}$ −12.49 ± 0.04 43.76 ± 0.04 27.79 −0.70

Note. Column (1): identification number assigned in this paper and Dong et al. (2012b), respectively. Column (2): ROSAT count rate (count s−1). Column (3): Galactic column density ${N}_{{\rm{H}}}$ (1020 cm−2). Column (4): photon index of X-ray spectrum. Column (5): X-ray flux in the 0.5–2.0 keV in log-scale (erg s−1 cm−2). Column (6): X-ray luminosity in the 0.5–2.0 keV in log-scale (erg s−1). Column (7): monochromatic luminosity at 2500 Å (erg s−1 Hz−1). Column (8): optical-to-X-ray effective spectral index αOX.

A machine-readable version of the table is available.

Download table as:  DataTypeset images: 1 1

For those undetected by the ROSAT, upper limits on the X-ray fluxes are estimated using the method in Yuan et al. (2008). An upper limit of 12 source counts is adopted for each undetected object, and thus the count rate upper limit is calculated using the corresponding effective exposure time from the RASS exposure map. Then, the flux upper limits are estimated using the method described above assuming the Galactic NH and the mean Γ of 2.36.

4.2.2. X-Ray versus [O iii] λ5007 Luminosity

The soft X-ray emission is susceptible to absorption, while the [O iii] luminosity is suggested to be an isotropic indicator of the intrinsic AGN power since the [O iii] line originates from the narrow-line region and should be unaffected by obscuration. A strong correlation between the X-ray and [O iii] λ5007 luminosity has been found in unobscured AGNs (e.g., Panessa et al. 2006). This can help to discriminate whether our low-mass AGNs are heavily obscured.

The relation between the 2--10 keV luminosity and the [O iii] λ5007 luminosity for the X-ray-detected sources is shown in Figure 6; overplotted are the low-mass AGN sample with ROSAT detections from GH07, LMBHs observed by Chandra from Dong et al. (2012a), more massive AGNs in Jin et al. (2012), and low-mass AGNs with log (Lbol/LEdd) < −1.5 observed using Chandra in Yuan et al. (2014) and Plotkin et al. (2016). Luminosities in the 2–10 keV band are calculated by extrapolating the spectra of ROSAT to 10 keV assuming a photon index of 2.5. The ROSAT-detected sources in our total sample are roughly consistent with the correlation between the X-ray and [O iii] λ5007 luminosity derived by more massive AGNs (Panessa et al. 2006). It indicates that our low-mass AGNs with X-ray detections suffer little or no obscuration in X-rays.

Figure 6.

Figure 6. X-ray luminosity in 2–10 keV vs. [O iii] λ5007 luminosity for our total LMBH sample detected with ROSAT (blue filled star symbols), GH07 sample with ROSAT detections (purple open inverted triangles), LMBHs with Chandra detections in Dong et al. (2012a) (black filled circles), low-mass active galaxies with low Eddington ratios from Yuan et al. (2014, orange red open diamonds) and Plotkin et al. (2016, cyan open triangles), and more massive AGNs from Jin et al. (2012, red open squares). The black solid line represents the relation for Seyfert galaxies and QSOs given by Panessa et al. (2006). Arrows denote upper limits.

Standard image High-resolution image

4.2.3. Optical–X-Ray Spectral Index

The optical/X-ray effective spectral index, αOX ≡ −0.384 log(f2500 Å/f2 keV), where f2500 Å and f2 keV are the rest-frame flux densities at 2500 Å and 2 keV, respectively, is commonly used to describe the relative dominance of the optical and X-ray emission (Tananbaum et al. 1979).14 The statistical properties of αOX have been well-studied for the classical AGNs (e.g., Avni & Tananbaum 1986; Yuan et al. 1998; Vignali et al. 2003; Steffen et al. 2006). Recently, Dong et al. (2012a) studied the αOX properties of 49 low-mass AGNs in the GH07 sample with Chandra observations. In this work, we briefly investigate the αOX properties of LMBHs based on ROSAT observations.

The αOX indices are calculated as follows. f2 keV is derived from the ROSAT observations as described in Section 4.2.1, and ${f}_{2500\mathring{\rm A} }$ is calculated from ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ using the relation between ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ and λLλ(5100 Å) given by Greene & Ho (2005) assuming a spectral index of 1.56 (fλ ∝ λ−1.56; Vanden Berk et al. 2001). The αOX values are listed in Table 4, and their distribution is shown in Figure 7. The αOX distribution has a large scatter, ranging from −1.58 to −0.70, with a median of $\langle {\alpha }_{\mathrm{OX}}\rangle =-1.12$, which is systematically flatter than those of AGNs with more massive BHs (e.g., αOX ∼ −1.5; Yuan et al. 1998). The αOX distribution of our low-mass BHs is similar to that of the ROSAT-detected sources in GH07, but slightly flatter than that of LMBHs with Chandra observations in Dong et al. (2012a). These LMBHs are roughly in accordance with the relation between αOX and the monochromatic luminosity at 2500 Å defined by those more massive AGNs (Steffen et al. 2006), albeit with a large scatter (see Figure 8). In addition, Figure 8 shows that low-mass AGNs have systematically flatter αOX and a larger scatter compared to those of more massive BHs, which is consistent with Dong et al. (2012a).

Figure 7.

Figure 7. Distributions of αOX for objects detected with ROSAT in the total LMBH sample (black shaded histogram). Chandra-detected sources in Dong et al. (2012a) (red dotted shaded histogram) and upper limits for those undetected in ROSAT (blue dotted–dashed histogram) in our sample are also plotted for comparison. The vertical lines represent the corresponding medians.

Standard image High-resolution image
Figure 8.

Figure 8. Optical-to-X-ray spectral index αOX vs. the monochromatic luminosity at 2500 Å for the total LMBH sample detected with ROSAT (blue filled star symbols). The 2500 Å monochromatic luminosities are derived from ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ using the scaling relation between ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ and λLλ(5100 Å) assuming a spectral shape of 1.56 (fλ ∝ λ−1.56). LMBHs in Dong et al. (2012a) (black filled circles), Yuan et al. (2014, orange red open diamonds), and Plotkin et al. (2016, cyan open triangles), as well as more luminous Seyfert galaxies and QSOs in Jin et al. (2012, red open squares), Wu et al. (2012, purple filled inverted triangles), and Just et al. (2007, orange filled triangles), are also plotted for comparison. The solid and dashed lines represent the relation and 1σ scatter given by Steffen et al. (2006), respectively. Arrows represent upper limits.

Standard image High-resolution image

As can be seen from Figure 8, although the distribution of our LMBH AGNs are mostly confined within the 1σ scatter of the αOXL2500 Å relationship, there is a much larger scatter when the Chandra samples of Dong et al. (2012a) and Yuan et al. (2014) are taken into account, which contains some ;X-ray-weak objects detected by Chandra owing to its high sensitivity. To explore further the large scatter of αOX for LMBH AGNs, we plot in Figure 9 αOX versus the Eddington ratio and BH mass for our sample as well as for the other AGN samples as in Figure 8. No significant correlation is found between αOX and Lbol/LEdd over a wide range for the LMBHs spanning nearly two orders of magnitudes. This is consistent with previous results for LMBH AGNs (e.g., Greene & Ho 2007c; Dong et al. 2012a). Regarding the αOXMBH relationship, no correlation is found for LMBHs only given the narrow dynamic range of MBH; however, a weak correlation (Spearman's rank correlation coefficient is −0.25 with a null probability of ∼10−7) appears to be present when taking into account AGNs with more massive BHs from the other samples, albeit with a large dispersion for the LMBHs. A similar weak correlation was also suggested by Dong et al. (2012a). Thus, the wide αOX distribution for LMBH AGNs is not caused by the distributions of the Eddington ratio or BH mass.

Figure 9.

Figure 9. Dependence of the optical-to-X-ray spectral index αOX on Eddington ratio (a) and BH mass (b). The blue filled stars represent the total LMBH sample. Overplotted symbols are the low-mass AGNs with the Eddington ratios log (Lbol/LEdd) < −1 from Yuan et al. (2014, orange red open diamonds) and Plotkin et al. (2016, cyan open triangles), LMBHs observed by Chandra from Dong et al. (2012a) (black filled circles), and more massive AGNs from Wu et al. (2012, purple filled inverted triangles) and Jin et al. (2012, red open squares), respectively. Arrows denote upper limits on αOX.

Standard image High-resolution image

Three possibilities might account for the large scatter in the αOX distribution. The first is X-ray variability. It has been demonstrated that the amplitude of short-timescale X-ray variability is anticorrelated with the BH mass for AGNs with MBH > 106 M (e.g., Ponti et al. 2012; Kelly et al. 2013), below which the relation flattens (Pan et al. 2015; Ludlam et al. 2015). This means that LMBH AGNs show the largest amplitude of X-ray variability among all AGNs. Moreover, LMBH AGNs appear to show strong variations in the X-ray spectral shape, as manifested by their wide range of spectral indices (Γ = 1.2–4.0). A typical example is NGC 4395, which was found to exhibit rapid and strong X-ray variability (Iwasawa et al. 2000; Dewangan et al. 2008), and its X-ray spectral slope varied from Γ < 1.25 to Γ > 1.7 over about one year (Moran et al. 2005). Similar behavior was also found in POX 52 (Thornton et al. 2008).

Second, LMBHs may have a wide distribution in the intrinsic X-ray luminosities. As commonly believed, the X-ray emission is mainly produced in the hot corona by inverse-Compton scattering of optical/UV photons from the disk. The relative dominance of the optical/UV and X-ray emission is determined by the fraction of energy deposited into the hot corona out of the total viscous energy produced in the accretion disk. For instance, in models where the corona is heated by magnetic reconnection (Liu et al. 2003), the fraction is mainly determined by the magnetic field strength in the corona, which may vary from one object to another.

Third, X-ray absorption may play a role as well. Although no strong effect of X-ray obscuration is present for most of the objects as suggested by the L2–10 keVL[O iii] relation above, moderate absorption cannot be ruled out in some of the objects. It is probable that the environment of LMBHs may be more gas-rich compared to normal AGNs powered by more massive BHs, whose stronger radiation makes it easier to expel circumnuclear materials.

4.3. Radio Properties

We explore the radio properties for our total LMBH sample using radio data at 20 cm from the VLA FIRST survey15 (Becker et al. 1995). There are 26 objects detected, leading to a low detection fraction of 5%. The low incidence of radio activity is consistent with the result in GH07. The radio powers at 20 cm of these sources range from 2 × 1021 to 4 × 1023 W Hz−1, with a median of 4 × 1022 W Hz−1 (see Table 5). Figure 10 shows their radio power versus the [O iii] luminosity, which appears broadly consistent with that found in GH07.

Figure 10.

Figure 10. Relation of the radio power at 6 cm and the [O iii] λ5007 luminosity for our total LMBH sample (blue stars), low-mass black holes from GH07 (red squares), and more massive Seyfert galaxies from Ho & Peng (2001, black circles). The filled and open symbols represent radio-loud and radio-quiet sources, respectively, and arrows denote upper limits. The solid and dashed lines represent the relation derived from more massive radio-loud and radio-quiet AGNs, respectively (Ho & Peng 2001).

Standard image High-resolution image

Table 5.  FIRST Detections

ID ${S}_{20\mathrm{cm}}$ log P20 cm log R
(1) (2) (3) (4)
Our new sample
23 1.14 ± 0.16 22.80 ± 0.06 1.36
28 0.85 ± 0.15 22.15 ± 0.07 1.22
55 2.86 ± 0.19 22.32 ± 0.03 1.49
89 1.22 ± 0.16 22.89 ± 0.06 1.51
95 16.75 ± 0.86 23.59 ± 0.02 2.49
123 2.72 ± 0.19 22.76 ± 0.03 1.05
143 1.97 ± 0.17 22.20 ± 0.04 1.51
199 2.67 ± 0.18 22.27 ± 0.03 0.80
201 1.79 ± 0.16 23.07 ± 0.04 1.34
Dong et al. (2012b)
22 1.82 ± 0.14 22.32 ± 0.03 1.57
30 3.18 ± 0.27 23.12 ± 0.04 2.03
36 4.53 ± 0.26 23.32 ± 0.03 2.06
54 2.07 ± 0.18 22.75 ± 0.04 1.43
78 1.44 ± 0.17 22.88 ± 0.05 1.44
118 1.04 ± 0.14 21.87 ± 0.06 0.84
140 2.38 ± 0.19 22.05 ± 0.03 1.57
141 2.55 ± 0.18 21.83 ± 0.03 1.58
170 1.03 ± 0.16 21.95 ± 0.07 1.58
181 1.29 ± 0.16 22.32 ± 0.06 0.91
206 2.33 ± 0.18 21.39 ± 0.03 1.49
208 8.81 ± 0.46 23.20 ± 0.02 2.21
222 2.69 ± 0.20 21.58 ± 0.03 1.15
257 2.52 ± 0.19 22.56 ± 0.03 1.50
273 1.43 ± 0.16 22.53 ± 0.05 1.09
277 0.68 ± 0.14 23.02 ± 0.09 1.36
305 3.78 ± 0.23 22.79 ± 0.03 1.98

Note. Column (1): identification number assigned in this paper and Dong et al. (2012b), respectively. Column (2): flux density at 20 cm from FIRST (mJy). Uncertainties include the 5% systematic uncertainty recommended by White et al. (1997). Column (3): corresponding radio power (W Hz−1) at 20 cm. Column (4): radio loudness in log-scale. R $\equiv \,{f}_{6\mathrm{cm}}/{f}_{4400\mathring{\rm A} }$ assuming a spectral index of αr  = 0.46 (radio; Ulvestad & Ho 2001) and αo = 0.44 (optical; Vanden Berk et al. 2001), where ${f}_{\nu }\propto {\nu }^{-{\alpha }_{o}}$.

A machine-readable version of the table is available.

Download table as:  DataTypeset image

As a common practice, the radio-loudness parameter, defined as the radio to optical flux ratio (R ≡ f6 cm/f4400 Å), is used to separate radio-loud AGNs from radio-quiet ones, with an operational dividing value of 10. The radio loudness is calculated for the radio-detected sources as follows: f6 cm is derived assuming a radio spectral index of αr = 0.46 (${f}_{\nu }\propto {\nu }^{-{\alpha }_{r}};$ Ulvestad & Ho 2001) from the flux density at 20 cm, while f4400 Å is estimated in the same way as f2500 Å. Among the 26 sources with FIRST detections, 23 are radio-loud, corresponding to a radio-loud fraction of 4%. For those sources undetected in FIRST,16 upper limits on radio loudness are estimated assuming a radio flux density of 1 mJy at 20 cm, which is the detection limit of the FIRST survey (see Figure 11).

Figure 11.

Figure 11. Distribution of the radio loudness of the total LMBH sample objects. Blue shaded histogram represents sources detected in FIRST, while the black histogram represents upper limits derived assuming a flux density limit of 1 mJy at 20 cm for the FIRST survey. The conventional demarcation line between radio-loud and radio-quiet AGNs is shown by the dashed line. The dotted–dashed lines denote corresponding medians.

Standard image High-resolution image

The radio-loud fraction in the LMBH AGN sample is lower than that in classical AGNs (∼15%; Ivezić et al. 2002). This seems to be consistent with the finding that radio sources may prefer to reside in more massive galaxies (Best et al. 2005). For our LMBHs, their host galaxies are indeed at the lower stellar mass end of the parent sample from SDSS DR7 (see Section 5.3 for details). However, it should be noted that the above radio-loud fraction should only be considered as a lower limit, since the FIRST flux sensitivity is not deep enough to give stringent constraints on the radio loudness. Further deep and high spatial resolution radio surveys are needed to determine the true radio-loud fraction for LMBHs.

As a natural interpretation, the radio emission of an active galaxy is a mixture of contribution from nuclear AGNs and star formation in the host galaxy. This may be particularly true for the fluxes listed in the FIRST catalog we used here, which were obtained by fitting the sources usually as extended radio sources (cf., the "Deconv.MajAx" column in the catalog). However, it is difficult to seek clues from the radio morphology since most of these sources have very weak radio emissions (S20 cm < ∼5 mJy, close to the detection limit of 1 mJy), and their FIRST images are unresolved. Only one brighter object (S20 cm > 10 mJy), SDSS J122412.51+160012.1, probably shows marginally extended structure with integrated to peak flux ratio of 1.5 and deconvolved major axis = 5farcs45. On the other hand, radio luminosities at 6 GHz may provide some information on the radio origins. It is suggested that for AGNs at low redshifts, low-luminosity radio emission with L6 GHz < 1023 W Hz−1 seems to be powered by star formation, while that with L6 GHz ≥ 1023 W Hz−1 is more likely AGN-dominated (e.g., Kimball et al. 2011; Condon et al. 2013; Kellermann et al. 2016). The luminosity at 6 GHz is derived assuming a spectral slope of 0.46. All the objects except one (SDSS J122412.51+160012.1) have 6 GHz luminosities located in the range of 1021 W Hz−1 ≤ L6 GHz ≤ 1023 W Hz−1, and thus their radio fluxes may also be contributed by host stellar processes.

4.4. Comparison with NLS1s

Most low-mass AGNs tend to have broad-line widths (FWHMs) narrower than ∼2000 km s−1, which is the conventional criterion for narrow-line Seyfert 1 galaxies (NLS1s; Osterbrock & Pogge 1985). It is thus interesting to compare these two AGN subclasses. In fact, some objects in our total sample have also been classified as NLS1s in a comprehensive study of NLS1 AGNs from SDSS DR3 by Zhou et al. (2006). NLS1s are characterized by some peculiar properties, including strong Fe ii lines, weak [O iii]/Hβ line ratios, high Eddington ratios, low radio-loud fraction, significant soft X-ray excess (in some), strong X-ray variability, and steep X-ray spectra (see, e.g., Laor 2000; Komossa 2008 for reviews). Based simply on the BH mass estimation method from the line width and luminosity, a relatively narrow broad-line width can result from either a low BH mass, or/and a high Eddington ratio. With their peculiar multiwavelength properties, NLS1s are found to cluster at one end of the so-called EV1, and are suggested to be driven by high accretion rates (Boroson & Green 1992). Thus, NLS1s are generally suggested to harbor relatively less massive BHs radiating near their Eddington limits.

However, LMBHs are selected only from BH masses regardless of the Eddington ratios. They are simply the low-BH-mass counterparts of classical Seyfert galaxies, and are expected to exhibit diverse properties depending on the Eddington ratios, which span a wide dynamical range. Greene & Ho (2004) found that their low-mass AGNs span a larger range in both the Fe ii and the [O iii] strengths relative to Hβ than classical NLS1s. There is a wide distribution of Eddington ratios in our sample, some as low as two orders of magnitudes below their Eddington limits. LMBHs show diverse properties in X-ray, with a wide range of X-ray photon indices (e.g., Γ = 1.0–2.7 in Desroches et al. 2009; Γ = 1.5–2.8 in Dong et al. (2012b), and Γ = 1.2–4.0 in the soft X-ray band for our sample). Moreover, some LMBHs, especially those with low Eddington ratios, do not show a soft X-ray excess, as found in the spectrum of NGC 4395. On the other hand, some typical NLS1s, such as those with extremely narrow widths of the broad lines as studied by Ai et al. (2011), do have BH masses as low as 106 MBH or lower.

It is clear that NLS1s, when solely selected from the line width, are a heterogeneous class of AGNs, which may include ordinary Seyfert galaxies having simply LMBHs, regardless of the Eddington ratios. Classical NLS1s were selected by also considering the strong Fe ii and weak [O iii] emissions. These tend to select objects with high Eddington ratios given the observed eigenvector i correlations. We thus suggest that a simple criterion based solely on the line width may not be a good approach for the selection of NLS1s and other criteria have to be incorporated. Similar remarks can also be found in Ai et al. (2011).

One interesting feature is that both types have relatively low fractions of radio-loud objects. This might be related to the suggested observed statistical relations among radio loudness, BH mass, and Eddington ratio. However, the underlying physical mechanisms driving these relations are not known.

5. Host Galaxy Properties

In this section, we briefly discuss the properties of the host galaxies for the LMBHs and the possible co-evolution with BHs by making use of SDSS data.

5.1. Luminosities, Colors, and Morphologies

The host galaxy luminosities are calculated following Dong et al. (2012b) and GH07. First, the AGN luminosities are estimated from the broad Hα luminosities17 that are free from starlight contamination, and are then subtracted from the SDSS photometric Petrosian g-band magnitude. The host galaxy luminosities are then derived after K-corrections using the routine of Blanton & Roweis (2007). The distributions of the g-band absolute magnitudes of the AGN, host, and the total are shown in Figure 12, with those from Dong et al. (2012b) and GH07 plotted for comparison. The distributions of these luminosities for the DR7 LMBH sample are similar to those of Dong et al. (2012b). The current sample has a median AGN luminosity of Mg = −17.80 mag and a median host galaxy luminosity of Mg = −20.14 mag (these values are −17.70 and −20.22, respectively, in Dong et al. 2012b). The median host galaxy luminosity is slightly brighter than the characteristic luminosity of ${M}_{g}^{* }=-20.1\,\mathrm{mag}$ at z = 0.1 (Blanton et al. 2003). Only six objects have host galaxy luminosities falling into the dwarf galaxy of Mg > −18.0 mag. It indicates that although LMBHs usually indeed reside in small stellar systems, their hosts may not necessarily be dominantly dwarf galaxies. Note that the host galaxy luminosities are possibly overestimated in our AGN-host separation method since fiber loss may result in the underestimation of the nuclear luminosity, though the systematic overestimation is smaller than 0.3 mag according to GH07. Furthermore, the AGN luminosities estimated from the Hα luminosities assuming a power-law continuum may be subject to an uncertainty of about 0.1 mag (Dong et al. 2012b).

Figure 12.

Figure 12. Distributions of the absolute g-band magnitudes of the AGN (top), host galaxy (middle), and the total (AGN plus host galaxy) (bottom) for the sample in this work (black solid line), Dong et al. (2012b) (blue dotted–dashed line), and GH07 (red dotted line; see text for details of the estimation of the AGN and host galaxy luminosities.)

Standard image High-resolution image

In general, it is not practical to visually classify the galaxy morphology for most of our sample objects using the SDSS images due to their limited depth and spatial resolution. Instead, we try to gain information on the morphologies from the host galaxy colors. The u − g colors are found to have a mean value of 1.28 mag with a standard deviation of 0.38 mag, and the g − r colors have a mean of 0.61 mag and a standard deviation of 0.16 mag. They correspond to the typical colors of Sbc galaxies according to Fukugita et al. (1995).

5.2. Stellar Populations

Two stellar indices, the 4000 Å break strength (D4000) and Balmer absorption line index (HδA), are used to infer the stellar populations of the host galaxies. In general, D4000 is small for young stellar populations (D4000 < 1.5 for ages < 1 Gyr) and large for old, metal-rich galaxies. On the other hand, galaxies that experienced starburst activity that has ended ∼0.1–1 Gyr ago tend to show strong Hδ absorption lines (Kauffmann et al. 2003b). Hence, these two indices can be used to constrain the mean stellar ages of the host galaxies of LMBHs and to diagnose the star formation history over the past few gigayears.

The two indices are calculated from the decomposed stellar component as described in Section 2.1 based on the definition and calculation method in Kauffmann et al. (2003b). Note that only 157 objects with AGN contribution less than 75% at 4000 Å in the SDSS spectra are used. The resulting typical errors (1σ) of the two indices for the LMBH sample are 0.1 and 1.6 Å, respectively (see Section 4.2.3 in Dong et al. 2012b for the details of error estimation).

The distribution of 419 low-mass AGNs (262 objects from Dong et al. 2012b also included) on the D4000–HδA plane is shown in Figure 13, and the distribution of ∼492,000 inactive galaxies is plotted as contours for comparison. These inactive galaxies are derived from SDSS DR7 spectra classified as "galaxies" by the SDSS pipeline and with S/N > 5 pixel−1 around 4000 Å in their spectra. In addition, broad-line AGNs and those Seyfert 2 galaxies defined by the BPT diagram are excluded. Note that the negative values of HδA are due to the definition and corresponding calculation method of the line index (Worthey & Ottaviani 1997; Kauffmann et al. 2003b). Most of our LMBH hosts (112 out of 157) have D4000 < 1.5, indicating that their mean stellar ages tend to be less than 1 Gyr. In general, for those galaxies with D4000 < 1.5, if they have experienced an instantaneous burst of star formation in the last few gigayears, they tend to have HδA > 5 Å. However, most objects with D4000 < 1.5 in our LMBH sample do not follow this tendency. It indicates that the distribution of our LMBHs on the D4000–HδA plane is more likely overlapping the locus of continuous star formation, which is in broad agreement with that of normal galaxies.

Figure 13.

Figure 13. Distribution of 423 low-mass BH host galaxies (blue filled circles) from the total LMBH sample on the plane of the 4000 Å break strength (D4000) vs. equivalent width of the Hδ absorption (HδA). As a comparison, the distribution of ∼492,000 non-AGN galaxies in SDSS DR7 is also plotted as contours. The cross at the upper-right corner represents the typical size of the 1σ errors.

Standard image High-resolution image

Studies of optical images and spectra indicated that there are two distinct populations of host galaxies of the LMBHs (Greene et al. 2008; Jiang et al. 2011). The majority are gas-rich, late-type galaxies similar to the prototype NGC 4395. The rest are spheroidal galaxies like POX 52, which are gas-deficient and have red color. Kormendy & Bender (2012) suggested that such spheroidal galaxies were transformed from irregular and late-type disk galaxies at an early epoch by both internal and secular processes. The low fraction of LMBH host galaxies with nearby companions also imply that these galaxies mainly evolve through secular processes (Jiang et al. 2011). This is also supported by the above finding that most of the host galaxies of our LMBHs have undergone continuous star formation in the past few gigayears.

5.3. Stellar Masses

The stellar masses of the host galaxies are estimated using three methods. The first is by using near-infrared (NIR) photometry, which can be considered as a tracer of stellar mass. In the NIR band, the luminosities of the LMBH AGNs are dominated by host galaxies, and AGNs have negligible contribution (∼7%).18 Moreover, the mass-to-light ratio in the NIR band is sensitive neither to dust absorption nor to star formation history. Of the total sample, 180 objects are detected in the Two Micron All-Sky Survey (2MASS19 ) and among the rest, 129 objects are detected in the UKIRT Infrared Deep Sky Surveys (UKIDSS; Lawrence et al. 2007). Here we adopt the mass-to-light ratios provided by Into & Portinari (2013; the scatter is about 0.1 dex), which are dependent on the galaxy colors r − i,

Equation (3)

where LKs is the Ks-band luminosity from 2MASS in units of L, and r − i is the color derived from the SDSS Petrosian magnitudes after K-correction and subtraction of the AGN continuum. The K-correction is calculated following Blanton & Roweis (2007). The K-band magnitudes of the UKIDSS-detected sources are transformed to the 2MASS Ks band using the color equations given by Hewett et al. (2006). For the sources not detected in 2MASS and UKIDSS, we derive their stellar masses from the MPA–JHU catalog20 (Kauffmann et al. 2003a; Brinchmann et al. 2004; Tremonti et al. 2004; Salim et al. 2007) if available, or otherwise using the relation between M* and the 3.4 μm luminosity from the Wide-field Infrared Survey Explorer (WISE) All-Sky survey calibrated using the MPA–JHU catalog (Wen et al. 2013). We compare the stellar masses given by the MPA–JHU catalog with those derived from the NIR luminosities and r − i colors for those objects detected in both surveys, and found that the two estimators are statistically in good agreement. The stellar masses for the total sample range from 6.0 × 108 M to 2.5 × 1012 M, with a median of 3.7 × 1010 M, slightly lower than that of the MPA–JHU sample (5.7 × 1010 M; see Figure 14 for their distributions). The stellar masses are mostly greater than 109.5 M, which is the mass of the Large Magellanic Cloud (LMC), indicating again that LMBHs do not necessarily reside in dwarf galaxies.

Figure 14.

Figure 14. Distributions of the stellar masses (M*) of the host galaxies of our LMBH sample (black solid histogram) and ∼925,000 sources in SDSS DR7, the stellar masses of which are taken from the MPA–JHU catalog (purple dotted–dashed histograms). The vertical lines represent the corresponding medians. The black shaded histogram represents radio-loud AGNs in our LMBH sample. (The MPA–JHU histogram is normalized to have a peak value of 100 for ease of comparison.)

Standard image High-resolution image

The 512 sources with reliable stellar mass estimations are plotted on the color–mass diagram (see Figure 15). The green lines indicate the location of the so-called green valley defined by Schawinski et al. (2014), which is a superposition of red and blue galaxies with the same intermediate optical colors. It has been suggested to be a special stage in the evolution of massive galaxies after star formation is quenched, and perhaps as evidence of AGN feedback. Our LMBH sources span almost the entire u − r color range, while most are located in the green valley or blue sequence, which is consistent with the suggestion that most low-mass AGNs in GH07 reside in gas-rich, disk galaxies (Jiang et al. 2011).

Figure 15.

Figure 15. Distribution of the total LMBH sample on the the u − r color vs. stellar mass plane. The green dashed lines denote the green valley defined by Schawinski et al. (2014). The red square, black filled circle, and blue star represent sources with their Eddington ratios in the ranges of log Lbol/LEdd < −1, −1 ≤ log Lbol/LEdd < −0.5, and log Lbol/LEdd ≥ −0.5, respectively. Open circles denote radio-loud sources. Our LMBH sample objects span almost the entire u − r color range, and most of the objects are located in the green valley or blue sequence, while about 11% of the sample are in the red sequence. The cross at the upper-right corner represents the typical size of 1σ errors.

Standard image High-resolution image

We note that there still exist about 11% of our sample objects in the red sequence. An interesting question arises: how were these LMBH AGNs triggered and fueled in such red and presumably gas-deficient galaxies? We know that AGN activity requires both available feeding fuel (gas or stars) and a process to get out of their angular momentum. As discussed in Kormendy & Ho (2013), the feeding rates of low-mass AGNs are very modest. Assuming these LMBHs accrete at their Eddington limits, the mass accretion rate, $\dot{M}=2.2\times {10}^{-8}(\eta /0.1)({M}_{\mathrm{BH}}/{M}_{\odot }){M}_{\odot }$, is only ∼0.02 M yr−1 for BHs with MBH ∼ 106 M (Kormendy & Ho 2013). These values are tiny even for red galaxies, and hence there may always be enough gas to feed the nuclear low-mass AGNs. On the other hand, various physical process may break the "angular momentum barrier" and push the gas into the nuclear regions of galaxies, such as wet major and minor mergers of the galaxies (e.g., Silk & Rees 1998; Hopkins et al. 2008), inflow along spiral arms or bar (for instance, NGC 1097; e.g., Fathi et al. 2006; Davies et al. 2009), and disk instabilities (e.g., cold flows; Dekel et al. 2009; Bournaud et al. 2011). Current observations give evidence that the triggering of low-mass AGNs are indeed dominated by secular processes. For instance, about 90% of LMBH AGNs in the low-redshift universe reside in galaxies with so-called pseudobulges, which are formed mainly by slow process without major mergers involved (Jiang et al. 2011). Our result on the D4000–HδA distribution also support this idea. The other processes, however, may not be completely ruled out. For example, SDSS J083803.68+540642.0, a red and gas-poor LINER in the sample of Dong et al. 2012b, shows a peculiar circumgalactic ring, which may be caused by collision with a gas-rich galaxy, suggesting a history of violent galaxy interaction (Liu et al. 2017). Such a process may also play a role in triggering low-mass AGNs. As discussed in Liu et al. (2017), major mergers of low-mass galaxies in the local universe may be more common than massive galaxies. Further studies are needed to fully understand this question.

5.4. Co-evolution of LMBHs and Their Host Galaxies?

As discussed above, LMBHs seem to evolve with secular processes and have not experienced major mergers. Thus, they can be used to trace the initial stage of evolution of BHs and their host galaxies. We briefly investigate the co-evolution of BHs and host galaxies in the low-mass regime by using the relation between BH mass accretion rates and host star formation rates (SFRs). The mass accretion rates are derived from the bolometric luminosities assuming the efficiency factor η = 0.1, while the SFRs are given by the MPA–JHU catalog with no aperture corrections. Figure 16 shows the distribution of 282 low-mass AGNs with both quantities available on the diagram of mass accretion rates versus SFRs, revealing a strong correlation between them.

Figure 16.

Figure 16. Mass accretion rate vs. star formation rate for 282 sources in the total LMBH sample with AGN contribution lower than 75%. The mass accretion rates are derived assuming an efficiency of η = 0.1. The star formation rates are obtained from the MPA–JHU catalog estimated using a technique described in Brinchmann et al. (2004). Blue stars, red squares, and black circles represent sources with redshift z ≥ 0.1, 0.05 ≤ z < 0.1, and z < 0.05, respectively.

Standard image High-resolution image

This suggests that there may exist a connection between AGN activity and star formation. However, more compelling and direct evidence is still needed to support the co-evolution of LMBHs and their hosts. In fact, there is no strong evidence that LMBHs have direct feedback on their host galaxies in the local universe. Furthermore, the link between BH accretion and star formation may also be explained by the fact that they both depend on the gas from the same reservoir (Kormendy & Kennicutt 2004). Nevertheless, our result hints at a possible co-evolution scenario between LMBHs and their host galaxies.

6. Summary

Using the optical spectrometric data from SDSS DR7, we obtain a sample of 204 new AGNs with LMBHs, and expand the SDSS LMBH sample from DR4 (Dong et al. 2012b) to a total of 513 objects. This is the largest optically selected, broad-line low-mass AGN sample so far. The BH masses, estimated using the virial method, are in the range of 1.1 × 105 to 2.0 × 106 M, with a median of 1.3 × 106 M. The Eddington ratios range from 0.01 to 2, with a median of 0.26. The properties and distributions of the new sample are statistically consistent with those of Dong et al. (2012b).

We present some statistical properties of the combined LMBH AGN sample from this work and Dong et al. (2012b), focusing on the emission lines and multiwavelength properties including X-ray and radio, as well as their host galaxies. Most of the LMBHs are located in the Seyfert galaxy and composite regions on the narrow-line diagnostic diagrams, confirming their AGN nature. A total of 102 sources were detected by ROSAT. LMBHs with X-ray detections tend to follow the correlation between the X-ray luminosities (2–10 keV) and [O iii] λ5007 luminosities derived from more massive AGNs. The optical/X-ray effective indices αOX of X-ray-detected AGNs show a large scatter ranging from −1.58 to −0.70, and are systematically flatter than more massive AGNs. In general, they are broadly consistent with the extrapolation of the αOXL2500 Å relation to the low-luminosity end. No dependence of αOX is found on Lbol/LEdd, whereas a weak correlation between αOX and MBH is suggested, which is consistent with Dong et al. (2012a). Only 5% of the sources are detected in the FIRST survey, which are mostly radio-loud. Thus, we suggest that LMBHs are predominantly radio-quiet, though further deep radio observations are needed to confirm this result.

The host galaxies of LMBHs have g-band magnitudes (Mg) ranging from −22.2 to −15.9 mag, with a median comparable to the characteristic luminosity of ${M}_{{\text{}}g}^{* }=-20.1\,\mathrm{mag}$ at z = 0.1. The colors of the galaxies suggest mostly a type of typical Sbc in general. The galaxies have stellar masses (M*) ranging from 108.8 to 1012.4 M, with a median of 1010.6 M, which is slightly lower than that of the SDSS DR7 sample from the MPA–JHU catalog. Only a few tens have M* < 1010 M or Mg > −18 mag. Thus, low-mass BHs may live in lower-mass stellar systems, but do not necessarily reside in the dwarf galaxies. Most of the galaxies have mean stellar ages younger than 1 Gyr from their D4000 values. The locus on the D4000–HδA diagram indicates that they tend to have experienced continuous star formation over the past few gigayears, which is consistent with the suggestion that the host galaxies of LMBHs mainly evolve via secular processes (e.g., Jiang et al. 2011; Kormendy & Ho 2013). Their distribution on the color versus M* diagram shows that most of these are blue, late-type galaxies.

With homogeneous selection and accurate measurements of the spectral parameters, our SDSS LMBH sample provides a useful database to further explore the properties of low-mass BHs and their host galaxies, as well as to study the BH mass function in the low-mass regime.

This work is supported by the National Natural Science Foundation of China (grant No. 11473035 and No. 11473062) and the National Program on Key Research and Development Project (grant No. 2016YFA0400804). W.L. acknowledges support from the Natural Science Foundation of China grant (NSFC 11703079) and the "Light of West China" Program of Chinese Academy of Sciences (CAS). H.L. thanks NAOC for providing the computing resources on the Zen cluster. We are grateful to the anonymous referee for constructive comments that improved the paper. This work is mainly based on the observations obtained by the SDSS. We acknowledge the entire SDSS team for providing the data that made this work possible. We have made use of the ROSAT Data Archive of the Max-Planck-Institut für extraterrestrische Physik (MPE) at Garching, Germany. This publication also makes use of data products from the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All-Sky Survey (2MASS). WISE is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, which are funded by the National Aeronautics and Space Administration (NASA). 2MASS is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the NASA and the National Science Foundation (NSF). We also use the data from the UKIDSS, which uses the UKIRT Wide Field Camera.

Footnotes

  • The Laser Interferometer Space Antenna (LISA) is a space-based gravitational-wave observatory led by the European Space Agency (ESA), designed to detect and accurately measure gravitational waves. The new LISA mission (based on the 2017 L3 competition) is a collaboration of ESA and NASA.

  • According to Dong et al. (2012b), the broad Hα widths (FWHMs) of LMBHs range from 500 to 2200 km s−1, with a median of 1000 km s−1, which is much lower than that of the entire parent broad-line AGN sample (≈3000 km s−1) and even slightly lower than the traditional demarcation value between AGN broad and narrow lines (1200 km s−1; cf. Hao et al. 2005).

  • Note that there exists a population of spectra dominated by AGNs, the starlight components of which are negligible, and the broad lines are broad and significant. For this case we fit simultaneously the AGN power-law continuum together with the emission lines, including the Fe ii multiplets, and the forbidden and the Balmer lines (details see Dong et al. 2008).

  • 10 

    For details, see Section 3 of Dong et al. (2012b).

  • 11 

    In fact, Greene & Ho (2005) give the correlation between λLλ(5100 Å) and the combined Hα luminosity of both broad and narrow components. However, as described in Greene & Ho (2005), "the best-fit parameters are virtually unchanged when only the broad component is considered." In addition, the narrow Hα component may be contaminated by the host galaxies for our LMBHs. Thus, in this study, λLλ(5100 Å) is estimated using merely broad Hα component.

  • 12 

    We adopt the same upper limit of LMBHs as in GH07 and Dong et al. (2012b) for consistency.

  • 13 

    The two hardness ratios are defined as HR1 = (B − A)/(B + A), HR2 = (DC)/(D + C), where A, B, C, and D are the number of source counts in the energy channels of 11–41, 52–201, 52–90, and 100–201, respectively.

  • 14 

    The definition is different from the original one in Tananbaum et al. (1979) by a negative sign.

  • 15 

    VLA FIRST aims to produce the Faint Image of the Radio Sky at Twenty centimeters using the Very Large Array, which is operated by the National Radio Astronomy Observatory.

  • 16 

    Note that ∼30 sources are not covered by the FIRST, thus their upper limits are not calculated.

  • 17 

    Assuming an optical spectral shape of 1.56 (fλ ∝ λ−1.56; Vanden Berk et al. 2001), we obtain the AGN luminosities using the scaling relation between ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ and λLλ(5100 Å) described in Greene & Ho (2005).

  • 18 

    We simply estimate the AGN contribution in the K band using the relative flux ratio between 5100 Å and 2.2 μm from the mean nuclear spectral energy distribution (SED) of the nearest Seyfert 1 galaxies (Prieto et al. 2010) that includes the AGN continuum and torus emission in the NIR band. The flux at 5100 Å is derived from ${L}_{{\rm{H}}{\alpha }^{{\rm{B}}}}$ (Greene & Ho 2005).

  • 19 

    Our LMBHs tend to be extended sources since they are mostly located in the low-redshift universe. Thus, the result is derived by cross-matching the LMBH sample with the 2MASS All-Sky Extended Source Catalog (XSC).

  • 20 

    The MPA–JHU catalog, which contains the properties of millions of galaxies from SDSS DR7, was produced by a collaboration of researchers from the Max Planck Institute for Astrophysics (MPA) and the Johns Hopkins University (JHU). It contains two subsets, the raw data and the derived data. The raw data contain spectral parameters, including the line fluxes, equivalent widths, and continuum indices, as well as basic information about the objects such as their redshifts, velocity dispersions and plate, fiber, and MJD. The derived data include gas-phase metallicities, star formation rates, and stellar masses. The stellar masses are calculated from fits to the photometry with population synthesis models following the philosophy of Kauffmann et al. (2003a) and Salim et al. (2007). In this study, we use the improved version of the stellar mass catalog, which can be obtained from http://home.strw.leidenuniv.nl/~jarle/SDSS/.

Please wait… references are loading.
10.3847/1538-4365/aab88e