Shear Viscosity of the Ionic Liquid 1-n-Butyl 3-Methylimidazolium Hexafluorophosphate [bmim][PF6] Computed by Reverse Nonequilibrium Molecular Dynamics

Confined ionic liquids films under shear: The importance of the chemical nature of the solid surface

Ionic liquids have generated interest in applications as lubricants and as additives to conventional lubricants due to their unique physical properties. In these applications, the liquid thin film can be subjected simultaneously to extremely high shear and loads in addition to nanoconfinement effects. Here, we use molecular dynamics simulations with a coarse-grained model to study a nanometric film of an ionic liquid confined between two planar solid surfaces both at equilibrium and at several shear rates. The strength of the interaction between the solid surface and the ions was changed by simulating three different surfaces with enhanced interactions with different ions. The increase in the interaction with either the cation or the anion leads to the formation of a solid-like layer that moves alongside the substrates; however, this layer can exhibit different structures and stability. An increase in the interaction with the high symmetry anion produces a more regular structure that is more resistant to the effects of shear and viscous heating. Two definitions were proposed and used for the calculation of the viscosity: a local definition based on the microscopic characteristics of the liquid and an engineering definition based on the forces measured at the solid surfaces, with the former displaying a correlation with the layered structure induced by the surfaces. Because of the shear thinning behavior of the ionic liquids as well as the temperature rise brought on by viscous heating, both the engineering and the local viscosities decrease as the shear rate increases.

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.

Nanoconfinement between Graphene Walls Suppresses the Near-Wall Diffusion of the Ionic Liquid [BMIM][PF6]

We identify two distinct regimes for the diffusion of the ionic liquid [BMIM][PF₆] confined between parallel graphene walls using molecular dynamics simulations. Within 2 nm of the wall, the cations and anions form a well-defined layered structure. In this region, the in-plane diffusion coefficients are suppressed when compared to their bulk values and increase monotonically with the distance away from the wall. Beyond 2 nm from the wall, the density profile and in-plane diffusion coefficients recover their bulk values. The channel-averaged in-plane diffusion coefficients increase monotonically with wall separation and recover the bulk values at a separation of 15 nm. A simple semianalytical model is proposed that mirrors this trend. The results also highlight the importance of applying a finite-size correction to molecular dynamics-predicted diffusion coefficients of confined liquids, which may otherwise be unusually larger than their bulk values.

Atomistic reverse nonequilibrium molecular dynamics simulation of the viscosity of ionic liquid 1-n-butyl 3-methylimidazolium bis(trifluoromethylsulfonyl)imide [bmim][Tf2N]

The shear viscosity of room-temperature ionic liquid (IL) 1-n-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [bmim][Tf2N] is calculated over a temperature range 298-353 K, using the reverse nonequilibrium molecular dynamics simulation technique. The results of this work show that while the use of equilibrium molecular dynamics simulation techniques might be inefficient for viscosity calculations of ILs, the reverse nonequilibrium molecular dynamics technique is an efficient tool for this purpose. Our findings indicate that the shear rate for crossover from the Newtonian plateau to the shear thinning regime, corresponds to the relaxation time for the slowest microscopic scale motions, i.e., exchange of counterions in an ion's solvation shell (ion-pair relaxation time). The closeness of the time scales and activation energies for zero-shear-rate viscosities to the relaxation times and the corresponding activation energies for ion-pair formation/rupture, connects macroscopic dynamic properties with local atomic-level motions of the IL. The calculated viscosity coefficients and relaxation times for reorientations of the cation and anion as well as their corresponding activation energies are in very good agreement with experimental data.

Rheology of phosphonium ionic liquids: a molecular dynamics and experimental study

Comparison between the theoretical and experimental viscosity of an ionic liquid.

Unique orientations and rotational dynamics of a 1-butyl-3-methyl-imidazolium hexafluorophosphate ionic liquid at the gas–liquid interface: the effects of the hydrogen bond and hydrophobic interactions

Two-stage rotational motions of the interfacial [BMIM]⁺ cations are essentially determined by both hydrophobic and hydrogen-bonding interactions.

Neural network modeling as an efficient approach to predict the density of ionic liquids/ethanol binary systems

