Supermassive Black Hole Growth in Massive Galaxies at Cosmic Dawn

기술 요약: 우주 초기 거대 은하 내 초대질량 블랙홀의 성장

문제 제기
제임스 웹 우주 망원경(JWST)의 최근 관측은 높은 별 질량($M_\star \gtrsim 10^8 M_\odot$)과 별 형성률(SFR)을 가진 우주 초기($z \gtrsim 9$)의 거대하고 밝은 은하 집단을 발견했다. 이들 중 여러 은하는 활동성 은하핵(AGN)을 시사하는 스펙트럼 특징을 보이며, 이는 국부적인 스케일링 관계에 비해 과도하게 큰 질량을 가진 초대질량 블랙홀(SMBH)의 존재를 암시한다. 표준 은하 형성 모델은 이러한 SMBH의 급격한 성장과 관찰된 높은 별 형성 효율을 설명하는 데 어려움을 겪고 있다. 이론적 제안들은 무거운 씨앗(heavy seeds)이나 초에딩(super-Eddington) 강착을 호출하지만, 초기 거대 은하의 난류적이고 다상적인 성간 매질(ISM) 내에서 SMBH 성장, 별 피드백, 그리고 AGN 피드백 사이의 상호작용은 여전히 제대로 이해되지 않고 있다.

방법론
저자들은 ramses 적응형 격자 세분화(AMR) 코드를 사용한 11개의 우주론적 유체역학 "줌인(zoom-in)" 시뮬레이션 세트를 수행하였다. 시뮬레이션은 $z \sim 9$까지 $\sim 10^{11} M_\odot$ 할로로 붕괴하는 특정 과밀 영역에 초점을 맞추었으며, 여기에는 $M_\star \gtrsim 10^8 M_\odot$인 거대 은하(Massive Dawn Galaxy, MDG)가 포함된다.

주요 방법론적 구성 요소는 다음과 같다:

• 별 형성 모델: Kennicutt-Schmidt 법칙을 사용하는 표준 시뮬레이션과 달리, 본 연구는 물리적으로 동기화된 난류 기반의 서브그리드 다중 자유 낙하(multi-freefall) 모델을 채택하였다. 이를 통해 ISM의 난류 운동 에너지에 의해 구동되어 $\gtrsim 10\%$에 도달할 수 있는 가변적인 국지적 별 형성 효율(SFE)을 구현하였다.
• SMBH 구현: SMBH는 $10^8 M_\odot$ 할로 내에서 국지적 가스 및 별 질량이 $10^5 M_\odot$를 초과할 때 $z \sim 15$에서 생성되는 싱크 입자(sink particles)로 모델링되었다. 씨앗 질량으로는 $10^4 M_\odot$와 $10^5 M_\odot$가 테스트되었다.
• 강착 물리: 강착은 해결되지 않은 밀도 요동을 고려하기 위한 부스트 인자(boost factor)가 포함된 수정된 Bondi-Hoyle-Lyttleton 공식을 따른다. 강착률은 $\lambda_{\text{edd}}$(값: 1, 3, $\infty$ [한계 없음])를 통해 시뮬레이션 전반에 걸쳐 변화하는 에딩 한계에 의해 제한된다.
• 피드백 메커니즘:
  • 별 피드백: 초기형의 열적 광이온화(thermal photoionization)와 후기형의 제2형 초신성(Type II supernovae)을 포함한다. ISM 구조에 미치는 영향을 격리하기 위해 광이온화를 제거한 "낮은 별 피드백" 변형 모델도 테스트되었다.
  • AGN 피드백: 싱크 구체(sink sphere) 내로 등방성 열 에너지를 주입하는 방식으로 모델링되었으며, 결합 효율 $\epsilon_c$를 테스트하였다 (0% 및 15%).
• 시뮬레이션 세트: 11개의 실행 모델은 AGN 피드백 강도, 씨앗 질량, 최대 강착 한계, 그리고 별 피드백 효율을 체계적으로 변화시켜 SMBH와 은하의 공진화에 미치는 영향을 격리하였다.

주요 결과

1. ISM 주도적 확률적 강착 (포식과 기아):
   표준 MDG 환경에서 별 피드백은 차갑고 밀도가 높은($n_H \sim 10^3 \text{ cm}^{-3}$) 구역과 뜨겁고 희박한($n_H \sim 10 \text{ cm}^{-3}$) 구역으로 구성된 난류적 다상 ISM을 생성한다. SMBH 싱크 구체는 이러한 상들을 확률적으로 샘플링한다.

   • 포식 모드 (Feast Mode): 싱크가 차갑고 밀도가 높은 구역에 진입하면 음속이 감소하여 강착이 (초-)에딩 한계에 도달할 수 있게 된다.
   • 기아 모드 (Starvation Mode): 싱크가 뜨겁고 희박한 구역에 진입하면 음속이 상승하고, Bondi 강착률이 크게 떨어진다.
   • 표본 모델에서 SMBH는 별 형성이 시작된 후($t \sim 300$ Myr) 약 50%의 시간을 기아 모드($\alpha \approx 0.5$)로 보내며, 이는 매끄러운 ISM에 비해 평균 성장률을 절반으로 억제하는 효과를 낸다.

