Premium
Molecular dynamics simulation of structures, bulk moduli, and volume thermal expansivities of silicate liquids in the system CaO‐MgO‐Al 2 O 3 ‐SiO 2
Author(s) -
Matsui Masanori
Publication year - 1996
Publication title -
geophysical research letters
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.007
H-Index - 273
eISSN - 1944-8007
pISSN - 0094-8276
DOI - 10.1029/96gl00260
Subject(s) - diopside , forsterite , enstatite , anorthite , silicate , bulk modulus , mineralogy , wollastonite , thermodynamics , materials science , mantle (geology) , solidus , geology , chemistry , meteorite , geophysics , composite material , chondrite , raw material , physics , organic chemistry , alloy , astronomy
Materials in the system CaO‐MgO‐Al 2 O 3 ‐SiO 2 are important constituents of the Earth's lower crust and mantle. Silicate liquids in this geophysically important system have been studied using molecular dynamics(MD) simulation with an empirical interatomic potential (CMAS94). MD simulations are quite satisfactory in reproducing well the observed structure and pressure‐volume‐temperature equation‐of‐state parameters of molten enstatite (MgSiO 3 ), wollastonite (CaSiO 3 ), diopside (CaMgSi 2 O 6 ), and anorthite (CaAl 2 Si 2 O 8 ) at 1900 K and 0 GPa. However, the MD simulated bulk modulus of molten forsterite (Mg 2 SiO 4 ) at 2300 K and 0 GPa, K 0 =18.0(6) GPa, is found to be much smaller than the value, K 0 =∼60 GPa at similar temperature and pressure conditions, estimated previously based on melting curve analyses of forsterite. In an attempt to investigate the possible occurrence of the density inversion between magmatic liquids and residual crystals in the upper mantle conditions, as proposed by Stolper et al . [1981], we have further applied the MD technique with the CMAS94 potential to the diopside system at 1900 K as an example, and have found that the density inversion between crystal and liquid in this system actually occurs at approximately 11 GPa.