Ionic liquids (ILs) especially their mixtures are of high interest within the different scientific societies due to their amazing properties. In this regard, a number of attempts have been made to measure, correlate, estimate and calculate the properties of ILs in the neat or mixture forms. Among the different possible predictive methods, artificial neural networks (ANNs) are widely used because of their unique and amazing capabilities for prediction of different parameters. With respect to this paper, a feed-forward ANN model is proposed to model the densities of different binary mixtures of ILs/ethanol. The proposed network is trained and tested with 1078 binary data points gathered by mining into the different published literatures. The data gathered from previously published literatures are separated into two different subsets namely training and testing. The statistical error analysis has shown that the proposed neural network correlated the binary densities with the overall mean absolute percentage error (MAPE), average relative deviation percentage error (ARD%), minimum relative deviation percent (RDmin%), maximum relative deviation Percent (RDmax%) and correlation coefficient ([Formula: see text] of 1.5%, [Formula: see text]0.1%, [Formula: see text]13.0%, 15.0% and 0.9712, respectively.

A molecular dynamics simulations study on the relations between dynamical heterogeneity, structural relaxation, and self-diffusion in viscous liquids

We perform molecular dynamics simulations for viscous liquids to study the relations between dynamical heterogeneity, structural (α) relaxation, and self-diffusion. For atomistic models of supercooled water, polymer melts, and an ionic liquid, we characterize the space-time characteristics of dynamical heterogeneity by the degree of deviations from Gaussian displacement statistics (α2), the size of clusters comprising highly mobile particles (S(w)), and the length of strings consisting of cooperatively moving particles (L(w)). Comparison of our findings with previous simulation results for a large variety of viscous liquids, ranging from monoatomic liquids to silica melt, reveals a nearly universal decoupling between the time scales of maximum non-Gaussian parameter (τ(α2)) and the time constant of the α relaxation (τ(α)) upon cooling, explicitly, τ(α2) ∝τ(α)(3/4). Such uniform relation was not observed between the peak times of S(w) or L(w) and τ(α). On the other hand, the temperature-dependent time scale of maximum string length (τ(L)) follows the inverse of the self-diffusion coefficient (D) for various systems at sufficiently low temperatures, i.e., τ(L) ∝ D(-1). These observations are discussed in view of a breakdown of the Stokes-Einstein relation for the studied systems. It is found that the degree of deviation from this relation is correlated with the stretching of the α relaxation.

Computational studies of ionic liquids: size does matter and time too

Taking the molecular ionic liquid 1-ethyl-3-methylimidazolium triflate as a reference system, the size and time dependence of molecular dynamics simulation studies is analyzed in a systematic way. Based on an all atom force field, trajectories of 70 ns length, covering samples of 8-2000 ion pairs, were generated and analyzed in terms of structure as well as single particle and collective dynamics. Although 50 ion pairs seemed sufficient for structure, at least 500 ion pairs were needed for the correct handling of dynamics. For larger systems a linear regime is found, i.e., the respective dynamical properties are a linear function of the inverse box length. In case of translational diffusion coefficients, this linear relation can be rationalised in hydrodynamic terms. The respective formula is essentially determined by viscosity and the inverse box length. Concerning the time dependence, consistent dynamical properties required a time period of 20-30 ns. Nevertheless, size dependence dominates time dependence and has to be primarily addressed.

Molecular dynamics simulations of nanoparticle self-assembly at ionic liquid-water and ionic liquid-oil interfaces

We have studied the self-assembly of hydrophobic nanoparticles at ionic liquid (IL)-water and IL-oil (hexane) interfaces using molecular dynamics (MD) simulations. For the 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF(6)])/water system, the nanoparticles rapidly approached the IL-water interface and equilibrated more into the IL phase although they were initially in the water phase. In contrast, when the nanoparticles were dispersed in the hexane phase, they slowly approached the IL-hexane interface but remained primarily in the hexane phase. Consequently, the IL-hexane interface was rather undisturbed by the nanoparticles whereas the IL-water interface changed significantly in width and morphology to accommodate the presence of the nanoparticles. The equilibrium positions of the nanoparticles were also supported and explained by potential of mean force (PMF) calculations. Interesting ordering and charge distributions were observed at the IL-liquid interfaces. At the IL-hexane interface, the [BMIM] cations preferentially oriented themselves so that they were immersed more in the hexane phase and packed efficiently to reduce steric hindrance. The ordering likely contributed to a heightened IL density and a slightly positive charge at the IL-hexane interface. In contrast, the cations at the IL-water interface were oriented isotropically unless in the presence of nanoparticles, where the cations aligned across the nanoparticle surfaces.

Correlations in the Structure and Dynamics of Ionic Liquids

A very great deal of the experimental work on room-temperature ionic liquids was done after high-performance computing became readily available for quantum-chemical or molecular-dynamic calculations. I explore the use of modern computational methods to guide or aid laboratory work, and the importance of ‘old-fashioned’ theory, from before the age of fast computers. Debye and Hückel published the first really important theoretical work on correlations between charged particles, and the Nernst–Einstein formula is still used to understand electrical conductivities. I assess the usefulness of all these theoretical methods and ideas, and discuss the particular difficulties presented by ionic liquids.