Why are ionic liquid ions mainly associated in water? A Car–Parrinello study of 1-ethyl-3-methyl-imidazolium chloride water mixture

Comparative assessment of the performance of density functionals and dispersion correction on different properties of dicationic ionic liquids - an ab initio molecular dynamics (AIMD) study

In this study, we investigated the effect of DFT density functionals and dispersion correction on an imidazolium-based dicationic ionic liquid (DIL) using ab initio molecular dynamics simulations. To achieve this purpose, the electronic structures, as well as the structural and dynamical properties of [C3(mim)2][NTF2]2 DIL, were obtained using the BLYP and PBE functionals, both with and without D3-correction, and the results were compared with experimental values. Radial distribution functions and structure factors revealed that applying D3-correction increases the interaction between the anion and hydrogen atoms of the rings and side chains. The simulation of the studied DIL with the BLYP-D3 functional depicted lower structural heterogeneity compared to the other functionals. Analysis of Voronoi tessellation and linkage chain conformations showed a reduction in the aggregation of the linkage alkyl chains in the presence of D3-correction, which is more pronounced in the BLYP functional than in PBE. Additionally, it was observed that the probability of forming a hydrogen-bond network depends on both the type of used density functionals and applying dispersion correction. The results of dynamical properties, such as the self-diffusion coefficients, velocity autocorrelation function, and van Hove correlation function, as well as ion pair, ion cage, and hydrogen bond dynamics, indicated that applying D3-correction in both density functionals leads to an increase in the dynamics of the studied DIL. Additionally, the ratio of self-diffusion coefficients of the anion to the cation in the BLYP functional is closer to experimental values compared to the PBE functional. Furthermore, the electronic structure, including dipole moment distribution, and also infrared (IR) and power spectra were studied. Applying D3-correction and the type of density functionals have a significant effect on the dipole moment distribution of ions. Moreover, the results of IR and power spectra demonstrated that only in the BLYP functional, by applying D3-correction, the hydrogen bonding between the anion and the hydrogen atoms of the cation is strengthened at high wavenumbers. Thus, we conclude that applying D3 correction to both the BLYP and PBE density functionals improves the accuracy in describing the various properties of the studied system. Overall, the evaluation of different structural, dynamical, and vibrational properties of [C3(mim)2][NTF2]2 DIL suggests that the BLYP-D3 density functional may be the best choice among the studied density functionals.

Ionic Pairing and Selective Solvation of Butylmethylimidazolium Chloride Ion Pairs in DMSO-Water Mixtures: A Comprehensive Examination via Molecular Dynamics Simulations and Potentials of Mean Force Analysis

Ionic liquids (ILs) with dimethyl sulfoxide (DMSO) and water act as a promising solvent medium for the dissolution of cellulose in an efficient manner. To develop a proper solvent system, it is really important to understand the thermodynamics of the molecular solutions consisting of ILs, DMSO, and water. The ion-pairing propensity of the ILs in the presence of DMSO and water plays a crucial role in governing the property of the solvent mixtures. Employing all-atom molecular dynamics simulations, we estimate the potentials of mean force between BMIM+ and Cl- ions in DMSO-water mixtures. Analysis reveals a significant increase in the thermodynamic stability of both contact ion pair (CIP) and solvent-assisted ion pair (SAIP) states with a rising DMSO mole fraction. Thermodynamic assessments highlight the entropic stabilization of CIP states and SAIP states in pure water, in DMSO-water mixtures, and in pure DMSO. The structural analysis reveals that in comparison to the DMSO local density, the local water density is relatively very high around ion pairs, more specifically in the solvation shell of a chloride ion. Preferential binding coefficients also consistently indicate exclusion of DMSO from the ion pair in DMSO-water mixtures. To enhance our understanding regarding the solvent molecules kinetics around the ion pairs, the survival probabilities of DMSO and water are computed. The calculations reveal that the water molecules prefer a prolonged stay in the solvation shell of Cl- ions.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Computational study and in vitro alpha-glucosidase inhibitory effects of medicinal plants from a Thai folk remedy

