|
|
 |
|||||||||||||||||
|
|||||||||||||||||
 |
|||||||||||||||||
| JOURNAL HOME | HELP | CONTACT PUBLISHER | SUBSCRIBE | ARCHIVE | SEARCH | TABLE OF CONTENTS |
1 Department of Earth and Planetary Sciences, Faculty of Science, Kyushu University, Hakozaki, Fukuoka 812-8581, Japan
2 School of Chemistry, University of Bath, Claverton Down, Bath BA2 7AY, U.K.
3 Daresbury Laboratory, Warrington WA4 4AD, Cheshire, U.K.
Correspondence: * E-mail: matsui{at}geo.kyushu-u.ac.jp
Molecular dynamics (MD) simulation is used to calculate the elastic constants and their temperature and pressure derivatives, and the T-P-V equation of state of MgO. The interionic potential is taken to be the sum of pairwise additive Coulomb, van der Waals attraction, and repulsive interactions. In addition, to account for the observed large Cauchy violation of the elastic constants of MgO, the breathing shell model (BSM) is introduced in MD simulation, in which the repulsive radii of O ions are allowed to deform isotropically under the effects of other ions in the crystal. Quantum correction to the MD pressure is made using the Wigner-Kirkwood expansion of the free energy. Required energy parameters, including oxygen breathing parameters, were derived empirically to reproduce the observed molar volume and elastic constants of MgO, and their measured temperature and pressure derivatives as accurately as possible. The MD simulation with BSM is found to be very successful in reproducing accurately the measured molar volumes and individual elastic constants of MgO over a wide temperature and pressure range. The errors in the simulated molar volumes are within 0.3% over the temperature range between 300 and 3000 K at 0 GPa, and within 0.1% over the pressure range from 0 up to 50 GPa at 300 K. The simulated bulk modulus is found to be correct to within 0.7% between 300 and 1800 K at 0 GPa. Here we present the MD simulated T-P-V equation of state of MgO as an accurate internal pressure calibration standard at high temperatures and high pressures.
This article has been cited by other articles:
|
|
 |
|
  M. MATSUI, T. KATSURA, A. KUWATA, K. HAGIYA, N. TOMIOKA, M. SUGITA, S. YOKOSHI, A. NOZAWA, and K.-I. FUNAKOSHI Equation of state of (Mg0.8,Fe0.2)2SiO4 ringwoodite from synchrotron X-ray diffraction up to 20 GPa and 1700 K European Journal of Mineralogy, October 1, 2006; 18(5): 523 - 528. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 E. Ohtani and K. D. Litasov The Effect of Water on Mantle Phase Transitions Reviews in Mineralogy and Geochemistry, January 1, 2006; 62(1): 397 - 420. [Full Text] [PDF]
|
|
|
|
 |
|
 D. Tinker, C. E. Lesher, G. M. Baxter, T. Uchida, and Y. Wang High-pressure viscometry of polymerized silicate melts and limitations of the Eyring equation American Mineralogist, November 1, 2004; 89(11-12): 1701 - 1708. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 T. Irifune Application of synchrotron radiation and Kawai-type apparatus to various studies in high-pressure mineral physics Mineralogical Magazine, October 1, 2002; 66(5): 769 - 790. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Shirasaka, E. Takahashi, Y. Nishihara, K. Matsukage, and T. Kikegawa In situ X-ray observation of the reaction dolomite = aragonite + magnesite at 900-1300 K American Mineralogist, July 1, 2002; 87(7): 922 - 930. [Abstract] [Full Text] [PDF]
|
|
| JOURNAL HOME | HELP | CONTACT PUBLISHER | SUBSCRIBE | ARCHIVE | SEARCH | TABLE OF CONTENTS |
