Molecular simulation of the shear viscosity and the self-diffusion coefficient of mercury along the vapor-liquid coexistence curve

In earlier work [G. Raabe and R. J. Sadus, J. Chem. Phys. 119, 6691 (2003)] we reported that the combination of an accurate two-body ab initio potential with an empirically determined multibody contribution enables the prediction of the phase coexistence properties, the heats of vaporization, and the pair distribution functions of mercury with reasonable accuracy. In this work we present molecular dynamics simulation results for the shear viscosity and self-diffusion coefficient of mercury along the vapor-liquid coexistence curve using our empirical effective potential. The comparison with experiment and calculations based on a modified Enskog theory shows that our multibody contribution yields reliable predictions of the self-diffusion coefficient at all densities. Good results are also obtained for the shear viscosity of mercury at low to moderate densities. Increasing deviations between the simulation and experimental viscosity data at high densities suggest that not only a temperature-dependent but also a density-dependent multibody contribution is necessary to account for the effect of intermolecular interactions in liquid metals. An analysis of our simulation data near the critical point yields a critical exponent of β=0.39, which is identical to the value obtained from the analysis of the experimental saturation densities.