프로디모: 가스·먼지·화학이 동시에 움직이는 원시행성계 원반 모델

초록

ProDiMo 코드를 이용해 방사선, 화학, 열·냉각, 수압 평형을 2차원적으로 결합한 원시행성계 원반 모델을 제시한다. T Tauri 별 주변 원반을 0.5 AU–500 AU 범위로 계산한 결과, 내부 가장자리(puffed‑up rim)가 부풀어 오르고, 0.1 < z/r < 0.5, 0.5–10 AU 구역에 온도 5000 K, 밀도 10⁷–10⁸ cm⁻³ 수준의 얇은 원자 가스 “헤일로”가 형성된다. 이 층은 별의 UV와 PAH가 열을 공급하고, Fe II와 O I 630 nm 라인으로 냉각한다. 먼 외곽(>100 AU)에서는 가스 온도가 먼지 온도와 분리돼 Tg≈2Td가 된다. 가스 열평형을 고려하지 않은 모델은 원반이 현저히 평평해짐을 보여, 가스 온도와 구조의 상호작용이 핵심임을 강조한다.

상세 분석

본 논문은 원시행성계 원반의 물리·화학 구조를 일관되게 해석하기 위한 새로운 수치 코드 ProDiMo(프로디모)를 소개한다. 기존 모델들은 대개 먼지 연속복사와 가스 화학을 별도로 다루거나, 가스 온도를 먼지 온도와 동일시하는 가정을 사용했다. ProDiMo는 2차원(반구면) 방사선 전달을 주파수 의존적으로 계산하고, UV와 X‑ray에 의한 광화학 반응망을 포함한다. 특히 Fe II와 CO의 로-진동 전이(rovibrational) 라인을 열·냉각 메커니즘에 추가함으로써, 별에 가까운 고밀도 영역(ρ ≈ 10¹⁰ cm⁻³)에서도 정확한 에너지 균형을 잡는다.

핵심은 가스와 먼지의 온도(Tg, Td)를 독립적으로 계산하고, 이를 수압 평형 방정식에 피드백한다는 점이다. 온도 차이가 큰 경우(예: 외곽 원반)에는 가스 압력이 증가해 원반의 부피가 팽창하고, 이는 다시 방사선 전파와 화학 반응에 영향을 미친다. 논문은 이러한 상호작용을 반복적인 이터레이션으로 수렴시켜, ‘열·화학·구조’가 서로 얽힌 자기일관적인 해를 얻는다.

모델 적용 사례는 T Tauri 별 주변 원반이며, 반경 0.5 AU에서 500 AU까지 확장한다. 결과는 다음과 같다. 첫째, 내측 가장자리(‘puffed‑up rim’)가 별빛을 직접 받으며 부풀어 오르고, 이는 내부 그림자 효과를 완화해 원반 전체에 더 많은 UV가 도달하게 만든다. 둘째, z/r ≈ 0.1–0.5, r ≈ 0.5–10 AU 구역에 온도 5 000 K, 밀도 10⁷–10⁸ cm⁻³ 수준의 얇은 원자 가스 층이 형성된다. 이 ‘halo’는 PAH와 작은 먼지에 의한 광전효과와 포톤 흡수에 의해 가열되고, Fe II 반감소선(세미포비든)과 O I 630 nm 전이로 효율적으로 냉각한다. 특히 Fe II 라인은 고밀도(>10⁸ cm⁻³)에서 비LTE 탈포화가 일어나며, 열 손실을 지배한다는 점이 새롭게 제시된다.

셋째, 먼 외곽(>100 AU)에서는 가스와 먼지 온도가 크게 분리된다. 방사선 차폐가 약해 별 UV가 직접 가스에 흡수되면서 Tg≈2Td까지 상승한다. 이때도 가스 압력 상승이 원반의 플레어(flare) 구조를 강화시켜, 원반이 더 넓게 퍼지는 형태를 만든다.

마지막으로, 가스 온도를 Td와 동일하게 가정한 ‘단순 모델’과 비교했을 때, 후자는 원반이 현저히 평평해져, 관측되는 IR·mm 파장 스펙트럼과 크게 불일치한다. 따라서 가스 열평형을 정확히 계산하는 것이 원반의 수직 구조와 방사선 전파, 나아가 관측 해석에 필수적임을 강력히 증명한다.

이러한 결과는 원반 내에서 행성 형성 초기 단계의 물질 이동, 화학적 풍부도, 그리고 가스·먼지 상호작용을 이해하는 데 중요한 기반을 제공한다. 특히, 고온 원자 ‘halo’는 원반 표면에서의 원소 이온화와 라인 방출을 설명하는 데 핵심적인 역할을 하며, 향후 ALMA·JWST 관측과 직접 연결될 수 있다.