Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics

A method is proposed for calculating the shear viscosity of a liquid from finite-size effects of self-diffusion coefficients in Molecular Dynamics simulations. This method uses the difference in the self-diffusivities, computed from at least two system sizes, and an analytic equation to calculate the shear viscosity. To enable the efficient use of this method, a set of guidelines is developed. The most efficient number of system sizes is two and the large system is at least four times the small system. The number of independent simulations for each system size should be assigned in such a way that 50%–70% of the total available computational resources are allocated to the large system. We verified the method for 250 binary and 26 ternary Lennard-Jones systems, pure water, and an ionic liquid ([Bmim][Tf2N]). The computed shear viscosities are in good agreement with viscosities obtained from equilibrium Molecular Dynamics simulations for all liquid systems far from the critical point. Our results indicate that the proposed method is suitable for multicomponent mixtures and highly viscous liquids. This may enable the systematic screening of the viscosities of ionic liquids and deep eutectic solvents.

[1]  Akihito Otani,et al.  Reliable Viscosity Calculation from Equilibrium Molecular Dynamics Simulations: A Time Decomposition Method. , 2015, Journal of chemical theory and computation.

[2]  J. García,et al.  General friction theory viscosity model for the PC‐SAFT equation of state , 2006 .

[3]  Wolfgang Nowak,et al.  Bayesian Model Selection Helps To Choose Objectively between Thermodynamic Models: A Demonstration of Selecting a Viscosity Model Based on Entropy Scaling , 2016 .

[4]  Edward J Maginn,et al.  Direct Correlation between Ionic Liquid Transport Properties and Ion Pair Lifetimes: A Molecular Dynamics Study. , 2015, The journal of physical chemistry letters.

[5]  Kurt Kremer,et al.  Molecular dynamics simulation of a polymer chain in solution , 1993 .

[6]  T. Vlugt,et al.  Maxwell-Stefan diffusivities in binary mixtures of ionic liquids with dimethyl sulfoxide (DMSO) and H2O. , 2011, The journal of physical chemistry. B.

[7]  Gaurav Pranami,et al.  Estimating error in diffusion coefficients derived from molecular dynamics simulations. , 2015, Journal of chemical theory and computation.

[8]  C. Brooks.  Computer simulation of liquids , 1989 .

[9]  Emma L. Smith,et al.  Deep eutectic solvents (DESs) and their applications. , 2014, Chemical reviews.

[10]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[11]  K. Hermansson,et al.  Effective three-body potentials for Li+(aq) and Mg2+(aq) , 2003 .

[12]  A. Panagiotopoulos,et al.  Molecular simulation of thermodynamic and transport properties for the H2O+NaCl system. , 2014, The Journal of chemical physics.

[13]  Craig M. Tenney,et al.  Limitations and recommendations for the calculation of shear viscosity using reverse nonequilibrium molecular dynamics. , 2010, The Journal of chemical physics.

[14]  Ioannis G. Economou,et al.  Thermodynamic and Transport Properties of H2O + NaCl from Polarizable Force Fields. , 2015, Journal of chemical theory and computation.

[15]  Aage Fredenslund,et al.  An improved corresponding states model for the prediction of oil and gas viscosities and thermal conductivities , 1987 .

[16]  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.

[17]  K. Kobe.  The properties of gases and liquids , 1959 .

[18]  Y. Rosenfeld,et al.  A quasi-universal scaling law for atomic transport in simple fluids , 1999 .

[19]  B. Rotenberg,et al.  Diffusion in bulk liquids: finite-size effects in anisotropic systems , 2015 .

[20]  G. Hummer,et al.  Divergent Diffusion Coefficients in Simulations of Fluids and Lipid Membranes. , 2016, The journal of physical chemistry. B.

[21]  R. Zwanzig,et al.  Time-Correlation Functions and Transport Coefficients in Statistical Mechanics , 1965 .

