Magnetic Ionic Liquids: Current Achievements and Future Perspectives with a Focus on Computational Approaches

Magnetic ionic liquids (MILs) stand out as a remarkable subclass of ionic liquids (ILs), combining the desirable features of traditional ILs with the unique ability to respond to external magnetic fields. The incorporation of paramagnetic species into their structures endows them with additional attractive features, including thermochromic behavior and luminescence. These exceptional properties position MILs as highly promising materials for diverse applications, such as gas capture, DNA extractions, and sensing technologies. The present Review synthesizes key experimental findings, offering insights into the structural, thermal, magnetic, and optical properties across various MIL families. Special emphasis is placed on unraveling the influence of different paramagnetic species on MILs' behavior and functionality. Additionally, the Review highlights recent advancements in computational approaches applied to MIL research. By leveraging molecular dynamics (MD) simulations and density functional theory (DFT) calculations, these computational techniques have provided invaluable insights into the underlying mechanisms governing MILs' behavior, facilitating accurate property predictions. In conclusion, this Review provides a comprehensive overview of the current state of research on MILs, showcasing their special properties and potential applications while highlighting the indispensable role of computational methods in unraveling the complexities of these intriguing materials. The Review concludes with a forward-looking perspective on the future directions of research in the field of magnetic ionic liquids.

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.

Synthesis and Physical Properties of Tunable Aryl Alkyl Ionic Liquids (TAAILs) Comprising Imidazolium Cations Blocked with Methyl-, Propyl- and Phenyl-Groups at the C2 Position

Imidazolium-based ionic liquids are very popular for different applications because of their low viscosity and melting point. However, the hydrogen atom at the C2 position of the imidazolium cation can easily be deprotonated by a base, resulting in a reactive carbene. If an inert ionic liquid is needed, it is necessary to introduce an unreactive alkyl or aryl group at the C2 position to prevent deprotonation. Tunable aryl alkyl ionic liquids (TAAILs) were first introduced by our group in 2009 and are characterized by a phenyl group at the N1 position, which offers the possibility to fine-tune the physicochemical properties by using different electron-donating or -withdrawing substituents. In this work, we present a new series of TAAILs where the C2 position is blocked by a methyl, propyl or phenyl group. For each of the blocking groups, the phenyl and three different phenyl derivatives (2-Me, 4-OMe, 2,4-F₂ ) are compared with respect to melting point, viscosity, conductivity and electrochemical window. In addition, the differences between blocked and unblocked TAAILs with regard to their electrochemical reduction potentials are investigated by quantum chemical methods.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

Effect of the cation structure on the properties of homobaric imidazolium ionic liquids

In this work we investigate the structure-property relationships in a series of alkylimidazolium ionic liquids with almost identical molecular weight. Using a combination of theoretical calculations and experimental measurements, we have shown that re-arranging the alkyl side chain or adding functional groups results in quite distinct features in the resultant ILs. The synthesised ILs, although structurally very similar, cover a wide spectrum of properties ranging from highly fluid, glass forming liquids to high melting point crystalline salts. Theoretical ab initio calculations provide insight on minimum energy orientations for the cations, which then are compared to experimental X-ray crystallography measurements to extract information on hydrogen bonding and to verify our understanding of the studied structures. Molecular dynamics simulations of the simplest (core) ionic liquids are used in order to help us interpret our experimental results and understand better why methylation of C² position of the imidazolium ring results in ILs with such different properties compared to their non-methylated analogues.

Unravelling the Interactions of Magnetic Ionic Liquids by Energy Decomposition Schemes: Towards a Transferable Polarizable Force Field

This work aims at unravelling the interactions in magnetic ionic liquids (MILs) by applying Symmetry-Adapted Perturbation Theory (SAPT) calculations, as well as based on those to set-up a polarisable force field model for these liquids. The targeted MILs comprise two different cations, namely: 1-butyl-3-methylimidazolium ([Bmim]⁺) and 1-ethyl-3-methylimidazolium ([Emim]⁺), along with several metal halides anions such as [FeCl₄]⁻, [FeBr₄]⁻, [ZnCl₃]⁻ and [SnCl₄]²⁻ To begin with, DFT geometry optimisations of such MILs were performed, which in turn revealed that the metallic anions prefer to stay close to the region of the carbon atom between the nitrogen atoms in the imidazolium fragment. Then, a SAPT study was carried out to find the optimal separation of the monomers and the different contributions for their interaction energy. It was found that the main contribution to the interaction energy is the electrostatic interaction component, followed by the dispersion one in most of the cases. The SAPT results were compared with those obtained by employing the local energy decomposition scheme based on the DLPNO-CCSD(T) method, the latter showing slightly lower values for the interaction energy as well as an increase of the distance between the minima centres of mass. Finally, the calculated SAPT interaction energies were found to correlate well with the melting points experimentally measured for these MILs.

