Computer simulations of liquid water: treatment of long-range interactions

Computer simulations by the molecular-dynamics method are used to study the physical properties of liquid water. Two three-point-charge models (SPC and TIPS) are analysed and compared using the Ewald-Kornfeld summation method to calculate long-range electrostatic interactions. Although these two models are not very different considering their geometry and energy parameters, they lead to different physical properties of liquid water. Simulations with the SPC water model and the Ewald-summation method are then used as a reference against which different termination functions for Coulombic interactions are compared. The largest discrepancies between simulations performed using the truncated potentials combined with termination functions and those performed using the non-truncated one occur in the radial distribution functions and in the diffusion coefficients. The origin of these discrepancies is analysed and it is shown that the termination function of the form 1 - 2(r/r c) n + (r/r c)2n with n = 1 provides...

[1]  L. Verlet Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[2]  Michiel Sprik,et al.  A polarizable model for water using distributed charge sites , 1988 .

[3]  Graham Hills,et al.  The computer simulation of polar liquids , 1979 .

[4]  M. Klein,et al.  A molecular dynamics study of the effect of pressure on the properties of water and ice , 1988 .

[5]  I. R. Mcdonald,et al.  Molecular dynamics studies of the structure of water at high temperatures and density , 1981 .

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

[7]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[8]  A. Barker Monte Carlo calculations of the radial distribution functions for a proton-electron plasma , 1965 .

[9]  F. Stillinger,et al.  Molecular Dynamics Study of Liquid Water , 1971 .

[10]  O. Matsuoka,et al.  CI study of the water dimer potential surface , 1976 .

[11]  D. Lévesque,et al.  A Monte Carlo study of dipolar hard spheres The pair distribution function and the dielectric constant , 1977 .

[12]  P. P. Ewald Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[13]  E. Clementi,et al.  Study of the structure of molecular complexes. XIII. Monte Carlo simulation of liquid water with a configuration interaction pair potential , 1976 .

[14]  B. Montgomery Pettitt,et al.  Integral equation predictions of liquid state structure for waterlike intermolecular potentials , 1982 .

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

[16]  B. Montgomery Pettitt,et al.  Structural and energetic effects of truncating long ranged interactions in ionic and polar fluids , 1985 .

[17]  Martin Karplus,et al.  SOLVATION. A MOLECULAR DYNAMICS STUDY OF A DIPEPTIDE IN WATER. , 1979 .

[18]  H. A. Levy,et al.  Liquid Water: Molecular Correlation Functions from X‐Ray Diffraction , 1971 .

[19]  J. A. Barker,et al.  Structure of water; A Monte Carlo calculation , 1969 .

[20]  L. V. Woodcock,et al.  Thermodynamic and structural properties of liquid ionic salts obtained by Monte Carlo computation. Part 1.—Potassium chloride , 1971 .

[21]  J. A. C. Rullmann,et al.  A polarizable water model for calculation of hydration energies , 1988 .

[22]  W. L. Jorgensen Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water , 2002 .

[23]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[24]  R. O. Watts,et al.  Monte Carlo studies of liquid water , 1974 .

[25]  J. E. Quinn,et al.  Cooperative effects in simulated water , 1979, Nature.

[26]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[27]  Bruce J. Berne,et al.  A Monte Carlo study of structural and thermodynamic properties of water: dependence on the system size and on the boundary conditions , 1980 .

[28]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

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

[30]  J. A. Barker,et al.  Monte Carlo studies of the dielectric properties of water-like models , 1973 .

[31]  F. Stillinger,et al.  Improved simulation of liquid water by molecular dynamics , 1974 .

[32]  E. Teller,et al.  Monte Carlo Study of a One‐Component Plasma. I , 1966 .

[33]  Roger Impey,et al.  Spectroscopic and transport properties of water , 1982 .