Prediction of macroscopic properties of protic ionic liquids by ab initio calculations

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.

Phase Transitions and Electrochemical Properties of Ionic Liquids and Ionic Liquid-Solvent Mixtures

Recent advances in studies of ionic liquids (IL) and ionic liquid-solvent mixtures are reviewed. Selected experimental, simulation, and theoretical results for electrochemical, thermodynamical, and structural properties of IL and IL-solvent mixtures are described. Special attention is paid to phenomena that are not predicted by the classical theories of the electrical double layer or disagree strongly with these theories. We focus on structural properties, especially on distribution of ions near electrodes, on electrical double layer capacitance, on effects of confinement, including decay length of a dissjoining pressure between confinig plates, and on demixing phase transition. In particular, effects of the demixing phase transition on electrochemical properties of ionic liquid-solvent mixtures for different degrees of confinement are presented.

Hydrogen Bond Kinetics, Ionic Dynamics, and Voids in the Binary Mixtures of Protic Ionic Liquids with Alkanolamines

Classical molecular dynamics simulations were used to investigate the structural and dynamical properties of the mixtures of ionic liquids (ILs) with the conjugate forms of the cation in a 1:1 molar ratio. The experimental studies suggested the combination of ethanolamines and ILs as novel absorbents for acidic gases such as CO₂ and H₂S, which provide the advantage of efficient absorption of gases at low pressures. However, the microscopic properties of the ionic mixtures are not studied. From our computational investigations, the densities of mixtures are reported and compared with the experimental results. The structural evolution of mixtures is reported by radial distribution functions, coordination numbers, void analysis, and spatial distribution functions. The mixtures' dynamic properties were studied by analyzing the hydrogen bond, ion-pair, and ion-cage lifetimes of the system. Monoethanolammonium and triethanolammonium ILs show different types of spatial distribution functions. The cations have lesser effect on dynamics compared with anions. The charge on the anion greatly affects the dynamics of mixtures. The dianion mixtures show slower dynamics than the monoanionic mixtures. The hydrogen bonding between cations and anions is stronger than that between cations and neutral molecules due to strong coulombic attractive forces. The cations spend more time around the dianions as compared to monoanions. The distributions of voids show that the void sizes are smaller in triethanolamine-based mixtures. The sulfobenzoate-based mixtures show voids smaller than those of pyridine-3-carboxylate-based mixtures due to more available free space between the entities, which facilitates the overall dynamics.

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.

Review of Chitosan-Based Polymers as Proton Exchange Membranes and Roles of Chitosan-Supported Ionic Liquids

Perfluorosulphonic acid-based membranes such as Nafion are widely used in fuel cell applications. However, these membranes have several drawbacks, including high expense, non-eco-friendliness, and low proton conductivity under anhydrous conditions. Biopolymer-based membranes, such as chitosan (CS), cellulose, and carrageenan, are popular. They have been introduced and are being studied as alternative materials for enhancing fuel cell performance, because they are environmentally friendly and economical. Modifications that will enhance the proton conductivity of biopolymer-based membranes have been performed. Ionic liquids, which are good electrolytes, are studied for their potential to improve the ionic conductivity and thermal stability of fuel cell applications. This review summarizes the development and evolution of CS biopolymer-based membranes and ionic liquids in fuel cell applications over the past decade. It also focuses on the improved performances of fuel cell applications using biopolymer-based membranes and ionic liquids as promising clean energy.

Modeling protic dicationic ionic liquids based on quaternary ammonium, imidazolium or pyrrolidinium cations and bis(trifluoromethanesulfonyl)imide anion: Structure and spectral characteristics

Protic dicationic ionic liquids (PDILs) have attracted growing attention owing to their applications in domains of electrochemistry, proton conducting materials and other diverse areas. In the present work protic dicationic ionic liquids (PDILs) comprising of quaternary ammonium-, imidazolium- or pyrrolidinium-dications and bis(trifluoromethanesulfonyl)imide (Tf₂N‾) anion have been modelled as the dication-(Tf₂N)₂ complexes. Electronic structure, vibrational and ¹H NMR spectra of these complexes have been derived employing the M06-2x density functional theory. Theoretical calculations have shown that the strength of cation-anion binding follows the order: methylpyrrolidinum > quaternary ammonium > butylpyrrolidinium > imidazolium, which can be attributed to number and strength of N-H---O and C-H---O interactions. The dication-(Tf₂N)₂ complexes emerge with signature in frequency up-shift of the characteristic N-H stretching in their infrared spectra. Underlying molecular interactions are unveiled through natural bond orbital analyses, Quantum theory of atoms in molecules (QTAIM) and noncovalent interaction reduced density gradient method. The calculations have shown that cation-anion binding energies increase linearly with kinetic energy density component G(r) in QTAIM analysis and proton affinities in the PDILs. A correlation between change in free energies accompanying the dication-(Tf₂N)₂ complexes and proton affinities has also been established.