Targeted modifications in ionic liquids - from understanding to design

Ionic liquids are extremely versatile and continue to find new applications in academia as well as industry. This versatility is rooted in the manifold of possible ion types, ion combinations, and ion variations. However, to fully exploit this versatility, it is imperative to understand how the properties of ionic liquids arise from their constituents. In this work, we discuss targeted modifications as a powerful tool to provide understanding and to enable design. A 'targeted modification' is a deliberate change in the structure of an ionic liquid. This includes chemical changes in an experiment as well as changes to the parameterisation in a computer simulation. In any case, such a change must be purposeful to isolate what is of interest, studying, as far as is possible, only one concept at a time. The concepts can then be used as design elements. However, it is often found that several design elements interact with each other - sometimes synergistically, and other times antagonistically. Targeted modifications are a systematic way of navigating these overlaps. We hope this paper shows that understanding ionic liquids requires experimentalists and theoreticians to join forces and provides a tool to tackle the difficult transition from understanding to design.

Developing design tools for introducing and tuning structural order in ionic liquids

Ionic liquids – ionic crystals – ionic liquid crystals? Structural order in imidazolium-based ILs, a series of asymmetrical 1-dodecyl-2-methyl-3-alkylimidazolium bromides, [C₁₂C₁C_nim][Br] with n = 0–12.

Gas-Phase Binding Energies and Dissociation Dynamics of 1-Alkyl-3-Methylimidazolium Tetrafluoroborate Ionic Liquid Clusters

Ionic liquids (ILs) have become increasingly popular due to their useful and unique properties, yet there are still many unanswered questions regarding their fundamental interactions. In particular, details regarding the nature and strength of the intrinsic cation-anion interactions and how they influence the macroscopic properties of ILs are still largely unknown. Elucidating the molecular-level details of these interactions is essential to the development of better models for describing ILs and enabling the purposeful design of ILs with properties tailored for specific applications. Current uses of ILs are widespread and diverse and include applications for energy storage, electrochemistry, designer/green solvents, separations, and space propulsion. To advance the understanding of the energetics, conformations, and dynamics of gas-phase IL clustering relevant to space propulsion, threshold collision-induced dissociation approaches are used to measure the bond dissociation energies (BDEs) of the 2:1 clusters of 1-alkyl-3-methylimidazolium cations and tetrafluoroborate, [2C_nmim:BF₄]⁺. The cation, [C_nmim]⁺, is varied across the series, 1-ethyl-3-methylimidazolium [C₂mim]⁺, 1-butyl-3-methylimidazolium [C₄mim]⁺, 1-hexyl-3-methylimidazolium [C₆mim]⁺, and 1-octyl-3-methylimidazolium [C₈mim]⁺, to examine the structural and energetic effects of the size of the 1-alkyl substituent on binding. Complementary electronic structure calculations are performed to determine the structures and energetics of the [C_nmim]⁺ and [BF₄]^- ions and their binding preferences in the (C_nmim:BF₄) ion pairs and [2C_nmim:BF₄]⁺ clusters. Several levels of theory, B3LYP, B3LYP-GD3BJ, and M06-2X, using the 6-311+G(d,p) basis set for geometry optimizations and frequency analyses and the 6-311+G(2d,2p) basis set for energetics, are benchmarked to examine their abilities to properly describe the nature of the binding interactions and to reproduce the measured BDEs. The modest structural variation among these [C_nmim]⁺ cations produces only minor structural changes and variation in the measured BDEs of the [2C_nmim:BF₄]⁺ clusters. Present findings indicate that the dominant cation-anion interactions involve the 3-methylimidazolium moieties and that these clusters are sufficiently small that differences in packing effects associated with the variable length of the 1-alkyl substituents are not yet significant.

