New force field for molecular simulation of guanidinium-based ionic liquids

An all-atom force field was proposed for a new class of room temperature ionic liquids (RTILs), N,N,N',N'-tetramethylguanidinium (TMG) RTILs. The model is based on the AMBER force field with modifications on several parameters. The refinements include (1) fitting the vibration frequencies for obtaining force coefficients of bonds and angles against the data obtained by ab initio calculations and/or by experiments and (2) fitting the torsion energy profiles of dihedral angles for obtaining torsion parameters against the data obtained by ab initio calculations. To validate the force field, molecular dynamics (MD) simulations at different temperatures were performed for five kinds of RTILs, where TMG acts as a cation and formate, lactate, perchlorate, trifluoroacetate, and trifluoromethylsulfonate act as anions. The predicted densities were in good agreement with the experimental data. Radial distribution functions (RDFs) and spatial distribution functions (SDFs) were investigated to depict the microscopic structures of the RTILs.     

A model for self-diffusion of guanidinium-based ionic liquids: a molecular simulation study

We propose a novel self-diffusion model for ionic liquids on an atomic level of detail. The model is derived from molecular dynamics simulations of guanidinium-based ionic liquids (GILs) as a model case. The simulations are based on an empirical molecular mechanical force field, which has been developed in our preceding work, and it relies on the charge distribution in the actual liquid. The simulated GILs consist of acyclic and cyclic cations that were paired with nitrate and perchlorate anions. Self-diffusion coefficients are calculated at different temperatures from which diffusive activation energies between 32-40 kJ/mol are derived. Vaporization enthalpies between 174-212 kJ/mol are calculated, and their strong connection with diffusive activation energies is demonstrated. An observed formation of cavities in GILs of up to 6.5% of the total volume does not facilitate self-diffusion. Instead, the diffusion of ions is found to be determined primarily by interactions with their immediate environment via electrostatic attraction between cation hydrogen and anion oxygen atoms. The calculated average time between single diffusive transitions varies between 58-107 ps and determines the speed of diffusion, in contrast to diffusive displacement distances, which were found to be similar in all simulated GILs. All simulations indicate that ions diffuse by using a brachiation type of movement: a diffusive transition is initiated by cleaving close contacts to a coordinated counterion, after which the ion diffuses only about 2 A until new close contacts are formed with another counterion in its vicinity. The proposed diffusion model links all calculated energetic and dynamic properties of GILs consistently and explains their molecular origin. The validity of the model is confirmed by providing an explanation for the variation of measured ratios of self-diffusion coefficients of cations and paired anions over a wide range of values, encompassing various ionic liquid classes as well as the simulated GILs. The proposed diffusion model facilitates the qualitative a priori prediction of the impact of ion modifications on the diffusive characteristics of new ionic liquids.     

All-atom and united-atom simulations of guanidinium-based ionic liquids

Ionic liquids (ILs) have been widely used in separation, catalysis, electrochemistry, etc., and one of the most outstanding characteristics is that ILs can be tailored and tuned for specific tasks. In order to design and make better use of ionic liquids, the structures and properties relationship is indispensable. Both molecular dynamics and Monte Carlo simulations have been proved useful to understand the behavior of molecules at the microscale and the properties of the system. However, the quality of such simulations depends on force field parameters describing the interactions between atoms. All-atom (AA) or the united-atom (UA) force fields will be chosen because of the demand for more exact results or the lower computational cost, respectively. In order to make a systematic comparison of the two force fields, molecular simulations for four kinds of acyclic guanidinium-based ionic liquids (cations: (R2N)2C=N+<, anion: nitric or perchloric acid) were performed based on the AA and the UA force fields in this work. AA force field parameters were derived from our previous work (Fluid Phase Equilib., 2008, 272: 1-7), and the UA parameters were proposed in this work. Molecular dynamics simulation results for the AA and UA force fields were compared. Simulation densities are very similar to each other. Center of mass radial distribution functions (RDFs), site to site RDFs and spatial distribution functions (SDFs) were also investigated to depict the microscopic structures of the ILs.     

