A molecular dynamics investigation of the structural and dynamic properties of the ionic liquid 1-n-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide

Abstract

Molecular dynamics simulations have been performed to investigate the structure and dynamics of the ionic liquid, 1-n-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([C(4)mim][Tf(2)N]) in the temperature range of 283-460 K. Extensive analysis was carried out to characterize a number of structural and dynamic features. Transport properties were computed using a variety of equilibrium methods that employed the Green-Kubo and Einstein formulations. Nonequilibrium techniques were also used. In general, different methods mostly yielded consistent results, although some differences were observed. Computed self-diffusivities and ionic conductivities tended to be slightly lower than experimental values, while computed viscosities were significantly higher than experiment. Computed thermal conductivities agreed reasonably well with experimental data. Despite these discrepancies, the simulations capture the experimental temperature-dependent trends for all these transport properties. Single ion dynamics were studied by examining diffusional anisotropy, the self-part of the van Hove function, non-Gaussian parameters, and incoherent intermediate scattering functions. It is found that cations diffuse faster than anions and are more dynamically heterogeneous. A clear anisotropy is revealed in cation displacement, with the motion normal to the imidazolium ring plane being the most hindered and the motion along the alkyl chain in the plane of the ring being the most facile. Cations structurally relax faster than anions but they rotationally relax slower than anions. There is a pronounced temperature dependence to the rotational anisotropy of the cations, but only a weak temperature dependence for the anions. The ionic conductivity deviates from the Nernst-Einstein relation due to the correlated motion of cations and anions. The results suggest that the dynamical behavior of this and related ionic liquids is extremely complex and consists of many different modes with widely varying timescales, making the prediction of dynamical trends extremely difficult.