On the force bias Monte Carlo simulation of water: methodology, optimization and comparison with molecular dynamics

Abstract

An optimized version of the force bias scheme is presented for the Monte Carlo simulation of water in which the random walk of the individual molecule is biased in the direction of the force and torque acting on the molecule. A new criterion is developed to judge the efficiency of sampling the configuration space by studying the translational and rotational diffusion of individual molecules during the random walk. The force bias method is compared with the uniform sampling method using ST2 water as an example, and with molecular dynamics. Similarities and differences between molecular dynamics and Monte Carlo are discussed with respect to pair correlations and energy distributions.

Additional information

Notes on contributors

C. Pangali

In partial fulfilment of the Ph.D. in Chemistry at Columbia University.