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

Transferability and Accuracy of Ionic Liquid Simulations with Equivariant Machine Learning Interatomic Potentials

Ionic liquids (ILs) are an exciting class of electrolytes finding applications in many areas from energy storage to solvents, where they have been touted as "designer solvents" as they can be mixed to precisely tailor the physiochemical properties. As using machine learning interatomic potentials (MLIPs) to simulate ILs is still relatively unexplored, several questions need to be answered to see if MLIPs can be transformative for ILs. Since ILs are often not pure, but are either mixed together or contain additives, we first demonstrate that a MLIP can be trained to be compositionally transferable; i.e., the MLIP can be applied to mixtures of ions not directly trained on, while only being trained on a few mixtures of the same ions. We also investigated the accuracy of MLIPs for a novel IL, which we experimentally synthesize and characterize. Our MLIP trained on ∼200 DFT frames is in reasonable agreement with our experiments and DFT.

A mechanistic study on conversion of carbon dioxide into formic acid promoted by 1-ethyl-2, 3-dimethyl-imidazolium nitrite

CONTEXT

The conversion of carbon dioxide (CO2) to formic acid (FA) through hydrogenation using 1-ethyl-2,3- dimethyl imidazolium nitrite (EDIN) ionic liquid was studied to understand the catalytic roles within EDIN. CO2 hydrogenation in various solvents has been explored, but achieving high efficiency and selectivity remains challenging due to the thermodynamic stability and kinetic inertness of CO2. This study explored two mechanistic pathways through theoretical calculations, revealing that the nitrite (NO2-) group is the most active site. The oxygen site on nitrite favorably activates H2, while the nitrogen site shows a minor activation barrier of 108.90 kJ/mol. The Gibbs energy variation indicates stable FA formation via EDIN, suggesting effective hydrogen (H2) activation and subsequent CO2 conversion. These insights are crucial for developing improved catalytic sites and processes in ionic liquid catalysts for CO2 hydrogenation.

METHODS

Quantum chemical calculations were conducted using the ORCA software package at the Restricted Hartree-Fock (RHF) and density functional theory (DFT) levels. The RHF method, known for its predictive abilities in simpler systems, provided a baseline description of electronic structures. In contrast, DFT was employed for its effectiveness in complex interactions involving significant electron correlation. A valence triple-zeta polarization (def2-TZVPP) basis set was employed for both RHF and DFT, ensuring accurate and correlated calculations. The B3LYP functional was utilized for its rapid convergence and cost-efficiency in larger molecules. Dispersion corrected functionals (DFT-D) addressed significant dispersion forces in ionic liquids, incorporating Grimme's D2, D3, and D4 corrections. Geometry optimizations, kinetics, and thermodynamic calculations were performed in the gas phase. The Nudged Elastic Band Transition State (NEB-TS) approach, combining Climbing Image-NEB (CINEB) and Eigenvector-Following (EF) methods, was used to find the minimum energy path (MEP) between reactants and products. Thermochemical analyses based on vibrational frequency calculations evaluated properties such as Enthalpy, Entropy, and Gibbs energy using ideal gas statistical mechanics.

Evolution of Surface Tension and Hansen Parameters of a Homologous Series of Imidazolium-Based Ionic Liquids

In this study, we systematically analyze the surface tension and Hansen solubility parameters (HSPs) of imidazolium-based ionic liquids (ILs) with different anions ([NTf2]-, [PF6]-, [I]-, and [Br]-). These anions are combined with the classical 1-alkyl-3-methyl-substituted imidazolium cations ([CnC1Im]+) and a group of oligoether-functionalized imidazolium cations ([(mPEGn)2Im]+) based on methylated polyethylene glycol (mPEGn). In detail, the influences of the length of the alkyl- and the mPEGn-chain, the anion size, and the water content are investigated experimentally. For [CnC1Im]+-based ILs, the surface tension decreases with increasing alkyl chain length in all cases, but the magnitude of this decrease depends on the size of the anion ([NTf2]- < [PF6]- < [Br]- ≤ [I]-). Molecular dynamics (MD) simulations on [CnC1Im]+-based ILs indicate that these differences are caused by the interplay of charged and uncharged domains, in particular in the different anions, which affects the ability of the alkyl chains of the cation to orient toward the liquid-gas interface. An increase in the mPEGn-chain length of the [(mPEGn)2Im][A] ILs does not significantly influence the surface tension. These changes upon variation of the cation/anion combination do not correlate with the evolution of the HSPs for the two sets of ILs. Finally, our data suggest that significant water contents up to water mole fractions of x(H2O) = 0.25 do not significantly affect the surface tension of the studied binary IL-water mixtures.

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.

