An ab initio Study of the Structure and Energetics of Hydrogen Bonding in Ionic Liquids

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.

Energy Applications of Ionic Liquids: Recent Developments and Future Prospects

Ionic liquids (ILs) consisting entirely of ions exhibit many fascinating and tunable properties, making them promising functional materials for a large number of energy-related applications. For example, ILs have been employed as electrolytes for electrochemical energy storage and conversion, as heat transfer fluids and phase-change materials for thermal energy transfer and storage, as solvents and/or catalysts for CO2 capture, CO2 conversion, biomass treatment and biofuel extraction, and as high-energy propellants for aerospace applications. This paper provides an extensive overview on the various energy applications of ILs and offers some thinking and viewpoints on the current challenges and emerging opportunities in each area. The basic fundamentals (structures and properties) of ILs are first introduced. Then, motivations and successful applications of ILs in the energy field are concisely outlined. Later, a detailed review of recent representative works in each area is provided. For each application, the role of ILs and their associated benefits are elaborated. Research trends and insights into the selection of ILs to achieve improved performance are analyzed as well. Challenges and future opportunities are pointed out before the paper is concluded.

Proton transfer between sulfonic acids and various propylamines by density functional theory calculations

CONTEXT

Proton transfer in acid-base systems is not well understood. Some acid-base reactions do not proceed to the extent that is expected from the difference in the pK_a values between the base and acid in aqueous solutions, yet some do. In that regard, we have computationally studied the process of proton transfer from the acids of varying strength (benzenesulfonic acid (BSu), methansulfonic acid (MsO), and sulfuric acid (SA)) to the amines with different numbers of propyl substituents on the nitrogen atom (propylamine (PrA), dipropylamine (DPrA), and tripropylamine (TPrA)) upon complexation. Density functional theory calculations were used to thoroughly examine the energetic and structural aspects of the molecular complexes and/or ionic pairs resulting from the acid-base interaction. The potential energy curves along the proton transfer coordinate in these acid-amine systems were analyzed. The change in free energies accompanying the molecular complexes and ionic pair formations was calculated, and the relationship between the energy values and the ΔРА parameter (difference in proton affinity of the acid anion and amine) was established. The larger ΔРА values were found to be unfavorable for the formation of ionic pairs. Using structural, energy, QTAIM, and NBO analyses, we determined that the hydrogen bonds in the molecular complexes PrA-MsO and PrA-BSu are stronger than those in their corresponding ionic pairs. The ionic pairs with the TPrA cation possess the strongest hydrogen bonds of all the ionic pairs being studied, regardless of the anion. The results showed that hydrogen bonding interactions in the molecular complexes contribute significantly to the energies of the acid-base interaction, while in the ionic pairs, the most important energy contribution comes from Coulomb interactions, followed by hydrogen bonding and dispersion forces. The ionic pairs with propylammonium, dipropylammonium, and tripropylammonium cations have stronger ion-ion interactions than tetrapropylammonium (TetPrA)-containing ionic pairs with the same anions. This effect rises with the order of the cation: TetPrA → TPrA → DPrA → PrA, and the sequence of anions is SA → BSu → MsO. The results obtained here expand the concept of acid-base interaction and provide an alternative to experimental searches for suitable acids and bases to obtain new types of protic ionic liquids.

METHODS

All quantum-chemical calculations were carried out by using the DFT/B3LYP-GD3/6-31++G(d,p) level as implemented in the Gaussian 09 software package. For the resulting structures, the electron density distribution was analyzed by the "atoms in molecules" (QTAIM) and the natural bond orbital (NBO) methods on the wave functions obtained at the same level of theory by AIMAll Version 10.05.04 and Gaussian NBO Version 3.1 programs, respectively.

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.

Understanding X-ray Photoelectron Spectra of Ionic Liquids: Experiments and Simulations of 1-Butyl-3-methylimidazolium Thiocyanate

We demonstrate a combined experimental and computational approach to probe the electronic structure and atomic environment of an ionic liquid, based on core level binding energies. The 1-butyl-3-methylimidazolium thiocyanate [C₄C₁Im][SCN] ionic liquid was studied using ab initio molecular dynamics, and results were compared against previously published and new experimental X-ray photoelectron spectroscopy (XPS) data. The long-held assumption that initial-state effects in XPS dominate the measured binding energies is proven correct, which validates the established premise that the ground state electronic structure of the ionic liquid can be inferred directly from XPS measurements. A regression model based upon site electrostatic potentials and intramolecular bond lengths is shown to account accurately for variations in core-level binding energies within the ionic liquid, demonstrating the important effect of long-range interactions on the core levels and throwing into question the validity of traditional single ion pair ionic liquid calculations for interpreting XPS data.

