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

Capturing the Effect of Anion Type on the Intermolecular Interactions between Water and Imidazolium-Based Ionic Liquids: A Comparative DFT Study

The studies on ionic liquids (ILs) and their interaction with different solvents have always been an interesting topic for experimental and computational chemists. Recently, however, deep insights on the molecular structures of the IL-water binary mixtures have been mainly performed through classical simulations. Here, a comprehensive quantum mechanical study is presented on seven 1-butyl-3-methylimidazolium-based ILs in the absence and presence of water. As the most important intermolecular interaction between ionic moieties of ILs and water molecules, hydrogen bonding is studied through different bonding analyses. The effect of different anions, [NO3]-, [HSO4]-, [SCN]-, [DCA]-, [BF4]-, [PF6]-, and [NTf2]-, on the behavior of ILs interacting with a sample of water molecules is investigated. Comparing the implicit and explicit approaches to consider water solvent indicated that the structure of ILs in the solvent depends on the selected solvent model. By considering explicit water molecules, we analyzed the intermolecular interactions between ILs and the water sample. The energy decomposition analysis indicated that the stability of the IL···water systems is mainly due to the electrostatic component of the total interaction energy. The interaction region indicator (IRI) analysis discovered that chemical bond and van der Waals (vdW) interactions are important in the IL···water systems. Indeed, investigation of each ion/ion pair surrounded by ten nearest neighbor water molecules discovered that the vdW interactions are responsible for the cation···anion and the cation···water interactions, while chemical bonding is important in the anion···water and the water···water interactions. Therefore, the anion···water interaction requires further analysis. The quantum theory of atoms in molecules verified the ionic nature of the H-bond in the anion···water interaction. The IRI analysis showed that the interaction between water molecules and cyano-based anions, [SCN]- and [DCA]-, is only due to chemical bonding, while in the oxygenated anions, [NO3]- and [HSO4]-, the vdW forces are also important. For the other anions, [BF4]-, [PF6]-, and [NTf2]-, the vdW forces have the main contribution in the anion···water interaction. Natural bond orbital analysis indicated that these intermolecular interactions originate from nanion → σO-H* electron transfer. Finally, the law of matching water affinity (LMWA) using energy-based parameters was used to predict the hydrophilicity of ILs as follows: [BMIM][NO3] > [BMIM][SCN] > [BMIM][DCA] > [BMIM][HSO4] > [BMIM][BF4] > [BMIM][NTf2] > [BMIM][PF6]. Results obtained in the current work give insights into the electronic nature of intermolecular interactions between ILs and water molecules, which is necessary due to importance of water in modifying properties of ILs in various applications.

Single Particle Dynamics at the Free Surface of Imidazolium-Based Ionic Liquids

In this work, we carry out a systematic computer simulation investigation of the single particle dynamics at the free surface of imidazolium-based room temperature ionic liquids by applying intrinsic surface analysis. Besides assessing the effect of the potential model and temperature, we focus in particular on the effect of changing the anion type, and, hence, their shape and size. Further, we also address the role of the length of the cation alkyl chains, known to protrude into the vapor phase, on the surface dynamics of the ions. We observe that the surface dynamics of ionic liquids, being dominated by strong electrostatic interactions, is about 2 orders of magnitude slower than that for common molecular liquids. Furthermore, the free energy driving force for exposing apolar chains to the vapor phase "pins" the cations at the surface layer for much longer than anions, allowing them to perform noticeable lateral diffusion at the liquid surface during their stay there. On the other hand, anions, accumulated in the second layer beneath the liquid surface, stay considerably longer here than in the surface layer. The ratio of the mean surface residence time of the cations and anions depends on the relative size of the two ions: larger size asymmetry typically corresponds to larger values of this ratio. We also find, in a clear contrast with the bulk liquid phase behavior, that anions typically diffuse faster at the liquid surface than cations. Finally, our results show that the surface dynamics of the ions is largely determined by the apolar layer of the cation alkyl chains at the liquid surface, as in the absence of such a layer, cations and anions are found to behave similarly with respect to their single particle dynamics.

