The Speciation and Coordination of a Deep Earth Carbonate-Silicate-Metal Melt.

Ab initio molecular dynamics calculations on a carbonate-silicate-metal melt were performed to study speciation and coordination changes as a function of pressure and temperature. We examine in detail the bond abundances of specific element pairs and the distribution of coordination environments over conditions spanning Earth's present-day mantle. Average coordination numbers increase continuously from 4 to 8 for Fe and Mg, from 4 to 6 for Si, and from 2 to 4 for C from 1 to 148 GPa (4,000 K). Speciation across all pressure and temperature conditions is complex due to the unusual bonding of carbon. With the increasing pressure, C-C and C-Fe bonding increase significantly, resulting in the formation of carbon polymers, C-Fe clusters, and the loss of carbonate groups. The increased bonding of carbon with elements other than oxygen indicates that carbon begins to replace oxygen as an anion in the melt network. We evaluate our results in the context of diamond formation and of metal-silicate partitioning behavior of carbon. Our work has implications for properties of carbon and metal-bearing silicate melts, such as viscosity, electrical conductivity, and reactivity with surrounding phases.

Pressure-Induced Coordination Changes in a Pyrolitic Silicate Melt From Ab Initio Molecular Dynamics Simulations.

With ab initio molecular dynamics simulations on a Na-, Ca-, Fe-, Mg-, and Al-bearing silicate melt of pyrolite composition, we examine the detailed changes in elemental coordination as a function of pressure and temperature. We consider the average coordination as well as the proportion and distribution of coordination environments at pressures and temperatures encompassing the conditions at which molten silicates may exist in present-day Earth and those of the Early Earth's magma ocean. At ambient pressure and 2,000 K, we find that the average coordination of cations with respect to oxygen is 4.0 for Si-O, 4.0 for Al-O, 3.7 for Fe-O, 4.6 for Mg-O, 5.9 for Na-O, and 6.2 for Ca-O. Although the coordination for iron with respect to oxygen may be underestimated, the coordination number for all other cations are consistent with experiments. By 15 GPa (2,000 K), the average coordination for Si-O remains at 4.0 but increases to 4.1 for Al-O, 4.2 for Fe-O, 4.9 for Mg-O, 8.0 for Na-O, and 6.8 for Ca-O. The coordination environment for Na-O remains approximately constant up to core-mantle boundary conditions (135 GPa and 4000 K) but increases to about 6 for Si-O, 6.5 for Al-O, 6.5 for Fe-O, 8 for Mg-O, and 9.5 for Ca-O. We discuss our results in the context of the metal-silicate partitioning behavior of siderophile elements and the viscosity changes of silicate melts at upper mantle conditions. Our results have implications for melt properties, such as viscosity, transport coefficients, thermal conductivities, and electrical conductivities, and will help interpret experimental results on silicate glasses.

Anharmonic contribution to the stabilization of Mg(OH)2 from first principles.

Geometrical and vibrational characterization of magnesium hydroxide was performed using density functional theory. Four possible crystal symmetries were explored: P3[combining macron] (No. 147, point group -3), C2/m (No. 12, point group 2), P3m1 (No. 156, point group 3m) and P3[combining macron]m1 (No. 164, point group -3m) which are the currently accepted geometries found in the literature. While a lot of work has been performed on Mg(OH)2, in particular for the P3[combining macron]m1 phase, there is still a debate on the observed ground state crystal structure and the anharmonic effects of the OH vibrations on the stabilization of the crystal structure. In particular, the stable positions of hydrogen are not yet defined precisely, which have implications in the crystal symmetry, the vibrational excitations, and the thermal stability. Previous work has assigned the P3[combining macron]m1 polymorph as the low energy phase, but it has also proposed that hydrogens are disordered and they could move from their symmetric position in the P3[combining macron]m1 structure towards P3[combining macron]. In this paper, we examine the stability of the proposed phases by using different descriptors. We compare the XRD patterns with reported experimental results, and a fair agreement is found. While harmonic vibrational analysis shows that most phases have imaginary modes at 0 K, anharmonic vibrational analysis indicates that at room temperature only the C2/m phase is stabilized, whereas at higher temperatures, other phases become thermally competitive.

Proton dynamics and the phase diagram of dense water ice.

All the different phases of water ice between 2 GPa and several megabars are based on a single body-centered cubic sub-lattice of oxygen atoms. They differ only by the behavior of the hydrogen atoms. In this study, we investigate the dynamics of the H atoms at high pressures and temperatures in water ice from first-principles molecular dynamics simulations. We provide a detailed analysis of the O-H⋯O bonding dynamics over the entire stability domain of the body-centered cubic (bcc) water ices and compute transport properties and vibrational density-of-states. We report the first ab initio evidence for a plastic phase of water and we propose a coherent phase diagram for bcc water ices compatible with the two groups of melting curves and with the multiple anomalies reported in ice VII around 15 GPa.

Pressure-dependent isotopic composition of iron alloys.

Our current understanding of Earth's core formation is limited by the fact that this profound event is far removed from us physically and temporally. The composition of the iron metal in the core was a result of the conditions of its formation, which has important implications for our planet's geochemical evolution and physical history. We present experimental and theoretical evidence for the effect of pressure on iron isotopic composition, which we found to vary according to the alloy tested (FeO, FeH(x), or Fe3C versus pure Fe). These results suggest that hydrogen or carbon is not the major light-element component in the core. The pressure dependence of iron isotopic composition provides an independent constraint on Earth's core composition.

Ab initio determination of the valence electron distribution in the average structure of the incommensurately modulated calaverite AuTe2.

The valence-electron density distribution of the average structure of incommensurately modulated calaverite, AuTe2, has been computed using density-functional theory. High-density regions, centered around the Au and Te atoms, are not spheric, but present charge concentrations along the Au-Te and Te-Te bonds. The electronic band structure and its corresponding density of states reveal the presence of three electronic band groups, constituted mainly by Te 5s, Au 5d and hybrids of Te 6p + Au 6s + Au 5d orbitals. The electrons belonging to the last block are responsible for the chemical bonds.