A systematic search for hidden type 1 AGNs: gas kinematics and scaling relations

TL;DR

This study identifies hidden type 1 AGNs within emission-line galaxies initially classified as type 2, revealing their properties, gas kinematics, and scaling relations, and highlighting differences in ionized gas outflows compared to type 2 AGNs.

Contribution

It presents a systematic method to uncover low-luminosity hidden type 1 AGNs and analyzes their gas kinematics and scaling relations, expanding understanding of AGN diversity.

Findings

Hidden type 1 AGNs have low luminosity and Eddington ratios.

Ionized gas in type 1 AGNs shows higher velocity dispersion and shifts.

Type 1 AGNs follow the black hole mass - stellar velocity dispersion relation.

Abstract

We search type 1 AGNs among emission-line galaxies, that are typically classified as type 2 AGNs based on emission line flux ratios if a broad component in the H$\alpha$ line profile is not properly investigated. Using ~24,000 type 2 AGNs at z $<$0.1 initially selected from Sloan Digital Sky Survey Data Release 7 by Bae, et al. 2014, we identify a sample of 611 type 1 AGNs based on the spectral fitting results and visual inspection. These hidden type 1 AGNs have relatively low luminosity with a mean broad \Ha\ luminosity, log L$_{\rm H\alpha }$ $=$ 40.73$\pm$0.32 \ergs\, and low Eddington ratio with a mean log L$_{bol}$/L$_{\rm Edd}$ $=$ -2.04$\pm$0.34, while they do follow the black hole mass - stellar velocity dispersion relation defined by the inactive galaxies and the reverberation-mapped type 1 AGNs. We investigate ionized gas outflows based on the [OIII] $\lambda$5007 kinematics, which show relatively high velocity dispersion and velocity shift, indicating that the line-of-sight velocity and velocity dispersion of the ionized gas in type 1 AGNs is on average larger than that of type 2 AGNs.

Equations

$${\lambda }_0=\frac{\int \lambda f_{\lambda }d\lambda }{\int f_{\lambda }d\lambda }$$

$${\sigma }^2=\frac{\int {\lambda }^2{f}_{\lambda }d\lambda }{\int f_{\lambda }d\lambda }-{\lambda }_0^2$$

$$M_{\mathrm{BH}}=f\times 10^{6.819}\times {\left(\frac{{\sigma }_{\mathrm{H}\beta }}{10^3\mathrm{km}{\mathrm{s}}^{-1}}\right)}^2\times$$

$${\left(\frac{\lambda L_{5100}}{10^{44}\mathrm{erg}{\mathrm{s}}^{-1}}\right)}^{0.533}{M}_{\odot }$$

$$M_{BH}=f\times 10^{6.544}{\left(\frac{L_{\mathrm{H}\alpha }}{10^{42}\mathrm{erg}{\mathrm{s}}^{-1}}\right)}^{0.46}$$

$$\times {\left(\frac{\mathrm{FWH}{\mathrm{M}}_{\mathrm{H}\alpha }}{10^3\mathrm{km}{\mathrm{s}}^{-1}}\right)}^{2.06}{M}_{\odot }$$

$$M_{BH}=f\times 10^{6.561}{\left(\frac{L_{\mathrm{H}\alpha }}{10^{42}\mathrm{erg}{\mathrm{s}}^{-1}}\right)}^{0.46}$$

$$\times {\left(\frac{{\sigma }_{\mathrm{H}\alpha }}{10^3\mathrm{km}{\mathrm{s}}^{-1}}\right)}^{2.06}{M}_{\odot }$$

$$L_{\mathrm{bol}}=2.21\times 10^{44}(\frac{L_{\mathrm{H}\alpha }}{10^{42}\mathrm{erg}{\mathrm{s}}^{-1}}{)}^{0.86}\mathrm{erg}{\mathrm{s}}^{-1}.$$

$${\chi }^2=\sum _{i=1}^N\frac{{({\mu }_i-\alpha -\beta s_i)}^2}{{{\sigma }_{\mu ,i}}^2+{\beta }^2{{\sigma }_{s,i}}^2+{{ϵ}_0}^2}+$$

$$\sum _{j=1}^M\frac{{({\mu }_{\mathrm{VP},\mathrm{j}}+\log \mathrm{f}-\alpha -\beta {\mathrm{s}}_{\mathrm{j}})}^2}{{\sigma }_{{\mu }_{\mathrm{VP},\mathrm{j}}}^2+{\beta }^2{{\sigma }_{\mathrm{s},\mathrm{j}}}^2+{{ϵ}_0}^2}$$

Full text

1. Introduction

The spectral features of AGNs (active galactic nuclei) are observed in different ways since the nucleus can be obscured by an optically thick dust torus (Antonucci, 1993; Urry & Padovani, 1995). According to the simplest AGN unification model, the observed physical properties show differences due to the orientation effect, which depends on the angle between the line-of-sight and the axis of the dust torus (Antonucci, 1993). If we directly observe the inner part of an AGN, a blue AGN continuum from an accretion disk and broad emission lines originated from the broad-line region (BLR) are present in the optical spectral range. In contrast, if the dust torus obscures the inner part of an AGN, we observe narrow emission lines originated from the photoionized narrow-line region (NLR), without the direct detection of a blue AGN continuum and broad emission lines.

If the FWHM of the broad emission line, e.g., H$\beta$ and H$\alpha$, is larger than 1000 km s*-1*, AGNs are typically classified as type 1 AGNs (e.g., Seyfert, 1943), for which the kinematics of gas in the BLR can be used to trace the gravitational potential of the central supermassive black hole for determining black hole mass (e.g., Blanford & Mckee, 1982; Kaspi et al., 2000; Peterson et al., 2004; Bentz et al., 2013). In addition, the relativistic effect can also produce gravitational redshift causing velocity shift of the broad emission lines, as the BLR gas has orbital velocity corresponding to a few percent of the speed of light (e.g., Zheng & Sulentic, 1990; Kollatschny, 2003; Tremaine et al., 2014).

