Insights into structure-property relationships in ionic liquids using cyclic perfluoroalkylsulfonylimides

Room temperature ionic liquids of cyclic sulfonimide anions ncPFSI (ring size: n = 4-6) with the cations [EMIm]+ (1-ethyl-3-methylimidazolium), [BMIm]+ (1-butyl-3-methylimidazolium) and [BMPL]+ (BMPL = 1-butyl-1-methylpyrrolidinium) have been synthesized. Their solid-state structures have been elucidated by single-crystal X-ray diffraction and their physicochemical properties (thermal behaviour and stability, dynamic viscosity and specific conductivity) have been assessed. In addition, the ion diffusion was studied by pulsed field gradient stimulated echo (PFGSTE) NMR spectroscopy. The decisive influence of the ring size of the cyclic sulfonimide anions on the physicochemical properties of the ILs has been revealed. All ILs show different properties compared to those of the non-cyclic TFSI anion. While these differences are especially distinct for ILs with the very rigid 6cPFSI anion, the 5-membered ring anion 5cPFSI was found to result in ILs with relatively similar properties. The difference between the properties of the TFSI anion and the cyclic sulfonimide anions has been rationalized by the rigidity (conformational lock) of the cyclic sulfonimide anions. The comparison of selected IL properties was augmented by MD simulations. These highlight the importance of π+-π+ interactions between pairs of [EMIm]+ cations in the liquid phase. The π+-π+ interactions are evident for the solid state from the molecular structures of the [EMIm]+-ILs with the three cyclic imide anions determined by single-crystal X-ray diffraction.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

What Is Special about Aromatic-Aromatic Interactions? Significant Attraction at Large Horizontal Displacement

High-level ab initio calculations show that the most stable stacking for benzene-cyclohexane is 17% stronger than that for benzene-benzene. However, as these systems are displaced horizontally the benzene-benzene attraction retains its strength. At a displacement of 5.0 Å, the benzene-benzene attraction is still ∼70% of its maximum strength, while benzene-cyclohexane attraction has fallen to ∼40% of its maximum strength. Alternatively, the radius of attraction (>2.0 kcal/mol) for benzene-benzene is 250% larger than that for benzene-cyclohexane. Thus, at relatively large distances aromatic rings can recognize each other, a phenomenon that helps explain their importance in protein folding and supramolecular structures.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Ion pairing controls rheological properties of "processionary" polyelectrolyte hydrogels

We demonstrated recently that polyelectrolytes with cationic moieties along the chain and a single anionic head are able to form physical hydrogels due to the reversible nature of the head-to-body ionic bond. Here we generate a variety of such polyelectrolytes with various cationic moieties and counterion combinations starting from a common polymeric platform. We show that the rheological properties (shear modulus, critical strain) of the final hydrogels can be modulated over three orders of magnitude depending on the cation/anion pair. Our data fit remarkably well within a scaling model involving a supramolecular head-to-tail single file between cross-links, akin to the behaviour of pine-processionary caterpillar. This model allows the quantitative measure of the amount of counterion condensation from standard rheology procedure.

Electronic Structure, NMR, Spin-Spin Coupling, and Noncovalent Interactions in Aromatic Amino Acid Based Ionic Liquids

