MD Study of Stokes Shifts in Ionic Liquids: Temperature Dependence

Reorganization Energy of Electron Transfer in Ionic Liquids

Bandshape analysis of charge-transfer optical bands in room-temperature ionic liquids (ILs) was performed to extract the reorganization energy of electron transfer. Remarkably, the reorganization energies in ILs are close to those in cyclohexane. This result runs against common wisdom in the field since conducting ILs, which are characterized by an infinite static dielectric constant, and nonpolar cyclohexane fall to the opposite ends of the polarity scale based on their dielectric constants. Theoretical calculations employing structure factors of ILs from molecular dynamics simulations support the low values of the reorganization energy. Standard dielectric arguments do not apply to solvation in ILs, and nonergodic reorganization energies are required for a quantitative analysis.

Polarizability in ionic liquid simulations causes hidden breakdown of linear response theory

The validity of linear response theory (LRT) in computer simulations of solvation dynamics, i.e. the time-dependent Stokes shift, has been debated widely during the last decades. Since the use of LRT is computationally less expensive than the calculation of the true nonequilibrium response, it is often invoked for large systems exhibiting a particularly slow solvation response, e.g. ionic liquids. In the case of ionic liquids, LRT does not only need to capture the correct overall dynamics of the system, but also the contributions and timescales of the respective cation and anion movement. We show by large scale computer simulations that the contribution of the permanent dipoles to the solvation response obeys LRT to some extent, whereas the induced contributions in polarizable simulations lead to a failure of LRT for the respective ion contributions.

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.

Solvation dynamics in polar solvents and imidazolium ionic liquids: failure of linear response approximations

This study presents the large scale computer simulations of two common fluorophores, N-methyl-6-oxyquinolinium betaine and coumarin 153, in five polar or ionic solvents. The validity of linear response approximations to calculate the time-dependent Stokes shift is evaluated in each system. In most studied systems linear response theory fails. In ionic liquids the magnitude of the overall response is largely overestimated, and linear response theory is not able to capture the individual contributions of cations and anions. In polar liquids, the timescales of solvation dynamics are often not correctly reproduced. These observations are complemented by a detailed analysis of Gaussian statistics including higher order correlation functions, variance of the energy gap distribution and its time evolution. The analysis of higher order correlation functions was found to be not suitable to predict a failure of linear response theory. Further analysis of radial distribution functions and hydrogen bonds in the ground and excited state, as well as the time evolution of the number of hydrogen bonds after solute excitation reveal an influence of solvent structure in some of the studied systems.

CO2 capture in ionic liquid 1-alkyl-3-methylimidazolium acetate: a concerted mechanism without carbene

Ionic liquids (ILs) provide a promising medium for CO₂ capture. Recently, the family of ILs comprising imidazolium-based cations and acetate anions, such as 1-ethyl-3-methylimidazolium acetate (EMI⁺OAc^-), has been found to react with CO₂ and form carboxylate compounds. N-Heterocyclic carbene (NHC) is widely assumed to be responsible by directly reacting with CO₂ though NHC has not been detected in these ILs. Herein, a computational analysis of CO₂ capture in EMI⁺OAc^- is presented. Quantum chemistry calculations predict that NHC is unstable in a polar environment, suggesting that NHC is not formed in EMI⁺OAc^-. Ab initio molecular dynamics simulations indicate that an EMI⁺ ion "activated" by the approach of a CO₂ molecule can donate its acidic proton to a neighboring OAc^- anion and form a carboxylate compound with the CO₂ molecule. Analysis of this termolecular process indicates that the EMI⁺-to-OAc^- proton transfer and the formation of 1-ethyl-3-methylimidazolium-2-carboxylate occur essentially concurrently. Based on these findings, a novel concerted mechanism that does not involve NHC is proposed for CO₂ capture.

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.