Spatial organization of the ions at the free surface of imidazolium-based ionic liquids

HYPOTHESIS

Experimental information on the molecular scale structure of ionic liquid interfaces is controversial, giving rise to two competing scenarios, namely the double layer-like and "chessboard"-like structures. This issue can be resolved by computer simulation methods, at least for the underlying molecular model. Systematically changing the anion type can elucidate the relative roles of electrostatic interactions, hydrophobic (or, strictly speaking, apolar) effects and steric restrictions on the interfacial properties.

SIMULATIONS

Molecular dynamics simulation is combined with intrinsic analysis methods both at the molecular and atomic levels, supplemented by Voronoi analysis of self-association.

FINDINGS

We see no evidence for the existence of a double-layer-type arrangement of the ions, or for their self-association at the surface of the liquid. Instead, our results show that cation chains associate into apolar domains that protrude into the vapour phase, while charged groups form domains that are embedded in this apolar environment at the surface. However, the apolar chains largely obscure the cation groups, to which they are bound, while the smaller and more mobile anions can more easily access the free surface, leading to a somewhat counterintuitive net excess of negative charge at the interface. Importantly, this excess charge could only be identified by applying intrinsic analysis.

Molecular dynamics simulation study of sodium ion structure & dynamics in water in ionic liquids electrolytes using 1-butyl-3-methyl imidazolium tetrafluoroborate and 1-butyl-3-methyl imidazolium hexafluorophosphate

In this paper, we have performed an all-atom molecular dynamics simulation to understand the structure and dynamics of Na+ ions in water mixed Ionic liquids (Water in Ionic liquid). Two ionic liquid (IL) systems consist of (1) 1-butyl-3-methylimidazolium [BMIM] tetrafluoroborate [BF4] and (2) 1-butyl-3-methylimidazolium [BMIM] hexafluorophosphate [PF6] were considered in this work. We understand various inter-molecular structures and dynamic and thermodynamic behaviours of Na+ ions in the water-mixed IL systems. The water (H2O) mole fractions (x) varied from 0.33 to 0.71. The neat ILs [BMIM][BF4] and [BMIM][PF6] pairwise radial distribution functions show a decrease with an increase in x. The [BMIM][PF6] exhibits a strong coordination structure with Na+ ions across the entire range of x values. The rdf between the pairs of Na+-[PF6] presents a significant interaction compared to Na+ and [BF4]. The Na + ions manifested greater coordination with H2O In H2O-[BMIM][PF6] compared to H2O-[BMIM][BF4]. The self-diffusion coefficient (D) values of Na + ions increase with the rise in x in both ILs. The D values of Na + ions are 10-fold higher in [BMIM][BF4] than [BMIM][PF6]. The ionic conductivity values are higher for [BMIM][BF4]. Overall, this paper unveils molecular-level insights for understanding the behavior of Na+ ions in the water in ionic liquid systems.

Ionic liquid mixtures as energy storage materials: a preliminary and comparative study based on thermal storage density

Fifteen equimolar binary mixtures are synthesized and thermophysically evaluated in this study. These mixtures are derived from six ionic liquids (ILs) based on methylimidazolium and 2,3-dimethylimidazolium cations with butyl chains. The objective is to compare and elucidate the impact of small structural changes on the thermal properties. The preliminary results are compared to previously obtained results with mixtures containing longer eight-carbon chains. The study demonstrates that certain mixtures exhibit an increase in their heat capacity. Additionally, due to their higher densities, these mixtures achieve a thermal storage density equivalent to that of mixtures with longer chains. Moreover, their thermal storage density surpasses that of some conventional materials commonly used for energy storage.

Anomalous Dependence of Translational Diffusion on the Water Mole Fraction for Solute Molecules Dissolved in a 1-Butyl-3-methylimidazolium Tetrafluoroborate/Water Mixture