Strategies to Improve the Quantum Computation Accuracy for Electrochemical Windows of Ionic Liquids

Rational design of ionic liquids (ILs) with wide electrochemical windows (ECWs) for high-voltage cathodes is essential to elevating the energy density of current rechargeable batteries. It is significant to determine appropriate strategies for accurately predicting the ECWs of ILs. In this study, we compare the calculated ECWs based on three quantum methods, including the highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) method, the ionization potential-electron affinity (IP-EA) method, and the thermodynamic method, under four unique combinations of simulation environments, and assess the discrepancies between the calculated and the experimental results of ECWs. For the three quantum methods, although HOMO-LUMO and IP-EA methods show limited prediction accuracy compared to the experimental values, they can qualitatively rank the anodic limits of anions and the cathodic limits of cations. For the thermodynamic method, we demonstrate that the highest accuracy can be achieved by considering the most rational redox reaction process. By varying the calculation environments, the calculated ECWs tend to be underestimated by considering separate cations and anions of ILs under gas phase, whereas they always exhibit overestimated results when calculating based on the cation-anion pairs. When the computation considers isolated ions with the solvent model plus proper thermodynamic corrections, we observe improved consistency with the experimental results. Though all methods have limitations to achieving perfect predictions of ECWs, we believe rational selection of calculation methods based on application-oriented scenarios can balance the efficiency and accuracy of the results for the development of a high-throughput and accurate database of ECWs for ILs.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.

Intermolecular interactions in clusters of ethylammonium nitrate and 1-amino-1,2,3-triazole

The intermolecular interaction energies, including hydrogen bonds (H-bonds), of clusters of the ionic liquid ethylammonium nitrate (EAN) and 1-amino-1,2,3-triazole (1-AT) based deep eutectic propellants (DeEP) are examined. 1-AT is introduced as a neutral hydrogen bond donor (HBD) to EAN in order to form a eutectic mixture. The effective fragment potential (EFP) is used to examine the bonding interactions in the DeEP clusters. The resolution of the Identity (RI) approximated second order Møller-Plesset perturbation theory (RI-MP2) and coupled cluster theory (RI-CCSD(T)) are used to validate the EFP results. The EFP method predicts that there are significant polarization and charge transfer effects in the EAN:1-AT complexes, along with Coulombic, dispersion and exchange repulsion interactions. The EFP interaction energies are in good agreement with the RI-MP2 and RI-CCSD(T) results. The quasi-atomic orbital (QUAO) bonding and kinetic bond order (KBO) analyses are additionally used to develop a conceptual and semi-quantitative understanding of the H-bonding interactions as a function of the size of the system. The QUAO and KBO analyses suggest that the H-bonds in the examined clusters follow the characteristic hydrogen bonding three-center four electron interactions. The strongest H-bonding interactions between the (EAN)1:(1-AT)n and (EAN)2:(1-AT)n (n = 1-5) complexes are observed internally within EAN; that is, between the ethylammonium cation [EA]+ and the nitrate anion ([NO3]-). The weakest H-bonding interactions occur between [NO3]- and 1-AT. Consequently, the average strengths of the H-bonds within a given (EAN)x:(1-AT)n complex decrease as more 1-AT molecules are introduced into the EAN monomer and EAN dimer. The QUAO bonding analysis suggests that 1-AT in (EAN)x:(1-AT)n can act as both a HBD and a hydrogen bond acceptor simultaneously. It is observed that two 1-AT molecules can form H-bonds to each other. Although the KBOs that correspond to H-bonding interactions in [EA]+:1-AT, [NO3]-:1-AT and between two 1-AT molecules are weaker than the H-bonds in EAN, those weak H-bond networks with 1-AT could be important to form a stable DeEP.

Strongly Bound Polynuclear Anions Comprising Scandium Fluoride Building Blocks