A Systematic Study of DFT Performance for Geometry Optimizations of Ionic Liquid Clusters

Clusters of two ion pairs of imidazolium-based ionic liquids were optimized with 43 different levels of theory, including DFT functionals and MP2-based methods combined with varying Dunning's basis sets, and added dispersion corrections. Better preforming DFT functionals were then applied to clusters consisting of four ion pairs. Excellent performance of some DFT functionals for the two ion pair clusters did not always match that of the four ion-paired clusters despite interionic distances remaining constant between the optimized two and four ion-paired clusters of the same ionic liquid. Combinations of DFT functional and basis set such as ωB97X-D/cc-pVDZ, M06-2X/aug-cc-pVDZ, B3LYP-D3/cc-pVTZ, and TPSS-D3/cc-pVTZ gave excellent results for geometry optimization of two ion-paired clusters of imidazolium ionic liquids but gave larger deviations when applied to the four ion-paired clusters of varying ionic liquids. Empirical dispersion corrections were seen to be crucial in correctly capturing correlation effects in the studied ionic liquid clusters, becoming more important in larger clusters. Dunning's double-ζ basis set, cc-pVDZ, is associated with the smallest root mean squared deviations for geometries; however, it also produces the largest deviations in total electronic energies. ωB97X-D and M06-2X produced the best performance with the augmented version of this basis set. The triple-ζ basis set, cc-pVTZ, leads to the best performance of most of the DFT functionals (especially the dispersion-corrected ones) used, whereas its augmented version, aug-cc-pVTZ, was not seen to improve results. The combinations of functional and basis set that gave the best geometry and energetics in both two and four ion-paired clusters were PBE-D3/cc-pVTZ, ωB97X-D/aug-cc-pVDZ, and BLYP-D3/cc-pVTZ. All three combinations are recommended for geometry optimizations of larger clusters of ionic liquids. PBE-D3/cc-pVTZ performed the best with an average deviation of 2.3 kJ mol^-1 and a standard deviation of 3.4 kJ mol^-1 for total electronic energy when applied to four ion-paired clusters. Geometries optimized with FMO2-SRS-MP2/cc-pVTZ produced total energy within 2.0 kJ mol^-1 off the benchmark in two ion-paired clusters, with the cc-pVDZ basis set performing unsurprisingly poorly with the same method. The error increased to 4.8 kJ mol^-1 on average in four ion-paired clusters, with the smallest RMSD deviations in geometries when compared to the benchmark ones. This study is the first report that investigated the performance of DFT functionals for two and four ion-paired clusters of a wide range of ionic liquids consisting of commonly used cations such as pyrrolidinium, imidazolium, pyridinium, and ammonium. It also identified the importance of assessing the performance of quantum chemical methods for ionic liquids on a variety of cation-anion combinations.

Low temperature glass/crystal transition in ionic liquids determined by H-bond vs. coulombic strength

Self-assembled ionic liquid crystals are anisotropic ionic conductors, with potential applications in areas as important as solar cells, battery electrolytes and catalysis. However, many of these applications are still limited by the lack of precise control over the variety of phases that can be formed (nematic, smectic, or semi/fully crystalline), determined by a complex pattern of different intermolecular interactions. Here we report the results of a systematic study of crystallization of several imidazolium salts in which the relative contribution of isotropic coulombic and directional H-bond interactions is carefully tuned. Our results demonstrate that the relative strength of directional H-bonds with respect to the isotropic Coulomb interaction determines the formation of a crystalline, semi-crystalline or glassy phase at low temperature. The possibility of pinpointing H-bonding directionality in ionic liquids make them model systems to study the crystallization of an ionic solid under a perturbed Coulomb potential.

Effect of alkyl-group flexibility on the melting point of imidazolium-based ionic liquids

