Limitations and recommendations for the calculation of shear viscosity using reverse nonequilibrium molecular dynamics.

The reverse nonequilibrium molecular dynamics (RNEMD) method calculates the shear viscosity of a fluid by imposing a nonphysical exchange of momentum and measuring the resulting shear velocity gradient. In this study we investigate the range of momentum flux values over which RNEMD yields usable (linear) velocity gradients. We find that nonlinear velocity profiles result primarily from gradients in fluid temperature and density. The temperature gradient results from conversion of heat into bulk kinetic energy, which is transformed back into heat elsewhere via viscous heating. An expression is derived to predict the temperature profile resulting from a specified momentum flux for a given fluid and simulation cell. Although primarily bounded above, we also describe milder low-flux limitations. RNEMD results for a Lennard-Jones fluid agree with equilibrium molecular dynamics and conventional nonequilibrium molecular dynamics calculations at low shear, but RNEMD underpredicts viscosity relative to conventional NEMD at high shear.

[1]  B. Hess.  Determining the shear viscosity of model liquids from molecular dynamics simulations , 2002 .

[2]  O. Borodin,et al.  Viscosity of a room temperature ionic liquid: predictions from nonequilibrium and equilibrium molecular dynamics simulations. , 2009, Journal of Physical Chemistry B.

[3]  F. Müller-Plathe,et al.  Shear viscosity of the ionic liquid 1-n-butyl 3-methylimidazolium hexafluorophosphate [bmim][pf6] computed by reverse nonequilibrium molecular dynamics. , 2008, The journal of physical chemistry. B.

[4]  On the response of an ionic liquid to external perturbations and the calculation of shear viscosity. , 2007, The journal of physical chemistry. B.

[5]  F. Müller-Plathe.  A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity , 1997 .

[6]  Florian Müller-Plathe,et al.  The shear viscosity of molecular fluids: A calculation by reverse nonequilibrium molecular dynamics , 2002 .

[7]  Palmer.  Transverse-current autocorrelation-function calculations of the shear viscosity for molecular liquids. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[8]  Christopher P Calderon,et al.  Comment on "reversing the perturbation in nonequilibrium molecular dynamics: an easy way to calculate the shear viscosity of fluids". , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  F Müller-Plathe,et al.  Reversing the perturbation in nonequilibrium molecular dynamics: an easy way to calculate the shear viscosity of fluids. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[10]  Gary P. Morriss,et al.  Statistical Mechanics of Nonequilibrium Liquids , 2008 .

[11]  E. Maginn,et al.  Effect of temperature and water content on the shear viscosity of the ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide as studied by atomistic simulations. , 2007, The journal of physical chemistry. B.

[12]  Hsueh-Chia Chang,et al.  Efficient viscosity estimation from molecular dynamics simulation via momentum impulse relaxation , 2000 .

[13]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[14]  E. Maginn,et al.  Prediction of viscosities and vapor–liquid equilibria for five polyhydric alcohols by molecular simulation , 2007 .

[15]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[16]  P. Daivis,et al.  Comparison of constant pressure and constant volume nonequilibrium simulations of sheared model decane , 1994 .

[17]  H. H. Chiang,et al.  Shear‐rate‐dependence modeling of polymer melt viscosity , 1992 .