First-principles calculations of elasticity of minerals at high temperature and pressure
Abstract The elasticity of minerals at high temperature and pressure ( PT ) is critical for constraining the composition and temperature of the Earth’s interior and understand better the deep water cycle and the dynamic Earth. First-principles calculations without introducing any adjustable parameters, whose results can be comparable to experimental data, play a more and more important role in investigating the elasticity of minerals at high PT mainly because of (1) the quick increasing of computational powers and (2) advances in method. For example, the new method reduces the computation loads to one-tenth of the traditional method with the comparable precise as the traditional method. This is extraordinarily helpful because first-principles calculations of the elasticity of minerals at high PT are extremely time-consuming. So far the elasticity of most of lower mantle minerals has been investigated in detail. We have good idea on the effect of temperature, pressure, and iron concentration on elasticity of main minerals of the lower mantle and the unusual softening in bulk modulus by the spin crossover of iron in ferropericlase. With these elastic data the lower mantle has been constrained to have 10–15 wt% ferropericlase, which is sufficient to generate some visible effects of spin crossover in seismic tomography. For example, the spin crossover causes that the temperature sensitivity of P wave at the depth of ~1700 km is only a fraction of that at the depth of ~2300 km. The disruptions of global P wave structure and of P wave image below hotspots such as Hawaii and Iceland at similar depth are in consistence with the spin crossover effect of iron in ferropericlase. The spin crossover, which causes anomalous thermodynamic properties of ferropericlase, has also been found to play a control role for the two features of the large low shear velocity provinces (LLSVPs): the sharp edge and high elevation up to 1000 km above core-mantle boundary. All these results clearly suggest the spin crossover of iron in the lower mantle. The theoretical investigations for the elasticity of minerals at the upper mantle and water effect on elasticity of minerals at the mantle transition zone and subducting slab have also been conducted extensively. These researches are critical for understanding better the composition of the upper mantle and water distribution and transport in the Earth’s mantle. Most of these were static calculations, which did not include the vibrational (temperature) effect on elasticity, although temperature effect on elasticity is basic because of high temperature at the Earth’s interior and huge temperature difference between the ambient mantle and the subducting slab. Including temperature effect on elasticity of minerals should be important future work. New method developed is helpful for these directions. The elasticity of iron and iron-alloy with various light elements has also been calculated extensively. However, more work is necessary in order to meet the demand for constraining the types and amount of light elements at the Earth’s core.