Binary Mixtures of Imidazolium-Based Protic Ionic Liquids. Extended Temperature Range of the Liquid State Keeping High Ionic Conductivities

Binary mixtures based on the two protic ionic liquids 1-ethylimidazolium triflate ([C2HIm][TfO]) and 1-ethylimidazolium bis(trifluoromethanesulfonyl)imide ([C2HIm][TFSI]) have been investigated, with focus on phase behavior, ionic conductivity, and intermolecular interactions as a function of composition (χTFSI indicating the mole fraction of the added compound). It is found that on addition of [C2HIm][TFSI] to [C2HIm][TfO], the melting temperature is first decreased (0 <χ≤ 0.3) and then suppressed (0.3 <χ≤ 0.8) resulting in mixtures with no phase transitions. These mixtures display a wide temperature range of the liquid state and should be interesting for use in devices operating at extreme temperatures. The ionic conductivity does not vary significantly across the composition range analyzed, as evidenced in the comparative Arrhenius plot. The activation energy, Ea, estimated by fitting with the Arrhenius relation in a limited temperature range (between 60 and 140 °C) varies marginally and keeps values between 0.17 and 0.21 eV. These marginal differences can be rationalized by the initially very similar values of the two neat protic ionic liquids. Vibrational spectroscopy, including both Raman and infrared spectroscopies, reveals weakening of the cation–anion interactions for increasing content of [C2HIm][TFSI], which is reflected by the blue shift of the average N-H stretching mode and the red shift of the S-O stretching mode in the TfO anion. These trends correlate with the higher disorder in the mixtures observed by DSC and are evidenced by the decrease and suppression of the melting temperature as the amount of [C2HIm][TFSI] is increased.

Identification of Entry Inhibitors against Delta and Omicron Variants of SARS-CoV-2

Entry inhibitors against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are urgently needed to control the outbreak of coronavirus disease 2019 (COVID-19). This study developed a robust and straightforward assay that detected the molecular interaction between the receptor-binding domain (RBD) of viral spike protein and the angiotensin-converting enzyme 2 (ACE2) receptor in just 10 min. A drug library of 1068 approved compounds was used to screen for SARS-CoV2 entry inhibition, and 9 active drugs were identified as specific pseudovirus entry inhibitors. A plaque reduction neutralization test using authentic SARS-CoV-2 virus in Vero E6 cells confirmed that 2 of these drugs (Etravirine and Dolutegravir) significantly inhibited the infection of SARS-CoV-2. With molecular docking, we showed that both Etravirine and Dolutegravir are preferentially bound to primary ACE2-interacting residues on the RBD domain, implying that these two drug blocks may prohibit the viral attachment of SARS-CoV-2. We compared the neutralizing activities of these entry inhibitors against different pseudoviruses carrying spike proteins from alpha, beta, gamma, and delta variants. Both Etravirine and Dolutegravir showed similar neutralizing activities against different variants, with EC50 values between 4.5 to 5.8 nM for Etravirine and 10.2 to 22.9 nM for Dolutegravir. These data implied that Etravirine and Dolutegravir may serve as general spike inhibitors against dominant viral variants of SARS-CoV-2.

Low-temperature Hg0 abatement by ionic liquid based on weak interaction

Low-temperature gaseous elemental mercury (Hg⁰) abatement is an objective demand in industrial flue gas treatment. In this work, we proposed a new approach for Hg⁰ capture via weak interaction of ionic liquids. Ionic liquids with varied anions (1-butyl-3-methylimidazolium thioacetate ([Bmim][ThAc]), 1-butyl-3-methylimidazolium diethyldithiocarbamate ([Bmim][DTCR]), and 1-butyl-3-methylimidazolium ethylxanthate ([Bmim][EX])) were designed and synthesized. The interaction energies between ionic liquids and elemental mercury were proved to be positively related to mercury removal efficiency, revealing that the electrostatic interaction derived physical adsorption from anions is the dominant factor affecting mercury removal performance. [Bmim][ThAc] with the largest anionic electrostatic interaction energy showed the best mercury abatement performance, achieving a Hg⁰ removal efficiency of over 98% and an adsorption capacity of 10.66 mg/g at 50 °C. The influence of temperature and the results of mercury temperature-programmed desorption (Hg-TPD), X-ray photoelectron spectroscopy (XPS) further confirmed that the ionic liquid combines with elemental mercury through physical adsorption. The work provides a new perspective on designing high-efficiency sorbents for mercury removal at low temperature.

