The range of meta stability of ice-water melting for two simple models of water

Abstract

A number of crystal structures of water have been `superheated' in Monte Carlo simulations. Two well known models for water were considered; namely the TIP4P model and the SPC/E model. By comparing the fluid-solid coexistence temperature to the temperature at which the solid becomes mechanically unstable and melts it is possible to determine the typical range of temperatures over which is possible to superheat the ice phases in conventional simulation studies. It is found that the ice phases can be superheated to approximately 90K beyond the fluid-solid coexistence temperature. Beyond this limit they spontaneously melt. This limit appears to depend weakly both on the type of ice phase considered and on the chosen model. Obviously only rigorous free energy calculations can determine the equilibrium fluid-solid coexistence of a model. However, a "rule of thumb" is that, by subtracting 90K from the mechanically stability limit of the the ice phase one is provided with a first guess as to the equilibrium fluid-solid coexistence temperature.

The preprint of the paper published in Molecular Physics, volume 103 pp. 1-5 (2005)

Cited by in corpus (6)

• 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
• The phase diagram of water at high pressures as obtained by computer simulations of the TIP4P/2005 model: the appearance of a plastic crystal phase
• Absence of superheating for ice Ih with a free surface : a new method of determining the melting point of different water models
• Non-Markovian melting: a novel procedure to generate initial liquid like phases for small molecules for use in computer simulation studies