The ice-vapor interface and the melting point of ice Ih for the polarizable POL3 water model

Eva Muchová, Ivan Gladich, Sylvain Picaud, Paul N.M. Hoang, Martina Roeselová

Research output: Contribution to journalArticle

14 Citations (Scopus)

Abstract

We use molecular dynamics simulations to determine the melting point of ice Ih for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys. 1992, 97, 2659). Simulations are performed on a slab of ice Ih with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields themelting point of the POL3 model to be Tm = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.

Original languageEnglish
Pages (from-to)5973-5982
Number of pages10
JournalJournal of Physical Chemistry A
Volume115
Issue number23
DOIs
Publication statusPublished - 16 Jun 2011
Externally publishedYes

Fingerprint

Ice
Freezing
melting points
Melting point
ice
Vapors
vapors
Water
water
field theory (physics)
simulation
Temperature
Undercooling
supercooling
Molecular Dynamics Simulation
yield point
Phase diagrams
temperature
Molecular dynamics
Hydrogen

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

The ice-vapor interface and the melting point of ice Ih for the polarizable POL3 water model. / Muchová, Eva; Gladich, Ivan; Picaud, Sylvain; Hoang, Paul N.M.; Roeselová, Martina.

In: Journal of Physical Chemistry A, Vol. 115, No. 23, 16.06.2011, p. 5973-5982.

Research output: Contribution to journalArticle

Muchová, Eva ; Gladich, Ivan ; Picaud, Sylvain ; Hoang, Paul N.M. ; Roeselová, Martina. / The ice-vapor interface and the melting point of ice Ih for the polarizable POL3 water model. In: Journal of Physical Chemistry A. 2011 ; Vol. 115, No. 23. pp. 5973-5982.
@article{15717d0472df4b0aa699485ff2ff56ce,
title = "The ice-vapor interface and the melting point of ice Ih for the polarizable POL3 water model",
abstract = "We use molecular dynamics simulations to determine the melting point of ice Ih for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys. 1992, 97, 2659). Simulations are performed on a slab of ice Ih with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields themelting point of the POL3 model to be Tm = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.",
author = "Eva Muchov{\'a} and Ivan Gladich and Sylvain Picaud and Hoang, {Paul N.M.} and Martina Roeselov{\'a}",
year = "2011",
month = "6",
day = "16",
doi = "10.1021/jp110391q",
language = "English",
volume = "115",
pages = "5973--5982",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "23",

}

TY - JOUR

T1 - The ice-vapor interface and the melting point of ice Ih for the polarizable POL3 water model

AU - Muchová, Eva

AU - Gladich, Ivan

AU - Picaud, Sylvain

AU - Hoang, Paul N.M.

AU - Roeselová, Martina

PY - 2011/6/16

Y1 - 2011/6/16

N2 - We use molecular dynamics simulations to determine the melting point of ice Ih for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys. 1992, 97, 2659). Simulations are performed on a slab of ice Ih with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields themelting point of the POL3 model to be Tm = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.

AB - We use molecular dynamics simulations to determine the melting point of ice Ih for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys. 1992, 97, 2659). Simulations are performed on a slab of ice Ih with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields themelting point of the POL3 model to be Tm = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.

UR - http://www.scopus.com/inward/record.url?scp=79959929008&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=79959929008&partnerID=8YFLogxK

U2 - 10.1021/jp110391q

DO - 10.1021/jp110391q

M3 - Article

C2 - 21452834

AN - SCOPUS:79959929008

VL - 115

SP - 5973

EP - 5982

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 23

ER -