Experimental and Quantum Chemical Studies of Nicotinamide-Oxalic Acid Salt: Hydrogen Bonding, AIM and NBO Analysis

The computational modeling supported with experimental results can explain the overall structural packing by predicting the hydrogen bond interactions present in any cocrystals (active pharmaceutical ingredients + coformer) as well as salts. In this context, the hydrogen bonding synthons, physiochemical properties (chemical reactivity and stability), and drug-likeliness behavior of proposed nicotinamide–oxalic acid (NIC–OXA) salt have been reported by using vibrational spectroscopic signatures (IR and Raman spectra) and quantum chemical calculations. The NIC–OXA salt was prepared by reactive crystallization method. X-ray powder diffraction (XRPD) and differential scanning calorimetry (DSC) techniques were used for the characterization and validation of NIC–OXA salt. The spectroscopic signatures revealed that (N7–H8)/(N23–H24) of the pyridine ring of NIC, (C═O), and (C–O) groups of OXA were forming the intermolecular hydrogen bonding (N–H⋯O–C), (C–H⋯O═C), and (N–H⋯O═C), respectively, in NIC–OXA salt. Additionally, the quantum theory of atoms in molecules (QTAIM) showed that (C10–H22⋯O1) and (C26–H38⋯O4) are two unconventional hydrogen bonds present in NIC–OXA salt. Also, the natural bond orbital analysis was performed to find the charge transfer interactions and revealed the strongest hydrogen bonds (N7–H8⋯O5)/(N23–H24⋯O2) in NIC–OXA salt. The frontier molecular orbital (FMO) analysis suggested more reactivity and less stability of NIC–OXA salt in comparison to NIC–CA cocrystal and NIC. The global and local reactivity descriptors calculated and predicted that NIC–OXA salt is softer than NIC–CA cocrystal and NIC. From MESP of NIC–OXA salt, it is clear that electrophilic (N7–H8)/(N23–H24), (C6═O4)/(C3═O1) and nucleophilic (C10–H22)/(C26–H38), (C6–O5)/(C3–O2) reactive groups in NIC and OXA, respectively, neutralize after the formation of NIC–OXA salt, confirming the presence of hydrogen bonding interactions (N7–H8⋯O5–C6) and (N23–H24⋯O2–C3). Lipinski’s rule was applied to check the activeness of salt as an orally active form. The results shed light on several features of NIC–OXA salt that can further lead to the improvement in the physicochemical properties of NIC.

Inclusion of More Physics Leads to Less Data: Learning the Interaction Energy as a Function of Electron Deformation Density with Limited Training Data

Machine learning (ML) approaches to predicting quantum mechanical (QM) properties have made great strides toward achieving the computational chemist's holy grail of structure-based property prediction. In contrast to direct ML methods, which encode a molecule with only structural information, in this work, we show that QM descriptors improve ML predictions of dimer interaction energy, both in terms of accuracy and data efficiency, by incorporating electronic information into the descriptor. We present the electron deformation density interaction energy machine learning (EDDIE-ML) model, which predicts the interaction energy as a function of Hartree-Fock electron deformation density. We compare its performance with leading direct ML schemes and modern DFT methods for the prediction of interaction energies for dimers of varying charge type, size, and intermolecular separation. Under a low-data regime, EDDIE-ML outperforms other direct ML schemes and is the only model readily transferrable to larger, more complex systems including base pair trimers and porous cages. The underlying physical connection between the density and interaction energy enables EDDIE-ML to reach an accuracy comparable to modern DFT functionals in fewer training data points compared to other ML methods.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

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.

Alkylimidazolium Protic Ionic Liquids: Structural Features and Physicochemical Properties

We focus on a series of protic ionic liquids (PILs) with imidazolium and alkylimidazolium (1R3HIm, R = methyl, ethyl, propyl, and butyl) cations. Using the literature data and our experimental results on the thermal and transport properties, we analyze the effects of the anion nature and the alkyl radical length in the cation structure on the above properties. DFT calculations in gas and solvent phase have resulted in microscopic insights into the structure and cation-anion binding in these PILs. We show that the higher thermodynamic stability of an ion pair raises the PIL decomposition temperature. The melting points of the salts with the same cation decrease as the hydrocarbon radical in the cation becomes longer, which correlates with the weaker ion-ion interaction in the ion pairs. A comparative analysis of the protic ILs and corresponding ILs (1R3MeIm) with the same radical (R) in the cation structure and the same anion has been performed. The lower melting points of the ILs with 1R3MeIm cations are assumed to result from the weakening both of the ion-ion interaction and hydrogen bond.

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.