When AGN luminosity is relatively low, however, the emission from host galaxy dilutes AGN emission, consequently dominating in the observed spectra. Thus, either a blue AGN continuum or broad emission lines can be missed from detection if no careful analysis is performed. In particular, a weak broad H$\alpha$ emission line from the BLR can be easily undetected if the narrow H$\alpha$ emission originated from the NLR and/or star forming region is very strong. Consequently, these low luminosity AGNs may be mis-classified as type 2 AGNs or star forming galaxies instead of type 1 AGNs. Previous studies by Woo et al. (2014) and Oh et al. (2015), for example, showed that by decomposing the broad and narrow components in the H$\alpha$ line profile, a hidden population of type 1 AGNs can be identified among emission line galaxies (see also Greene & Ho, 2005). The hidden population of type 1 AGNs is important to constrain the AGN unification model and to characterize dust torus since the empirical number ratio between type 1 and type 2 AGNs is a key to understand the orientation effect (Willott et al., 2000; Simpson, 2005; Oh et al., 2015) although it is difficult to overcome the different bias and systematic difference in selecting type 1 and type 2 AGNs. Searching for intermediate mass black holes, other studies examined the H$\alpha$ line profile to investigate the presence of a broad H$\alpha$ component and identified type 1 AGNs (Reines et al., 2013; Baldassare et al., 2016).

Hidden type 1 AGNs are characterized with low luminosity and low Eddington ratio (e.g., Woo et al., 2014; Oh et al., 2015). Thus, the optical continuum mainly represents stellar population while narrow emission lines, e.g., [O iii] $\lambda$5007, are similar to those of type 2 AGNs. Compared to typical type 1 AGNs, these hidden type 1 AGNs have advantage in studying their host galaxies and investigating the coevolution of black holes and galaxies by constraining the black hole mass scaling relations (Ferrarese & Merritt, 2000; Gültekin et al., 2009; Kormendy & Ho, 2013). For example, while it is very difficult to measure stellar velocity dispersion or bulge mass of type 1 AGN host galaxies due to the high flux ratios between AGNs and stars (e.g., Woo et al., 2006; Park et al., 2015) it is easy to measure the properties of host galaxies if the AGN continuum is weak (Reines & Volonteri, 2015; Woo et al., 2015). Thus these AGNs can be utilized to investigate the black hole-galaxy scaling relations and their evolution (Barth et al., 2005; Woo et al., 2013).

In our pilot study (Woo et al., 2014), we searched hidden type 1 AGNs among emission line galaxies at z $<$0.05, by carefully investigating the presence of a broad component in the H$\alpha$ emission line profile, using a large sample of $\sim$24,000 type 2 AGNs based on Sloan Digital Sky Survey (SDSS) Data Release (DR) 7 (Bae & Woo, 2014). In this paper we extend our previous work by enlarging the survey volume out to z=0.1, and improving the detection scheme for the broad component in H$\alpha$. Based on the newly identified sample of type 1 AGNs with a broad H$\alpha$ component, we investigate the MBH-$\sigma_{*}$ relation, the kinematics of ionized gas outflows, and the effect of gravitational redshift. We present sample selection and analysis methods in Section 2 and 3, respectively. In Section 4, we present the results on the MBH-$\sigma_{*}$ relation. In Section 5, we investigate the kinematics of [O iii] and the gravitational redshift of the broad H$\alpha$. Discussion is given in Section 6 and summary and conclusions follow in Section 7. Throughout the paper, we use the cosmological parameters as $H_{0}=70$ km s*-1* Mpc*-1*, $\Omega_{\rm m}=0.30$, and $\Omega_{\Lambda}=0.70$.

2. Sample Selection

We utilized the type 2 AGN catalog by Bae & Woo (2014) (see also Woo et al., 2016, 2017), which identified type 2 AGNs at 0.02$<$z$<$0.1, in order to search hidden type 1 AGNs. The sample is selected based on a couple of criteria: signal-to-noise ratio (S/N)$>$3 for the four emission lines, namely, H$\beta$, [O iii] $\lambda 5007$, H$\alpha$, and [N ii] $\lambda 6583$ from the MPA-JHU catalog for SDSS DR7 galaxies (Abazajian et al., 2009). Using the flux ratios of these four emission lines (Baldwin et al., 1981) and the demarcation line for separating AGNs from star-forming galaxies (Kauffmann et al., 2003), 23,517 type 2 AGNs (pure AGNs and composite objects)were identified with the amplitude-to-noise ratio >5 for [O iii] and H$\alpha$. More details of the selection procedure can be found in Bae & Woo (2014) and Woo et al. (2016).

To identify hidden type 1 AGNs among type 2 AGNs, we first subtract stellar continuum by fitting SDSS spectra using the penalized pixel-fitting code (pPXF) (Cappellari & Emsellem., 2004), which find the best-fit stellar population model of the given galaxy spectrum. We use MILES simple stellar population models with solar metallicity (Sánchez-Blázquez et al., 2006). In addition to the optical emission lines, we mask the continuum around the H$\alpha$ region (6300 $\sim$ 6900 Å) to prevent the stellar population model from fitting a potential broad component as a stellar continuum.

