Assessing foundational atomistic models for iron alloys under Earth's core conditions
Abstract
We assess the capability of recently developed foundational atomistic models (FAMs) to simulate iron alloys under the extreme pressures and temperatures of Earth's core. Static equations of state for hexagonal close-packed (hcp) and body-centered cubic (bcc) iron, computed using 17 FAMs, are benchmarked against ab initio calculations. Two representative models, MatterSim and MACE, are further evaluated for their ability to reproduce phonon spectra, liquid structure, and melting relations of iron at core conditions. While both models capture several key properties, MACE substantially overestimates the stability of bcc iron and fails to correctly describe the stability of hcp iron. Their performance is also examined for binary liquids, superionic phases, and a seven-component Fe-Ni-Si-S-O-H-C liquid. Although these FAMs were not explicitly trained on data from core conditions, they can reproduce several structural and dynamical properties across a wide range of compositions. However, none of the tested models consistently reproduces all first-principles benchmarks. By analyzing the origins of these discrepancies, we identify several limitations of current FAMs, particularly the lack of an explicit treatment of thermal electronic excitations, which significantly affect phase stability and thermodynamic properties under core conditions. We further discuss directions for improving FAMs to enable predictive simulations of core-forming materials under extreme conditions.
Related
- A domain-specific atomistic model for Earth’s core materials
- Early Hadean Earth conditions were ideal for prebiotic chemistry
- Earth observation and causal inference in assessing protected area effectiveness
- Van der Waals interactions in supercritical water under Earth's mantle conditions
- Artificial upwelling via Stommel pipes under present and future climate conditions: an Earth system model assessment
BibTeX
@misc{wan2026assessing,
title = {{Assessing foundational atomistic models for iron alloys under Earth's core conditions}},
author = {Tianqi Wan and Liangrui Wei and Zepeng Wu and Renata M. Wentzcovitch and Yang Sun},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2605.13594},
eprint = {2605.13594},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2605.13594},
}
RIS
TY - UNPB
TI - Assessing foundational atomistic models for iron alloys under Earth's core conditions
AU - Tianqi Wan
AU - Liangrui Wei
AU - Zepeng Wu
AU - Renata M. Wentzcovitch
AU - Yang Sun
PY - 2026
DA - 2026/09/23
PB - arXiv
DO - 10.48550/arxiv.2605.13594
UR - https://arxiv.org/abs/2605.13594
AB - We assess the capability of recently developed foundational atomistic models (FAMs) to simulate iron alloys under the extreme pressures and temperatures of Earth's core. Static equations of state for hexagonal close-packed (hcp) and body-centered cubic (bcc) iron, computed using 17 FAMs, are benchmarked against ab initio calculations. Two representative models, MatterSim and MACE, are further evaluated for their ability to reproduce phonon spectra, liquid structure, and melting relations of iron at core conditions. While both models capture several key properties, MACE substantially overestimates the stability of bcc iron and fails to correctly describe the stability of hcp iron. Their performance is also examined for binary liquids, superionic phases, and a seven-component Fe-Ni-Si-S-O-H-C liquid. Although these FAMs were not explicitly trained on data from core conditions, they can reproduce several structural and dynamical properties across a wide range of compositions. However, none of the tested models consistently reproduces all first-principles benchmarks. By analyzing the origins of these discrepancies, we identify several limitations of current FAMs, particularly the lack of an explicit treatment of thermal electronic excitations, which significantly affect phase stability and thermodynamic properties under core conditions. We further discuss directions for improving FAMs to enable predictive simulations of core-forming materials under extreme conditions.
ER -