Mass transport coefficients play an important role in process design and in compositional grading of oil reservoirs. As experimental measurements of these properties can be costly and hazardous, Molecular Dynamics simulations emerge as an alternative approach. In this work, we used Molecular Dynamics to calculate the self-diffusion coefficients of methane/n-hexane mixtures at different conditions, in both liquid and supercritical phases. We evaluated how the finite box size and the choice of the force field affect the calculated properties at high pressures. Results show a strong dependency between self-diffusion and the simulation box size. The Yeh-Hummer analytical correction [J. Phys. Chem. B, 108, 15873 (2004)] can attenuate this effect, but sometimes makes the results depart from experimental data due to issues concerning the force fields. We have also found that different all-atom and united-atom models can produce biased results due to caging effects and to different dihedral configurations of the n-alkane.
Optimal estimates of self-diffusion coefficients from molecular dynamics simulations
