vix.ing · top · new · best · stats · spec

Pressure‐volume‐temperature equation of state of MgSiO3 perovskite from molecular dynamics and constraints on lower mantle composition

2000/10/31 by Frederic C. Marton, Joel Ita, Ronald E. Cohen
Earth and Planetary Sciences · Physics and Astronomy · #Geological and Geochemical Analysis #High-pressure geophysics and materials #earthquake and tectonic studies #physics.geo-ph

paper · pdf · doi:10.1029/2000jb900457

published as J. Geophys. Res., 106, 8615-8627, 2001 · 26 pages, 11 figures, 3 tables. Uses fixed version of agupp.sty

arxiv created 2001/03/14 · openalex publication_date 2001/05/10 · arxiv updated 2009/12/01 · openalex created_date 2016/06/24 · openalex updated_date 2026/07/28

Abstract

The composition of the lower mantle can be investigated by examining densities and seismic velocities of compositional models as functions of depth. In order to do this it is necessary to know the volumes and thermoelastic properties of the compositional constituents under lower mantle conditions. We determined the thermal equation of state (EOS) of MgSiO3 perovskite using the nonempirical variational induced breathing (VIB) interatomic potential with molecular dynamics simulations at pressures and temperatures of the lower mantle. We fit our pressure‐volume‐temperature results to a thermal EOS of the form P ( V , T ) = P 0( V , T 0) + Δ P th( T ), where T 0 = 300 K and P 0 is the isothermal Universal EOS. The thermal pressure Δ P th can be represented by a linear relationship Δ P th = a + bT . We find V 0 = 165.40 Å 3 , K T 0 = 273 GPa, K T 0′ = 3.86, a = −1.99 GPa, and b = 0.00664 GPa K −1 for pressures of 0–140 GPa and temperatures of 300–3000 K. By fixing V 0 to the experimentally determined value of 162.49 Å 3 and calculating density and bulk sound velocity profiles along a lower mantle geotherm we find that the lower mantle cannot consist solely of (Mg,Fe)SiO3 perovskite with XMg ranging from 0.9–1.0. Using pyrolitic compositions of 67 vol % perovskite (XMg = 0.93–0.96) and 33 vol % magnesiowüstite (XMg = 0.82–0.86), however, we obtained density and velocity profiles that are in excellent agreement with seismological models for a reasonable geotherm.

Related