Thermodynamics of hexagonal-close-packed iron under Earth’s core conditions

The free energy and other thermodynamic properties of hexagonal-close-packed iron are calculated by direct ab initio methods over a wide range of pressures and temperatures relevant to the Earth’s core. The ab initio calculations are based on density-functional theory in the generalized-gradient approximation, and are performed using the projector augmented wave 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üneisen parameter, and comparisons are made with values obtained from diamond-anvil-cell and shock experiments

The earth’s core: an approach from first principles

The Earth’s core is largely composed of iron (Fe), alloyed with less dense elements such as sulphur, silicon and/or oxygen. The phase relations and physical properties of both solid and liquid Fe-alloys are therefore of great geophysical importance. As a result, over the past fifty years the properties of Fe and its alloys have been extensively studied experimentally. However, achieving the extreme pressures (up to 360 GPa) and temperatures (~6000K) found in the core provide a major experimental challenge, and it is not surprising that there are still considerable discrepancies in the results obtained by using different experimental techniques. In the past fifteen years quantum mechanical techniques have been applied to predict the properties of Fe. Here we review the progress that has been made in the use of first principles methods to study Fe and its alloys, and as a result of these studies we conclude: (i) that pure Fe adopts an hexagonal close packed structure under core conditions and melts at ~6200 K at 360 GPa, (ii) that thermodynamic equilibrium and observed seismic data are satisfied by a liquid Fe alloy outer core with a composition of ~10 mole% S (or Si) and 8 mole% O crystallising at ~ 5500 K to give an Fe alloy inner core with ~8 mole% S (or Si) and 0.2 mole % O, and (iii) that with such concentrations of S (or Si), an Fe alloy might adopt a body centred cubic structure in all or part of the inner core. In the future the roles of Ni, C, H and K in the core need to be studied, and techniques to predict the transport and rheological properties of Fe alloys need to be developed

Ab-initio simulation of high-temperature liquid selenium

Ab initio molecular dynamics simulation is used to investigate the structure and dynamics of liquid Se at temperatures of 870 and 1370~K. The calculated static structure factor is in excellent agreement with experimental data. The calculated radial distribution function gives a mean coordination number close to 2, but we find a significant fraction of one-fold and three-fold atoms, particularly at 1370~K, so that the chain structure is considerably disrupted. The self-diffusion coefficient has values ($\sim 1 \times 10^{-8}$~m~s$^{-1}$) typical of liquid metals.Comment: 10 pages, 4 Poscript figures, uses REVTE

The Ab-Initio Simulation of the Liquid Ga-Se System

