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.

The Partial Molar Volume of Fe2O3 in Alkali Silicate Melts: Evidence for an Average Fe3+ Coordination Number Near Five

The density and compressibility of CaO–FeO–SiO 2 liquids at one bar: Evidence for four-coordinated Fe 2+ in the CaFeO 2 component

The Partial Molar Volume and Thermal Expansivity of TiO2 in Alkali Silicate Melts: Systematic Variation with Ti Coordination

Density and Sound Speed Measurements on Model Basalt (an–di–hd) Liquids at One Bar: New Constraints on the Partial Molar Volume and Compressibility of the FeO Component

The structure of CaO–MgO–Al2O3–SiO2 melts and glasses doped with FeOX–NiO

Multi-Technique Equation of State for Fe2sio4 Melt and the Density of Fe-Bearing Silicate Melts from 0 to 161 Gpa

Thermal Expansion of Metal Oxides in Silicate Melts: A Simple Theoretical Analysis

Revision 1 1 ! Recoil-free Fractions of Iron in Aluminous Bridgmanite from Temperature 2 ! Dependent Mössbauer Spectra 3 ! 4 !

Solution behavior of C-O-H volatiles in FeO-Na2O-Al2O3-SiO2 melts in equilibrium with liquid iron alloy and graphite at 4 GPa and 1550°C

Effect of Coordination Structure of Iron Ions on Iron Oxide Activity Coefficients in the CaO‐SiO2‐FeO‐Fe2O3 Slags at 1573 K Under Oxygen Partial Pressures Between 10−9 and 10−4 atm

Influence of oxygen fugacity on the solubility of carbon and hydrogen in FeO-Na2O-SiO2-Al2O3 melts in equilibrium with liquid iron at 1.5 GPa and 1400°C

Acoustic velocity measurements on Na2O–TiO2–SiO2 liquids: Evidence for a highly compressible TiO2 component related to five-coordinated Ti

Origin of Fe 3+ in Fe-containing, Al-free Mantle Silicate Perovskite

Estimation of the proportions of aluminum-fluorine species in NaF-AlF 3 mixtures by in situ Raman spectroscopy and theoretical simulations

Density, viscosity and surface tension of high-silicate CaO–SiO2 and CaO–SiO2–Fe2O3 slags derived by aerodynamic levitation. The behavior of Fe3+ in high-silicate melts

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

Crystal Structure of Calcium-Ferrite Type NaAlSiO4 Up to 45 GPa

Quantitative studies on the microstructure of ternary CaO–Al2O3–SiO2 glasses by Raman spectroscopy, 27Al MAS NMR and quantum chemistry ab initio calculation

Micro-structure studies of the molten binary K3AlF6–Al2O3 system by in situ high temperature Raman spectroscopy and theoretical simulation

Temperature dependence of magnetoelastic properties of Fe 100 − xSix

Recoil-free Fractions of Iron in Aluminous Bridgmanite from Temperature-Dependent Mössbauer Spectra