Heterogeneous Orientational Relaxations and Translation–Rotation Decoupling in (Choline Chloride + Urea) Deep Eutectic Solvents: Investigation through Molecular Dynamics Simulations and Dielectric Relaxation Measurements

Molecular dynamics insights into the dynamical behavior of structurally modified water in aqueous deep eutectic solvents (ADES)

Recent studies have demonstrated that the presence of water in deep eutectic solvents (DESs) significantly affects their dynamics, structure, and physical properties. Although the structural changes due to the addition of water are well understood, the microscopic dynamics of these changes have been rarely studied. Here, we performed molecular dynamics simulation of 30% (v/v) (∼0.57 molar fraction) water mixture of DES containing CH3CONH2 and NaSCN/KSCN at various salt fractions to understand the microscopic structure and dynamics of water. The simulated results reveal a heterogeneous environment for water molecules in aqueous DES (ADES), which is influenced by the nature of the cation. The diffusion coefficients of water in ADESs are significantly lower than that in neat water and concentrated aqueous NaSCN/KSCN solution. When Na+ ions are replaced by K+ ions in the ADES system, the diffusion coefficient increases, which is consistent with the measured nuclear magnetic resonance data. Self-dynamic structure factor for water and other simulated dynamic quantities, such as reorientation, hydrogen-bond, and residence time correlation functions, show markedly slower dynamics inside ADES than in the neat water and aqueous salt solution. Moreover, these dynamics become faster when Na+ ions in ADES are replaced by K+ ions. The results suggest that the structural environment of water in Na+-rich ADES is rigid due to the presence of cation-bound water and geometrically constrained water. The medium becomes less rigid as the KSCN fraction increases due to the relatively weaker interaction of K+ ions with water than Na+ ions, which accelerates the dynamical processes.

Temperature-dependent dielectric relaxation measurements of (acetamide + K/Na SCN) deep eutectic solvents: Decoding the impact of cation identity via computer simulations

The impact of successive replacement of K+ by Na+ on the megahertz-gigahertz polarization response of 0.25[fKSCN + (1 - f)NaSCN] + 0.75CH3CONH2 deep eutectic solvents (DESs) was explored via temperature-dependent (303 ≤ T/K ≤ 343) dielectric relaxation (DR) measurements and computer simulations. Both the DR measurements (0.2 ≤ ν/GHz ≤ 50) and the simulations revealed multi-Debye relaxations accompanied by a decrease in the solution static dielectric constant (ɛs) upon the replacement of K+ by Na+. Accurate measurements of the DR response of DESs below 100 MHz were limited by the well-known one-over-frequency divergence for conducting solutions. This problem was tackled in simulations by removing the zero frequency contributions arising from the ion current to the total simulated DR response. The temperature-dependent measurements revealed a much stronger viscosity decoupling of DR times for Na+-containing DES than for the corresponding K+ system. The differential scanning calorimetry measurements indicated a higher glass transition temperature for Na+-DES (∼220 K) than K+-DES (∼200 K), implying more fragility and cooperativity for the former (Na+-DES) than the latter. The computer simulations revealed a gradual decrease in the average number of H bonds (⟨nHB⟩) per acetamide molecule and increased frustrations in the average orientational order upon the replacement of K+ by Na+. Both the measured and simulated ɛs values were found to decrease linearly with ⟨nHB⟩. Decompositions of the simulated DR spectra revealed that the cation-dependent cross interaction (dipole-ion) term contributes negligibly to ɛs and appears in the terahertz regime. Finally, the simulated collective single-particle reorientational relaxations and the structural H-bond fluctuation dynamics revealed the microscopic origin of the cation identity dependence shown by the measured DR relaxation times.

Absorption of Volatile Organic Compounds Toluene and Acetaldehyde in Choline Chloride-Based Deep Eutectic Solvents

Unrestricted emission of volatile organic compounds (VOCs)─a threat to human health and the environment─can be controlled to a large extent by the capturing mechanism. Few recent experimental studies explored the efficacy of the deep eutectic solvent (DES), a designer solvent with some fascinating properties, as a VOC-capturing medium. Through the partition coefficient measurement, it was found that the choline chloride-based DESs exhibit excellent VOC-capturing potencies. However, a molecular picture of the above absorption process is still lacking. Here, we study the molecular mechanism of the absorption of two commonly occurring VOCs, toluene and acetaldehyde, in two different choline chloride-based DESs with varying donor molecules, urea, and levulinic acid via the molecular dynamics simulation technique. Strong absorption of the VOCs is observed in both the DESs. The free energy profile for the absorption process has been explored using the umbrella sampling method. The VOCs are preferentially solvated near the liquid/vapor interface. The simulated partition coefficients for the VOCs from the vapor to the liquid phase show good agreement with the experimental results. Detailed analyses of the spatial and orientational structure of the VOCs and different components of DESs are performed to elucidate the interaction among them. The above analyses have indicated that DES is a better VOC-capturing medium compared to a room-temperature ionic liquid, which is more extensively studied in the literature.

Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

Structural properties of supercooled deep eutectic solvents: choline chloride–thiourea compared to reline

Deep eutectic solvents (DESs) are eutectic mixtures of hydrogen bond acceptors and hydrogen bond donors which melt at much lower temperatures than the individual components. DESs attract growing interest because...