Ab-initio dynamical simulation is used to study the liquid Ga-Se system at the three concentrations Ga$_2$Se, GaSe and Ga$_2$Se$_3$ at the temperature 1300~K. The simulations are based on the density functional pseudopotential technique, with the system maintained on the Born-Oppenheimer surface by conjugate gradients minimization. We present results for the partial structure factors and radial distribution functions, which reveal how the liquid structure depends on the composition. Our calculations of the electrical conductivity $\sigma$ using the Kubo-Greenwood approximation show that $\sigma$ depends very strongly on the composition. We show how this variation of $\sigma$ is related to the calculated electronic density of states. Comparisons with recent experimental determinations of the structure and conductivity are also presented.Comment: REVTEX, 8 pages, 4 uuencoded poscript figures, ([email protected]

Stability of the iterative solutions of integral equations as one phase freezing criterion

A recently proposed connection between the threshold for the stability of the iterative solution of integral equations for the pair correlation functions of a classical fluid and the structural instability of the corresponding real fluid is carefully analyzed. Direct calculation of the Lyapunov exponent of the standard iterative solution of HNC and PY integral equations for the 1D hard rods fluid shows the same behavior observed in 3D systems. Since no phase transition is allowed in such 1D system, our analysis shows that the proposed one phase criterion, at least in this case, fails. We argue that the observed proximity between the numerical and the structural instability in 3D originates from the enhanced structure present in the fluid but, in view of the arbitrary dependence on the iteration scheme, it seems uneasy to relate the numerical stability analysis to a robust one-phase criterion for predicting a thermodynamic phase transition.Comment: 11 pages, 2 figure

The particle-in-cell model for ab initio thermodynamics: implications for the elastic anisotropy of the Earth's inner core

We assess the quantitative accuracy of the particle-in-cell (PIC) approximation used in recent ab initio predictions of the thermodynamic properties of hexagonal-close-packed iron at the conditions of the Earth's inner core. The assessment is made by comparing PIC predictions for a range of thermodynamic properties with the results of more exact calculations that avoid the PIC approximation. It is shown that PIC gives very accurate results for some properties, but that it gives an incorrect treatment of anharmonic lattice vibrations. In addition, our assessment does not support recent PIC-based predictions that the hexagonal c/a ratio increases strongly with increasing temperature, and we point out that this casts doubt on a proposed re-interpretation of the elastic anisotropy of the inner core.Comment: 25 pages, 9 figures, submitted to Physics of the Earth and Planetary Interior

A Self-Consistent First-Principles Technique Having Linear Scaling

An algorithm for first-principles electronic structure calculations having a computational cost which scales linearly with the system size is presented. Our method exploits the real-space localization of the density matrix, and in this respect it is related to the technique of Li, Nunes and Vanderbilt. The density matrix is expressed in terms of localized support functions, and a matrix of variational parameters, L, having a finite spatial range. The total energy is minimized with respect to both the support functions and the elements of the L matrix. The method is variational, and becomes exact as the ranges of the support functions and the L matrix are increased. We have tested the method on crystalline silicon systems containing up to 216 atoms, and we discuss some of these results.Comment: 12 pages, REVTeX, 2 figure

Towards a Linear-Scaling DFT Technique: The Density Matrix Approach

A recently proposed linear-scaling scheme for density-functional pseudopotential calculations is described in detail. The method is based on a formulation of density functional theory in which the ground state energy is determined by minimization with respect to the density matrix, subject to the condition that the eigenvalues of the latter lie in the range [0,1]. Linear-scaling behavior is achieved by requiring that the density matrix should vanish when the separation of its arguments exceeds a chosen cutoff. The limitation on the eigenvalue range is imposed by the method of Li, Nunes and Vanderbilt. The scheme is implemented by calculating all terms in the energy on a uniform real-space grid, and minimization is performed using the conjugate-gradient method. Tests on a 512-atom Si system show that the total energy converges rapidly as the range of the density matrix is increased. A discussion of the relation between the present method and other linear-scaling methods is given, and some problems that still require solution are indicated.Comment: REVTeX file, 27 pages with 4 uuencoded postscript figure

A first principles study of sub-monolayer Ge on Si(001)

Experimental observations of heteroepitaxial growth of Ge on Si(001) show a (2xn) reconstruction for sub-monolayer coverages, with dimer rows crossed by missing-dimer trenches. We present first-principles density-functional calculations designed to elucidate the energetics and relaxed geometries associated with this reconstruction. We also address the problem of how the formation energies of reconstructions having different stoichiometries should be compared. The calculations reveal a strong dependence of the formation energy of the missing-dimer trenches on spacing n, and demonstrate that this dependence stems almost entirely from elastic relaxation. The results provide a natural explanation for the experimentally observed spacings in the region of n \~ 8.Comment: 13 pages, 4 figures, submitted to Surface Scienc

The Adsorption of H2O on TiO2 and SnO2(110) Studied by First-Principles Calculations

First-principles calculations based on density functional theory and the pseudopotential method have been used to investigate the energetics of H$_2$O adsorption on the (110) surface of TiO$_2$ and SnO$_2$. Full relaxation of all atomic positions is performed on slab systems with periodic boundary conditions, and the cases of full and half coverage are studied. Both molecular and dissociative (H$_2$O $\rightarrow$ OH$^-$ + H$^+$) adsorption are treated, and allowance is made for relaxation of the adsorbed species to unsymmetrical configurations. It is found that for both TiO$_2$ and SnO$_2$ an unsymmetrical dissociated configuration is the most stable. The symmetrical molecularly adsorbed configuration is unstable with respect to lowering of symmetry, and is separated from the fully dissociated configuration by at most a very small energy barrier. The calculated dissociative adsorption energies for TiO$_2$ and SnO$_2$ are in reasonable agreement with the results of thermal desorption experiments. Calculated total and local electronic densities of states for dissociatively and molecularly adsorbed configurations are presented and their relation with experimental UPS spectra is discussed