A molecular dynamics computer simulation study of room-temperature ionic liquids. I. Equilibrium solvation structure and free energetics

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts - often in collaboration with experimental groups - to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO2 capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.

A Computer Simulation Study of Thermal and Mechanical Properties of Poly(Ionic Liquid)s

Thermal and mechanical properties of poly(ionic liquid)s (PILs), an epoxidized ionic liquid-amine network, are studied via molecular dynamics simulations. The poly(ionic liquid)s are designed with two different ionic liquid monomers, 3-[2-(Oxiran-2-yl)ethyl]-1-{4-[(2-oxiran-2-yl)ethoxy]phenyl}imidazolium (EIM2) and 1-{4-[2-(Oxiran-2-yl)ethyl]phenyl}-3-{4-[2-(oxiran-2-yl)ethoxy]benzyl}imidazolium (EIM1), each of which is networked with tris(2-aminoethyl)amine, paired with different anions, bis(trifluoromethanesulfonyl)imide (TFSI⁻) and chloride (Cl⁻). We investigate how ionic liquid monomers with high ionic strength affect structures of the cross-linked polymer networks and their thermomechanical properties such as glass transition temperature (T_g) and elastic moduli, varying the degree of cross-linking. Strong electrostatic interactions between the cationic polymer backbone and anions build up their strong structures of which the strength depends on their molecular structures and anion size. As the anion size decreases from TFSI⁻ to Cl⁻, both T_g and elastic moduli of the PIL increase due to stronger electrostatic interactions present between their ionic moieties, making it favorable for the PIL to organize with stronger bindings. Compared to the EIM2 monomer, the EIM1 monomers and TFSI⁻ ions generate a PIL with higher T_g and elastic moduli. This attributes to the less flexible structure of the EIM1 monomer for the chain rotation, in which steric hindrance by ring moieties in the EIM1-based PIL enhances their structural rigidity. The π-π stacking structures between the rings are found to increase in EIM1-based PIL compared to the EIM2-based one, which becomes stronger with smaller Cl⁻ ion rather than TFSI⁻. The effect of the degree of the cross-linking on thermal and mechanical properties is also examined. As the degree of cross-linking decreases from 100% to 60%, T_g also decreases by a factor of 10–20%, where the difference among the given PILs becomes decreased with a lower degree of cross-linking. Both the Young’s (E) and shear (G) moduli of all the PILs decrease with degree of cross-linking, which the reduction is more significant for the PIL generated with EIM2 monomers. Transport properties of anions in PILs are also studied. Anions are almost immobilized globally with very small structural fluctuations, in which Cl⁻ presents lower diffusivity by a factor of ~2 compared to TFSI⁻ due to their stronger binding to the cationic polymer backbone.

To unravel the connection between the non-equilibrium and equilibrium solvation dynamics of tryptophan: success and failure of the linear response theory of fluorescence Stokes shift

The connections between the non-equilibrium solvation dynamics upon optical transitions and the system's equilibrium fluctuations are explored in aqueous liquid. Linear response theory correlates time-dependent fluorescence with the equilibrium time correlation functions. In the previous work [T. Li, J. Chem. Theory Comput., 2017, 13, 1867], Stokes shift was explicitly decomposed into the contributions of various order time correlation functions on the excited state surface. Gaussian fluctuations of the solute-solvent interactions validate linear response theory. Correspondingly, the deviation of the Gaussian statistics causes the inefficiency of linear response evaluation. The above mechanism is thoroughly tested in this study. By employing molecular simulations, multiple non-equilibrium processes, not necessarily initiated from the ground state equilibrium minimum, were examined for tryptophan. Both the success and failure of linear response theory are found for this simple system and the mechanism is analyzed. These observations, assisted by the width dynamics, the initial state linear response approach, and the variation of the solvation structures, integrally verify the virtue of the excited state Gaussian statistics on the dynamics of Stokes shift.

CS2 capture in the ionic liquid 1-alkyl-3-methylimidazolium acetate: reaction mechanism and free energetics

