Black Hole Masses, Scaling Relations, and The Co-Evolution of Black Holes and Galaxies

Abstract

초거대 블랙홀의 성장과 그 모은하 간의 공동 진화를 이해하고자 다음과 같은 일련의 연구들을 수행하였다. 먼저 이 연구에서 핵심적인 물리량인 블랙홀 질량을 정확히 측정하기 위해, 수소 방출선의 한 시점 분광 관측을 통해 측정되는 비리얼 블랙홀 질량의 불확실성에 대해 체계적으로 연구하고 이러한 질량 측정법에 중요한 비리얼 계수를 추정하여 보다 향상된 질량 측정법을 제시하였다. 높은 적색편이에 존재하는 활동성 은하의 블랙홀 질량 측정에 유용한 탄소 방출선 기반 질량 추정법의 정확도를 향상시켰다. 이러한 경험과 기술을 바탕으로 먼 우주에 존재하는 블랙홀의 질량과 그 모은하의 팽대부 밝기를 측정하여 이들의 물리적 상관 관계에 대한 시간 진화를 연구하였다.

가까운 우주에 존재하는 9개의 활동성 은하에 대해 빛 반향 분포 측정법을 통해 얻어진 균질하고 고품질의 분광 자료에 다성분 분광 분해 기법을 적용하여, 은하의 밝기 변화에서 기인하는 임의적 오차(12\%)와 방출선의 선폭 차이에서 기인하는 체계적 오차(60\%)를 측정하고 이를 보정하는 새로운 질량 측정기법을 개발하였다. 또한 이 연구에 기반한 질량 측정치가 이론적인 비리얼 정리와 크기-밝기 관계에서 예측되는 값과 잘 부합함을 보였다.

최근 연구 문헌에서 정립된 활동성/비활동성 은하 자료에 다양한 선형회귀분석법을 적용하여 블랙홀질량-별속도분산 상관 관계식과 비리얼 계수를 유도하였다. 이러한 결과로부터 샘플 선택에서 기인하는 차이와 선형회귀분석의 독립변수 선택에서 기인하는 차이로 인해 현재 비리얼 계수 결정법에는 두배 가량의 체계적 오차가 존재함을 보였다. 또한 선형회귀분석법 중 FITEXY와 베이지안 방법론이 블랙홀 척도관계 연구에 가장 적합함을 보였다.

빛 반향 분포 측정법이 적용된 활동성 은하의 최신 샘플과 자료를 바탕으로 다성분 분광 분해 기법을 적용해 탄소 방출선 기반 블랙홀 질량 측정법의 정확도를 향상시켰다. 이러한 새 측정법이 기존의 측정법에 비해 얼마나 차이를 주는지 슬로안 분광전천탐사 자료에서 얻어진 퀘이사 목록에 적용하여 연구하였다.

지금으로부터 대략 40억년 전과 60억년 전에 존재하는 52개의 활동성 은하 샘플에 대해 Keck 지상 망원경과 허블 우주 망원경을 통해 관측하여 얻어진 고품질의 분광 자료와 고분해능의 측광 자료에 정교한 스펙트럼 분석 기법과 이미지 분석 기법을 적용하여 정확한 블랙홀 질량과 은하 팽대부 밝기 측정치를 얻었다. 이 결과에 기반해 먼 우주에 존재하는 블랙홀이 현재 우주의 블랙홀에 비해 상대적으로 작은 은하의 중심에 위치하고 있다는 관측적 증거를 제시하였다. 샘플 선택 효과를 보정하도록 설계된 몬테 칼로 시뮬레이션 기법을 적용해 블랙홀과 은하의 공동 진화 양상을 통계적으로 유도하여 블랙홀의 성장이 은하 진화에 선행한다는 것을 시사하였다.

To better understand the cosmic growth of supermassive black hole (BH) population and the co-evolution with their host galaxies, we first perform three studies for improving determination of AGN BH masses, i.e., 1) investigating the uncertainties of \Hb\ single-epoch (SE) virial BH mass estimates, 2) constraining the virial factor of the BH mass estimators, and 3) improving the calibration of the \CIV\ SE virial BH mass estimators for high-redshift AGNs. Then we present a study for direct observational constraints on the redshift evolution of the black hole mass $-$ bulge luminosity relation.

