Shear stress relaxation and diffusion in simple liquids by molecular dynamics simulations: Analytic expressions and...

The results are reported of an equilibrium molecular dynamics simulation study of the shear viscosity, η, and self-diffusion coefficient, D, of the Lennard-Jones liquid using the Green-Kubo (GK) method. Semiempirical analytic expressions for both GK time correlation functions were fitted to the simulation data and used to derive analytic expressions for the time dependent diffusion coefficient and shear viscosity, and also the correlation function frequency transforms. In the case of the shear viscosity for a state point near the triple point, a sech function was found to fit the correlation function significantly better than a gaussian in the ballistic short time regime. A reformulation of the shear GK formula in terms of a time series of time integrals ("viscuits") and contributions to the viscosity from components based on the initial stress ("visclets") enable the GK expressions to be recast in terms of probability distributions which could be used in coarse grained stochastic models of nanoscale flow. The visclet treatment shows that stress relaxation is statistically independent of the initial stress for equilibrium and metastable liquids, suggesting that shear stress relaxation in liquids is diffusion controlled. By contrast, the velocity autocorrelation function is sensitive to the initial velocity. Weak oscillations and a plateau at intermediate times originate to a greater extent from the high velocity tail of the Maxwell-Boltzmann velocity distribution. Simple approximate analytic expressions for the mean square displacement and the self Van Hove correlation function are also derived.