After subtracting the best-fit stellar population model, we investigate whether a broad H$\alpha$ component is present around the H$\alpha$ region by visually inspecting the residual spectra (i.e., emission-line spectra). In this process we initially checked that $\sim$1000 objects show a broad feature around the narrow H$\alpha$ and [N ii] lines. To confirm the presence of a broad H$\alpha$ component, we need a reliable emission line decomposition process since the H$\alpha$ region is complex due to the blending of the narrow H$\alpha$ and the [N ii] doublet as well as the broad H$\alpha$ component if present. Therefore, we use three different fitting approaches for the H$\alpha$ region as follows: (a) If the narrow emission lines do not show particular wing components and if there is no broad H$\alpha$ component, we use a single Gaussian model for each of the H$\alpha$ and two [N ii] lines. In case of the [N ii] doublet, we fix their centers and flux ratio at their theoretical value while we use the same line width; (b) If wing components are present in the narrow emission lines (H$\alpha$ and [N ii] doublet), we utilize a double Gaussian model (one for a core component and the other for a wing component) to account for outflow kinematics in the NLR as we used for typical type 2 AGNs (Woo et al., 2016) (c) If a very broad H$\alpha$ component with FWHM > 1000 km s*-1* is present, we add a single Gaussian component or a Gauss-Hermite component depending on whether this component is symmetric or asymmetric. For the broad and narrow H$\alpha$ components, we use free parameters for the line centers, widths and amplitudes.

For the total sample, we first apply model (a) and identify type 1 AGN candidates by visually inspecting whether the model-subtracted spectra show a significantly high residual around the H$\alpha$ region. Note that the residual can be present due to a very broad H$\alpha$ component or wing components of the narrow H$\alpha$ and [N ii] lines. For these type 1 AGN candidates, we then apply either model (b) or (c) depending on the line profile of the [S ii] doublet. If the [S ii] doublet shows wing components based on a single Gaussian fit and visual inspection, then we apply model (b) by including a wing component for each of H$\alpha$ and [N ii], while we use model (c) only if both model (a) and (b) do not provide an acceptable fit (i.e., a significant residual is present based on visual inspection). In this process, we conservatively identify type 1 AGNs only for the case with model (c) in order to avoid false detection of the broad H$\alpha$ component. In total, we obtained 611 type 1 AGNs with a broad H$\alpha$ lines with FWHM > 1000 km s*-1*. Note that we use FWHM $=$ 1000 km s*-1as a lower limit of the broad H$\alpha$ line width as a conventional definiton of type 1 AGNs (Vanden Berk et al., 2006) since a H$\alpha$ line narrower than 1000km s-1* may not be originated from the BLR. In Figure 1, we present examples with the best-fit models in the H$\alpha$ region.

3. Analysis

3.1. Measurement of gas kinematics

To investigate gas kinematics, we measure the velocity shift with respect to the systemic velocity based on stellar lines, velocity dispersion and luminosity of the [O iii] $\lambda$5007. We fit [O iii] using single or double Gaussian models to investigate the ionized gas outflows. First, we fit the total sample with a single Gaussian model. If the single Gaussian fit is not acceptable and a wing component is present based on visual inspection, we use a double Gaussian model. Since in some cases a noise is fit with a second Gaussian component, we accept the double Gaussian fit only if the peak of the wing component is a factor of three larger than the noise level as we adopted in Woo et al. (2016)(see Figure 1). Based on the best-fit model, we calculate the first and second moment of the total line profile as:

[TABLE]

[TABLE]

where $f_{\lambda}$ is the flux density and $\sigma$ refers to the second moment line dispersion. For the given $\lambda_{0}$ , the velocity shift is calculated with respect to the systemic velocity. Note that we follow the same procedure to analyze the [O iii] kinematics presented in our previous study of type 2 AGNs (Woo et al., 2016).

3.2. Morphology classification

We use the Galaxy Zoo human classification data to obtain the host galaxy morphology (Lintott et al., 2008, 2011). The host galaxies of the hidden type 1 AGNs are classified as 221 spirals ($\sim$36%) and 381 ellipticals ($\sim$62%). Nine ($\sim$1.5%) of them are not classified due to the low image quality. Spiral host galaxies are further divided into face-on and edge-on galaxies based on the minor-to-major axis ratio (i.e., b$/$a) obtained from the SDSS DR7. If the ratio is larger than 0.5, we classify them as edge-on galaxies. Otherwise, we classify them as face-on galaxies. According to this criteria, 221 spirals are divided into 44 edge-on spirals and 177 face-on spirals.

3.3. Black Hole Mass and Eddington ratio

To calculate black hole mass, we use the single-epoch virial mass estimator calibrated by Woo et al. (2015) as:

[TABLE]

where $\sigma_{\rm H\beta}$ is the line dispersion of the H$\beta$ line and L5100 represents AGN continuum luminosity at 5100Å. However, since we do not detect AGN continuum and barely see the H$\beta$ line due to their low flux in the large fraction of our sample, we use the width (i.e., FWHM and line dispersion) and the luminosity of the broad H$\alpha$ line as a proxy of H$\beta$. When we use the FWHM of H$\alpha$ broad component, we calculate the H$\alpha$-based MBH following Woo et al. (2015) as:

[TABLE]

where $L_{\rm H\alpha}$ is the luminosity of the broad H$\alpha$ and the FWHMHα is the width of the broad H$\alpha$ component. When we instead use the line dispersion of the broad H$\alpha$ ($\sigma_{H\alpha}$), we calculate MBH following Woo et al. (2015) as:

[TABLE]

We utilize log f=0.05 for the FWHM-based MBH estimator and log f=0.65 for $\sigma$-based MBH, which were obtained from the recent calibration based on the MBH-$\sigma_{*}$ relation of quiescent galaxies and H$\beta$ reverberation-mapped AGNs (Woo et al., 2015). We use the luminosity of the broad H$\alpha$ component ($L_{\rm H\alpha}$) as a proxy for the bolometric luminosity ($L_{\rm bol}$), by utilizing the relation between H$\alpha$ and the continuum luminosity at 5100Å (Greene & Ho, 2005), and the bolometric correction 9.26 for $L_{5100}$ (Richards et al, 2006) as:

