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

The crucial role of transient tri-coordinated oxygen in the flow of silicate melts

The transport properties of high-temperature silicate melts control magma flow and are crucial for a wide variety of industrial processes involving minerals. However, anomalous melt properties have been observed that cannot be explained by the traditional polymerization degree theory, which was derived based on quenched melts. Ab initio molecular dynamics (AIMD) simulations were conducted to investigate the flow mechanism of CaO-Al2O3-SiO2 melts under high temperature atmospheric conditions. By analyzing the dynamic structure of melted silicates and employing molecular orbital theory, we gained a fundamental understanding of the flow mechanism from a chemistry perspective. Transient tri-coordinated oxygen (TO) bonded with one Si and two Al atoms (SiOAl2) was found to be a pivotal intermediate in melt flow and atomic diffusion processes. Frequent chemical transition between TO in SiOAl2 and bridging oxygen (BO) dominated the fluidity of melted silicates. The presence of such transitions is facilitated by the unstable nature of [SiAlO2] 4-membered rings, which are susceptible to instability due to the intense repulsion between the O 2p lone pairs and the excessively bent O-Al-O angle. Additionally, the density of SiOAl2 type TO motif could serve as an indicator to determine the relationship between structure and fluidity. Our results challenge the traditional polymerization degree theory and suggest the need to reassess high-temperature liquid properties that govern processes in the Earth and industry by monitoring transient motifs.

The structure of molten calcium ferrite under various redox conditions

Laser-heated melts based on the 43CaO-57Fe2O3-x eutectic, close to the calcium ferrite (CF) composition, were measured with high-energy X-ray diffraction using aerodynamic levitation over a range of redox states controlled by CO/CO2 gas atmospheres. The iron-oxygen coordination number was found to rise from 4.4 ± 0.3 at 15% Fe3+ to 5.3 ± 0.3 at 87% Fe3+. Empirical potential structure refinement modelling was used to obtain the ferric and ferrous partial pair distribution functions. It was found that the Fe2+ iron-oxygen coordination number is consistently approximately 10% higher in CF than in pure iron oxide, while Fe3+ is essentially identical in all but the most oxygen-rich environments (where it is higher in CF compared with FeOx). The model also shows calcium octahedra to be the dominant species across all redox environments, although the population of CaO7 increases with the availability of oxygen at the expense of CaO4 and CaO5. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

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.

Buoyancy and Structure of Volatile-Rich Silicate Melts

The early Earth was marked by at least one global magma ocean. Melt buoyancy played a major role for its evolution. Here we model the composition of the magma ocean using a six-component pyrolite melt, to which we add volatiles in the form of carbon as molecular CO or CO2 and hydrogen as molecular H2O or through substitution for magnesium. We compute the density relations from first-principles molecular dynamics simulations. We find that the addition of volatiles renders all the melts more buoyant compared to the reference volatile-free pyrolite. The effect is pressure dependent, largest at the surface, decreasing to about 20 GPa, and remaining roughly constant to 135 GPa. The increased buoyancy would have enhanced convection and turbulence, and thus promoted the chemical exchanges of the magma ocean with the early atmosphere. We determine the partial molar volume of both H2O and CO2 throughout the magma ocean conditions. We find a pronounced dependence with temperature at low pressures, whereas at megabar pressures the partial molar volumes are independent of temperature. At all pressures, the polymerization of the silicate melt is strongly affected by the amount of oxygen added to the system while being only weakly affected by the specific type of volatile added.

The Critical Point and the Supercritical State of Alkali Feldspars: Implications for the Behavior of the Crust During Impacts

The position of the vapor-liquid dome and of the critical point determine the evolution of the outermost parts of the protolunar disk during cooling and condensation after the Giant Impact. The parts of the disk in supercritical or liquid state evolve as a single thermodynamic phase; when the thermal trajectory of the disk reaches the liquid-vapor dome, gas and melt separate leading to heterogeneous convection and phase separation due to friction. Different layers of the proto-Earth behaved differently during the Giant Impact depending on their constituent materials and initial thermodynamic conditions. Here we use first-principles molecular dynamics to determine the position of the critical point for NaAlSi3O8 and KAlSi3O8 feldspars, major minerals of the Earth and Moon crusts. The variations of the pressure calculated at various volumes along isotherms yield the position of the critical points: 0.5-0.8 g cm-3 and 5500-6000 K range for the Na-feldspar, 0.5-0.9 g cm-3 and 5000-5500 K range for the K-feldspar. The simulations suggest that the vaporization is incongruent, with a degassing of O2 starting at 4000 K and gas component made mostly of free Na and K cations, O2, SiO and SiO2 species for densities below 1.5 g cm-3. The Hugoniot equations of state imply that low-velocity impactors (<8.3 km s-1) would at most melt a cold feldspathic crust, whereas large impacts in molten crust would see temperatures raise up to 30000 K.