Improvement in molecule exchange efficiency in Gibbs ensemble Monte Carlo: development and implementation of the continuous fractional component move.
The continuous fractional component Monte Carlo (CFC MC) move (J Chem Theory Comput, 2007, 3, 1451) is extended to the Gibbs ensemble. The algorithm is validated against conventional simulations for the Lennard Jones fluid and a flexible water model. The method is also used to compute the vapor-liquid coexistence densities of a model for SO(2). The CFC molecule exchange move relies on the gradual insertion and deletion of molecules in conjunction with a self-adapting bias potential. As a result, the method does not require the formation of spontaneous voids in the dense fluid phase to be successful, leading to molecule exchange acceptance probabilities that are nearly independent of temperature. For example, over 1% of the vapor-liquid molecule exchange moves are successful for water at 280 K, whereas advanced rotational and configurational bias methods have success rates of less than 0.03%. The CFC move can be combined with other Monte Carlo moves to enable efficient simulation of dense strongly associating fluids that are to this point infeasible to model with standard methods.
Shi
Wei
W
Department of Chemical and Biomolecular Engineering, University of Notre Dame, 182 Fitzpatrick Hall, Notre Dame, IN 46556-5637, USA.
Maginn
Edward J
EJ