Seismic velocities of CaSiO3 perovskite can explain LLSVPs in Earth's lower mantle.

Seismology records the presence of various heterogeneities throughout the lower mantle1,2, but the origins of these signals-whether thermal or chemical-remain uncertain, and therefore much of the information that they hold about the nature of the deep Earth is obscured. Accurate interpretation of observed seismic velocities requires knowledge of the seismic properties of all of Earth's possible mineral components. Calcium silicate (CaSiO3) perovskite is believed to be the third most abundant mineral throughout the lower mantle. Here we simultaneously measure the crystal structure and the shear-wave and compressional-wave velocities of samples of CaSiO3 perovskite, and provide direct constraints on the adiabatic bulk and shear moduli of this material. We observe that incorporation of titanium into CaSiO3 perovskite stabilizes the tetragonal structure at higher temperatures, and that the material's shear modulus is substantially lower than is predicted by computations3-5 or thermodynamic datasets6. When combined with literature data and extrapolated, our results suggest that subducted oceanic crust will be visible as low-seismic-velocity anomalies throughout the lower mantle. In particular, we show that large low-shear-velocity provinces (LLSVPs) are consistent with moderate enrichment of recycled oceanic crust, and mid-mantle discontinuities can be explained by a tetragonal-cubic phase transition in Ti-bearing CaSiO3 perovskite.

First-principles constraints on diffusion in lower-mantle minerals and a weak D'' layer.

Post-perovskite MgSiO(3) is believed to be present in the D'' region of the Earth's lowermost mantle. Its existence has been used to explain a number of seismic observations, such as the D'' reflector and the high degree of seismic anisotropy within the D'' layer. Ionic diffusion in post-perovskite controls its viscosity, which in turn controls the thermal and chemical coupling between the core and the mantle, the development of plumes and the stability of deep chemical reservoirs. Here we report the use of first-principles methods to calculate absolute diffusion rates in post-perovskite under the conditions found in the Earth's lower mantle. We find that the diffusion of Mg(2+) and Si(4+) in post-perovskite is extremely anisotropic, with almost eight orders of magnitude difference between the fast and slow directions. If post-perovskite in the D'' layer shows significant lattice-preferred orientation, the fast diffusion direction will render post-perovskite up to four orders of magnitude weaker than perovskite. The presence of weak post-perovskite strongly increases the heat flux across the core-mantle boundary and alters the geotherm. It also provides an explanation for laterally varying viscosity in the lowermost mantle, as required by long-period geoid models. Moreover, the behaviour of very weak post-perovskite can reconcile seismic observation of a D'' reflector with recent experiments showing that the width of the perovskite-to-post-perovskite transition is too wide to cause sharp reflectors. We suggest that the observed sharp D'' reflector is caused by a rapid change in seismic anisotropy. Once sufficient perovskite has transformed into post-perovskite, post-perovskite becomes interconnected and strain is partitioned into this weaker phase. At this point, the weaker post-perovskite will start to deform rapidly, thereby developing a strong crystallographic texture. We show that the expected seismic contrast between the deformed perovskite-plus-post-perovskite assemblage and the overlying isotropic perovskite-plus-post-perovskite assemblage is consistent with seismic observations.

Possible thermal and chemical stabilization of body-centred-cubic iron in the Earth's core.

The nature of the stable phase of iron in the Earth's solid inner core is still highly controversial. Laboratory experiments suggest the possibility of an uncharacterized phase transformation in iron at core conditions and seismological observations have indicated the possible presence of complex, inner-core layering. Theoretical studies currently suggest that the hexagonal close packed (h.c.p.) phase of iron is stable at core pressures and that the body centred cubic (b.c.c.) phase of iron becomes elastically unstable at high pressure. In other h.c.p. metals, however, a high-pressure b.c.c. form has been found to become stabilized at high temperature. We report here a quantum mechanical study of b.c.c.-iron able to model its behaviour at core temperatures as well as pressures, using ab initio molecular dynamics free-energy calculations. We find that b.c.c.-iron indeed becomes entropically stabilized at core temperatures, but in its pure state h.c.p.-iron still remains thermodynamically more favourable. The inner core, however, is not pure iron, and our calculations indicate that the b.c.c. phase will be stabilized with respect to the h.c.p. phase by sulphur or silicon impurities in the core. Consequently, a b.c.c.-structured alloy may be a strong candidate for explaining the observed seismic complexity of the inner core.

The ab initio simulation of the Earth's core.

The Earth has a liquid outer and solid inner core. It is predominantly composed of Fe, alloyed with small amounts of light elements, such as S, O and Si. The detailed chemical and thermal structure of the core is poorly constrained, and it is difficult to perform experiments to establish the properties of core-forming phases at the pressures (ca. 300 GPa) and temperatures (ca. 5000-6000 K) to be found in the core. Here we present some major advances that have been made in using quantum mechanical methods to simulate the high-P/T properties of Fe alloys, which have been made possible by recent developments in high-performance computing. Specifically, we outline how we have calculated the Gibbs free energies of the crystalline and liquid forms of Fe alloys, and so conclude that the inner core of the Earth is composed of hexagonal close packed Fe containing ca. 8.5% S (or Si) and 0.2% O in equilibrium at 5600 K at the boundary between the inner and outer cores with a liquid Fe containing ca. 10% S (or Si) and 8% O.

Theoretical investigation of metastable Al2SiO5 polymorphs.

Using theoretical simulations based on density functional theory within the generalized gradient approximation, a series of metastable phase transitions occurring in low-pressure Al2SiO5 polymorphs (andalusite and sillimanite) are predicted; similar results were obtained using semiclassical interatomic potentials within the ionic shell model. Soft lattice modes as well as related structural changes are analysed. For sillimanite, an isosymmetric phase transition at ca 35 GPa is predicted; an incommensurately modulated form of sillimanite can also be obtained at low temperatures and high pressures. The high-pressure isosymmetric phase contains five-coordinate Si and Al atoms. The origin of the fivefold coordination is discussed in detail. Andalusite was found to transform directly into an amorphous phase at ca 50 GPa. This study provides an insight into the nature of metastable modifications of crystal structures and the ways in which they are formed. Present results indicate the existence of a critical bonding distance, above which interatomic interactions cannot be considered as bonding. The critical distance for the Si-O bond is 2.25 A.

The elastic constants of MgSiO3 perovskite at pressures and temperatures of the Earth's mantle.

The temperature anomalies in the Earth's mantle associated with thermal convection can be inferred from seismic tomography, provided that the elastic properties of mantle minerals are known as a function of temperature at mantle pressures. At present, however, such information is difficult to obtain directly through laboratory experiments. We have therefore taken advantage of recent advances in computer technology, and have performed finite-temperature ab initio molecular dynamics simulations of the elastic properties of MgSiO3 perovskite, the major mineral of the lower mantle, at relevant thermodynamic conditions. When combined with the results from tomographic images of the mantle, our results indicate that the lower mantle is either significantly anelastic or compositionally heterogeneous on large scales. We found the temperature contrast between the coldest and hottest regions of the mantle, at a given depth, to be about 800 K at 1,000 km, 1,500 K at 2,000 km, and possibly over 2,000 K at the core-mantle boundary.