Parametric dependence on shear viscosity of SPC/E water from equilibrium and non-equilibrium molecular dynamics simulations

Abstract

A parametric dependent study is crucial for the accurate determination of transport coefficients such as shear viscosity. In this study, we calculate the shear viscosity of extended simple point charge water using a transverse current auto-correlation function (TCAF) from equilibrium molecular dynamics (EMD) and the periodic perturbation method from non-equilibrium molecular dynamics (NEMD) simulations for varying coupling time and system sizes. Results show that the shear viscosity calculated using EMD simulations with different thermostats varies significantly with coupling times and system size. The use of Berendsen and velocity-rescale thermostats in NEMD simulations generates a significant drift from the target temperature and results in an inconsistent shear viscosity with coupling time and system size. The use of Nosé–Hoover thermostat in NEMD simulations offers thermodynamic stability which results in a consistent shear viscosity for various coupling times and system sizes.