[TABLE]

4. Mass distribution and the scaling relations

4.1. Comparison with typical type 1 AGNs

In this section, we investigate the properties of the newly identified type 1 AGNs. The most prominent feature of these type 1 AGNs is that they have very weak AGN continuum, hence, the observed continuum represents the stellar continuum. To compare the hidden type 1 AGNs with typical type 1 AGNs, we choose 882 type 1 AGNs at the matched redshift range z$<$0.1 from the catalogue by Oh et al. (2015). Using the luminosity and FWHM of the broad H$\alpha$ emission lines of the typical type 1 AGNs, we obtain bolometric luminosity, black hole mass and the Eddington ratio in the same way as we described in Section 3.3. Note that we use FWHMHα based MBH in this section.

The hidden type 1 AGNs have relatively low broad H$\alpha$ luminosity while they have similar FWHM of the broad H$\alpha$ compared with typical type 1 AGNs (Figure 2). The mean H$\alpha$ luminosity of the hidden type 1 AGNs, $L_{\rm H\alpha}$, is $10^{40.73\pm 0.32}$ erg s*-1* while it is $10^{41.61\pm 0.38}$ erg s*-1* for the typical type 1 AGNs. In contrast, there is no significant difference in the distribution of FWHM of H$\alpha$ (Figure 2 right panel). Figure 3 presents the distribution of the hidden type 1 AGNs (red dots) and the typical type 1 AGNs (blue dots) in the $L_{\rm bol}$ - $M_{\rm BH}$ plane. The mean bolometric luminosity of the hidden type 1 AGNs is $10^{43.25\pm 0.28}$ erg s*-1* while that of typical type 1 AGNs is $10^{44.01\pm 0.33}$ erg s*-1*, indicating that the hidden type 1 AGNs have lower bolometric luminosity than the typical type 1 AGNs. In the case of black hole mass, the mean mass is $10^{7.20\pm 0.42}$M*⊙* and $10^{7.56\pm 0.44}$M*⊙*, respectively for the hidden type 1 AGNs and the typical type 1 AGNs. Combining mass and luminosity, the mean Eddington ratio of the hidden type 1 AGNs (i.e., 0.01 $\pm$ 0.01) is slightly lower than that of typical type 1 AGNs (i.e., 0.03 $\pm$ 0.02). Overall, the newly identified type 1 AGNs have relatively low luminosity while the width of H$\alpha$ is comparable to that of typical type 1 AGNs.

4.2. MBH-$\sigma_{*}$ relation of the hidden type 1 AGNs

We investigate whether the newly identified type 1 AGNs follow the black hole mass scaling relations as inactive galaxies and reverberation-mapped type 1 AGNs, using stellar velocity dispersion ($\sigma_{*}$) measurements from the SDSS DR7. We adopt the most recent calibration by Woo et al. (2015) for $\rm\sigma_{\rm H\beta}$-based MBH estimator as log (M{}_{\rm BH}/M_{\odot})$$=8.34 + 4.97 log ($\sigma_{*}$/200 km s*-1*) and for FWHMHβ-based MBH estimator as log (M${}_{\rm BH}/M_{\odot})$ $=$ 8.34 + 5.04 log ($\sigma_{*}$/200 km s*-1*). We find that the hidden type 1 AGNs locate slightly below the MBH-$\sigma_{*}$ relation (Figure 4 left panels) with an average offset of 0.28 $\pm$ 0.40 dex and 0.12 $\pm$ 0.41 dex, respectively for the $\rm\sigma_{H\beta}$ and $\rm FWHM_{H\beta}$ based $M_{\rm BH}$ estimator.

The offset may be interpreted in different ways. First these hidden type 1 AGNs may not follow the same MBH-$\sigma_{*}$ relation compared to typical type 1 AGNs. Second, since we use the H$\alpha$ velocity and luminosity as a proxy for H$\beta$ velocity and 5100Å luminosity in estimating black hole masses, systematic effects in the calibration of mass estimators may cause the offset, particularly at low mass scales where our targets are mainly located in the MBH-$\sigma_{*}$ plane. Note that the empirical correlation between $L_{5100}$ and $L_{\rm H\alpha}$ has not been directly calibrated for the hidden type 1 AGNs. Third, there is a possibility that the virial factor may not be same for the hidden type 1 AGNs. Four, stellar velocity dispersion may be overestimated particularly for the late type galaxies since the significant rotational velocity may broaden stellar lines observed through the 3" aperture of the SDSS spectroscopy (see e.g., Kang et al., 2013; Woo et al., 2013, 2015). Without more detailed data, e.g., dynamical black hole mass measurements and spatially resolved kinematics measurements, it is difficult to evaluate these different scenarios.

Instead, we further investigate whether the hidden type 1 AGNs have a different MBH-$\sigma_{*}$ relation compared to inactive and reverberation-mapped AGNs by following the joint-fit analysis presented by Woo et al. (2015). By combining the compilation of black hole mass and stellar velocity dispersion of inactive galaxies and reverberation-mapped AGNs from Woo et al. (2015) with those of our type 1 AGNs, we fit the MBH-$\sigma_{*}$ relation log (M{}_{BH}/M_{\odot})$$=$$\alpha + $\beta$ log ($\sigma_{*}$/200 km s*-1*), in order to determine the virial factor, intrinsic scatter, intercept($\alpha$) and slope($\beta$) based on the $\chi^{2}$ minimization as

[TABLE]

