ALMA로 본 은하 중심의 억태양질량 블랙홀 측정

초록

NGC 4061의 중심에 위치한 초대질량 블랙홀(MBH≈1.2×10⁹ M⊙)을 고해상도 ALMA ¹²CO(2‑1) 데이터와 KinMS 전방 모델링으로 정밀 측정하였다. 0’’.16 × 0’’.13의 빔으로 블랙홀 영향반경을 직접 해상도에 포함시켰으며, 은하의 I‑밴드 질량‑광도비(M/L≈3.5)와도 동시에 추정했다. 이전에 보고된 서로 상이한 속도분산값을 해소하고, 고질량 영역에서의 MBH‑σ 관계 검증에 중요한 사례를 제공한다.

상세 분석

본 연구는 초대질량 블랙홀(MBH) 측정을 위해 ALMA Band 6에서 ¹²CO(2‑1) 전이선을 이용한 고해상도 분자 가스 동역학 분석을 수행하였다. 기존에 제시된 두 개의 서로 다른 별속도분산(σ≈459 km s⁻¹와 σ≈290 km s⁻¹)으로 인해 NGC 4061의 블랙홀 질량 추정에 큰 불확실성이 존재했으나, 본 논문은 0’’.16 × 0’’.13(≈80 pc × 70 pc)의 합성빔을 확보해 블랙홀 영향반경(rₛₚₕ≈0’’.15)을 충분히 해상도에 포함시켰다.

데이터는 Cycle 6(2018.1.00397.S)와 Cycle 7(2019.1.00036.S)의 12 m 배열을 결합해 총 6 시간 이상의 온소스 시간을 확보했으며, 10 km s⁻¹ 채널폭으로 ³⁰ km s⁻¹ 이하의 속도 해상도를 달성하였다. 연속파는 uv‑plane에서 선형 파워‑로우 모델로 빼고, ¹²CO(2‑1) 라인만을 포함한 3‑차원 데이터 큐브를 생성하였다. MOM0, MOM1, MOM2 맵은 모두 원형 회전 디스크 형태를 보이며, 고정된 PA≈175°와 거의 엣지‑온(i≈80°)인 기울기를 갖는다.

은하의 별질량 분포는 HST F814W 이미지에서 Multi‑Gaussian Expansion(MGE) 기법을 적용해 2‑D 광도 프로파일을 추출하고, I‑밴드 질량‑광도비(M/L_F814W)를 자유 파라미터로 두어 중력 포텐셜을 구성하였다. 가스 표면밀도는 두 가지 방식(데이터‑드리븐 비대칭 모델과 분석적 퐁크션)으로 모델링했으며, 두 경우 모두 KinMS(Kinematic Molecular Simulation) 프레임워크를 이용해 전체 데이터 큐브를 전방 시뮬레이션했다.

베이즈 추정에서는 블랙홀 질량(MBH), M/L, 디스크 기울기, 중심 좌표 등을 파라미터로 두고, MCMC 샘플링을 통해 사후 확률분포를 도출하였다. 통계적 불확실성은 ΔM_BH≈+0.08/−0.10 × 10⁹ M⊙이며, 시스템적 오차(모델 선택, M/L 기울기, 거리 불확실성 등)를 포함하면 총 ±0.43 × 10⁹ M⊙가 된다. 최종 결과는 M_BH = (1.17 ± 0.08_stat ± 0.43_sys) × 10⁹ M⊙, M/L_F814W = 3.46 ± 0.07_stat ± 0.10_sys M⊙/L⊙이다.

주요 검증 절차로는(1) 가스 디스크의 중심 좌표를 자유롭게 두어도 M_BH가 변하지 않음, (2) M/L 프로파일에 방사형 변화를 허용해도 M_BH 편차가 <5%, (3) 비대칭 가스 분포를 반영한 모델과 단순 원형 가우시안 모델이 일관된 M_BH를 제공함을 확인했다. 또한, 블랙홀 질량이 σ≈459 km s⁻¹에서 예상되는 1.2×10¹⁰ M⊙와는 크게 차이 나며, 이는 이전 σ 측정이 과대평가된 결과임을 시사한다.

이 연구는 고해상도 ALMA 분자 가스 동역학이 수백 메가파섹스 거리에서도 블랙홀 영향반경을 직접 해상도에 포함시켜 정확한 M_BH 측정이 가능함을 입증한다. 특히, 고질량(>10⁹ M⊙) 영역에서 기존 별속도분산 기반 추정치와의 불일치를 해결하고, M_BH‑σ 관계의 상단을 보다 견고히 하는 데 기여한다. 향후 비슷한 구조를 가진 초대질량 타원은 ALMA를 활용한 가스 동역학 측정으로 대규모 샘플을 구축하면, 블랙홀 질량 함수와 은하 진화 모델에 중요한 제약을 제공할 것으로 기대된다.