[165] Influence of simulation protocols on the efficiency of Gibbs ensemble Monte Carlo simulations

A.D. Cortes-Morales, I.G. Econonmou, C.J. Peters, and J.I. Siepmann

Molec. Simul. 39, 1135-1142 (2013)

Publication Abstract

The Gibbs ensemble Monte Carlo (GEMC) method is a versatile approach for the prediction of fluid phase equilibria from particle-based simulations. For a one-component system, a GEMC simulation utilises two separate simulation boxes for the vapour and the liquid phases and a significant fraction of the computational effort is expended on special trial moves that transfer (swap) particles and exchange volume between the two boxes. The user needs to specify the frequency of swap and volume moves and the overall volume that controls the phase ratio. In this study, the efficiency of GEMC simulation protocols that yield three different frequencies of accepted swap and volume moves and three different phase ratios is assessed for the computation of the saturated vapour pressure and liquid density of n-octane and water at three reduced temperatures. Differences in the simulation efficiency of up to an order of magnitude are observed, and recommendations are made for suitable GEMC simulation protocols.

Influence of simulation protocols on the efficiency of Gibbs ensemble Monte Carlo simulations