Dynamics of solvation of an ion in a dense dipolar liquid

A molecular theory of the dynamics of solvation of an ion in a dense dipolar liquid is presented. The theory is based on an extended hydrodynamic approach that properly includes the interparticle correlations that are present at molecular length scales. The effects of the solvent inertial and viscoelastic responses are also included consistently. Numerical studies reveal rich relaxation behaviour such as short-time oscillations followed by a slow long-time decay. The results are in semi-quantitative agreement with recent computer simulation studies.     

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.     

Dynamics of solvent and rotational relaxation of Coumarin 153 in a room temperature ionic liquid, 1-butyl-3-methylimidazolium octyl sulfate, forming micellar structure

The solvent and rotational relaxation of Coumarin 153 (C-153) was investigated by picosecond time-resolved fluorescence spectroscopy in a room temperature ionic liquid (RTIL), 1-butyl-3-methylimidazolium octyl sulfate ([C4mim][C8SO4]). This is a typical RTIL, which form micellar structure above certain concentration of the RTIL (0.031 M). Dynamic light scattering (DLS) measurements show that the average hydrodynamic diameter ( Dh) of a [C4mim][C8SO4]-water micelle is 2.8 (+/-0.2) nm. Both the solvent and rotational relaxation of C-153 are retarded in this micelle compared to the solvation time of a similar type of dye in neat water. However, the solvent relaxation in this ionic liquid surfactant is different from that of a conventional ionic surfactant. The slow component of the solvation dynamics in C8H17SO4Na or TX-100 micelle is on the nanoseconds time scale, whereas in [C4mim][C8SO4] micelle the same component is on the subnanoseconds time scale. The different molecular motions with different time scale is the main reason behind this difference in the solvation time in micelles composed of RTIL with other conventional micelles.     

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

The solvation dynamics of ionic liquids have been the subject of intense experimental study but remain poorly understood. We present the results of molecular dynamics simulations of the solvation dynamics of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate in response to photoexcitation of the fluorescent dye coumarin-153. We reproduce the time-resolved fluorescence Stokes shift using linear response theory, then use novel statistical techniques to analyze cation and anion contributions to the signal. We find that the solvation dynamics are dominated by collective ionic motion and characterize the time scale for various features of the collective response. Further, we use the Steele analysis [Mol. Phys. 61, 1031 (1987)] to characterize the contributions to the observed Stokes shift made by translational and rovibrational degrees of freedom. Our results indicate that in contrast to molecular liquids, the rovibrational response is trivial and the observed fluorescence response arises almost entirely from ionic translation. Our results resolve previously open questions in the literature about the nature of the rapid dynamics in room-temperature ionic liquids and offer insight into the physical principles governing ionic liquid behavior on longer time scales.     

Solute rotation and solvation dynamics in an alcohol-functionalized room temperature ionic liquid

Steady-state and time-resolved fluorescence behaviors of two dipolar solutes, coumarin 153 and 4-aminophthalimide, have been studied in an alcohol-functionalized room-temperature ionic liquid, 1-(hydroxyethyl)-3-methylimidazolium bis(trifluoromethanesulfonyl)imide. The steady-state fluorescence parameters have been exploited for the estimation of the polarity of this ionic liquid and to obtain information on the hydrogen bonding interaction between the ionic liquid and the probe molecules. The time-resolved measurements have been focused on the dynamics of solvation by studying the dynamic Stokes shift in the ps-ns time scale and solute rotation by measuring the time dependence of the fluorescence anisotropy. The time-resolved anisotropy studies reveal a significant slow down of the rotational motion of one of the probe molecules. The time-dependent fluorescence Stokes shift measurements suggest that the time-resolvable part of the dynamics is biphasic in nature, highly dependent on the probe molecule and the ultrafast component is comparatively less than that in other ionic liquids. The influence of the hydrogen bonding interaction between the probe molecules and the ionic liquids on the solute rotation and the various components of the solvation dynamics is carefully analyzed in an attempt to obtain further insight into the mechanism of solvation in these novel media.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

