2025-11-21T14:19:15.500377

Electrical and Thermal Conductivity of Earth's Iron-enriched Basal Magma Ocean

Dragulet, Stixrude

The Earth's earliest magnetic field may have originated in a basal magma ocean, a layer of silicate melt surround the core that could have persisted for billions of years. Recent studies show that the electrical conductivity of liquid with a bulk silicate Earth composition exceeds 10000 S/m at basal magma ocean conditions, potentially surprising the threshold for dynamo activity. Over most of its history however, the basal magma ocean is more enriched in iron than the bulk silicate Earth, due to iron's incompatibility in the mineral assemblages of the lower mantle. Using ab-initio molecular dynamics calculations, we examine how iron content affects the silicate dynamo hypothesis. We investigate how the electrical conductivity of silicate liquid changes with iron enrichment, at pressures and temperatures relevant for Earth's basal magma ocean. We also compute the electronic contribution to the thermal conductivity , to evaluate convective instability of basal magma oceans. Finally, we apply our results to model the thermal and magnetic evolution of Earth's basal magma ocean over time.

academic

地球の鉄富化基底マグマオーシャンの電気伝導率と熱伝導率

基本情報

論文ID: 2508.00791
タイトル: Electrical and Thermal Conductivity of Earth's Iron-enriched Basal Magma Ocean
著者: Francis Dragulet, Lars Stixrude (UCLA)
分類: astro-ph.EP cond-mat.mtrl-sci physics.geo-ph
発表日時: 2025年8月17日(ドラフト版)
投稿先: Proceedings of the National Academy of Sciences
論文リンク: https://arxiv.org/abs/2508.00791

要旨

地球の初期磁場は、地核を取り囲むケイ酸塩融体の層である基底マグマオーシャン(BMO)に由来する可能性がある。BMOは数十億年間存在していた可能性がある。最新の研究によると、基底マグマオーシャンの条件下では、地球全体のケイ酸塩組成を持つ液体の電気伝導率は10⁴ S/mを超え、発電機活動の閾値に達する可能性がある。しかし、下部マントル鉱物組合における鉄の非相溶性のため、基底マグマオーシャンはその歴史の大部分において、地球全体のケイ酸塩よりも鉄に富んでいた。本研究は、第一原理分子動力学計算を用いて、鉄含有量がケイ酸塩発電機仮説にどのように影響するかを検証する。

研究背景と動機

科学的問題

地球初期磁場起源の謎：地球磁場は少なくとも35億年間活動しており、地球の居住可能性の維持に不可欠であるが、現在の外核液体鉄対流によって駆動される発電機は10億年未満の歴史しか持たない可能性がある
基底マグマオーシャン発電機仮説：BMOが初期磁場の可能な起源として、磁気レイノルズ数Rm > 40の条件を満たす必要がある
鉄富化効果の未知性：既存の電気伝導率計算は主に簡略化されたシステムに基づいており、BMO結晶化過程における鉄富化効果を考慮していない

研究の重要性

地球初期磁場が大気圏を保護し、大気流失を防ぐメカニズムの説明
岩石質系外惑星の磁場進化に対する理論的基礎の提供
地球内部構造と熱進化の理解の深化

既存手法の限界

従来の電気伝導率計算はMgO、MgSiO₃、SiO₂などの簡略化されたシステムに集中している
BMO結晶化過程における鉄富化の体系的研究が不足している
熱伝導率と電気伝導率の結合関係が十分に検討されていない

核心的貢献

初めての体系的研究：鉄富化がBMOの電気伝導率と熱伝導率に与える影響
定量的関係の確立：電気伝導率と鉄含有量は二次関係を示し、Mott-Ziman理論と一致
重要メカニズムの発見：鉄3d電子とO 2p状態のハイブリッド化により形成される広いエネルギーバンドが電気伝導率増強の根本原因
進化モデルの構築：熱進化モデルと組み合わせてBMO発電機の時間進化を予測
実現可能性の証明：鉄富化によってケイ酸塩発電機の寿命が33億年まで大幅に延長される

方法の詳細

タスク定義

異なる鉄含有量(XFe = Fe/(Fe+Mg) = 0.12, 0.5, 1.0)のケイ酸塩液体のBMO条件下(100-140 GPa, 4000-6000 K)における電気伝導率σelと熱伝導率kelを計算する。

計算方法アーキテクチャ

1. 分子動力学シミュレーション

理論的基礎：密度汎関数理論(DFT) + PBEsol近似 + Hubbard U補正(U-J = 2.5 eV)
実装プラットフォーム：VASPコード、投影増強平面波法
シミュレーションパラメータ：
- システム：149原子、6種類の元素(Mg, Fe, Si, O, Ca, Al)
- タイムステップ：1 fs、総時間：10-15 ps
- エネルギーカットオフ：500 eV
- K点サンプリング：ガンマ点
- アンサンブル：NVT、Nosé-Hoover温度制御

2. 電子輸送特性の計算

Chester-Thellung形式のKubo-Greenwood法を採用：

電気伝導率： $\sigma_{el}(\omega) = L_{11}(\omega)$

熱伝導率： $k_{el}(\omega) = \frac{1}{e^2T}\left(L_{22}(\omega) - \frac{L_{12}(\omega)L_{21}(\omega)}{L_{11}(\omega)}\right)$

ここで動力学係数： $L_{\alpha\beta}(\omega) = \frac{2\pi e^2\hbar^2}{3m_e^2\omega\Omega}\sum_{i,j}(f_i-f_j)\delta(\epsilon_i-\epsilon_j-\hbar\omega)|\langle\psi_i|\nabla|\psi_j\rangle|^2(-1)^{\alpha+\beta}(\epsilon_i-\mu)^{\alpha-1}(\epsilon_j-\mu)^{\beta-1}$