Novel guanidinium-based ionic liquids as stationary phases for capillary gas chromatography

<正>The present study describes guanidinium-based ionic liquids(GBILs) as stationary phases for capillary gas chromatography (CGC) and to the best of our knowledge,no related reports are available up to now.In this study,a hexaalkylguanidinium ionic liquid(DOTMG-NTf_2) was synthesized and coated statically onto capillary columns.Selectivity of the stationary phase was evaluated by separating Grob test mixture,test mixture,alcohols mixture,and fatty acid methyl esters mixture,and thermal stability was investigated as well.The present study demonstrates that GBILs as CGC stationary phases exhibit satisfactory selectivity and thermal stability and have a great potential as new candidates for CGC stationary phases.     

Molecular dynamics simulation of dual amino-functionalized imidazolium-based ionic liquids

In this work, the all-atom (AA) force fields were set up for three kinds of dual amino-functionalized imidazolium-based ionic liquids (ILs), composed by cations with different alkyl chain length and amino acid anion [Gly]⁻. The force field was based on our previous work and the default parameters were developed in this study. Molecular dynamics simulations were performed. Validation was carried out by comparing simulation densities with experimental data, and good agreement was obtained. Molar volume and heat capacity at constant pressure were predicted. Mean square displacements for these ILs were computed and these ILs were proved to move very slowly. It may be caused by hydrogen-bonded network between ions and the terminal azyl. To depict the microscopic structures of the ILs, many types of radial distribution functions were investigated. It is interesting to find that not only the cation and anion, but also the anions themselves will form hydrogen bonds.     

A force field for guanidinium-based ionic liquids that utilizes the electron charge distribution of the actual liquid: a molecular simulation study

We propose a new all-atom force field for guanidinium-based ionic liquids (GILs) which is based on the charge distribution in the actual liquid. It comprises all cations that can be built by attaching alkyl chains of variable length to an acyclic or cyclic guanidinium compound and that are paired with nitrate or perchlorate anions. We based the parametrization of the force field on liquid-phase charge distributions to improve the prediction of energetic and dynamic properties of GILs. The impact of electron charge transfer and polarization on various properties of GILs is systematically assessed. A significant average electron charge transfer between -0.12 and -0.06 e from anions to the central guanidinium group of the cations and a strong polarization of acyclic cations are observed by applying a combined quantum mechanical/molecular mechanical (QM/MM) approach. Molecular dynamics simulations of GILs are performed, utilizing the proposed force field. Derived structures approach the accuracy of QM/MM structures, and a previously reported crystal structure remains stable throughout the simulations. Mass densities are reproduced with a deviation of only 1.4% from experimental data. The calculated melting point of a GIL crystal deviates only 8% from the measured value. Self-diffusion coefficients of various GILs are reported, and a comparison with a diffusion coefficient derived from experimental data indicates that the values are within a reasonable range. We observe that the melting point of a GIL crystal was lowered up to 60 K and that diffusion coefficients are substantially increased by a factor of up to 3.5 upon consideration of charge transfer and polarization. The results demonstrate that liquid-phase partial charges are capable of improving the quality of ionic liquid force fields substantially and that their utilization led to a model that can be applied to predict structural, energetic, and dynamic properties of GILs.     

Molecular simulation study of some thermophysical and transport properties of triazolium-based ionic liquids

Results of a molecular dynamics study of several triazolium-based ionic liquids are reported. Triazolium cations include 1,2,4-triazolium, 1,2,3-triazolium, 4-amino-1,2,4-triazolium, and 1-methyl-4-amino-1,2,4-triazolium. Each cation was paired with a nitrate or perchlorate anion. These materials are part of a class of ionic compounds that have been synthesized recently but for which little physical property data are available. Properties of the more common ionic liquid, 1-n-butyl-3-methylimidazolium nitrate, are also computed and compared with the properties of the triazolium-based compounds. A molecular mechanics force field was developed for these materials using a mix of ab initio calculations and parameter fitting using the molecular compound 1H-1,2,4-triazole as a basis for the triazolium cations. Liquid-phase properties that were computed include heat capacities, cohesive energy densities, gravimetric densities/molar volumes as a function of temperature and pressure, self-diffusivities, rotational time constants, and various pair correlation functions. In the solid phase, heat capacities and lattice parameters were computed. Of all of these properties, only lattice parameters have been measured experimentally (and only for four of the triazolium compounds). The agreement with the experimental crystal structures was good. When compared with that of the imidazolium-based ionic liquid, the triazolium-based materials have much smaller molar volumes, higher cohesive energy densities, and larger specific heat capacities. They also tend to be less compressible, have a higher gravimetric density, and have faster rotational dynamics but similar translational dynamics.     

