Dissecting the Vaporization Enthalpies of Ionic Liquids by Exclusively Experimental Methods: Coulomb Interaction,Hydrogen Bonding,and Dispersion Forces

The quantification of hydrogen bonding and dispersion energies from vaporization enthalpies is a great challenge. Dissecting interaction energies is particularly difficult for ionic liquids (ILs), for which the composition of the different types of interactions is known neither for the liquid nor for the gas phase. In this study, we demonstrate the existence of ion pairs in the gas phase and dissect the interaction energies exclusively from measured vaporization enthalpies of different alkylated protic ILs (PILs) and aprotic ILs (AILs) and the molecular analogues of their cations. We demonstrate that the evaporated ion pairs are characterized by H‐bond‐enhanced Coulomb interaction. The overall interaction energy for the ILs in the bulk phase is composed of Coulomb interaction (76 kJ mol^?1), hydrogen bonding (38 kJ mol^?1), and minor dispersion interaction (10 kJ mol^?1). Thus, hydrogen bonding prominently contributes to the overall interaction energy of PILs, which is reflected in the properties of this class of liquids.     

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 provide further 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 inthe 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 of both the ion-ion interaction and the hydrogen bond.     

Dispersion and Hydrogen Bonding Rule: Why the Vaporization Enthalpies of Aprotic Ionic Liquids Are Significantly Larger than those of Protic Ionic liquids

It is well known that gas‐phase experiments and computational methods point to the dominance of dispersion forces in the molecular association of hydrocarbons. Estimates or even quantification of these weak forces are complicated due to solvent effects in solution. The dissection of interaction energies and quantification of dispersion interactions is particularly challenging for polar systems such as ionic liquids (ILs) which are characterized by a subtle balance between Coulomb interactions, hydrogen bonding, and dispersion forces. Here, we have used vaporization enthalpies, far‐infrared spectroscopy, and dispersion‐corrected calculations to dissect the interaction energies between cations and anions in aprotic (AILs), and protic (PILs) ionic liquids. It was found that the higher total interaction energy in PILs results from the strong and directional hydrogen bonds between cation and anion, whereas the larger vaporization enthalpies of AILs clearly arise from increasing dispersion forces between ion pairs.     

Computer‐Aided Design of Ionic Liquids as CO2 Absorbents

Ionic liquids (ILs), vary strongly in their interaction with CO₂. We suggest simple theoretical approach to predict the CO₂ absorption behavior of ILs. Strong interaction of the CO₂ with the IL anions corresponds to chemical absorption whereas weak interaction indicates physical absorption. A predictive estimate with a clear distinction between physical and chemical absorption can be simply obtained according to geometries optimized in the presence of a solvation model instead of optimizing it only in gas phase as has been done to date. The resulting Gibbs free energies compare very well with experimental values and the energies were correlated with experimental capacities. Promising anions, for ionic liquids with reversible CO₂ absorption properties can be defined by a reaction Gibbs free energy of absorption in the range of ?30 to 16 kJ mol^?1.     

Corrected small basis set Hartree‐Fock method for large systems

A quantum chemical method based on a Hartree‐Fock calculation with a small Gaussian AO basis set is presented. Its main area of application is the computation of structures, vibrational frequencies, and noncovalent interaction energies in huge molecular systems. The method is suggested as a partial replacement of semiempirical approaches or density functional theory (DFT) in particular when self‐interaction errors are acute. In order to get accurate results three physically plausible atom pair‐wise correction terms are applied for London dispersion interactions (D3 scheme), basis set superposition error (gCP scheme), and short‐ranged basis set incompleteness effects. In total nine global empirical parameters are used. This so‐called Hartee‐Fock‐3c (HF‐3c) method is tested for geometries of small organic molecules, interaction energies and geometries of noncovalently bound complexes, for supramolecular systems, and protein structures. In the majority of realistic test cases good results approaching large basis set DFT quality are obtained at a tiny fraction of computational cost. © 2013 Wiley Periodicals, Inc.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Ab initio calculations for absorption and emission energies of alkali halide crystals doped with thallium

We calculate transition energies associated with optical properties of thallium doping in alkali halide crystals via an atomic cluster of minimal size where an sp‐valence‐shell impurity enters as a substitutional defect in the model crystal. Hartree–Fock (HF), density functional theory (DFT), and configuration interaction (CI) [CIS (CI with single excitation) and QCISD (single plus double and quadruple excitation)] calculations are performed to theoretically obtain the absorption and emission energies as vertical transitions evaluated at the ground and first excited‐state optimized geometries, respectively, where the optimization is carried out separately with the HF and DFT methods. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 785–790, 2000     

