Calculation of the mutual diffusion coefficient by equilibrium and nonequilibrium molecular dynamics

A nonequilibrium molecular dynamics method for the calculation of the mutual diffusion coefficient for a mixture of hard spheres is described. The method is applied to a 50-50 mixture of equidiameter particles having a mass ratio of 0.1 for the two species, at a volume of three times close-packing. By extrapolating the results to the limit of vanishing concentration gradient and infinite system size, we obtain a value in statistical agreement with the result obtained using a Green-Kubo molecular dynamics procedure which is also described. The nonequilibrium calculation yields a mutual diffusion coefficient which decreases slightly with increasing concentration gradient. The Green-Kubo time correlation function for mutual diffusion displays a slow decay with time, qualitatively similar to the long-time tail which has been predicted by the hydrodynamic theory of Pomeau