Molecular simulation study of dynamical properties of room temperature ionic liquids with carbon pieces

Room temperature ionic liquids (RTILs) with dispersed carbon pieces exhibit distinctive physiochemical properties. To explore the molecular mechanism, RTILs/carbon pieces mixture was investigated by molecular dynamics (MD) simulation in this work. Rigid and flexible carbon pieces in the form of graphene with different thicknesses and carbon nanotubes in different sizes were dispersed in a representative RTIL 1-butyl-3-methyl-imidazolium dicyanamide ([Bmim][DCA]). This study demonstrated that the diffusion coefficients of RTILs in the presence of flexible carbons are similar to those of bulk RTILs at varying temperatures, which is in contrast to the decreased diffusion of RTILs in the presence of rigid carbons. In addition, interfacial ion number density at rigid carbon surfaces was higher than that at flexible ones, which is correlated with the accessible external surface area of carbon pieces. The life time of cation-anion pair in the presence of carbon pieces also exhibited a dependence on carbon flexibility. RTILs with dispersed rigid carbon pieces showed longer ion pair life time than those with flexible ones, in consistence with the observation in diffusion coefficients. This work highlights the necessity of including the carbon flexibility when performing MD simulation of RTILs in the presence of dispersed carbon pieces in order to obtain the reliable dynamical and interfacial structural properties.     

Ab initio molecular dynamics simulation of ionic liquids

Ab initio Car-Parinnello molecular dynamics is used to simulate the structure and the dynamics of 1-butyl-3-methylimidazolium iodide ([bmim]I) ionic liquid at 300 K. Site-site pair correlation functions reveal that the anion has a strong interaction with any three C-H's of the imidazolium ring. The ring bends over and wraps around the anion such that the two nitrogen atoms take a distance to the anion. Electron donating butyl group contributes the electronic polarization in addition to geometrical (out-of-plane) polarization of the ring due to the liquid environment. This facilitates bending of the ring along the axis passing through nitrogen atoms. The average bending angle depends largely on the alkyl chain length and slightly on the anion type. Redistribution of electron density over the ring caused by the electron donating alkyl group provides additional independent evidence to the instability of lattice structure, hence the low melting point of the ionic liquid. Simulated viscosity and diffusion coefficients of [bmim]I are in quite agreement with the experiments.     

Molecular ions of ionic liquids in the gas phase

Ionic liquids form neutral ion pairs (CA) upon evaporation. The softness of the gas-phase ionization of field ionization has been used to generate "molecular ions," CA(+*), of ionic liquids, most probably by neutralization of the anion. In detail, 1-ethyl-3-methylimidazolium-thiocyanate, [C(6)H(11)N(2)](+) [SCN](-), 1-butyl-3-methylimidazolium-tricyanomethide, [C(8)H(15)N(2)](+) [C(4)N(3)](-), N-butyl-3-methylpyridinium-dicyanamide, [C(10)H(16)N](+) [C(2)N(3)](-), and 1-butyl-1-methylpyrrolidinium-bis[(trifluormethyl)sulfonyl]amide, [C(9)H(20)N](+) [C(2)F(6)NO(4)S(2)](-) were used. The assignment as CA(+*) ions, which has been confirmed by accurate mass measurements and misassignments due to thermal decomposition of the ionic liquids, has been ruled out by field desorption and electrospray ionization mass spectrometry of the residues.     

Molecular dynamics simulations of the surface tension of ionic liquids

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.     

