Molecular Mechanism of CO2 and SO2 Molecules Binding to the Air/Liquid Interface of 1-Butyl-3-methylimidazolium Tetrafluoroborate Ionic Liquid: A Molecular Dynamics Study with Polarizable Potential Models

Structure, Molecular Interactions, and Dynamics of Aqueous [BMIM][BF4] Mixtures: A Molecular Dynamics Study

Molecular dynamics simulations with many-body polarizable force fields were carried out to investigate the thermodynamic, structural, and dynamic properties of aqueous solutions of 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF₄]). The radial distribution functions exhibit well-defined features, revealing favored structural correlations between [bmim]⁺, [BF₄]^-, and H₂O. The addition of water is shown to alter ionic liquid structural organizations by replacing counterions in the coordination shells and disrupt the cation-anion network. At low water concentration, the majority of water molecules are isolated from each other and have lower average dipole moment than that in pure water. With increasing hydration level, while [bmim][BF₄] ionic network breaks up and becomes isolated ion pairs or free ions in the dilute limit, water begins to form clusters of increasing sizes and eventually forms a percolating network. As a result, the average water dipole moment increases and approaches its bulk value. Water is also observed to have a substantial influence on the dynamics of ionic liquids. At low water content, the cation and anion have similar diffusion coefficients due to the correlated ionic motion of long-lived ion pairs. As the water concentration increases, both ions exhibit greater mobility and faster rotations from the breakup of ionic network. Consequently, the ionic conductivity of [bmim][BF₄] aqueous solutions rises with increasing water composition.

Predicting gas selectivity in organic ionic plastic crystals by free energy calculations

Organic ionic plastic crystals (OIPCs) are molecularly disordered solids, and their potential for the development of gas separation membranes has recently been demonstrated. Here, the gas absorption capability of the OIPC, diethyl(methyl)(isobutyl)phosphonium hexafluorophosphate ([P_122i4][PF₆]), for four gases is predicted through potential of mean force (PMF) calculations based on two methods – average force method and adaptive biasing force method. Both methods correctly predicted the different trends of adsorption and absorption of these gases across the OIPC–gas interface. The distinct energy barriers of the PMF profiles of CO₂ and N₂ near the interface directly reflect the good selectivity of OIPC to these two gases. However, the selectivity of CH₄ and O₂ cannot be accurately reflected by the PMF curve near the interface, because the relative energy varies greatly at different positions inside the OIPC. Thus the average free energy change should be calculated over the entire OIPC box to evaluate the difference in selectivity between the two gases. This also suggests that gas absorption in OIPCs is greatly affected by the structural order and chemical environment. The adaptive biasing force method overall outperforms the average force method. The method should be able to provide a prediction of gas selectivity for a wider range of organic ionic plastic crystals and other solid materials.

The free energy calculation shows the different free energy changes of the adsorption and absorption of gas molecules into an organic ionic plastic crystal, successfully predicting the gas selectivity of this new type of gas separation material.

Evidence of Direct Dissolution of CO2 into the Ionic Liquid [C4min] [NTf2] during Their Initial Interaction

Ionic liquids (ILs) are known for their high ability to capture CO₂. However, the mechanism of CO₂ solubility into ILs during their initial interaction remains controversial. In this study, we analyzed the initial dissolution of CO₂ into an IL 1-butyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)imide ([C₄min] [NTf₂]) by measuring its solubility using a combination of a molecular beam and a flowing liquid jet sheet beam (FJSB) and the King and Wells method (KW method). The temperature dependence of the initial dissolution probability indicates that the solubility of CO₂ in the IL [C₄min] [NTf₂] increases with increasing temperature. This result is not consistent with what has been reported in an equilibrium state. The initial dissolution probability was well-fitted by the Vogel-Fulcher-Tammann (VFT) equation, which describes the dynamical cage structure in IL systems. We also find that the initial dissolution probability was correlated to the cage lifetime and correlation length. The simple model of CO₂ dissolution into an IL with the cage model was implemented to explain the experimental results in this study. Our results indicate that the initial dissolution of CO₂ into the IL corresponds to a direct solution and not an uptake process.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Mass Transfer Thermodynamics through a Gas-Liquid Interface

Molecular level information about thermodynamic variations (enthalpy, entropy, and free energy) of a gas molecule as it crosses a gas-liquid interface is strongly lacking from an experimental perspective under equilibrium conditions. Herein, we perform in situ measurements of water interacting with the ionic liquid (IL) 1-butyl-3-methylimidazolium acetate, [C₄mim][Ace], using ambient pressure X-ray photoelectron spectroscopy in order to assess the interfacial uptake of water quantitatively as a function of temperature, pressure, and water mole fraction ( x_w). The surface spectroscopy results are compared to existing bulk water absorption experiments, showing that the amount of water in the interfacial region is consistently greater than that in the bulk. The enthalpy and entropy of water sorption vary significantly between the gas-liquid interface and the bulk as a function of x_w, with a crossover that occurs near x_w = 0.6 where the water-IL mixture converts from being homogeneous ( x_w < 0.6) to nanostructured ( x_w > 0.6). Free energy results reveal that water at the gas-IL interface is thermodynamically more favorable than that in the bulk, consistent with the enhanced water concentration in the interfacial region. The results herein show that the efficacy for an ionic liquid to absorb a gas phase molecule is not merely a function of bulk solvation parameters but also is significantly influenced by the thermodynamics occurring across the gas-IL interface during the mass transfer process.