where $\mu=\log(M_{\rm BH}/M_{\odot})$ of inactive galaxies, $\mu_{VP}$ is the virial product of AGNs, and s represents log ($\sigma_{*}$/200 km s*-1*). For our type 1 AGNs, we calculate the virial product by diving MBH by the virial factor f in Eqs. 4 and 5. As a result of joint-fit analysis, we obtained a slope of 4.26 and log$f$ $=$ 1.03 for the $\sigma_{\rm H\alpha}$-based MBH estimator and a slope of 4.30 and log$f$ $=$ 0.30 for the $\rm FWHM_{H\alpha}$-based MBH estimator (Figure 4 right panels). The offset of the hidden type 1 AGNs from the MBH-$\sigma_{*}$ relation becomes negligible (i.e., 0.01 and 0.03 dex, respectively for $\sigma_{\rm H\alpha}$-based MBH and FWHMHα-based MBH), and the scatter of the sample (0.37 dex) is comparable to that of inactive galaxies and reverberation-mapped AGNs. Based on the joint fit analysis we find no strong evidence that the hidden type 1 AGNs have a different MBH-$\sigma_{*}$ relation compared to inactive galaxies or the reverberation-mapped AGNs.

We examine whether the measured stellar velocity dispersion may be affected by the rotational velocity by comparing face-on and edge-on late type galaxies in the MBH-$\sigma_{*}$ plane (see Figure 5). As expected, we find that the majority of edge-on galaxies are located below the MBH-$\sigma_{*}$ relation, suggesting that stellar velocity dispersion was overestimated due to the contribution of the rotational velocity, while face-on galaxies are more randomly distributed since the line-of-sight rotational velocity becomes much lower than that of edge-on galaxies.

In Figure 6, we compare the b/a ratio as a proxy for the inclination angle with the offset from the MBH-$\sigma_{*}$ relation using late-type galaxies. We find a weak trend that the average offset increases with increasing b/a ratio, indicating that more inclined late type galaxies have larger offset from the MBH-$\sigma_{*}$ relation. These results suggest that spatially resolved kinematics are required to overcome the inclination and rotation effect and to better investigate the MBH-$\sigma_{*}$ relation of the hidden type 1 AGNs. Nevertheless, we find that the hidden type 1 AGNs follow the same MBH-$\sigma_{*}$ relation as inactive galaxies and reverberation-mapped AGNs within the limitation of the current data.

5. Gas kinematics

5.1. Narrow Line Region gas kinematics

In this section, we investigate the gas kinematics in the NLR using [O iii] and [N ii] emission lines compared to stellar kinematics. If the kinematics of NLR gas is governed by the gravitational potential, the velocity dispersions of emission lines and stellar absorption lines should be comparable (Komossa et al., 2008). However, the ionized gas, e.g., [O iii] of AGNs typically shows larger velocity dispersion than stars (Woo et al., 2016, 2017), suggesting the presence of non-gravitational component.

First, we compare the velocity dispersion of [N ii] with stellar velocity dispersion, finding one-to-one relationship with a mean ratio, 0.99$\pm$0.20. This suggests that [N ii] gas does not show any non-gravitational component and the broadening of [N ii] gas can be entirely explained by the gravitational potential. In the case of [O iii], however, we find that the mean ratio between gas and stellar velocity dispersion is 1.22$\pm$0.50, indicating the presence of outflow component in [O iii]. When we separately calculate the velocity dispersion of narrow core and wing components, the core component shows a consistent velocity dispersion compared to stellar velocity dispersion with a mean ratio 0.96$\pm$0.29 while the velocity dispersion of the wing component of [O iii] is a factor of 2.11$\pm$0.80 higher than stellar velocity dispersion, representing non-gravitational kinematics, i.e., outflows (e.g., Boroson, 2005; Crenshaw et al., 2010; Woo et al., 2016, 2017).

Second, we investigate the velocity shift of each emission line with respect to the systemic velocity (Figure 7). When we use the total profile of [O iii], the mean velocity shift is -23$\pm$ 44 km s*-1* with a skewed distribution toward blueshift. Compared to the distribution of [O iii] velocity shift of type 2 AGNs (Woo et al., 2016), type 1 AGNs show a strong asymmetric distribution with a much larger number of AGNs with blueshifted [O iii] than AGNs with redshifted [O iii], presumably due to the fact that the direction of outflows is close to the line-of-sight (see the discussion on the number ratio between blueshifted and redshifted [O iii] in Bae & Woo (2016)). Once we use the narrow core and broad wing components separately, the mean velocity shift is -10$\pm$38 km s*-1and -102$\pm$108 km s-1*, respectively. These measurements indicate that the core component of [O iii] reflects the gravitational component while the wing component represents the outflow kinematics as discussed in the previous studies (e.g., Boroson, 2005; Bae & Woo, 2014; Woo et al., 2016). We also measure the velocity shift of narrow H$\alpha$, finding no significant velocity shift with a mean velocity shift 1$\pm$ 22 km s*-1* and a similar distribution compared to the [O iii] core component. In comparing velocity shifts of [O iii] and H$\alpha$, we find that the high ionization line is mainly blueshifted, which is consistent with the previous studies of type 1 AGNs (e.g., Boroson, 2005; Komossa et al., 2008; Zamanov et al., 2002). For instance, Boroson (2005) reported that 77% of the type 1 QSOs have blueshifted [O iii], while low ionization lines (e.g., [S ii], [N ii], [O ii]) do not show any significant velocity shift. Based on these results, we will use the wing component in [O iii] as a tracer of gas outflows to further investigate AGN-driven gas outflow in the next section.

5.2. AGN-driven outflow

5.2.1 Velocity shift and velocity dispersion of [O iii]

