Phase transition and elasticity of enstatite under pressure from experiments and first-principles studies

We have investigated the thermodynamic stability, crystal structure, elastic constants, and sound velocities of MgSiO3-enstatite using data from X-ray diffraction and ultrasonic measurements up to 16.8 GPa and first-principles calculations up to 30 GPa. The calculated enthalpies provide theoretical support for the phase transition from Pbca to P21/c between 9 and 14 GPa previously observed in natural orthoensatite and MgSiO3 enstatite. A density increase of 1.4–1.5% for the Pbca → P21/c transition is obtained from both first-principles and experimental studies. Elastic constants of Pbca, P21/c, C2/c, P21ca and Pbca-II are all calculated, and a softening in the shear constant C55 is predicted for Pbca and Pbca-II phases. C55 of Pbca is found to be closely correlated with the A-site SiO4 tetrahedra chain angle while C44 and C66 are correlated with the B-site chain angle. Pbca, P21/c and C2/c all exhibit similar volumetric compressibilities at all pressures. The calculated velocities of the P21/c phase at 12 GPa are equal to those of Pbca for P and 1.3% higher for S waves. The experimentally observed P and S wave velocity anomalies can be qualitatively described by the transformation from Pbca to P21/c; however, the magnitudes of the velocity decreases between 10 and 14 GPa remain to be verified by future single crystal data or polycrystalline measurements at high pressures. The predicted velocity jumps of 2.8% and 4.5% for P and S waves, respectively, between Pbca and C2/c in the pressure range of 5–12 GPa are in excellent agreement with the values of ∼3(1)% and ∼5(1)% obtained from the directly measured data, thereby making it a plausible candidate for the seismic X-discontinuity at depths of 250–300 km in the Earth’s upper mantle.