DIFFUSION IN LENNARD-JONES FLUIDS USING DUAL CONTROL-VOLUME GRAND-CANONICAL MOLECULAR-DYNAMICS SIMULATION (DCV-GCMD)

A new approach to calculating diffusivities, both transport as well as equilibrium, is presented. The dual control volume grand canonical mo lecular dynamics (or DCV-GCMD) method employs two local control volume s for chemical potential control via particle creation/destruction as in grand canonical Monte Carlo (GCMC) simulations. The control volumes are inserted in a standard NVT molecular dynamics simulation yielding a simulation with stochastic chemical potential control that may be t hought of as a hybrid GCMC-MD approach. The geometrical control of the chemical potential enables a steady state chemical potential gradient to be established in the system. By measuring the density profile and flux, Fick's law is used to determine the diffusivity. An example cal culation is presented for a simple Lennard-Jones system.