Testing recent charge-on-spring type polarizable water models. I. Melting temperature and ice properties

Péter T. Kiss, Péter Bertsyk, A. Baranyai

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

We determined the freezing point of eight molecular models of water. All models use the charge-on-spring (COS) method to express polarization. The studied models were the COSG2, COSG3 H. Yu and W. F. van Gunsteren, J. Chem. Phys. 121, 9549 (2004)10.1063/1.1805516, the COSB2 H. Yu, T. Hansson, and W. F. van Gunsteren, J. Chem. Phys. 118, 221 (2003)10.1063/1.1523915, the SWM4-DP G. Lamoureux, A. D. MacKerell, Jr., and B. Roux, J. Chem. Phys. 119, 5185 (2003)10.1063/1.1598191, the SWM4-NDP G. Lamoureux, E. Harder, I. V. Vorobyov, B. Roux, and A. D. MacKerell, Jr., Chem. Phys. Lett. 418, 245 (2006)10.1016/j.cplett.2005.10.135, and three versions of our model, the BKd1, BKd2, and BKd3. The BKd1 is the original Gaussian model P. T. Kiss, M. Darvas, A. Baranyai, and P. Jedlovszky, J. Chem. Phys. 136, 114706 (2012)10.1063/1. 3692602 with constant polarization and with a simple exponential repulsion. The BKd2 applies field-dependent polarizability A. Baranyai and P. T. Kiss, J. Chem. Phys. 135, 234110 (2011)10.1063/1.3670962, while the BKd3 model has variable size to approximate the temperature-density (T-ρ) curve of water P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)10.1063/1.4746419. We used the thermodynamic integration (TI) and the Gibbs-Helmholtz equation to determine the equality of the free energy for liquid water and hexagonal ice (Ih) at 1 bar. We used the TIP4P and the SPCE models as reference systems of the TI. The studied models severely underestimate the experimental melting point of ice Ih. The calculated freezing points of the models are the following: COSG2, 215 K; COSG3, 149 K; SWM4-DP, 186 K; BKd1, 207 K; BKd2, 213 K; BKd3, 233 K. The freezing temperature of the SWM4-NDP system is certainly below 120 K. It might even be that the water phase is more stable than the ice Ih at 1 bar in the full temperature range. The COSB2 model melts below 100 K. The best result was obtained for the BKd3 model which indicates that correct description of the (T-ρ) curve improves the estimation of the freezing point. We also determined and compared the densities of high-pressure polymorphs of ice for these models.

Original languageEnglish
Article number194102
JournalThe Journal of Chemical Physics
Volume137
Issue number19
DOIs
Publication statusPublished - Nov 21 2012

Fingerprint

Ice
Melting point
ice
melting
Water
Testing
water
Freezing
melting points
temperature
Gibbs-Helmholtz equations
Thermodynamics
Polarization
thermodynamics
Helmholtz equation
reference systems
curves
polarization
Polymorphism
Temperature

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

Testing recent charge-on-spring type polarizable water models. I. Melting temperature and ice properties. / Kiss, Péter T.; Bertsyk, Péter; Baranyai, A.

In: The Journal of Chemical Physics, Vol. 137, No. 19, 194102, 21.11.2012.

Research output: Contribution to journalArticle

@article{e87a3021929a4a7da5dff177de7dd769,
title = "Testing recent charge-on-spring type polarizable water models. I. Melting temperature and ice properties",
abstract = "We determined the freezing point of eight molecular models of water. All models use the charge-on-spring (COS) method to express polarization. The studied models were the COSG2, COSG3 H. Yu and W. F. van Gunsteren, J. Chem. Phys. 121, 9549 (2004)10.1063/1.1805516, the COSB2 H. Yu, T. Hansson, and W. F. van Gunsteren, J. Chem. Phys. 118, 221 (2003)10.1063/1.1523915, the SWM4-DP G. Lamoureux, A. D. MacKerell, Jr., and B. Roux, J. Chem. Phys. 119, 5185 (2003)10.1063/1.1598191, the SWM4-NDP G. Lamoureux, E. Harder, I. V. Vorobyov, B. Roux, and A. D. MacKerell, Jr., Chem. Phys. Lett. 418, 245 (2006)10.1016/j.cplett.2005.10.135, and three versions of our model, the BKd1, BKd2, and BKd3. The BKd1 is the original Gaussian model P. T. Kiss, M. Darvas, A. Baranyai, and P. Jedlovszky, J. Chem. Phys. 136, 114706 (2012)10.1063/1. 3692602 with constant polarization and with a simple exponential repulsion. The BKd2 applies field-dependent polarizability A. Baranyai and P. T. Kiss, J. Chem. Phys. 135, 234110 (2011)10.1063/1.3670962, while the BKd3 model has variable size to approximate the temperature-density (T-ρ) curve of water P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)10.1063/1.4746419. We used the thermodynamic integration (TI) and the Gibbs-Helmholtz equation to determine the equality of the free energy for liquid water and hexagonal ice (Ih) at 1 bar. We used the TIP4P and the SPCE models as reference systems of the TI. The studied models severely underestimate the experimental melting point of ice Ih. The calculated freezing points of the models are the following: COSG2, 215 K; COSG3, 149 K; SWM4-DP, 186 K; BKd1, 207 K; BKd2, 213 K; BKd3, 233 K. The freezing temperature of the SWM4-NDP system is certainly below 120 K. It might even be that the water phase is more stable than the ice Ih at 1 bar in the full temperature range. The COSB2 model melts below 100 K. The best result was obtained for the BKd3 model which indicates that correct description of the (T-ρ) curve improves the estimation of the freezing point. We also determined and compared the densities of high-pressure polymorphs of ice for these models.",
author = "Kiss, {P{\'e}ter T.} and P{\'e}ter Bertsyk and A. Baranyai",
year = "2012",
month = "11",
day = "21",
doi = "10.1063/1.4767063",
language = "English",
volume = "137",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "19",

}