New approach to free energy of solvation applying continuum models to molecular dynamics simulation

A new approach to the calculation of the free energy of solvation from trajectories obtained by molecular dynamics simulation is presented. The free energy of solvation is computed as the sum of three contributions originated at the cavitation of the solute by the solvent, the solute-solvent nonpolar (repulsion and dispersion) interactions, and the electrostatic solvation of the solute. The electrostatic term is calculated based on ideas developed for the broadly used continuum models, the cavitational contribution from the excluded volume by the Claverie-Pierotti model, and the Van der Waals term directly from the molecular dynamics simulation. The proposed model is tested for diluted aqueous solutions of simple molecules containing a variety of chemically important functions: methanol, methylamine, water, methanethiol, and dichloromethane. These solutions were treated by molecular dynamics simulations using SPC/E water and the OPLS force field for the organic molecules. Obtained free energies of solvation are in very good agreement with experimental data.     

A noniterative mean-field QM/MM-type approach with a linear response approximation toward an efficient free-energy evaluation

Mean-field treatment of solvent provides an efficient technique to investigate chemical processes in solution in quantum mechanics/molecular mechanics (QM/MM) framework. In the algorithm, an iterative calculation is required to obtain the self-consistency between QM and MM regions, which is a time-consuming step. In the present study, we have proposed a noniterative approach by introducing a linear response approximation (LRA) into the solvation term in the one-electron part of Fock matrix in a hybrid approach between molecular-orbital calculations and a three-dimensional (3D) integral equation theory for molecular liquids (multicenter molecular Ornstein–Zernike self-consistent field [MC-MOZ-SCF]; Kido et al., J. Chem. Phys. 2015, 143, 014103). To save the computational time, we have also developed a fast method to generate electrostatic potential map near solute and the solvation term in Fock matrix, using Fourier transformation (FT) and real spherical harmonics expansion (RSHE). To numerically validate the LRA and FT-RSHE method, we applied the present approach to water, carbonic acid, and their ionic species in aqueous solution. Molecular properties of the solutes were evaluated by the present approach with four different types of initial wave functions and compared with those by the original (MC-MOZ-SCF). We found that an initial wave function considering solvation effects is needed to appropriately reproduce the properties by MC-MOZ-SCF. Furthermore, a benchmark test for 32 solute molecules was performed to evaluate the accuracy of the present approach for solvation free energy (SFE) and measure the speedup ratio for MC-MOZ-SCF. The error of SFE for MC-MOZ-SCF does not correlate with the SFE but increases in proportion to the electronic reorganization energy. Similar to water and carbonic acid, an initial wave function with solvation effects is also important to make the error small. From the averaged speed up ratio, the present approach is 13.5 times faster than MC-MOZ-SCF. © 2019 Wiley Periodicals, Inc.     

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.     

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.     

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.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

离子液中溶剂化电子的结构及其演化动力学

本文综述了离子型液体介质中过剩电子的结构、存在状态及其时间演化动力学特征。基于从头算分子动力学模拟及计算结果,重点阐述了咪唑型、吡啶型、碱金属离子型熔盐氯化物离子液中与过剩电子溶剂化密切相关的溶剂化能量学、结构特征、可能的存在状态以及态-态转化稳态动力学机制,分析了此类离子型介质中电子高效传导的内在本质及离子液组成离子的重要作用。阳离子的最低未占轨道组成的导带结构是离子液中过剩电子的溶剂化态及其稳定性的决定因素,任何能影响或改变其导带结构的因素均能显著影响过剩电子溶剂化。但快速的态-态转化及电子迁移并不明显取决于其组成离子扩散动力学,而是敏感地受离子液涨落所控制。这种基于溶剂化电子的迁移模式构成了此类离子型介质甚至其它液态介质中电子转移的新途径。     

Insights into Tris‐(2‐Hydroxylethyl)methylammonium Methylsulfate Aqueous Solutions

