Phase transitions of ice VIII-VII-X: A potential energy landscape perspective

Abstract

Water ice, an archetypal molecular system, exhibits a complex phase diagram characterized by numerous phase transitions under varying pressure-temperature conditions. However, it remains a significant challenge to theoretical modeling of these transitions owing to the high dimensionality in representing the potential energy landscapes (PEL) of ice crystalline phases. In an attempt to deeply understand the mechanisms of phase transitions between ice VII, VIII, and X, a graph-based stationary-point network was employed to simplify the representation of the PEL and locally constructed the PEL of ice at different pressures by our developed high-index saddle dynamic scheme. The stationary-point networks at 40 and 80 GPa indicate a pronounced flattened PEL of ice over a wide range of pressure (up to at least 80 GPa), where multiple structures (i.e., local minima) are energetically degenerate with phase VIII and separated by low-energy barriers (∼10 meV/atom at 40 GPa). These features indicate the disordered VII phase is more favorable at finite temperatures, in agreement with previous experimental observations. Additionally, we demonstrated that the PEL topology undergoes a fundamental alteration at higher pressures, with a single funnel emerging at 120 GPa, facilitating the crystallization of phase X. This work serves as a proof of concept for using the stationary-point network to explicitly represent the PEL, highlighting their potential in elucidating complex phase transitions.