마그마 해양과 금속성에 따른 실리콘‑메탄 경쟁: TOI‑421b 대기 진단

초록

이 논문은 실리콘수소화물(SiH₄)과 메탄(CH₄)의 경쟁을 통해 서브넵튠의 대기 금속성 및 내부 마그마 해양 존재 여부를 진단한다. 고압·고온 마그마‑대기 경계에서 실제 기체 상태방정식을 적용해 화학 평형을 계산하고, 용해도와 실리콘 함량이 온도·압력·용융률에 따라 어떻게 변하는지 보여준다. 결과는 낮은 금속성(1×태양)에서는 SiH₄가 우세하고, 높은 금속성(100×태양)에서는 CH₄가 지배한다는 점을 강조한다. 마그마 용융률이 클수록 SiH₄/CH₄ 비율이 높아지며, 이는 관측 가능한 대기 스펙트럼에서 금속성 및 마그마 해양 존재 여부를 추정하는 새로운 지표가 된다.

상세 분석

본 연구는 서브넵튠의 마그마‑대기 경계(MEB, >10 kbar)에서 기체가 이상기체 가정을 벗어나 실제 기체 거동을 보이는 점에 주목한다. 저자들은 H‑He‑C‑N‑O‑Si 6원소 시스템을 대상으로, 최신 실험 데이터와 고압·고온 충격 실험 결과를 이용해 SiO와 SiH₄의 상태방정식을 직접 피팅하였다. SiO는 Redlich‑Kwong EOS와 대응상태법을 결합해 Z₍c₎≈0.33을 재현했고, SiH₄는 44–138 GPa, 1200–4100 K 범위의 부피‑압력 데이터를 2차 Virial EOS( Z = 1 + A(T)P + B(T)P² )에 비선형 최소제곱법으로 적합해 A₀≈3.86 GPa⁻¹, B₀≈‑0.019 GPa⁻² 등 6개의 파라미터를 도출하였다. 이는 기존에 CH₄와 동일한 압축성을 가정했던 선행 연구와 달리, 실리콘 수소화물의 비이상성을 정량적으로 반영한다는 점에서 의미가 크다.

모델 설정은 TOI‑421b( Rₚ=2.64 R⊕, Mₚ=6.7 M⊕, T_eq≈920 K)를 기준으로, MEB 온도 3000 K, 압력 10–30 kbar를 가정하고, 다양한 H 질량분율(0.1–1 wt %)과 금속성(1×, 100× 태양) 및 마그마 용융률(1 %–100 %)을 탐색한다. 원소 질량 예산은 Lodders(2009)와 Bulk Silicate Earth(45.4 wt % SiO₂) 값을 기반으로, 금속성에 따라 He, C, N, Si, O의 비중을 조정하였다.

화학 평형 계산 결과, H와 N이 마그마에 가장 높은 용해도를 보이며 He와 C가 그 뒤를 잇는다. 완전 용융된 마그마(100 % melt)에서는 Si가 풍부해 SiH₄가 MEB에서 지배적인 실리콘‑수소 화합물이 된다. 그러나 금속성이 100× 태양으로 증가하면 C와 N의 절대량이 급증해 CH₄와 NH₃ 등 탄소·질소 화합물이 우세해지며, SiH₄는 상대적으로 억제된다. 마그마 용융률이 감소하면 (10 %·1 %) Si와 O의 총량이 크게 줄어들어 SiH₄의 기체 압력이 감소하고, 동시에 H와 He의 마그마 용해도가 낮아져 대기 상부에 H₂·He‑풍부한 조성이 남는다.

관측 가능한 대기(1 mbar–100 bar)로 확장했을 때, Si‑함유 응결 클라우드(예: SiO₂, SiC) 형성으로 SiH₄와 SiO의 기체 농도가 크게 감소한다. 그럼에도 불구하고, 낮은 금속성·고용융률 경우 SiH₄는 여전히 검출 가능한 수준을 유지한다. 결과적으로 SiH₄/CH₄ 비율과 Si/C 원소 비율은 (i) 마그마 용융률이 클수록 상승, (ii) 금속성이 높을수록 감소하는 경향을 보이며, 이는 JWST와 같은 전이천문학 관측에서 서브넵튠의 내부 상태를 역추정하는 새로운 진단 도구가 된다.

이 논문의 주요 기여는: (1) 고압·고온에서 SiH₄의 비이상성을 정량화한 EOS 제공, (2) 마그마‑대기 상호작용을 포함한 전체 원소 질량 보존 모델 구축, (3) SiH₄와 CH₄의 경쟁을 금속성·마그마 용융률의 함수로 제시함으로써 관측 가능 스펙트럼에 직접 연결한 점이다. 향후 연구에서는 다른 실리콘‑산소 화합물(SiO, Si₂O₃ 등)과 광학적 클라우드 모델을 결합해, 전파·적외선 전이선의 상세한 전이 확률을 계산하는 것이 필요하다.