Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published July 15, 2015 | Accepted Version
Journal Article Open

A self-consistent optimization of multicomponent solution properties: ab initio molecular dynamic simulations and the MgO-SiO_2 miscibility gap under pressure

Abstract

We propose a new approach to parameterizing the Gibbs energy of a multicomponent solution as a function of temperature, pressure and composition. It uses the quasichemical model in the second nearest neighbour approximation and considers both a polynomial representation (for low pressure) and an exponential decay representation (for moderate-to-high pressure) of the excess molar volume v^(xs) to extend thermodynamic behaviour to elevated pressure. This approach differs from previous configuration-independent regular or associated solution-type models of multicomponent silicate liquids at elevated pressure and can account for any structural or short-range order data that may be available. A simultaneous least squares fit of the molar volume and the molar enthalpy of mixing data obtained from First Principles Molecular Dynamics (FPMD) simulations at various pressures enables complete parameterization of the excess thermodynamic properties of the solution. Together with consistently optimized properties of coexisting solids, this enables prediction of pressure-temperature-composition phase diagrams associated with melting. Although the method is extensible to natural multicomponent systems, we apply the procedure as a first test case to the important planetary model system MgO-SiO_2 using FPMD data found in the literature. One key result of this optimization, which depends only on the derived excess properties of the liquid phase, is that the consolute temperature of the SiO_2-rich miscibility gap is predicted to decrease with increasing pressure. This appears to be in disagreement with available experimental constraints and suggests possible thermodynamic inconsistency between FPMD data and experimental phase equilibrium data in the 0-5 GPa pressure range. We propose a new thermodynamic consistency criterion relating the signs of v^(xs) and other excess properties and discuss the need for precise calculations of derivatives of excess properties. Finally, the potential reappearance of the miscibility gap in the MgO-SiO_2 system above 5 GPa is discussed in light of this work.

Additional Information

© 2015 Elsevier B.V. Received 15 August 2014, Accepted 2 April 2015, Available online 13 April 2015. We would like to thank Dr. James Sangster and Dr. Pierre Hudon for their constructive criticisms of the present work as well as the editors and reviewers at GCA for pushing us to improve the manuscript. This work was supported in part by Defence Research and Development Canada and by the US NSF through award EAR-1426526.

Attached Files

Accepted Version - 1-s2.0-S0016703715001970-main.pdf

Files

1-s2.0-S0016703715001970-main.pdf
Files (1.0 MB)
Name Size Download all
md5:70880a589f89528f18ea6c1924770cbe
1.0 MB Preview Download

Additional details

Created:
August 22, 2023
Modified:
October 23, 2023