버스티 은하가 작은 붉은 점을 품을 수 있는가

초록

JWST가 발견한 고‑z 작은 붉은 점(LRD)들이 저질량, 버스티 은하에서 활발히 성장하는 초대질량 블랙홀을 포함한다는 가설을 검증하기 위해 FIRE‑2 시뮬레이션에 두 가지 흡수 모델(중력 토크 구동 모델과 자유 낙하 모델)을 적용하였다. 두 모델 모두 중앙 가스 공급이 1 % 이하의 효율로 자유 낙하한다면 관측된 LRD 수와 밝기 함수를 충분히 재현하지만, 저광도 AGN를 과다 예측한다. 초과를 완화하려면 초극초광속(Eddington 제한)으로 성장하는 $M_{\rm BH}\gtrsim2\times10^{5},M_\odot$ 블랙홀이 $M_\star\gtrsim2\times10^{7},M_\odot$ 은하에 존재해야 한다는 시나리오가 제시된다.

상세 분석

본 연구는 JWST가 보고한 $z\gtrsim4$의 풍부한 작은 붉은 점(LRD) 집단이 실제로는 저질량, 버스티(high‑z bursty) 은하 내에서 활발히 성장하는 블랙홀에 의해 구동된 AGN임을 검증하고자 한다. 이를 위해 FIRE‑2 프로젝트의 고‑z 줌‑인 시뮬레이션(z = 5–7, $M_{\rm halo}=10^{10-12}M_\odot$)을 사용했으며, 시뮬레이션 자체는 블랙홀 물리와 AGN 피드백을 포함하지 않는다. 저자들은 가상의 블랙홀을 각 은하 중심(최대 별 밀도 지점)에 삽입하고, 두 가지 포스트‑프로세싱 흡수 모델을 적용했다. 첫 번째는 Hopkins & Quataert(2011)의 중력 토크 구동 모델(GTDA)으로, 은하 내 디스크 질량, 가스 분포, 토크 효율 $\epsilon_T$ 등을 이용해 $\dot M_{\rm torque}= \epsilon_T f_d (M_{\rm tot,b}/10^9M_\odot)^{5/2}(R/100\ {\rm pc})^{-3/2}\dots$ 형태의 식을 사용한다. 두 번째는 단순 자유 낙하 모델로, 반경 $R$ 안의 가스 질량 $M_{\rm gas}(<R)$와 자유 낙하 시간 $t_{\rm ff}$를 이용해 $\dot M_{\rm ff}= \epsilon_{\rm ff} M_{\rm gas}(<R)/t_{\rm ff}$ 로 정의한다. 여기서 $\epsilon_{\rm ff}$는 가스가 자유 낙하당 얼마나 효율적으로 블랙홀에 공급되는지를 나타내는 파라미터이며, 0–1 사이의 값을 가진다.

핵심 결과는 두 모델 모두 $\epsilon_{\rm ff}\lesssim0.01$ 혹은 $\epsilon_T\sim5$(GTDA의 높은 정규화)일 때, $L_{\rm bol}\sim10^{45},{\rm erg,s^{-1}}$ 수준의 고광도 AGN를 충분히 생산한다는 점이다. 그러나 저광도 구간($L_{\rm bol}\lesssim10^{44},{\rm erg,s^{-1}}$)에서는 관측된 LRD 수보다 과다하게 예측한다. 이는 시뮬레이션이 중앙 가스 공급을 과도하게 유지하거나, 실제 은하에서 블랙홀 피드백이 가스 흐름을 억제하는 효과를 무시했기 때문일 가능성이 있다.

저자들은 과잉 예측을 완화하기 위한 두 가지 물리적 시나리오를 제시한다. 첫째, 블랙홀 질량이 $M_{\rm BH}\gtrsim2\times10^{5}M_\odot$ 수준으로 비교적 크게 시작하면, 동일한 가스 공급에서도 Eddington 한계에 도달해 초극초광속(초과) 성장하게 된다. 이 경우 실제 방출되는 복사량은 Eddington 제한에 의해 억제되므로, 관측된 LRD의 밝기 분포와 일치한다. 둘째, 은하의 UV 광도 분포와 비교했을 때, LRD가 호스트 은하보다 약간 더 밝은 경우가 많아, 블랙홀 성장 단계가 은하 성장보다 앞서 진행된다는 ‘조기 과성장’ 시나리오가 가능하다.

또한, 저자들은 모델 민감도 검증을 위해 반경 $R$(50–200 pc), 블랙홀–별 질량 비율($M_{\rm BH}=0.01M_\star$ vs. 다른 스케일링), 그리고 $\epsilon_T$, $\epsilon_{\rm ff}$의 로그 정규산포(σ=0.5) 등을 변형했다. 결과는 전반적으로 고광도 AGN 예측에 크게 영향을 주지 않으며, 저광도 과잉은 주로 가스 공급 지속성에 기인한다는 결론을 내렸다.

마지막으로, AGN 피드백이 포함되지 않은 점을 인정하고, 향후 고해상도 ‘하이퍼‑리파인먼트’ 시뮬레이션을 통해 중앙 가스 흐름과 피드백 상호작용을 직접 모델링할 필요성을 강조한다.