Structure and intermolecular interactions in ionic liquid 1-ethyl-3-methylimidazolium bromide and its aqueous solutions investigated by vibrational spectroscopy and quantum chemical computations

The recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions (S. A. Katsyuba, S. Spicher, T. P. Gerasimova, S. Grimme, J. Phys. Chem. B 2020, 124, 6664) is applied to ionic liquid (IL) 1-ethyl-3-methylimidazolium bromide (EmimBr), its C2-deuterated analog [Emim-d]Br and its aqueous solutions. It is shown that the solvation strongly modifies frequencies and IR intensities of the CH/CD stretching vibrations (νCH/νCD) of the imidazolium ring. The main vibrational spectroscopic features of the neat IL are reproduced by the simulations for a cluster (EmimBr)9, in which all three imidazolium CH moieties of the solvated cation form short contacts with three Br- anions, and another two Br- anions are located on top and bottom of imidazolium ring. Cluster models of aqueous solutions reproduce the experimental vibrational frequencies of actual solutions, provided that the Br- anion of solvated contact ion pair (CIP) is situated on top of imidazolium ring, and CH/CD moieties of the latter participate in short contacts with surrounding water molecules. Both structural and spectroscopic analysis allow to interpret the short contacts CH/CD⋯Br- and CH/CD⋯OH2 as hydrogen bonds of approximately equal strength. Enthalpies of bonding of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to ca. 1.4 kcal⋅mol-1, while the analogous estimates obtained for the gas-phase charged species [Emim]2Br+ increase to 5.6 kcal⋅mol-1. It is shown that formation of solvent-shared ion pair (SIP) in aqueous solution, where the counterions of IL are separated by two water molecules H-bonded to a Br- anion, produces frequency shifts ΔνCH/CD, strongly different from the case of CIP formation. This difference can be used for IR/Raman spectroscopic differentiation of the type of solvated ion pairs of EmimBr or other related ILs.

What quantum chemical simulations tell us about the infrared spectra, structure and interionic interactions of a bulk ionic liquid

The recently developed efficient protocol to explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is applied to ionic liquid 1-ethyl-3-methyl-imidazolium tetrafluoroborate [Emim][BF₄], and its C2-deuterated analog [Emim-d][BF₄]. It is shown that the solvation strongly modifies the frequencies and IR intensities of both cationic and anionic components of the ionic liquids. The main features of the bulk spectra are reproduced by the simulations for cluster ([Emim][BF₄])₈, representing an ion pair solvated by the first solvation shell. The geometry of the cluster closely resembles the solid-state structure of the actual ionic liquid and is characterized by short contacts of all CH moieties of the imidazolium ring with [BF₄]^- anions. Both structural and spectroscopic analyses allow the contacts to be interpreted as hydrogen bonds of approximately equal strength. The enthalpies of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to 1.2-1.5 kcal mol^-1, while the analogous estimates obtained for the gas-phase charged species [Emim][BF₄]₂^- and [Emim]₂[BF₄]⁺ increase to 3.6-3.9 kcal mol^-1.

A density functional theory based approach for predicting melting points of ionic liquids

Accurate prediction of melting points of ILs is important both from the fundamental point of view and from the practical perspective for screening ILs with low melting points and broadening their utilization in a wider temperature range. In this work, we present an ab initio approach to calculate melting points of ILs with known crystal structures and illustrate its application for a series of 11 ILs containing imidazolium/pyrrolidinium cations and halide/polyatomic fluoro-containing anions. The melting point is determined as a temperature at which the Gibbs free energy of fusion is zero. The Gibbs free energy of fusion can be expressed through the use of the Born-Fajans-Haber cycle via the lattice free energy of forming a solid IL from gaseous phase ions and the sum of the solvation free energies of ions comprising IL. Dispersion-corrected density functional theory (DFT) involving (semi)local (PBE-D3) and hybrid exchange-correlation (HSE06-D3) functionals is applied to estimate the lattice enthalpy, entropy, and free energy. The ions solvation free energies are calculated with the SMD-generic-IL solvation model at the M06-2X/6-31+G(d) level of theory under standard conditions. The melting points of ILs computed with the HSE06-D3 functional are in good agreement with the experimental data, with a mean absolute error of 30.5 K and a mean relative error of 8.5%. The model is capable of accurately reproducing the trends in melting points upon variation of alkyl substituents in organic cations and replacement one anion by another. The results verify that the lattice energies of ILs containing polyatomic fluoro-containing anions can be approximated reasonably well using the volume-based thermodynamic approach. However, there is no correlation of the computed lattice energies with molecular volume for ILs containing halide anions. Moreover, entropies of solid ILs follow two different linear relationships with molecular volume for halides and polyatomic fluoro-containing anions. Continuous progress in predicting crystal structures of organic salts with halide anions will be a key factor for successful prediction of melting points with no prior knowledge of the crystal structure.

