Dynamics of a supercooled ionic liquid studied by optical and dielectric spectroscopy

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming ionic liquid dibutylammonium formate near its glass transition temperature Tg=153 K. The Stokes-shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical line width changes systematically along the solvation coordinate. The solvent dynamics in the viscous regime is compared with the rotational behavior of the solute and with the dielectric relaxation of the ionic liquid. Among the different quantities derived from the dielectric experiments, the relaxation of the macroscopic electric field, i.e., the modulus Mt, matches best the solvent response Ct regarding time scale and stretching exponent. Many other properties are reminiscent of the behavior of polar molecular liquids which lack the ionic character.     

Dipolar solvation dynamics in room temperature ionic liquids: an effective medium calculation using dielectric relaxation data

Solvation dynamics in four imidazolium cation based room temperature ionic liquids (RTIL) have been calculated by using the recently measured dielectric relaxation data [ J. Phys. Chem. B 2008, 112, 4854 ] as an input in a molecular hydrodynamic theory developed earlier for studying solvation energy relaxation in polar solvents. Coumarin 153 (C153), 4-aminophthalimide (4-AP), and trans-4-dimethylamino-4'-cyanostilbene (DCS) have been used as probe molecules for this purpose. The medium response to a laser-excited probe molecule in an ionic liquid is approximated by that in an effective dipolar medium. The calculated decays of the solvent response function for these RTILs have been found to be biphasic and the decay time constants agree well with the available experimental and computer simulation results. Also, no probe dependence has been found for the average solvation times in these ionic liquids. In addition, dipolar solvation dynamics have been predicted for two other RTILs for which experimental results are not available yet. These predictions should be tested against experiments and/or simulation studies.     

Kinetic study of thermal Z to E isomerization reactions of azobenzene and 4-dimethylamino-4'-nitroazobenzene in ionic liquids [1-R-3-methylimidazolium bis(trifluoromethylsulfonyl)imide with R = butyl, pentyl, and hexyl

Thermal Z to E isomerization reactions of azobenzene and 4-dimethylamino-4'-nitroazobenzene were examined in three ionic liquids of general formula 1-R-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (R = butyl, pentyl, and hexyl). The first-order rate constants and activation energies for the reactions of azobenzene measured in these ionic liquids were consistent with those measured in ordinary organic solvents, which indicated that the slow isomerization through the inversion mechanism with a nonpolar transition state was little influenced by the solvent properties, such as the viscosity and dielectric constant, of ionic liquids. On the other hand, the rate constants and the corresponding frequency factors of the Arrhenius plot were significantly reduced for the isomerization of 4-dimethylamino-4'-nitroazobenzene in ionic liquids compared with those for the isomerization in ordinary organic molecular solvents with similar dielectric properties. Although these ionic liquids are viscous, the apparent viscosity dependence of the rate constant could not be explained either by the Kramers-Grote-Hynes model or by the Agmon-Hopfield model for solution reactions. It is proposed that the positive and the negative charge centers of a highly polar rotational transition state are stabilized by the surrounding anions and cations, respectively, and that the ions must be rearranged so as to form highly ordered solvation shells around the charge centers of the reactant in the transition state. This requirement for the orderly solvation in the transition state results in unusually small frequency factors of 10(4)-10(7) s(-1).     

Solvation dynamics and electric field relaxation in an imidazolium-PF6 ionic liquid: from room temperature to the glass transition

Time-resolved phosphorescence spectra and anisotropy of quinoxaline were measured in an ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate (BMIM-HFP), in its supercooled state near the glass-transition temperature. The solvation dynamics results are compared with the rotational motion of the probe and with the dielectric behavior of the neat ionic liquid. The dynamics in the viscous state are highly dispersive and show a super-Arrhenius temperature dependence, as typical for glass-forming materials. Combined with room-temperature results, solvation dynamics is observed to follow the structural relaxation times in terms of eta/T for more than 10 decades, from subnanoseconds at room temperature to seconds near the glass-transition temperature T(g). The dielectric modulus relaxation follows this trend only for temperatures T > 1.2T(g) and departs significantly from eta/T in the 1.1T(g) > T > T(g) range. This deviation is reminiscent of the enhanced translational diffusion or fractional Stokes-Einstein behavior observed in many fragile supercooled liquids. Because the electric field relaxation in BMIM-HFP includes dc conductivity, this correlation function involves translational motion and thus displays the effect of enhanced diffusivity. A microscopic model is required for rationalizing the decoupling of solvation dynamics from the longitudinal time scales and the limitation of this effect to the viscous regime with T < 1.2T(g).     

