An equation of state for silicate melts. III. Analysis of stoichiometric liquids at elevated pressure: Shock compression data, molecular dynamics simulations and mineral fusion curves

Experimental data and molecular dynamics simulations that constrain the densities of liquids of "mineral-like" stoichiometry at elevated pressure are evaluated using the liquid Equation of State (EOS) of Ghiorso (2004a). Phase equilibrium constraints on melt density are developed from experimental brackets on mineral fusion curves, specifically congruent melting of albite, anorthite, cristobalite/quartz/stishovite, diopside, enstatite, fayalite, forsterite, jadeite, nepheline, pyrope, sanidine, and titanite. A self consistent thermodynamic analysis of shock compression is applied to experimentally determined brackets on the Hugoniot for the liquid compositions CaMgSi2O6, CaAl2Si2O8, Fe2SiO4, and (CaMgSi2O6)(0.64)(CaAl2Si2O8)(0.36). Molecular dynamics (MD) simulations of both melt density and melt structure are analyzed utilizing the liquid EOS of Ghiorso (2004a) and an ideal associated solution model that accounts explicitly for the effect of pressure and temperature on the oxygen coordination environment of silicon and aluminum. MD data on Mg2Si2O6, CaMgSi2O6, SiO2, Na2Si4O9, CaAl2Si2O8, NaAlSiO4, NaAlSi2O6, and NaAlSi3O8 liquids are considered in conjunction with the fusion curve and shock compression experimental data sets. An internally consistent assessment of liquid EOS parameters from all available data sources is attempted for each liquid composition considered. Shock compression experiments on more chemically complex liquids of komatiite and MORB bulk composition are also examined. In support of the analysis of mineral fusion curves, volumetric properties as a function of temperature and pressure of diopside, enstatite (both Pbca and C2/c), fayalite, forsterite, and pyrope are evaluated from literature data. Experimental observations are parameterized to a Universal EOS (Vinet and others, 1986, 1987, 1989). It is found that the reference pressure properties derived by Ghiorso and Kress (2004) are internally consistent with the majority of high-pressure constraints on melt density. The notable exceptions are molten SiO2 and melts in the system MgO-SiO2. For the latter system, the calibration of Ghiorso and Kress (2004) fails to recover reference pressure volumes. The effect of configurational collapse (the volumetric response associated with changes in melt structure, such as changing oxygen coordination number about a cation) is systematically assessed from molecular dynamics simulation data. Fully depolymerized melts exhibit an similar to11 percent density increase associated with Si and Al transforming from IV- to V-fold coordination with respect to oxygen. Alkali aluminosilicate melts demonstrate similar to27 percent increase in melt density and SiO2 and Na2Si4O9 liquids show a density increase of similar to17 percent.