지구 내핵에서 초이온성 수소의 방사형 농도 구배

초록

본 연구는 고압·고온 조건에서 Fe‑H 시스템의 초이온성 고체와 액체 상의 자유에너지를 직접 계산하고, 이를 이용해 내핵 전역에 걸친 초이온성‑액체 평형을 구축한다. 온도를 순수 철의 융점으로 정규화하면 압력에 대한 의존성이 거의 사라지는 스케일링 관계를 발견했으며, 이 관계를 통해 기존 이론 간의 분배계수 차이를 해소한다. 열역학적 평형 조건을 적용하면 내핵 중심에서 외부 경계(ICB)까지 수소 함량이 약 7 at.% 감소하는 방사형 농도 구배가 자연스럽게 발생함을 보인다.

상세 분석

이 논문은 지구 내핵 깊이(323–360 GPa, 5000–6800 K)에서 Fe‑H 합금의 초이온성(hcp, bcc) 및 액체 상에 대한 ab initio Gibbs 자유에너지를 직접 계산한다. 먼저 AIMD 시뮬레이션을 통해 Fe‑H 결합 특성을 확인하고, Fe‑H 전용 머신러닝 포텐셜을 구축해 대규모 초이온‑액체 공존 시뮬레이션과 고정‑점 자유에너지 적분을 수행하였다. 결과적으로 bcc 초이온상은 hcp보다 자유에너지가 높아 메타안정 상태이며, 실제 내핵에서는 hcp 초이온상이 유일하게 액체와 공존한다는 점을 확인했다.

압력에 따른 상도표를 그린 뒤, 온도를 순수 Fe의 융점 (T_m^{Fe}(P))으로 나누는 축소 온도 (\Theta = T/T_m^{Fe})를 도입하였다. (\Theta)를 기준으로 하면 323 GPa와 360 GPa에서의 고체‑액체 경계선이 거의 겹치며, 압력은 단순히 에너지 스케일을 재조정할 뿐 수소 혼합열역학에는 큰 영향을 미치지 않음을 보여준다. 이 스케일링은 기존 연구들(예: Yuan & Steinle‑Neumann, Liu & Jing, Zhang 등)에서 보고된 서로 다른 분배계수 (K_{H}^{s/l}) 차이를 통합적으로 설명한다.

열역학적 평형 조건 (\mu_H^{solid}(P,T,x_s)=\mu_H^{liquid}(P,T,x_l))과 내핵이 거의 등온(≈5500 K)이라는 가정 하에, 자유에너지 데이터베이스를 이용해 각 깊이에서 수소 화학퍼텐셜을 계산한다. 화학퍼텐셜이 방사형으로 일정해야 하므로, 압력(깊이) 변화에 따라 수소 함량 (x)가 조정된다. 구체적으로 ICB(≈330 GPa)에서는 약 16 at.%의 수소가 용해되고, 중심부(≈360 GPa)에서는 약 12 at.%로 감소한다. 온도 변동(5000 K–6000 K)이나 중력 퍼텐셜을 포함한 보정에도 구배 형태는 유지되며, 중력 효과는 화학퍼텐셜 변화에 비해 미미하다.

마지막으로, 계산된 분배계수 (K_{H}^{s/l})는 압력·온도에 따라 0.57–0.70 사이이며, 이는 Fe‑H 시스템에서 수소가 액체보다 고체에 약 60 % 정도 더 많이 농축된다는 의미다. 기존 실험(1900 K, 45 GPa 등)과도 정량적으로 일치한다. 전체적으로, 초이온성 수소의 존재와 그 열역학적 거동이 내핵 내부의 방사형 화학조성을 자연스럽게 설명한다는 점이 가장 큰 과학적 기여이다.