Mw. Mahoney et Wl. Jorgensen, Rapid estimation of electronic degrees of freedom in Monte Carlo calculations for polarizable models of liquid water, J CHEM PHYS, 114(21), 2001, pp. 9337-9349
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.