고대 은하의 금속성 진화와 질소 풍부 현상의 새로운 해석

초록

THESAN‑ZOOM 고해상도 시뮬레이션을 이용해 3 < z < 12 구간의 은하 질량‑금속성 관계(MZR)와 기본 금속성 관계(FMR)를 조사했다. 가스와 별의 MZR는 z≈12부터 이미 완만한 곡선을 이루며, z=3까지 금속 함량이 각각 ≈0.2 dex, ≈0.4 dex 상승한다. 금속 보유 효율이 주요 조절인자이며, 저질량 은하( M★ ≲ 10⁹ M⊙ )에서는 가스‑FMR이 붕괴·역전한다. 질소‑산소 비(N/O)와 C/O, Fe/O는 관측과 일치하고, 질소‑풍부 은하(NRG, log (N/O) > ‑0.6)는 폭발적인 별 형성 후 AGB 별이 주도하는 질소 주입과 차등적인 은하풍에 의해 100 Myr 이내에 급격히 상승한다. NRG는 주계열 아래에 위치하고 EW(Hβ) ≲ 40 Å를 보이지만, 관측된 고 EW NRG와의 차이는 새로운 별 폭발 단계에서 기존 질소‑풍부 가스를 재조명함으로써 해소된다.

상세 분석

THESAN‑ZOOM은 기존 THESAN 대규모 시뮬레이션을 기반으로, 4×, 8×, 16× 세 가지 해상도(최고 1.4 × 10² M⊙ 입자 질량)에서 방사선‑유체역학과 비평형 열화학을 동시에 구현한다. 특히 SMUGGLE 기반의 다채널 피드백(광이온화, 방사압, 항성풍, 초신성, 초기 피드백)을 포함해 별 형성 효율을 파괴적인 수준으로 억제하고, 은하 규모의 가스 흐름을 정밀히 추적한다. 이러한 물리적 세팅은 은하 내 금속 보유율을 직접 계산하게 하며, 금속 손실이 주로 은하풍에 의해 질량‑의존적으로 조절된다는 결론을 뒷받침한다.

MZR 분석에서는 가스와 별의 금속성이 모두 ‘곡선형’으로 진화한다는 점이 핵심이다. 저질량 은하에서는 가스‑MZR의 기울기가 완만해지며, 이는 외부 저금속 가스 유입이 지속적으로 금속을 희석하기 때문이다. 반면 별‑MZR는 금속 보유율이 시간이 지남에 따라 증가하면서 더 뚜렷한 상승을 보인다. 이러한 차이는 금속이 별 형성 후 장기적으로 은하 내에 머무는 비율이 가스와 별 각각에서 다르게 작용함을 의미한다.

FMR에 대해서는, 기존 연구에서 제시된 ‘금속‑질량‑SFR 삼중 관계’가 별에 대해서는 유지되지만, 가스에 대해서는 저질량 은하에서 급격히 붕괴하고 심지어 역전한다는 새로운 현상을 발견했다. 이는 가스‑FMR이 SFR에 의해 금속이 희석되는 정도를 반영하지만, 저질량 은하에서는 유입 가스가 지속적으로 저금속을 공급해 SFR과 무관하게 금속을 낮추기 때문이다.

핵심적인 화학적 발견은 N/O 비율의 급격한 변동이다. 시뮬레이션은 별 폭발 직후 초신성 풍으로 인해 α‑원소(O 등)가 대량 배출되는 반면, AGB 별이 풍부한 질소를 방출하기까지는 ≈30–100 Myr가 소요된다. 이 시차가 ‘질소‑풍부 은하(NRG)’를 만들며, N/O가 ‑0.6 dex 이상으로 급등한다. 차등적인 은하풍(초신성 풍은 주로 가벼운 가스를, AGB 풍은 무거운 가스를 남김)도 N/O 상승에 기여한다. NRG는 주계열 아래에 위치하고, Hβ 등가폭(EW) ≲ 40 Å를 보이는데, 이는 별 폭발 직후 가스가 거의 비워진 상태에서 남은 AGB 질소가 희소 가스에 섞이기 때문이다. 관측된 고 EW(NRG)와의 차이는, 실제 은하가 새로운 별 폭발을 시작하면서 기존 질소‑풍부 가스를 재조명하고, 그 결과 EW가 급격히 상승한다는 시나리오로 설명된다.

시뮬레이션은 또한 C/O와 Fe/O 비율이 관측된 고‑z 은하와 일치함을 보여, 금속 생산 메커니즘(핵융합, 초신성, AGB)과 피드백이 일관된 방식으로 작동함을 확인한다. 전체적으로, 금속 보유 효율, 가스 유입·유출, 그리고 별 형성 주기의 복합적인 상호작용이 초기 은하의 화학 진화를 좌우한다는 결론에 도달한다.