Equations of state of CaSiO3 Perovskite: a molecular dynamics study

The molar volumes and bulk moduli of CaSiO₃ perovskite are calculated in the temperature range from 300 to 2,800 K and the pressure range from 0 to 143 GPa using molecular dynamics simulations that employ the breathing shell model for oxygen and the quantum correction in addition to the conventional pairwise interatomic potential models. The performance of five equations of state, i.e., the Keane, the generalized-Rydberg, the Holzapfel, the Stacey–Rydberg, and the third-order Birch–Murnaghan equations of state are examined using these data. The third-order Birch–Murnaghan equation of state is found to have a clear tendency to overestimate the bulk modulus at very high pressures. The Stacey–Rydberg equation of state degrades slightly at very high pressures along the low-temperature isotherms. In comparison, the Keane and the Holzapfel equations of state remain accurate in the whole temperature and pressure range considered in the present study. K ₀′ derived from the Holzapfel equation of state also agrees best with that calculated independently from molecular dynamics simulations. The adiabatic bulk moduli of CaSiO₃ perovskite along lower mantle geotherms are further calculated using the Keane and the Mie-Grüneisen–Debye equations of state. They are found to be constantly higher than those of the PREM by ~5%, and also very similar to those of the MgSiO₃ perovskite. Our results support the view that CaSiO₃ perovskite remains invisible in the Earth’s lower mantle.