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

Basic Information

• Paper ID: 2508.00791
• Title: Electrical and Thermal Conductivity of Earth's Iron-enriched Basal Magma Ocean
• Authors: Francis Dragulet, Lars Stixrude (UCLA)
• Classification: astro-ph.EP cond-mat.mtrl-sci physics.geo-ph
• Publication Date: August 17, 2025 (Draft version)
• Target Journal: Proceedings of the National Academy of Sciences
• Paper Link: https://arxiv.org/abs/2508.00791

Abstract

Earth's earliest magnetic field likely originated from the basal magma ocean (BMO)—a layer of silicate melt surrounding the core that may have persisted for billions of years. Recent research indicates that under BMO conditions, liquids with Earth's bulk silicate composition exhibit electrical conductivity exceeding 10⁴ S/m, potentially reaching the threshold for dynamo activity. However, due to iron's incompatibility in lower mantle mineral assemblages, the basal magma ocean was likely iron-enriched relative to bulk silicate Earth throughout most of its history. This study employs first-principles molecular dynamics calculations to examine how iron content affects the silicate dynamo hypothesis.

Research Background and Motivation

Scientific Questions

1. The Mystery of Early Earth's Magnetic Field: Earth's magnetic field has been active for at least 3.5 billion years, crucial for maintaining planetary habitability, yet the current dynamo driven by outer core liquid iron convection may have existed for less than 1 billion years
2. The Basal Magma Ocean Dynamo Hypothesis: The BMO as a possible source of early magnetic fields requires satisfying the condition of magnetic Reynolds number Rm > 40
3. Unknown Iron Enrichment Effects: Existing conductivity calculations primarily rely on simplified systems, neglecting iron enrichment effects during BMO crystallization

Research Significance

• Explains the mechanism by which Earth's early magnetic field protected the atmosphere and prevented atmospheric loss
• Provides theoretical foundation for magnetic field evolution in rocky exoplanets
• Deepens understanding of Earth's internal structure and thermal evolution

Limitations of Existing Methods

• Previous conductivity calculations focused on simplified systems such as MgO, MgSiO₃, and SiO₂
• Lack of systematic investigation of iron enrichment during BMO crystallization
• Insufficient exploration of coupling relationships between thermal and electrical conductivity

Core Contributions

1. First systematic study of iron enrichment effects on BMO electrical and thermal conductivity
2. Establishes quantitative relationships: Electrical conductivity exhibits quadratic dependence on iron content, consistent with Mott-Ziman theory
3. Discovers key mechanisms: Wide energy bands formed by hybridization of iron 3d electrons with oxygen 2p states are the fundamental cause of conductivity enhancement
4. Constructs evolutionary models: Predicts temporal evolution of BMO dynamo by combining thermal evolution models
5. Demonstrates feasibility: Iron enrichment significantly extends silicate dynamo lifetime to 3.3 billion years

Methodology Details

Task Definition

Calculate electrical conductivity σel and thermal conductivity kel of silicate liquids with varying iron content (XFe = Fe/(Fe+Mg) = 0.12, 0.5, 1.0) under BMO conditions (100-140 GPa, 4000-6000 K).

Computational Architecture

1. Molecular Dynamics Simulations

• Theoretical Foundation: Density Functional Theory (DFT) + PBEsol approximation + Hubbard U correction (U-J = 2.5 eV)
• Implementation Platform: VASP code, projector augmented-wave method
• Simulation Parameters:
  • System: 149 atoms, 6 elements (Mg, Fe, Si, O, Ca, Al)
  • Time step: 1 fs, total duration: 10-15 ps
  • Energy cutoff: 500 eV
  • K-point sampling: Gamma point
  • Ensemble: NVT, Nosé-Hoover thermostat

2. Electronic Transport Property Calculations

Employs the Kubo-Greenwood method in Chester-Thellung form:

Electrical Conductivity: $\sigma_{el}(\omega) = L_{11}(\omega)$

Thermal Conductivity: $k_{el}(\omega) = \frac{1}{e^2T}\left(L_{22}(\omega) - \frac{L_{12}(\omega)L_{21}(\omega)}{L_{11}(\omega)}\right)$

where the kinetic coefficients are: $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}$

3. Spin State Treatment

Accounts for high-spin to low-spin transition of iron through thermodynamic integration of free energy difference: $\Delta F_{HS-LS} = \int_0^1 \langle\Delta U\rangle_\lambda d\lambda$

Equilibrium electrical conductivity: $\frac{1}{\sigma_{el}^{eq}} = \frac{f_{eq}}{\sigma_{el}^{HS}} + \frac{1-f_{eq}}{\sigma_{el}^{LS}}$

Technical Innovations

1. Multi-scale Coupling: Combines atomic-scale DFT calculations with macroscopic thermal evolution models
2. Spin State Treatment: Accurately accounts for spin transition effects of iron at high pressure
3. Compositional Evolution: Temporal evolution of BMO composition based on elemental partition coefficients
4. Nonlinear Phase Diagram: Employs nonlinear liquidus temperature model, more realistically reflecting Fe-Mg system behavior

Experimental Setup

Computational Parameters

• Pressure Range: 100-140 GPa (BMO-relevant pressures)
• Temperature Range: 4000-6000 K (BMO-relevant temperatures)
• Iron Content: XFe = 0.12 (olivine), 0.5, 1.0 (iron-rich end member)
• Convergence Criteria: Energy 3 meV/atom, pressure 0.2 GPa

Evaluation Metrics

• Electrical Conductivity: σel (S/m), compared with dynamo threshold (~10⁴ S/m)
• Thermal Conductivity: kel (W/m/K), assesses convection stability
• Magnetic Reynolds Number: Rm = μ₀vlσ, dynamo criterion (Rm > 40)