Understanding SO2 Capture by Ionic Liquids

Ionic liquids have generated interest for efficient SO2 absorption due to their low vapor pressure and versatility. In this work, a systematic investigation of the structure, thermodynamics, and dynamics of SO2 absorption by ionic liquids has been carried out through quantum chemical calculations and molecular dynamics (MD) simulations. MP2 level calculations of several ion pairs complexed with SO2 reveal its preferential interaction with the anion. Results of condensed phase MD simulations of SO2-IL mixtures manifested the essential role of both cations and anions in the solvation of SO2, where the solute is surrounded by the "cage" formed by the cations (primarily its alkyl tail) through dispersion interactions. These structural effects of gas absorption are substantiated by calculated Gibbs free energy of solvation; the dissolution is demonstrated to be enthalpy driven. The entropic loss of SO2 absorption in ionic liquids with a larger anion such as [NTf2](-) has been quantified and has been attributed to the conformational restriction of the anion imposed by its interaction with SO2. SO2 loading IL decreases its shear viscosity and enhances the electrical conductivity. This systematic study provides a molecular level understanding which can aid the design of task-specific ILs as electrolytes for efficient SO2 absorption.

A Rational Approach to CO2 Capture by Imidazolium Ionic Liquids: Tuning CO2 Solubility by Cation Alkyl Branching

Branching at the alkyl side chain of the imidazolium cation in ionic liquids (ILs) was evaluated towards its effect on carbon dioxide (CO2 ) solubilization at 10 and 80 bar (1 bar=1×10(5)  Pa). By combining high-pressure NMR spectroscopy measurements with molecular dynamics simulations, a full description of the molecular interactions that take place in the IL-CO2 mixtures can be obtained. The introduction of a methyl group has a significant effect on CO2 solubility in comparison with linear or fluorinated analogues. The differences in CO2 solubility arise from differences in liquid organization caused by structural changes in the cation. ILs with branched cations have similar short-range cation-anion orientations as those in ILs with linear side chains, but present differences in the long-range order. The introduction of CO2 does not cause perturbations in the former and benefits from the differences in the latter. Branching at the cation results in sponge-like ILs with enhanced capabilities for CO2 capture.

Water Effect on Acid-Gas Capture Using Choline Lactate: A DFT Insight beyond Molecule-Molecule Pair Simulations

The suitability of CO2 and SO2 capture by using choline lactate ionic liquid as a sorbent and the effect of water content for acid-gas absorption were investigated through density functional theory (DFT) simulations in this work. Simulations that contain model systems considering up to four molecules (cholinium, lactate, water, and CO2/SO2) have been analyzed, and compositional effects on small cluster(s) formed by four ionic pairs and variable number of water molecules have been studied in this work. Assessment of the effect of water content on acid-gas capture that uses exotic ionic liquids is a rare study, and our results showed that water presence hinders CO2/SO2 affinity and solubility dramatically, mainly due to the dominated affinity between the ionic pair and water molecule rather than the CO2/SO2 molecule. Moreover, our studies also showed that affinity between ionic liquid and CO2 is hindered by more than ionic liquid and SO2 rich system with the presence of water in the environment.

SO2 Solvation in the 1-Ethyl-3-Methylimidazolium Thiocyanate Ionic Liquid by Incorporation into the Extended Cation-Anion Network

We have carried out an ab initio molecular dynamics study on the sulfur dioxide (SO₂) solvation in 1-ethyl-3-methylimidazolium thiocyanate for which we have observed that both cations and anions play an essential role in the solvation of SO₂. Whereas, the anions tend to form a thiocyanate- and much less often an isothiocyanate-SO₂ adduct, the cations create a “cage” around SO₂ with those groups of atoms that donate weak interactions like the alkyl hydrogen atoms as well as the heavy atoms of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi $$\end{document}π-system. Despite these similarities between the solvation of SO₂ and CO₂ in ionic liquids, an essential difference was observed with respect to the acidic protons. Whereas CO₂ avoids accepting hydrogen bonds form the acidic hydrogen atoms of the cations, SO₂ can from O(SO₂)–H(cation) hydrogen bonds and thus together with the strong anion-adduct it actively integrates in the hydrogen bond network of this particular ionic liquid. The fact that SO₂ acts in this way was termed a linker effect by us, because the SO₂ can be situated between cation and anion operating as a linker between them. The particular contacts are the H(cation)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\cdots $$\end{document}⋯O(SO₂) hydrogen bond and a S(anion)–S(SO₂) sulfur bridge. Clearly, this observation provides a possible explanation for the question of why the SO₂ solubility in these ionic liquids is so high.

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.