The stability of polynuclear anions composed of ScF3 building blocks was studied by using ab initio and density functional theory electronic structure methods and flexible basis sets. Thorough exploration of ground state potential energy surfaces of (Sc2F7)-, (Sc3F10)-, and (Sc4F13)- anions which may be viewed as comprising ScF3 fragments and the additional fluorine atom led to determining the isomeric structures thereof. It was found that the most stable isomers which are predicted to dominate at room temperature correspond to the compact structures enabling the formation of a large number of Sc-F-Sc bridging linkages rather than to the chain-like structures. The vertical electron detachment energies of the (ScnF3n+1)- anions were found to be very large (spanning the 10.85-12.29 eV range) and increasing with the increasing number of scandium atoms (n) and thus the ScF3 building blocks involved in the structure. Thermodynamic stability of (ScnF3n+1)- anions (i.e., their susceptibility to fragmentation) was also verified and discussed.

Estimation of the Thermophysical Properties of Pentaalkylguanidinium-Based Magnetic Ionic Liquids with Unusual Thermal Expansion Coefficient

The thermal expansion coefficient (αexp. ), the molecular volume (Vm ), the entropy of surface formation (Sa ), and the Gibbs energy of surface formation (Ea ) of four pentaalkylguanidinium-based MILs [Cn TMG][FeCl3 Br] (n=2, 4, 6, 8) were calculated based on the density and surface tension data determined from 278.15 to 323.15 K. In terms of classical semiempirical methods, the standard molar entropy (S0 ), the lattice energy (UPOT ), the molar enthalpy of evaporation (Δ l g H m 0 ( T b ) ${{\rm{\Delta }}_{\rm{l}}^{\rm{g}} H_{\rm{m}}^0 (T_{\rm{b}} )}$ ,Δ l g H m 0 ${{\rm{\Delta }}_{\rm{l}}^{\rm{g}} H_{\rm{m}}^0 }$ (298 K)), and the thermal expansion coefficient (αest. ) of the MILs were further estimated. The estimation results indicate that the classical semiempirical methods are suitable for estimating the thermophysical properties of the MILs, except the unusual αexp. , which were extremely larger than those of representative non-magnetic ionic liquids (ILs). We further optimized the estimation methods and discussed the potential reasons for the unusual thermal expansion coefficient of the MILs.

Toward an extreme-scale electronic structure system

Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree-Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency.

Structural elucidation of polydopamine facilitated by ionic liquid solvation

Minimal understanding of the formation mechanism and structure of polydopamine (pDA) and its natural analogue, eumelanin, impedes the practical application of these versatile polymers and limits our knowledge of the origin of melanoma. The lack of conclusive structural evidence stems from the insolubility of these materials, which has spawned significantly diverse suggestions of pDA's structure in the literature. We discovered that pDA is soluble in certain ionic liquids. Using these ionic liquids (ILs) as solvents, we present an experimental methodology to solvate pDA, enabling us to identify pDA's chemical structure. The resolved pDA structure consists of self-assembled supramolecular aggregates that contribute to the increasing complexity of the polymer. The underlying molecular energetics of pDA solvation and a macroscopic picture of the disruption of the aggregates using IL solvents have been investigated, along with studies of the aggregation mechanism in water.

Purely Predicting the Pharmaceutical Solubility: What to Expect from PC-SAFT and COSMO-RS?

A pair of popular thermodynamic models for pharmaceutical applications, namely, the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state and the conductor-like screening model for real solvents (COSMO-RS) are thoroughly benchmarked for their performance in predicting the solubility of active pharmaceutical ingredients (APIs) in pure solvents. The ultimate goal is to provide an illustration of what to expect from these progressive frameworks when applied to the thermodynamic solubility of APIs based on activity coefficients in a purely predictive regime without specific experimental solubility data (the fusion properties of pure APIs were taken from experiments). While this kind of prediction represents the typical modus operandi of the first-principles-aided COSMO-RS, PC-SAFT is a relatively highly parametrized model that relies on experimental data, against which its pure-substance and binary interaction parameters (k_ij) are fitted. Therefore, to make this benchmark as fair as possible, we omitted any binary parameters of PC-SAFT (i.e., k_ij = 0 in all cases) and preferred pure-substance parameter sets for APIs not trained to experimental solubility data. This computational approach, together with a detailed assessment of the obtained solubility predictions against a large experimental data set, revealed that COSMO-RS convincingly outperformed PC-SAFT both qualitatively (i.e., COSMO-RS was better in solvent ranking) and quantitatively, even though the former is independent of both substance- and mixture-specific experimental data. Regarding quantitative comparison, COSMO-RS outperformed PC-SAFT for 9 of the 10 APIs and for 63% of the API-solvent systems, with root-mean-square deviations of the predicted data from the entire experimental data set being 0.82 and 1.44 log units, respectively. The results were further analyzed to expand the picture of the performance of both models with respect to the individual APIs and solvents. Interestingly, in many cases, both models were found to qualitatively incorrectly predict the direction of deviations from ideality. Furthermore, we examined how the solubility predictions from both models are sensitive to different API parametrizations.