Along with [O iii] line dispersion, [O iii] velocity shift is also an outflow signature (Woo et al., 2016). As presented in the Figure 7, the velocity shift of the [O iii] total profile shows a blueward tail in its distribution. To estimate the velocity shift error in [O iii] of each galaxy, we performed Monte Carlo simulation to produce 100 mock spectra. We measure the velocity shift in each mock spectra and consider the 1$\sigma$ dispersion as the uncertainty of velocity shift. The mean measurement uncertainty of the [O iii] velocity shift is 17.0 \pm\ 27.2 km s*-1*. If we restrict the samples whose [O iii] velocity shifts are larger than 1$\sigma$ measurement uncertainty, 414 AGNs ($\sim$ 68%) show velocity shifts. Among the 414 AGNs, 318 AGNs ($\sim$ 77%) are blueshifted. In type 2 AGNs, the blueshifted fraction is 56%, when adopting the same measurement uncertainty using Monte Carlo Simulation (Woo et al., 2016). This result indicates that the blueshift is more common in type 1 AGNs than the type 2 AGNs (see also Crenshaw et al., 2010). The relatively common blueshifts in type 1 AGNs can be explained by the combined model of biconical outflows and dust extinction (Crenshaw et al., 2010; Bae & Woo, 2016). When we observe type 1 AGNs, it is more likely that the receding cone is obscured by the dusty stellar disk, resulting in blueshift of the observed [O iii] line profile. Therefore, [O iii] of type 1 AGNs are more likely to be blueshifted than type 2 AGNs due to the orientation effect.

To further investigate kinematics of the AGNs with gas outflow in the NLR, we separate [O iii] with a wing component (i.e., fitted with double Gaussian models; Group A) from [O iii] without a wing component (i.e., fitted with single Gaussian; Group B). Group A consists of 226 AGNs (about 34.2 % of total sample) while Group B is composed of 385 AGNs. The mean [O iii] velocity shift is -41 $\pm$ 54 km s*-1* and -12.1 $\pm$ 31.5 km s*-1*, respectively for Group A and Group B. Approximately 77% of AGNs in Group A shows velocity shift larger than 1$\sigma$ measurement uncertainty, while 63% of AGNs in Group B has velocity shift with >1$\sigma$ measurement uncertainty. If we count only these measurements (i.e., $>$ 1$\sigma$ error), the mean [O iii] velocity shift is -49.5 $\pm$ 56.6 km s*-1* and -18.9 $\pm$ 37.7 km s*-1*, respectively for Group A and B, indicating that Group A has on average larger velocity shift than Group B (see Figure 8).

In Figure 8 we compare the velocity shift and velocity dispersion of [O iii]. The normalized [O iii] line dispersions (log $\sigma_{\rm OIII}/\sigma_{*}$) of Group B are centered around 0, indicating that [O iii] line dispersion and stellar velocity dispersion are comparable. In contrast, Group A shows on average larger velocity shift and velocity dispersion than Group B. Compared to stellar velocity dispersion, [O iii] velocity dispersion is much larger, indicating the presence of non-gravitational kinematics, i.e., outflows. There is a trend that velocity shift becomes larger with increasing velocity dispersion as similarly found in type 2 AGNs (Woo et al., 2016), which is a characteristic feature of biconical outflows as demonstrated by Bae & Woo (2016).

5.2.2 Outflow fraction

We investigate the outflow fraction as a function of [O iii] luminosity and the Eddington ratio. Since the wing component of [O iii] indicates gas outflows, we use the presence of the wing in [O iii] as an evidence of outflows. The outflow fraction (i.e., the fraction of double Gaussian [O iii]) rapidly increases with increasing luminosity (Figure 9). At high luminosity, e.g., $L_{[\mbox{O\,{iii}}]}$ > $10^{41.0}$ erg s*-1*, the outflow fraction is over 85%, suggesting that most high-luminosity AGNs have strong outflows, which is consistent with the trend found in type 2 AGNs (Woo et al., 2016). On the other hand, the outflow fraction in low luminosity AGNs is lower than that of type 2 AGNs. For example, the outflow fraction of AGNs at $10^{40.0}$ erg s*-1* is 40% in the case of type 2 AGNs, while it is $\sim$20% in the hidden type 1 AGNs. It is not clear what is the origin of this discrepancy. A systematic comparison of the outflow fraction between type 1 and type 2 AGNs is yet to be available. The detection of a wing component in [O iii] is more difficult for lower luminosity AGNs since the wing component can be easily diluted by noise. Thus, the outflow fraction of the low luminosity AGNs should be regarded as a minimum (see also Woo et al., 2017). We note that our hidden type 1 AGN sample is far from complete at low luminosity, since a broad H$\alpha$ component will be much weaker and easily missed in identifying type 1 AGNs. Therefore, we interpret the difference of the outflow fraction at low luminosity is insignificant.

The outflow fraction also increases as Eddington ratio increases albeit with a shallower slope. Over the Eddington ratio range covered by our sample, i.e., log $(L_{\rm bol}/L_{\rm Edd})$ = -2.5 to -1.5, the outflow fraction increases from $\sim$30% to $\sim$50%. These results indicate that ionized gas outflows are directly connected to AGN activity and the outflow properties are qualitatively similar between type 1 and type 2 AGNs.

5.3. Gravitational Redshift

The velocity of BLR gas is relatively high up to several percent of the speed of light, suggesting that the the broad lines can be redshifted due to relativistic effects (e.g., Kollatschny, 2003). To investigate the gravitational redshift using the newly identified type 1 AGNs, we calculate the velocity shift of the broad H$\alpha$ with respect to the systemic velocity. In contrast to typical type 1 AGNs, for which systemic velocity cannot be accurately measured (Tremaine et al., 2014), we are able to use the luminosity-weighted stellar absorption lines to determine the systemic velocity since AGN continuum is too weak to dilute stellar absorption lines. Thus, we can decrease the systematic effect in measuring velocity shift of broad emission lines caused by the large uncertainty in systemic velocity.