The number of patients with type 2 diabetes mellitus (T2DM) has increased worldwide. Although an instant cure was achieved with the standard treatment acabose, unsatisfactory symptoms associated with cardiovascular disease after acabose administration have been reported. Therefore, it is important to explore new treatments. A Thai folk recipe has long been used for T2DM treatment, and it effectively decreases blood glucose. However, the mechanism of this recipe has never been proven. Therefore, the potential anti-T2DM effect of this recipe, which is used in Thai hospitals, was determined to involve alpha-glucosidase (AG) inhibition with a half maximal inhibitory concentration (IC50). In vitro experiments showed that crude Cinnamomum verum extract (IC50 = 0.35 ± 0.12 mg/mL) offered excellent inhibitory activity, followed by extracts from Tinospora crispa (IC50 = 0.69 ± 0.39 mg/mL), Stephania suberosa (IC50 = 1.50 ± 0.17 mg/mL), Andrographis paniculate (IC50 = 1.78 ± 0.35 mg/mL), and Thunbergia laurifolia (IC50 = 4.66 ± 0.27 mg/mL). However, the potencies of these extracts were lower than that of acabose (IC50 = 0.55 ± 0.11 mg/mL). Therefore, this study investigated and developed a formulation of this recipe using computational docking. Among 61 compounds, 7 effectively inhibited AG, including chlorogenic acid (IC50 = 819.07 pM) through 5 hydrogen bonds (HBs) and 2 hydrophobic interactions (HIs); β-sitosterol (IC50 = 4.46 nM, 6 HIs); ergosterol peroxide (IC50 = 4.18 nM, 6 HIs); borapetoside D (IC50 = 508.63 pM, 7 HBs and 2 HIs); borapetoside A (IC50 = 1.09 nM, 2 HBs and 2 His), stephasubimine (IC50 = 285.37 pM, 6 HIs); and stephasubine (IC50 = 1.09 nM, 3 HBs and 4 HIs). These compounds bind with high affinity to different binding pockets, leading to additive effects. Moreover, the pharmacokinetics of six of these seven compounds (except ergosterol peroxide) showed poor absorption in the gastrointestinal tract, which would allow for competitive binding to AG in the small intestine. These results indicate that the development of these 6 compounds into oral antidiabetic agents is promising.

Hydration and dynamics of L-glutamate ion in aqueous solution

Aqueous solutions of sodium l-glutamate (NaGlu) in the concentration range 0 < c/M ≤ 1.90 at 25 °C were investigated by dielectric relaxation spectroscopy (DRS) and statistical mechanics (1D-RISM and 3D-RISM calculations) to study the hydration and dynamics of the l-glutamate (Glu-) anion. Although at c → 0 water molecules beyond the first hydration shell are dynamically affected, Glu- hydration is rather fragile and for c ⪆ 0.3 M apparently restricted to H2O molecules hydrogen bonding to the carboxylate groups. These hydrating dipoles are roughly parallel to the anion moment, leading to a significantly enhanced effective dipole moment of Glu-. However, l-glutamate dynamics is determined by the rotational diffusion of individual anions under hydrodynamic slip boundary conditions. Thus, the lifetime of the hydrate complexes, as well as of possibly formed [Na+Glu-]0 ionpairs and l-glutamate aggregates, cannot exceed the characteristic timescale for Glu- rotation.

Computational NMR Study of Ion Pairing of 1-Decyl-3-methyl-imidazolium Chloride in Molecular Solvents

