Pressure‐volume‐temperature equation of state of MgSiO3 perovskite from molecular dynamics and constraints on lower mantle composition

The composition of the lower mantle can be investigated by examining densities and seismic velocities of compositional models as functions of depth. In order to do this it is necessary to know the volumes and thermoelastic properties of the compositional constituents under lower mantle conditions. We determined the thermal equation of state (EOS) of MgSiO3 perovskite using the nonempirical variational induced breathing (VIB) interatomic potential with molecular dynamics simulations at pressures and temperatures of the lower mantle. We fit our pressure‐volume‐temperature results to a thermal EOS of the form P ( V , T ) = P 0( V , T 0) + Δ P th( T ), where T 0 = 300 K and P 0 is the isothermal Universal EOS. The thermal pressure Δ P th can be represented by a linear relationship Δ P th = a + bT . We find V 0 = 165.40 Å 3 , K T 0 = 273 GPa, K T 0′ = 3.86, a = −1.99 GPa, and b = 0.00664 GPa K −1 for pressures of 0–140 GPa and temperatures of 300–3000 K. By fixing V 0 to the experimentally determined value of 162.49 Å 3 and calculating density and bulk sound velocity profiles along a lower mantle geotherm we find that the lower mantle cannot consist solely of (Mg,Fe)SiO3 perovskite with XMg ranging from 0.9–1.0. Using pyrolitic compositions of 67 vol % perovskite (XMg = 0.93–0.96) and 33 vol % magnesiowüstite (XMg = 0.82–0.86), however, we obtained density and velocity profiles that are in excellent agreement with seismological models for a reasonable geotherm.