Adiabatic Nuclear and Electronic Sampling Monte Carlo Simulations in the Gibbs Ensemble: Application to Polarizable Force Fields for Water

The adiabatic nuclear and electronic sampling Monte Carlo algorithm (ANES-MC) is extended to simulations in the Gibbs ensemble. Whereas the maximum displacements used for translational, rotational, and volume trial moves can be adjusted to foster efficient sampling in the adiabatic limit, the transfer (swap) of particles always causes a major disturbance of the electronic structures of the two phases (supplying and receiving the particle). To reequilibrate the electronic structures requires additional sampling of the electronic degrees of freedom. A simple, distance-dependent criterion for the preferential selection of the electronic degrees of freedom, for which a move is to be attempted, is shown to improve the efficiency of the particle swap move. The ANES-MC algorithm is applied to the polarizable simple point charge-fluctuating charge (SPC−FQ) and transferable intermolecular potential 4 point-fluctuating charge (TIP4P−FQ) models proposed by Rick et al. (J. Chem. Phys. 1994, 101, 6141). For both model...