Crystal Structure and Properties of Imidazo-Pyridine Ionic Liquids

Computational studies were performed on novel protic ionic liquids imidazolium-[1,2-a]-pyridine trifluoroacetate [ImPr][TFA] synthesized by the reaction of imidazo-[1,2a]-pyridine (ImPr) with trifluoroacetic acid (TFA), and on fused salt imidazolium-[1,2-a]-pyridine maleamic carbonate [ImPr][Mal] synthesized by reaction of ImPr with maleamic acid (Mal). Synthesis was performed as one-pot reactions, which applies green chemistry tenets. Both these compounds begin to decompose at 180°C. Our computational studies suggest another thermal reaction channel, in which [ImPr][Mal] can also thermally polymerizes to polyacrylamide which then cyclizes. This is thermal product remains stable up to 700 degrees, consistent with our thermogravimetric studies. [ImPr][TFA] exhibited good conductivity and ideal ionic behavior, as evaluated by a Walden plot. X-ray crystallography of [ImPr][TFA] revealed a tightly packed system for the crystals as a result of strong ionic interaction, pi-stacking, and fluorine-CH interactions. Both synthesized compounds exhibited some CO₂ absorptivity, with [ImPr][Mal] outperforming [ImPr][TFA] in this regard. The quantum chemistry based computational methods can shed light on many properties of these ionic liquids, but they are challenged in fully describing their ionic nature. © 2017 Wiley Periodicals, Inc.

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.

Amphoteric water as acid and base for protic ionic liquids and their electrochemical activity when used as fuel cell electrolytes

Amphoteric water was mixed with equimolar amounts of a super-strong acid, trifluoromethanesulfonic acid (TfOH), and a super-strong base, 1,8-diazabicyclo[5.4.0]-7-undecene (DBU). Bulk physicochemical and electrochemical properties of the mixtures were compared with those of the best ever reported protic ionic liquid (PIL), diethylmethylammonium trifluoromethanesulfonate ([dema][TfO]), which has excellent physicochemical properties as a fuel cell electrolyte. The acidic mixture ([H₃O][TfO]) behaved as a protic ionic liquid, while the basic mixture ([DBU]OH) showed incomplete proton transfer. The Walden plot indicated that [H₃O][TfO] behaves as a good PIL, similar to [dema][TfO], whereas [DBU]OH behaves as a poor PIL. [H₃O][TfO] showed excellent H₂/O₂ fuel cell performance at 80 °C; however, the performance deteriorated as the bulk water content increased, because of the retardation of the electrode kinetics due to the oxidation of Pt in the presence of bulk water. On the other hand, [DBU]OH exhibited very poor performance possibly because of the existence of neutral species in the system.

Protic ammonium carboxylate ionic liquids: insight into structure, dynamics and thermophysical properties by alkyl group functionalization

This study is aimed at characterising the structure, dynamics and thermophysical properties of five alkylammonium carboxylate ionic liquids (ILs) from classical molecular dynamics simulations.

A Reinvestigation of the Ionic Liquid Diisopropylethylammonium Formate by NMR and DFT Methods

The complex between diisopropylethylamine (DIPEA) and formic acid has been reinvestigated. Mixing the compounds in the ratio 1:1 leads to a phase separation in which the upper phase is DIPEA and the lower phase is the "ionic liquid" named DIPEF. A combined NMR and DFT study shows that the lower phase primarily is formic acid:formate and diisopropylammonium ions in the ratio 2:1 (acid:base) plus the formic acid dimer. Addition of more acid leads to more and more of the acid dimer. The proton transfer in the system is 65-80%. The structural picture presented in this paper is very different from that presented elsewhere. However, the present picture should be considered using acids and bases with a pK_a difference less than 8. The formic acid content in the DIPEF ionic liquid causes desorption of the dye-sensitized solar cell (DSC) dye N719 from the photo anode, and DIPEF is therefore not a suitable electrolyte for DSCs.

Brønsted acidity of protic ionic liquids: a modern ab initio valence bond theory perspective

Modern ab initio valence bond theory calculations of representative primary and tertiary ammonium protic ionic liquids indicate that modern ab initio valence bond theory can be employed to assess the acidity and ionicity of protic ionic liquids a priori.