Infrared and Terahertz Spectroscopic Investigation of Imidazolium, Pyridinium, and Tetraalkylammonium Tetrafluoroborate Ionic Liquids

We performed terahertz time-domain spectroscopy and infrared spectroscopy of imidazolium-based, pyridinium-based, and tetraalkylammonium-based tetrafluoroborate ionic liquids to study their characteristic intermolecular and intramolecular vibrational modes to clarify interactions between various cations and the tetrafluoroborate anion. It was found that the central frequency of the intermolecular vibrational band for these ionic liquids has a relatively high frequency, ranging from 90 to 100 cm^-1. In the 900-1150 cm^-1 range, the intramolecular vibrational absorption band of the 3-fold degenerate mode of tetrafluoroborate anions in the ionic liquids was observed. Although the tetrafluoroborate anion is attributable to one of the weakly coordinated anions, the spectroscopic splitting behavior of the 3-fold degenerate mode differs depending on the cation species. It was revealed that the degenerate mode is very sensitive to local interactions between the tetrafluoroborate anion and each cation.

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.

Predicting the Thermodynamics of Ionic Liquids: What to Expect from PC-SAFT and COSMO-RS?

Two popular thermodynamic modeling frameworks, namely, the PC-SAFT equation of state and the COSMO-RS model, are benchmarked for their performance in predicting the thermodynamic properties of pure ionic liquids (ILs) and the solubility of CO₂ in ILs. The ultimate goal is to provide an illustration of what to expect from these frameworks when applied to ILs in a purely predictive way with established parametrization approaches, since the literature generally lacks their mutual comparisons. Two different modeling approaches with respect to the description of the molecular structure of ILs are tested within both models: a cation-anion pair as (i) a single electroneutral supermolecule and (ii) a pair of separately modeled counterions (ion-based approach). In general, we illustrate that special attention should be paid when estimating unknown thermodynamic data of ILs even with these two progressive thermodynamic frameworks. For both PC-SAFT and COSMO-RS, the supermolecule approach generally yields better results for the vapor pressure and the vaporization enthalpy of pure ILs, while the ion-based approach is found to be more suitable for the solubility of CO₂. In spite of some shortcomings, COSMO-RS with the supermolecule approach shows the best overall predictive capabilities for the studied properties. The ion-based strategy within both models has significant limitations in the case of the vaporization properties of ILs. In COSMO-RS, these limitations can, to a certain extent, be surpassed by additional quantum mechanical calculations of the ion pairing in the gas phase, while the ion-based PC-SAFT approach still needs a sophisticated improvement to be developed. As an initiating point, we explore one possible and simple route considering a high degree of cross associations between the counterions in the gas phase.

Structures and properties of ionic crystals and condensed phase ionic liquids predicted with the generalized energy-based fragmentation method

The generalized energy-based fragmentation (GEBF) approach is extended to facilitate ab initio investigations of structures, lattice energies, vibrational spectra and ¹ H NMR chemical shifts of ionic crystals and condensed-phase ionic liquids (ILs) with the periodic boundary conditions (PBC). For selected periodic systems, our results demonstrate that the so-called PBC-GEBF approach can provide satisfactory descriptions on ground-state energies, structures, and vibrational spectra of ionic crystals and IL crystals. The PBC-GEBF approach is then applied to three realistic condensed phase systems. For three ionic crystals (LiCl, NaCl, and KCl), we apply the PBC-GEBF approach with MP2 theory as well as some popular DFT methods to investigate their crystal structures and lattice energies. Our calculations indicate that the crystal structures obtained with PBC-GEBF-MP2/6-311 + G** are very close to the corresponding X-ray structures, while PBC-GEBF-ωB97X-D/6-311 + G** provides satisfactory prediction for crystal structures and lattice energies. For two polymorphs of [n-C₄ mim][Cl] crystals, we find that the PBC-GEBF approach at the M06-2X/6-311 + G** level can give a satisfactory descriptions on structures and Raman spectra of these two crystals. Furthermore, for [C₂ mim][BF₄ ] ILs, we demonstrate that their ¹ H NMR chemical shifts can be estimated from averaging over 5 typical snapshots (extracted from MD simulations) with the PBC-GEBF approach at the B97-2/pcSseg-2 level. The calculated results account for the observed experimental data quite well. Therefore, we expect that the PBC-GEBF approach, combined with various quantum chemistry methods, will become an effective tool in predicting structures and properties of ionic crystals and condensed-phase ILs.