Basicity Characterization of Imidazolyl Ionic Liquids and Their Application for Biomass Dissolution

Alkalinity determination is of crucial significance for the applications of basic ionic liquids with imidazolyl. In this work, the ionization constant pKb value and acid function H- values of ionic liquids synthesized were calculated by pH method and UV spectrum-Hammett method. The dissolution ratio of biomass in these ionic liquids was measured at different temperatures. Finally, the relationship between the alkalinity and structure of these ionic liquids was discussed, and the relationship between the alkalinity of ionic liquid and the dissolution mechanism biomass was also discussed. The results show that the basicity of carboxylate ionic liquids is determined mainly by their anions, whereas cations take some finely tuned roles. Furthermore, cations and anions are equally important and are involved in dissolution mechanisms.

Ionic liquids as solvents of polar and non-polar solutes: affinity and coordination

Evolution of H₂O and CO₂ interactions with an ionic liquid (IL) from gas phase to IL phase is described. Affinity is lost and coordination patterns vary in the process, favouring H₂O–anion and CO₂–cation interactions.

New di-ferrocenyl-ethynylpyridinyl triphenylphosphine copper halide complexes and related di-ferricenyl electro-crystallized materials

Three di-ferrocenyl-ethynylpyridinyl copper complexes have been synthesised and CV measurements made.

How strong is hydrogen bonding in ionic liquids? Combined X-ray crystallographic, infrared/Raman spectroscopic, and density functional theory study

Hydrogen bonding in ionic liquids based on the 1-(2'-hydroxylethyl)-3-methylimidazolium cation ([C₂OHmim](+)) and various anions ([A](-)) of differing H-bond acceptor strength, viz. hexafluorophosphate [PF6](-), tetrafluoroborate [BF₄](-), bis(trifluoromethanesulfonimide) [Tf₂N](-), trifluoromethylsulfonate [OTf](-), and trifluoroacetate [TFA](-), was studied by a range of spectroscopic and computational techniques and, in the case of [C₂OHmim][PF6], by single crystal X-ray diffraction. The first quantitative estimates of the energy (E(HB)) and the enthalpy (-ΔH(HB)) of H-bonds in bulk ILs were obtained from a theoretical analysis of the solid-state electron-density map of crystalline [C₂OHmim][PF6] and an analysis of the IR spectra in crystal and liquid samples. E(HB) for OH···[PF6](-) H-bonds amounts to ~3.4-3.8 kcal·mol(-1), whereas weaker H-bonds (2.8-3.1 kcal·mol(-1)) are formed between aromatic C2H group of imidazolium ring and the [PF6](-) anion. The enthalpy of the OH···[A](-) H-bonds follows the order: [PF6] (2.4 kcal·mol(-1)) < [BF₄] (3.3 kcal·mol(-1)) < [Tf₂N] (3.4 kcal·mol(-1)) < [OTf] (4.7 kcal·mol(-1)l) < [TFA] (6.2 kcal·mol(-1)). The formation of aggregates of self-associated [C₂OHmim](+) cations is present in liquid [C₂OHmim][PF6], [C₂OHmim][BF₄], and [C₂OHmim][Tf₂N], with the energy of the OH···OH H-bonds amounting to ~6 kcal·mol(-1). Multiple secondary interactions in the bulk ILs influence their structure, vibrational spectra, and H-bond strength. In particular, these interactions can blue-shift the stretching frequencies of the CH groups of the imidazolium ring in spite of red-shifting CH···[A](-) H-bonds. They also weaken the H-bonding in the IL relative to the isolated ion pairs, with these anticooperative effects amounting to ca. 50% of the E(HB) value.