Nanostructure of mixtures of protic ionic liquids and lithium salts: effect of alkyl chain length

The bulk structure of mixtures of two protic ionic liquids, propylammonium nitrate and butylammonium nitrate, with a salt with a common anion, is analyzed using small angle X-ray scattering and classical molecular dynamics simulations.

Solvation of lithium salts in protic ionic liquids: a molecular dynamics study

The structure of solutions of lithium nitrate in a protic ionic liquid with a common anion, ethylammonium nitrate, at room temperature is investigated by means of molecular dynamics simulations. Several structural properties, such as density, radial distribution functions, hydrogen bonds, spatial distribution functions, and coordination numbers, are analyzed in order to get a picture of the solvation of lithium cations in this hydrogen-bonded, amphiphilically nanostructured environment. The results reveal that the ionic liquid mainly retains its structure upon salt addition, the interaction between the ammonium group of the cation and the nitrate anion being only slightly perturbed by the addition of the salt. Lithium cations are solvated by embedding them in the polar nanodomains of the solution formed by the anions, where they coordinate with the latter in a solid-like fashion reminiscent of a pseudolattice structure. Furthermore, it is shown that the average coordination number of [Li](+) with the anions is 4, nitrate coordinating [Li](+) in both monodentate and bidentate ways, and that in the second coordination layer both ethylammonium cations and other lithiums are also found. Additionally, the rattling motion of lithium ions inside the cages formed by their neighboring anions, indicative of the so-called caging effect, is confirmed by the analysis of the [Li](+) velocity autocorrelation functions. The overall picture indicates that the solvation of [Li](+) cations in this amphiphilically nanostructured environment takes place by means of a sort of inhomogeneous nanostructural solvation, which we could refer to as nanostructured solvation, and which could be a universal solvation mechanism in ionic liquids.

Probing molecular interaction in ionic liquids by low frequency spectroscopy: Coulomb energy, hydrogen bonding and dispersion forces

Low vibrational spectroscopy provides detailed information on the strength and type of interaction and their influence on the properties of ionic liquids.

Dual mode of extraction for Cs+ and Na+ ions with dicyclohexano-18-crown-6 and bis(2-propyloxy)calix[4]crown-6 in ionic liquids: density functional theoretical investigation

Experimentally observed distribution constant and theoretically predicted values of ΔG_ext for Cs⁺ and Na⁺ ions with DCH18C6 and BPC6 ligand in ionic liquid and octanol.

Structural heterogeneity and unique distorted hydrogen bonding in primary ammonium nitrate ionic liquids studied by high-energy X-ray diffraction experiments and MD simulations

Liquid structure and the closest ion-ion interactions in a series of primary alkylammonium nitrate ionic liquids [C(n)Am(+)][NO(3)(-)] (n = 2, 3, and 4) were studied by means of high-energy X-ray diffraction (HEXRD) experiments with the aid of molecular dynamics (MD) simulations. Experimental density and X-ray structure factors are in good accordance with those evaluated with MD simulations. With regard to liquid structure, characteristic peaks appeared in the low Q (Q: a scattering vector) region of X-ray structure factors S(Q)'s for all ionic liquids studied here, and they increased in intensity with a peak position shift toward the lower Q side by increasing the alkyl chain length. Experimentally evaluated S(Q(peak))(r(max)) functions, which represent the S(Q) intensity at a peak position of maximum intensity Q(peak) as a function of distance (actually a integration range r(max)), revealed that characteristic peaks in the low Q region are related to the intermolecular anion-anion correlation decrease in the r range of 10-12 Å. Appearance of the peak in the low Q region is probably related to the exclusion of the correlations among ions of the same sign in this r range by the alkyl chain aggregation. From MD simulations, we found unique and rather distorted NH···O hydrogen bonding between C(n)Am(+) (n = 2, 3, and 4) and NO(3)(-) in these ionic liquids regardless of the alkyl chain length. Subsequent ab initio calculations for both a molecular complex C(2)H(5)NH(2)···HONO(2) and an ion pair C(2)H(5)NH(3)(+)···ONO(2)(-) revealed that such distorted hydrogen bonding is specific in a liquid state of this family of ionic liquids, though the linear orientation is preferred for both the N···HO hydrogen bonding in a molecular complex and the NH···O one in an ion pair. Finally, we propose our interpretation of structural heterogeneity in PILs and also in APILs.

The influence of hydrogen bonding on the physical properties of ionic liquids

