Premium
Comparison between the lattice dynamics and molecular dynamics methods: Calculation results for MgSiO 3 perovskite
Author(s) -
Matsui Masanori,
Price Geoffrey D.,
Patel Atul
Publication year - 1994
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/94gl01370
Subject(s) - anharmonicity , molecular dynamics , bulk modulus , thermodynamics , materials science , molar volume , volume (thermodynamics) , lattice (music) , atmospheric temperature range , thermal expansion , perovskite (structure) , mineralogy , chemistry , condensed matter physics , physics , computational chemistry , crystallography , acoustics
The lattice dynamics (LD) and molecular dynamics (MD) methods have been used to calculate the structure, bulk modulus, and volume thermal expansivity of MgSiO 3 perovskite, in order to investigate the reliability of the two simulation techniques over a wide range of temperature and pressure conditions. At an intermediate temperature of 500 K and zero pressure, the LD and MD values are in exellent agreement for both the structure and bulk modulus of MgSiO 3 perovskite. At high temperatures and zero pressure, however, the LD method, which is based on the quasi‐harmonic approximation, increasingly overestimates the molar volume of MgSiO 3 perovskite because of the neglect of higher‐order anharmonic terms. At the high temperatures and high pressures prevailing in the lower mantle, the errors in the LD values for both the molar volume and bulk modulus, relative to the MD values, are generally small or negligible. However, since anharmonicity decreases substantially with pressure but increases rapidly with temperature, the error in the LD simulated volume thermal expansivity is serious, especially in the lower pressure region.