The low melting point of room temperature ionic liquids is usually explained in terms of the presence of bulky, low-symmetry, and flexible ions, with the first two factors related to the lattice energy while an entropic effect is attributed to the latter. By means of molecular dynamics simulations, the melting points of 1-ethyl-3-methyl-imidazolium hexafluorophosphate and 1-decyl-3-methyl-imidazolium hexafluorophosphate were determined, and the effect of the molecular flexibility over the melting point was explicitly computed by restraining the rotation of dihedral angles in both the solid and the liquid phases. The rotational flexibility over the bond between the ring and the alkyl chain affects the relative ordering of the anions around the cations and results in substantial effects over both the enthalpy and the entropy of melting. For the other dihedral angles of the alkyl group, the contributions are predominantly entropic and an alternating behavior was found. The flexibility of some dihedral angles has negligible effects on the melting point, while others can lead to differences in the melting point as large as 20 K. This alternating behavior is rationalized by the different probabilities of conformation defects in the crystal.

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.

Heterogeneity in the microstructure and dynamics of tetraalkylammonium hydroxide ionic liquids: insight from classical molecular dynamics simulations and Voronoi tessellation analysis

Microscopic structural and dynamic heterogeneities were investigated for three ionic liquids (ILs), tetraethylammonium hydroxide, tetrapropylammonium hydroxide, and tetrabutylammonium hydroxide employing classical molecular dynamics (MD) simulations.

Distance Angle Descriptors of the Interionic and Ion-Solvent Interactions in Imidazolium-Based Ionic Liquid Mixtures with Aprotic Solvents: A Molecular Dynamics Simulation Study

The aim of this paper is to quantify the changes of the interionic and ion-solvent interactions in mixtures of imidazolium-based ionic liquids, having tetrafluoroborate (BmimBF₄), hexafluorophosphate (BmimPF₆), trifluoromethylsulfonate (BmimTFO), or bis(trifluoromethanesulfonyl)imide (BmimTFSI), anions, and polar aprotic molecular solvents, such as acetonitrile (AN), γ-butyrolactone (GBL), and propylene carbonate (PC). For this purpose, we calculate, using the nearest-neighbor approach, the average distance between the imidazolium ring H atom in positions 2, 4, and 5 (H^2,4,5) and the nearest high-electronegativity atom of the solvent or anion (X) as distance descriptors, and the mean angle formed by the C^2,4,5-H^2,4,5 bond and the H^2,4,5···X axis around the H^2,4,5 atom as angular descriptors of the cation-anion and cation-solvent interactions around the ring C-H groups. The behavior of these descriptors as a function of the ionic liquid mole fraction is analyzed in detail. The obtained results show that the extent of the change of these descriptors with respect to their values in the neat ionic liquid depends both on the nature of the anion and on the mixture composition. Thus, in the case of the mixtures of the molecular solvents with BmimBF₄ and BmimTFO, a small change of the distance and a drastic increase of the angular descriptor corresponding to the cation-anion interactions are observed with decreasing mole fraction of the ionic liquid, indicating that the anion moves from the above/below position (with respect to the imidazolium ring plane) to a position that is nearly linearly aligned with the C²-H² bond and hinders the possible interaction between the C²-H² group and the solvent molecules. On the other hand, in the case of mixtures of BmimTFSI and BmimPF₆ with the molecular solvents, both the observed increase of the distance descriptor and the slight change of the angular descriptor with decreasing ionic liquid mole fraction are compatible with the direct interactions of the solvent with the C²-H² group. The behavior of these descriptors is correlated with the experimentally observed ¹H chemical shift of the C²-H² group and the red shift of the C²-H² vibrational mode, particularly at low ionic liquid mole fractions. The present results are thus of great help in interpreting these experimental observations.

Computational study of halogen-free Boron based dicationic ionic liquids of [bis-Mim][BMB]2 and [bis-Mim][BScB]2

In this paper, the structures and energetics of 1,3-bis[3-methylimidazolium-yl] pentane ([bis-Mim]²⁺) dication and bis(mandelato) borate [BMB]^- and bis(salicylato) borate [BScB]^-anions in isolated forms, ion pairs (IPs) and dicationic ionic liquids (DILs) were studied by Density Functional Theory (DFT) at the M06-2X/6-31G(d,p) level of theory. According to the IUPAC criteria, the hydrogen bonds between anion and cation were characterized and classified in the optimized geometries of the isolated ions, IPs and DILs. Inspection of the optimized structures revealed that the interionic hydrogen binding has important effect on the ions structures. The interaction energies between a dication and anions have been described in terms of NBO charge distribution, the stabilization energy E⁽²⁾ values, changes of vibrational frequencies and the reduced density gradient (RDG) analysis. Also, the reactivity and interactions between chemical species were interpreted in terms of global electronic properties. Electrostatic potential surfaces (ESP) have been applied for visualizing the charge related properties and characterization of the most energetic sites of isolated ions and ion complexes.

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.