Reaction pathways for CS2 and COS in the ionic liquid, 1-ethyl-3-methylimidazolium (EMI+) acetate (OAc-), are studied using the ab initio self-consistent reaction field theory (SCRF) and molecular dynamics (MD) computer simulations. It is found that while CS2 converts to COS nearly at the 100% level through S/O exchange with acetate, both conversion and capture processes are kinetically possible for COS, yielding CO2/thioacetate and 1-ethyl-3-methylimidazole-2-thiocarboxylate (EMI-COS)/acetic acid as reaction products, respectively. These findings are in excellent agreement with recent experimental observations in the closely related 1-butyl-3-methylimidazolium acetate (BMI+OAc-) ionic liquid system. Constrained ab initio MD indicates that the capture reaction of COS (and CS2 if allowed) proceeds in a concerted fashion; viz., proton transfer from EMI+ to OAc- and carboxylation of EMI+ by COS (and CS2) occur concurrently, analogous to the concerted pathway proposed recently for CO2 capture in the imidazolium acetate ionic liquid family. As N-heterocyclic carbene (NHC) is not required, the concerted mechanism is fully consistent with the experimental fact that NHC has not been detected directly in this ionic liquid family. Computational analysis further predicts that if NHC would be present in the ionic liquid, it would react with CS2 and produce 1-ethyl-3-imidazole-2-dithiocarboxylate, prior to the conversion of CS2 to COS. Since such a dithiocarboxylate compound was not detected experimentally, the present analysis lends support to the view that NHC is not formed in the pure imidazolium acetate ionic liquid family.

A molecular dynamics study of the ionic liquid, choline acetate

Structural and dynamic properties of the ionic liquid (IL) choline acetate are studied using molecular dynamics (MD) simulations. The hydroxyl group of choline shows significant hydrogen-bonding interactions with the oxygen atoms of acetate. Nearly all choline cations are found to form a hydrogen bond with acetate anions at 400 K, while about 67% of cations participate in hydrogen-bonding interactions at 600 K. At 400 K, subdiffusive and prominent non-Gaussian behavior persist for t > 10 ns. At 600 K, the usual diffusion regime is obtained after a few hundred ps of subdiffusive behavior. Analysis of reorientational motions of acetate ions, particularly those of their short axes, indicates a high degree of dynamic heterogeneity, in agreement with previous work on different IL systems.

Effect of a Tertiary Butyl Group on Polar Solvation Dynamics in Aqueous Solution: A Computational Approach

The current computational study investigates the changes in solvation dynamics of water when introducing hydrophobic side chains to the molecular probe N-methyl-6-oxyquinolinium betaine. High-precision transient fluorescence and absorption measurements published in the companion article (10.1021/acs.jpcb.7b05031) revealed an influence of hydrophobic side chain alterations on the observed solvation dynamics of a chromophore in water. As the influence of shape, size, and structure of chromophores on the time-dependent Stokes shift was so far thought to play a role only in slowly rotating solvents compared to the solute or if the hydrogen bonding ability of the solute changes, this finding is quite unexpected. Analysis of the time-dependent Stokes shift obtained from nonequilibrium simulations corroborates experimental retardation factors and activation energies, and indicates that solute motion, namely vibration, is mainly responsible for the observed retardation of solvation dynamics. The faster dynamics around the smaller chromophore is in fact achieved by some normal modes located at the pyridinium part of the chromophore. Rotation also contributes to a very small extent to hydration dynamics, but for small and large derivatives alike. Local residence times furthermore reveal slight retardations in the first solvent shell around the chromophores. The current picture of the solute acting as a passive molecular probe therefore needs to be revised even for solvents like water.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Excitation-energy dependence of solvation dynamics in room-temperature ionic liquids

