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 Mg2SiO4. 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.

[1]  Y. Takéuchi,et al.  Electron-density distributions of three orthopyroxenes, Mg2Si2O6, Co2Si2O6, and Fe2Si2O6 , 1982 .

[2]  M. Catti.  The lattice energy of forsterite. Charge distribution and formation enthalpy of the SiO44− ion , 1981 .

[3]  A. B. Lidiard,et al.  Point Defects in Ionic Crystals , 1972 .

[4]  K. Fujino,et al.  On the estimation of atomic charges by the X-ray method for some oxides and silicates , 1980 .

[5]  G. R. Barsch,et al.  Elastic constants of single‐crystal forsterite as a function of temperature and pressure , 1969 .

[6]  M. Catti.  Atomic charges in Mg2SiO4 (forsterite), fitted to thermoelastic and structural properties , 1982 .

[7]  William R. Busing,et al.  WMIN: A computer program to model molecules and crystals in terms of potential energy functions , 1981 .

[8]  L. Pauling.  The Nature Of The Chemical Bond , 1939 .

[9]  John C. Slater,et al.  Quantum Theory of Molecules and Solids , 1951 .

[10]  Seishiro Mizukami,et al.  High-pressure X-ray diffraction studies on - and ?-Mg2SiO4 , 1975 .

[11]  Robert M. Hazen,et al.  Effects of temperature and pressure on the crystal structure of forsterite , 1976 .

[12]  Gerald V. Gibbs,et al.  Molecules as models for bonding in silicates , 1982 .

[13]  Hiroyasu Shimizu,et al.  Pressure Dependence of Elastic Moduli of Forsterite by Brillouin Scattering in the Diamond Cell , 1982 .

[14]  H. Takeda,et al.  Atomic diffusion coefficients calculated for transition metals in olivine , 1983, Nature.

[15]  C. Prewitt,et al.  Size and Compressibility of Ions at High Pressure , 1982 .

[16]  H. Takeda,et al.  An interpretation of the structures of mantle minerals at high pressure in terms of interionic forces , 1980 .

[17]  W. D. Kingery,et al.  Introduction to Ceramics , 1976 .

[18]  G. V. Gibbs,et al.  Molecular orbital studies of geometries and spectra of minerals and inorganic compounds , 1977 .

[19]  C. Prewitt,et al.  Single‐crystal X ray study of γ Mg2SiO4 , 1982 .

[20]  Y. Ida.  Interionic repulsive force and compressibility of ions , 1976 .

[21]  G. Herzberg,et al.  Molecular Spectra and Molecular Structure: I. Spectra of Diatomic Molecules , 1944 .

[22]  R. Hazen.  Comparative crystal chemistry , 1982 .

[23]  R. Hazen,et al.  The crystal structures of forsterite and hortonolite at several temperatures up to 900°C , 1973 .

[24]  T. Gilbert,et al.  Soft‐Sphere Model for Closed‐Shell Atoms and Ions , 1968 .

[25]  G. V. Gibbs,et al.  A molecular orbital study of shared-edge distortions in linked polyhedra , 1976 .

[26]  J. A. Mandarino.  The Gladstone-Dale relationship; Part I, Derivation of new constants , 1976 .

[27]  A. Das,et al.  Preparation and Characterization of Forsterite and Measurement of its Dielectric Constant and Loss Factor in the Frequency Range of 100 Kc/S to 25 Mc/S , 1979 .

[28]  J. M. Thomas,et al.  Simulating silicate structures and the structural chemistry of pyroxenoids , 1982, Nature.

[29]  K. Fujino,et al.  The ionic compressibilities and radii estimated for some transition metals in olivine structures , 1982 .

[30]  Stephen C. Parker,et al.  Prediction of mineral crystal structures , 1983 .

[31]  H. Sawamoto,et al.  ORTHORHOMBIC PEROVSKITE (Mg, Fe)SiO3 AND CONSTITUTION OF THE LOWER MANTLE , 1977 .

[32]  R. D. Shannon,et al.  Effective ionic radii in oxides and fluorides , 1969 .

[33]  Mineo Kumazawa,et al.  Elastic moduli, pressure derivatives, and temperature derivatives of single‐crystal olivine and single‐crystal forsterite , 1969 .

[34]  M. Born,et al.  Dynamical Theory of Crystal Lattices , 1954 .

[35]  P. A. Reynolds,et al.  X-ray determination of electron distributions in forsterite, fayalite and tephroite , 1981 .

[36]  C. Catlow,et al.  Point defect and electronic properties of uranium dioxide , 1977, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[37]  J. Tossell.  A comparison of silicon-oxygen bonding in quartz and magnesian olivine from X-ray spectra and molecular orbital calculations , 1977 .