Ab initio melting temperatures of bcc and hcp iron under the Earth's inner core condition

Abstract

There has been a long debate on the stable phase of iron under the Earth's inner core conditions. Because of the solid-liquid coexistence at the inner core boundary, the thermodynamic stability of solid phases directly relates to their melting temperatures, which remains considerable uncertainty. In the present study, we utilized a semi-empirical potential fitted to high-temperature ab initio data to perform a thermodynamic integration from classical systems described by this potential to ab initio systems. This method provides a smooth path for thermodynamic integration and significantly reduces the uncertainty caused by the finite-size effect. Our results suggest the hcp phase is the stable phase of pure iron under the inner core conditions, while the free energy difference between the hcp and bcc phases is tiny, on the order of 10s meV/atom near the melting temperature.

4 figures

References in corpus (3)

• Temperature of the inner-core boundary of the Earth: Melting of iron at high pressure from first-principles coexistence simulations
• Two-step nucleation of the Earth's inner core
• Automated free energy calculation from atomistic simulations

Cited by in corpus (9)

• Searching for iron nanoparticles with a general-purpose Gaussian approximation potential
• Deep-learning-based prediction of the tetragonal$\to$cubic transition in davemaoite
• Melting Temperature of Iron Under the Earth's Inner Core Condition from Deep Machine Learning
• The Fe-Ni phase diagram and the Earth's inner core structure
• Ab initio superionic-liquid phase diagram of Fe1-xOx under Earth's inner core conditions
• Understanding the two-step nucleation of iron at Earth's inner core conditions: a comparative molecular dynamics study
• Thermal Induced Structural Competitiveness and Metastability of Body-centered Cubic Iron under Non-Equilibrium Conditions
• Structure and Melting of Fe, MgO, SiO2, and MgSiO3 in Planets: Database, Inversion, and Phase Diagram
• Incorporating Gibbs free energy into interatomic potential fitting