Thermal and Catalytic Decomposition of 2-Hydroxyethylhydrazine and 2-Hydroxyethylhydrazinium Nitrate Ionic Liquid

To develop chemical kinetics models for the combustion of ionic liquid-based monopropellants, identification of the elementary steps in the thermal and catalytic decomposition of components such as 2-hydroxyethylhydrazinium nitrate (HEHN) is needed but is currently not well understood. The first decomposition step in protic ionic liquids such as HEHN is typically the proton transfer from the cation to the anion, resulting in the formation of 2-hydroxyethylhydrazine (HEH) and HNO₃. In the first part of this investigation, the high-temperature thermal decomposition of HEH is probed with flash pyrolysis (<1400 K) and vacuum ultraviolet (10.45 eV) photoionization time-of-flight mass spectrometry (VUV-PI-TOFMS). Next, the investigation into the thermal and catalytic decomposition of HEHN includes two mass spectrometric techniques: (1) tunable VUV-PI-TOFMS (7.4-15 eV) and (2) ambient ionization mass spectrometry utilizing both plasma and laser ionization techniques whereby HEHN is introduced onto a heated inert or iridium catalytic surface and the products are probed. The products can be identified by their masses, their ionization energies, and their collision-induced fragmentation patterns. Formation of product species indicates that catalytic surface recombination is an important reaction process in the decomposition mechanism of HEHN. The products and their possible elementary reaction mechanisms are discussed.

A DLPNO-CCSD(T) benchmarking study of intermolecular interactions of ionic liquids

The accuracy of correlation energy recovered by coupled cluster single-, double-, and perturbative triple-excitations, CCSD(T), has led to the method being considered the gold standard of computational chemistry. The application of CCSD(T) has been limited to medium-sized molecular systems due to its steep scaling with molecular size. The recent development of alternative domain-based local pair natural orbital coupled-cluster method, DLPNO-CCSD(T), has significantly broadened the range of chemical systems to which CCSD(T) level calculations can be applied. Condensed systems such as ionic liquids (ILs) have a large contribution from London dispersion forces of up to 150 kJ mol^-1 in large-scale clusters. Ionic liquids show appreciable charge transfer effects that result in the increased valence orbital delocalization over the entire ionic network, raising the question whether the application of methods based on localized orbitals is reliable for these semi-Coulombic materials. Here the performance of DLPNO-CCSD(T) is validated for the prediction of correlation interaction energies of two data sets incorporating single-ion pairs of protic and aprotic ILs. DLPNO-CCSD(T) produced results within chemical accuracy with tight parameter settings and a non-iterative treatment of triple excitations. To achieve spectroscopic accuracy of 1 kJ mol^-1 , especially for hydrogen-bonded ILs and those containing halides, the DLPNO settings had to be increased by two orders of magnitude and include the iterative treatment of triple excitations, resulting in a 2.5-fold increase in computational cost. Two new sets of parameters are put forward to produce the performance of DLPNO-CCSD(T) within chemical and spectroscopic accuracy.

Origin of low melting point of ionic liquids: dominant role of entropy

Large structural entropy makes salts liquid at room temperature.

Establishing the accuracy of density functional approaches for the description of noncovalent interactions in ionic liquids

We test a number of dispersion corrected versatile Generalized Gradient Approximation (GGA) and meta-GGA functionals for their ability to predict the interactions of ionic liquids, and show that most can achieve energies within 1 kcal mol^-1 of benchmarks. This compares favorably with an accurate dispersion corrected hybrid, ωB97X-V. Our tests also reveal that PBE (Perdew-Burke-Ernzerhof GGA) calculations using the plane-wave projector augmented wave method and Gaussian Type Orbitals (GTOs) differ by less than 0.6 kJ mol^-1 for ionic liquids, despite ions being difficult to evaluate in periodic cells - thus revealing that GTO benchmarks may be used also for plane-wave codes. Finally, the relatively high success of explicit van der Waals density functionals, compared to elemental and ionic dispersion models, suggests that improvements are required for low-cost dispersion correction models of ions.

