자기중력과 클러핑이 만든 은하계 분자수소 부족 현상 해소

초록

관측에 비해 태양권 근처에서 시뮬레이션이 H₂ 비율을 2‑4배 낮게 예측한다는 문제를, GHOSDT 고해상도 MHD 시뮬레이션과 다양한 물리 모델(광이온화, 초신성 에너지, 자기장, 서브그리드 클러핑)을 통해 조사하였다. 해상도 1 M⊙ 이하에서도 Rₘₒₗ이 수렴하지 않으며, 클러핑 모델을 도입했을 때만 관측값 Rₘₒₗ≈0.25–0.42를 재현한다는 결론을 얻었다.

상세 분석

본 논문은 은하계 패치 시뮬레이션에서 관측된 HI‑to‑H₂ 전이 압력(P_DE/k_B≈10⁴ K cm⁻³)과 실제 분자수소 비율(R_mol) 사이의 불일치를 정량적으로 분석한다. GHOSDT 시뮬레이션은 질량 해상도 100 M⊙부터 0.25 M⊙까지, 공간 해상도 20 pc에서 0.05 pc까지 500 Myr 동안 진행되었다. 주요 변수로는 (1) 광이온화 방사선 포함 여부, (2) 초신성(SN) 에너지(10⁵¹ erg 기본값 대비 변동), (3) 자기장(초기 B≈5 µG) 포함 여부, (4) 해상도, (5) 서브그리드 클러핑 모델(밀도 분포의 로그정규형태를 가정한 보정)이다.

해상도 의존성 결과는 10 M⊙→1 M⊙로 갈수록 R_mol이 평균 2배 증가함을 보여, 현재 1 M⊙ 수준에서도 완전 수렴이 이루어지지 않음을 시사한다. 이는 작은 스케일의 밀도 구조가 H₂ 형성에 결정적 역할을 함을 의미한다. 광이온화와 자기장을 각각 제외하면 R_mol이 약 2배씩 감소한다. 광이온화는 FUV 차폐를 감소시켜 H₂ 파괴를 가속화하고, 자기장은 가스 압축을 억제해 평균 밀도를 낮춘다.

가장 중요한 결과는 서브그리드 클러핑 모델을 적용했을 때 R_mol이 약 3배 상승한다는 점이다. 클러핑은 해상도 한계 이하의 고밀도 코어를 통계적으로 보정함으로써, 실제 ISM에서 기대되는 작은 스케일 클라우드의 존재를 반영한다. 이 모델을 포함한 경우, 시간 평균 R_mol≈0.25가 되며, 이는 태양권 관측값(0.19–0.26)과 PHANGS‑ALMA 샘플의 중위값(0.42) 사이에 위치한다.

또한 SN 에너지 변화를 시험했을 때, 에너지를 2배 증가시키면 R_mol이 약 30% 감소하고, 0.5배로 감소시키면 약 20% 증가한다. 이는 SN 피드백이 가스 압축·팽창을 통해 밀도 분포를 재조정하고, 결과적으로 H₂ 형성 효율에 직접적인 영향을 미친다.

전체적으로, 논문은 현재 가장 정교한 은하 패치 시뮬레이션조차도 미세한 밀도 구조와 광이온화·자기장 상호작용을 완전히 포착하지 못한다는 점을 강조한다. 서브그리드 클러핑과 같은 통계적 보정이 없으면, 시뮬레이션은 관측된 H₂ 함량을 지속적으로 과소평가한다는 것이 핵심 결론이다.