The Role of Hydrogen Bonding in the Preferential Solvation of 5-Aminoquinoline in Binary Solvent Mixtures

5-Aminoquinoline (5AQ) has been used as a fluorescent probe of preferential solvation (PS) in binary solvent mixtures in which the nonpolar component is diethyl ether and the polar component is protic (methanol) or aprotic (acetonitrile). Hence, the roles of solvent polarity and solute-solvent hydrogen bonding have been delineated. Positive deviations of spectral shifts from a linear dependence on the concentration of the polar component, signifying PS, are markedly more pronounced in case of the protic solvent. Solvation dynamics on a nanosecond time scale mark the formation of the solvation shell around the fluorescent probe. Time-resolved area-normalized emission spectra indicate the occurrence of the continuous solvation of the excited state when the polar component is acetonitrile. In contrast, two distinct states were observed when the polar component was methanol, the second state being the hydrogen bonded one. Translational diffusion is the rate-determining step for formation of the solvation shell. The time constant associated with it has been estimated from rise times observed in fluorescence transients monitored at the red end of the fluorescence spectra and also from the time evolution of the spectral width of time-resolved emission spectra.

Temperature-Dependent Dielectric Relaxation in Ionic Acetamide Deep Eutectics: Partial Viscosity Decoupling and Explanations from the Simulated Single-Particle Reorientation Dynamics and Hydrogen-Bond Fluctuations

We report here temperature-dependent (293 ≤ T (K) ≤ 336) dielectric relaxation (DR) measurements of (acetamide + LiBr/NO₃^-/ClO₄^-) deep eutectic solvents (DESs) in the frequency window of 0.2 ≤ ν (GHz) ≤ 50 and explore, via molecular dynamics simulations, the relative roles for the collective single-particle reorientational relaxations and the H-bond dynamics of acetamide in the measured DR response. In addition, DR measurements of neat molten acetamide were performed. Recorded DR spectra of these DESs require multi-Debye fits and produce well-separated DR time scales that are spread over several picoseconds to ∼1 ns. Simulations suggest DR time scales derive contributions from both the collective reorientational (Cl(t)) relaxation and structural H-bond (C_HB(t)) dynamics of acetamide. A good correlation between the measured and simulated activation energies further reveals a strong connection between the measured DR and the simulated Cl(t) and C_HB(t). Average DR times exhibit a strong fractional viscosity dependence, suggesting substantial microheterogeneity in these media. Simulations of Cl(t) and C_HB(t) reveal strong stretched exponential relaxations with a stretching exponent, 0.4 ≤ β ≤ 0.7. The ratio between the average reorientational correlation times of first and second ranks, ⟨τ⟩l=1/⟨τ⟩l=2, deviates appreciably from Debye's l(l+1) law for homogeneous media. Importantly, a pronounced translation-rotation decoupling between the simulated reorientation and center-of-mass diffusion times was observed.

Assessing rotation and solvation dynamics in ethaline deep eutectic solvent and its solutions with methanol

Steady-state and time-resolved fluorescence were used to investigate the solvation of coumarin 153 (C153) and coumarin 343 (C343) in methanol + ethaline binary solutions, a deep eutectic solvent composed of a 1:2 molar ratio choline chloride + ethylene glycol. In addition, time-resolved anisotropy decays were used to determine the solute's rotational reorientation time as a function of viscosity. Measurements were made in solutions covering the entire range of mole fraction. Viscosity measurements were used to characterize the bulk solvent properties, and as expected, addition of methanol resulted in an decreased viscosity, showing an exponential decrease with mole fraction, up to ∼50-fold at x_MeOH = 1.0. Probe rotational reorientation times were found to be biexponential at x_MeOH < 0.3 for C153 and x_MeOH < 0.5 for C343 and monoexponential at richer methanol content. In proportion to viscosity, C153 and C343 average rotation times decreased ∼30-fold from x_MeOH = 0 to 0.9 and showed a power law dependence of ∼η^0.85. Rotation times approached the stick boundary limit on dilution with methanol. Time-resolved Stokes shifts quantified the solvation dynamics and were nearly single exponential for C153 but were clearly biexponential for C343. Solvation times also tracked with viscosity according to a power law dependence, with exponents of 0.3 and 0.4 for C153 and C343, respectively. The dilution effect of methanol was not linear in proportion to the viscosity change and alone cannot account for the change in solvation. Dilution also showed a different correlation to solvation than did temperature variations to govern the viscosity change.

Dynamic heterogeneity and viscosity decoupling: origin and analytical prediction

The molecular-level structure and dynamics decide the functionality of solvent media. Therefore, a significant amount of effort is being dedicated continually over time in understanding their structural and dynamical features. One intriguing aspect of solvent structure and dynamics is heterogeneity. In these systems, the dynamics follow , where p is the measure of viscosity decoupling. We analytically predicted that in such cases, the Stokes-Einstein relationship is modified to due to microdomain formation, and the second term on the right-hand side leads to viscosity decoupling. We validated our prediction by estimating the p values of a few solvents, and they matched well with the literature. Overall, we believe that our approach gives a simple yet unique physical picture to help us understand the heterogeneity of solvent media.