Characterization of the solvation dynamics of an ionic liquid via molecular dynamics simulation

Solvation dynamics on the diffusion timescale elucidated using energy-represented dynamics theory

Photoexcitation of a solute alters the solute-solvent interaction, resulting in the nonequilibrium relaxation of the solvation structure, often called a dynamic Stokes shift or solvation dynamics. Thanks to the local nature of the solute-solvent interaction, the characteristics of the local solvent environment dissolving the solute can be captured by the observation of this process. Recently, we derived the energy-represented Smoluchowski-Vlasov (ERSV) equation, a diffusion equation for molecular liquids, which can be used to analyze the solvation dynamics on the diffusion timescale. This equation expresses the time development for the solvent distribution on the solute-solvent pair interaction energy (energy coordinate). Since the energy coordinate can effectively treat the solvent flexibility in addition to the position and orientation, the ERSV equation can be utilized in various solvent systems. Here, we apply the ERSV equation to the solvation dynamics of 6-propionyl-2-dimethylamino naphthalene (Prodan) in water and different alcohol solvents (methanol, ethanol, and 1-propanol) for clarifying the differences of the relaxation processes among these solvents. Prodan is a solvent-sensitive fluorescent probe and is thus widely utilized for investigating heterogeneous environments. On the long timescale, the ERSV equation satisfactorily reproduces the relaxation time correlation functions obtained from the molecular dynamics (MD) simulations for these solvents. We reveal that the relaxation time coefficient on the diffusion timescale linearly correlates with the inverse of the translational diffusion coefficients for the alcohol solvents because of the Prodan-solvent energy distributions among the alcohols. In the case of water, the time coefficient deviates from the linear relationship for the alcohols due to the difference in the extent of importance of the collective motion between the water and alcohol solvents.

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.

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.

Understanding Structure-Property Correlation in Monocationic and Dicationic Ionic Liquids through Combined Fluorescence and Pulsed-Field Gradient (PFG) and Relaxation NMR Experiments

Steady state, time-resolved fluorescence and NMR experiments are carried out to gain deeper insights into the structure-property correlation in structurally similar monocationic and dicationic room-temperature ionic liquids (RTILs). The excitation wavelength dependent fluorescence response of fluorophore in 1-methy-3-propyllimidazolium bis(trifluoromethylsulfonyl)amide [C3MIm][NTf2] is found to be different from that of 1,6-bis(3-methylimidazolium-1-yl)hexane bis(trifluoromethylsulfonyl)amide [C6(MIm)2][NTf2]2 and 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide [C6MIm][NTf2]. The outcomes of the present solvent dynamics study in [C3MIm][NTf2] when compared with those in [C6(MIm)2][NTf2]2 and in [C6MIm][NTf2] from our previous studies (Phys. Chem. Chem. Phys. 2014, 16, 12918-12928) indicate the involvement of dipolar rotation of imidazolium cation during solvation. To correlate the findings of solvation dynamics study with the dipolar rotation of the imidazolium ring, pulsed-field gradient (PFG)-NMR technique for translational diffusion coefficient measurement and (1)H as well as (19)F spin-lattice relaxation measurements are employed. NMR investigation reveals that an ultrafast component of solvation can be related to the dipolar rotation of imidazolium cation; hence, the role of dipolar rotation of cations in governing the dynamics of solvation in ILs cannot be ignored. Analysis of the rotational relaxation dynamics data by the Stokes-Einstein-Debye hydrodynamic theory unveils distinctive features of solute-solvent interaction in [C3MIm][NTf2] and [C6(MIm)2][NTf2]2.

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.

Molecular dynamics investigation of the vibrational spectroscopy of isolated water in an ionic liquid

Experimental studies examining the structure and dynamics of water in ionic liquids (ILs) have revealed local ion rearrangements that occur an order of magnitude faster than complete randomization of the liquid structure. Simulations of an isolated water molecule embedded in 1-butyl-3-methyl imidazolium hexafluorophosphate, [bmim][PF6], were performed to shed insight into the nature of these coupled water-ion dynamics. The theoretical calculations of the spectral diffusion dynamics and the infrared absorption spectra of the OD stretch of isolated HOD in [bmim][PF6] agree well with experiment. The infrared absorption line shape of the OD stretch is narrower and blue-shifted in the IL compared to those in aqueous solution. Decomposition of the OD frequency time correlation function revealed that translational motions of the anions dominate the spectral diffusion dynamics.

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.

Fluorescence response of a dipolar organic solute in a dicationic ionic liquid (IL): is the behavior of dicationic IL different from that of usual monocationic IL?

The solvation dynamics of coumarin 153 are compared in di and in mono-cationic-based ionic liquids for the first time.

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.

Aggregation behavior of Triton X-100 with a mixture of two room-temperature ionic liquids: can we identify the mutual penetration of ionic liquids in ionic liquid containing micellar aggregates?

In this manuscript, we have characterized two different micellar aggregates containing all nonvolatile components. We have shown (i) the effect of ethylammonium nitrate (EAN) addition on the properties of micellar solution of Triton X-100 in 1-butyl-3-methylimidazolium hexafluorophosphate (bmimPF(6)) and (ii) the effect of bmimPF(6) addition on the properties of micellar solution of Triton X-100 in EAN. To investigate the effect, we have used (1)H NMR, pulsed-field gradient spin-echo NMR (PFGSE NMR), and methyl orange (MO) and coumarin 153 (C-153) as absorption and emission probes, respectively. The penetration of added EAN inside the Triton X-100/bmimPF(6) micellar aggregates is indicated by (i) red shift in both the absorption spectra of MO and emission spectra of C-153 and (ii) downfield shift of proton signals of ethylene oxide units in Triton X-100. On the other hand, (1)H NMR and PFGSE NMR indicates the penetration of added bmimPF(6) inside the Triton X-100/EAN micellar aggregates. However, the constancy of both the absorption spectra of MO and emission spectra of C-153 indicates that the microenvironment around the probe molecules remains unaffected. We have also investigated the effect of micelle formation and the effect of penetration of ionic liquids (ILs) in micellar aggregates, on the solvation dynamics of C-153. The solvent relaxation around C-153 gets retarded on going from neat ILs to the micellar solution of Triton X-100 in ILs. In addition to this, we have also observed that with the addition of EAN in Triton X-100/bmimPF(6) micellar aggregates the solvation dynamics becomes faster, whereas with the addition of bmimPF(6) in Triton X-100/EAN micellar aggregates we did not observe any notable change in solvation dynamics. This observation further supports the conclusions drawn from UV-visible and NMR studies.

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.