Molecular simulation of chemical reaction equilibria by Kinetic Monte Carlo

ABSTRACT

We describe a new algorithm for the molecular simulation of chemical reaction equilibria, which we call the Reactive Kinetic Monte Carlo (ReKMC) algorithm. It is based on the use of the equilibrium Kinetic Monte Carlo (eKMC) method (Ustinov et al., J. Colloid Interface Sci., 2012, 366, 216–223) to generate configurations in the underlying nonreacting system and to calculate the species chemical potentials at essentially zero marginal computational cost. We consider in detail the typical case of specified temperature, T and pressure, P, but extensions to other thermodynamic constraints are straightforward in principle. In the course of this work, we also demonstrate an alternative method for calculating simulation box volume changes in NPT ensemble simulations to achieve the specified P. We consider two sets of example reacting systems previously considered in the literature, and compare the ReKMC results and computational efficiencies with those of different implementations of the REMC algorithm (Turner et al., Molec. Simulation, 2008, 34, 119–146).