The 1H NMR spectra of 10-5 mole fraction solutions of 1-decyl-3-methyl-imidazolium chloride ionic liquid in water, acetonitrile, and dichloromethane have been measured. The chemical shift of the proton at position 2 in the imidazolium ring of 1-decyl-3-methyl-imidazolium (H2) is rather different for all three samples, reflecting the shifting equilibrium between the contact pairs and free fully solvated ions. Classical molecular dynamics simulations of the 1-decyl-3-methyl-imidazolium chloride contact ion pair as well as of free ions in water, acetonitrile, and dichloromethane have been conducted, and the quantum mechanics/molecular mechanics methods have been applied to predict NMR chemical shifts for the H2 proton. The chemical shift of the H2 proton was found to be primarily modulated by hydrogen bonding with the chloride anion, while the influence of the solvents-though differing in polarity and capabilities for hydrogen bonding-is less important. By comparing experimental and computational results, we deduce that complete disruption of the ionic liquid into free ions takes place in an aqueous solution. Around 23% of contact ion pairs were found to persist in acetonitrile. Ion-pair breaking into free ions was predicted not to occur in dichloromethane.

Aqueous solutions of hydroxyl-functionalized ionic liquids: Molecular dynamics studies

A series of aqueous solutions of 1-(n-hydroxyalkyl)-3-(n-hydroxyalkyl) imidazolium bromide ([HOC_nC_mOHIm][Br], with n and m = 2, 6,10 and 14) were studied by atomistic molecular dynamics simulations. Structural properties were characterized by the radial distribution functions between different pairs, angular distributions and aggregation numbers. Dynamics of the system has been investigated by computing the diffusion of the ions and molecules. Structures of the aggregates formed depend upon the length of the hydroxyalkyl chains. The long-distance spatial correlations observed in solutions with cations having long chain substituent are arising due to the formation of intercalated structures. A thin film like structure is formed in solutions having longer hydroxyalkyl chains, with the structure stabilized by the dispersion interactions between the interdigitated alkyl chains and the hydrogen bond formation between the hydroxyl group of a cation with head group of a different cation. Anions are dispersed near the surface of the film.

Molecular Dynamics Insights into the Nanoscale Structural Organization and Local Interaction of Aqueous Solutions of Ionic Liquid 1-Butyl-3-methylimidazolium Nitrate

Considering the growing number of applications of the aqueous ionic liquids (ILs), atomistic molecular dynamics (MD) simulations were used to probe the effect of water molar fraction, x_w, ranging from 0.00 to 0.90, on the nanoscale local structure of 1-butyl-3-methylimidazolium nitrate, [bmim][NO₃], IL. The results prove that, with water addition, the cation-anion, cation-cation, and anion-anion structural correlations are weakened, while strong anion-water and unconventional cation-water hydrogen bonds are formed in the solutions. Water molecules were detected as bridges between nitrate anions, and the water cluster size distribution at different x_w's was investigated. Simulation shows a similar pattern of probability densities for water and anion around the acidic hydrogen atoms of the reference cation ring, while both species move away from the cation butyl chain. Increasing the water concentration to x_w = 0.90 causes decreasing of the local arrangement of the nearest-neighboring cations, because of the weakening of cation-cation π-π stacking. In addition, this dilution reduces the probability of the in-plane cation-anion conformation, disrupts both the polar ionic network and nonpolar domains, and diminishes the nanoaggregation of the cation butyl chains compared to those of the neat IL. These results can rationalize the origins of the fluidity enhancements and transport property trends upon adding water to the imidazolium-based ILs. The current study proposes a deep atomistic-level insight into the complex coupling between water concentration, microscopic structure, and local interactions of aqueous imidazolium-based ILs with hydrophilic anions.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Glucose in dry and moist ionic liquid: vibrational circular dichroism, IR, and possible mechanisms

