A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics

Liquid Structure of Ionic Liquids with [NTf2]- Anions, Derived from Neutron Scattering

The liquid structure of three common ionic liquids (ILs) was investigated by neutron scattering for the first time. The ILs were based on the bis(trifluoromethanesulfonyl)imide anion, abbreviated in the literature as [NTf2]- or [TFSI]-, and on the following cations: 1-ethyl-3-methylimidazolium, [C2mim]+; 1-decyl-3-methylimidazolium, [C10mim]+; and trihexyl(tetradecyl)phosphonium, [P666,14]+. Comparative analysis of the three ILs confirmed increased size of nonpolar nanodomains with increasing bulk of alkyl chains. It also sheds light on the cation-anion interactions, providing experimental insight into strength, directionality, and angle of hydrogen bonds between protons on the imidazolium ring, as well as H-C-P protons in [P666,14]+, to oxygen and nitrogen atoms in the [NTf2]-. The new Dissolve data analysis package enabled, for the first time, the analysis of neutron scattering data of ILs with long alkyl chains, in particular, of [P666,14][NTf2]. Results generated with Dissolve were validated by comparing outputs from three different models, starting from three different sets of cation charges, for each of the three ILs, which gave convergent outcomes. Finally, a modified method for the synthesis of perdeuterated [P666,14][NTf2] has been reported, with the aim of reporting a complete set of synthetic and data processing approaches, laying robust foundations that enable the study of the phosphonium ILs family by neutron scattering.

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.

Linear Response Theory for Decomposition Energies of Stokes Shift in Proteins

In aqueous solution, fluorescence Stokes shift experiments monitor the relaxation of the solute-solvent interactions upon photon excitation of the solute chromophore. Linear response (LR) theory expects the identical dynamics between the Stokes shift and the system's spontaneous fluctuations. Whether this identity guarantees similar dynamics between the nonequilibrated and equilibrium processes for the decomposition energy of the Stokes shift is the main focus of this study. In our previous work [Li, T. J. Chem. Theory Comput. 2017, 13, 1867-1873], Stokes shift is properly correlated with various order time-correlation functions. As a continuation, its decomposition energy from the subsystem is further represented as the full summation of all of the cross-time correlation functions between the decomposition energy and the total solute-solvent interactions. Gaussian statistics of the total solute-solvent interactions ensure the same decay rates among the odd orders not only for the time-correlation functions but also for the cross-time correlation functions, validating the LR of the Stokes shift and the decompositions, respectively. The above mechanism is verified by molecular dynamics simulations in the protein Staphylococcus nuclease and is robust even as the decomposed energy associated with an individual residue exhibits typical non-Gaussian properties. Further examinations reveal the consistent molecular motions for a specific residue over the nonequilibrium and equilibrium processes, which are responsible for the nonequilibrium dynamics of the associated decomposed energy. Our results show the appropriateness of LR on finer molecular scales.

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.

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.

Combining non-equilibrium simulations and coarse-grained modelling allows for a fine-grained decomposition of solvation dynamics

The time-dependent Stokes shift is shown to be a localized and short-ranged effect in ionic liquids.

Dielectric spectra of ionic liquids and their conversion to solvation dynamics: a detailed computational analysis of polarizable systems

For the three molecular ionic liquids 1-ethyl-3-methylimidazolium tetrafluoroborate, 1-ethyl-3-methylimidazolium trifluoromethanesulfonate and 1-butyl-3-methylimidazolium tetrafluoroborate, dielectric spectra were calculated from molecular dynamics simulations based on polarizable force fields. Using the reaction field continuum model the dielectric spectra were converted to the solvation dynamics of coumarin 153. It is shown in detail that the inclusion of the static conductivity in this model is essential. When simplifying the dielectric spectrum to the static conductivity hyperbola, the solvation response function becomes mono-exponential. Taking into account the frequency dependence of the conductivity, the typical two time-regimes of the solvation response function in ionic liquids are already obtained. However, the mean relaxation time remains the same. When converting the complete dielectric spectrum, i.e. also including frequency-dependent dielectric permittivity, quantitative changes are observed, but the qualitative shape is conserved. In accordance with previous experimental studies, solvation dynamics in ionic liquids predicted by the reaction field continuum model is too fast for longer times. This correlates with the suppression of the fine structure of the dielectric spectrum at low frequencies by the static conductivity hyperbola. By scaling down the static conductivity this effect can be partially amended. In addition to the impact of the solvent dielectric spectrum on solvation dynamics, solute-specific effects, i.e. anisotropy in shape and charge distribution as well as polarizability, were also studied.

Observations of probe dependence of the solvation dynamics in ionic liquids

Solvation and rotational dynamics of 4-aminophthalimide (4AP) in four ionic liquids (ILs) are measured using a combination of fluorescence upconversion spectroscopy and time-correlated single photon counting.

Ultrafast solvation dynamics and charge transfer reactions in room temperature ionic liquids

In this perspective, we review the recent studies concerning the liquid structure and solvation dynamics of ionic liquids, and their influence upon electron transfer and charge transfer reactions.

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.

Solvent and rotational relaxation of coumarin-153 and coumarin-480 in ionic liquid (1-butyl-3-methylimidazolium tetrafluoroborate) modified sodium 1,4-bis(2-ethylhexyl) sulfosuccinate (NaAOT) micelle

Understanding ion transport dynamics, structure of surfactant aggregates in ionic liquids or ionic liquid/water solutions are quite interesting and potentially important due to widespread applications of surfactant-based systems. In this manuscript we have investigated the effect of 1-butyl-3-methylimidazolium tetrafluoroborate (bmimBF(4)) addition on solvent and rotational relaxation of coumarin-153 (C-153) and coumarin-480 (C-480) in aqueous solution of sodium 1,4-bis(2-ethylhexyl) sulfosuccinate (NaAOT) using steady state and picosecond time resolved fluorescence spectroscopy. The strong adsorption of the bmim(+) at the interface and the role of the ionic liquid particularly the cation bmim(+) in the modification of the interfacial geometry were probed by the analysis of decay parameters and the rotational relaxation parameters. Since the addition of the NaAOT in water-bmimBF(4) mixture above critical micellar concentration (48 mM, obtained from observing pyrene fluorescence) causes strong adsorption of the ionic liquid particularly the cation bmim(+), the average solvation time, particularly the slow component increases significantly. More importantly we have found the probe dependent solvation dynamics due to the different location of the probe molecules, C-153 and C-480. C-153 being hydrophobic in nature resides in the stern layer and the adsorption of the bmim(+) at the interface modifies stern layer more effectively. So we have observed more pronounced change in solvation dynamics in case of C-153 compared to that in case of C-480. The fluorescence anisotropy decays of the probe molecules were found to be biexponential in nature. The anisotropy decay was interpreted by using a model which consists of the wobbling (rotational) and translational diffusion of the dye coupled with the rotational motion of the micelle as a whole.