Kinetics of Hydrogen Bonding between Ions with Opposite and Like Charges in Hydroxyl-Functionalized Ionic Liquids

The Coiling Effect in Ether Ionic Liquids: Exploiting Acetate as a Probe for Transport Properties and Microenvironment Analysis

The inherently high viscosity of ionic liquids (ILs) can limit their potential applications. One approach to address this drawback is to modify the cation side chain with ether groups. Herein, we assessed the structure-property relationship by focusing on acetate (OAc), a strongly coordinating anion, with 1,3-dialkylimidazolium cations with different side chains, including alkyl, ether, and hydroxyl functionalized, as well as their combinations. We evaluated their viscosity, thermal stabilities, and microstructure using Raman and infrared (IR) spectroscopies, allied to density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. The viscosity data showed that the ether insertion significantly enhances the fluidity of the ILs, consistent with the coiling effect of the cation chain. Through a combined experimental and theoretical approach, we analyzed how the OAc anion interacts with ether ILs, revealing a characteristic bidentate coordination, particularly in hydroxyl functionalized ILs due to specific hydrogen bonding with the OH group. IR spectroscopy showed subtle shifts in the acidic hydrogens of imidazolium ring C(2)-H and C(4,5)-H, suggesting weaker interactions between OAc and the imidazolium ring in ether-functionalized ILs. Additionally, spatial distribution functions (SDF) and dihedral angle distribution obtained via AIMD confirmed the intramolecular hydrogen bonding due to the coiling effect of the ether side chain.

Ion Mobility in Hydroxy-Functionalized Ionic Liquids Depends on Cationic Clustering: Tracking the Alkyl Chain Length Behavior with Deuteron NMR Relaxation

Observing and quantifying the like-charge attraction in liquids and solutions is still challenging. However, we showed that elusive cation-cation hydrogen bonding may govern the structure and interaction in hydroxyl-functionalized ionic liquids. Therefore, cationic cluster formation depends on the shape, charge distribution, and functionality of the ions. We demonstrated by means of solid-state 2H NMR spectroscopy that cationic clusters change the structure and dynamics of ionic liquids. With increasing alkyl chain length, we observed two deuteron quadrupole coupling constants for the OD groups, differing by about 30 kHz. The lower value was assigned to the cation-cation interaction, indicating that the average (c-c) hydrogen bonds are stronger than the (c-a) hydrogen bonds between the cation and the anion despite the repulsive and attractive Coulomb interaction in the first and latter cases. Ion mobility could be studied by 2H NMR spectroscopy, although the deuterons in the hydrogen-bonded clusters underwent fast exchange. Our results also showed that simple relaxation models are not applicable anymore and that anisotropic motion must be considered.

How Like-Charge Attraction Influences the Mobility of Cations in Hydroxyl-Functionalized Ionic Liquids

Attractive interactions between ions of like charge remain an elusive concept. Observing and quantifying this type of interaction in liquids and solutions is still a major challenge. Recently, we have shown that cation-cation interactions are present in hydroxyl-functionalized ionic liquids and that they can be controlled by the shape, charge distribution and functionality of the ions. In the present study, we demonstrate that cationic cluster formation does not only change the local structures of the ionic liquids but also influences the dynamics of the cations in a characteristic way. We show that solid-state ²H NMR spectroscopy is well suited for the study of molecular motion, even if the hydrogen bonded species of interest are indistinguishable due to fast deuteron exchange. We also provide valuable information about the applicability of well-accepted relaxation models.

Quantification and Distribution of Three Types of Hydrogen Bonds in Mixtures of an Ionic Liquid with the Hydrogen-Bond-Accepting Molecular Solvent DMSO Explored by Neutron Diffraction and Molecular Dynamics Simulations