Molecular dynamic simulation methods for anisotropic liquids

Methods of molecular dynamics simulations for anisotropic molecules are presented. The new methods, with an anisotropic factor in the cell dynamics, dramatically reduce the artifacts related to cell shapes and overcome the difficulties of simulating anisotropic molecules under constant hydrostatic pressure or constant volume. The methods are especially effective for anisotropic liquids, such as smectic liquid crystals and membranes, of which the stacks of layers are compressible (elastic in direction perpendicular to the layers) while the layer itself is liquid and only elastic under uniform compressive force. The methods can also be used for crystals and isotropic liquids as well.     

Raman spectra of ionic liquids: interpretation via computer simulation

Theoretical Raman spectra of the complex-forming ionic liquids LaCl3 and ScCl3, derived from molecular dynamics computer simulations, are presented. These simulations, which use polarizable ion interaction models, have previously been shown to predict structural properties in excellent agreement with diffraction experiments. The dependence of the polarizability of the melt on the ionic positions, which determines the Raman spectrum through the time dependence of the polarizability correlation function, is modeled on the basis of ab initio electronic structure calculations carried out on alkali chlorides. New simulation techniques are introduced in order to allow the spectrum to be calculated with acceptable statistics. The calculated spectra are in semiquantitative agreement with experimental data. The distinctive bands which appear in the spectra of such complex melts are linked to the vibrations of the transient coordination complexes which form in these melts and new interpretations for the origin of several well-known features are proposed. The simulations thus enable a link between the structure of a melt as perceived through Raman spectroscopy and through diffraction experiments to be made.     

Molecular simulation and experimental insights into the separation of dimethyl carbonate and isopropyl alcohol using imidazolium ionic liquids

In the process of polymerizing vinylidene fluoride, a mixture of isopropyl alcohol (IPA) and dimethyl carbonate (DMC) was generated. DMC is hard to separate from IPA since the two substances can form an azeotrope. Liquid-liquid extraction is often used to separate azeotropic systems due to their advantages of energy-saving and low cost. In this work, ionic liquids were selected as the extractants to separate the azeotrope. To obtain the suitable ILs, we calculated the σ-profiles of DMC, IPA, [H₂PO₄]⁻, [HSO₄]⁻, [NO₃]⁻, [PF₆]⁻, [BF₄]⁻, and [OTF]⁻by the COSMO model, and further revealed the extraction mechanism through the electrostatic potential (ESP), reduced density gradient (RDG) and interaction energies analysis. To verify the theoretical calculation results, the liquid-liquid equilibrium (LLE) data of DMC-IPA-[Emim][H₂PO₄]/[Emim][HSO₄]/[Emim][NO₃] was measured at 298.15 K and 101.32 kPa, the experimental results are consistent with the theoretical calculation results. Meanwhile, the NRTL interaction parameters were regressed using the least squares method, and the RMSD values were all less than 0.009, which indicated the NRTL model was suitable for this system. The combination of experiment and theoretical calculation provides a new idea for separating the DMC-IPA system using ILs.     

Molecular solutions of cellulose in mixtures of ionic liquids with pyridine

New systems suitable for determination of molecular characteristics of cellulose, mixtures of ionic liquids based on 1-n-alkyl-3-methylimidazolium with pyridine, were found. In ionic liquid-pyridine mixtures, cellulose is dispersed on the molecular level. The cellulose-ionic liquid-pyridine systems with 1-ethyl-3-methylimidazolium acetate and 1-ethyl-3-methylimidazolium diethyl phosphate are stable in time. The dynamic viscosity and refractive index of the mixtures can be controlled by varying the ionic liquid to pyridine ratio. The viscometric and dynamooptical properties of cellulose in these mixtures were compared with those in Cadoxen.     

Molecular dynamics simulation studies of the influence of imidazolium structure on the properties of imidazolium/azide ionic liquids

