Unveiling the population of massive quenched galaxies at $z\ge2$ in the COLIBRE simulations -- I. Galaxy demographics

기술 요약: COLIBRE 시뮬레이션에서의 $z \ge 2$ 고적색편이 거대 정지 은하(Massive Quenched Galaxies) 인구 규명 – I. 은하 인구 통계

문제 제로 (Problem Statement)
최근 제임스 웹 우주 망원경(JWST)의 관측은 높은 적색편이($z \gtrsim 2$)에서 거대 정지 은하(Massive Quenched Galaxies; MQGs; $M_\star \gtrsim 10^{10} M_\odot$)의 상당한 인구를 발견했다. 추정된 수밀도와 급격한 형성 시간대는 계층적 $\Lambda$CDM 모델 및 기존 은하 형성 시뮬레이션과 상당한 긴장 관계를 형성했는데, 기존 모델들은 우주 역사 초기 단계에서 이토록 거대한 수동적(passive) 시스템을 재현하는 데 어려움을 겪어왔다. 이후의 관측 분석을 통해 수밀도가 하향 조정되었으나, 현재의 모델들과 최신 JWST 데이터 사이에는 MQG의 풍부도, 질량 함수, 그리고 물리적 특성에 관한 불일치가 여전히 존재한다.

방법론 (Methodology)
본 연구는 이전의 대규모 부피 시뮬레이션들과 비교하여 서브그리드 물리(sub-grid physics)의 핵심적인 혁신을 도입한 새로운 COLIBRE 우주론적 유체역학 시뮬레이션을 활용한다. 분석은 피듀셜(fiducial) L200m6 시뮬레이션(평균 가스 입자 질량 $m_g \sim 1.84 \times 10^6 M_\odot$인 $200 \, \text{cMpc}^3$ 부피)을 중심으로 수행되며, 부피, 해상도, 그리고 활동 은하핵(AGN) 피드백 모델(열적 피드백 vs. 하이브리드 열적+제트 피드백)을 변화시킨 비교 실행들을 함께 다룬다.

주요 방법론적 특징은 다음과 같다:

• 서브그리드 물리 (Sub-grid Physics): 차가운 분자 상(phase)과 먼지 입자의 명시적 추적 (HYBRID-CHIMES 모델을 통해), 자기 일관적 화학적 농축, 그리고 업데이트된 AGN 피드백 처방(초에딩(super-Eddington) 강착 및 제트 피드백 포함)을 도입하였다.
• 샘플 선택: MQG는 $z \ge 2$에서 $M_\star > 10^{10} M_\odot$이고 비정규 별 형성률(sSFR)이 $< 0.2/t_{\text{age}}$인 은하로 정의된다.
• 포워드 모델링 (Forward Modeling): 관측 불확실성을 해결하기 위해, 저자들은 시뮬레이션의 별 질량과 별 형성률(SFR)을 가우시안 분포($\sigma = 0.3$ dex)로 컨볼루션하여 SED 피팅 및 광도 적색편이에서 발생하는 오차를 모사하였다.
• 특성 분석: 본 연구는 수밀도, 별 질량 함수(SMFs), 별 형성사(SFHs), 형성 및 정지 시간($t_{50}$ 및 $t_q$), 먼지 및 분자 수소($H_2$) 분율, 그리고 운동학($v_\star/\sigma_\star$)을 조사한다.

주요 기여 및 결과 (Key Contributions and Results)

1. 수밀도 및 해상도 의존성:

   • COLIBRE 내 MQG의 수밀도는 질량 해상도와 AGN 피드백 구현 방식 모두에 의존한다. 더 높은 해상도의 시뮬레이션(예: L200m7)은 더 적은 수의 정지 은하를 생성하는데, 이는 높은 해상도가 블랙홀 주변의 가스를 더 밀도 있게 만들고, 이로 인해 정지를 위해 더 강력한 피드백이 필요하게 되며, AGN 피드백이 덜 돌발적(less bursty)으로 변하기 때문이다.
   • 하이브리드 AGN 피드백(열적 + 제트)은 피듀셜 열적 모델에 비해 $z \gtrsim -2$에서 유의미하게 적은 수의 정지 은하를 생성하는데, 이는 하이브리드 시나리오에서 블랙홀이 초기 우주에서 더 느리게 성장하기 때문이다.
   • 관측 불확실성(0.3 dex)을 포워드 모델링했을 때, 피듀셜 L200m6 시뮬레이션이 예측하는 수밀도는 최신 JWST 분광 및 광도 제약 조건과 대체로 일치하며, 이는 이전 연구들에서 나타났던 긴장 관계를 크게 완화시킨다.

