Structural and elastic anisotropy of crystals at high pressures and temperatures from quantum mechanical methods: The case of Mg2SiO4 forsterite

Abstract

We report accurate ab initio theoretical predictions of the elastic, seismic, and structural anisotropy of the orthorhombic Mg2SiO4 forsterite crystal at high pressures (up to 20 GPa) and temperatures (up to its melting point, 2163 K), which constitute earth’s upper mantle conditions. Single-crystal elastic stiffness constants are evaluated up to 20 GPa and their first- and second-order pressure derivatives reported. Christoffel’s equation is solved at several pressures: directional seismic wave velocities and related properties (azimuthal and polarization seismic anisotropies) discussed. Thermal structural and average elastic properties, as computed within the quasi-harmonic approximation of the lattice potential, are predicted at high pressures and temperatures: directional thermal expansion coefficients, first- and second-order pressure derivatives of the isothermal bulk modulus, and P-V-T equation-of-state. The effect on computed properties of five different functionals, belonging to three different classes of approximations, of the density functional theory is explicitly investigated.