Elasticity of CaSiO(3) perovskite at high pressure and high temperature

Ab initio molecular dynamic (AIMD) simulations were performed to calculate the equation of state (EOS) of CaSiO(3) perovskite at mantle pressure-temperature conditions. At temperatures above 2000 K, even though the hydrostatic crystal structure is metrically tetragonal in the pressure range of 13-123 GPa, the symmetry of the elastic moduli is consistent with cubic symmetry. Our results show that elastic constants and velocities are independent of temperature at constant volume. Referenced to room pressure and 2000 K, we find: Gruneisen parameter is gamma(V) = gamma(0)(V/V(0))(q) with gamma(0) = 1.53 and q = 1.02(5), and the Anderson Gruneisen parameter is given by (alpha/alpha(0)) = (V/V(0))(delta T) in which alpha(0) = 2.89 x 10(-5) K(-1) and delta(T) = 4.09(5). Using the third order Birch Murnaghan equation of state to fit our data, we have for ambient P and T, K(0) = 236.6(8) GPa, K(0)' = 3.99(3), and V(0) = 729.0(6) angstrom(3). Calculated acoustic velocities show the following P-T dependence: (partial derivative ln V(P)/partial derivative V)(T) (or P) = 1.9 x 10(-3); (partial derivative ln V(S)/partial derivative V)(T or P) = -1.5 x 10(-3); (partial derivative ln V(Phi)/partial derivative V)(T or) (P) = -2.4 x 10(-3); (partial derivative ln V(S)/partial derivative ln V(P))(T or P) = 0.79; (partial derivative ln V(S)/partial derivative ln V(Phi))(T or P) = 0.63, indicating that the variations in bulk modulus overpower the variations in shear modulus. The bulk modulus of CaSiO(3) perovskite is up to 10% lower than MgSiO(3) perovskite under lower mantle conditions. The difference diminishes with pressure and temperature. The shear modulus of CaSiO(3) perovskite is almost 25% lower compared with MgSiO(3) perovskite for shallow lower mantle pressures and temperatures and about 3% lower at the base of the lower mantle. The difference in density of these two perovskite is about 3-4% for all conditions. Both the density and bulk modulus differ from PREM by less than 2% throughout the lower mantle. The shear modulus is similar to 10% lower at shallow depths grading to similar to 5% by the core-mantle boundary. Thus, the seismic velocity of CaSiO(3) perovskite will be lower (0-6%) than PREM. (c) 2006 Elsevier B.V. All rights reserved