2. 별 질량 함수 (SMF):

   • COLIBRE의 MQG 별 질량 함수는 $M_\star \approx 10^{10.5} M_\odot$ 부근에서 정점을 찍으며, 이는 Baker et al. (2025c) 및 Stevenson et al. (2026)의 관측 데이터와 일치한다.
   • 저질량 및 고질량 끝단을 동시에 맞추는 데 어려움을 겪는 다른 모델들(예: EAGLE, IllustrisTNG, Simba)과 달리, COLIBRE는 관측 오차를 고려할 때 전체 질량 범위에 걸쳐 SMF를 재현한다. 이는 차가운 가스 상을 명시적으로 모델링함으로써 더 높은 가스 밀도, 초에딩 블랙홀 강착, 그리고 결과적으로 더 효율적인 초기 정지를 가능하게 했기 때문이다.

3. 별 형성사 (SFHs):

   • 형성 시간 ($t_{50}$): COLIBRE의 MQG는 빠르게 형성되며, 중앙값 형성 시간은 $t_{50} \approx 0.5 - 1.5$ Gyr이다. 이는 연장된 형성과 그 뒤를 잇는 급격한 정지를 나타낸다.
   • 정지 시간 ($t_q$): 시뮬레이션은 주로 강력한 스타버스트 에피소드에 의해 구동되는 급격한 정지($t_q \lesssim 0.6$ Gyr)를 예측한다.
   • SFH는 종종 최종 정지 전 여러 번의 스타버스트 에피소드를 보여주는데, 이는 관측적 복원을 복잡하게 만드는 특징이지만 JWST 데이터로부터 얻은 추론들과 질적으로 일치한다.

4. 먼지 및 분자 가스 함량:

   • 이 모델은 MQG가 거대 별 형성 은하들보다 먼지 및 $H_2$ 분율이 1 dex 이상 낮을 것이라고 예측한다.
   • 이러한 낮은 분율은 가용한 희박한 관측 추정치 및 상한선과 일치하며, 이는 피드백 과정이 이러한 시스템에서 가스 저장량을 효율적으로 고갈시키고 먼지 성장을 억제한다는 견해를 뒷받침한다.

5. 크기 및 운동학:

   • COLIBRE의 MQG와 거대 별 형성 시스템은 전반적으로 유사한 크기와 운동학($v_\star/\sigma_\star$)을 보인다. 이는 $z \ge 2$에서 유의미한 크기 또는 형태적 변형이 정지 상태의 전제 조건이라기보다는, 정지 이후에 발생함을 시사한다.

의의 (Significance)
본 논문은 최첨단 서브그리드 물리(특히 먼지, 차가운 가스, 그리고 AGN 피드백)를 활용한 COLIBRE 시뮬레이션이 높은 적색편이에서의 거대 정지 은하의 인구 통계를 견고하게 재현할 수 있음을 확립한다. 이 연구는 현실적인 관측 불확실성을 포워드 모델링할 때 시뮬레이션과 JWST 관측 사이의 외견상 긴장이 대부분 완화된다는 것을 보여준다. MQG의 풍부도, 질량 조립, 그리고 정지 시간에 대한 일관된 프레임워크를 제공함으로써, 본 연구는 이 초기 수동적 시스템들의 기원을 이해하기 위한 실행 가능한 경로를 제시하며, 이들의 급격한 형성과 정지를 구동하는 구체적인 물리적 메커니즘을 조사할 Paper II를 위한 토대를 마련한다.