Translational diffusion coefficients of carbon monoxide (CO), diphenylacetylene (DPA), and diphenylcyclopropenone (DPCP) were determined in mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate ([C₄mim]BF₄) and water using transient grating spectroscopy at different mole fractions of water (x_w). While DPA exhibited a larger diffusion coefficient than DPCP at low water mole fractions (x_w < 0.7), as observed for conventional liquids and ionic liquids (ILs), it was smaller at high mole fractions (x_w > 0.9). The apparent molecular radius of DPA determined using the Stokes-Einstein equation at x_w > 0.9 is close to the radius of an IL cluster in a water pool as determined from small-angle neutron scattering experiments (J. Bowers et al., Langmuir, 2004, 20, 2192-2198), suggesting that the DPA molecules are trapped in IL clusters in the water pool and move together. The solvation state of DPCP in the mixture was studied using Raman spectroscopy. Dramatically strong water/DPCP hydrogen bonding was observed at higher water mole fractions, suggesting that DPCP is located near the cluster interfaces. The large diffusion coefficient of DPCP suggests that hopping of DPCP between IL clusters occurs through hydrogen bonding with water.

Quantitative Relation between Ionic Diffusivity and Ionic Association in Ionic Liquid-Water Mixtures

Molecular dynamic simulations of aqueous mixtures of imidazolium ionic liquids (ILs) were performed to elucidate the dependence of the ionic diffusivity on the microscopic structures changed by water. Two distinct regimes of the average ionic diffusivity (D_ave) were identified with the increased water concentrations: the jam regime with slowly increased D_ave and the exponential regime with rapidly increased D_ave, which are found to be directly correlated to the ionic association. Further analysis leads to two general relationships independent of IL species between D_ave and the degree of ionic association: (i) a consistent linear relationship between D_ave and the inverse of ion-pair lifetimes (1/τ_IP) in the two regimes and (ii) an exponential relationship between normalized diffusivities (D̃_ave) and short-ranged interactions between cations and anions (Ẽ_ions), with different interdependent strengths in the two regimes. These findings revealed and quantified the direct correlation between dynamic properties and ionic association in IL-water mixtures.

Water-Controlled Structural Transition and Charge Transfer of Interfacial Ionic Liquids

Clarification of the water-induced structural transitions and electron transfer between ionic liquids (ILs) and a solid surface allows for establishing a unified view of the electrical properties of interfacial ILs via a hitherto unexplored pathway. Here, we propose a simple and effective method to quantitatively identify and extract the transferred electrons between ILs and a solid surface, while demonstrating the critical structural transition of interfacial ILs from ordered stripe structures to disordered aggregation structures. The formation of hydrated anions, rooted in the hydrogen bonds of O-H···O between the anion and water, lies at the tipping point where electron transfer ends and aggregation structure begins. In addition, it is discovered to what extent the hydrophilicity of substrates can affect electron transfer, and a regulation method based on the electric field is explored. These experimental findings may refresh our knowledge of interfacial ILs and provide an effective method for evaluating the intrinsic electrical features of the ILs-solid surface.

Interfacial Water and Microheterogeneity in Aqueous Solutions of Ionic Liquids

In this work, aqueous solutions of two prototypical ionic liquids (ILs), [BMIM][BF4] and [BMIM][TfO], were investigated by UV Raman spectroscopy and small-angle neutron scattering (SANS) in the water-rich domain, where strong heterogeneities at mesoscopic length scales (microheterogeneity) were expected. Analyzing Raman data by a differential method, the solute-correlated (SC) spectrum was extracted from the OH stretching profiles, emphasizing specific hydration features of the anions. SC-UV Raman spectra pointed out the molecular structuring of the interfacial water in these microheterogeneous IL/water mixtures, in which IL aggregates coexist with bulk water domains. The organization of the interfacial water differs for the [BMIM][BF4] and [BMIM][TfO] solutions, being affected by specific anion-water interactions. In particular, in the case of [BMIM][BF4], which forms weaker H-bonds with water, the aggregation properties clearly depend on concentration, as reflected by local changes in the interfacial water. On the other hand, stronger water-anion hydrogen bonds and more persistent hydration layers were observed for [BMIM][TfO], which likely prevent changes in IL aggregates. The modeling of SANS profiles, extended to [BPy][BF4] and [BPy][TfO], evidences the occurrence of significant concentration fluctuations for all of the systems: this appears as a rather general phenomenon that can be ascribed to the presence of IL aggregation, mainly induced by (cation-driven) hydrophobic interactions. Nevertheless, larger concentration fluctuations were observed for [BMIM][BF4], suggesting that anion-water interactions are relevant in modulating the microheterogeneity of the mixture.