We report herein a combined experimental–computational study on tris‐(2‐hydroxylethyl)methylammonium methylsulfate in water solutions, as a representative ionic liquid of the aqueous‐solution behavior of hydroxylammonium‐based ionic liquids. Relevant thermophysical properties were measured as a function of mixture composition and temperature. Classical molecular dynamics simulations were performed to infer microscopic structural features. The reported results for ionic liquid in water‐rich solutions show that it behaves as isolated non‐interacting ions solvated by water molecules, through well‐defined solvation shells, exerting a disrupting effect on the water hydrogen bonding network. Nevertheless, as ionic liquid concentration increase, interionic association increases, even for diluted water solutions, evolving from the typical behavior of strong electrolytes in solution toward large interacting structures. For ionic‐liquid‐rich mixtures, water exerts a minor disrupting effect on the fluid’s structuring because it occupies regions around each ion (developing water–ion hydrogen bonds) but without significantly weakening anion–cation interactions.     

Solvent response and dielectric relaxation in supercooled butyronitrile

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming dipolar liquid butyronitrile near its glass transition temperature T(g)=95 K. The Stokes shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical linewidth changes along the solvation coordinate in a nonmonotonic fashion. These features are characteristic of solvation in viscous solvents and clearly indicate heterogeneous dynamics, i.e., spatially distinct solvent response times. Using the dielectric relaxation data of viscous butyronitrile as input, a microscopic model of dipolar solvation captures the relaxation time, the relaxation dispersion, and the amplitude of the dynamical Stokes shift remarkably well.     

Dynamics of solvent and rotational relaxation of coumarin 153 in room-temperature ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate confined in Brij-35 micelles: a picosecond time-resolved fluorescence spectroscopic study

The dynamics of solvent and rotational relaxation of Coumarin 153 (C-153) in ionic liquid (IL) 1-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6]) and in the ionic liquid confined in Brij-35 micellar aggregates have been investigated using steady-state and time-resolved fluorescence spectroscopy. We observed slower dynamics in the presence of micellar aggregates as compared to the pure IL. However, the slowing down in the solvation time on going from neat IL to IL-confined micelles is much smaller compared to that on going from water to water-confined micellar aggregates. The increase in solvation and rotational time in micelles is attributed to the increase in viscosity of the medium. The slow component is assumed to be dependent on the viscosity of the solution and involves large-scale rearrangement of the anions and cations while fast component is assumed to originate from the initial response of the anions during excitation. The slow component increases due to the increase in the viscosity of the medium and increase in fast component is probably due to the hydrogen bonding between the anions and polar headgroup of the surfactant. The dynamics of solvent relaxation was affected to a small extent due to the micelle formation.     

Femtosecond solvation dynamics in a neat ionic liquid and ionic liquid microemulsion: excitation wavelength dependence

Solvation dynamics in a neat ionic liquid, 1-pentyl-3-methyl-imidazolium tetra-flouroborate ([pmim][BF4]) and its microemulsion in Triton X-100 (TX-100)/benzene is studied using femtosecond up-conversion. In both the neat ionic liquid and the microemulsion, the solvation dynamics is found to depend on excitation wavelength (lambda(ex)). The lambda(ex) dependence is attributed to structural heterogeneity in neat ionic liquid (IL) and in IL microemulsion. In neat IL, the heterogeneity arises from clustering of the pentyl groups which are surrounded by a network of cation and anions. Such a nanostructural organization is predicted in many recent simulations and observed recently in an X-ray diffraction study. In an IL microemulsion, the surfactant (TX-100) molecules aggregate in form of a nonpolar peripheral shell around the polar pool of IL. The micro-environment in such an assembly varies drastically over a short distance. The dynamic solvent shift (and average solvation time) in neat IL as well as in IL microemulsions decreases markedly as lambda(ex) increases from 375 to 435 nm. In a [pmim][BF4]/water/TX-100/benzene quaternary microemulsion, the solvation dynamics is slower than that in a microemulsion without water. This is ascribed to the smaller size of the water containing microemulsion. The anisotropy decay in an IL microemulsion is found to be faster than that in neat IL.