Noncovalent interactions accompanying phenylalanine (Phe), tryptophan (Trp), and tyrosine (Tyr) amino acids based ionic liquids (AAILs) composed of 1-methyl-3-butyl-imidazole and its methyl-substituted derivative as cations have been analyzed employing the dispersion corrected density functional theory. It has been shown that cation-anion binding in these bioionic ILs is primarily facilitated through hydrogen bonding in addition to lp---π and CH---π interactions those arising from aromatic moieties which can be probed through (1)H and (13)C NMR spectra calculated from the gauge independent atomic orbital method. Characteristic NMR spin-spin coupling constants across hydrogen bonds of ion pair structures viz., Fermi contact, spin-orbit and spin-dipole terms show strong dependence on mutual orientation of cation with the amino acid anion. The spin-spin coupling mechanism transmits spin polarization via electric field effect originating from lp---π interactions whereas the electron delocalization from lone pair on the carbonyl oxygen to antibonding C-H orbital is facilitated by hydrogen bonding. It has been demonstrated that indirect spin-spin coupling constants across the hydrogen bonds correlate linearly with hydrogen bond distances. The binding energies and dissected nucleus independent chemical shifts (NICS) document mutual reduction of aromaticity of hydrogen bonded ion pairs consequent to localization of π-character. Moreover the nature and type of such noncovalent interactions governing the in-plane and out-of-plane NICS components provide a measure of diatropic and paratropic currents for the aromatic rings of varying size in AAILs. Besides the direction of frequency shifts of characteristic C═O and NH stretching vibrations in the calculated vibrational spectra has been rationalized.

Hydrogen bonding and π–π interactions in imidazolium-chloride ionic liquid clusters

The importance of 1° and 2° hydrogen-bonding and anion–π⁺ interactions for ionic liquid structuring.

Competitive pi interactions and hydrogen bonding within imidazolium ionic liquids

In this paper we have explored the structural and energetic landscape of potential π(+)-π(+) stacked motifs, hydrogen-bonding arrangements and anion-π(+) interactions for gas-phase ion pair (IP) conformers and IP-dimers of 1,3-dimethylimidazolium chloride, [C1C1im]Cl. We classify cation-cation ring stacking as an electron deficient π(+)-π(+) interaction, and a competitive anion on-top IP motif as an anion-donor π(+)-acceptor interaction. 21 stable IP-dimers have been obtained within an energy range of 0-126 kJ mol(-1). The structures have been found to exhibit a complex interplay of structural features. We have found that low energy IP-dimers are not necessarily formed from the lowest energy IP conformers. The sampled range of IP-dimers exhibits new structural forms that cannot be recovered by examining the ion-pairs alone, moreover the IP-dimers are recovering additional key features of the local liquid structure. Including dispersion is shown to impact both the relative energy ordering and the geometry of the IPs and IP-dimers, however the impact is found to be subtle and dependent on the underlying functional.

Nature of noncovalent interactions in catenane supramolecular complexes: calibrating the MM3 force field with ab initio, DFT, and SAPT methods

The design and assembly of mechanically interlocked molecules, such as catenanes and rotaxanes, are dictated by various types of noncovalent interactions. In particular, [C-H···O] hydrogen-bonding and π-π stacking interactions in these supramolecular complexes have been identified as important noncovalent interactions. With this in mind, we examined the [3]catenane 2·4PF6 using molecular mechanics (MM3), ab initio methods (HF, MP2), several versions of density functional theory (DFT) (B3LYP, M0X), and the dispersion-corrected method DFT-D3. Symmetry adapted perturbation theory (DFT-SAPT) provides the highest level of theory considered, and we use the DFT-SAPT results both to calibrate the other electronic structure methods, and the empirical potential MM3 force field that is often used to describe larger catenane and rotaxane structures where [C-H···O] hydrogen-bonding and π-π stacking interactions play a role. Our results indicate that the MM3 calculated complexation energies agree qualitatively with the energetic ordering from DFT-SAPT calculations with an aug-cc-pVTZ basis, both for structures dominated by [C-H···O] hydrogen-bonding and π-π stacking interactions. When the DFT-SAPT energies are decomposed into components, we find that electrostatic interactions dominate the [C-H···O] hydrogen-bonding interactions, while dispersion makes a significant contribution to π-π stacking. Another important conclusion is that DFT-D3 based on M06 or M06-2X provides interaction energies that are in near-quantitative agreement with DFT-SAPT. DFT results without the D3 correction have important differences compared to DFT-SAPT, while HF and even MP2 results are in poor agreement with DFT-SAPT.