Structure of Sulfuric Acid Solutions Using Pair Distribution Function Analysis

Solvation and mesoscale ordering of sulfuric acid and other strong acid solutions leads to suppressed freezing points and strong rheological changes with varying concentration. While the solid-state structures are well-understood, studies focused on the evolving solvation structure in the solution phase have probed a limited concentration range (∼1-6 M). This study applies a total scattering approach in both the wide-angle X-ray scattering (WAXS) and pair distribution function (PDF) regimes to elucidate the evolving solvation structure over its full range of acid concentration (0-18 M). The emergence of a prepeak in the WAXS regime at intermediate concentrations indicates a transition from noninteracting sulfate molecules in the dilute limit to sterically limited sulfates at concentrations near its deep eutectic point. Fits to the PDF data quantify this trend, showing a transition from octahedrally hydrated sulfates up to 6-7 M concentrations, followed by gradual dehydration, and eventually reaching a solution structure similar to that of water-in-salt electrolyte systems at high acid concentrations.

Computer Simulation of the Surface of Aqueous Ionic and Surfactant Solutions

The surface of aqueous solutions of simple salts was not the main focus of scientific attention for a long while. Considerable interest in studying such systems has only emerged in the past two decades, following the pioneering finding that large halide ions, such as I^-, exhibit considerable surface affinity. Since then, a number of issues have been clarified; however, there are still several unresolved points (e.g., the effect of various salts on lateral water diffusion at the surface) in this respect. Computer simulation studies of the field have largely benefited from the appearance of intrinsic surface analysis methods, by which the particles staying right at the boundary of the two phases can be unambiguously identified. Considering complex ions instead of simple ones opens a number of interesting questions, both from the theoretical point of view and from that of the applications. Besides reviewing the state-of-the-art of intrinsic surface analysis methods as well as the most important advances and open questions concerning the surface of simple ionic solutions, we focus on two such systems in this Perspective, namely, the surface of aqueous mixtures of room temperature ionic liquids and that of ionic surfactants. In the case of the former systems, for which computer simulation studies have still scarcely been reported, we summarize the theoretical advances that could trigger such investigations, which might well be of importance also from the point of view of industrial applications. Computer simulation methods are, on the other hand, widely used in studies of the surface of surfactant solutions. Here we review the most important theoretical advances and issues to be addressed and discuss two areas of applications, namely, the inclusion of information gathered from such simulations in large scale atmospheric models and the better understanding of the airborne transmission of viruses, such as SARS-CoV-2.

Special Mixing Behavior of Chelate-based Ionic Liquid with Methanol

Compared to the general ionic liquids (ILs), a significant deviation of the binary mixtures of 1-decyl-3-methylimidazolium tri(hexafluoroacetylaceto)-copper(II) ([C₁₀ mim][Cu(hfacac)₃ ]) with methanol was found, indicating the way methanol interacts with ILs might be governed by the special structure of the chelating anion. IR results showed that the v (C2-H) of 1-decyl-3-methylimidazolium hexafluoroacetylacetonate ([C₁₀ mim][hfacac]) blue-shifted more significantly than that of [C₁₀ mim][Cu(hfacac)₃ ], meanwhile the v (C=O) red-shifted in [C₁₀ mim][Cu(hfacac)₃ ], which is contrast with that in [C₁₀ mim][hfacac]. Two-dimensional correlation analysis of the FTIR spectra indicated that the chelating cavity has little effect on the sequence of the ILs sites that interact with methanol. Combined with small angle X-ray scattering (SAXS) results, the picture of mixing processes in these two systems were proposed. Methanol interacts directly with the anion followed by the cation in [C₁₀ mim][hfacac], while methanol preferentially enters the chelating cavity and enhances the packing effect in the [C₁₀ mim][Cu(hfacac)₃ ] system.