Molecular dynamics simulation of the density and surface tension of water by particle-particle particle-mesh method.

In this work, molecular dynamics simulation is performed to study the density and surface tension of water for a range of temperatures from 300 to 600 K. The extended simple point charge interaction potential for water is used. The particle-particle particle-mesh method, which automatically includes untruncated long-range terms, is used for the Lennard-Jones and the Coulombic terms. The results show that the long-range correction for the Lennard-Jones term is very important for the calculation of surface tension. It is found that the calculated density and surface tension of water fit well with experimental data for temperatures less than 500 K. Near the critical temperature, the simulation results are off from the experimental data.

[1]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[2]  A. Pohorille,et al.  Molecular dynamics of the water liquid-vapor interface. , 1987, The Journal of physical chemistry.

[3]  Peter T. Cummings,et al.  Molecular simulation of water along the liquid–vapor coexistence curve from 25 °C to the critical point , 1990 .

[4]  M. Matsumoto,et al.  Study on liquid–vapor interface of water. I. Simulational results of thermodynamic properties and orientational structure , 1988 .

[5]  D. Tildesley,et al.  Molecular dynamics simulation of the orthobaric densities and surface tension of water , 1995 .

[6]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[7]  J. D. Bernal,et al.  A Theory of Water and Ionic Solution, with Particular Reference to Hydrogen and Hydroxyl Ions , 1933 .

[8]  B. Guillot,et al.  A computer simulation study of the liquid–vapor coexistence curve of water , 1993 .

[9]  W. Reynolds.  Thermodynamic properties in SI , 1979 .

[10]  M. Deserno,et al.  HOW TO MESH UP EWALD SUMS. II. AN ACCURATE ERROR ESTIMATE FOR THE PARTICLE-PARTICLE-PARTICLE-MESH ALGORITHM , 1998, cond-mat/9807100.

[11]  Bertrand Guillot,et al.  A reappraisal of what we have learnt during three decades of computer simulations on water , 2002 .

[12]  A. D. McLean,et al.  Monte Carlo simulation of the liquid-vapor interface of water using an ab initio potential , 1993 .

[13]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[14]  A. Laaksonen,et al.  Surface tension of water droplets: A molecular dynamics study of model and size dependencies , 1997 .

[15]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[16]  G. W. Robinson,et al.  Flexible simple point-charge water in a self-supporting thin film , 1991 .

[17]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[18]  Christian Holm,et al.  How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines , 1998 .

[19]  W. V. van Gunsteren,et al.  A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations , 2001 .