A gentler approach to RNEMD: nonisotropic velocity scaling for computing thermal conductivity and shear viscosity

We present a new method for introducing stable nonequilibrium velocity and temperature gradients in molecular dynamics simulations of heterogeneous systems. This method extends earlier reverse nonequilibrium molecular dynamics (RNEMD) methods which use momentum exchange swapping moves. The standard swapping moves can create nonthermal velocity distributions and are difficult to use for interfacial calculations. By using nonisotropic velocity scaling (NIVS) on the molecules in specific regions of a system, it is possible to impose momentum or thermal flux between regions of a simulation while conserving the linear momentum and total energy of the system. To test the method, we have computed the thermal conductivity of model liquid and solid systems as well as the interfacial thermal conductivity of a metal-water interface. We find that the NIVS-RNEMD improves the problematic velocity distributions that develop in other RNEMD methods.