Density scaling of the transport properties of molecular and ionic liquids

Casalini and Roland [Phys. Rev. E 69, 062501 (2004); J. Non-Cryst. Solids 353, 3936 (2007)] and other authors have found that both the dielectric relaxation times and the viscosity, η, of liquids can be expressed solely as functions of the group (TV?(γ)), where T is the temperature, V is the molar volume, and γ a state-independent scaling exponent. Here we report scaling exponents γ, for the viscosities of 46 compounds, including 11 ionic liquids. A generalization of this thermodynamic scaling to other transport properties, namely, the self-diffusion coefficients for ionic and molecular liquids and the electrical conductivity for ionic liquids is examined. Scaling exponents, γ, for the electrical conductivities of six ionic liquids for which viscosity data are available, are found to be quite close to those obtained from viscosities. Using the scaling exponents obtained from viscosities it was possible to correlate molar conductivity over broad ranges of temperature and pressure. However, application of the same procedures to the self-diffusion coefficients, D, of six ionic and 13 molecular liquids leads to superpositioning of poorer quality, as the scaling yields different exponents from those obtained with viscosities and, in the case of the ionic liquids, slightly different values for the anion and the cation. This situation can be improved by using the ratio (D∕T), consistent with the Stokes-Einstein relation, yielding γ values closer to those of viscosity.     

Sound velocity dispersion in room temperature ionic liquids studied using the transient grating method

Sound velocity is determined by the transient grating method in a range from 10(6) to 10(10) Hz in three room temperature ionic liquids, 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide, 1-butyl-3-methylimidazolium hexafluorophosphate, and N,N,N-trimethyl-N-propylammonium bis(trifluoromethanesulfonyl)imide. In all room temperature ionic liquids studied, the sound velocity increased with increasing frequency. The cause of this change is posited to be structural relaxation in the room temperature ionic liquids. Frequency dependence of the sound velocity is not reproduced by a simple Debye relaxation model. The sound velocity dispersion relation in 1-butyl-3-methylimidazolium hexafluorophosphate matches a Cole-Davidson function with parameters determined by a dielectric relaxation [C. Daguenet et al., J. Phys. Chem. B 110, 12682 (2006)], indicating that structural and reorientational relaxations are strongly coupled. Conversely, the sound velocity dispersions of the other two ionic liquids measured do not match those measured for dielectric relaxation, implying that structural relaxation is much faster than the reorientational relaxation. This difference is discussed in relation to the motilities of anions and cations.     

Fluorescence studies in a pyrrolidinium ionic liquid: polarity of the medium and solvation dynamics

While the imidazolium ionic liquids have been studied for some time, little is known about the pyrrolidinium ionic liquids. In this work, steady-state and picosecond time-resolved fluorescence behavior of three electron donor-acceptor molecules, coumarin-153 (C153), 4-aminophthalimide (AP), and 6-propionyl-2-dimethylaminonaphthalene (PRODAN), has been studied in a pyrrolidinium ionic liquid, N-butyl-N-methylpyrrolidinium bis(trifluoromethanesulfonyl)imide, abbreviated here as [bmpy][Tf2N]. The steady-state fluorescence data of the systems suggest that the microenvironment around these probe molecules, which is measured in terms of the solvent polarity parameter, E(T)(30), is similar to that in 1-decanol and that the polarity of this ionic liquid is comparable to that of the imidazolium ionic liquids. All three systems exhibit wavelength-dependent fluorescence decay behavior, and the time-resolved fluorescence spectra show a progressive shift of the fluorescence maximum toward the longer wavelength with time. This behavior is attributed to solvent-mediated relaxation of the fluorescent state of these systems. The dynamics of solvation, which is studied from the time-dependent shift of the fluorescence spectra, suggests that approximately 45% of the relaxation is too rapid to be measured in the present setup having a time resolution of 25 ps. The remaining observable components of the dynamics consist of a short component of 115-440 ps (with smaller amplitude) and a long component of 610-1395 ps (with higher amplitude). The average solvation time is consistent with the viscosity of this ionic liquid. The dynamics of solvation is dependent on the probe molecule, and nearly 2-fold variation of the solvation time depending on the probe molecule could be observed. No correlation of the solvation time with the probe molecule could, however, be observed.     