The concept of hydrogen bonding is celebrating its 100th birthday. Hydrogen bonds (H-bonds) play a key role in the structure and function of biological molecules, the strength of materials, and molecular binding. Herein, we study H-bonding in mixtures of a hydroxyl-functionalized ionic liquid with the neutral, H-bond-accepting molecular liquid dimethylsulfoxide (DMSO) using neutron diffraction experiments and molecular dynamics simulations. We report the geometry, strength, and distribution of three different types of H-bond OH···O, formed between the hydroxyl group of the cation and either the oxygen atom of another cation, the counteranion, or the neutral molecule. Such a variety of different strengths and distributions of H-bonds in one single mixture could hold the promise of providing solvents with potential applications in H-bond-related chemistry, for example, to alter the natural selectivity patterns of catalytic reactions or the conformation of catalysts.

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.

An exact a posteriori correction for hydrogen bond population correlation functions and other reversible geminate recombinations obtained from simulations with periodic boundary conditions. Liquid water as a test case

The kinetics of breaking and re-formation of hydrogen bonds (HBs) in liquid water is a prototype of reversible geminate recombination. HB population correlation functions (HBPCFs) are a means to study the HB kinetics. The long-time limiting behavior of HBPCFs is controlled by translatoric diffusion and shows a t^-3/2 time-dependence, which can be described by analytical expressions based on the HB acceptor density and the donor-acceptor inter-diffusion coefficient. If the trajectories are not properly "unwrapped," the presence of periodic boundary conditions (PBCs) can perturb this long-time limiting behavior. Keeping the trajectories "wrapped," however, allows for a more efficient calculation of HBPCFs. We discuss the consequences of PBCs in combination with "wrapped" trajectories following from the approximations according to Luzar-Chandler and according to Starr, each deviating in a different fashion from the true long-time limiting behavior, but enveloping the unperturbed function. A simple expression is given for estimating the maximum time up to which the computed HBPCFs reliably describe the long-time limiting behavior. In addition, an exact a posteriori correction for systems with PBCs for "wrapped" trajectories is derived, which can be easily computed and which is able to fully recover the true t^-3/2 long-time behavior. For comparison, HBPCFs are computed from MD simulations of TIP4P/2005 model water for varying system sizes and temperatures of 273 and 298 K using this newly introduced correction. Implications for the computations of HB lifetimes and the effect of the system-size are discussed.

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

Hydroxyl-functionalized ionic liquids (ILs) represent a new interesting class of ILs where hydrogen bonds (HBs) play an important role: here, "typical" HBs between cations and anions (ca) are competing with "atypical" HBs connecting pairs of cations (cc). We study the equilibrium and kinetics of (cc) and (ca) HBs in 1-(n-hydroxyalkyl)-pyridinium bis(trifluoromethlysulfonyl)imide [HOC_nPy][NTf₂] ILs by means of molecular dynamics simulations. (cc) HBs are found to be between 0.96 and 3.76 kJ mol^-1 stronger than their (ca) counterparts, depending on the alkyl chain length. HB lifetimes and kinetics are analyzed by means of HB population and reactive flux correlation functions. Essentially, four different HB lifetimes have to be considered, spanning about 3 orders of magnitude, each valid in its own right and each associated with different aspects of HB breaking and HB reformation. The long-time limiting behavior of the HB population correlation function is controlled by diffusion of the ions and can be quantitatively described by analytical expressions. The short-time HB behavior is tied to the localized dynamics of the hydroxyl group exploring its local solvation environment. A minimalist kinetic two-domain model is introduced to realistically describe the time evolution of the HB population correlation function for both (ca) and (cc) HBs over 5 orders of magnitude. By employing the reactive flux method, we determine the kinetics of HB breaking, unaffected by diffusion processes. We determine both, the ultrafast upper boundary and the average rate of HB breaking, allowing recrossing-events during the transient relaxation time period. For sufficiently long alkyl chains, all those computed HB lifetimes indicate a higher kinetic stability of (cc) HBs over (ca) HBs; for short chains, it is vice-versa.