SMBH와 은하의 인과 관계 전환: M·–σ₀ 관계의 새로운 해석

초록

이 논문은 NIHAO 시뮬레이션을 이용해 초대질량 블랙홀(M·)과 은하 중심 속도분산(σ₀) 사이의 인과 관계가 별 형성 피크 전후에 반전된다는 사실을 밝힌다. 결과는 고‑z에서 기존 스케일링 관계를 그대로 적용하면 블랙홀 질량이 과대평가될 수 있음을 시사한다.

상세 분석

본 연구는 관측에서 도출된 인과 관계를 시간적 차원에서 검증하기 위해 NIHAO 프로젝트의 고해상도 줌‑인 수치 시뮬레이션을 활용하였다. 시뮬레이션은 Gasoline2 코드를 기반으로 하며, ΛCDM 우주론, 블랙홀 시드, Bondi‑Hoyle‑Lyttleton 급착 및 Eddington 제한을 포함한 피드백 모델을 적용한다. 각 은하에 대해 M·, σ₀, Rₑ, M* 및 sSFR 다섯 개 변수를 매 스냅샷마다 추출하였다.

인과 탐지는 베이지안 점수 기반 정확한 사후 분포를 계산하는 방식으로, 5개의 변수 조합에 대해 29,281개의 가능한 DAG를 전부 평가하였다. 관측 연구(Paper I)와 달리 시뮬레이션은 시간 순서를 알 수 있기 때문에, 별 형성 피크(T_peak) 전후의 인과 흐름을 직접 비교할 수 있다.

주요 결과는 다음과 같다. 별 형성 피크 이전, 즉 활발한 가스 공급이 존재하는 시기에는 σ₀ → M· 보다 M· → σ₀ 가 우세하였다. 즉, 블랙홀 활동이 은하 중심의 동역학을 주도한다는 의미다. 반대로 피크 이후, 별 형성이 급감하고 가스가 고갈되는 단계에서는 σ₀ → M· 가 지배적으로 나타나, 은하 구조가 블랙홀 성장의 주요 원천이 된다. 이는 관측에서 ‘spiral 은하는 SMBH가 원인, elliptical 은하는 SMBH가 결과’라는 결론과 일맥상통하지만, 시계열 관점에서 인과 방향이 전환되는 구체적 메커니즘을 제시한다.

이러한 전환은 고‑z에서 SMBH 질량을 추정할 때 기존의 M·–σ₀ 스케일링을 그대로 적용하면, 실제보다 1–2 dex 과대평가될 위험이 있음을 보여준다. 저자들은 시뮬레이션 기반 인과 정보를 반영한 새로운 스케일링 관계를 도출했으며, 이를 적용하면 고‑z 블랙홀 질량이 현지(z = 0) 관계와 비교해 약 100배 낮게 추정된다.

논문은 또한 시뮬레이션의 한계—예를 들어, ISM의 다상성 미해상도, 별 형성 버스트의 과소평가, 그리고 샘플이 질량에 따라 비균등하게 선택된 점—을 인정하고, 향후 고해상도 ISM 모델과 더 큰 샘플을 이용한 인과 탐색이 필요함을 제언한다.