Solvation and rotational dynamics of coumarin 153 in ionic liquids: comparisons to conventional solvents

Steady-state and time-resolved emission spectroscopy with 25 ps resolution are used to measure equilibrium and dynamic aspects of the solvation of coumarin 153 (C153) in a diverse collection of 21 room-temperature ionic liquids. The ionic liquids studied here include several phosphonium and imidazolium liquids previously reported as well as 12 new ionic liquids that incorporate two homologous series of ammonium and pyrrolidinium cations. Steady-state absorption and emission spectra are used to extract solvation free energies and reorganization energies associated with the S0 <--> S1 transition of C153. These quantities, especially the solvation free energy, vary relatively little in ionic liquids compared to conventional solvents. Some correlation is found between these quantities and the mean separation between ions (or molar volume). Time-resolved anisotropies are used to observe solute rotation. Rotation times measured in ionic liquids correlate with solvent viscosity in much the same way that they do in conventional polar solvents. No special frictional coupling between the C153 and the ionic liquid solvents is indicated by these times. But, in contrast to what is observed in most low-viscosity conventional solvents, rotational correlation functions in ionic liquids are nonexponential. Time-resolved Stokes shift measurements are used to characterize solvation dynamics. The solvation response functions in ionic liquids are also nonexponential and can be reasonably represented by stretched-exponential functions of time. The solvation times observed are correlated with the solvent viscosity, and the much slower solvation in ionic liquids compared to dipolar solvents can be attributed to their much larger viscosities. Solvation times of the majority of ionic liquids studied appear to follow a single correlation with solvent viscosity. Only liquids incorporating the largest phosphonium cation appear to follow a distinctly different correlation.     

Chromatographic and spectroscopic methods for the determination of solvent properties of room temperature ionic liquids

Room temperature ionic liquids are novel solvents with favorable environmental and technical features. Synthetic routes to over 200 room temperature ionic liquids are known but for most ionic liquids physicochemical data are generally lacking or incomplete. Chromatographic and spectroscopic methods afford suitable tools for the study of solvation properties under conditions that approximate infinite dilution. Gas-liquid chromatography is suitable for the determination of gas-liquid partition coefficients and activity coefficients as well as thermodynamic constants derived from either of these parameters and their variation with temperature. The solvation parameter model can be used to define the contribution from individual intermolecular interactions to the gas-liquid partition coefficient. Application of chemometric procedures to a large database of system constants for ionic liquids indicates their unique solvent properties: low cohesion for ionic liquids with weakly associated ions compared with non-ionic liquids of similar polarity; greater hydrogen-bond basicity than typical polar non-ionic solvents; and a range of dipolarity/polarizability that encompasses the same range as occupied by the most polar non-ionic liquids. These properties can be crudely related to ion structures but further work is required to develop a comprehensive approach for the design of ionic liquids for specific applications. Data for liquid-liquid partition coefficients is scarce by comparison with gas-liquid partition coefficients. Preliminary studies indicate the possibility of using the solvation parameter model for interpretation of liquid-liquid partition coefficients determined by shake-flask procedures as well as the feasibility of using liquid-liquid chromatography for the convenient and rapid determination of liquid-liquid partition coefficients. Spectroscopic measurements of solvatochromic and fluorescent probe molecules in room temperature ionic liquids provide insights into solvent intermolecular interactions although interpretation of the different and generally uncorrelated "polarity" scales is sometimes ambiguous. All evidence points to the ionic liquids as a unique class of polar solvents suitable for technical development. In terms of designer solvents, however, further work is needed to fill the gaps in our knowledge of the relationship between ion structures and physicochemical properties.     

