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

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 Cp , while the erros for ϰ and α would still be elusive.