Slow solvation in ionic liquids: Connections to non-Gaussian moves and multi-point correlations

Solute dynamics of a hydrophobic molecule in a menthol-thymol based type-V deep eutectic solvent: effect of composition of the components

Type-V deep eutectic solvents (DESs) are a newly emerging unique class of solvents obtained by physical mixing and heating of non-ionic components. These solvents show deviation from the thermodynamic ideality. Compared to type-I to IV DESs, type-V DESs are less explored and their physical chemistry is in its nascent stage. In this work, we have chosen a type-V DES based on menthol-thymol (MT) for our working media. Solvent and rotational dynamics were studied with varying temperature using a well-known solvatochromic probe, Coumarin 153 (C153). We prepared the MT-based DES using a reported procedure at three molar ratios: 1 : 1 (M1T1), 1 : 1.5 (M1T1.5), and 2 : 1 (M2T1) of menthol (M) and thymol (T). Time-resolved emission spectra (TRES) were constructed with varying temperature. Utilizing TRES, the decay of the solvent correlation function (C(t)) was plotted. We have correlated the solvent relaxation time in these DESs as a function of viscosity. The time-resolved anisotropy decays were also collected to perceive the rotational relaxation dynamics of C153 as a function of temperature. The decay of solvent relaxation was found to be bi-exponential, and the average solvation time (〈τs〉) in M2T1 was found to be longer than those of M1T1.5 and M1T1. The rotational reorientation times (〈τrot〉) also follow the same trend. We have analysed the rotational dynamics of C153 in type-V DESs employing the Stokes-Einstein-Debye (SED) hydrodynamic model. The rotational dynamics in DESs demonstrate a good correlation with the SED model with a little deviation. In MT-based DESs, the solute's rotational relaxation times approach hydrodynamic stick boundary condition at low viscosity (or at high temperatures) for all molar compositions. Using the Arrhenius-type equations, we have correlated the activation energies for the rotational motion of C153, along with the viscous flow and non-radiative pathways for all the DESs.

Difference in "Supercooling" Affinity between (Acetamide + Na/KSCN) Deep Eutectics: Reflections in the Simulated Anomalous Motions of the Constituents and Solution Microheterogeneity Features

Deep depression of freezing points of ionic amide deep eutectic solvents (DESs) is known to exhibit a significant dependence on the identity of ions present in those systems and the nature of the functional group attached to the host amide. This deep depression of the freezing point is sometimes termed as "supercooling". For (acetamide + electrolyte) DESs, experiments have revealed signatures of ion-dependent spatiotemporal heterogeneity features. The focus of this work is to provide microscopic explanations of these experimentally observed macroscopic system properties in terms of particle jumps and insights about the origin of the cation dependence. For this purpose, extensive molecular dynamics simulations have been performed employing (acetamide + Na/KSCN) deep eutectics as representative ionic systems at 303, 318, 333, and 348 K. The individual translational motions of acetamide and the ions are followed, and their connections to solution heterogeneity are explored. The center-of-mass motion for Na⁺ has been found to be more anomalous than that for K⁺. This difference corroborates well with experimental reports on heterogeneous relaxations in these systems. Simulated viscosity coefficients and dynamic heterogeneity features also reflect this difference. Moreover, simulated reorientational relaxations of acetamide molecules in these ionic DESs suggest that a Na⁺-containing DES is more heterogeneous than the corresponding K⁺-containing system. Estimated void and neck distributions for acetamide molecules differ as the alkali metal ions differ. In brief, this study provides a detailed microscopic view of the cation dependence of the microheterogeneous relaxation dynamics of these DESs reported repeatedly by different experiments.

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

Molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. In this work, we develop and assess two multiscale solvation models for computing solvation thermodynamic properties. The new SLIC/CDC model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with new statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). Given the structures of 500 solutes, the SLIC/CDC model predicts Gibbs energies of solvation in water with an average accuracy better than 1 kcal/mol, when compared to experimental measurements, and better than 0.8 kcal/mol, when compared to explicit-solvent molecular dynamics simulations. The individual SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The previously published SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent-accessible surface area (SASA) nonpolar energy mode. With our new, improved parametrization method, the SLIC/SASA model now predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems, compared to experimental and explicit-solvent molecular dynamics, and better than 1.6 kcal/mol average accuracy in ionic liquids, compared to explicit-solvent molecular dynamics. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

Solvent-Assisted Li-Ion Transport and Structural Heterogeneity in Fluorinated Battery Electrolytes

The electrolytes diluted with fluorinated solvents show promising properties toward better battery technology than existing ones. The transport of Li ions in these fluorinated electrolytes is essential to access the performance of a battery. It is believed that the transport of the Li ion in these electrolytes occurs through polar solvents in the matrix of nonpolar solvent molecules. The atomistic details of this mechanism are yet to be proved using the dynamics of these mixtures. In this study, we performed classical molecular dynamics simulations at various temperatures to probe this mechanism through the structure and dynamics of electrolytes at the atomic level. Here, we have shown that the polar fluorinated solvents assist the Li-ion transport in a region of nonpolar solvent. Highly polar molecules also solvate the Li ion at a lower temperature. The nonpolar solvent solvates the Li ion weakly as compared to others. The calculated values of the ionic conductivity from the Green-Kubo relation provide a better match than that from an experimental conductivity meter. Furthermore, we probed the heterogeneity in the dynamics of the electrolytes by calculating the non-Gaussian parameter. We also show that the transport mechanism of the Li ion in diluted concentrated electrolytes is different than a few of the other reported electrolytes. We have also calculated the ion-pair and ion-cage lifetimes to see the most and least lived ion/ion-solvent pairs. The mechanism given from the present study may help to design the fluorinated electrolytes for Li-ion batteries.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

