13억년 동안 은하의 케넌벳‑슈미트 관계를 밝히는 COLIBRE 시뮬레이션

초록

COLIBRE 우주론적 수치 시뮬레이션은 비평형 화학, 먼지 진화, 10 K까지의 냉각을 직접 구현해 원자수소(H I)와 분자수소(H₂) 면밀한 KS 관계를 kpc 규모에서 0 ≤ z ≤ 8까지 추적한다. z≈0에서 관측된 H I·H₂ KS 곡선을 재현하고, 산란이 별표면밀도, sSFR, 금속성에 강하게 연관됨을 보인다. H I KS는 저질량 은하에서 가팔라지고, H₂ KS는 sSFR이 높은 은하에서 정상화가 상승한다. H₂ 고갈시간은 z=0에서 z=8까지 약 20배 감소하며, 이는 주로 가스 금속성 감소에 기인한다. 결과는 관측과 일치한다.

상세 분석

본 연구는 최신 우주론적 수치 시뮬레이션인 COLIBRE를 이용해 Kennicutt‑Schmidt(KS) 관계의 물리적 기원을 다각도로 검증한다. 가장 큰 강점은 비평형 H/He 화학망과 실시간 먼지 입자 성장·소멸 모델을 포함해 온도 10 K까지 냉각을 직접 계산한다는 점이다. 이로써 H I와 H₂의 질량 분포를 사후 처리 없이 시뮬레이션 자체에서 예측할 수 있다. 별 형성은 자유 낙하 시간당 1 %의 고정 효율을 갖는 Schmidt 법칙을 따르며, 불안정성 기준(열·난류 압력 대비 중력)으로 가스 입자를 선별한다. 해상도 의존성을 최소화하기 위해 불안정성 기준을 적용하고, 고해상도 실행에서는 더 높은 밀도에서 별이 형성되는 것을 확인한다.

z≈0에서 시뮬레이션은 관측된 H I KS(Σ_SFR ∝ Σ_HI^N, N≈1.5)와 H₂ KS(N≈1) 모두를 재현한다. 특히 산란이 별표면밀도(Σ_), 지역 sSFR, 가스 금속성(Z_gas)과 강하게 상관함을 보여, 기존 연구에서 제시된 “두 번째 매개변수” 역할을 물리적으로 입증한다. 저질량 은하(M_ < 10^9 M_⊙)에서는 H I KS의 기울기가 급격히 가팔라져, 원자수소가 별 형성의 주요 연료가 됨을 시사한다. 반면, sSFR이 높은 은하에서는 H₂ KS의 정규화가 상승해, 동일한 Σ_H₂에 대해 더 높은 Σ_SFR을 보인다.

시간 진화 측면에서 H₂ 고갈시간(t_dep = M_H₂/Ṁ_*)은 z=0에서 약 2 Gyr에서 z=8에서는 ≈0.1 Gyr까지 감소한다. 저금속성 환경에서는 H₂ 형성 효율이 감소하고, 같은 SFR을 유지하기 위해 더 많은 H I가 필요하므로 H I/Σ_SFR 비율이 증가한다. 이는 관측된 고‑z 은하들의 높은 sSFR과 낮은 금속성에 대한 자연스러운 설명을 제공한다. 또한, sSFR와 H₂ 고갈시간 사이의 반비례 관계가 시뮬레이션 전 범위(0 ≤ z ≤ 5)에서 관측과 일치한다는 점은 모델의 신뢰성을 높인다.

시뮬레이션은 또한 Σ_*와 Z_gas가 KS 산란에 미치는 영향을 정량화한다. Σ_*가 높을수록 중력적 안정성이 강화돼 같은 가스 밀도에서 별 형성 효율이 낮아지며, 이는 KS 곡선의 아래쪽으로 이동한다. 반대로, 금속성이 높을수록 H₂ 형성이 촉진돼 Σ_SFR가 상승한다. 이러한 다변량 관계는 미래 관측(예: JWST, ALMA)에서 KS 관계를 해석할 때 필수적인 교정 인자로 활용될 수 있다.

마지막으로, 해상도 수렴 테스트와 다양한 볼륨·해상도 조합을 통해 결과의 견고함을 검증하였다. 특히, kpc 스케일에서의 KS 관계는 해상도 1 kpc 이하에서도 수렴함을 확인했으며, 이는 대규모 우주 시뮬레이션에서도 신뢰할 수 있는 결과임을 의미한다.