Relating the structure and dynamics of ionic liquids under shear by means of reverse non-equilibrium molecular dynamics simulations

The effect of the shear rate on the viscosity and the structure of 1-ethyl-3-methylimidazolium based ionic liquids with three different anions (tetrafluoroborate, dicyanamide, and bis(trifluoromethylsulfonyl)imide) was studied by means of reverse non-equilibrium molecular dynamics (RNEMD) simulations using a polarizable force field. The three liquids display a Newtonian plateau followed by a shear thinning regime at shear rates of the order of GHz. Even though the main features of the liquid structure remains under shear, systematic changes were noticed at the GHz rates, with coordination shells becoming more diffuse as noticed by the reduction in the difference between consecutive maxima and minima in the radial distribution function. Interestingly, these structural changes with the shear rate can be precisely fitted using the Carreau equation, which is a well-known expression for the shear rate dependence of the viscosity. The fitting parameters for different distributions can be used to explain qualitatively the shear thinning behavior of these liquids. In the GHz range, the cations and, in a minor extension, some anions, tend to assume preferentially a parallel orientation with the flux, which contributes to the shear thinning behavior and may have consequences for adhesion in applications as lubricants.