Absence of superheating for ice Ih with a free surface : a new method of determining the melting point of different water models

Abstract

Molecular dynamic simulations were performed for ice Ih with a free surface. The simulations were carried out at several temperatures and each run lasted more than 7ns. At high temperatures the ice melts. It is demonstrated that the melting process starts at the surface and propagates to the bulk of the ice block. Already at the temperatures below the melting point, we observe a thin liquid layer at the ice surface, but the block of ice remains stable along the run. As soon as the temperature reaches the melting point the entire ice block melts. Our results demonstrate that, unlike in the case of conventional simulations in the NpT ensemble, overheating of the ice Ih with a free surface does not occur. That allows to estimate the melting point of ice at zero pressure. We applied the method to the following models of water: SPC/E, TIP4P, TIP4P/Ew, TIP4P/Ice and TIP4P/2005, and found good agreement between the melting temperatures obtained by this procedure and the values obtained either from free energy calculations or from direct simulations of the ice/water interface.

References in corpus (2)

• Melting and nonmelting of solid surfaces and nanosystems
• The range of meta stability of ice-water melting for two simple models of water

Cited by in corpus (10)

• What ice can teach us about water interactions: a critical comparison of the performance of different water models
• Determination of phase diagrams via computer simulation: Methodology and applications to water, electrolytes and proteins
• The thickness of a liquid layer on the free surface of ice as obtained from computer simulation
• Intermolecular forces at ice and water interfaces: premelting, surface freezing and regelation
• The kinetics of the ice-water interface from ab initio machine learning simulations
• Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations
• Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate
• Effect of Sodium Chloride Adsorption on the Surface Premelting of Ice
• Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation
• Is it possible to overheat ice? The activated melting of TIP4P/Ice at solid-vapor coexistence