Atomistic molecular dynamics simulations were performed on 1-butyl-3-methyl-imidazolium azide [bmim][N(3)], 1-butyl-2,3-dimethylimidazolium azide [bmmim][N(3)], and 1-butynyl-3-methyl-imidazolium azide [bumim][N(3)] ionic liquids. The many-body polarizable APPLE&P force field was augmented with parameters for the azide anion and the bumim cation. Good agreement between the experimentally determined and simulated crystal structure of [bumim][N(3)] as well as the liquid-state density and ionic conductivity of [bmmim][N(3)] were found. Methylation of bmim (yielding bmmim) resulted in dramatic changes in ion structuring in the liquid and slowing of ion motion. Conversely, replacing the butyl group of bmim with the smaller 2-butynyl group resulted in an increase of ion dynamics.     

Microwave-assisted separation of ionic liquids from aqueous solution of ionic liquids

Microwave-assisted separation has been applied to recover ionic liquid (IL) from its aqueous solution as an efficient method with respect to time and energy compared to the conventional vacuum distillation. Hydrophilic ILs such as 1-butyl-3-methylimidazolium tetrafluoroborate ([Bmim][BF(4)]), 1-butyl-3-methylimidazolium trifluoromethanesulfonate ([Bmim][TfO]) and 1-ethyl-3-methylimidazolium methylsulfate ([Emim][MS]) could be recovered in 6 min from the mixture of ILs and water (1:1, w/w) under microwave irradiation at constant power of 10 W while it took at least 240 min to obtain ILs containing same water content (less than 0.5 wt%) by conventional vacuum oven at 363.15 K with 90 kPa of vacuum pressure. Energy consumptions per gram of evaporated water from the homogeneous mixture of hydrophilic ILs and water (1:1, w/w) by microwave-assisted separation were at least 52 times more efficient than those in conventional vacuum oven. It demonstrated that microwave-assisted separation could be used for complete recovery of ILs in sense of time and energy as well as relevant purity.     

Molecular aspects of glucose dehydration by chromium chlorides in ionic liquids

A combined experimental and computational study of the ionic‐liquid‐mediated dehydration of glucose and fructose by Cr^II and Cr^III chlorides has been performed. The ability of chromium to selectively dehydrate glucose to 5‐hydroxymethylfurfural (HMF) in the ionic liquid 1‐ethyl‐3‐methyl imidazolium chloride does not depend on the oxidation state of chromium. Nevertheless, Cr^III exhibits higher activity and selectivity to HMF than Cr^II. Anhydrous CrCl₂ and CrCl₃?6 H₂O readily catalyze glucose dehydration with HMF yields of 60 and 72 %, respectively, after 3 h. Anhydrous CrCl₃ has a lower activity, because it only slowly dissolves in the reaction mixture. The transformation of glucose to HMF involves the formation of fructose as an intermediate. The exceptional catalytic performance of the chromium catalysts is explained by their unique ability to catalyze glucose to fructose isomerization and fructose to HMF dehydration with high selectivity. Side reactions leading to humins by means of condensation reactions take predominantly place during fructose dehydration. The higher HMF selectivity for Cr^III is tentatively explained by the higher activity in fructose dehydration compared to Cr^II. This limits the concentration of intermediates that are involved in bimolecular condensation reactions. Model DFT calculations indicate a substantially lower activation barrier for glucose isomerization by Cr^III compared to Cr^II. Qualitatively, glucose isomerization follows a similar mechanism for Cr^II and Cr^III. The mechanism involves ring opening of D ‐glucopyranose coordinated to a single Cr ion, followed by a transient self‐organization of catalytic chromium complexes that promotes the rate‐determining hydrogen‐shift step.     

Molecular solvent effect on the acidity constant of protic ionic liquids

In this work how the microscopic properties of a molecular solvent affect the chemical environment of the protic ionic liquids (PILs) was analyzed. Using Reichardt’s dye as indicator of acidity, new acidity constant values for eight PILs (pK_a^PILs) were determined by spectrophotometric titration. Modifying the character hydrogen bonding donor of the molecular solvent it is possible to handle the PIL acid strength. Thus, we can turn basic PILs into acidic ones thereby the molecular solvent could be used as ‘additive’ for PILs, which allowed us to tune PILs design.