We investigate the calibration and uncertainties of the black hole mass estimates based on the single-epoch method, using homogeneous and high-quality 428 multi-epoch spectra obtained by the Lick AGN Monitoring Project 2008 for 9 local Seyfert~1 galaxies with black hole masses $<10^{8}~{\rm M}_{\odot}$. By decomposing the spectra into their AGN and stellar components, we study the variability of the single-epoch H$\beta$ line width (full width at half-maximum intensity, \FWHM

or dispersion, \sigmaline) and of the AGN continuum luminosity at 5100~\AA~($L_{\rm 5100}$). From the distribution of the ``virial products'' ($\propto$~\FWHM$^2$ $L_{\rm 5100}^{0.5}$ or \sigmaline$^2$ $L_{\rm 5100}^{0.5}$) measured from SE spectra, we estimate the uncertainty due to the combined variability as $\sim0.05$ dex (12\%). This is subdominant with respect to the total uncertainty in SE mass estimates, which is dominated by uncertainties in the size-luminosity relation and virial coefficient, and is estimated to be $\sim 0.46$ dex (factor of $\sim3$). By comparing the \Hb~line profile of the SE, mean, and root-mean-square (rms) spectra, we find that the \Hb\ line is broader in the mean (and SE) spectra than in the rms spectra by $\sim0.1$ dex (25\%) for our sample with \FWHM\ $<3000$ \kms. This result is at variance with larger mass black holes where the difference is typically found to be much less than $0.1$ dex. To correct for this systematic difference of the \Hb\ line profile, we introduce a line-width dependent virial factor, resulting in a recalibration of SE black hole mass estimators for low-mass AGNs.

Determining the virial factor of the broad-line region gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the \msigma~relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studies. We investigate the origin of the difference by measuring the \msigma~relation using an updated galaxy sample from the literature, and explore the dependence of the virial factor on various fitting methods. We find that the discrepancy is primarily caused by the sample selection, while the difference stemming from the various regression methods is marginal. However, we generally prefer the \texttt{FITEXY} and \texttt{Bayesian} estimators based on Monte Carlo simulations for the \msigma~relation. In addition, the choice of independent variable in the regression leads to $\sim 0.2$ dex variation in the virial factor inferred from the calibration process. Based on the determined virial factor, we present the updated \msigma~relation of local active galaxies.

We present the single-epoch black hole mass estimators based on the \CIV~$\lambda$1549 broad emission line, using the updated sample of the reverberation-mapped AGNs and high-quality UV spectra. By performing multi-component spectral fitting analysis, we measure the \CIV\ line widths (FWHM$_{\rm CIV}$ and line dispersion, $\sigma_{\rm CIV}$) and the continuum luminosity at 1350 \AA~($L_{1350}$) to calibrate the \CIV-based mass estimators. By comparing with the \Hb\ reverberation-based masses, we provide new mass estimators with the best-fit relationships, i.e., $M_{\rm BH} \propto L_{1350}^{0.50\pm0.07} \sigma_{\rm CIV}^{2}$ and $M_{\rm BH} \propto L_{1350}^{0.52\pm0.09}{\rm FWHM}_{\rm CIV}^{0.56\pm0.48}$. The new \CIV-based mass estimators show significant mass-dependent systematic difference compared to the estimators commonly used in the literature. Using the published Sloan Digital Sky Survey QSO catalog, we show that the black hole mass of high-redshift QSOs decreases on average by $\sim 0.25$ dex if our recipe is adopted.

We investigate the cosmic evolution of the black hole mass -- bulge luminosity relation with a sample of 52 moderate luminosity galaxies at $z \sim 0.36$ and $z \sim 0.57$, corresponding look-back times of 4 and 6 Gyrs. By employing multicomponent spectral and structural decomposition methods to the high-quality Keck spectra and high-resolution \HST\ images, black hole masses (\mbh) are estimated using the \Hb\ broad emission line with the 5100\AA\ nuclear luminosity, and bulge luminosities (\Lbul) are derived from the surface photometry. Taking into account selection effects, we obtain the evolution trend of the form $M_{\rm BH} / L_{\rm bul} \propto (1+z)^{1.9 \pm 0.7}$, which is consistent with previous studies based on the bulge luminosity, as well as on the bulge mass and stellar velocity dispersion. This result indicates that black holes in the non-local universe live in smaller bulges than today, thus implying that black holes grow first and then the host galaxies follow up in the context of the co-evolution of black holes and galaxies.