Rapid estimation of electronic degrees of freedom in Monte Carlo calculations for polarizable models of liquid water

The inclusion of electronic polarization within Monte Carlo calculations of simple models of molecular liquids is hampered, relative to its inclusion within molecular dynamics calculations, by the need to fully determine the variables that specify the electronic configuration every time each molecul e is moved, i.e., N times per cycle, rather than once per cycle. Classical statistical mechanical Monte Carlo calculations on two models of liquid wat er have been performed. For each of the models, electronic degrees of freed om are modeled by polarizable sites; thus it is the components of the induc ed dipole vector that must be determined at each step. Commonly used approx imation methods have been characterized and found to be inadequate. Efficie nt procedures have been devised to estimate the dipole vector and have been tested on reproducing electronic, thermodynamic, and structural properties of the two polarizable water models. The most promising procedure, conside ring both computational time saved and accuracy at reproducing pure liquid properties, involves approximating the induced dipoles at each step by an i nitial perturbative modification of the dipoles from the previous step, fol lowed by an iteration of the induced dipoles on only the moved molecule. Wi th this procedure, the CPU time is dramatically reduced, and the thermodyna mic and structural properties are estimated correctly to within a few perce nt. They are reproduced more rapidly and with greater accuracy than in calc ulations in which the dipole vector is estimated by a single iterative cycl e starting with the dipoles from the previous Monte Carlo step. (C) 2001 Am erican Institute of Physics.