Dielectric and shear mechanical alpha and beta relaxations in seven glass-forming liquids

We present shear mechanical and dielectric measurements taken on seven liquids: triphenylethylene, tetramethyltetra-phenyltrisiloxane (Dow Corning 704 diffusion pump fluid), polyphenyl ether (Santovac 5 vacuum pump fluid), perhydrosqualene, polybutadiene, decahydroisoquinoline (DHIQ), and tripropylene glycol. The shear mechanical and dielectric measurements are for each liquid performed under identical thermal conditions close to the glass transition temperature. The liquids span four orders of magnitude in dielectric relaxation strength and include liquids with and without Johari-Goldstein beta relaxation. The shear mechanical data are obtained by the piezoelectric shear modulus gauge method giving a large frequency span (10(-3)-10(4.5) Hz). This allows us to resolve the shear mechanical Johari-Goldstein beta peak in the equilibrium DHIQ liquid. We moreover report a signature (a pronounced rise in the shear mechanical loss at frequencies above the alpha relaxation) of a Johari-Goldstein beta relaxation in the shear mechanical spectra for all the liquids which show a beta relaxation in the dielectric spectrum. It is found that both the alpha and beta loss peaks are shifted to higher frequencies in the shear mechanical spectrum compared to the dielectric spectrum. It is in both the shear and dielectric responses found that liquids obeying time-temperature superposition also have a high-frequency power law with exponent close to -12. It is moreover seen that the less temperature dependent the spectral shape is, the closer it is to the universal -12 power-law behavior. The deviation from this universal power-law behavior and the temperature dependencies of the spectral shape are rationalized as coming from interactions between the alpha and beta relaxations.     

Why are ionic liquids liquid? A simple explanation based on lattice and solvation energies

We have developed a simple and quantitative explanation for the relatively low melting temperatures of ionic liquids (ILs). The basic concept was to assess the Gibbs free energy of fusion (Delta(fus)G) for the process IL(s) --> IL(l), which relates to the melting point of the IL. This was done using a suitable Born-Fajans-Haber cycle that was closed by the lattice (i.e., IL(s) --> IL(g)) Gibbs energy and the solvation (i.e., IL(g) --> IL(l)) Gibbs energies of the constituent ions in the molten salt. As part of this project we synthesized and determined accurate melting points (by DSC) and dielectric constants (by dielectric spectroscopy) for 14 ionic liquids based on four common anions and nine common cations. Lattice free energies (Delta(latt)G) were estimated using a combination of Volume Based Thermodynamics (VBT) and quantum chemical calculations. Free energies of solvation (Delta(solv)G) of each ion in the bulk molten salt were calculated using the COSMO solvation model and the experimental dielectric constants. Under standard ambient conditions (298.15 K and 10(5) Pa) Delta(fus)G degrees was found to be negative for all the ILs studied, as expected for liquid samples. Thus, these ILs are liquid under standard ambient conditions because the liquid state is thermodynamically favorable, due to the large size and conformational flexibility of the ions involved, which leads to small lattice enthalpies and large entropy changes that favor melting. This model can be used to predict the melting temperatures and dielectric constants of ILs with good accuracy. A comparison of the predicted vs experimental melting points for nine of the ILs (excluding those where no melting transition was observed and two outliers that were not well described by the model) gave a standard error of the estimate (s(est)) of 8 degrees C. A similar comparison for dielectric constant predictions gave s(est) as 2.5 units. Thus, from very little experimental and computational data it is possible to predict fundamental properties such as melting points and dielectric constants of ionic liquids.     

Fluorescence probing of temperature-dependent dynamics and friction in ionic liquid local environments

