Black Hole Feedback, Galaxy Quenching and Outflows at Cosmic Dawn: Analysis of the SEEDZ Simulations

기술 요약: 우주 초기의 블랙홀 피드백, 은하 퀘칭(Quenching), 그리고 유출(Outflows)

문제 제기
제임스 웹 우주 망원경(JWST)의 최근 관측은 초기 우주의 10억 년 이내($z \gtrsim 6$, $z \sim 12.3$까지 확장)에 존재하는 $10^6 M_\odot$을 초과하는 질량을 가진 거대 질량 블랙홀(MBH)의 상당한 인구 집단을 밝혀냈다. 이러한 발견은 압축 천체의 형성과 성장에 관한 현재의 이론적 모델에 도전 과제를 던진다. 구체적으로, 블랙홀이 주어진 시간 내에 어떻게 초기 씨앗 질량으로부터 그토록 큰 규모로 성장할 수 있었는지는 여전히 불분명하다. 특히, 고적색편이 관측은 국부 우주의 관계($M_{BH}/M_{star} \simeq 10^{-4}$)와 비교했을 때 과도하게 큰 블랙홀 대비 항성 질량 비($M_{BH}/M_{star} \gtrsim 10^{-2}$)를 나타내며, 이는 블랙홀 씨앗의 형성 이후에 항성 성분이 구축되는 독특한 진화 경로를 시사한다. SEEDZ 시뮬레이션 제품군은 초기 은하 및 MBH 진화 경로를 조사하기 위해 설계되었으며, 특히 우주 초기에 MBH 성장을 조절하거나 중단시키는 물리적 메커니즘에 초점을 맞추었다.

방법론
본 연구는 움직이는 격자 코드인 Aremo2를 사용하여 수행된 완전한 3D 유체 역학적 우주론적 시뮬레이션 세트인 SEEDZ 시뮬레이션 제품군을 분석한다. 시뮬레이션은 $40 \, \text{cMpc} \, h^{-1}$의 부모 박스 내에서 세 개의 서로 다른 영역(Normal1, Normal2, Rarepeak)에 집중하는 "줌인(zoom-in)" 접근 방식을 사용한다. 고해해상도 영역의 암흑 물질 입자 질량은 약 $\approx 6.9 \times 10^4 M_\odot h^{-1}$이며, 최소 셀 크기는 $z=15$에서 $\approx 6 \, \text{pc}$(물리적 단위)이다.

주요 수치 및 물리적 구성 요소는 다음과 같다:

• 씨앗 형성(Seeding): 시뮬레이션은 "가벼운 씨앗"(Population III 별의 잔해)과 "무거운 씨앗"(조밀한 클러스터에서의 폭주 충돌, 초거대 별 또는 준성(quasi-star)을 통해 형성됨)을 모두 모델링한다. 무거운 씨앗은 멱함수 분포에 따라 $5 \times 10^3 M_\odot$에서 $10^5 M_\odot$ 사이의 질량으로 형성된다.
• 강착(Accretion): 블랙홀 강착은 커널 반경을 통해 시뮬레이션 해상도에 맞춰 조정된 Bondi-Hoyle-Lyttleton 정식화(formalism)를 사용하여 모델링되며, 고각운동량 가스의 과도한 강착을 방지하기 위해 가스 와도(vorticity)에 의해 수정된다.
• 피드백(Feedback): 열 에너지는 강착 반경 내의 가스 셀에 등방적으로 주입된다. 복사 효율($\epsilon$)은 동적으로 계산된다. 이는 표준적인 아디아바틱(sub-Eddington) 값을 따르지만($\sim 0.1$), 초에딩겐(super-Eddington) 강착 중에는 광자 트래핑(photon trapping)으로 인해 유효 효율이 감소하고 피드백이 약화되는 "슬림 디스크(slim disc)" 모델로 전환된다.
• 역학(Dynamics): 소프트닝 길이(softening lengths)에 의한 중력 상호작作用의 억제를 설명하기 위한 역학적 마찰 보정과 근접성 및 상대 속도에 기반한 병합 기준을 포함한다.
• 범위: 분석은 $z=12.5$까지의 제품군 내 가장 질량이 큰 블랙홀들의 성장 이력에 초점을 맞춘다.

주요 결과

