Search

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., Comment: 19 pages, two figures

{\em Ab initio} techniques based on density functional theory in the projector-augmented-wave implementation are used to calculate the free energy and a range of other thermodynamic properties of liquid iron at high pressures and temperatures relevant to the Earth's core. The {\em ab initio} free energy is obtained by using thermodynamic integration to calculate the change of free energy on going from a simple reference system to the {\em ab initio} system, with thermal averages computed by {\em ab initio} molecular dynamics simulation. The reference system consists of the inverse-power pair-potential model used in previous work. The liquid-state free energy is combined with the free energy of hexagonal close packed Fe calculated earlier using identical {\em ab initio} techniques to obtain the melting curve and volume and entropy of melting. Comparisons of the calculated melting properties with experimental measurement and with other recent {\em ab initio} predictions are presented. Experiment-theory comparisons are also presented for the pressures at which the solid and liquid Hugoniot curves cross the melting line, and the sound speed and Gr\"{u}neisen parameter along the Hugoniot. Additional comparisons are made with a commonly used equation of state for high-pressure/high-temperature Fe based on experimental data., Comment: 16 pages including 6 figures and 5 tables

We present a new technique which allows the fully {\em ab initio} calculation of the chemical potential of a substitutional impurity in a high-temperature crystal, including harmonic and anharmonic lattice vibrations. The technique uses the combination of thermodynamic integration and reference models developed recently for the {\em ab initio} calculation of the free energy of liquids and anharmonic solids. We apply the technique to the case of the substitutional oxygen impurity in h.c.p. iron under Earth's core conditions, which earlier static {\em ab initio} calculations indicated to be thermodynamically very unstable. Our results show that entropic effects arising from the large vibrational amplitude of the oxygen impurity give a major reduction of the oxygen chemical potential, so that oxygen dissolved in h.c.p. iron may be stabilised at concentrations up a few mol % under core conditions.

The free energy and other thermodynamic properties of hexagonal-close-packed iron are calculated by direct {\em ab initio} methods over a wide range of pressures and temperatures relevant to the Earth's core. The {\em ab initio} calculations are based on density-functional theory in the generalised-gradient approximation, and are performed using the projector augmented wave (PAW) approach. Thermal excitation of electrons is fully included. The Helmholtz free energy consists of three parts, associated with the rigid perfect lattice, harmonic lattice vibrations, and anharmonic contributions, and the technical problems of calculating these parts to high precision are investigated. The harmonic part is obtained by computing the phonon frequencies over the entire Brillouin zone, and by summation of the free-energy contributions associated with the phonon modes. The anharmonic part is computed by the technique of thermodynamic integration using carefully designed reference systems. Detailed results are presented for the pressure, specific heat, bulk modulus, expansion coefficient and Gr\"{u}neisen parameter, and comparisons are made with values obtained from diamond-anvil-cell and shock experiments., Comment: 18 pages, 14 figures

First principles electronic structure calculations based on density functional theory have been used to study the thermodynamic, structural and transport properties of solid solutions and liquid alloys of iron and oxygen at Earth's core conditions. Aims of the work are to determine the oxygen concentration needed to account for the inferred density in the outer core, to probe the stability of the liquid against phase separation, to interpret the bonding in the liquid, and to find out whether the viscosity differs significantly from that of pure liquid iron at the same conditions. It is shown that the required concentration of oxygen is in the region 25-30 mol percent, and evidence is presented for phase stability at these conditions. The Fe-O bonding is partly ionic, but with a strong covalent component. The viscosity is lower than that of pure liquid iron at Earth's core conditions. It is shown that earlier first-principles calculations indicating very large enthalpies of formation of solid solutions may need reinterpretation, since the assumed crystal structures are not the most stable at the oxygen concentration of interest., Comment: 21 pages, 12 figures