Improved carbon dioxide absorption in double-charged ionic liquids

Four divalent ionic liquids based on imidazolium cations with alkyl or ether functionalized side-chains were synthesised and characterized: 3,3'-(tetraethyleneglycol-1,11-diyl)bis(1-methyl-1H-imidazolium)bromide, [tetraEG(mim)₂][Br]₂, 3,3'-(tetraethyleneglycol-1,11-diyl)bis(1-methyl-1H-imidazolium)acetate, [tetraEG(mim)₂][OAc]₂, 1-butyl-3-methylimidazolium malonate, [C₄mim]₂[Mal], and 3-butyl-1-methylimidazolium glutarate, [C₄mim]₂[Glut]. Their densities vary between 1.1 and 1.5 g cm^-3 and their viscosities between 0.2 and 4 Pa s at 353 K. We found that the molar volumes are not additive, especially in the case of the divalent ionic liquids based on the double-charged imidazolium cations, meaning that they cannot be predicted using common group contribution methods. The reason for this behaviour could be explained by the structure of the cations, which is dominated by intramolecular hydrogen bonding. The carboxylate-based divalent ionic liquids absorb reversibly large quantities of carbon dioxide following a chemical mechanism described before. An improved 1 : 1 stoichiometry is achieved both in a double-charged imidazolium acetate ionic liquid and in imidazolium carboxylate salts with double charged anions. This behaviour places these ionic liquids amongst the best performing for carbon dioxide absorption.

Does Explicit Polarizability Improve Simulations of Phase Behavior of Ionic Liquids?

Molecular dynamics simulations are performed for a test set of 20 aprotic ionic liquids to investigate whether including an explicit polarizability model in the force field leads to higher accuracy and reliability of the calculated phase behavior properties, especially the enthalpy of fusion. A classical nonpolarizable all-atom optimized potentials for liquid simulations (OPLS) force-field model developed by Canongia Lopes and Pádua (CL&P) serves as a reference level of theory. Polarizability is included either in the form of Drude oscillators, resulting in the CL&P-D models, or in the framework of the atomic multipole optimized energetics for biomolecular application (AMOEBA) force field with polarizable atomic sites. Benchmarking of the calculated fusion enthalpy values against the experimental data reveals that overall the nonpolarizable CL&P model and polarizable CL&P-D models perform similarly with average deviations of about 30%. However, fusion enthalpies from the CL&P-D models exhibit a stronger correlation with their experimental counterparts. The least successful predictions are interestingly obtained from AMOEBA (deviation ca. 60%), which may indicate that a reparametrization of this force-field model is needed to achieve improved predictions of the fusion enthalpy. In general, all FF models tend to underestimate the fusion enthalpies. In addition, quantum chemical calculations are used to compute the electronic cohesive energies of the crystalline phases of the ionic liquids and of the interaction energies within the ion pair. Significant positive correlations are found between the fusion enthalpy and the cohesive energies. The character of the present anions predetermines the magnitude of individual mechanistic components of the interaction energy and related enthalpic and cohesive properties.

Experimental measurement and prediction of ionic liquid ionisation energies

Ionic liquid (IL) valence electronic structure provides key descriptors for understanding and predicting IL properties. The ionisation energies of 60 ILs are measured and the most readily ionised valence state of each IL (the highest occupied molecular orbital, HOMO) is identified using a combination of X-ray photoelectron spectroscopy (XPS) and synchrotron resonant XPS. A structurally diverse range of cations and anions were studied. The cation gave rise to the HOMO for nine of the 60 ILs presented here, meaning it is energetically more favourable to remove an electron from the cation than the anion. The influence of the cation on the anion electronic structure (and vice versa) were established; the electrostatic effects are well understood and demonstrated to be consistently predictable. We used this knowledge to make predictions of both ionisation energy and HOMO identity for a further 516 ILs, providing a very valuable dataset for benchmarking electronic structure calculations and enabling the development of models linking experimental valence electronic structure descriptors to other IL properties, e.g. electrochemical stability. Furthermore, we provide design rules for the prediction of the electronic structure of ILs.

