A Force Field for Guanidinium-Based Ionic Liquids That Utilizes the Electron Charge Distribution of the Actual Liquid: A Molecular Simulation Study

Computational approaches to understanding reaction outcomes of organic processes in ionic liquids

The utility of using a combined experimental and computational approach for understanding ionic liquid media, and their effect on reaction outcome, is highlighted through a number of case studies.

How does the dehydration change the host-guest association under homogeneous and heterogeneous conditions?

In this study, the thermodynamic properties of association of some inorganic ions (ClO4(-) and SO4(2-)) with β-cyclodextrins (β-CD) in aqueous solution are determined under both free β-CD and surface confined β-CD conditions using atomistic simulations. The potential of mean force (PMF) is calculated as a function of the environment and the thermodynamic properties of association are deduced by integrating the free energy profiles. No inclusion complex between SO4(2-) and β-CD is detected. Nevertheless, the PMF curve obtained for gold-confined CD seems to evidence a small minimum at a larger separation distance that shows specific interactions such as hydrogen bonding outside the cavity. As concerns ClO4(-), our simulations reveal the formation of an inclusion complex with free β-CD in perfect agreement with the available experimental results. Nevertheless, we do not detect any formation of the host-guest inclusion complex under heterogeneous conditions. Finally, the differences observed as a function of the anions are interpreted through an atomistic description. The general trend of weaker complex stabilities with the increasing free energy of hydration of the anions is found in homogeneous systems.

Variation and decomposition of the partial molar volume of small gas molecules in different organic solvents derived from molecular dynamics simulations

The partial molar volumes, V(i), of the gas solutes H2, CO, and CO2, solvated in acetone, methanol, heptane, and diethylether are determined computationally in the limit of infinite dilution and standard conditions. Solutions are described with molecular dynamics simulations in combination with the OPLS-aa force field for solvents and customized force field for solutes. V(i) is determined with the direct method, while the composition of V(i) is studied with Kirkwood-Buff integrals (KBIs). Subsequently, the amount of unoccupied space and size of pre-formed cavities in pure solvents is determined. Additionally, the shape of individual solvent cages is analyzed. Calculated V(i) deviate only 3.4 cm(3) mol(-1) (7.1%) from experimental literature values. Experimental V(i) variations across solutions are reproduced qualitatively and also quantitatively in most cases. The KBI analysis identifies differences in solute induced solvent reorganization in the immediate vicinity of H2 (<0.7 nm) and solvent reorganization up to the third solvation shell of CO and CO2 (<1.6 nm) as the origin of V(i) variations. In all solutions, larger V(i) are found in solvents that exhibit weak internal interactions, low cohesive energy density and large compressibility. Weak internal interactions facilitate solvent displacement by thermal solute movement, which enhances the size of solvent cages and thus V(i). Additionally, attractive electrostatic interactions of CO2 and the solvents, which do not depend on internal solvent interactions only, partially reversed the V(i) trends observed in H2 and CO solutions where electrostatic interactions with the solvents are absent. More empty space and larger pre-formed cavities are found in solvents with weak internal interactions, however, no evidence is found that solutes in any considered solvent are accommodated in pre-formed cavities. Individual solvent cages are found to be elongated in the negative direction of solute movement. This wake behind the moving solute is more pronounced in case of mobile H2 and in solvents with weaker internal interactions. However, deviations from a spherical solvent cage shape do not influence solute-solvent radial distribution functions after averaging over all solvent cage orientations and hence do not change V(i). Overall, the applied methodology reproduces V(i) and its variations reliably and the used V(i) decompositions identify the underlying reasons behind observed V(i) variations.

Extraction of tryptophan with ionic liquids studied with molecular dynamics simulations

Extraction of amino acids from aqueous solutions with ionic liquids (ILs) in biphasic systems is analyzed with molecular dynamics (MD) simulations. Extraction of tryptophan (TRP) with the imidazolium-based ILs [C(4)mim][PF(6)], [C(8)mim][PF(6)], and [C(8)mim][BF(4)] are considered as model cases. Solvation free energies of TRP are calculated with MD simulations and thermodynamic integration in combination with an empirical force field, whose parametrization is based on the liquid-phase charge distribution of the ILs. Calculated solvation free energies reproduce successfully all observed experimental trends according to the previously reported partition of TRP between water and IL phases. Water is present in ILs as a cosolvent, due to direct contact with the aqueous phase during extraction, and is found to play a major role in the extraction of TRP. Water improves solvation of cationic TRP by 7.8 and 5.1 kcal/mol in [C(4)mim][PF(6)] and [C(8)mim][PF(6)], respectively, which is in the case of [C(4)mim][PF(6)] sufficient to extract TRP. Extraction in [C(8)mim][PF(6)] is not feasible, since the hydrophobic octyl groups of the cations limit the water concentration in the IL. The solvation of cationic TRP is 2.4 kcal/mol less favorable in [C(8)mim][PF(6)] than in [C(4)mim][PF(6)]. Water improves the solvation of TRP in ILs mostly through dipole-dipole interactions with the polar backbone of TRP. Extraction is most efficient with [C(8)mim][BF(4)], where hydrophilic BF(4)(-) anions substantially increase the water concentration in the IL. Additionally, stronger direct electrostatic interactions of TRP with BF(4)(-) anions improve its solvation in the IL further. The solvation of cationic TRP in [C(8)mim][BF(4)] is 3.4 kcal/mol more favorable than in [C(8)mim][PF(6)]. Overall, the extractive power of the ILs correlates with the water saturation concentration of the IL phase, which in turn is determined by the hydrophilicity of the constituting ions. The results of this work identify relations between the extraction performance of ILs and the basic chemical properties of the ions, which provide guidelines that could contribute to the design of improved novel ILs for amino acid extraction.

Free energy calculations in electroactive self-assembled monolayers (SAMs): impact of the chain length on the redox reaction

The free energy approach is used to study the effect of the relative chain length of the two constituents of electroactive self-assembled monolayers (SAMs) on gold. In this study, the ferrocene groups are exposed to the electrolyte solution. This situation is achieved by using shorter diluent alkanethiol chains. To this end, the mixed monolayers formed by the self-assembly of 11-ferrocenylundecanethiol and butanethiol FcC(11)S/C(4)S and of 6-ferrocenylhexanethiol and butanethiol FcC(6)S/C(4)S onto a gold surface are studied. Calculation of enthalpy and entropy differences are also performed using molecular simulations. Additionally, the electrochemical signatures of these systems are determined to allow a direct comparison with our calculations. The thermodynamic properties are discussed in terms of enthalpy and entropy changes. Two effects account for the thermodynamic behavior. The first one involves the ion pairing between the ferrocenium group and the perchlorate anion. The second one concerns the desolvation of the first hydration shell of the anions. Finally, this work is also completed with a microscopic description associated with an energy characterization of these SAMs as a function of the surface coverage under conditions close to experiments.

The structure of geminal imidazolium bis(trifluoromethylsulfonyl)amide ionic liquids: a theoretical study of the gas phase ionic complexes

In this work we report molecular mechanics and ab initio calculations on the geminal di-imidazolium bis(trifluoromethylsulfonyl)amide ionic liquid in the gas phase. We report the likely energetically preferred geometries of the ionic complex and its main features in terms of charge distribution, electronic density, structure, and energetics. We find that the gas phase structure of the ionic complex is quite compact and that the alkyl chain connecting the two imidazolium charged rings is strongly bent in order to maximize their electrostatic interactions with the two anions.