In Figure 10, we present the velocity shift of the broad H$\alpha$ line as a function of the line width of H$\alpha$. The mean velocity shift is 115 $\pm$ 389 km s*-1* while 427 out of 611 AGNs show a redshifted H$\alpha$ ($\sim$ 70%), which may be interpreted as the gravitational redshift although the velocity shift may be also caused by other mechanisms such as non-gravitational kinematics (i.e., inflows and outflows) or orbital motions (Zheng & Sulentic, 1990). Note that the other 1/3 of the sample shows a blueshifted H$\alpha$. We find several extreme objects with a very large blue shift ($<$ -1000 km s*-1*) and a large velocity dispersion, while the H$\alpha$ line profile is very asymmetric. The nature of these objects is unclear and beyond the scope of the current study.

Using the SDSS quasar sample, Tremaine et al. (2014) statistically investigated the velocity shift of H$\beta$ to compare with the prediction of their BLR geometry models. The [O iii] narrow emission line was used to calculate the systemic velocity since stellar absorption lines are not detected due to strong AGN continuum. Comparing with their result (Figure 1 in Tremaine et al. (2014)), we obtain a similar result within the observed range of the broad H$\alpha$ velocity dispersion $<$ $\sim$3000 km s*-1* in our sample.

The geometry of the BLR is often assumed with spherical or disk models, and these models predict that the effect of the gravitational redshift increases with the BLR gas velocity (Tremaine et al., 2014). By taking the expected velocity shift of the broad H$\alpha$ as VHα = 1.5$\sigma_{\rm H\alpha}^{2}/c$ based on a spherical model of BLR from Tremaine et al. (2014), we compare the observed mean velocity shift with the prediction as a function of the line dispersion of broad H$\alpha$. Although there is a large scatter of velocity shift for given velocity dispersion, the trend of the mean velocity shift is consistent with the model, suggesting that the observed velocity shift is consistent with the gravitational redshift effect.

6. Discussion : Comparison with previous works

Although the presence of broad emission lines is a characteristic of optical type 1 AGNs, in some cases it is difficult to detect broad lines due to obscuration in the optical spectral range (e.g., red AGNs) or intrinsically low luminosity diluted by host galaxy emission (e.g., Simpson, 2005; Gilkman et al., 2007, 2013; Woo et al., 2014; Oh et al., 2015; Kim et al., 2015). Several studies tried to find hidden type 1 AGNs using different methods. For example, using the SDSS DR2 catalogue, Simpson (2005) found 63 objects (1.8% of total galaxies), which have a relatively broad H$\alpha$ line based on visual inspection. Similar works by Woo et al. (2014) and Oh et al. (2015) reported the detection of hidden type 1 AGNs using SDSS DR7.

We compare the selection methods and results between our work and the previous studies. In our pilot study (Woo et al., 2014), we identified type 1 AGNs at 0.02$<$z$<$0.05 by detecting a broad component in the H$\alpha$ line profile. Using the same type 2 AGN catalogue, we extend our previous work out to z=0.1, enlarging the sample size from 142 to 611 objects. However, we note that the selection scheme is improved based on the revised spectral decomposition and visual inspection. For example, we model the H$\alpha$ region more carefully using 3 different cases as described in Section 2.2, while Woo et al. (2014) used a simple fitting scheme using one broad and one narrow components for H$\alpha$. Thus, readers are suggested to use the results in this extended work.

Oh et al. (2015) searched type 1 AGNs among emission line galaxies at z$<$0.2 using the SDSS DR7 catalogue. First, they used the flux ratio at two narrow band regions, respectively representing a broad H$\alpha$ (6523-6543Å) and continuum (6460-6480Å), to trace the presence of a broad H$\alpha$ component. This flux ratio $F_{6533}/F_{6470}$ is useful to find broad H$\alpha$ candidates since it increases as the broad H$\alpha$ line becomes more prominent. Combining the signal-to-statistical noise (defined from continuum spectral range) and the flux ratio $F_{6533}/F_{6470}$, they defined the 1$\sigma$ demarcation line to select type 1 AGN candidates (see Figure 1 in Oh et al., 2015). Second, they used the area flux ratio, which is defined by the area of the red wing of the broad H$\alpha$ component, after excluding the [N ii] $\lambda$6584 emission line (see the blue + grey area in Figure 12), divided by the noise area defined with the 1$\sigma$ noise level of the continuum (blue area in Figure 12), in order to avoid unreliable detection of a broad H$\alpha$. Only if the area flux ratio is larger than 2, then they classified a target as type 1 AGNs. Using these criteria, they identified 1611 type 1 AGNs at z < 0.1. Among these AGNs, 882 objects are SDSS specClass =3 objects, which are already classified as type 1 AGNs in the SDSS, while the other 729 objects are SDSS specClass =2 (i.e., galaxy) objects. Thus, they found 729 hidden type 1 AGNs among type 2 or star forming galaxies by detecting a broad H$\alpha$, while in our study we found 611 type 1 AGNs from our type 2 AGN catalogue.

Among 729 type 1 AGNs identified by Oh et al. (2015), 230 objects are overlapped with our sample of 611 type 1 AGNs, while the other 499 objects are not identified as type 1 AGNs in our study. When we examine these 499 objects, the majority of them (412 objects) are located in the star forming region in the BPT diagram. Since we only used the type 2 AGN catalogue to search hidden type 1 AGNs, we simply did not investigate these 412 AGNs. The remaining 87 objects were not classified as type 1 AGNs in our study since a broad H$\alpha$ component was not required to fit the H$\alpha$ region. Instead, we were able to fit H$\alpha$ and [N ii] with double Gaussian models for the majority of these objects (70 objects, see the top panel in Figure 11) or single Gaussian models provide a good fit for 7 objects (see the bottom panel in Figure 11). For some cases (10 objects), the noise in the vicinity of H$\alpha$ region is too large to decide whether single or double Gaussian models provide a better fit.