The explicit role of electron exchange in the hydrogen bonded molecular complexes

We applied a set of advanced bonding descriptors to establish the hidden electron density features and binding energy characteristics of intermolecular DH∙∙∙A hydrogen bonds (OH∙∙∙O, NH∙∙∙O and SH∙∙∙O) in 150 isolated and solvated molecular complexes. The exchange-correlation and Pauli potentials as well as corresponding local one-electron forces allowed us to explicitly ascertain how electron exchange defines the bonding picture in the proximity of the H-bond critical point. The electron density features of DH∙∙∙A interaction are governed by alterations in the electron localization in the H-bond region displaying itself in the exchange hole. At that, they do not depend on the variations in the exchange hole mobility. The electrostatic interaction mainly defines the energy of H-bonds of different types, whereas the strengthening/weakening of H-bonds in complexes with varying substituents depends on the barrier height of the exchange potential near the bond critical point. Energy variations between H-bonds in isolated and solvated systems are also caused the electron exchange peculiarities as follows from the corresponding potential and the interacting quantum atom analyses complemented by electron delocalization index calculations. Our approach is based on the bonding descriptors associated with the characteristics of the observable electron density and can be recommended for in-depth studies of non-covalent bonding.

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.

Cholinium amino acid-based ionic liquids

Boosted by the simplicity of their synthesis and low toxicity, cholinium and amino acid-based ionic liquids have attracted the attention of researchers in many different fields ranging from computational chemistry to electrochemistry and medicine. Among the uncountable IL variations, these substances occupy a space on their own due to their exceptional biocompatibility that stems from being entirely made by metabolic molecular components. These substances have undergone a rather intensive research activity because of the possibility of using them as greener replacements for traditional ionic liquids. We present here a short review in the attempt to provide a compendium of the state-of-the-art scientific research about this special class of ionic liquids based on the combination of amino acid anions and cholinium cations.

Can a gas phase contact ion pair containing a hydrocarbon carbocation be formed in the ground state?

So far, no conclusive evidence of a ground-state contact ion-pair containing a hydrocarbon carbocation has been given in the gas phase.

Multi-Interactions in Ionic Liquids for Natural Product Extraction

Natural products with a variety of pharmacological effects are important sources for commercial drugs, and it is very crucial to develop effective techniques to selectively extract and isolate bioactive natural components from the plants against the background of sustainable development. Ionic liquids (ILs) are a kind of designable material with unique physicochemical properties, including good thermal stability, negligible vapor pressure, good solvation ability, etc. ILs have already been used in pharmaceuticals for extraction, purification, drug delivery, etc. It has been reported that multi-interactions, like hydrogen bonding, hydrophobic interactions, play important roles in the extraction of bioactive components from the plants. In this review, recent progress in the understanding of scientific essence of hydrogen bonding, the special interaction, in ILs was summarized. The extraction of various natural products, one important area in pharmaceutical, by conventional and functional ILs as well as the specific roles of multi-interactions in this process were also reviewed. Moreover, problems existing in bioactive compound extraction by ILs and the future developing trends of this area are given, which might be helpful for scientists, especially beginners, in this field.

Insights into the Stabilization of Fluoride Ions in Ionic Liquids: Pointers to Better Fluorinating Agents

The fluorination efficiency of a fluorinating agent depends on the free availability of the fluoride ions, which in turn depends on its interaction with its solvation shell. A stable fluoride-based poor solvate ionic liquid (SIL) comprising 1-ethyl-3-methylimidazolium (EMIM) cation and ethylene glycol (EG) was recently reported and demonstrated as a fluorinating agent. Herein, we performed ab initio calculations and ab initio molecular dynamics simulations to gain a microscopic understanding of the intermolecular interactions in this SIL in gas, liquid, and crystalline phases. Ethylene glycol (EG), being capable of forming hydrogen bond(s) with the fluoride ion, prevents the latter from reacting with the EMIM cation. Fluoride forms hydrogen bonds with both the cation and the EG molecule, but it was found to have more affinity toward EG, forming a stronger hydrogen bond with its hydroxyl proton than with the acidic proton of the cation. An optimal concentration of EG in the SIL balances its contribution to stabilizing the fluoride ion and yet making fluoride available for fluorination.

