Computer simulations of the structural and physical properties of the olivine and spinel polymorphs of Mg2SiO4

The aim of the work presented is to develop a computer simulation technique which will predict the structure and physical properties of forsterite and ringwoodite, the major mantle-forming polymorphs of Mg₂SiO₄. The technique is based upon energy minimization, in which all structural parameters are varied until the configuration with the lowest energy is achieved. The lattice energy and physical properties (e.g. elasticity and dielectric constants) are calculated from interatomic potentials, which generally include electrostatic and short-range terms. We investigate several types of traditional potential models, and present a new type of model which includes partial ionic charges and a Morse potential to describe the effect of covalency on the Si-O bond. This new form of potential model is highly successful, and not only reproduces the zero-pressure structural, elastic and dielectric properties of forsterite and ringwoodite, but also accurately describes their pressure dependence.