1. 성장 및 질량 한계: 시뮬레이션의 가장 질량이 큰 블랙홀들은 $z=12.5$까지 $\sim 10^6 M_\odot$의 질량으로 성장한다. 이러한 성장은 5~30 Myr 동안 지속되는 짧은 초에딩겐 강착 폭 중에 발생한다. 최종 질량은 초기 씨앗 질량에 의해 결정되는 것이 아니라 환경 조건에 의해 결정된다.
2. 피드백에 의한 퀘칭(Quenching): MBH 성장을 중단시키는 주요 메커니즘은 초신성이나 가스 고갈이 아니라 블랙홀 자체로부터 오는 피드백이다. 가장 활발하게 강착하는 MBH의 경우, 주입된 열 에너지는 호스트 헤일로의 결합 에너지(binding energy)를 초과한다. 이는 호스트 헤일로로부터 은하간 매질(IGM)으로 가스를 완전히 배출하는 결과를 초래하여, 효과적으로 은하를 퀘칭시킨다.
3. 유출 특성: 강착 중단 시점에 가스 유출은 $\sim 10^3 \, \text{km s}^{-1}$의 속도와 $10^7 \, \text{K}$의 온도에 도달한다. 이러한 유출은 냉각되기 전까지 $\sim 20 \, \text{kpc}$ 규모까지 확장된다. 일부 헤일로는 20~100 Myr의 시간 척도 내에 가스 유입률이 회복되는 징후를 보이기도 하지만, 이 시뮬레이션의 블랙홀들은 시뮬레이션 기간 내에 두 번째 주요 성장 단계를 경험하지 않는다.
4. 헤일로 질량 의존성: 최대 블랙홀 질량은 피드백 모델과 호스트 헤일로의 결합 퍼텐셜(binding potential)의 조합에 의해 실질적으로 결정된다. 시뮬레이션된 헤일로(질량 최대 $\sim 10^{10} M_\odot$)에서, 가스를 결합 해제(unbind)하는 데 필요한 피드백 에너지는 블랙홀 성장을 약 $10^6 M_\odot$로 제한한다.
5. 항성 형성사: 가장 질량이 큰 블랙홀을 수용하는 헤일로들은 무거운 씨앗이 형성되기 전까지 매우 적은 항성 형성을 보이며, 이를 통해 효율적인 강착이 가능하다. 일단 강착이 시작되면, 결과적인 AGN 피드백은 추가적인 항성 형성을 억제한다. 반면, 지배적인 항성 인구를 가진 헤일로들은 지속적인 항성 형성을 보이지만, 질량이 큰 BH 호스트에서 보이는 급격하고 방해받지 않는 강착은 나타나지 않는다.

의의 및 결론
본 논문은 구현된 피드백 모델 하에서, 초기 $\Lambda$CDM 우주 내의 거대 질량 블랙홀들이 강착만으로는 $z=12.5$까지 $\sim 10^6 M_\odot$를 초과하기 어려울 것이라고 결론짓는다. 급격한 강착으로부터 발생하는 피드백은 호스트 헤일로의 가스 저장고를 결합 해제할 만큼 충분히 강력하며, 이로 인해 블랙홀의 성장이 자기 조절(self-regulation)된다.

저자들은 이 한계를 초과하는 고적색편이의 거대 질량 블랙홀 인구 집단을 설명하기 위해 두 가지 해석을 제안한다:

1. 더 약한 피드백: 초기 우주의 MBH 피드백은 모델링된 것보다 훨씬 덜 효율적일 수 있다 (예: 성장이 강착 중심이 아닌 병합 중심이거나, 피드백이 좁은 영역 내에 갇혀 있는 경우).
2. 2단계 형성: 관측된 고적색하이 은하들은 MBH가 처음에 가스 저장고를 배출함으로써 호스트 은하를 퀘칭시키고, 이후 가스가 풍부한 은하들과의 병합 또는 우주 거대 구조(cosmic web)로부터의 강착을 통해 바리온 저장고를 재구성하여 후속 성장을 허용하는 2단계 과정을 거쳐 형성되었을 수 있다.

저자들은 현재의 해상도가 블랙홀 병합을 완전히 해결하는 데 한계가 있음을 언급하며, 이는 후기 단계에서 지배적인 성장 메커니즘이 될 수 있다. 향-후 더 높은 해상도의 시뮬레이션은 블랙홀 병합을 더 잘 해결하고, 특히 JWST에서 관측된 "리틀 레드 닷(Little Red Dot)" 시스템과 관련하여 피드백 매개변수를 더 잘 규명하는 데 사용될 계획이다. 결과적으로, 피드백이 극도로 비효фектив적이거나 국소적으로 제한되지 않는 한, $z=12.5$ 이전의 블랙홀 최대 질량은 $10^6 M_\odot$를 크게 초과하지 않을 것임을 시사한다.