Gradient-Driven Diffusion Using Dual Control Volume Grand Canonical Molecular Dynamics (DCV-GCMD)

Recently we developed a new nonequilibrium molecular simulation method [1] that allows the direct study of interdiffusion in multicomponent mixtures. The method combines stochastic insertion and deletion moves characteristic of grand canonical (GC) simulations with molecular dynamics (MD) to control the chemical potential μi of a species i. Restricting the insertions and deletions to two separate control volumes (CV's) one can apply different μ's in distinct locations, and thus create chemical potential gradients. DCV-GCMD can be used to study transient phenomena such as the filling of micropores or used in steady-state mode to determine the diffusion coefficients in multicomponent fluid mixtures. We report on the effects of molecular interactions and demonstrate how in a sufficiently nonideal ternary mixture this can lead to up-hill or reverse diffusion. In addition we introduce a novel extension of DCV-GCMD that is specifically designed for the study of gradient-driven diffusion of molecules that are simply too large to be inserted and deleted.