2. 자기 조절(Self-Regulation)의 조건:
   AGN 피드백은 블랙홀이 그 가열이 싱크 구체 내의 냉각을 압도할 만큼 충분히 커질 때($M_{\text{sink}} > M_{\text{sink,crit}}^{\text{cool}}$)만 SMBH 성장을 자기 조절할 수 있다.

   • 이 임계 질량은 가스 밀도와 금속 함량에 따라 달라진다. 난류적인 MDG 환경에서 SMBH가 가스를 할로로부터 해방시키기 위한 폭주 가열(runaway heating)을 달성하려면, 뜨겁고 희박한 포켓에서 보내는 시간 중 $\gtrsim 50\%$ 동안 임계치를 초과해야 한다.
   • 에딩 제한 강착($\lambda_{\text{edd}}=1$)과 표준 씨앗 질량을 가진 시뮬레이션은 종종 자기 조절에 도달하지 못하는데, 이는 금속 함량이 높아짐에 따라 임계 질량 임계치가 SMBH가 성장하는 속도보다 빠르게 상승하거나, SMBH가 임계치 이상에서 머무는 시간이 불충분하기 때문이다.
   • 자기 조절은 더 높은 씨앗 질량($10^5 M_\odot$) 또는 초에딩 강착($\lambda_{\text{edd}} \ge 3$)이 제공될 때 달성되며, 관측된 SMBH 질량과 일치하려면 열적 결합 효율이 충분히 낮아야 한다 ($\epsilon_c \lesssim 0.15\%$).

3. 은하 규모의 퀜칭(Quenching) 결여:
   저적색편이 모델의 예상과는 반대로, 이 시뮬레이션에서의 AGN 피드백는 전역적인 별 형성률(SFR)을 억제하거나 호스트 은하의 별 질량 성장을 중단시키지 못한다.

   • 할로는 우주 필라멘트로부터 유입되는 차가운 가스를 유입하는 비율($\dot{M}_{\text{halo}}$)이 SMBH의 최대 성장률($\dot{M}_{\text{SMBH}}$)을 상당히 초과한다.
   • 할로의 성장 시간 규모($\tau_{\text{halo}}$)가 SMBH의 성장 시간 규모와 비슷하거나 더 짧기 때문에, 신선한 차가운 가스가 지속적으로 ISM을 보충하여 AGN 피드백으로부터 오는 가열을 압도한다. 따라서 은하는 SMBH 활동과 관계없이 지수적으로 계속 성장한다.

4. 관측과의 비교:
   우주 초기의 과도한 질량을 가진 SMBH 집단(예: UHZ-1, GN-z11)을 가장 잘 설명하는 시뮬레이션 실행 모델은 super_edd_no_AGN 런($M_{\text{seed}} \sim 10^4 M_\odot$, $\lambda_{\text{edd}} = 3$, AGN 피드백 없음)이다.

   • 자기 조절되는 SMBH를 통해 관측된 질량에 도달하려면 열적 결합 효율이 극도로 낮아야 한다 ($\epsilon_c \lesssim 0.15\%$).
   • 결과는 초에딩 강착이 $10^6 M_\odot$만큼 높은 씨앗 질량을 요구하지 않고도 SMBH가 자기 조절 단계에 진입할 수 있을 만큼 성장하는 데 필요할 가능성이 높음을 시사한다.

의의 및 결론
본 논문은 우주 초기 거대 은하 내 SMBH 성장이 AGN 피드백 단독이라기보다는, 별 피드백에 의해 유도된 난류적이고 다상적인 ISM의 확률적 특성에 의해 주로 조절된다고 주장한다. 저자들의 결론은 다음과 같다:

• 이러한 환경의 SMBH는 자주 뜨겁고 희박한 가스 포켓에 의해 "기아" 상태에 놓이며, 이는 초에딩 강착을 유지하지 못하는 한 성장을 제한한다.
• AGN 피드백는 할로의 급격한 가스 유입이 연료 공급을 지속적으로 보충하기 때문에 이 시대의 은하 규모 별 형성을 억제하는 데 효과적이지 않다.
• JWST에서 관찰된 거대한 SMBH를 달성하려면, 은하가 더 안정적인 저적색편이 상태로 진 evolution 하기 전의 상승하는 임계 냉각 질량 임계치를 극복하기 위해 초에딩 강착과 비효율적인 AGN 피드백 결합, 또는 매우 무거운 씨앗의 조합이 필요할 가능성이 높다.

저자들은 본 결과가 해상도 의존적임을 언급하며, 이러한 발견을 완전히 검증하고 피드백 증착의 물리적 척도를 정교화하기 위해서는 비평형 열화학, 복사 전달, 그리고 더 높은 공간 해상도($\sim 1$ pc)를 포함하는 향후 연구가 필요하다고 밝혔다.