DFT study of migration enthalpies in MgSiO3 perovskite

The effect of pressure on ionic diffusion in orthorhombic MgSiO3 perovskite has been investigated using density functional theory. An intensive investigation of possible silicon pathways revealed new positions of the saddle-points and an enthalpy of migration at 26.2 GPa of 4.7 eV that is in fair agreement with the experimental values of about 3.5 eV at 25 GPa. This is much lower than found in previous studies (~9 eV) and removes the need to explain silicon diffusion by a complicated process involving coupled oxygen vacancies, as has been previously proposed. Our migration enthalpies for oxygen and magnesium are in excellent agreement with experiments. We find that oxygen diffusion occurs via a chain of several inequivalent jumps along the octahedron edges, and that magnesium occurs via two inequivalent [110] jumps and one [001] jump. We also present activation volumes for all three species at 25 and 135 GPa.