Computational investigation of non-covalent interactions in 1-butyl 3-methylimidazolium/bis(trifluoromethylsulfonyl)imide [bmim][Tf2N] in EMD and NEMD.

Abstract

Non-covalent interactions (NCIs) play a crucial role in the behavior and properties of ionic liquids (ILs). These interactions are particularly important for non-equilibrium properties such as the change in viscosity due to shearing forces (shear viscosity). Therefore, a detailed understanding of these interactions can improve our understanding of these important classes of liquids. Here, we have employed quantum mechanical energy decomposition analysis (EDA) and NCI analysis to investigate a series of representative 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([bmim][Tf2N]) ion pairs extracted from classical equilibrium and non-equilibrium molecular dynamics simulations. EDA based on symmetry-adapted perturbation theory (SAPT) for the complete monomers, as well as fragment SAPT (FSAPT), for the functional fragments has been carried out. In general, the electrostatic component comprises ≈80% of the intermolecular interaction, and significant contributions from other components (induction and dispersion) are also observed, especially for interactions involving bifurcated hydrogen bonds. The FSAPT analysis suggests that caution is warranted when employing simplified assumptions for non-bonded interactions, e.g., focusing only on hydrogen bonds between functional fragments, since this view may not provide a complete picture of the complicated interactions between the ions. In non-equilibrium molecular dynamics, the total interaction energies of some fragments have a significant qualitative change as the shear rate increases. Our results indicate that the inter-fragment interactions play a fundamental role in the viscous behavior of ILs, suggesting that the exclusive use of geometric criteria to analyze inter-molecular interactions in these systems is not sufficient to investigate shear-thinning effects.