Rotational dynamics of a diatomic solute in the room-temperature ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate

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.

Coarse-grained simulation of the translational and rotational diffusion of globular proteins by dissipative particle dynamics

With simplified interactions and degrees of freedom, coarse-grained (CG) simulations have been successfully applied to study the translational and rotational diffusion of proteins in solution. However, in order to reach larger lengths and longer timescales, many CG simulations employ an oversimplified model for proteins or an implicit-solvent model in which the hydrodynamic interactions are ignored, and thus, the real kinetics are more or less unfaithful. In this work, we develop a CG model based on the dissipative particle dynamics (DPD) that can be universally applied to different types of proteins. The proteins are modeled as a group of rigid DPD beads without conformational changes. The fluids (including solvent and ions) are also modeled as DPD beads. The electrostatic interactions between charged species are explicitly considered by including charge distributions on DPD particles. Moreover, a surface friction between the protein and fluid beads is applied to control the slip boundary condition. With this model, we investigate the self-diffusion of a single globular protein in bulk solution. The translational and rotational diffusion coefficients of the protein can be tuned by the surface frictional constant to fit the predictions of the Stokes-Einstein (SE) relation. We find that both translational and rotational diffusion coefficients that meet with the prediction of the SE relation based on experimental results of the hydrodynamic radius are reached at almost the same frictional constant for different types of proteins. Such scaling behavior indicates that the model can be applied to simulate the translational and rotational diffusion together for various types of proteins.

Dynamics in the BMIM PF6/acetonitrile mixtures observed by femtosecond optical Kerr effect and molecular dynamics simulations

We have performed the measurements of the optical Kerr effect signal time evolution up to 4 ns for a mixture of 1-alkyl-3-methyl-imidazolium hexafluorophosphate (BMIM PF₆) ionic liquid and acetonitrile in the whole mole fractions range. The long delay line in our experimental setup allowed us to capture the complete reorientational dynamics of the ionic liquid. We have analysed the optical Kerr effect signal in the time and frequency domains with help of molecular dynamics simulations. In our approximation of the slow picosecond dynamics with a multi-exponential decay, we distinguish three relaxation times. The highest two are assigned to the reorientation of the cation and acetonitrile molecules that are in the vicinity of the imidazolium ring. The third one is recognized as originating from cation rotations and reorientation of acetonitrile molecules in the bulk or in the vicinity of the aliphatic chains of the cation. With help of the simulation we interpret the intermolecular band in the reduced spectral density, obtained from Kerr signal, as follows: its low-frequency side results from oscillations of one of the components in the cage formed by its neighbors, while the high-frequency side is attributed to the librations of the cation and acetonitrile molecule as well as the intermolecular oscillations of system components involved in specific interactions. We use this assignment and concentration dependence of the spectra obtained from velocity and angular velocity correlations to explain the mole fraction dependence of Kerr reduced spectral density.

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.

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.

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.

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.

Fluorescence, Phosphorescence, and Delayed Fluorescence of Benzil in Imidazolium Ionic Liquids

Temperature dependence of the emission behaviour of benzil has been studied in three imidazolium ionic liquids differing in their polarity and viscosity. Room temperature absorption and steady-state emission spectra suggest that the ground and excited state conformers of benzil in ionic liquids are similar to those in conventional organic solvents. The non-degassed solutions of benzil in ionic liquids show phosphorescence at room temperature in contrast to conventional solvents where phosphorescence is commonly observed in degassed conditions. This study reveals that a thermally activated reverse intersystem crossing (T1↝S1) process is responsible for the drastic change in phosphorescence intensity with temperature in ionic liquids. The rate constant () of this process is found to be dependent on the polarity of the media and is 5 times higher in most polar ionic liquids. The evidence of the presence of multiple conformers of benzil in frozen conditions is obtained from the excitation wavelength dependence of the phosphorescence spectra.

Free Energy and Dynamics of Electron-Transfer Reactions in a Room Temperature Ionic Liquid

Reaction free energetics and dynamics of unimolecular electron-transfer processes in ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI+PF6-) are investigated via molecular dynamics computer simulations employing a model diatomic solute and compared with those in aprotic acetonitrile. Using the free energy perturbation method, diabatic free energy curves relevant to charge separation and recombination processes are studied over a wide range of the reaction coordinate. The diabatic curves are found to vary with the solute charge distribution, especially in EMI+PF6-. Nevertheless, if the free energy of reaction is not that substantial, the Marcus free energy relationship holds reasonably well, provided that the reorganization free energy averaged between the reactant and product states is employed. The effective polarity, measured as solvation-induced stabilization of dipolar solutes, is higher for EMI+PF6- than for acetonitrile, consonant with many solvatochromic measurements. Thus, in the normal regime, activation barriers for charge separation and recombination reactions are, respectively, lower and higher in EMI+PF6- than in acetonitrile. The influence of solvent dynamics on reaction kinetics through modulations of activation, deactivation, and barrier crossing is analyzed. Even though overall solvent relaxation dynamics in EMI+PF6- are considerably slower than those in acetonitrile, the deviation of the rate constant from the transition state theory predictions is found to be small for both solvents. Implications of this finding for other reactions in ionic liquids are briefly discussed.

Measurements of the Complete Solvation Response in Ionic Liquids

Dynamic Stokes shift measurements of the solvatochromic probe trans-4-dimethylamino-4'-cyanostilbene were used to measure the solvation response of five imidazolium and one pyrrolidinium ionic liquid at 25 degrees C. The Kerr-gated emission and time-correlated single-photon-counting techniques were used to measure spectral dynamics occurring over the time ranges of 100 fs-200 ps and 50 ps-5 ns, respectively, and a combination of data sets from these two techniques enabled observation of the complete solvation response. Observed response functions were found to be biphasic, consisting of a sub-picosecond component of modest (10-20%) amplitude and a dominant slower component relaxing over times of a few picoseconds to several nanoseconds. The faster component could be correlated to inertial characteristics of the constituent ions, and the slower component to solvent viscosity. Dielectric continuum calculations of the sort previously used to predict solvation dynamics in dipolar liquids were shown to work poorly for predicting the response in these ionic liquids.

Effects of Solute Electronic Polarizability on Solvation in a Room-Temperature Ionic Liquid

The effects of solute polarizability on solvation and solute transport in the room-temperature ionic liquid 1-ethyl-3-methylimidazolium hexafluorophosphate (EMI+PF(6)-) are investigated via molecular dynamics simulations. A valence-bond description is employed to account for the instantaneous adjustment of the solute electronic charge distribution to the fluctuating solvent environment. It is found that the ultrafast inertial component of solvation dynamics becomes slower as the solute polarizability grows. Moreover, its contribution to overall solvent relaxation becomes reduced with increasing polarizability, especially in the case of nonequilibrium solvation dynamics. Overall, the inclusion of the solute electronic polarizability in the simulations improves the agreement with time-dependent Stokes shift measurements.