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.

New thermodynamic regularity for cesium over the whole liquid range

Characterizing Phase Transitions in Liquid Cesium by a Soft-core and Large Attractive Equation of State

Linear Exp6 Isotherm for Compressed Molten Cesium Over the Whole Liquid range Including Metal-nonmetls Transition and Tc

Common compression factor and bulk modulus quiescence points of liquid cesium

Lattice Dynamics and Thermal Conductivity of Cesium Chloride Via First-Principles Investigation

Thermodynamics of cesium complexes formation with 18-crown-6 in ionic liquids

Possible phase transition in liquid Cesium at ambient pressure

Thermodynamics of complexation of cesium ions with dibenzo-18-crown-6 in hydrophobic ionic liquids

Generalized linear isotherm regularity equation of state applied to metals

Theoretical and practical investigation of ion–ion association in electrolyte solutions

Thermodynamic Properties and Equations of State for Solid and Liquid Magnesium

A fresh look at the vibrational and thermodynamic properties of liquids within the soft potential model

Stability constants of cesium complexes with 18-crown-6 in ionic liquids

New Equations of State describing both the Dynamic Viscosity and Self-Diffusion Coefficient for Potassium and Thallium in their fluid phases

Molecular dynamics simulation of potassium along the liquid-vapor coexistence curve

The thermodynamic properties of liquid 3He-4He mixtures between 0 and 20 atm in the limit of absolute zero temperature

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

Prediction of Viscosity of Liquid Alkali Metals Using an Improved Mean Spherical Approximation

Modeling of the thermodynamic properties of aqueous ionic liquid solutions with an equation of state for square-well chain fluid with variable range

A new correlation for thermal conductivity of liquids