Nonlinear measurements of kinetics and generalized dynamical modes. II. Application to a simulation of solvation dynamics in an ionic liquid

Solvation dynamics in ionic liquids show features that are often associated with supercooled liquids, including "stretched" nonexponential relaxation. To better understand the mechanism behind the stretching, the nonlinear mode-correlation methods proposed in Paper I [S. R. Hodge and M. A. Berg, J. Chem. Phys. 155, 024122 (2021)] are applied to a simulation of a prototypical ionic liquid. A full Green's function is recovered. In addition, specific tests for non-Gaussian dynamics are made. No deviations from Gaussian dynamics are found. This finding is incompatible with rate heterogeneity as a cause of the nonexponential relaxation and appears to be in conflict with an earlier multidimensional analysis of the same data. Although this conflict is not resolved here, this work does demonstrate the practicality of mode-correlation analysis in the face of finite datasets and calculations.

Breakdown of the Stokes-Einstein Relation in Supercooled Water/Methanol Binary Mixtures: Explanation Using the Translational Jump-Diffusion Approach

A recent experiment has directly checked the validity of the Stokes-Einstein (SE) relation for pure water, pure methanol, and their binary mixtures of three different compositions at different temperatures. The effect of composition on the nature of breakdown of the SE relation is interesting. While in the majority of the systems, an increasing SE breakdown is observed with decreasing temperature, the breakdown is already significant at higher temperatures for the equimolar mixture. Violations of the SE relation in pure supercooled water at different temperatures and pressures have been previously explained using the translational jump-diffusion (TJD) approach, which provides a fundamental molecular basis, by directly connecting the SE breakdown with jump-diffusion of the molecules. We have used the same TJD approach for explaining the SE breakdown for the methanol/water binary mixtures of compositions studied in the experiment over a wide range of temperatures between 220 K and 300 K. We have understood that the jump-diffusion is the key responsible factor for the SE breakdown. The maximum jump-diffusion contribution gives rise to the early SE breakdown observed for the equimolar mixture observed in the experiment. This study, therefore, provides molecular insight into the SE breakdown for the supercooled water/methanol binary mixture, as found in the experiment.

Heterogeneity in the microstructure and dynamics of tetraalkylammonium hydroxide ionic liquids: insight from classical molecular dynamics simulations and Voronoi tessellation analysis

Microscopic structural and dynamic heterogeneities were investigated for three ionic liquids (ILs), tetraethylammonium hydroxide, tetrapropylammonium hydroxide, and tetrabutylammonium hydroxide employing classical molecular dynamics (MD) simulations.

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.

Molecular dynamics simulation of geminal dicationic ionic liquids [Cn(mim)2][NTf2]2 – structural and dynamical properties

The structural and dynamical properties of two dicationic ionic liquids, i.e. [C_n(mim)₂][NTf₂]₂ with n = 3 and 5, have been studied to obtain a fundamental understanding of the molecular basis of the macroscopic and microscopic properties of the bulk liquid phase.

Fluorescence Dynamics in the Endoplasmic Reticulum of a Live Cell: Time-Resolved Confocal Microscopy

Fluorescence dynamics in the endoplasmic reticulum (ER) of a live non-cancer lung cell (WI38) and a lung cancer cell (A549) are studied by using time-resolved confocal microscopy. To selectively study the organelle, ER, we have used an ER-Tracker dye. From the emission maximum (λmaxem) of the ER-Tracker dye, polarity (i.e. dielectric constant, ϵ) in the ER region of the cells (≈500 nm in WI38 and ≈510 nm in A549) is estimated to be similar to that of chloroform (λmaxem =506 nm, ϵ≈5). The red shift by 10 nm in λmaxem in the cancer cell (A549) suggests a slightly higher polarity compared to the non-cancer cell (WI38). The fluorescence intensity of the ER-Tracker dye exhibits prolonged intermittent oscillations on a timescale of 2-6 seconds for the cancer cell (A549). For the non-cancer cell (WI38), such fluorescence oscillations are much less prominent. The marked fluorescence intensity oscillations in the cancer cell are attributed to enhanced calcium oscillations. The average solvent relaxation time (<τs >) of the ER region in the lung cancer cell (A549, 250±50 ps) is about four times faster than that in the non-cancer cell (WI38, 1000±50 ps).

Heterogeneous dynamics and its length scale in simple ionic liquid models: a computational study

We numerically investigate the dynamic heterogeneity and its length scale found in coarse-grained ionic liquid model systems.

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.

Time scale of dynamic heterogeneity in model ionic liquids and its relation to static length scale and charge distribution

We find a general power-law behavior: , where ζ_dh ≈ 1.2 for all the ionic liquid models, regardless of charges and the length scale of structural relaxation.

Study on the temperature-dependent coupling among viscosity, conductivity and structural relaxation of ionic liquids

The frequency-dependent viscosity and conductivity of three imidazolium-based ionic liquids were measured at several temperatures in the MHz region, and the results are compared with the intermediate scattering functions determined by neutron spin echo spectroscopy.

Solvent effects of 1-ethyl-3-methylimidazolium acetate: solvation and dynamic behavior of polar and apolar solutes

We study the solvation properties of the ionic liquid 1-ethyl-3-methylimidazolium acetate ([eMIM]⁺[ACE]⁻) and the resulting dynamic behavior for differently charged model solutes at room temperature via atomistic molecular dynamics (MD) simulations of 500 ns length.