[22]  M. Schoen,et al.  The shear viscosity of a Lennard-Jones fluid calculated by equilibrium molecular dynamics , 1985 .

[23]  Matthias Krack,et al.  Static and Dynamical Properties of Liquid Water from First Principles by a Novel Car-Parrinello-like Approach. , 2009, Journal of chemical theory and computation.

[24]  B. Rotenberg,et al.  Diffusion coefficient and shear viscosity of rigid water models , 2012, Journal of physics. Condensed matter : an Institute of Physics journal.

[25]  P. Tavan,et al.  Polarizable six-point water models from computational and empirical optimization. , 2014, The journal of physical chemistry. B.

[26]  G. Grest,et al.  Viscosity calculations of n-alkanes by equilibrium molecular dynamics , 1997 .

[27]  Stephan Kabelac,et al.  Transport coefficients of the Lennard-Jones model fluid. I. Viscosity. , 2004, The Journal of chemical physics.

[28]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[29]  Edward J Maginn,et al.  System-size corrections for self-diffusion coefficients calculated from molecular dynamics simulations: The case of CO2, n-alkanes, and poly(ethylene glycol) dimethyl ethers. , 2016, The Journal of chemical physics.

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

[31]  E. Maginn,et al.  Thermal and Transport Properties of Six Ionic Liquids: An Experimental and Molecular Dynamics Study , 2012 .

[32]  André Bardow,et al.  Finite-Size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics , 2018, Journal of chemical theory and computation.

[33]  D. Frenkel,et al.  Finite-size corrections to the chemical potential , 1992 .

[34]  Xiaojing Wang,et al.  A review on the transport properties of ionic liquids , 2014 .

[35]  C. Boned,et al.  Molecular Dynamics Study of the Lennard-Jones Fluid Viscosity: Application to Real Fluids , 2005 .

[36]  Randall Q. Snurr,et al.  A new perspective on the order-n algorithm for computing correlation functions , 2009 .

[37]  D. Bedeaux,et al.  Fick diffusion coefficients of liquid mixtures directly obtained from equilibrium molecular dynamics. , 2011, The journal of physical chemistry. B.

[38]  J. G. Powles,et al.  Self-diffusion coefficient of the hard-sphere fluid: system size dependence and empirical correlations. , 2007, The journal of physical chemistry. B.

[39]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[40]  Josefa Fernández,et al.  Scaling of the viscosity of the Lennard-Jones chain fluid model, argon, and some normal alkanes. , 2011, The Journal of chemical physics.

[41]  Amyn S. Teja,et al.  Generalized corresponding states method for the viscosities of liquid mixtures , 1981 .

[42]  J. Abascal,et al.  The shear viscosity of rigid water models. , 2010, The Journal of chemical physics.

[43]  José Mario Martínez,et al.  PACKMOL: A package for building initial configurations for molecular dynamics simulations , 2009, J. Comput. Chem..

[44]  Athanassios Z Panagiotopoulos,et al.  Atomistic Molecular Dynamics Simulations of Carbon Dioxide Diffusivity in n-Hexane, n-Decane, n-Hexadecane, Cyclohexane, and Squalane. , 2016, The journal of physical chemistry. B.

[45]  Hui Zhang,et al.  Corrections of the periodic boundary conditions with rectangular simulation boxes on the diffusion coefficient, general aspects , 2017 .

[46]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

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

[48]  André Bardow,et al.  Diffusion Coefficients from Molecular Dynamics Simulations in Binary and Ternary Mixtures , 2013 .

[49]  E. Maginn,et al.  Molecular Topology and Local Dynamics Govern the Viscosity of Imidazolium-Based Ionic Liquids. , 2015, The journal of physical chemistry. B.

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

[51]  F. Llovell,et al.  Free-volume theory coupled with soft-SAFT for viscosity calculations: comparison with molecular simulation and experimental data. , 2013, The journal of physical chemistry. B.

