The energetics of aluminum solubility into MgSiO3 perovskite at lower mantle conditions

Experiments [T. Irifune (1994) Nature 370, 131–133; E. Ito et al. (1998) Geophys. Res. Lett. 25, 821–824; A. Kubo, M. Akaogi (2000) Phys. Earth Planet. Int. 121, 85–102] indicate that (Mg,Fe)SiO3 perovskite, commonly believed to be the most abundant mineral in the Earth, is the preferred host phase of Al2O3 in the Earth’s lower mantle. Aiming to better understand the effects of Al2O3 on the thermoelastic properties of the lower mantle, we use atomistic models to examine the chemistry and elasticity of solid solutions within the MgSiO3(perovskite)–Al2O3(corundum)–MgO(periclase) mineral assemblage under conditions pertinent to the lower mantle: low Al cation concentrations, P=25–100 GPa, and T=1000–2000 K. We assess the relative stabilities of two likely substitution mechanisms of Al into MgSiO3 perovskite in terms of reactions involving MgSiO3, MgO, and Al2O3, in a manner similar to the 0 Kelvin calculations of Brodholt [J.P. Brodholt (2000) Nature 407, 620–622] and Yamamoto et al. [T. Yamamoto et al. (2003) Earth Planet. Sci. Lett. 206, 617–625]. We determine the equilibrium composition of the assemblage by examining the chemical potentials of the Al2O3 and MgO components in solid solution with MgSiO3, as functions of concentration. We find that charge coupled substitution dominates at lower mantle pressures and temperatures. Oxygen vacancy-forming substitution accounts for 3–4% of Al substitution at shallow lower mantle conditions, and less than 1% in the deep mantle. For these two pressure regimes, the corresponding adiabatic bulk moduli of aluminous perovskite are 2% and 1% lower than that of pure MgSiO3 perovskite.