초기 우주 은하 크기 조절은 중심 별폭발과 내부 소거에 달렸다

초록

THESAN‑ZOOM 고해상도 방사선‑수소역학 시뮬레이션을 이용해 3 < z < 13, 10⁶–10¹⁰ M⊙ 범위의 은하 크기 진화를 조사하였다. 별 형성 주계열 위에 있는 은하는 중심 별폭발(central starburst)로 인해 반경이 일시적으로 축소되고, 이후 내부 소거(inside‑out quenching) 과정에서 별 형성이 외곽으로 퍼지며 크기가 다시 증가한다. 이 순환적 성장으로 고‑z에서 양의 크기‑질량 관계가 유지되며, 정규화는 z≈6에서 전환되는 이중 거듭 제곱법칙을 따른다. Hα는 UV 연속보다 평균 1.7배 크게 보이며, 이는 LyC 방출이 중심 별폭발에서 급증해 스트뮐렌 구( Strömgren sphere) 크기가 확장되기 때문이라는 해석이 제시된다.

상세 분석

본 연구는 THESAN‑ZOOM 시뮬레이션의 뛰어난 공간·질량 해상도(최고 1.4 × 10² M⊙)와 다상성( multi‑phase) ISM 모델을 활용해, 기존 대용량 시뮬레이션에서 흔히 사용되는 유효 상태 방정식(effective equation‑of‑state) 접근법이 놓치는 ‘버스티(starburst‑driven)‘와 ‘내부 소거(inside‑out quenching)’ 메커니즘을 정량화한다. 은하가 주계열 위에 있을 때, 방사선 피드백과 초신성·방사압이 중앙 밀도 높은 가스 구역을 급격히 압축해 별 형성 효율을 순간적으로 10‑100배 상승시킨다. 이때 별 형성 반경(Rₛₓ)은 급격히 감소해 반경‑질량 관계가 일시적으로 아래로 이동한다. 별 형성률이 감소하면, 방사선과 초신성 피드백이 가스 흐름을 외부로 밀어내고, 남은 가스는 은하 외곽에 재분포한다. 결과적으로 새로운 별 형성이 외부에서 진행되면서 반경이 다시 팽창하고, 은하는 ‘크기‑질량 진동(oscillatory)’ 패턴을 보인다.

시뮬레이션은 z > 6에서 대부분의 은하가 상승하는 별 형성 역사를 보이며, 이때는 ‘컴팩션(compaction) 단계’에 머무른다. 이러한 단계는 크기‑질량 정규화가 두 개의 거듭 제곱법칙(α₁, α₂)으로 기술될 수 있음을 의미한다. z≈6에서 전환점이 나타나는 이유는, 이 시점부터 은하가 내부 소거를 겪으며 별 형성 영역이 점차 외부로 이동하기 때문이다.

관측적 측면에서는 Hα와 UV 연속의 크기 차이를 직접 측정했으며, Hα가 평균 1.7배 더 크게 나타난다. 기존 해석은 ‘inside‑out 성장’에 기인한다고 보았지만, 저자들은 LyC 방출이 중심 별폭발에서 급증해 스트뮐렌 구의 반경이 크게 늘어나는 현상을 제시한다. 스트뮐렌 반경 Rₛ ∝ (Qₕ / nₑ²)¹ᐟ³ (Qₕ: LyC 광자 방출률, nₑ: 전자 밀도) 를 이용해, 고‑z 은하의 높은 Qₕ와 낮은 nₑ가 Hα 확장을 자연스럽게 설명한다.

또한, 대용량 시뮬레이션이 예측하는 ‘음의 크기‑질량 기울기(negative slope)’와 관측(및 THESAN‑ZOOM)에서 보이는 ‘양의 기울기(positive slope)’ 사이의 차이를, 버스티 별 형성을 억제하는 유효 EOS와 물리적 피드백이 결합된 다상성 ISM 모델의 차이로 해석한다. 이는 차세대 관측(JWST, ELT)과 이론 모델 간의 격차를 해소하는 데 중요한 단서를 제공한다.