The theory of the correction for the difference in the fluctuations among the different ensembles was applied to derive the time correlation function formula required to compute the distinct diffusion coefficients in the barycentric reference frame by molecular dynamics. In the process it was found that the dynamical variables in the time correlation function of any mass transport coefficient of any reference frame were the same in the canonical as in the molecular dynamics ensemble. The result agreed with the accepted formulas for computing the kinetic part of the mutual diffusion coefficient in a binary mixture of non-electrolytes and for computing the electrical conductivity in molten salts and electrolyte solutions. The static velocity correlations in the molecular dynamics ensemble were analyzed without recourse to the other ensembles. In sharp contrast with the ensembles which were not constrained with respect to the total momentum, the laboratory-frame velocities of two different particles at the same time were found to be correlated. The result obtained was equivalent to the initial value of the canonical time correlation function of the barycentric distinct diffusion coefficients derived in the preceding paper.

Velocity Correlations in the Molecular Dynamics Ensemble: Computation of the Distinct Diffusion Coefficients. Raineri, F.O., Friedman, H.L.: The Journal of Chemical Physics, 1989, 91[9], 5642-7