Solvent-shift Monte Carlo: a cluster algorithm for solvated systems

Abstract

We present a cluster algorithm for the efficient simulation of solvated systems that we term solvent-shift Monte Carlo (SSMC). The algorithm involves a conformational change in a solvated solute molecule of interest, followed by a concerted movement of solvent particles about a rotation axis. The method satisfies detailed balance and can be applied to existing schemes to sample conformational space, where an axis or plane of rotation can be defined. We demonstrate that the algorithm significantly enhances the sampling of phase space in solvated systems, and may be easily combined with other advanced sampling techniques.