Ionic liquids and their mixtures with water show remarkable features in cellulose processing. For this reason, understanding the behavior of carbohydrates in ionic liquids is important. In the present study, we investigated three d-glucose isomers (α, β and open-chain) in 1-ethyl-3-methylimidazolium acetate in the presence and absence of water, through ab initio molecular dynamics simulations. In the complex hydrogen bonding network of these mixtures, the most interesting observation is that upon water addition every hydrogen bond elongates, except the glucose-glucose hydrogen bond for the open-chain and the α-form which shortens, clearly showing the beginning of the crystallization process. The ring glucose rearranges from on-top to in-plane and the open form changes from a coiled to a more linear arrangement when adding water which explains the contradiction that the center of mass distances of the glucose molecules with other glucose molecules grow while the hydrogen bonds shorten. The appearance of coiled open forms indicates that the previously suggested isomerization between these forms is possible and might play a role in the solubility of the related carbohydrates. The calculated IR and VCD spectra reveal insight into the intermolecular interactions, with good to excellent agreements with experimental spectra. Investigating the role of the cation, distances between the acidic carbon atom of the cation and the glucose carbon atom where ring closure and opening occurs are found, which are way shorter than dispersion-like interactions between aliphatic hydrocarbons.

Decelerated Liquid Dynamics Induced by Component-Dependent Supercooling in Hydrogen and Deuterium Quantum Mixtures

Isotopic mixtures of p-H₂ and o-D₂ molecules have been an attractive binary system because they include two kinds of purely isotopic molecules which possess the same electronic potential but the twice different mass inducing differently pronounced nuclear quantum effects (NQEs). Accessing details of structures and dynamics in such quantum mixtures combining complex molecular dynamics with NQEs of different strengths remains a challenging problem. Taking advantage of the nonempirical molecular dynamics method which describes p-H₂ and o-D₂ molecules, we found that the liquid dynamics slows down at a specific mixing ratio, which can be connected to the observed anomalous slowdown of crystallization in the quantum mixtures. We attributed the decelerated dynamics to the component-dependent supercooling of p-H₂ taking place in the mixtures, demonstrating that there is an optimal mixing ratio to hinder crystallization. The obtained physical insights will help in experimentally controlling and achieving unknown quantum mixtures including superfluid.

Current Status of AMOEBA-IL: A Multipolar/Polarizable Force Field for Ionic Liquids

Computational simulations of ionic liquid solutions have become a useful tool to investigate various physical, chemical and catalytic properties of systems involving these solvents. Classical molecular dynamics and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of IL systems have provided significant insights at the atomic level. Here, we present a review of the development and application of the multipolar and polarizable force field AMOEBA for ionic liquid systems, termed AMOEBA–IL. The parametrization approach for AMOEBA–IL relies on the reproduction of total quantum mechanical (QM) intermolecular interaction energies and QM energy decomposition analysis. This approach has been used to develop parameters for imidazolium– and pyrrolidinium–based ILs coupled with various inorganic anions. AMOEBA–IL has been used to investigate and predict the properties of a variety of systems including neat ILs and IL mixtures, water exchange reactions on lanthanide ions in IL mixtures, IL–based liquid–liquid extraction, and effects of ILs on an aniline protection reaction.

Distance Angle Descriptors of the Interionic and Ion-Solvent Interactions in Imidazolium-Based Ionic Liquid Mixtures with Aprotic Solvents: A Molecular Dynamics Simulation Study

