Heating- and pressure-induced transformations in amorphous and hexagonal ice: A computer simulation study using the TIP4P/2005 model.

Abstract

We characterize the phase behavior of glassy water by performing extensive out-of-equilibrium molecular dynamics simulations using the TIP4P/2005 water model. Specifically, we study (i) the pressure-induced transformations between low-density (LDA) and high-density amorphous ice (HDA), (ii) the pressure-induced amorphization (PIA) of hexagonal ice (Ih), (iii) the heating-induced LDA-to-HDA transformation at high pressures, (iv) the heating-induced HDA-to-LDA transformation at low and negative pressures, (v) the glass transition temperatures of LDA and HDA as a function of pressure, and (vi) the limit of stability of LDA upon isobaric heating and isothermal decompression (at negative pressures). These transformations are studied systematically, over a wide range of temperatures and pressures, allowing us to construct a P-T phase diagram for glassy TIP4P/2005 water. Our results are in qualitative agreement with experimental observations and with the P-T phase diagram obtained for glassy ST2 water that exhibits a liquid-liquid phase transition and critical point. We also discuss the mechanism for PIA of ice Ih and show that this is a two-step process where first, the hydrogen-bond network (HBN) is distorted and then the HBN abruptly collapses. Remarkably, the collapse of the HB in ice Ih occurs when the average molecular orientations order, a measure of the tetrahedrality of the HBN, is of the same order as in LDA, suggesting a common mechanism for the LDA-to-HDA and Ih-to-HDA transformations.