Influence of the excitation energy of a probe solute molecule on its solvation dynamics and emission spectrum in 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI(+)PF6 (-)) is studied via molecular dynamics simulations using a coarse-grained model description. By exciting the probe at different energies, each with an extremely narrow distribution, ensuing solvent relaxation and its dynamic variance are monitored using the isoconfigurational ensemble method. Resulting Stokes shift function, S(t), indicates that long-time solvent relaxation becomes slower with the decreasing excitation energy and approaches the equilibrium correlation function, C(t), of solvent fluctuations. This suggests that the system excited at the red-edge of the spectrum observes linear response better than that at the blue-edge. A detailed analysis of nonequilibrium trajectories shows that the effect of initial configurations on variance of relaxation dynamics is mainly confined to short times; it reaches a maximum around 0.1 ≲ t ≲ 1 ps and diminishes as time further increases. The influence of the initial velocity distribution, on the other hand, tends to grow with time and dominates the long-time variations of dynamics. The emission spectrum shows the red-edge effect in accord with previous studies.

Solvent effects of 1-ethyl-3-methylimidazolium acetate: solvation and dynamic behavior of polar and apolar solutes

We study the solvation properties of the ionic liquid 1-ethyl-3-methylimidazolium acetate ([eMIM]⁺[ACE]⁻) and the resulting dynamic behavior for differently charged model solutes at room temperature via atomistic molecular dynamics (MD) simulations of 500 ns length.

Dynamic propensity as an indicator of heterogeneity in room-temperature ionic liquids

Dynamic propensity of an RTIL system exhibits broad and asymmetric distributions, and spatial patterns of the dynamic propensity and mobility distribution are shown.

Double-layer formation of [Bmim][PF6] ionic liquid triggered by surface negative charge

Applications of ionic liquids (ILs) in electrified interfaces and electrochemical systems require insight into the molecular-level structure and properties of the interfacial ILs. Using atomistic molecular dynamics (MD) simulations, we show here that a new double-layer stacking formation of the [Bmim][PF(6)] IL can be triggered by the surface negative charge. We also found that the double-layer formation induced by the surface charge thoroughly extended into the bulk phase, implying a strong unscreened ion effect in our IL system. Further study indicated that the double-layer formation in the bulk phase was due to a rapid structural transition. Different IL formations, including the conventional adsorption layer and the double-layer formation, can be achieved in sequence by increasing the surface negative charge. Moreover, the diffusion ability of the new double-layer formation in the bulk phase is much lower when compared to that observed in its original uncharged condition. The structure and properties of the ILs formation may be attributed to the tail-tail aggregation hypothesis of the nonpolar domain in the IL.

Charge effect on the diffusion coefficient and the bimolecular reaction rate of diiodide anion radical in room temperature ionic liquids

The diffusion coefficients of diiodide anion radical, I(2)(-), in room temperature ionic liquids (RTILs) were determined by the transient grating (TG) method using the photochemical reaction of iodide. The diffusion coefficients we obtained were larger in RTILs than the theoretical predictions by the Stokes-Einstein relation, whereas both values are similar in conventional solvents. By comparison with the diffusion coefficients of neutral molecules, it was suggested that the Coulomb interaction between I(2)(-) and constituent ions of RTILs strongly affects the diffusion coefficients. The bimolecular reaction rates between I(2)(-) were calculated by the Debye-Smoluchowski equation using the experimentally determined diffusion coefficients. These calculated reaction rate were much smaller than the experimentally determined rates (Takahashi, K.; et al. J. Phys. Chem. B 2007, 111, 4807), indicating the charge screening effect of RTILs.

Ionic liquids: Solvation ability and polarity

The role of ionic liquids (ILs) as solvents in chemistry is limited by the poor understanding of the solvation phenomenon in these media. The usual classification criteria used for molecular solvents through various experimental measurements fail to insert ILs into a univocal classification for ILs. Here, we first discuss the unsuitability of the usual interpretative scheme for molecular liquids and elucidate schematically the mechanism of solvation in ILs, pointing out the peculiarities that differentiate them with respect to molecular liquids. Second, we focus on reactivity and reaction kinetics in ILs, underlining the many problems that the complexity of these media reflects on the interpretation of kinetic data and some possible approaches to understand qualitatively the (often not trivial) kinetic problems for reactions performed in ILs.