Calculating high-pressure PAO4 viscosity with equilibrium molecular dynamics simulations

The development of optimized lubricants is hindered by missing knowledge of fluid properties, in particular the viscosity, in the range of extreme pressures and temperatures relevant for application. Molecular dynamics simulations can be used to calculate viscosity, but the necessary computational effort imposes practical limits for high viscosities. In this study, the viscosity of PAO4 oil was extracted from equilibrium molecular dynamics simulations as a function of pressure and temperature reaching viscosities up to 20 Pas. Three calculation methods based on different microscopic expressions for the viscosity were used. The methods exhibit considerably different performance with respect to preciseness and computational efficiency. The highest viscosities were found to be calculated most efficiently via the Stokes–Einstein relation, by computing the diffusion coefficient from the velocity correlation function. This offers a new, more effective route to push viscosity calculations in equilibrium molecular dynamics simulations to higher pressure systems.