[52]  J. H. Dymond,et al.  Hard-sphere theories of transport properties , 1985 .

[53]  H. Shirota,et al.  Dicationic versus monocationic ionic liquids: distinctive ionic dynamics and dynamical heterogeneity. , 2013, The journal of physical chemistry. B.

[54]  K. Palanivelu,et al.  Carbon Dioxide Capture and Utilization by Alkanolamines in Deep Eutectic Solvent Medium , 2015 .

[55]  S. Passerini,et al.  From Nanoscale to Microscale: Crossover in the Diffusion Dynamics within Two Pyrrolidinium-Based Ionic Liquids. , 2017, The journal of physical chemistry letters.

[56]  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.

[57]  B. Rotenberg,et al.  Diffusion under Confinement: Hydrodynamic Finite-Size Effects in Simulation. , 2017, Journal of chemical theory and computation.

[58]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[59]  S. Aparicio,et al.  Systematic Study on the Viscosity of Ionic Liquids: Measurement and Prediction , 2015 .

[60]  Carlos Vega,et al.  Simulating water with rigid non-polarizable models: a general perspective. , 2011, Physical chemistry chemical physics : PCCP.

[61]  Sung Chul Bae,et al.  When Brownian diffusion is not Gaussian. , 2012, Nature Materials.

[62]  Thijs J. H. Vlugt,et al.  State-of-the-Art of CO2 Capture with Ionic Liquids , 2012 .

[63]  K. R. Harris.  Scaling the transport properties of molecular and ionic liquids , 2016 .

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

[65]  Berend Smit,et al.  Are pressure fluctuation-based equilibrium methods really worse than nonequilibrium methods for calculating viscosities? , 2009, The Journal of chemical physics.

[66]  Hans Hasse,et al.  Prediction of self-diffusion coefficient and shear viscosity of water and its binary mixtures with methanol and ethanol by molecular simulation. , 2011, The Journal of chemical physics.

[67]  Stephan Kabelac,et al.  Transport coefficients of the Lennard-Jones model fluid. II Self-diffusion. , 2004, The Journal of chemical physics.

[68]  Gerhard Hummer,et al.  System-Size Dependence of Diffusion Coefficients and Viscosities from Molecular Dynamics Simulations with Periodic Boundary Conditions , 2004 .

[69]  Rajamani Krishna,et al.  The Darken Relation for Multicomponent Diffusion in Liquid Mixtures of Linear Alkanes: An Investigation Using Molecular Dynamics (MD) Simulations , 2005 .

[70]  Sean Garrett-Roe,et al.  Ultrafast Structure and Dynamics in Ionic Liquids: 2D-IR Spectroscopy Probes the Molecular Origin of Viscosity. , 2014, The journal of physical chemistry letters.

[71]  Dick Bedeaux,et al.  Kirkwood-Buff Integrals for Finite Volumes. , 2013, The journal of physical chemistry letters.

[72]  F. Montel,et al.  Influence of the mass ratio on viscosity in Lennard–Jones mixtures: The one-fluid model revisited using nonequilibrium molecular dynamics , 2005 .

[73]  B. D. Todd,et al.  Transient-time correlation function applied to mixed shear and elongational flows. , 2012, The Journal of chemical physics.

[74]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[75]  Marcelle B. M. Spera,et al.  Self-diffusion coefficient of bulk and confined water: a critical review of classical molecular simulation studies , 2018, Molecular Simulation.

[76]  E. Maginn,et al.  Rapid shear viscosity calculation by momentum impulse relaxation molecular dynamics. , 2005, The Journal of chemical physics.

[77]  Howard J. M. Hanley,et al.  Prediction of transport properties. 1. Viscosity of fluids and mixtures , 1981 .

[78]  Jeffrey M. Young,et al.  System-Size Dependence of Electrolyte Activity Coefficients in Molecular Simulations. , 2018, The journal of physical chemistry. B.