Comparison Benchmarks

• Ionic thermal conductivity of pure MgSiO₃ system (~4-5 W/m/K)
• Electrical and thermal conductivity of liquid iron in Earth's core
• Wiedemann-Franz law predictions

Experimental Results

Main Results

1. Electrical Conductivity Variation with Iron Content

• XFe = 0.12: σel ~ 2×10⁴ S/m (6000K, 100GPa)
• XFe = 1.0: σel > 10⁵ S/m, approximately 5-fold increase
• Temperature Dependence: σel increases with temperature
• Pressure Dependence: Minimal effect within BMO pressure range

2. Thermal Conductivity Evolution

• XFe = 0.12: kel ~ 5 W/m/K
• XFe = 1.0: kel ~ 30 W/m/K
• Still Lower than Core: Significantly smaller than liquid iron outer core (>100 W/m/K)

3. Electronic Density of States Analysis

• Iron 3d electrons hybridize with oxygen 2p states forming wide energy bands
• Density of states at Fermi level g(EF) increases linearly with XFe
• σel ∝ g(EF)², consistent with Mott-Ziman theory

Magnetic Field Evolution Simulation

Temporal Evolution Results

• Constant Composition (XFe=0.12): Dynamo persists for 1.4 Gyr
• Considering Iron Enrichment: Dynamo persists for 3.3 Gyr (until 1.2 Ga ago)
• Critical Thickness: Minimum BMO thickness required to maintain dynamo decreases from 200 km to 100 km

Dynamo Regime Analysis

Via regime diagram determination:

• BMO consistently remains in the thermal convection dynamo regime
• Does not reach the regime where high thermal conductivity suppresses dynamo activity
• Compared to the core, BMO exhibits more favorable σ-k parameter combinations

Ablation Studies

• Spin State Effects: Low-spin state conductivity approximately twice that of high-spin state
• Velocity Scaling: CIA scaling produces 1.3 Gyr shorter dynamo lifetime than MLT scaling
• Ionic Contribution: Electronic conduction accounts for 70-90% of total conductivity

Related Work

Electrical Conductivity Calculations

• Holmström et al. (2018): Simple MgO, MgSiO₃ systems
• Scipioni et al. (2017): Earth's bulk silicate composition
• Soubiran & Militzer (2018): Silicate liquids at high pressure

BMO Evolution Models

• Labrosse et al. (2007): Classical BMO thermal evolution model
• Stixrude et al. (2020): Silicate dynamo theoretical framework
• Ziegler & Stegman (2013): BMO dynamo numerical simulations

Iron Spin Transitions

• Holmström & Stixrude (2016): Spin transitions in liquid oxides
• Shim et al. (2023): Shock compression experimental verification

Conclusions and Discussion

Main Conclusions

1. Iron Enrichment Significantly Enhances Conductivity: Electrical conductivity increases approximately one order of magnitude as XFe increases from 0.12 to 1.0
2. Extended Dynamo Lifetime: Considering iron enrichment, the silicate dynamo can persist for 3.3 billion years
3. Convection Remains Sustainable: Despite increased thermal conductivity, BMO maintains convective instability
4. Early Magnetic Field Feasibility: Provides a reasonable physical mechanism for Earth's early magnetic field

Limitations

1. Compositional Simplification: Does not account for SiO₂ content reduction during crystallization
2. Convection Patterns: Does not include contributions from compositional convection
3. Radiative Heat Transfer: Neglects radiative thermal conductivity contributions
4. Three-Dimensional Effects: Based on one-dimensional thermal evolution model

Future Directions

1. Super-Earth Applications: Investigate electrical conductivity properties at higher pressures
2. Compositional Convection: Quantify convection effects driven by iron settling
3. Multi-element Systems: Consider effects of Al, Ca, and other elements
4. Experimental Verification: High-temperature, high-pressure experimental measurements for validation

In-Depth Evaluation

Strengths

1. Advanced Methodology: Combines first-principles calculations with macroscopic models, multi-scale investigation
2. Clear Physics: Explains conduction mechanisms from electronic structure perspective
3. Important Results: Provides quantitative constraints on early Earth's magnetic field origin
4. Complete Model: Accounts for multiple physical processes including spin transitions and compositional evolution

Weaknesses

1. Computational Cost: DFT+U method parameter selection may affect results
2. Model Assumptions: One-dimensional thermal evolution model is overly simplified
3. Experimental Gaps: Lacks high-temperature, high-pressure experimental data for validation
4. Uncertainties: Certain parameters (e.g., partition coefficients) exhibit substantial uncertainty

Impact

1. Academic Value: Provides important theoretical foundation for geophysics and planetary science
2. Application Prospects: Extensible to exoplanet magnetic field research
3. Methodological Exemplar: Demonstrates research paradigm in computational geophysics
4. Policy Implications: Enhances understanding of Earth's habitability evolution

Applicable Scenarios

• Early Earth magnetic field evolution research
• Rocky exoplanet magnetic field prediction
• Silicate properties at high temperature and pressure
• Earth's internal structure and dynamics modeling

References

This paper cites important literature from geophysics, computational physics, and planetary science, including:

• Labrosse et al. (2007): Classical BMO thermal evolution model
• Christensen & Aubert (2006): Dynamo magnetic Reynolds number threshold
• Pozzo et al. (2012): Earth's core thermal conductivity calculations
• McDonough & Sun (1995): Earth's olivine composition model

This paper provides quantitative physical constraints on the major scientific question of early Earth's magnetic field origin through rigorous first-principles calculations and thermal evolution modeling, representing an important contribution to the field of computational geophysics.