Hydrogen bonding mediated ion pairs of some aprotic ionic liquids and their structural transition in aqueous solution

Ion pair speciation of ionic liquids(ILs) has an important effect on the physical and chemical properties of ILs and recognition of the structure of ion pairs in solution is essential. It has been reported that ion pairs of some ILs can be formed by hydrogen bonding interactions between cations and anions of them. Considering the fact that far-IR(FIR) spectroscopy is a powerful tool in indicating the intermolecular and intramolecular hydrogen bonding, in this work, this spectroscopic technique has been combined with molecular dynamic(MD) simulation and nuclear magnetic resonance hydrogen spectroscopy(~1H NMR) to investigate ion pairs of aprotic ILs [Bmim][NO_3], [BuPy][NO_3], [Pyr_(14)][NO_3], [PP_(14)][NO_3] and [Bu-choline][NO_3] in aqueous IL mixtures. The FIR spectra have been assigned with the aid of density functional theory(DFT) calculations, and the results are used to understand the effect of cationic nature on the structure of ion pairs. It is found that contact ion pairs formed in the neat aprotic ILs by hydrogen bonding interactions between cation and anion, were still maintained in aqueous solutions up to high water mole fraction(say 0.80 for [BuPy][NO3]). When water content was increased to a critical mole fraction of water(say 0.83 for [BuPy][NO3]), the contact ion pairs could be transformed into solvent-separated ion pairs due to the formation of the hydrogen bonding between ions and water. With the further dilution of the aqueous ILs solution, the solvent-separated ion pairs was finally turned into free cations and free anions(fully hydrated cations or anions). The concentrations of the ILs at which the contact ion pairs were transformed into solvent-separated ion pairs and solvent-separated ion pairs were transformed into free ions(fully hydrated ion) were dependent on the cationic structures. These information provides direct spectral evidence for ion pair structures of the aprotic ILs in aqueous solution. MD simulation and ~1H NMR results support the conclusion drawn from FIR spectra investigations.     

关于限域离子液体吸附CO_2的密度泛函理论(DFT)的研究

近年来,温室效应日趋严重,因此吸收CO_2的材料受到了广泛的关注.采用了密度泛函理论(DFT)研究以SiO_2为载体的限域离子液体对CO_2的吸附.对比纯净离子液体(ILs)以及限域离子液体与CO_2的相互作用情况,在这两种状态下两种体系的吸附情况大不相同.从几何结构、相互作用以及电荷分析等方面对ILs、 SiO_2以及ILs/SiO_2复合结构进行研究.计算结果表明,载体、离子液体和CO_2之间都存在较强的相互作用.离子液体的负载不仅改变了SiO_2载体的结构,而且受载体的影响阴阳离子之间的相互作用力也发生了改变.计算结果为进一步深入限域离子液体对CO_2的吸附打下了理论基础.     

The Double‐Faced Nature of Hydrogen Bonding in Hydroxy‐Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations

We characterize the double‐faced nature of hydrogen bonding in hydroxy‐functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H???O and O???O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation–cation interaction is stronger than the cation–anion interaction. We compare the hydrogen‐bond geometries of both “doubly charged hydrogen bonds” with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like‐charge attraction is almost competitive with conventional ion‐pair formation.     

C?H Bond dissociation of acetylene: Local density functional calculations

The C? H bond dissociation energy of acetylene was computed by both ab initio approaches and density functional theory in a local density approximation (DFT–LDA ). Structures and energies for acetylene and its dissociation products (the ethynyl and hydrogen radicals) are presented and compared. Using directly computed HCCH and HCC· energies and the exact H· value, the DFT–LDA calculations are found to yield C? H dissociation energies ranging from 129 to 131 kcal/mol, in good agreement with recent experimental and the highest level theoretical results. The DFT–LDA results show little dependence upon the computational procedure used to obtain geometries.     

咪唑氟硼酸类离子液体熔点与分子内相互作用能的关系