TY - JOUR

T1 - Testing recent charge-on-spring type polarizable water models. I. Melting temperature and ice properties

AU - Kiss, Péter T.

AU - Bertsyk, Péter

AU - Baranyai, A.

PY - 2012/11/21

Y1 - 2012/11/21

N2 - We determined the freezing point of eight molecular models of water. All models use the charge-on-spring (COS) method to express polarization. The studied models were the COSG2, COSG3 H. Yu and W. F. van Gunsteren, J. Chem. Phys. 121, 9549 (2004)10.1063/1.1805516, the COSB2 H. Yu, T. Hansson, and W. F. van Gunsteren, J. Chem. Phys. 118, 221 (2003)10.1063/1.1523915, the SWM4-DP G. Lamoureux, A. D. MacKerell, Jr., and B. Roux, J. Chem. Phys. 119, 5185 (2003)10.1063/1.1598191, the SWM4-NDP G. Lamoureux, E. Harder, I. V. Vorobyov, B. Roux, and A. D. MacKerell, Jr., Chem. Phys. Lett. 418, 245 (2006)10.1016/j.cplett.2005.10.135, and three versions of our model, the BKd1, BKd2, and BKd3. The BKd1 is the original Gaussian model P. T. Kiss, M. Darvas, A. Baranyai, and P. Jedlovszky, J. Chem. Phys. 136, 114706 (2012)10.1063/1. 3692602 with constant polarization and with a simple exponential repulsion. The BKd2 applies field-dependent polarizability A. Baranyai and P. T. Kiss, J. Chem. Phys. 135, 234110 (2011)10.1063/1.3670962, while the BKd3 model has variable size to approximate the temperature-density (T-ρ) curve of water P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)10.1063/1.4746419. We used the thermodynamic integration (TI) and the Gibbs-Helmholtz equation to determine the equality of the free energy for liquid water and hexagonal ice (Ih) at 1 bar. We used the TIP4P and the SPCE models as reference systems of the TI. The studied models severely underestimate the experimental melting point of ice Ih. The calculated freezing points of the models are the following: COSG2, 215 K; COSG3, 149 K; SWM4-DP, 186 K; BKd1, 207 K; BKd2, 213 K; BKd3, 233 K. The freezing temperature of the SWM4-NDP system is certainly below 120 K. It might even be that the water phase is more stable than the ice Ih at 1 bar in the full temperature range. The COSB2 model melts below 100 K. The best result was obtained for the BKd3 model which indicates that correct description of the (T-ρ) curve improves the estimation of the freezing point. We also determined and compared the densities of high-pressure polymorphs of ice for these models.

AB - We determined the freezing point of eight molecular models of water. All models use the charge-on-spring (COS) method to express polarization. The studied models were the COSG2, COSG3 H. Yu and W. F. van Gunsteren, J. Chem. Phys. 121, 9549 (2004)10.1063/1.1805516, the COSB2 H. Yu, T. Hansson, and W. F. van Gunsteren, J. Chem. Phys. 118, 221 (2003)10.1063/1.1523915, the SWM4-DP G. Lamoureux, A. D. MacKerell, Jr., and B. Roux, J. Chem. Phys. 119, 5185 (2003)10.1063/1.1598191, the SWM4-NDP G. Lamoureux, E. Harder, I. V. Vorobyov, B. Roux, and A. D. MacKerell, Jr., Chem. Phys. Lett. 418, 245 (2006)10.1016/j.cplett.2005.10.135, and three versions of our model, the BKd1, BKd2, and BKd3. The BKd1 is the original Gaussian model P. T. Kiss, M. Darvas, A. Baranyai, and P. Jedlovszky, J. Chem. Phys. 136, 114706 (2012)10.1063/1. 3692602 with constant polarization and with a simple exponential repulsion. The BKd2 applies field-dependent polarizability A. Baranyai and P. T. Kiss, J. Chem. Phys. 135, 234110 (2011)10.1063/1.3670962, while the BKd3 model has variable size to approximate the temperature-density (T-ρ) curve of water P. T. Kiss and A. Baranyai, J. Chem. Phys. 137, 084506 (2012)10.1063/1.4746419. We used the thermodynamic integration (TI) and the Gibbs-Helmholtz equation to determine the equality of the free energy for liquid water and hexagonal ice (Ih) at 1 bar. We used the TIP4P and the SPCE models as reference systems of the TI. The studied models severely underestimate the experimental melting point of ice Ih. The calculated freezing points of the models are the following: COSG2, 215 K; COSG3, 149 K; SWM4-DP, 186 K; BKd1, 207 K; BKd2, 213 K; BKd3, 233 K. The freezing temperature of the SWM4-NDP system is certainly below 120 K. It might even be that the water phase is more stable than the ice Ih at 1 bar in the full temperature range. The COSB2 model melts below 100 K. The best result was obtained for the BKd3 model which indicates that correct description of the (T-ρ) curve improves the estimation of the freezing point. We also determined and compared the densities of high-pressure polymorphs of ice for these models.

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

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

U2 - 10.1063/1.4767063

DO - 10.1063/1.4767063

M3 - Article

VL - 137

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 19

M1 - 194102

ER -