Structural Transformations of Ice at High Pressures Via Molecular Dynamics Simulations

Abstract

A simple classical model is used for the study of the structural transformations of ice under high pressures, such as ice VIII to VII and X, via classical molecular dynamics (MD) simulation. In the present MD simulation, pair potentials of a simple form between pair of atoms and a thee-body potential representing the H-O-H angle dependence, originally developed by Kawamura et al., were used. Starting with a stable ice VIII at low pressure and low temperature, we have carried out two different MD runs, one with increasing pressure keeping the temperature constant (simulation I) and the other with increasing temperature under constant pressure (simulation II). From these MD simulations we have obtained the structural transformations from ice VIII to VII for both simulations; the former was finally transformed into ice X for the simulation I. The present results are compatible with recent experiments on high pressure ices.