Ab-initio chemical potentials of solid and liquid solutions and the chemistry of the Earth's core

Abstract

A general set of methods is presented for calculating chemical potentials in solid and liquid mixtures using {\em ab initio} techniques based on density functional theory (DFT). The methods are designed to give an {\em ab initio} approach to treating chemical equilibrium between coexisting solid and liquid solutions, and particularly the partitioning ratio of solutes between such solutions. For the liquid phase, the methods are based on the general technique of thermodynamic integration, applied to calculate the change of free energy associated with the continuous interconversion of solvent and solute atoms, the required thermal averages being computed by DFT molecular dynamics simulation. For the solid phase, free energies and hence chemical potentials are obtained using DFT calculation of vibrational frequencies of systems containing substitutional solute atoms, with anharmonic contributions calculated, where needed, by thermodynamic integration. The practical use of the methods is illustrated by applying them to study chemical equilibrium between the outer liquid and inner solid parts of the Earth's core, modelled as solutions of S, Si and O in Fe. The calculations place strong constraints on the chemical composition of the core, and allow an estimate of the temperature at the inner-core/outer-core boundary.

19 pages, two figures

Cited by in corpus (10)

• The melting curve of MgO from first principles simulations
• Temperature of the inner-core boundary of the Earth: Melting of iron at high pressure from first-principles coexistence simulations
• First principles simulations of direct coexistence of solid and liquid aluminium
• On the existence and structure of a mush at the inner core boundary of the Earth
• The axial ratio of hcp iron at the conditions of the Earth's inner core
• Thermal convection in Earth's inner core with phase change at its boundary
• The melting curve of iron at extreme pressures: implications for planetary cores
• Implication of the lopsided growth for the viscosity of Earth's inner core
• Ab initio superionic-liquid phase diagram of Fe1-xOx under Earth's inner core conditions
• Prediction of crystal structures and motifs in the Fe-Mg-O system at Earth's core pressures