MOLECULAR DYNAMICS SIMULATIONS OF THE PHASE TRANSITION BETWEEN LOW-TEMPERATURE AND HIGH-TEMPERATURE CLINOENSTATITES

Abstract

Phase transition between low-temperature clinoenstatite (LT-CEn) and high-temperature clinoenstatite (HT-CEn) was studied by using molecular dynamics (MD) simulations, based on empirical potential parameters. Starting from LT-CEn, the MD calculations were carried out at atmospheric pressure and at elevated pressures (1-6 GPa). At elevated temperatures the transformation from the starting LT-CEn to HT-CEn occurred at any pressure. It was confirmed that the HT-CEn has the same space group C2/c as diopside but the M2 site is six-coordinated, unlike diopside. A significant difference in the MD-simulated cell volumes between LT-CEn and HT-CEn was also observed, showing a first-order transition. In addition, there were some temperature ranges where LT-CEN and HT-CEn would be coexistent and very small thermal hystereses between increasing and decreasing temperatures during the transition. These behaviors are consistent with the characteristic of a thermoelastic-martensitic transformation. The phase boundary between LT-CEn and HT-CEn was determined for the first time. Its positive dT/dP slope strongly shows that the high-pressure clinoenstatite is a significantly distinct phase from HT-CEn although the both phases have the same space group, C2/c.