The aim of this paper is to quantify the changes of the interionic and ion-solvent interactions in mixtures of imidazolium-based ionic liquids, having tetrafluoroborate (BmimBF₄), hexafluorophosphate (BmimPF₆), trifluoromethylsulfonate (BmimTFO), or bis(trifluoromethanesulfonyl)imide (BmimTFSI), anions, and polar aprotic molecular solvents, such as acetonitrile (AN), γ-butyrolactone (GBL), and propylene carbonate (PC). For this purpose, we calculate, using the nearest-neighbor approach, the average distance between the imidazolium ring H atom in positions 2, 4, and 5 (H^2,4,5) and the nearest high-electronegativity atom of the solvent or anion (X) as distance descriptors, and the mean angle formed by the C^2,4,5-H^2,4,5 bond and the H^2,4,5···X axis around the H^2,4,5 atom as angular descriptors of the cation-anion and cation-solvent interactions around the ring C-H groups. The behavior of these descriptors as a function of the ionic liquid mole fraction is analyzed in detail. The obtained results show that the extent of the change of these descriptors with respect to their values in the neat ionic liquid depends both on the nature of the anion and on the mixture composition. Thus, in the case of the mixtures of the molecular solvents with BmimBF₄ and BmimTFO, a small change of the distance and a drastic increase of the angular descriptor corresponding to the cation-anion interactions are observed with decreasing mole fraction of the ionic liquid, indicating that the anion moves from the above/below position (with respect to the imidazolium ring plane) to a position that is nearly linearly aligned with the C²-H² bond and hinders the possible interaction between the C²-H² group and the solvent molecules. On the other hand, in the case of mixtures of BmimTFSI and BmimPF₆ with the molecular solvents, both the observed increase of the distance descriptor and the slight change of the angular descriptor with decreasing ionic liquid mole fraction are compatible with the direct interactions of the solvent with the C²-H² group. The behavior of these descriptors is correlated with the experimentally observed ¹H chemical shift of the C²-H² group and the red shift of the C²-H² vibrational mode, particularly at low ionic liquid mole fractions. The present results are thus of great help in interpreting these experimental observations.

The peculiar effect of water on ionic liquids and deep eutectic solvents

Ionic liquids (ILs) and deep eutectic solvents (DESs) have been suggested as eco-friendly alternatives to organic solvents. A trace amount of water is often unavoidable as impurity, and water is also added on purpose to reduce their problematically high viscosity and lower their high price. Understanding the distinct effects of water on the properties of ILs/DESs is highly important. In this review, we collect published experimental and theoretical results for IL/DES-H2O systems at varied water concentrations and analyze them. Results from mechanistic studies, thermodynamic modelling and advanced experiments are collected and critically discussed. Six commonly studied IL/DES-H2O systems were selected to map experimental observations onto microscopic results obtained in mechanistic studies. A great variety of distinct contours of the excess properties can be observed over the entire compositional range, indicating that the properties of IL/DES-H2O systems are highly unpredictable. Mechanistic studies clearly demonstrate that the added H2O rapidly changes the heterogeneous 3D structures of pure ILs/DESs, leading to very different properties and behaviour. There are similarities between aqueous electrolytes and IL/DES solutions but the bulky and asymmetric organic cations in ILs/DESs do not conform to the standard salt dissolution and hydration concepts. Thermodynamic modelling previously assumes ILs/DESs to be either a neutral ion-pair or completely dissociated ions, neglecting specific ion hydration effects. A new conceptual framework is suggested for thermodynamic modelling of IL/DES-H2O binary systems to enable new technologies for their practical applications.

Hydrophilic Ionic Liquid Mixtures of Weakly and Strongly Coordinating Anions with and without Water

With the aid of ab initio molecular dynamics simulations, we investigate an ionic liquid (IL) mixture composed of three components 1-butyl-3-methylimidazolium [C₄C₁Im]⁺, tetrafluoroborate [BF₄]^-, and chloride [Cl]^- without and with water. In the pure IL mixture, we observe an already complex network of interactions between cations and anions, and addition of water to the system even extends the complexity. Observed number integrals show that the coordination number between cations and anions is reduced in the system with water compared to that in the pure system. Further studies show that the Coulombic network of the strongly coordinating anion [Cl]^- is disturbed by water, while that of the weakly coordinating anion [BF₄]^- is not. These observations can also be confirmed by the Voronoi polyhedra analysis, which shows that the polar network of microheterogeneous IL collapses by the introduction of water. Hydrogen-acceptor interactions revealed that the [Cl]^- anions are transferred from being situated in the IL to the water continuum, while [BF₄]^- is almost unperturbed; these effects mainly influence the interplay of the ionic liquid network.