An overview of the performance of the COSMO-RS approach in predicting the activity coefficients of molecular solutes in ionic liquids and derived properties at infinite dilution

An overview of performance of state-of-the-art thermodynamic model COSMO-RS in capturing various effects of structure on interactions in ionic liquid binary systems, expressed in terms of limiting activity coefficients of molecular solutes, is presented.

Nature of the C2-methylation effect on the properties of imidazolium ionic liquids

Methylation at the C2 position of 1,3-disubstituted imidazolium-based ionic liquids (ILs) is one of the structural features that has gained attention due to its drastic impact on thermophysical and transport properties.

Local Structure in Terms of Nearest-Neighbor Approach in 1-Butyl-3-methylimidazolium-Based Ionic Liquids: MD Simulations

Description of the local microscopic structure in ionic liquids (ILs) is a prerequisite to obtain a comprehensive understanding of the influence of the nature of ions on the properties of ILs. The local structure is mainly determined by the spatial arrangement of the nearest neighboring ions. Therefore, the main interaction patterns in ILs, such as cation-anion H-bond-like motifs, cation-cation alkyl tail aggregation, and ring stacking, were considered within the framework of the nearest-neighbor approach with respect to each particular interaction site. We employed classical molecular dynamics (MD) simulations to study in detail the spatial, radial, and orientational relative distribution of ions in a set of imidazolium-based ILs, in which the 1-butyl-3-methylimidazolium (C4mim(+)) cation is coupled with the acetate (OAc(-)), chloride (Cl(-)), tetrafluoroborate (BF4(-)), hexafluorophosphate (PF6(-)), trifluoromethanesulfonate (TfO(-)), or bis(trifluoromethanesulfonyl)amide (TFSA(-)) anion. It was established that several structural properties are strongly anion-specific, while some can be treated as universally applicable to ILs, regardless of the nature of the anion. Namely, strongly basic anions, such as OAc(-) and Cl(-), prefer to be located in the imidazolium ring plane next to the C-H(2/4-5) sites. By contrast, the other four bulky and weakly coordinating anions tend to occupy positions above/below the plane. Similarly, the H-bond-like interactions involving the H(2) site are found to be particularly enhanced in comparison with the ones at H(4-5) in the case of asymmetric and/or more basic anions (C4mimOAc, C4mimCl, C4mimTfO, and C4mimTFSA), in accordance with recent spectroscopic and theoretical findings. Other IL-specific details related to the multiple H-bond-like binding and cation stacking issues are also discussed in this paper. The secondary H-bonding of anions with the alkyl hydrogen atoms of cations as well as the cation-cation alkyl chain aggregation turned out to be poorly sensitive to the nature of the anion.

Protonic Ammonium Nitrate Ionic Liquids and Their Mixtures: Insights into Their Thermophysical Behavior

This study is centered on the thermophysical characterization of different families of alkylammonium nitrate ionic liquids and their binary mixtures, namely the determination at atmospheric pressure of densities, electric conductivities and viscosities in the 288.15 < T/K < 353.15 range. First, measurements focusing on ethylammonium, propylammonium and butylammonium nitrate systems, and their binary mixtures, were determined. These were followed by studies involving binary mixtures composed of ethylammonium nitrate (with three hydrogen bond donor groups) and different homologous ionic liquids with differing numbers of hydrogen bond donor groups: diethylammonium nitrate (two hydrogen bond donors), triethylammonium nitrate (one hydrogen bond donor) and tetraethylammonium nitrate (no hydrogen bond donors). Finally, the behavior of mixtures with different numbers of equivalent carbon atoms in the alkylammonium cations was analyzed. The results show a quasi-ideal behavior for all monoalkylammonium nitrate mixtures. In contrast, the other mixtures show deviations from ideality, namely when the difference in the number of carbon atoms present in the cations increases or the number of hydrogen bond donors present in the cation decreases. Overall, the results clearly show that, besides the length and distribution of alkyl chains present in a cation such as alkylammonium, there are other structural and interaction parameters that influence the thermophysical properties of both pure compounds and their mixtures.

