Open AccessSelf‐consistent thermodynamic description of silicate liquids, with application to shock melting of MgO periclase and MgSiO 3 perovskite
Author(s) -
De Koker Nico,
Stixrude Lars
Publication year - 2009
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2009.04142.x
Subject(s) - periclase , enstatite , thermodynamics , mantle (geology) , silicate , perovskite (structure) , melting curve analysis , melting temperature , materials science , stishovite , mineralogy , meteorite , geology , spinel , geophysics , high pressure , chemistry , physics , metallurgy , crystallography , chondrite , polymerase chain reaction , biochemistry , organic chemistry , composite material , astronomy , gene
SUMMARY We develop a self‐consistent thermodynamic description of silicate liquids applicable across the entire mantle pressure and temperature regime. The description combines the finite strain free energy expansion with an account of the temperature dependence of liquid properties into a single fundamental relation, while honouring the expected limiting behaviour at large volume and high temperature. We find that the fundamental relation describes well previous experimental and theoretical results for liquid MgO, MgSiO 3 , Mg 2 SiO 4 and SiO 2 . We apply the description to calculate melting curves and Hugoniots of solid and liquid MgO and MgSiO 3 . For periclase, we find a melting temperature at the core–mantle boundary (CMB) of 7810 ± 160 K , with the solid Hugoniot crossing the melting curve at 375 GPa, 9580 K , and the liquid Hugoniot crossing at 470 GPa, 9870 K . For complete shock melting of periclase we predict a density increase of 0.14 g cm −3 and a sound speed decrease of 2.2 km s −1 . For perovskite, we find a melting temperature at the CMB of 5100 ± 100 K with the perovskite section of the enstatite Hugoniot crossing the melting curve at 150 GPa, 5190 K , and the liquid Hugoniot crossing at 220 GPa, 5520 K . For complete shock melting of perovskite along the enstatite principal Hugoniot, we predict a density increase of 0.10 g cm −3 , with a sound speed decrease of 2.6 km s −1 .