The solvation dynamics and local orientational friction for a series of four ionic liquids have been probed using coumarin 153 (C153) as a function of temperature. These ionic liquids are comprised of nonaromatic organic cations paired with a common anion, bis(trifluoromethylsulfonyl)imide (NTf(2)-). The specific liquids are as follows: N-methyl-tri-N-butylammonium NTf(2)- (N(1444)+/NTf(2)-), N-hexyl-tri-N-butylammonium NTf(2)- (N(6444)+/NTf(2)-), N-methyl-N-butylpyrrolidinium NTf(2)- (Pyrr(14)+/NTf(2)-), and N-methyl-N-ethoxyethylpyrrolidinium NTf(2)- (Pyrr(1(2O2))+/NTf(2)-). The observed solvation dynamics and fluorescence depolarization dynamics occur over a broad range of time scales that can only be adequately fit by functions including three or more exponential components. Stretched exponential distributions cannot adequately fit our data. The solvation and reorientational dynamics of the C153 probe are studied over a range of temperatures from 278.2 to 353.2 K. For both the solvation dynamics and the probe reorientational dynamics, the observed temperature dependence is well fit by a Vogel-Tammann-Fulcher law. To correlate the observed microscopic dynamics with macroscopic physical properties, temperature-dependent viscosities are also measured. Differential scanning calorimetry is used to study the thermodynamics of the phase transitions from the liquid to supercooled liquid to glassy states. For the two tetraalkylammonium liquids, the observed melting transitions occur near 300 K, so we are able to study the dynamics in a clearly supercooled regime. Very long time scale orientational relaxation time constants dynamics on the order of 100 ns are observed in the C153 fluorescence anisotropy. These are interpreted to arise from long-lived local structures in the environment surrounding the C153 probe.     

Effects of ion atmosphere relaxation on dipole isomerization reactions in electrolyte solutions

A theoretical study of the effects of ion atmosphere relaxation on the rate of a model dipole isomerization reaction in electrolyte solutions is presented. The time-dependent ion atmosphere friction is calculated by using a molecular hydrodynamic theory which properly includes the static and dynamic interionic correlations through ionic structure factors and van Hove functions. The rate constant is determined by employing the well-known Grote-Hynes theory. Numerical results are obtained for the time-dependent ion atmosphere friction and for the rate of isomerization reaction in electrolyte solutions of varying ion concentration and dielectric constant. It is found that the ion atmosphere friction can have significant effects in reducing the rate of isomerization below the prediction of equilibrium solvation transition state theory.     

Solvation studies of a zinc finger protein in hydrated ionic liquids

The solvation of the zinc finger protein with the PDB-ID “5ZNF” in hydrated ionic liquids was studied at varying water content. 1-Ethyl-3-methylimidazolium and trifluoromethanesulfonate were the cation and anion, respectively. The protein stability as well as the solvation structure, the shell dynamics and the shell resolved dielectric properties were investigated by means of molecular dynamics simulations. The lengths of the respective trajectories extended up to 200 nanoseconds in order to cover the complete solvent dynamics. Considering the above mentioned properties as a function of the water content they all exhibit a maximum or minimum at the very same mole fraction. While the exact value x(H(2)O) = 0.927 depends on the underlying force field, its origin may be traced back to the competition between the van der Waals and the electrostatic energy of the protein as well as to the transition from aqueous dielectric screening to ionic charge screening with decreasing water content. The parameter-free Voronoi decomposition of space served as a basis for the analysis of most results. In particular, solvation shells were naturally inferred from this concept. In addition to the molecular analysis a mesoscopic view is given in terms of dielectric properties. Thereby, the net dielectric constant is decomposed into contributions from the protein, the first and second solvation shells as well as the bulk. Cross-terms between these components are given, too.     

Orientational and translational dynamics in room temperature ionic liquids

The authors investigate the dynamics of a series of room temperature ionic liquids, based on the same 1-butyl-3-methylimidazolium cation with different anions, by means of broadband (10(-6)-10(9) Hz) dielectric spectroscopy and depolarized light scattering in the temperature range from 400 K down to 35 K. Typical ionic conductivity is observed above the glass transition temperature Tg. Below Tg the authors detect relaxation processes that exhibit characteristics of secondary relaxations, as typically observed in molecular glasses. At high temperatures, the characteristic times of cation reorientation, deduced from the light scattering data, are approximately equal to the electric modulus relaxation times related to ionic conductivity. In the supercooled regime and close to Tg, the authors observe decoupling of conductivity from structural relaxation. Overall, room temperature ionic liquids exhibit typical glass transition dynamics, apparently unaltered by Coulomb interactions.     

