지구 외핵 액체 실리콘‑산소‑철 혼합물의 전도성 및 열전도성 연구

초록

본 논문은 지구 외핵의 대표 조성인 Fe₀.₈₂Si₀.₁₀O₀.₀₈과 Fe₀.₇₉Si₀.₀₈O₀.₁₃ 액체 혼합물에 대해 밀도범함수이론(DFT)과 Kubo‑Greenwood 공식으로 열전도도와 전기전도도를 계산하였다. 외핵 상부에서는 열전도도 k≈100 W m⁻¹ K⁻¹, 전기전도도 σ≈1.1 × 10⁶ Ω⁻¹ m⁻¹이며, 하부에서는 각각 k≈160 W m⁻¹ K⁻¹, σ≈1.3 × 10⁶ Ω⁻¹ m⁻¹을 얻었다. 이는 기존 추정치보다 2‑3배 높으며, 핵의 급속 냉각 혹은 방사성 원소 함량이 높음을 시사한다. 또한 부분 라디얼 분포함수, 평균제곱변위, 점도, 음속 등 구조·동역학적 특성도 제시한다.

상세 분석

이 연구는 지구 외핵의 물리‑화학적 상태를 보다 정밀하게 규정하기 위해 두 가지 실리콘‑산소‑철 합금을 선택하였다. Fe₀.₈₂Si₀.₁₀O₀.₀₈은 외핵 상부 압·온도(≈135 GPa, 4000 K)에, Fe₀.₇₉Si₀.₀₈O₀.₁₃은 하부(≈330 GPa, 6000 K)에 대응하도록 설계되었다. 전자구조는 퍼포먼스가 검증된 PBE‑GGA 교환‑상관 함수와 PAW 포텐셜을 사용해 200 fs 정도의 마이크로캐노니컬 MD 시뮬레이션으로 샘플링하였다. Kubo‑Greenwood 공식은 전자 전이 행렬 요소를 직접 계산함으로써 전기전도도 σ와 열전도도 k를 동시에 얻는 데 활용되었으며, 윌슨-라우어-스미스(WS) 규칙을 적용해 라우어-라우어 상호작용을 보정하였다. 결과적으로 상부에서는 σ≈1.1 × 10⁶ Ω⁻¹ m⁻¹, k≈100 W m⁻¹ K⁻¹, 하부에서는 σ≈1.3 × 10⁶ Ω⁻¹ m⁻¹, k≈160 W m⁻¹ K⁻¹를 도출했으며, 이는 기존 실험·모델 추정치(σ≈0.4‑0.6 × 10⁶, k≈30‑50 W m⁻¹ K⁻¹)보다 2‑3배 높은 값이다. 이러한 차이는 전자‑이온 상호작용 강도가 기존 모델에서 과소평가되었음을 의미한다. 구조적 분석에서는 부분 라디얼 분포함수(g_{Fe‑Fe}, g_{Fe‑Si}, g_{Fe‑O} 등)를 통해 Fe‑Fe 네트워크가 고밀도 액체 구조를 유지하면서도 Si와 O가 부분적으로 1차 근접 껍질에 삽입되는 모습을 확인했다. 평균제곱변위(MSD)와 그에 기반한 확산계수는 10⁻⁹ m² s⁻¹ 수준으로, 점도는 약 0.5‑0.8 mPa·s 정도이며, 이는 외핵이 매우 낮은 점성을 갖는 유동성 액체임을 뒷받침한다. 음속은 10‑12 km s⁻¹ 범위로, 지진학적 모델과 일치한다. 이러한 물리량들은 열전도도가 높아질 경우 외핵의 열손실이 가속화되어 지구 전체의 냉각 속도가 기존 추정보다 1.5‑2배 빨라질 수 있음을 시사한다. 또한, 높은 전도도는 지구 자기장의 지속에 필요한 전류 밀도를 감소시키며, 방사성 원소(예: ⁴⁰K, ⁵⁶U)의 함량이 기존보다 크게 존재할 가능성을 열어준다. 연구는 계산 불확실성을 5‑10 % 수준으로 추정했으며, 향후 고압·고온 실험과 비교해 모델을 정교화할 필요성을 강조한다.