Please cite this article as: B'ejina, F., Blanchard, M., Wright, K., Price, G.D., A computer simulation study of the effect of pressure on Mg diffusion in forsterite, Physics of the Earth and Planetary Interiors (2007), doi:10.1016/j.pepi.2008 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. A c c e p t e d M a n u s c r i p t
AbstractComputer simulation techniques were used to investigate the effect of pressure on magnesium diffusion in forsterite between 0 and 10 GPa. We studied the diffusion path along the c crystallographic axis (we always refer to the Pbnm system) via a vacancy mechanism. Using a Mott-Littleton approach within the code GULP, we were able to precisely map the diffusion path of a Mg vacancy and we found the activation energy, E = 3.97 eV at 0 GPa (with E f = 3.35 eV for the formation energy and E m = 0.62 eV for the migration) and E = 4.46 eV at 10 GPa (E f = 3.81 eV and E m = 0.65 eV). Preliminary results using the supercell technique gave the same saddle point coordinates and energies. This saddle point of the Mg vacancy diffusion found with GULP was then introduced in an ab initio code, confirming the values of the migration energy both at 0 and 10 GPa. We were therefore able to estimate the activation volume (∆V ) to be around 5 cm 3 /mol and d(∆V )/dP ≃ 0.The effect of pressure applies mostly on defect formation and little on migration. A c c e p t e d M a n u s c r i p t
AbstractComputer simulation techniques were used to investigate the effect of pressure on magnesium diffusion in forsterite between 0 and 10 GPa. We studied the diffusion path along the c crystallographic axis (we always refer to the Pbnm system) via a vacancy mechanism. Using a Mott-Littleton approach within the code GULP, we were able to precisely map the diffusion path of a Mg vacancy and we found the activation energy, E = 3.97 eV at 0 GPa (with E f = 3.35 eV for the formation energy and E m = 0.62 eV for the migration) and E = 4.46 eV at 10 GPa (E f = 3.81 eV and E m = 0.65 eV). Preliminary results using the supercell technique gave the same saddle point coordinates and energies. This saddle point of the Mg vacancy diffusion found with GULP was then introduced in an ab initio code, confirming the values of the migration energy both at 0 and 10 GPa. We were therefore able to estimate the activation volume (∆V ) to be around 5 cm 3 /mol and d(∆V )/dP ≃ 0.The effect of pressure applies mostly on defect formation and little on migration.