运用密度泛函理论B3LYP方法及6-311++G(d,p)基组对11种咪唑氟硼酸离子液体进行了研究.选择相应化合物的离子体系{[XIM][BF4]n}(n-1)-(n=2,3)作为研究对象,即研究体系由一个烷基咪唑阳离子XIM+和2-3个BF4-阴离子构成,对其进行结构优化.在优化得到的最低能量构型的基础上计算了分子内阳离子与阴离子间的相互作用能,同时考虑了基组重叠误差的修正.结果表明所研究离子体系的离子间相互作用能与离子液体的实验熔点之间存在明确的线性关系,并且所得到的线性方程与氨基酸阳离子型离子液体中存在的线性关系相近.我们的工作为今后借助量子化学方法设计功能化离子液体提供了一定的理论基础.     

气态[HPro（C1）]＋和[Sac]-相互作用及质子转移机制的理论研究

质子化功能离子液体在许多重要领域具有潜在的应用价值,然而人们对其相互作用模式、质子转移行为等方面的认识尚不清晰．本文利用DFT／B3LYP和MP2方法,以脯氨酸阳离子[HPro]＋、脯氨酸甲酯阳离子[HProC1]＋和糖精阴离子[Sac]=组成的质子化离子液体（PILs）为研究对象,探讨气态离子对、分子对、双聚体团簇中的结构单元及其作用模式．利用过渡态和内禀反应坐标（IRC）理论研究气态[HPro]＋[Sac]-和[HProCl]＋[Sac]-中的质子转移反应,AIM（atomsinmolecules）理论分析给出氢键相互作用本质等．计算结果表明,气态单聚体中氢转移能垒很小,体系中存在离子对和分子对的动态平衡．质子转移发生后,体系内部基本作用单元改变,作用强度下降,形成分子对的相互作用能量远远小于离子对的相互作用能量．双聚体团簇计算说明体系中没有质子转移反应发生,在[HPro]2＋[Sac]2-中相互作用的基本结构单元为离子、分子和两性离子,酯化后双聚体[HProC1]2＋[Sac]2-中全部为阴阳离子相互作用．质子转移反应、两性离子和酯化作用等的深入研究对于理解功能化PILs的性质、结构因素及其应用具有一定的实际意义．     

A force field for molecular simulation of tetrabutylphosphonium amino acid ionic liquids

An all-atom force field was set up for a new class of ionic liquids (ILs), tetrabutylphosphonium amino acid, on the basis of the AMBER force field with determining parameters related to the phosphorus atom and modifying several parameters. Ab initio quantum chemical calculations were employed to obtain molecular geometries, infrared frequencies, and torsion energy profiles. Atom partial charges were obtained by using the one-conformation, two-step restraint electrostatic potential approach. Molecular dynamics simulation was carried out in the isothermal-isobaric ensemble for 14 tetrabutylphosphonium amino acid ILs at two temperatures to validate the force field against the experimental densities and heat capacities at constant pressure. Computed thermodynamic properties are in good agreement with available experimental values. Moreover, radial distribution functions were investigated to depict the microscopic structures of these ILs.     

Intermolecular potentials based on symmetry-adapted perturbation theory with dispersion energies from time-dependent density-functional calculations

Recently, three of us have proposed a method [Phys. Rev. Lett. 91, 33201 (2003)] for an accurate calculation of the dispersion energy utilizing frequency-dependent density susceptibilities of monomers obtained from time-dependent density-functional theory (DFT). In the present paper, we report numerical calculations for the helium, neon, water, and carbon dioxide dimers and show that for a wide range of intermonomer separations, including the van der Waals and short-range repulsion regions, the method provides dispersion energies with accuracies comparable to those that can be achieved using the current most sophisticated wave-function methods. If the dispersion energy is combined with (i) the electrostatic and first-order exchange interaction energies as defined in symmetry-adapted perturbation theory (SAPT) but computed using monomer Kohn-Sham (KS) determinants, and (ii) the induction energy computed using the coupled KS static response theory, (iii) the exchange-induction and exchange-dispersion energies computed using KS orbitals and orbital energies, the resulting method, denoted by SAPT(DFT), produces very accurate total interaction potentials. For the helium dimer, the only system with nearly exact benchmark values, SAPT(DFT) reproduces the interaction energy to within about 2% at the minimum and to a similar accuracy for all other distances ranging from the strongly repulsive to the asymptotic region. For the remaining systems investigated by us, the quality of the SAPT(DFT) interaction energies is so high that these energies may actually be more accurate than the best available results obtained with wave-function techniques. At the same time, SAPT(DFT) is much more computationally efficient than any method previously used for calculating the dispersion and other interaction energy components at this level of accuracy.     