Ion Pair Structures and Hydrogen Bonding in RnNH4-n Alkylammonium Ionic Liquids with Hydrogen Sulfate and Mesylate Anions by DFT Computations

Density function theory calculations are employed to study the interaction of amines bearing different numbers of alkyl substituents of different sizes on the nitrogen atom with sulfuric and methanesulfonic acids. The proton affinities of the studied amines are calculated, and it is shown that the higher the value is, the more probable is its protonation. The most stable structures of the ion pairs resulting from the acid-base interaction are obtained and characterized. The geometric parameters of the ion pairs and the characteristics derived from the NBO and QTAIM analysis show that there are hydrogen bonding interactions between the cation and the anion. The hydrogen bonding character of the ion pairs and the strength of the interaction between the ions strongly depend on the nature of the cation itself. The interaction between the ions in the ion pairs weakens with the increase in the cation size. The trend of change in the structural parameters of the H-bonds and energetic characteristics in the cation series for the studied ion pairs is not dependent on the nature of the anion.

Structural Features of Triethylammonium Acetate through Molecular Dynamics

I have explored the structural features and the dynamics of triethylammonium acetate by means of semi-empirical (density functional tight binding, DFTB) molecular dynamics. I find that the results from the present simulations agree with recent experimental determinations with only few minor differences in the structural interpretation. A mixture of triethylamine and acetic acid does not form an ionic liquid, but gives rise to a very complex system where ionization is only a partial process affecting only few molecules (1 over 4 experimentally). I have also found that the few ionic couples are stable and remain mainly embedded inside the AcOH neutral moiety.

Ab Initio Molecular Dynamics Study of Phospho-Amino Acid-Based Ionic Liquids: Formation of Zwitterionic Anions in the Presence of Acidic Side Chains

We present a computational analysis of the complex proton-transfer processes in two protic ionic liquids based on phosphorylated amino acid anions. The structure and the short time dynamics have been analyzed via ab initio and semi-empirical molecular dynamics. Given the presence of mobile protons on the side chain, such ionic liquids may represent a viable prototype of highly conductive ionic mediums. The results of our simulations are not entirely satisfactory in this respect. Our results indicate that conduction in these liquids may be limited due to a quick quenching of the proton-transfer processes. In particular, we have found that, while proton migration does occur on very short timescales, the amino groups act as proton scavengers preventing an efficient proton migration. Despite their limits as conductive mediums, we show that these ionic liquids possess an unconventional microscopic structure, where the anionic component is made by amino acid anions that the aforementioned proton transfer has transformed into zwitterionic isomers. This unusual chemical structure is relevant because of the recent use of amino acid-based ionic liquids, such as CO₂ absorbent.

Unexpected Strong Acidity Enhancing the Effect in Protic Ionic Liquids Quantified by Equilibrium Acidity Studies: A Crucial Role of Cation Structures on Dictating the Solvation Properties

The equilibrium acidities for several series of structural and electronic different organic acids were measured in 3-pyrrolidinium-based aprotic and protic ionic liquid (IL) analogues, that is, [Bmpy][NTf₂], [BpyH][NTf₂], and [PyH₂][NTf₂], by the UV-Vis method. The acidities of neutral acids are found to be much stronger in the protic ILs (PILs) than aprotic ILs (AILs), and the acidifying effect in the two PILs roughly increases proportionally to the number of protons in the cation of the PIL. On the other hand, interestingly, the cationic N⁺-H acids exhibit similar acidities in [Bmpy][NTf₂] and [BpyH][NTf₂] but much weaker than those in [PyH₂][NTf₂]. The Hammett ρ values for the acidic dissociation of para-substituted benzoic acids in two PILs are about the same as that in water (1.00) but significantly smaller than that in the AIL [Bmpy][NTf₂] (2.66). The correlations between the acidities in the PILs and water show double-linear relationships with different slopes and intercepts for the neutral and cationic acids. These, together with previous observations in the PILs [DBUH][OTf] (DBU = 1,5-diazabicyclo(5.4.0)-5-undecene) and EAN (ethylammonium nitrate), clearly indicate that the structure of the cation plays a subtle but a crucial role on sensing the electronic nature of solutes and strongly affects the solvation behaviors of PILs.