JWST와 PHANGS가 밝힌 별 형성 임베디드 단계의 짧은 지속시간

초록

PHANGS‑JWST 37개 근거리 은하를 대상으로 21 µm 중적외선과 CO, Hα 데이터를 결합해 별 형성의 임베디드 단계(먼지에 가려진 단계) 지속시간을 측정했다. 결과는 대부분의 은하에서 이 단계가 4 Myr 이하, 경우에 따라 1 Myr 미만으로 매우 짧으며, 바가 있는 나선 은하에서는 다소 길고 저질량 불규칙 은하에서는 거의 없다는 점을 보여준다. 이는 초신성 이전의 방사선, 풍압, 풍 등 사전‑초신성 피드백이 구름을 빠르게 분산시킨다는 증거다.

상세 분석

본 연구는 PHANGS‑ALMA, PHANGS‑Hα, PHANGS‑JWST (MIRI 21 µm) 데이터를 30–250 pc 해상도로 통합한 뒤, Kruijssen & Longmore(2014)와 Kruijssen et al.(2018)에서 제시된 “spatial decorrelation” 방법을 적용하였다. 구체적으로, CO(2‑1) 강도가 최대인 위치와 21 µm 강도가 최대인 위치 사이의 공간적 겹침 정도를 통계적으로 측정하고, 이를 CO‑Hα, Hα‑21 µm, CO‑21 µm 삼중 트레이서 조합에 대해 각각 모델링하였다. 모델은 각 트레이서가 나타내는 물리적 단계(분자 구름, 임베디드 별 형성, 노출된 H II 영역)의 평균 지속시간과 겹침 시간을 파라미터화한다. 베이즈 추정과 MCMC 샘플링을 통해 각 은하별 최적 파라미터를 도출했으며, 은하 전체 평균값은 t_embedded ≈ 2–3 Myr, t_exposed ≈ 1–2 Myr, t_overlap ≈ 0.5–1 Myr 로 나타났다.

특히, 은하 형태와 질량에 따른 차이를 조사했을 때, 바가 있는 대형 나선 은하(NGC 1300, NGC 1365 등)에서는 평균 t_embedded이 3–4 Myr까지 늘어나는 반면, 질량이 10⁹ M⊙ 이하인 불규칙·왜소 은하에서는 1 Myr 이하로 급감한다. 이는 바 구조가 가스 흐름을 집중시켜 GMC의 질량과 난류를 증가시키고, 따라서 별 형성 전 단계에서 더 오래 머무르게 만든다는 물리적 해석과 일치한다. 또한, 금속량이 낮은 은하(12+log(O/H) < 8.3)에서는 t_embedded이 짧아지는 경향이 관측되었으며, 이는 낮은 금속성이 먼지 양을 감소시켜 방사선 피드백이 더 효율적으로 구름을 파괴한다는 가설을 뒷받침한다.

이와 동시에, 21 µm 밴드가 CO와 Hα와 동시에 보이는 “overlap” 구간이 짧은 이유는, 중간 적외선이 주로 매우 작은 입자(VSG)의 열에 의해 구동되며, 이 입자들은 별이 방출하는 UV/광학 광자에 즉각적으로 반응하기 때문이다. 따라서 21 µm은 실제로 “임베디드” 단계의 실시간 지표가 되며, 전통적인 24 µm(Spitzer)보다 더 정확한 타이밍을 제공한다는 점이 강조된다.

통계적 검증을 위해 각 은하별 잔차 분석과 부트스트랩 재샘플링을 수행했으며, 전체 샘플에 대한 평균 χ²는 1.2 수준으로 모델이 데이터에 잘 맞는 것을 확인했다. 또한, 해상도 제한(≤ 150 pc) 이하에서는 CO‑21 µm 겹침이 감지되지 않아 15개의 은하가 제외되었으며, 이는 현재 관측 능력으로는 구름 규모보다 작은 피드백 구조를 분리하기 어려운 한계점을 보여준다.

결론적으로, 본 연구는 JWST의 고해상도 중적외선이 사전‑초신성 피드백을 정량화하는 데 핵심적인 역할을 하며, 은하 규모에서 피드백 효율이 구조·질량·금속성에 따라 어떻게 변하는지를 최초로 통계적으로 입증했다.