Distributed polarizability models for imidazolium-based ionic liquids

Quantum chemical calculations are used to derive distributed polarizability models sufficiently accurate and compact to be used in classical molecular dynamics simulations of imidazolium-based room temperature ionic liquids. Two distributed polarizability models are fitted to reproduce the induction energy of three imidazolium cations (1,3-dimethyl-, 1-ethyl-3-methyl-, and 1-butyl-3-methylimidazolium) and four anions (tetrafluoroborate, hexafluorophosphate, nitrate, and thiocyanate) polarized by a point charge located successively on a grid of surrounding points. The first model includes charge-flow polarizabilities between first-neighbor atoms and isotropic dipolar polarizability on all atoms (except H), while the second model includes anisotropic dipolar polarizabilities on all atoms (except H). For the imidazolium cations, particular attention is given to the transferability of the distributed polarizability sets. The molecular polarizability and its anisotropy rebuilt by the distributed models are found to be in good agreement with the exact ab initio values for the three cations and 23 additional conformers of 1-ethyl-3-methyl-, 1-butyl-3-methyl-, 1-pentyl-3-methyl-, and 1-hexyl-3-methylimidazolium cations.

Calculation of the intrinsic solvation free energy profile of an ionic penetrant across a liquid-liquid interface with computer simulations

We introduce the novel concept of an intrinsic free energy profile, allowing one to remove the artificial smearing caused by thermal capillary waves, which renders difficulties for the calculation of free energy profiles across fluid interfaces in computer simulations. We apply this concept to the problem of a chloride ion crossing the interface between water and 1,2-dichloroethane and show that the present approach is able to reveal several important features of the free energy profile which are not detected with the usual, nonintrinsic calculations. Thus, in contrast to the nonintrinsic profile, a free energy barrier is found at the aqueous side of the (intrinsic) interface, which is attributed to the formation of a water "finger" the ion pulls with itself upon approaching the organic phase. Further, by the presence of a nonsampled region, the intrinsic free energy profile clearly indicates the coextraction of the first hydration shell water molecules of the ion when entering the organic phase.

Interface behaviors of acetylene and ethylene molecules with 1-butyl-3-methylimidazolium acetate ionic liquid: a combined quantum chemistry calculation and molecular dynamics simulation study

Although imidazolium-based ionic liquids (ILs) combined with oxygen-containing anions were proposed as the potential solvents for the selective separation of acetylene (C(2)H(2)) and ethylene (C(2)H(4)), the detailed mechanism at the molecular level is still not well understood. The present work focuses on a most effective IL for removing C(2)H(2) from a C(2)H(4) stream, 1-butyl-3-methylimidazolium acetate ([BMIM][OAc]), aiming at understanding the first steps of the adsorption process of the molecules at the IL surface. We present a combined quantum mechanical (QM) calculation and molecular dynamics (MD) simulation study on the structure and property of the IL as well as its interaction with C(2)H(2) and C(2)H(4) molecules. The calculated results indicate that C(2)H(2) presents a stronger interaction with the IL than C(2)H(4) and the anion of the IL is mainly responsible for the stronger interaction. QM calculations show a stronger hydrogen-binding linkage between an acidic proton of C(2)H(2)/C(2)H(4) and the basic oxygen atom in [OAc](-) anion, in contrast to the relative weaker association via the C-H···π interaction between C(2)H(2)/C(2)H(4) and the cation. From MD simulations, it is observed that in the interfacial region, the butyl chain of cations and methyl of anions point into the vapor phase. The coming molecules on the IL surface may be initially wrapped by the extensive butyl chain and then devolved to the interface or caught into the bulk by the anion of IL. The introduction of guest molecules significantly influences the anion distribution and orientation on the interface, but the cations are not disturbed because of their larger volume and relatively weaker interaction with the changes in the guest molecules. The theoretical results provide insight into the molecular mechanism of the observed selective separation of C(2)H(2) form a C(2)H(4) stream by ILs.

Influence of the ionic liquid/gas surface on ionic liquid chemistry

Applications such as gas storage, gas separation, NP synthesis and supported ionic liquid phase catalysis depend upon the interaction of different species with the ionic liquid/gas surface. Consequently, these applications cannot proceed to the full extent of their potential without a profound understanding of the surface structure and properties. As a whole, this perspective contains more questions than answers, which demonstrates the current state of the field. Throughout this perspective, crucial questions are posed and a roadmap is proposed to answer these questions. A critical analysis is made of the field of ionic liquid/gas surface structure and properties, and a number of design rules are mined. The effects of ionic additives on the ionic liquid/gas surface structure are presented. A possible driving force for surface formation is discussed that has, to the best of my knowledge, not been postulated in the literature to date. This driving force suggests that for systems composed solely of ions, the rules for surface formation of dilute electrolytes do not apply. The interaction of neutral additives with the ionic liquid/gas surface is discussed. Particular attention is focussed upon H(2)O and CO(2), vital additives for many applications of ionic liquids. Correlations between ionic liquid/gas surface structure and properties, ionic liquid surfaces plus additives, and ionic liquid applications are given.