Convergence of Monte Carlo simulations of liquid water in the NPT ensemble

Abstract Long Monte Carlo simulations of liquid water at 25° and 1 atm have been carried out to study the convergence characteristics of the calculations. The recently reported TIPS2 potential was employed with system sizes of 125 and 216 monomers. The enthalpy, volume and radial distribution functions converge rapidly and show little size dependence. However, the rates of convergence are much slower for the fluctuation properties and are in the order: heat capacity ( Cp ) ⪢ isothermal compressibility (ϰ) ⪢ coefficient of thermal expansion (α). In fact,the weak coupling of the enthalpy and volume allows only crude estimates for α. In addition, the estimation of statistical error bounds for the thermodynamic properties is analyzed. It is found that very long (2500-3500 k) simulations are needed to yield valid estimates of the errors for the enthalpy, volume and C p , while the erros for ϰ and α would still be elusive.

[1]  Mihaly Mezei,et al.  Convergence characteristics of Monte Carlo–Metropolis computer simulations on liquid water , 1979 .

[2]  W. S. Benedict,et al.  Rotation‐Vibration Spectra of Deuterated Water Vapor , 1956 .

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

[4]  Monte Carlo calculations in the isothermal-isobaric ensemble. 1. Liquid water , 1977 .

[5]  W. L. Jorgensen.  Transferable intermolecular potential functions. Application to liquid methanol including internal rotation , 1981 .

[6]  William L. Jorgensen.  Quantum and statistical mechanical studies of liquids. 12. Simulation of liquid ethanol including internal rotation , 1981 .

[7]  Structures and properties of organic liquids: n-alkyl ethers and their conformational equilibriums , 1981 .

[8]  W. L. Jorgensen,et al.  Structures and properties of organic liquids: n-butane and 1,2-dichloroethane and their conformation equilibriums , 1981 .

[9]  W. L. Jorgensen.  Pressure dependence of the structure and properties of liquid n-butane , 1981 .

[10]  W. L. Jorgensen,et al.  Pressure dependence of hydrogen bonding in liquid methanol , 1982 .

[11]  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 .

[12]  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 .