Viscosity of a room temperature ionic liquid: predictions from nonequilibrium and equilibrium molecular dynamics simulations.
Nonequilibrium molecular dynamics (NEMD) simulations have been performed on 1-methyl-3-ethyl-imidazolium bis(trifluoromethane)sulfonimide [emim][Ntf(2)] using Lees-Edwards boundary conditions. A range of inverse shear rates corresponding to a fraction of the rotational relaxation time for the slowest relaxing molecular axis of anion and cation to 20 rotational relaxation times (1/20 tau(rot) < gamma < 5/tau(rot)) has been investigated. An extrapolation of the shear-rate-dependent viscosity obtained from these simulations to zero shear rate using the empirical three-parameter Carreau equation yielded excellent agreement with the viscosity obtained from equilibrium MD simulations. Based upon the Carreau equation fit to the simulation data, shear-thinning behavior was observed in [emim][Ntf(2)] for all shear rates investigated, implying that Newtonian behavior is observed in [emim][Ntf(2)] only for shear rates significantly lower than the inverse rotational relaxation time. A close resemblance between the apparent time-dependent viscosity extracted from equilibrium MD simulations and the shear-rate-dependent viscosity extracted from NEMD simulations has been found and discussed. MD simulations accurately predicted [emim][Ntf(2)] density, self-diffusion coefficients, heat of vaporization, and lattice parameters for the crystalline phase.
Borodin
Oleg
O
Department of Materials Science & Engineering, 122 South Central Campus Drive, Room 304, University of Utah, Salt Lake City, Utah 84112-0560, USA. Oleg.Borodin@utah.edu
Smith
Grant D
GD
Kim
Hojin
H
D003198:Computer Simulation; D052578:Ionic Liquids / Q000737:chemistry*; D008956:Models, Chemical; D011789:Quantum Theory; D012399:Rotation; D013696:Temperature; D014783:Viscosity
D052578:Ionic Liquids