Controlling the Subtle Energy Balance in Protic Ionic Liquids: Dispersion Forces Compete with Hydrogen Bonds

The properties of ionic liquids are determined by the energy‐balance between Coulomb‐interaction, hydrogen‐bonding, and dispersion forces. Out of a set of protic ionic liquids (PILs), including trialkylammonium cations and methylsulfonate and triflate anions we could detect the transfer from hydrogen‐bonding to dispersion‐dominated interaction between cation and anion in the PIL [(C₆H₁₃)₃NH][CF₃SO₃]. The characteristic vibrational features for both ion‐pair species can be detected and assigned in the far‐infrared spectra. Our approach gives direct access to the relative strength of hydrogen‐bonding and dispersion forces in a Coulomb‐dominated system. Dispersion‐corrected density functional theory (DFT) calculations support the experimental findings. The dispersion forces could be quantified to contribute about 2.3 kJ mol^?1 per additional methylene group in the alkyl chains of the ammonium cation.     

A sobering assessment of small‐molecule force field methods for low energy conformer predictions

We have carried out a large scale computational investigation to assess the utility of common small‐molecule force fields for computational screening of low energy conformers of typical organic molecules. Using statistical analyses on the energies and relative rankings of up to 250 diverse conformers of 700 different molecular structures, we find that energies from widely used classical force fields (MMFF94, UFF, and GAFF) show unconditionally poor energy and rank correlation with semiempirical (PM7) and Kohn–Sham density functional theory (DFT) energies calculated at PM7 and DFT optimized geometries. In contrast, semiempirical PM7 calculations show significantly better correlation with DFT calculations and generally better geometries. With these results, we make recommendations to more reliably carry out conformer screening.     

Choline saccharinate and choline acesulfamate: ionic liquids with low toxicities

Choline saccharinate and choline acesulfamate are two examples of hydrophilic ionic liquids, which can be prepared from easily available starting materials (choline chloride and a non-nutritive sweetener). The (eco)toxicity of these ionic liquids in aqueous solution is very low in comparison to other types of ionic liquids. A general method for the synthesis and purification of hydrophilic ionic liquids is presented. The method consists of a silver-free metathesis reaction, followed by purification of the ionic liquid by ion-exchange chromatography. The crystal structures show a marked difference in hydrogen bonding between the two ionic liquids, although the saccharinate and the acesulfamate anions show structural similarities. The optimized structures, the energetics, and the charge distribution of cation-anion pairs in the ionic liquids were studied by density functional theory (DFT) and second-order (M?ller-Plesset) perturbation theory calculations. The occupation of the non-Lewis orbitals was considered to obtain a qualitative picture of the Lewis structures. The calculated interaction energies and the dipole moments for the ion pairs in the gas phase were discussed.     

Simulation of electronic and geometric degrees of freedom using a kink-based path integral formulation: application to molecular systems

A kink-based path integral method, previously applied to atomic systems, is modified and used to study molecular systems. The method allows the simultaneous evolution of atomic and electronic degrees of freedom. The results for CH4, NH3, and H2O demonstrate this method to be accurate for both geometries and energies. A comparison with density functional theory (DFT) and second-order Moller-Plesset (MP2) level calculations show the path integral approach to produce energies in close agreement with MP2 energies and geometries in close agreement with both DFT and MP2 results.     

Geometrical and electronic structures of the dication and ion pair in the geminal dicationic ionic liquid 1,3-bis[3-methylimidazolium-yl]propane bromide

Geminal dicationic ionic liquids (ILs), a new category of IL family, have been developed recently and found to possess unique properties compared to conventional monocationic ILs. To establish a basis for understanding their novel properties, we studied the geometrical and electronic structures of the dication ([(mim)C₃(mim)]²⁺) and the ion pair ([(mim)C₃(mim)]²⁺-2Br⁻) in the geminal dicationic IL 1,3-bis[3-methylimidazolium-yl]propane bromide by performing density functional theory calculations. The geometrical structures and relative stabilities for the dication and the ion pair are discussed, and their electronic properties are analyzed in detail. The intrinsic interaction between the dication and Br anions in the most stable conformer was investigated by performing the natural bond orbital analyses. Results for the dication and the ion pair are compared with those of the corresponding monocation ([C₄mim]⁺) and ion pair ([C₄mim]⁺-Br⁻). ¹H NMR spectroscopy for the most stable ion pair has been calculated and the general trend is found to be in fairly agreement with the experimental data.