Journal of Chemical Physics
The Journal of Chemical Physics
   
 
 
 
Previous Article
Direct calculation of solid-liquid coexistence points of a binary mixture by thermodynamic integration
We present a new thermodynamic integration method that directly connects the liquid and the solid phases of a binary mixture by a reversible path. The states along the path are simulated in the isothe...
Next Article
Solvent phase behavior and the interaction of uniform and patterned solutes
Isotropic and anisotropic hypernetted-chain (HNC) integral equation theories are used to obtain the interaction of solutes both near and far from the solvent liquid-vapor coexistence. Spherically symm...

Characterization of the TIP4P-Ew water model: Vapor pressure and boiling point

J. Chem. Phys. 123, 194504 (2005); doi:10.1063/1.2085031

Published 14 November 2005

You are not logged in to this journal. Log in

Hans W. Horn, William C. Swope, and Jed W. Pitera
International Business Machines Corporation (IBM) Almaden Research Center, 650 Harry Road, San Jose, California 95120
The liquid-vapor-phase equilibrium properties of the previously developed TIP4P-Ew water model have been studied using thermodynamic integration free-energy simulation techniques in the temperature range of 274–400  K. We stress that free-energy results from simulations need to be corrected in order to be compared to the experiment. This is due to the fact that the thermodynamic end states accessible through simulations correspond to fictitious substances (classical rigid liquids and classical rigid ideal gases) while experiments operate on real substances (liquids and real gases, with quantum effects). After applying analytical corrections the vapor pressure curve obtained from simulated free-energy changes is in excellent agreement with the experimental vapor pressure curve. The boiling point of TIP4P-Ew water under ambient pressure is found to be at 370.3±1.9  K, about 7  K higher than the boiling point of TIP4P water (363.7±5.1  K; from simulations that employ finite range treatment of electrostatic and Lennard-Jones interactions). This is in contrast to the approximately +15  K by which the temperature of the density maximum and the melting temperature of TIP4P-Ew are shifted relative to TIP4P, indicating that the temperature range over which the liquid phase of TIP4P-Ew is stable is narrower than that of TIP4P and resembles more that of real water. The quality of the vapor pressure results highlights the success of TIP4P-Ew in describing the energetic and entropic aspects of intermolecular interactions in liquid water. ©2005 American Institute of Physics
History: Received 1 July 2005; accepted 26 August 2005; published 14 November 2005
Permalink: http://link.aip.org/link/?JCPSA6/123/194504/1
BUY THIS ARTICLE   (US$28)
Download HTML Download Sectioned HTML Download PDF (412 kB) View Cart

KEYWORDS and PACS

Keywords
PACS
  • 65.20.+w
    Thermal properties of liquids: heat capacity, thermal expansion, etc
  • 64.70.Fx
    Liquid–vapor transitions
  • YEAR: 2005

RELATED DATABASES


To view database links for this article,
you need to log in.
To view database links for this article,
you need to log in.

PUBLICATION DATA

ISSN:
0021-9606 (print)   1089-7690 (online)
Publisher:
AIP is a member of CrossRef AIP

REFERENCES (36)