Spectroscopic and Computational Study of Boronium Ionic Liquids and Electrolytes

Boronium cation-based ionic liquids (ILs) have demonstrated high thermal stability and a >5.8 V electrochemical stability window. Additionally, IL-based electrolytes containing the salt LiTFSI have shown stable cycling against the Li metal anode, the "Holy grail" of rechargeable lithium batteries. However, the basic spectroscopic characterisation needed for further development and effective application is missing for these promising ILs and electrolytes. In this work, attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy and density functional theory (DFT) calculations are used in combination to characterise four ILs and electrolytes based on the [NNBH₂ ]⁺ and [(TMEDA)BH₂ ]⁺ boronium cations and the [FSI]^- and [TFSI]^- anions. By using this combined experimental and computational approach, proper understanding of the role of different ion-ion interactions for the Li cation coordination environment in the electrolytes was achieved. Furthermore, the calculated vibrational frequencies assisted in the proper mode assignments for the ILs and in providing insights into the spectroscopic features expected at the interface created when they are adsorbed on a Li(001) surface. A reproducible synthesis procedure for [(TMEDA)BH₂ ]⁺ is also reported. The fundamental findings presented in this work are beneficial for any future studies that utilise IL based electrolytes in next generation Li metal batteries.

A review on machine learning algorithms for the ionic liquid chemical space

There are thousands of papers published every year investigating the properties and possible applications of ionic liquids. Industrial use of these exceptional fluids requires adequate understanding of their physical properties, in order to create the ionic liquid that will optimally suit the application. Computational property prediction arose from the urgent need to minimise the time and cost that would be required to experimentally test different combinations of ions. This review discusses the use of machine learning algorithms as property prediction tools for ionic liquids (either as standalone methods or in conjunction with molecular dynamics simulations), presents common problems of training datasets and proposes ways that could lead to more accurate and efficient models.

Screening Ionic Liquids Based on Ionic Volume and Electrostatic Potential Analyses

Ionic liquids (ILs) are known to have tunable solvation properties, based on the pairing of different anions and cations, but the compositional landscape is vast and challenging to navigate efficiently. Some computational screening protocols are available, but they can be either time-consuming or difficult to implement. Herein, we perform a detailed investigation of the fundamental role of electrostatic interactions in these systems. We effectively develop a bridge between the previous volume-based approach with a quantum structure-property relationship approach to create fast, simple screening guidelines. We propose a new parameter that is applicable to both monovalent and multivalent ions, the ionic polarity index (IPI), which is defined as the ratio of the average electrostatic surface potential (V̅) of the ion to the net charge of the ion (q). The IPI correlation has been tested on a diverse data set of 121 ions, and reliable predictions can be obtained within a homologous series of IL compounds.

Solubility Behavior of CO2 in Ionic Liquids Based on Ionic Polarity Index Analyses

Ionic liquids (ILs) can serve as effective CO₂ solvents with an appropriate selection of different anions and cations. However, due to the large library of potential IL compositions, rapid screening methods are needed for characterizing and ranking the expected properties. We have recently proposed the ionic polarity index (IPI) parameter, effectively connecting volume-based approaches and electrostatic potential analyses and providing a single metric that can potentially be used to rapidly screen for desirable IL properties. In this work, the corresponding anion and cation IPIs are used to generate correlations with respect to the CO₂ volumetric solubility in ILs. The relationships are generally applicable to groups of ILs within a homologous ion series, and this can be particularly valuable for prescreening different ion pairings for maximizing gas solvation performance.

Infrared Spectroscopy in the Middle Frequency Range for Various Imidazolium Ionic Liquids-Common Spectroscopic Characteristics of Vibrational Modes with In-Plane +C(2)-H and +C(4,5)-H Bending Motions and Peak Splitting Behavior Due to Local Symmetry Breaking of Vibrational Modes of the Tetrafluoroborate Anion

