Hydrogen Bonds between Ions of Opposite and Like Charge in Hydroxyl-Functionalized Ionic Liquids: an Exhaustive Examination of the Interplay between Global and Local Motions and Intermolecular Hydrogen Bond Lifetimes and Kinetics

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach

In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.

Ionic liquids revolutionizing biomedicine: recent advances and emerging opportunities

Ionic liquids (ILs), due to their inherent structural tunability, outstanding miscibility behavior, and excellent electrochemical properties, have attracted significant research attention in the biomedical field. As the application of ILs in biomedicine is a rapidly emerging field, there is still a need for systematic analyses and summaries to further advance their development. This review presents a comprehensive survey on the utilization of ILs in the biomedical field. It specifically emphasizes the diverse structures and properties of ILs with their relevance in various biomedical applications. Subsequently, we summarize the mechanisms of ILs as potential drug candidates, exploring their effects on various organisms ranging from cell membranes to organelles, proteins, and nucleic acids. Furthermore, the application of ILs as extractants and catalysts in pharmaceutical engineering is introduced. In addition, we thoroughly review and analyze the applications of ILs in disease diagnosis and delivery systems. By offering an extensive analysis of recent research, our objective is to inspire new ideas and pathways for the design of innovative biomedical technologies based on ILs.

Quantitative Relation between Ionic Diffusivity and Ionic Association in Ionic Liquid-Water Mixtures

Molecular dynamic simulations of aqueous mixtures of imidazolium ionic liquids (ILs) were performed to elucidate the dependence of the ionic diffusivity on the microscopic structures changed by water. Two distinct regimes of the average ionic diffusivity (D_ave) were identified with the increased water concentrations: the jam regime with slowly increased D_ave and the exponential regime with rapidly increased D_ave, which are found to be directly correlated to the ionic association. Further analysis leads to two general relationships independent of IL species between D_ave and the degree of ionic association: (i) a consistent linear relationship between D_ave and the inverse of ion-pair lifetimes (1/τ_IP) in the two regimes and (ii) an exponential relationship between normalized diffusivities (D̃_ave) and short-ranged interactions between cations and anions (Ẽ_ions), with different interdependent strengths in the two regimes. These findings revealed and quantified the direct correlation between dynamic properties and ionic association in IL-water mixtures.

Perspectives in the Computational Modeling of New Generation, Biocompatible Ionic Liquids

In this Perspective, I review the current state of computational simulations on ionic liquids with an emphasis on the recent biocompatible variants. These materials are used here as an example of relatively complex systems that highlights the limits of some of the approaches commonly used to study their structure and dynamics. The source of these limits consists of the coexistence of nontrivial electrostatic, many-body quantum effects, strong hydrogen bonds, and chemical processes affecting the mutual protonation state of the constituent molecular ions. I also provide examples on how it is possible to overcome these problems using suitable simulation paradigms and recently improved techniques that, I expect, will be gradually introduced in the state-of-the-art of computational simulations of ionic liquids.

Thermodynamically Stable Cationic Dimers in Carboxyl-Functionalized Ionic Liquids: The Paradoxical Case of "Anti-Electrostatic" Hydrogen Bonding

We show that carboxyl-functionalized ionic liquids (ILs) form doubly hydrogen-bonded cationic dimers (c⁺=c⁺) despite the repulsive forces between ions of like charge and competing hydrogen bonds between cation and anion (c⁺–a⁻). This structural motif as known for formic acid, the archetype of double hydrogen bridges, is present in the solid state of the IL 1−(carboxymethyl)pyridinium bis(trifluoromethylsulfonyl)imide [HOOC−CH₂−py][NTf₂]. By means of quantum chemical calculations, we explored different hydrogen-bonded isomers of neutral (HOOC–(CH₂)_n–py⁺)₂(NTf₂⁻)₂, single-charged (HOOC–(CH₂)_n–py⁺)₂(NTf₂⁻), and double-charged (HOOC– (CH₂)_n−py⁺)₂ complexes for demonstrating the paradoxical case of “anti-electrostatic” hydrogen bonding (AEHB) between ions of like charge. For the pure doubly hydrogen-bonded cationic dimers (HOOC– (CH₂)_n−py⁺)₂, we report robust kinetic stability for n = 1–4. At n = 5, hydrogen bonding and dispersion fully compensate for the repulsive Coulomb forces between the cations, allowing for the quantification of the two equivalent hydrogen bonds and dispersion interaction in the order of 58.5 and 11 kJmol⁻¹, respectively. For n = 6–8, we calculated negative free energies for temperatures below 47, 80, and 114 K, respectively. Quantum cluster equilibrium (QCE) theory predicts the equilibria between cationic monomers and dimers by considering the intermolecular interaction between the species, leading to thermodynamic stability at even higher temperatures. We rationalize the H-bond characteristics of the cationic dimers by the natural bond orbital (NBO) approach, emphasizing the strong correlation between NBO-based and spectroscopic descriptors, such as NMR chemical shifts and vibrational frequencies.