Potential applications of ionic liquids depend on the properties of this class of liquid material. To a large extent the structure and properties of these Coulomb systems are determined by the intermolecular interactions among anions and cations. In particular the subtle balance between Coulomb forces, hydrogen bonds and dispersion forces is of great importance for the understanding of ionic liquids. The purpose of the present paper is to answer three questions: Do hydrogen bonds exist in these Coulomb fluids? To what extent do hydrogen bonds contribute to the overall interaction between anions and cations? And finally, are hydrogen bonds important for the physical properties of ionic liquids? All these questions are addressed by using a suitable combination of experimental and theoretical methods including newly synthesized imidazolium-based ionic liquids, far infrared spectroscopy, terahertz spectroscopy, DFT calculations, differential scanning calorimetry (DSC), viscometry and quartz-crystal-microbalance measurements. The key statement is that although ionic liquids consist solely of anions and cations and Coulomb forces are the dominating interaction, local and directional interaction such as hydrogen bonding has significant influence on the structure and properties of ionic liquids. This is demonstrated for the case of melting points, viscosities and enthalpies of vaporization. As a consequence, a variety of important properties can be tuned towards a larger working temperature range, finally expanding the range of potential applications.

Towards large-scale, fully ab initio calculations of ionic liquids

Ionic liquids have attracted a substantial amount of interest as replacement of traditional electrolytes in high efficiency electrochemical devices for generation and storage of energy due to their superior physical and chemical properties, especially low volatility and high electrochemical stability. For enhanced performance of the electrochemical devices ionic liquids are required to be highly conductive and low viscous. Long-range Coulomb and short-range dispersion interactions between ions affect physical and chemical properties of ionic liquids in a very complex way, thus preventing direct correlations to the chemical structure. Considering a vast combination of available cations and anions that can be used to synthesize ionic liquids, development of predictive theoretical approaches that allow for accurate tailoring of their physical properties has become crucial to further enhance the performance of electrochemical devices such as lithium batteries, fuel and solar cells. This perspective article gives a thorough overview of current theoretical approaches applied for studying thermodynamic (melting point and enthalpy of vapourisation) and transport (conductivity and viscosity) properties of ionic liquids, emphasizing their reliability and limitations. Strategies for improving predictive power and versatility of existing theoretical approaches are also outlined.

From salts to ionic liquids by systematic structural modifications: a rational approach towards the efficient modular synthesis of enantiopure imidazolium salts

This paper reports a simple and robust modular synthetic strategy that leads to a large variety of configurationally and structurally diverse imidazole-based chiral ionic liquids (CILs) by lipase-catalyzed resolution. The intimate microscopic interactions of the supramolecular ionic network of these imidazolium chiral salts at the molecular level are investigated both spectroscopically (NMR, FT-IR-ATR) and theoretically, and a topological analysis of the experimental electron densities obtained by X-ray diffraction of single crystals is performed. Our results support the key role played by the relative configuration of the -OR group on the hydrogen-bonding pattern and its strong influence on the final physical properties of the imidazolium salt. We also obtained a reasonable correlation between the observed melting point and the non-covalent interactions. The spectroscopic data and the topological analysis reflect the key role played by hydrogen bonds between the OH and imidazolium C2H groups in both cation-anion and cation-cation interactions, with the presence of an OH group leading to an additional inter-cation interaction. This interaction significantly affects the properties of stereoisomeric salts. Even more interestingly, we also studied the effect of the chirality by comparing enantiopure CILs with their racemic mixtures and found that, with the exception of trans-Cy6-OH-Im-Bn-Br, the melting points of the racemic mixtures are higher than those of the corresponding enantiomerically pure forms. For stereoisomeric examples, we have successfully explained the differences in melting temperatures in light of the corresponding structural data. Chirality should therefore be taken into account as a highly attractive design vector in the preparation of ILs with specifically desired properties.

Imidazolium-based ionic liquids. 1-methyl imidazolium nitrate: thermochemical measurements and ab initio calculations

In this work data of the molar enthalpies of formation of the ionic liquid 1-methylimidazolium nitrate [H-MIM][NO3] was measured by means of combustion calorimetry. The molar enthalpy of fusion of [H-MIM][NO3] was measured using DSC. Experiments to vaporize the ionic liquid into vacuum or nitrogen stream in order to obtain vaporization enthalpy have been performed. Ab initio calculations of the enthalpy of formation in the gaseous phase have been performed for the ionic species using the G3MP2 theory. The combination of traditional combustion calorimertry with modern high-level ab initio calculations allow the determination of the molar enthalpy of vaporization of the ionic liquid under study. The ab initio calculations indicate that [H-MIM][NO3] is most probably separated into the neutral species methyl-imidazole and HNO3 in the gaseous phase at conditions of the vaporization experiments.