Note that in our fitting procedure we often used double Gaussian models whenever a wing component in the narrow line profile is required. In particular, since we also examined the line profile of the [S ii] doublet, in order to check whether a wing component is necessary to fit the narrow component of H$\alpha$ and [N ii] lines, we tried to avoid a false detection of fitting the wing components of narrow lines as a broad H$\alpha$ component. In contrast, Oh et al. (2015) used single Gaussian models for narrow H$\alpha$ and [N ii] lines, hence they may have treated the wing components of the narrow lines as a broad H$\alpha$ component. In fact, the [S ii] line profile clearly shows the presence of a wing component for $\sim$a half of these objects (43/87) as demonstrated in Figure 11 (top panel). Since we used a double Gaussian profile for narrow lines (H$\alpha$, [N ii], and [S ii]), our fitting method for detecting a broad component of H$\alpha$ seems to be more conservative than that of Oh et al. (2015).

We also investigate why 315 objects among 611 objects in our hidden type 1 AGN sample were not classified as type 1 AGNs by Oh et al. (2015) . Many of them were excluded due to the area ratio criterion by Oh et al. (2015). For example, when we apply the area ratio cut defined by Oh et al. (2015) to the 315 AGNs, 136 objects do not satisfy the criterion while in the current study we conservatively concluded that a broad H$\alpha$ is present (see Figure 12). In particular, if the noise level is relatively large, the area flux ratio becomes below 2, hence, these AGNs will be rejected by the criterion defined by Oh et al. (2015). In Figure 12 we demonstrate two cases, each of which clearly shows a broad H$\alpha$ component in visual inspection while the area flux ratio does not satisfy the criterion (i.e., ratio$<$2). In addition, if the broad H$\alpha$ is blueshifted, the area ratio becomes also smaller since the area defined in the red side of the H$\alpha$ becomes smaller. For the remaining 179 objects among 315 objects, we find that the area flux ratio is slightly larger than 2 when we calculate it using our own noise estimates. Thus, it is not clear why these objects were rejected by Oh et al. (2015). It is possible that there is difference in noise calculation and in modeling continuum between these two studies.

In summary, although the selection method of Oh et al. (2015) provides consistent and quantitative criteria for identifying a hidden broad H$\alpha$ component, it is somewhat limited due to the fact that the wing components in the narrow lines are not included in their fitting, and that their area ratio cut may reject a weak and shallow broad H$\alpha$ component. More detailed studies on the detection scheme of a weak H$\alpha$ component is necessary to better detect hidden type 1 AGNs.

7. SUMMARY AND CONCLUSIONS

By detecting a broad component in the H$\alpha$ line profile, we conservatively identified a sample of 611 type 1 AGNs at 0.02 $<$ z $<$ 0.1, using the catalogue of Bae & Woo (2014), which provides a large sample of type 2 AGNs classified based on the emission line flux ratios. We increased the sample size of the hidden type 1 AGNs from 142 (Woo et al., 2014) to 611 by extending the redshift range out to z=0.1 and using a consistent and improved emission line analysis. The main findings are summarized as follows.

$\bullet$ The hidden type 1 AGNs have a similar range of the H$\alpha$ width, but on average lower luminosity compared to typical type 1 AGNs at the same redshift range, indicating that these hidden type 1 AGNs with a relatively weak AGN continuum is useful to study the properties of low luminosity AGNs and their host galaxies. The mean black hole mass estimated based on the H$\alpha$ FWHM and luminosity is log MBH =7.20$\pm$0.42, while the mean Eddington ratio is log Lbol/LEdd = -2.04$\pm$0.34.

$\bullet$ These AGNs seem to slightly offset from the MBH-$\sigma_{*}$ relation defined by inactive galaxies and reverberation-mapped AGNs (Woo et al., 2015), presumably due to the systematic difference in estimating black hole mass (i.e., depending on H$\beta$ or H$\alpha$) and potentially different virial factors between typical type 1 AGNs and the hidden type 1 AGNs with a weak AGN continuum. In contrast, when we perform a joint-fit analysis by combining quiescent galaxies, RM AGNs, and hidden type 1 AGNs, we find no significant difference between typical type 1 AGNs and hidden type 1 AGNs.

$\bullet$ By investigating the kinematics of ionized gas, we find that the velocity dispersion of [N ii] and the core component of [O iii] is roughly consistent with stellar velocity dispersion, indicating that the host galaxy gravitational potential is responsible for the broadening of these lines. In contrast, the wing component of [O iii] represents non-gravitational kinematics, i.e., outflows, which are consistent with the finding of type 2 AGNs (Woo et al., 2016, 2017).

$\bullet$ The velocity dispersion and velocity shift of [O iii] show strong non-gravitational kinematics, i.e., outflows. The fraction of AGNs with a wing component in [O iii] strongly increases with AGN luminosity, suggesting that the non-gravitational kinematics are directly connected to AGN activity.

$\bullet$ The line-of-sight velocity and velocity dispersion of the ionized gas in type 1 AGNs is on average larger than that of type 2 AGNs (e.g. Woo et al., 2016), which is consistent with the biconical outflow models and orientation effect (e.g., Bae & Woo, 2016). Although the uncertainty of the velocity shift of [O iii] is large, we find that there are more AGNs with blueshifted [O iii] than AGNs with redshifted [O iii].

Based on these results, we conclude that the hidden type 1 AGNs follow the general characteristics of typical broad line AGNs except for low luminosity and low Eddington ratio. More detailed studies on the detection scheme and the completeness of the broad line AGNs may provide a better understanding of these hidden type 1 AGNs, which can be used as a unique channel for studying AGN unification and black hole-galaxy connection.

We thank the anonymous referee for valuable comments. Support for this work was provided by the National Research Foundation of Korea grant funded by the Korea government (No. 2016R1A2B3011457 and No. 2010-0027910).