Quasiharmonic thermal elasticity of crystals: An analytical approach

First-principles quasiharmonic calculations play a very important role in mineral physics because they can predict the structural and thermodynamic properties of materials at pressure and temperature conditions of the Earth's interior that are still challenging for experiments. They also enable calculations of thermal elastic properties by providing second-order derivatives of free energies with respect to strain. The latter are essential to interpret seismic tomography of the mantle in terms of temperature, composition, and mineralogy, in the context of geophysical processes. However, these are exceedingly demanding computations requiring up to ∼103 parallel jobs running on tens or more processors each. Here we introduce an analytical and computationally simpler approach that requires only calculations of static elastic constants and phonon density of states for unstrained configurations. This approach, currently implemented for crystals with up to orthorhombic symmetry, decreases the computational effort, i.e., CPU time and human labor, by up to two orders of magnitude. Results for the major mantle phases periclase (MgO) and forsterite (α-Mg₂SiO₄) show excellent agreement with previous first-principles results and experimental data.