The local structure in the BmimPF6/acetonitrile mixture: the charge distribution effect

The effect of the charge distribution on the local structure in the binary mixture of 1-butyl-3-methylimidazolium hexafluorophosphate (BmimPF₆) ionic liquid and acetonitrile is investigated over the entire composition range.

Hydration and Counterion Binding of [C12MIM] Micelles

Surface-active ionic liquids based on imidazolium cations are promising targets for micellar catalysis in aqueous solution, yielding enhanced rate constants compared to surfactants based on n-alkyltrimethylammonium cations and exhibiting a pronounced counterion dependence ( Bica Chem. Commun. 2012 , 48 , 5013 - 5015 ; Cognigni Phys. Chem. Chem. Phys. 2016 , 18 , 13375 - 13384 ). Probably most relevant to these effects is the interplay between headgroup hydration and counterion binding. To obtain more detailed information on these effects, aqueous solutions of 1-dodecyl-3-methylimidazolium ([C₁₂MIM]) bromide, iodide, and triflate (TfO^-) were investigated at 45 °C using broadband dielectric spectroscopy, viscosity measurements, and small-angle X-ray scattering experiments. Effective hydration numbers were determined, and information on the locations and mobilities of the condensed counterions, X^-, was derived. It was found that [C₁₂MIM] halide micelles were less hydrated than the corresponding n-dodecyltrimethylammonium ([C₁₂TA]X) aggregates. Together with their somewhat weaker counterion condensation, this difference probably explains their higher catalytic activity. Whereas [C₁₂MIM]Br micelles remained roughly spherical in the studied concentration range, rodlike aggregates were formed at high concentrations of the iodide and, in particular, the triflate surfactants. It appears that the much lower mobility of condensed TfO^- counterions is the reason for the very low catalytic activity of [C₁₂MIM]TfO micelles.

Aqueous ionic liquids and their effects on protein structures: an overview on recent theoretical and experimental results

Ionic liquids (ILs) are used in a variety of technological and biological applications. Recent experimental and simulation results reveal the influence of aqueous ionic liquids on the stability of protein and enzyme structures. Depending on different parameters like the concentration and the ion composition, one can observe distinct stabilization or denaturation mechanisms for various ILs. In this review, we summarize the main findings and discuss the implications with regard to molecular theories of solutions and specific ion effects. A preferential binding model is introduced in order to discuss protein-IL effects from a statistical mechanics perspective. The value of the preferential binding coefficient determines the strength of the ion influence and indicates a shift of the chemical equilibrium either to the native or the denatured state of the protein. We highlight the role of water in order to explain the self-association behavior of the IL species and discuss recent experimental and simulation results in the light of the observed binding effects.

Proteins in Ionic Liquids: Current Status of Experiments and Simulations

In the last two decades, while searching for interesting applications of ionic liquids as potent solvents, their solvation properties and their general impact on biomolecules, and in particular on proteins, gained interest. It turned out that ionic liquids are excellent solvents for protein refolding and crystallization. Biomolecules showed increased solubilities and stabilities, both operational and thermal, in ionic liquids, which also seem to prevent self-aggregation during solubilization. Biomolecules can be immobilized, e.g. in highly viscous ionic liquids, for particular biochemical processes and can be designed to some extent by the proper choice of the ionic liquid cations and anions, which can be characterized by the Hofmeister series.

Bio ionic liquids and water mixtures: a structural study

In this study we have explored, by means of ab initio molecular dynamics, a subset of three different water/cho⁺–phe⁻ mixtures.

The properties of residual water molecules in ionic liquids: a comparison between direct and inverse Kirkwood–Buff approaches