Significance of weak interactions in imidazolium picrate ionic liquids: spectroscopic and theoretical studies for molecular level understanding

The effects of interionic hydrogen bonding and π-π stacking interactions on the physical properties of a new series of picrate anion based ionic liquids (ILs) have been investigated experimentally and theoretically. The existence of aromatic (C2-HO) and aliphatic (C7-HO-N22 and C6-HO-N20) hydrogen bonding and π-π stacking interactions in these ILs has been observed using various spectroscopic techniques. The aromatic and aliphatic C-HO hydrogen bonding interactions are found to have a crucial role in binding the imidazolium cation and picrate anion together. However, the π-π stacking interactions between two successive layers are found to play a decisive role in tight packing in ILs leading to differences in physical properties. The drastic difference in the melting points of the methyl and propyl derivatives (mmimPic and pmimPic respectively) have been found to be primarily due to the difference in the strength and varieties of π-π stacking interactions. While in mmimPic, several different types of π-π stacking interactions between the aromatic rings (such as picrate-picrate, picrate-imidazole and imidazolium-imidazolium cation rings) are observed, only one type of π-π stacking interaction (picrate-picrate rings) is found to exist in the pmimPic IL. NMR spectroscopic studies reveal that the interaction of these ILs with solvent molecules is different and depends on the dielectric constant of the solvent. While an ion solvation model explains the solvation in high dielectric solvents, an ion-pair solvation model is found to be more appropriate for low dielectric constant solvents. The enhanced stability of these investigated picrate ILs compared with that of inorganic picrate salts under high doses of γ radiation clearly indicates the importance of weak interionic interactions in ILs, and also opens up the possibility of the application of picrate ILs as prospective diluents in nuclear separation for advanced fuel cycling process.

Electrical conductivity and glass formation in nitrile-functionalized pyrrolidinium bis(trifluoromethylsulfonyl)imide ionic liquids: chain length and odd-even effects of the alkyl spacer between the pyrrolidinium ring and the nitrile group

The electrical conductivity of a series of pyrrolidinium bis(trifluoromethylsulfonyl)imide ionic liquids, functionalized with a nitrile (cyano) group at the end of an alkyl chain attached to the cation, was studied in the temperature range between 173 K and 393 K. The glass formation of the ionic liquids is influenced by the length of the alkyl spacer separating the nitrile function from the pyrrolidinium ring. The electrical conductivity and the viscosity do not show a monotonic dependence on the alkyl spacer length, but rather an odd-even effect. An explanation for this behavior is given, including the potential energy landscape picture for the glass transition.

Molecular dynamics study of the effect of alkyl chain length on melting points of [CnMIM][PF6] ionic liquids

Based on molecular dynamics simulations, the melting points Tm of a series of 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids [CnMIM][PF6] with n = 2, 4, 10, 12, and 14 were studied using the free energy-based pseudosupercritical path (PSCP) method. The experimental trend that the Tm decreases with increasing alkyl chain length for ILs with short alkyl chains and increases for the ones with long alkyl chains was correctly captured. Further analysis revealed that the different trends are the results of the balance between fusion enthalpy and fusion entropy. For the ILs with short alkyl chains (ethyl and butyl groups), fusion entropy plays the dominant role so that [C4MIM][PF6], which has a larger fusion entropy due to its higher liquid phase entropy has the lower melting temperature. As for the ILs with long alkyl chains, due to the enhanced van der Waals interactions brought about by the long non-polar alkyl chains, enthalpy becomes the deciding factor and the melting points increase when the alkyl chain goes from C10 to C14. While the melting points for [C2MIM][PF6] and [C4MIM][PF6] were quantitatively predicted and the trends for the long chain ILs were captured correctly, the absolute melting points for [C10MIM][PF6], [C12MIM][PF6] and [C14MIM][PF6] were systematically overestimated in the simulations. Three possible reasons for the overestimation were studied but all ruled out. Further simulation or experimental studies are needed to explain the difference.

Molecular features contributing to the lower viscosity of phosphonium ionic liquids compared to their ammonium analogues

Classical molecular dynamics simulations identify the size of the central atom and the angle flexibility as important molecular features contributing to the lower viscosity of phosphonium based ionic liquids compared to their ammonium analogues.