Heterogeneous thermal excitation and relaxation in supercooled liquids

We investigate a phenomenological model which rationalizes the effects of dielectric hole burning on the basis of heterogeneous dielectric and specific heat relaxation in supercooled liquids. The quantitative agreement between model predictions and dielectric hole-burning observations is lost if the assumption of correlated dielectric and thermal relaxation times is removed from the model. This suggests that dynamically distinct domains in real liquids are associated with a time constant which characterizes both the structural and thermal relaxation behaviors. The calculations demonstrate that the observed burn-induced modifications reflect the spectral selectivity and persistence time of the fictive temperatures within these domains, and that 100 or more cycles of the sinusoidal burn field can be required to saturate the heat accumulated in the slow degrees of freedom. It is also shown that the recovery of dielectric holes is entirely accounted for by the model, and that the persistence times do not provide direct insight into rate exchange processes. Additionally, the model predicts that the heating effects considered here are a significant source of nonlinear dielectric behavior, even in the absence of deliberate frequency selective hole burning.     

Unravelling free volume in branched-cation ionic liquids based on silicon

The branching of ionic liquid cation sidechains utilizing silicon as the backbone was explored and it was found that this structural feature leads to fluids with remarkably low density and viscosity. The relatively low liquid densities suggest a large free volume in these liquids. Argon solubility was measured using a precise saturation method to probe the relative free volumes. Argon molar solubilities were slightly higher in ionic liquids with alkylsilane and siloxane groups within the cation, compared to carbon-based branched groups. The anion size, however, showed by far the dominant effect on argon solubility. Thermodynamic solvation parameters were derived from the solubility data and the argon solvation environment was modelled utilizing the polarizable CL&Pol force field. Semiquantitative analysis was in agreement with trends established from the experimental data. The results of this investigation demonstrate design principles for targeted ionic liquids when optimisation for the free volume is required, and demonstrate the utility of argon as a simple, noninteracting probe. As more ionic liquids find their way into industrial processes of scale, these findings are important for their utilisation in the capture of any gaseous solute, gas separation, or in processes involving the transformation of gases or small molecules.

The branching of ionic liquid cation sidechains utilizing silicon as the backbone was explored and it was found that this structural feature leads to fluids with remarkably low density and viscosity.     

On the influence of hydrated ionic liquids on the dynamical structure of model proteins: a computational study

The solvation of the protein ubiquitin (PDB entry "1UBQ") in hydrated molecular ionic liquids was studied for varying water content or, equivalently, a diversity of ionic strengths. The cations and anions were 1-ethyl-3-methylimidazolium and trifluoromethanesulfonate, respectively. The protein's shape and stability as well as the solvation structure, the shell dynamics and the shell resolved dielectric properties were investigated by means of molecular dynamics simulations. The respective simulation trajectories covered 200 nanoseconds. Besides the characteristic point already found for the zinc finger motif at the transition from the pure aqueous environment to the ionic solution an even more pronounced state is found where several properties show extremal behaviour (maximum or minimum). This second characteristic point occurs at the transition from the ionic solution to the hydrated ionic melt where water changes its role from a solvent to a co-solvent. Most of the data analysis presented here is based on the Voronoi decomposition of space.     

A molecular Debye-Hückel theory and its applications to electrolyte solutions

In this report, a molecular Debye-Hu?ckel theory for ionic fluids is developed. Starting from the macroscopic Maxwell equations for bulk systems, the dispersion relation leads to a generalized Debye-Hu?ckel theory which is related to the dressed ion theory in the static case. Due to the multi-pole structure of dielectric function of ionic fluids, the electric potential around a single ion has a multi-Yukawa form. Given the dielectric function, the multi-Yukawa potential can be determined from our molecular Debye-Hu?ckel theory, hence, the electrostatic contributions to thermodynamic properties of ionic fluids can be obtained. Applications to binary as well as multi-component primitive models of electrolyte solutions demonstrated the accuracy of our approach. More importantly, for electrolyte solution models with soft short-ranged interactions, it is shown that the traditional perturbation theory can be extended to ionic fluids successfully just as the perturbation theory has been successfully used for short-ranged systems.     

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

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.