Atomistic molecular dynamics simulations of aqueous ionic liquid mixtures were performed in order to compare the resulting Kirkwood–Buff integrals with experimental data and the corresponding integrals derived by an inverse Kirkwood–Buff approach.

ILPC: simple chemometric tool supporting the design of ionic liquids

Background

Ionic liquids (ILs) found a variety of applications in today’s chemistry. Since their properties depend on the ions constituting particular ionic liquid, it is possible to synthetize IL with desired specification, dependently on its further function. However, this task is not trivial, since knowledge regarding the influence of particular ion on the property of concern is crucial. Therefore, there is a strong need for new, fast and inexpensive methods supporting the process of ionic liquids’ design, making it possible to predefine IL’s properties even before the synthesis.

Results

We have developed a simple tool (called Ionic Liquid PhysicoChemical predictor: ILPC) that allows for the simultaneous qualitative prediction of four physicochemical properties of ionic liquids: viscosity, n-octanol–water partition coefficient, solubility and enthalpy of fusion. By the means of Principal Component Analysis, we studied 172 ILs and defined distribution trends of those four properties, dependently on the ILs structures. We proved that the qualitative prediction of mentioned properties could be performed on the basis of most simple information we can deliver about ILs, which are their molecular formulas.

Conclusions

Created tool presented in this paper allows fast, pre-synthesis screening of ILs, with the omission of any experimental steps. It can be helpful in the process of designing ILs with preferred properties. We proved that the information encrypted in molecular formula of ionic liquid could be a valuable source of knowledge regarding the IL’s viscosity, n-octanol–water partition coefficient, solubility and enthalpy of fusion. Moreover, we proved that the influence of both ions, constituting the IL, on each of those four properties indicates same, additive trend.Graphical Abstract

Water Nanocluster Formation in the Ionic Liquid 1-Butyl-3-methylimidazolium Tetrafluoroborate ([C4mim][BF4])-D2O Mixtures

The microstructure of mixtures of deuterated water in the ionic liquid [C4mim][BF4] is investigated by small-angle neutron scattering (SANS) measurement. In the salt-rich region, water dissolves in the ionic liquid up to 0.7 mole fraction, whereupon distinct, nanometer-sized water clusters are observed. These water nanoclusters increase in size with increasing water addition while the mole ratio of water dissolved into the ionic liquid nanostructure increases from 2 to 4. These results provide direct confirmation for recent simulations as well insight into the source of nonidealities in some thermophysical and transport properties (e.g., density and viscosity) of salt-rich aqueous mixtures reported in the literature.

Filling environmental data gaps with QSPR for ionic liquids: Modeling n-octanol/water coefficient

Ionic liquids (ILs) form a wide group of compounds characterized by specific properties that allow using ILs in different fields of science and industry. Regarding that the growing production and use of ionic liquids increase probability of their emission to the environment, it is important to estimate the ability of these compounds to spread in the environment. One of the most important parameters that allow evaluating environmental mobility of compound is n-octanol/water partition coefficient (KOW). Experimental measuring of the KOW values for a large number of compounds could be time consuming and costly. Instead, computational predictions are nowadays being used more often. The paper presents new Quantitative Structure-Property Relationship (QSPR) model that allows predicting the logarithmic values of KOW for 335 ILs, for which the experimentally measured values had been unavailable. We also estimated bioaccumulation potential and point out which group of ILs could have negative impact on environment.

Confined water in imidazolium based ionic liquids: a supramolecular guest@host complex case

Traces of water in some ionic liquids can be regarded as a guest@host supramolecular structure even when diluted in solvents with high dielectric constants.

Mobility and association of ions in aqueous solutions: the case of imidazolium based ionic liquids

Combining MD simulations and RISM calculations with experiments, we demonstrated that the interionic interactions of investigated ILs as model 1,1 electrolytes in water solutions are weak but evidently dependent on the molecular structure.

Dilute or Concentrated Electrolyte Solutions? Insight from Ionic Liquid/Water Electrolytes