For access to fully linked references, you need to log in. For access to fully linked references, you need to Log in.
  1. H. W. Horn, W. C. Swope, J. W. Pitera, J. D. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon, J. Chem. Phys. 120, 9665 (2004).
  2. W. L. Jorgensen, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983).
  3. H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. 91, 6269 (1987).
  4. J.P. M. Postma, Ph.D. thesis, University of Groningen (1985).
  5. W. L. Jorgensen and J. D. Madura, Mol. Phys. 56, 1381 (1985).
  6. M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 115, 10758 (2001).
  7. C. Vega, E. Sanz, and J. L. F. Abascal, J. Chem. Phys. 122, 114507 (2005).
  8. M.S. Gulam Razul (PhD thesis, Dalhousie University, Halifax, NS (Canada), Aug. 2005).
  9. T. Dick (PhD thesis, Duquesne University, Duquesne, PA, 2006).
  10. J. G. Kirkwood, J. Chem. Phys. 3, 300 (1935).
  11. J.E. Rice, H. Horn, B.H. Lengsfield III et al., MULLIKEN, a computational chemistry program, Version 2.0 (Almaden Research Center, IBM Research Division, San Jose, CA, 1996).
  12. W. C. Swope, H. C. Andersen, H. P. Berens, and K. R. Wilson, J. Chem. Phys. 76, 637 (1982).
  13. V. Kraeutler, W. F. van Gunsteren, and P. H. Huenenberger, J. Comput. Chem. 22, 501 (2001).
  14. H. C. Andersen, J. Chem. Phys. 72, 2384 (1980).
  15. amu=atomic mass unit.
  16. M. Born and T. V. Karman, Phys. Z. 13, 297 (1912).
  17. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, J. Chem. Phys. 21, 1087 (1953).
  18. L. Verlet, Phys. Rev. 159, 98 (1967).
  19. P. P. Ewald, Ann. Phys. 64, 253 (1921).
  20. T. C. Beutler, A. E. Mark, R. C. van Schaik, P. R. Gerber, and W. F. van Gunsteren, Chem. Phys. Lett. 222, 529 (1994).
  21. A. Delle Site, J. Phys. Chem. Ref. Data 26, 157 (1997).
  22. M. R. Shirts, J. W. Pitera, W. C. Swope, and V. S. Pande, J. Chem. Phys. 119, 5740 (2003).
  23. J.R. Taylor, An Introduction to Error Analysis, 2nd ed. (University Science Books, Sausalito, 1997), Chap. 8, pp. 201, 204.
  24. W. Wagner and A. Pruss, J. Phys. Chem. Ref. Data 31, 387 (2003).
  25. O. A. Karim and A. D. J. Haymet, J. Chem. Phys. 89, 6889 (1988).
  26. See also the detailed discussion on this in the appendix of Ref. 22.
  27. H. R. Zelsmann, J. Mol. Struct. 350, 95 (1995).
  28. P. F. Bernath, Phys. Chem. Chem. Phys. 4, 1501 (2002).
  29. J. K. Vij, D. R. J. Simpson, and O. E. Panarina, J. Mol. Liq. 112, 125 (2003).
  30. G. Seeley and T. Keyes, J. Chem. Phys. 91, 5581 (1989).
  31. F. Sciortino and S. Sastry, J. Chem. Phys. 100, 3881 (1994).
  32. H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. 91, 6269 (1987).
  33. The experimental values used for µgas and alphagas are: µgas=1.854989  D and alphagas=1.6633  10–40  F  m2, as recommended by IAPWS;
  34. Guideline on the Use of Fundamental Physical Constants and Basic Constants of Water, http://www.iapws.org/relguide/fundam.pdf (2001).
  35. Fundamental vibrational frequencies nu  (cm–1); from IR spectra of water gas and liquid from Ref. 4 that have been used for the present study: [bullet] intramolecular(gas):         {3755.7, 3656.6, 1594.6};[bullet] intramolecular(liquid):      {3490, 3280, 1645};[bullet] intermolecular(liquid):      {800, 500, 200, 50}.More recent experiments (Refs. 27,28,29) have yielded frequencies of the following:[bullet] intramolecular(gas):         {3755.93, 3657.05, 1594.75};[bullet] intramolecular(liquid):      {3490, 3277, 1645};[bullet] intermolecular(liquid):      {683.3, 395.5, 183.4, 50}.Using these "newer" frequencies for correcting DeltaGvap would cause the following shifts in the properties of the models studied: (a) Tboil by +4  K, (b) DeltaHvap by +0.2  kcal/mol, and (c) DeltaSvap by +0.3  cal/mol. For the sake of consistency, we have chosen to use the "older" frequencies from Ref. 4 because they were also used in the parametrization in the TIP4P-Ew model (see Ref. 1). A comprehensive, excellent (and well-maintained) source of theoretical and experimental information about water on the web is http://www.lsbu.ac.uk/water/index.html by M. Chaplin.
  36. M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys. 112, 8910 (2000).
  37. H.J. C. Berendsen, J.P. M. Postma, W.F. van Gunsteren, and J. Hermans, in Intermolecular Forces, edited by B. Pullmann (Reidel, Dordrecht, 1981), p. 331.

CITING ARTICLES

For access to citing articles, you need to log in.
For access to citing articles, you need to Log in.