Various alkyl-methylimidazolium ionic liquids (ILs) were inspected using infrared spectroscopy in the middle frequency range. In the 1050-1200 cm^-1 range, there is a skeletal vibrational mode accompanied with a large in-plane ⁺C(2)-H bending motion and ⁺C(4)-H and ⁺C(5)-H motions, and in the 1500-1650 cm^-1 range, there are two skeletal vibrational modes with in-plane ⁺C(4,5)-H bending motions. Interestingly, in both ranges, we found that skeletal vibrational modes with a large in-plane ⁺C(2)-H bending motion and in-plane ⁺C(4,5)-H bending motions are insensitive to increases in the basicity of anions or the strengthening of hydrogen bond-type interactions, and the behaviors are completely different from those in the ⁺C-H stretching vibrational modes in the 3000-3200 cm^-1 range and the skeletal vibrational modes with large out-of-plane ⁺C-H motions in the 700-950 cm^-1 range. Furthermore, in alkyl-methylimidazolium tetrafluoroborate [C _n mim⁺][BF₄ ^-] ILs, we found that absorption due to the (threefold) degenerate vibrational mode of [BF₄ ^-] was observed as a broad absorption band with three splitting peaks in the 900-1150 cm^-1 range as a result of local symmetry breaking due to the cation-anion interactions.

Machine learning approaches to understand and predict rate constants for organic processes in mixtures containing ionic liquids

The ability to tailor the constituent ions in ionic liquids (ILs) is highly advantageous as it provides access to solvents with a range of physicochemical properties. However, this benefit also leads to large compositional spaces that need to be explored to optimise systems, often involving time consuming experimental work. The use of machine learning methods is an effective way to gain insight based on existing data, to develop structure-property relationships and to allow the prediction of ionic liquid properties. Here we have applied machine learning models to experimentally determined rate constants of a representative organic process (the reaction of pyridine with benzyl bromide) in IL-acetonitrile mixtures. Multiple linear regression (MLREM) and artificial neural networks (BRANNLP) were both able to model the data well. The MLREM model was able to identify the structural features on the cations and anions that had the greatest effect on the rate constant. Secondly, predictive MLREM and BRANNLP models were developed from the full initial set of rate constant data. From these models, a large number of predictions (>9000) of rate constant were made for mixtures of different ionic liquids, at different proportions of ionic liquid and molecular solvent, at different temperatures. A selection of these predictions were tested experimentally, including through the preparation of novel ionic liquids, with overall good agreement between the predicted and experimental data. This study highlights the benefits of using machine learning methods on kinetic data in ionic liquid mixtures to enable the development of rigorous structure-property relationships across multiple variables simultaneously, and to predict properties of new ILs and experimental conditions.

Predicting Entropic Effects of Water Mixing with Ionic Liquids Containing Anions of Strong Hydrogen Bonding Ability: Role of the Cation

Ionic liquids (ILs) such as choline dihydrogen phosphate exhibit an extraordinary solubilizing ability for proteins such as cytochrome C when mixed with 20 wt % water. Most widely used imidazolium-based ionic liquids coupled with dihydrogen phosphate do not exhibit the same solubilizing properties, suggesting that a multifunctional cation such as choline might play a key role in enhancing these properties of ionic liquid mixtures with water. In this theoretical work, we compare intermolecular interactions between the water molecule and ionic liquid ions in two ion-paired clusters of choline- and 1-butyl-3-methyl-imidazolium-based ionic liquids coupled with acetate, dihydrogen phosphate, and mesylate. Gibbs free energy (GFE) of solvation of water in these ionic liquids was calculated. Incorporation of a water molecule into ionic liquid clusters was accompanied by negative GFEs of solvation in both types of cations. These results were in good agreement with previously reported experimental GFEs of solvation of water in ILs. Compared to imidazolium-based clusters, strong interionic interactions of choline ionic liquids resulted in more negative GFEs due to their smaller deformation upon the addition of a water molecule, with dihydrogen phosphate and mesylate predicting the lowest GFEs of -30.1 and -43.5 kJ/mol^-1, respectively. Lower GFEs of solvation of water in choline-based clusters were also accompanied with smaller entropic penalties, suggesting that water easily incorporates itself into the existing ionic network. Analysis of the intramolecular bonds within the water molecule showed that the choline hydroxyl group donates electron density to the neighboring water molecule, leading to additional polarization. The predicted infrared spectra of clusters of ionic liquids with water showed a pronounced red shift due to strongly polarized O-H bonds, in excellent agreement with the experimentally measured infrared spectra of ionic liquid mixtures with water. Increased polarization of water in choline-based ionic liquids undoubtedly creates more effective solvents for stabilizing biological molecules such as proteins.