When room-temperature ionic liquids (IL) are used as an electrolyte, their transport behaviors are still under heavy debate due to their complicated ion-associations. In this article, we conducted molecular dynamics simulations to study the molecular scale ion associations from the very dilute 1-butyl-3-methylimidazolium iodide/water solution to the pure IL. It revealed that ions are localized in a multicoordinated ion cage structure with nanoseconds in concentrated IL solutions. Dynamics analyses indicate that the transport of this solution can be depicted by the Debye-Hückel model only in dilute IL/water electrolyte. The velocity and rotational correlation functions showed that the lifetime of translational and rotational motions are at the level of picoseconds and nanoseconds, respectively, because of the ion cage effect. The lifetime of ion association demonstrated that the recombination of association ions was prevalent in IL solutions. It means that the dipolar or stable contact ion-pairs model may not be suitable for depicting the IL transport.

Aqueous Brønsted-Lowry Chemistry of Ionic Liquid Ions

Ionic liquids have become commonplace materials found in research laboratories the world over, and are increasingly utilised in studies featuring water as co-solvent. It is reported herein that proton activities, aH (+) , originating from auto-protolysis of H2O molecules, are significantly altered in mixtures with common ionic liquids comprised of Cl(-), [HSO4 ](-), [CH3SO4 ](-), [CH3COO](-), [BF4](-), relative to pure water. paH (+) values, recorded in partially aqueous media as -log(aH (+)), are observed over a wide range (∼0-13) as a result of hydrolysis (or acid dissociation) of liquid salt ions to their associated parent molecules (or conjugate bases). Brønsted-Lowry acid-base character of ionic liquid ions observed is rooted in equilibria known to govern the highly developed aqueous chemistry of classical organic and inorganic salts, as their well-known aqueous pKs dictate. Classical salt behaviour observed for both protic and aprotic ions in the presence of water suggests appropriate attention need be given to relevant chemical systems in order to exploit, or avoid, the nature of the medium formed.

Dissolution of cellulose in ionic liquids: an ab initio molecular dynamics simulation study

Interactions determining the dissolution of a monomer of β-cellulose, i.e., cellobiose in a room temperature ionic liquid, [Emim][OAc], have been studied using ab initio molecular dynamics simulations. Although anions are the predominant species in the first coordination shell of cellobiose, cations too are present to a minor extent around it. The presence of low concentration of water in the solution does not significantly alter the nature of the coordination environment of cellobiose. All intra-molecular hydrogen bonds of anti-syn cellobiose are replaced by inter-molecular hydrogen bonds formed with the anions, whereas the anti-anti conformer retains an intramolecular hydrogen bond.

On the collective network of ionic liquid/water mixtures. IV. Kinetic and rotational depolarization

Dielectric spectroscopy is a measure of the collective Coulomb interaction in liquid systems. Adding ionic liquids to an aqueous solution results in a decrease of the static value of the generalized dielectric constant which cannot be attributed to kinetic depolarization models characterized by the static conductivity and rotational relaxation constant. However, a dipolar Poisson-Boltzmann model computing the water depolarization in the proximity of ions is not only successful for simple electrolytes but also in case of molecular ionic liquids. Moreover, our simple geometric hydration model is also capable to explain the dielectric depolarization. Both models compute the dielectric constant of water and obtain the overall dielectric constant by averaging the values of its components, water and the ionic liquid, weighted by their volume occupancies. In this sense, aqueous ionic liquid mixtures seem to behave like polar mixtures.

A comparative study of the effects of the Hofmeister series anions of the ionic salts and ionic liquids on the stability of α-chymotrypsin

Direct interactions between the anion and the catalytic amino acid residues lead to denaturation of CT.

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.

Multiscale modeling of the trihexyltetradecylphosphonium chloride ionic liquid

Hierarchical trihexyltetradecylphosphonium cationic and chloride anionic models.