Abstract:Molecular dynamics simulations for liquid CaAl2Si2O8 have been carried out at 72 state points spanning ranges in density (ρ: 2398–4327kg/m3), temperature (T: 3490–6100K) and pressure (P: 0.84–120GPa) relevant to geosystems. The atomic scale structure of the melt is determined by analysis of nearest neighbor coordination statistics as a function of T and P. Dramatic structural change occurs as pressure increases especially for 0< P <20GPa at all temperatures. Changes in structure are encapsulated by examining the coordination of Si, Al, Ca and O around oxygen and vice versa. Si and Al change from predominantly fourfold at low-P to dominantly sixfold for P >∼20GPa. Pentahedrally coordinated Si and Al in distorted trigonal bipyramids attain abundance maxima corresponding to ∼60% of total (Si, Al)O n at 3–5GPa and weakly depend on T. The coordination of Ca by oxygen increases from 7 to 10 for 0< P <20GPa and changes slowly for P >20GPa at 3500K. Similar behavior is seen at 6000K except that the interval of rapid changes occurs at higher pressure. Oxygen with only one nearest Si or Al neighbor (i.e., non-bridging oxygen, NBO) decreases whereas oxygen with two or three nearest neighbors of Si, Al or Ca increases as pressure increases. Changes in melt structure are reflected in the variation of thermodynamic and transport properties of the liquid. Values of the self-diffusivities of Ca, Al, Si and O are fit to a modified Arrhenian expression and compare well to limited laboratory data. Self-diffusivities are best fit using ‘low P’ and ‘high-P’ expressions, identical in form but with different parameters, with activation energies and activation volumes in the range 150–200kJ/mol and +5 to −1cm3/mol, respectively. Green–Kubo calculations for liquid shear viscosity are presented and compare well with limited laboratory results. Application of the Eyring model to determine the characteristic size and number of atoms in the activated cluster based on independently computed D and η suggests that the activated cluster decreases from ∼8 to ∼3 atoms from low to high pressure while its characteristic size shrinks from ∼14Å to ∼3Å providing insight into dynamics of atom mobility and possible cooperative behavior. The equation of state and variation of internal energy with T and V are used in Part II (Ghiorso et al., 2009) to derive a comprehensive thermodynamic description of liquid CaAl2Si2O8. This is best accomplished by allowing for EOS expressions broken into high and low pressure intervals consistent with coordination statistics and MD-derived transport properties.

Critical Vaporization of MgSiO3.

EXPERIMENTAL STUDY ON THE PHASE RELATIONS OF THE SYSTEM Mg2SiO4-MgAl2O4 AT HIGH PRESSURE AND HIGH TEMPERATURE

Probing the critical point of MgSiO3 using deep potential simulation

First principles study of density, viscosity, and diffusion coefficients of liquid MgSiO3 at conditions of the Earth's deep mantle

The Crucial Role of Transient Tri-Coordinated Oxygen in the Flow of Silicate Melts.

Structures and Transport Properties of Supercritical SiO2-H2O and NaAlSi3O8-H2O Fluids

Shock Response and Phase Transitions of MgO at Planetary Impact Conditions

VapoRock: Thermodynamics of Vaporized Silicate Melts for Modeling Volcanic Outgassing and Magma Ocean Atmospheres

Properties of high-temperature phase diagram and critical point parameters in silica

Room temperature Raman spectroscopy and 29Si MAS NMR combined with high temperature Raman spectroscopy and DFT calculation of xMgO-(1-x)CaO–SiO2 glasses and melts

Phase diagram and density of SiO2–H2O fluid across critical conditions

Structural properties of molten silicates from ab initio molecular-dynamics simulations: comparison between CaO-Al$_2$O$_3$-SiO$_2$ and SiO$_2$

X-ray Raman Scattering Study of MgSiO3 Glass at High Pressure: Implication for Triclustered MgSiO3 Melt in Earth's Mantle

Structure, thermodynamic and transport properties of CaAl2Si2O8 liquid. Part I: Molecular dynamics simulations

In Situ Observations of Liquid–liquid Phase Separation in Aqueous MgSO4 Solutions: Geological and Geochemical Implications

Thermodynamic and transport properties of meteor melt constituents from ab initio simulations: MgSiO 3, SiO 2, and MgO

Thermal and pressure ionization in warm, dense MgSiO3 studied with first-principles computer simulations

Decaying shock studies of phase transitions in MgOSiO2 systems: implications for the Super-Earths interiors

Molecular Dynamics Simulation on Liquid Nucleation in Ideal Silicon Crystal upon Melting

Raman Spectroscopic Observations of the Ion Association between Mg(2+) and SO4(2-) in MgSO4-Saturated Droplets at Temperatures of ≤380 °C.

Silica solubility and molecular speciation in water vapor at 400–800 °C