OPLS Force Field for Choline Chloride-Based Deep Eutectic Solvents

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst-Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

Molecular Understanding of CO2 Absorption by Choline Chloride/Urea Confined within Nanoslits

Clarifying the potential relationship between the microstructure of nanoconfined choline chloride/urea (ChClU) and CO2 absorption performance is key to understanding the abnormal increase in CO2 under nanoconfinement. In this study, we used molecular dynamics simulations and grand canonical Monte Carlo (GCMC) to systematically study the mechanism underlying the absorption of CO2 by ChClU within nanoslits. According to the spatial distribution, ChClU can form two different laminar regions within nanoslits, namely, the interfacial region (region I) and beyond region I (region II). In region II, the interface induces rearrangement of ChClU, resulting in an increase in free volume and subsequent increase in CO2 solubility. In region I, changing the interface from hydrophobic to hydrophilic (e.g., S_I to S_IV) by setting the appropriate charge patterns, the urea molecules gradually change from "disordered" to "ordered standing" relative to the solid surface. The preferential orientation of the urea molecules causes competition between the ChClU's free volume and urea molecules, resulting in a non-monotonic change in CO2 solubility. Specifically, from S_I to S_III, the increase in urea molecules enhances the CO2 solubility. In S_IV, space for CO2 absorption is insufficient due to the accumulation of urea molecules, and thus CO2 solubility decreases.

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. The simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

Influence of Deep Eutectic Solvent Composition on Micelle Properties: A Molecular Dynamics Study

This study investigates the structural and transport properties of SDS, CTAB, and SB3-12 micelles in three deep eutectic solvents (DESs), Ethaline, Glyceline, and Reline, using molecular dynamics (MD) simulations. The influence of solvent composition on micelle morphology, interactions, and dynamics was explored, revealing key differences driven by the DES environment. Structural analyses, including eccentricity and radius of gyration, demonstrated that micelle shape and compactness vary significantly depending on the solvent. In Ethaline and Reline, larger micelles showed significant deviations from spherical shapes, while micelles in Glyceline became more spherical and compact, particularly those formed by SB3-12. Radial distribution functions highlighted different levels of micelle-solvent interactions, with SDS showing strong interactions with HBD components and SB3-12 exhibiting prominent self-interaction. According to hydrogen bonding analysis, micelles slightly disrupt the DES hydrogen bond network, with SB3-12 establishing the most significant hydrogen bond connections. The transport property analysis revealed that larger micelles have lower diffusion coefficients, whereas smaller micelles enhance DESs' component mobility. These findings advance the understanding of micelle behavior in DESs and also help in the optimization of DES-surfactant systems for applications such as electrodeposition, nanomaterial templating, and drug delivery. Future research will focus on surfactant interactions with surfaces to further improve these applications.

Temperature Dependence of Intermolecular Dynamics and Liquid Properties of Deep Eutectic Solvent, Reline

We investigated the temperature dependence of the intermolecular dynamics, including intermolecular vibrations and collective orientational relaxation, of one of the most typical deep eutectic solvents, reline, using femtosecond Raman-induced Kerr effect spectroscopy (fs-RIKES), subpicosecond optical Kerr effect spectroscopy (ps-OKES), and molecular dynamics (MD) simulations. According to fs-RIKES results, the temperature-dependent intermolecular vibrational band peak at ∼90 cm-1 exhibited a redshift with increasing temperature. The density-of-state (DOS) spectrum of reline by MD simulations reproduced this fs-RIKES spectral feature. The decomposition analysis of the DOS spectra showed that all constituent components, including urea, cholinium cation, and chloride anion, also exhibited similar magnitudes of redshifts upon heating, indicating that the three species intermolecularly interact one another. The temperature sensitivity of the intermolecular vibrational frequency was high compared to that of ionic liquids. According to ps-OKES results, the slow orientational relaxation rate increased with increasing temperature; however, this phenomenon was not well explained by the Stokes-Einstein-Debye hydrodynamic model. Analysis of the orientational relaxation time based on the Stokes-Einstein-Debye model indicates that the decoupling between the orientational relaxation time and viscosity occurs at temperatures below ∼330 K. Quantum chemistry calculations of urea and the cholinium cation based on the MP2/6-311++G(d,p) level of theory confirmed that the contribution of intramolecular vibrational bands to low-frequency bands below 200 cm-1 was minimal. The densities, viscosities, electrical conductivities, and surface tensions of reline at various temperatures were also estimated and compared with the dynamics data.

Deep Eutectic Solvents Meet Non-Aqueous Cosolvents: A Modeling and Simulation Perspective-A Tutorial Review

Deep eutectic solvents (DESs) have recently gained attention due to their tailorable properties and versatile applications in several fields, including green chemistry, pharmaceuticals, and energy storage. Their tunable properties can be enhanced by mixing DESs with cosolvents such as ethanol, acetonitrile, and water. DESs are structurally complex, and molecular modeling techniques, including quantum mechanical calculations and molecular dynamics simulations, play a crucial role in understanding their intricate behavior when mixed with cosolvents. While the most studied cosolvent is water, in some applications, even a small content of water is considered a contaminant, for example, when the processes of interest require dry conditions. Only quite recently have modeling studies begun to focus on DES mixed with cosolvents other than water. This tutorial provides the first comprehensive overview of these studies. It highlights how modern molecular modeling increases our understanding of their structural organization, transport properties, phase behavior, and thermodynamic properties. Additionally, case studies and recent developments in the field are discussed along with the challenges and future directions in molecular modeling of DES in cosolvent mixtures. Overall, this review offers valuable insights into the molecular-level understanding of DES-cosolvent systems and their implications for designing novel solvent mixtures with tailored properties for various applications.

Solvation Dynamics and Microheterogeneity in Deep Eutectic Solvents

Deep eutectic solvents have attracted considerable attention due to their unique properties and their potential to replace conventional solvents in diverse applications, such as catalysis, energy storage, and green chemistry. However, despite their broad use, the microscopic mechanisms governing solvation dynamics and the role of hydrogen bonding in deep eutectic solvents remain insufficiently understood. In this article, we present our contributions toward unravelling the micro heterogeneity within deep eutectic solvents by combining vibrational Stark spectroscopy and two-dimensional infrared spectroscopy with molecular dynamics simulations. Our findings demonstrate how the composition, constituents, and addition of water significantly influence the heterogeneous hydrogen bonding network and solvent dynamics within these systems. These insights provide valuable guidance for the design of next-generation solvents tailored to specific applications. By integrating experimental and computational approaches, this work sheds light on the intricate relationship between solvation dynamics and nanostructure in deep eutectic solvents, ultimately paving the way for innovative advances in solvent design.

Assessing Melting and Solid-Solid Transition Properties of Choline Chloride via Molecular Dynamics Simulations

Choline chloride (ChCl) is used extensively as a hydrogen bond donor in deep eutectic solvents (DESs). However, determining its melting properties experimentally is challenging due to decomposition upon melting, leading to widely varying literature values. Accurate melting properties are crucial for understanding the solid-liquid phase behavior of ChCl-containing DESs. Here, we employ molecular dynamics simulations to compute the phase transitions of ChCl, testing a variety of atomistic force fields. We find that the results are sensitive to the choice of force field, but a melting temperature of 627 K and a melting enthalpy of 7.8 kJ/mol seem most reasonable, in good agreement with some literature values. We suggest these as the likely melting properties of ChCl, though the results are tentative due to limited experimental data for the liquid ChCl phase.

Extraction mechanism of phenolic compounds by a choline chloride/glycerol solvent: DFT and molecular dynamics studies

The mechanism of a solvent consisting of choline chloride and glycerol (ChCl/GLY) for extracting phenolic compounds from coal tar was theoretically studied using density functional theory calculations and molecular dynamics simulations. The thermodynamic properties, interaction essence, and molecular dynamics properties of the extraction system were investigated, as well as the influence of ChCl/GLY on the vibration spectra of phenolic compounds. The results show that the solvation free energy of phenolic compounds in ChCl/GLY is more negative than that in coal tar, leading to the spontaneous transfer of phenolic compounds from coal tar to ChCl/GLY. The electrostatic and dispersion interactions between phenolic compounds and ChCl/GLY have similar significance in the extraction process, with interaction energies ranging from -46 to -53 kJ mol-1. The mixing of phenolic compounds with ChCl/GLY has minimal impact on their internal molecular structure, however, it does reduce the diffusion coefficients of each component in ChCl/GLY and shortens the lifetime of hydrogen bonds in both phenolic compounds and ChCl/GLY. The first shell of each phenolic compound is surrounded by 1.15 chloride ions. Following dissolving in ChCl/GLY, the stretching vibration peaks of phenolic compounds, namely the -OH and C-H/-CH3 regions, undergo a shift. The results enhance comprehension of the extraction process of phenolic compounds by DES.

Force Fields for Deep Eutectic Mixtures: Application to Structure, Thermodynamics and 2D-Infrared Spectroscopy

Parametrizing energy functions for ionic systems can be challenging. Here, the total energy function for an eutectic system consisting of water, SCN-, K+ and acetamide is improved vis-a-vis experimentally measured properties. Given the importance of electrostatic interactions, two different types of models are considered: the first (model M0) uses atom-centered multipole whereas the other two (models M1 and M2) are based on fluctuating minimal distributed charges (fMDCM) that respond to geometrical changes of SCN-. The Lennard-Jones parameters of the anion are adjusted to best reproduce experimentally known hydration free energies and densities, which are matched to within a few percent for the final models irrespective of the electrostatic model. Molecular dynamics simulations of the eutectic mixtures with varying water content (between 0 and 100%) yield radial distribution functions and frequency correlation functions for the CN-stretch vibration. Comparison with experiments indicates that models based on fMDCM are considerably more consistent than those using multipoles. Computed viscosities from models M1 and M2 are within 30% of measured values and their change with increasing water content is consistent with experiments. This is not the case for model M0.

Osmolyte-induced conformational stabilization of a hydrophobic polymer

Elucidating the mechanistic role of osmolytes on conformations of hydrophobic prototypical macromolecules in principle is the stepping stone towards understanding the effect of osmolytes on proteins. Motivated by this, we use equilibrium simulations and umbrella sampling techniques to dissect the underlying mechanism of osmolyte-induced conformational stability of a hydrophobic polymer. Our results unveil a remarkable osmolyte-dependent conformational stabilization of the polymer. In an aqueous solution of 4 M choline chloride (ChCl), the polymer has an even more compact structure than in water. On the other hand, an aqueous solution of 8 M urea stabilizes the extended state of the polymer. Interestingly, the polymer adopts an intermediate hairpin conformation in a mixed osmolyte solution of 4 M ChCl and 8 M urea in water due to the interplay of ChCl and urea. Our simulations identify the relative accumulation of water and the hydrophilic part of choline or preferential binding of urea near the collapsed and the extended states, respectively. Analyses split out the enthalpic and entropic contributions to the overall free energy. This decides the stabilization of the preferred conformation in the chosen osmolyte solution. Our simulations show that in an aqueous solution of ChCl, the hairpin state is stabilized by entropy gain. In contrast, the enthalpic contribution stabilizes the hairpin state in mixed environments. However, a collapsed state is energetically not favored in the presence of urea. In brief, via employing an in silico approach, the current findings indicate the importance of osmolytes in stabilizing the conformational states of hydrophobic polymers.

Combined experimental and computational investigation of tetrabutylammonium bromide-carboxylic acid-based deep eutectic solvents

Aiming at shedding light on the molecular interactions in deep eutectic solvents (DESs), the DESs based on tetrabutylammonium bromide (TBAB) as hydrogen bond acceptor (HBA) and carboxylic acids (CAs) (formic acid (FA), oxalic acid (OA), and malonic acid (MA)) as hydrogen bond donor (HBD) were investigated by both experimental and theoretical techniques. The thermal behaviors of the prepared DESs were investigated by differential scanning calorimetry (DSC) method. In order to study the hydrogen bond formation between the DESs constituents, the FT-IR analysis was carried out. The large positive deviations of the iso solvent activity lines of ternary HBA + HBD + 2-propanol mixtures determined by the isopiestic technique from the semi-ideal behavior indicate that CAs interact strongly with TBAB and therefore they can form DESs. Molecular dynamics (MD) simulations were performed to present an atomic-scale image of the components and describe the microstructure of DESs. From the MD simulations, the radial distribution functions (RDFs), coordination numbers (CNs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) were calculated to investigate the interaction between the components and three-dimensional visualization of the DESs. The obtained results confirmed the importance of hydrogen bonds in the formation of TBAB/CAs DESs.

A case study using spectroscopy and computational modelling for Co speciation in a deep eutectic solvent

Cobalt has a vital role in the manufacturing of reliable and sustainable clean energy technologies. However, the forecasted supply deficit for cobalt is likely to reach values of 150 kT by 2030. Therefore, it is paramount to consider end-of-life devices as secondary resources for cobalt. Electrorecovery of cobalt from leached solutions has attracted attention due to the sustainability of the recovery process over solvent extraction followed by chemical precipitation. Recently, we reported Co electrorecovery from two different cobalt sources (CoCl2·6H2O and CoSO4·7H2O) using ethylene glycol : choline chloride (EG : ChCl) in a 4.5 : 1 molar ratio, leading to higher purity and easier electrodeposition when sulfate was present as an additive. Here, Co2+ speciation is reported for the two EG : ChCl systems depending on the cobalt source using several spectroscopic techniques (e.g. NMR, EPR, FTIR) in combination with molecular dynamics simulations. Monodentate coordination of SO42- to Co2+, forming the tetrahedral [CoCl3(SO4)]3- was observed as the dominant structure in the system containing CoSO4·7H2O, whereas the system comprising CoCl2·6H2O shows a homoleptic tetrahedral [CoCl4]2- as the dominant structure. This resulted in knowledge being gained regarding Co2+ speciation and the correlation with electrochemistry will contribute to the science required for designing safe electrolytes for efficient electrorecovery.

Temperature-Dependent Dielectric Relaxation Measurements of (Betaine + Urea + Water) Deep Eutectic Solvent in Hz-GHz Frequency Window: Microscopic Insights into Constituent Contributions and Relaxation Mechanisms

A combined experimental and simulation study of dielectric relaxation (DR) of a deep eutectic solvent (DES) composed of betaine, urea, and water with the composition [Betaine:Urea:Water = 11.7:12:1 (weight ratio) and 9:18:5 (molar ratio)] was performed to explore and understand the interaction and dynamics of this system. Temperature-dependent (303 ≤ T/K ≤ 343) measurements were performed over 9 decades of frequency, combining three different measurement setups. Measured DR, comprising four distinct steps with relaxation times spreading over a few picoseconds to several nanoseconds, was found to agree well with simulations. The simulated total DR spectra, upon dissection into three self (intraspecies) and three cross (interspecies) interaction contributions, revealed that the betaine-betaine self-term dominated (∼65%) the relaxation, while the urea-urea and water-water interactions contributed only ∼7% and ∼1%, respectively. The cross-terms (betaine-urea, betaine-water, and urea-water) together accounted for <30% of the total DR. The slowest DR component with a time constant of ∼1-10 ns derived dominant contribution from betaine-betaine interactions, where betaine-water and urea-water interactions also contributed. The subnanosecond (0.1-0.6 ns) time scale originated from all interactions except betaine-water interaction. An extensive interaction of water with betaine and urea severely reduced the average number of water-water H-bonds (∼0.7) and heavily decreased the static dielectric constant of water in this DES (εs ∼ 2). Furthermore, simulated first rank collective single particle reorientational relaxations (C1(t)) and the structural H-bond fluctuation dynamics (CHB (t)) exhibited multiexponential kinetics with time scales that corresponded well with those found both in the simulated and measured DR.

Viscosity effects on the dynamics of diols and diol-based deep eutectic solvents

Diols, characterized by the presence of two hydroxyl groups, form extended hydrogen-bonded networks. Increasing hydrocarbon chain length is known to elevate the viscosity of diols. Given the established influence of viscosity on solvent dynamics, it becomes imperative to comprehend the impact of viscosity on the fluctuation dynamics within diols and establish connections with hydrogen bond formation and breaking dynamics. In this study, we employ two-dimensional infrared spectroscopy to investigate the viscosity dependence of the structural evolution dynamics in three diols with varying chain lengths. Complementing our experimental approach, molecular dynamics simulations are conducted to extract hydrogen bond lifetimes. Our findings reveal a linear correlation between bulk viscosity, solvent fluctuation timescales, and hydrogen bond lifetimes. Notably, the selected diols exhibit the capability to form deep eutectic solvents upon mixing with choline chloride at specific molar ratios. In contrast to molecular solvents like diols, deep eutectic solvents are characterized by the formation of heterogeneous nanodomains, comprising various intercomponent hydrogen-bonded networks. Interestingly, our observations indicate that while the fluctuation dynamics decelerate with increasing bulk viscosity in diol-based deep eutectic solvents, the relationship between viscosity and dynamics is not linear, in contrast to the observed linearity in diols. This nuanced understanding contributes to the broader comprehension of the interplay between viscosity and dynamics in both molecular and deep eutectic solvents.

Molecular dynamics simulations of choline chloride and ascorbic acid deep eutectic solvents: Investigation of structural and dynamics properties

Deep eutectic solvents (DESs) composed of choline chloride (ChCl) and ascorbic acid (AA) were investigated using the molecular dynamics (MD) simulations. The analyses of the configuration, radial distribution function (RDFs), coordination number, spatial distribution function (SDFs), interaction energies, hydrogen bond number, and self-diffusion coefficient of the ChCl/AA binary systems of different concentrations showed that the stability of the hydrogen bond network and the mutual attraction between systems were the strongest at the experimental eutectic concentration (molar ratio of 2:1). In our simulated temperature range from 303.15 to 353.15 K, the hydrogen bonding network of ChCl/AA DES does not undergo considerable alterations, indicating that its stability was insensitive to temperature. In addition, the influence of the water content on the ChCl/AA DES system was further investigated. The simulated results revealed that the water molecules could disrupt the formation of the hydrogen bonding network by occupyin positions that are essential for the formation of hydrogen bonds within the DES system.

Molecular dynamics simulation study of sodium ion structure & dynamics in water in ionic liquids electrolytes using 1-butyl-3-methyl imidazolium tetrafluoroborate and 1-butyl-3-methyl imidazolium hexafluorophosphate

In this paper, we have performed an all-atom molecular dynamics simulation to understand the structure and dynamics of Na+ ions in water mixed Ionic liquids (Water in Ionic liquid). Two ionic liquid (IL) systems consist of (1) 1-butyl-3-methylimidazolium [BMIM] tetrafluoroborate [BF4] and (2) 1-butyl-3-methylimidazolium [BMIM] hexafluorophosphate [PF6] were considered in this work. We understand various inter-molecular structures and dynamic and thermodynamic behaviours of Na+ ions in the water-mixed IL systems. The water (H2O) mole fractions (x) varied from 0.33 to 0.71. The neat ILs [BMIM][BF4] and [BMIM][PF6] pairwise radial distribution functions show a decrease with an increase in x. The [BMIM][PF6] exhibits a strong coordination structure with Na+ ions across the entire range of x values. The rdf between the pairs of Na+-[PF6] presents a significant interaction compared to Na+ and [BF4]. The Na + ions manifested greater coordination with H2O In H2O-[BMIM][PF6] compared to H2O-[BMIM][BF4]. The self-diffusion coefficient (D) values of Na + ions increase with the rise in x in both ILs. The D values of Na + ions are 10-fold higher in [BMIM][BF4] than [BMIM][PF6]. The ionic conductivity values are higher for [BMIM][BF4]. Overall, this paper unveils molecular-level insights for understanding the behavior of Na+ ions in the water in ionic liquid systems.

Importance of the Electrostatic Correlations in Surface Tension of Hydrated Reline Deep Eutectic Solvent from Combined Experiments and Molecular Dynamics Simulations

In this study, the surface tension and the structure of hydrated reline are investigated by using diverse methods. Initially, the surface tension displays a nonlinear pattern as water content increases, decreasing until reaching 45 wt %, then gradually matching that of pure water. This fluctuation is associated with strong electrostatic correlations present in pure reline, which decrease as more water is added. Changes in surface tension reflect a shift from charge layering in pure reline to an increased interfacial hydrogen bonding as the water content rises. This shift causes the segregation of urea molecules into the bulk phase and a gradual anchoring of water molecules to the air-reline interface. An interesting observation is the antisurfactant effect, where heightened interfacial anchoring results in an unexpected increase in real contribution of surface tension. This, along with weakened electrostatic correlations beyond 45 wt % due to reinforced interfacial hydrogen bonding, contributes to the complex behavior of surface tension observed in this study.

Effect of K+ Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion-ion and ion-solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K+ ion. Diffusion-based Nernst-Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

Investigating Biomolecules in Deep Eutectic Solvents with Molecular Dynamics Simulations: Current State, Challenges and Future Perspectives

Deep eutectic solvents (DESs) have recently gained increased attention for their potential in biotechnological applications. DESs are binary mixtures often consisting of a hydrogen bond acceptor and a hydrogen bond donor, which allows for tailoring their properties for particular applications. If produced from sustainable resources, they can provide a greener alternative to many traditional organic solvents for usage in various applications (e.g., as reaction environment, crystallization agent, or storage medium). To navigate this large design space, it is crucial to comprehend the behavior of biomolecules (e.g., enzymes, proteins, cofactors, and DNA) in DESs and the impact of their individual components. Molecular dynamics (MD) simulations offer a powerful tool for understanding thermodynamic and transport processes at the atomic level and offer insights into their fundamental phenomena, which may not be accessible through experiments. While the experimental investigation of DESs for various biotechnological applications is well progressed, a thorough investigation of biomolecules in DESs via MD simulations has only gained popularity in recent years. Within this work, we aim to provide an overview of the current state of modeling biomolecules with MD simulations in DESs and discuss future directions with a focus for optimizing the molecular simulations and increasing our fundamental knowledge.

Development of a QM/MM(ABEEM) method for the deprotonation of neutral and cation radicals in the G-tetrad and GGX(8-oxo-G) tetrad

The rapid deprotonation of G˙+ in the DNA strand impedes positive charge (hole) transfer, whereas the slow deprotonation rate of G˙+ in the G-tetrad makes it a more suitable carrier for hole conduction. The QM/MM(ABEEM) combined method, which involves the integration of QM and the ABEEM polarizable force field (ABEEM PFF), was developed to investigate the deprotonation of neutral and cation free radicals in the G-tetrad and GGX(8-oxo-G) tetrad (xanthine and 8-oxoguanine dual substituted G-tetrad). By incorporating valence-state electronegativity piecewise functions χ*(r) and implementing charge local conservation conditions, QM/MM(ABEEM) possesses the advantage of accurately simulating charge transfer and polarization effect during deprotonation. The activation energy calculated by the QM method of X˙ is the lowest among other bases in the GGX(8-oxo-G) tetrad, which is supported by the computation of the average electronegativity calculated by ABEEM PFF. By utilizing QM/MM(ABEEM) with a two-way free energy perturbation method, the deprotonation activation energy of X˙ in the GGX(8-oxo-G) tetrad is determined to be 33.0 ± 2.1 kJ mol-1, while that of G˙+ in the G-tetrad is 20.7 ± 0.6 kJ mol-1, consistent with the experimental measurement of 20 ± 1.0 kJ mol-1. These results manifest that X˙ in the GGX(8-oxo-G) tetrad exhibits a slower deprotonation rate than G˙+ in the G-tetrad, suggesting that the GGX(8-oxo-G) tetrad may serve as a more favorable hole transport carrier. Furthermore, the unequal average electronegativities of bases in the GGX(8-oxo-G) tetrad impede the deprotonation rate. This study provides a potential foundation for investigating the microscopic mechanism of DNA electronic devices.

Efficient Molecular Dynamics Simulations of Deep Eutectic Solvents with First-Principles Accuracy Using Machine Learning Interatomic Potentials

In recent years, deep eutectic solvents emerged as highly tunable and ecofriendly alternatives to common organic solvents and liquid electrolytes. In the present work, the ability of machine learning (ML) interatomic potentials for molecular dynamics (MD) simulations of these liquids is explored, showcasing a trained neural network potential for a 1:2 ratio mixture of choline chloride and urea (reline). Using the ML potentials trained on density functional theory data, MD simulations for large systems of thousands of atoms and nanosecond-long time scales are feasible at a fraction of the computational cost of the target first-principles simulations. The obtained structural and dynamical properties of reline from MD simulations using our machine learning models are in good agreement with the first-principles MD simulations and experimental results. Running a single MD simulation is highlighted as a general shortcoming of typical first-principles studies if the dynamic properties are investigated. Furthermore, velocity cross-correlation functions are employed to study the collective dynamics of the molecular components in reline.

Nature of a Tetrabutylammonium Chloride-Levulinic Acid Deep Eutectic Solvent

A deep eutectic solvent was formed by considering the mixtures of tetrabutylammonium chloride and levulinic acid, and it is studied via a combined theoretical and experimental approach. Physicochemical properties were measured as a function of temperature, providing a macroscopic characterization of the fluid. Quantum chemistry and classical molecular dynamics simulations were carried out for the nanoscopic characterization, providing attention to the nature, extension, and dynamics of the hydrogen bonding network, which is at the root of the properties of the fluid. The reported study allows multiscale characterization of this fluid as an archetypical example of a natural, low-cost, and sustainable fluid.

BILFF: All-Atom Force Field for Modeling Triazolium- and Benzoate-Based Ionic Liquids

We present an extension of our previously developed all-atom force field BILFF (Bio-polymers in Ionic Liquids Force Field) to three different ionic liquids: 1-ethyl-3-methyl-1,2,3-triazolium acetate ([EMTr][OAc]), 1-ethyl-3-methyl-1,2,3-triazolium benzoate ([EMTr][OBz]), and 1-ethyl-3-methylimidazolium benzoate ([EMIm][OBz]). These ionic liquids are of practical importance as they have the ability to dissolve significant amounts of cellulose even at room temperature. Our force field is optimized to accurately reproduce the strong hydrogen bonding in the system with nearly quantum chemical accuracy. A very good agreement between the microstructure of the quantum chemical simulations over a wide temperature range and experimental density data with the results of BILFF were observed. Non-trivial effects, such as the solvation shell structure and π-π stacking of the cations, are also accurately reproduced. Our force field enables accurate simulations of larger systems, such as solvated cellulose in different (aqueous) ionic liquids, and is the first to present the optimized parameters for mixtures of these solvents and water.

Lysozyme stability in various deep eutectic solvents using molecular dynamics simulations

The ability of neat deep eutectic solvents (DESs) to influence protein structure and function has gained considerable interest due to the unstable nature of enzymes or therapeutic proteins, which are often exposed to thermal, chemical, or mechanical stresses when handled at an industrial scale. In this study, we simulated a model globular protein, lysozyme, in water and six choline chloride-based DES using molecular dynamics simulations, to investigate the structural changes in various solvent environments, giving insights into the overall stability of lysozyme. Root mean square deviation (RMSD) and root mean square fluctuations (RMSF) of the C-α backbone indicated that most DESs induced a less flexible and rigid lysozyme structure compared to water. The radius of gyration and end-to-end distance calculations pointed towards higher structural compactness in reline and levuline, while the structure of lysozyme considerably expanded in oxaline. Protein-solvent interactions were further analysed by hydrogen bonding interactions and radial distribution functions (RDF), which indicated a higher degree of lysozyme-hydrogen bond donor (HBD) interactions compared to lysozyme-choline hydrogen bonding. Surface area analysis revealed an overall % increase in total positive, negative, donor, and acceptor surface areas in malicine and oxaline compared to water and other DESs, indicating the exposure of a larger number of residues to interactions with the solvent. Reline, levuline, and polyol-based DESs comparatively stabilized lysozyme, even though changes in the secondary/tertiary structures were observed.Communicated by Ramaswamy H. Sarma.

Anion dynamics and motional decoupling in a glycerol-choline chloride deep eutectic solvent studied by one- and two-dimensional 35Cl NMR

Chlorine-35 is among the few nuclides that provide an experimental handle on the anion dynamics in choline based deep eutectic solvents. By combining several nuclear magnetic resonance (NMR) techniques, the present work examines the Cl- motions within glyceline, a glycerol : choline chloride 2 : 1 solution, in a large temperature range down to the glass transition temperature Tg. The applied methods include spin relaxometry, second-order line shape analysis, as well as two-dimensional central-transition exchange and stimulated-echo spectroscopy. The finding of unstructured central-transition NMR spectra characterized by a relatively small average quadrupolar coupling attests to a highly disordered, essentially nondirectional anionic coordination in glyceline. For temperatures larger than about 1.3Tg the chlorine motions are well coupled to those of the glycerol and the choline moieties. At lower temperatures the local translational anion dynamics become Arrhenian and increasingly faster than the motion of glyceline's matrix molecules. Upon further cooling, the overall ionic conductivity continues to display a super-Arrhenius behavior, implying that the choline cations rather than the Cl anions dominate the long-range charge transport also near Tg.

Supercritical CO2/Deep Eutectic Solvent Biphasic System as a New Green and Sustainable Solvent System for Different Applications: Insights from Molecular Dynamics Simulations

Deep eutectic solvents (DESs) are one of the most interesting research subjects in green chemistry nowadays. Due to their low toxicity, simple synthesis, and lower prices, they have gradually taken the place of other green solvents such as ionic liquids (ILs) in sustainable processes. However, problems such as high viscosity and high polarity limit the applications of DESs in areas such as extraction, catalysis, and biocatalysis. In this work, we introduce and evaluate the potential application of scCO2/DES for the first time. Molecular dynamics simulations were used to examine the phase behavior, polarity, molecular mobilities, and microstructure of this system. Results show that CO2 molecules can significantly diffuse to the DES phase, while DES components do not appear in the scCO2 phase. The diffused CO2 molecules significantly enhanced the molecular mobility of the DES components. The presence of CO2 molecules changes the DES polarity so that hexane can be solubilized and dispersed in the DES phase. Radial distribution functions show that the solubilized CO2 molecules have negligible effects on the microstructure of DES. It was shown that chloride and urea are the main interaction sites of CO2 in DES. The results of this study show that scCO2/DES as a new class of green and versatile solvents can open a new promising window for research in sustainable chemistry and engineering.

Phenol, the simplest aromatic monohydroxy alcohol, displays a faint Debye-like process when mixed with a nonassociating liquid

Solvated in propylene carbonate, viscous phenol is studied using dielectric spectroscopy and shear rheology. In addition, several oxygen-17 and deuteron nuclear magnetic resonance (NMR) techniques are applied to specifically isotope labeled equimolar mixtures. Quantum chemical calculations are used to check the electrical field gradient at phenol's oxygen site. The chosen combination of NMR methods facilitates the selective examination of potentially hydrogen-bond related contributions as well as those dominated by the structural relaxation. Taken together the present results for phenol in equimolar mixtures with the van der Waals liquid propylene carbonate provide evidence for the existence of a very weak Debye-like process that originates from ringlike supramolecular associates.

A Systematic Study on the Role of Hydrogen Bond Donors in Dictating the Dynamics of Choline-Based Deep Eutectic Solvents

Deep eutectic solvents, promising green alternatives to conventional solvents, consist of a hydrogen bond donor and a hydrogen bond acceptor. The hydrogen bonding components in deep eutectic solvents form an extended hydrogen bonding network, which can be tuned to specific applications by changing the hydrogen bond donors. In this work, we have changed the hydrogen bond donor from a diol to a dicarboxylic acid by systematically replacing a hydroxyl group with an acid group one at a time to investigate the solvation structure and dynamics of the deep eutectic systems. Using a combination of ultrafast vibrational spectroscopy and molecular dynamics simulations, we compared the spectral diffusion and orientational relaxation dynamics of three deep eutectic systems using the vibrational responses of a dissolved anion. Our results indicate that although the solvation structures are marginally different across the systems, distinct differences are present in the solvent fluctuation and solute reorientation dynamics. This work provides a detailed molecular understanding of carboxylic-acid-based deep eutectic systems and how they differ from alcohol-based deep eutectic systems.

Counteraction Effects of Ammonium-Based Ionic Liquids on Urea-Induced Denaturation of α-Lactalbumin: A Comprehensive Molecular Simulation Study

Ionic liquids (ILs) are known to stabilize protein conformations in aqueous medium. Importantly, ILs can also act as refolding additives in urea-driven denaturation of proteins. However, despite the importance of the problem, detailed microscopic understanding of the counteraction effects of ILs on urea-induced protein denaturation remains elusive. In this work, atomistic molecular dynamics (MD) simulations of the protein α-lactalbumin have been carried out in pure aqueous medium, in 8 M binary urea-water solution and in ternary urea-IL-water solutions containing ammonium-based ethyl ammonium acetate (EAA) as the IL at different concentrations (1-4 M). Attempts have been made to quantify detailed molecular-level understanding of the origin behind the counteraction effects of the IL on urea-induced partial unfolding of the protein. The calculations revealed significant conformational changes of the protein with multiple free energy minima due to its partial unfolding in binary urea-water solution. The counteraction effect of the IL was evident from the enhanced structural rigidity of the protein with propensity to transform into a single native free energy minimum state in ternary urea-IL-water solutions. Such an effect has been found to be associated with preferential direct binding of the IL components with the protein and simultaneous expulsion of urea from the interface, thereby providing additional stabilization of the protein in ternary solutions. Most importantly, modified rearrangement of the hydrogen bond network at the interface due to the formation of stronger protein-cation (PC) and protein-anion (PA) hydrogen bonds by breaking relatively weaker protein-urea (PU) and protein-water (PW) hydrogen bonds has been recognized as the microscopic origin behind the counteraction effects of EAA on urea-induced partial unfolding of the protein.

Deep Eutectic Solvents: Properties and Applications in CO2 Separation

Nowadays, many researchers are focused on finding a solution to the problem of global warming. Carbon dioxide is considered to be responsible for the "greenhouse" effect. The largest global emission of industrial CO₂ comes from fossil fuel combustion, which makes power plants the perfect point source targets for immediate CO₂ emission reductions. A state-of-the-art method for capturing carbon dioxide is chemical absorption using an aqueous solution of alkanolamines, most frequently a 30% wt. solution of monoethanolamine (MEA). Unfortunately, the usage of alkanolamines has a number of drawbacks, such as the corrosive nature of the reaction environment, the loss of the solvent due to its volatility, and a high energy demand at the regeneration step. These problems have driven the search for alternatives to that method, and deep eutectic solvents (DESs) might be a very good substitute. Many types of DESs have thus far been investigated for efficient CO₂ capture, and various hydrogen bond donors and acceptors have been used. Deep eutectic solvents that are capable of absorbing carbon dioxide physically and chemically have been reported. Strategies for further CO₂ absorption improvement, such as the addition of water, other co-solvents, or metal salts, have been proposed. Within this review, the physical properties of DESs are presented, and their effects on CO₂ absorption capacity are discussed in conjunction with the types of HBAs and HBDs and their molar ratios. The practical issues of using DESs for CO₂ separation are also described.

Carbon dioxide solubility in choline chloride-based deep eutectic solvents under diverse conditions

CONTEXT

Global warming is a severe problem experiencing the climate crisis due to rising CO₂ emissions. Deep eutectic solvents (DESs) have recently attracted a lot of attention as potential absorbents to mitigate carbon dioxide CO₂ emissions because of their large CO₂ capacities and stability under diverse conditions. Designing a potent DES requires knowledge of molecular-level understanding including structure, dynamics, and interfacial properties in DESs. In this study, we investigate the CO₂ sorption and diffusion in different DESs at different temperatures and pressure using molecular dynamics (MD) simulations. Our results demonstrate that CO₂ molecules preferentially concentrate at the CO₂-DES interface, and the diffusion of CO₂ in bulk DESs increases with increasing pressure and temperature. The solubility of CO₂ in the three DESs increases as ChCL-ethylene glycol < ChCL-urea < ChCL-glycerol at high pressure (58.6 bar).

METHODS

The initial configuration for MD simulations included DES and CO₂ and produced the solvation box using PACKMOL software. The geometries are optimized in the Gaussian 09 software at the theoretical level of B3LYP/6-311 + G*. The partial atomic charges were fitted to an electrostatic surface potential using the CHELPG method. MD simulations were carried out by using the NAMD version 2.13 software. VMD software was used to take the snapshots. TRAVIS software is used to determine spatial distribution functions.

Uncovering the Critical Factors that Enable Extractive Desulfurization of Fuels in Ionic Liquids and Deep Eutectic Solvents from Simulations

Environmental regulatory agencies have implemented stringent restrictions on the permissible levels of sulfur compounds in fuel to reduce harmful emissions and improve air quality. Problematically, traditional desulfurization methods have shown low effectiveness in the removal of refractory sulfur compounds, e.g., thiophene (TS), dibenzothiophene (DBT), and 4-methyldibenzothiophene (MDBT). In this work, molecular dynamics (MD) simulations and free energy perturbation (FEP) have been applied to investigate the use of ionic liquids (ILs) and deep eutectic solvents (DESs) as efficient TS/DBT/MDBT extractants. For the IL simulations, the selected cation was 1-butyl-3-methylimidazolium [BMIM] and the anions included chloride [Cl], thiocyanate [SCN], tetrafluoroborate [BF₄], hexafluorophosphate [PF₆], and bis(trifluoromethylsulfonyl)amide [NTf₂]. The DESs were composed of choline chloride with ethylene glycol (CCEtg) or with glycerol (CCGly). Calculation of excess chemical potentials predicted the ILs to be more promising extractants with energies lower by 1-3 kcal/mol compared to DESs. Increasing IL anion size was positively correlated to enhanced solvation of S-compounds, which was influenced by energetically dominant solute-anion interactions and favorable solute-[BMIM] π-π stacking. For the DESs, the solvent components offered a range of synergistic, yet comparatively weaker, electrostatic interactions that included hydrogen bonding and cation-π interactions. An in-depth analysis of the structure of IL and DES systems is presented, along with a discussion of the critical factors behind experimental trends of S-compound extraction efficiency.

Liquid Structure Scenario of the Archetypal Supramolecular Deep Eutectic Solvent: Heptakis(2,6-di-O-methyl)-β-cyclodextrin/levulinic Acid

The concept of supramolecular solvents has been recently introduced, and the extended liquid-state window accessible for mixtures of functionalized cyclodextrins (CDs) with hydrogen bond (HB) donor species, e.g., levulinic acid, led to the debut of supramolecular deep eutectic solvents (SUPRA-DES). These solvents retain CD's inclusion ability and complement it with enhanced solvation effectiveness due to an extended HB network. However, so far, these promising features were not rationalized in terms of a microscopic description, thus hindering a more complete capitalization. This is the first joint experimental and computational study on the archetypal SUPRA-DES: heptakis(2,6-di-O-methyl)-β-CD/levulinic acid (1:27). We used X-ray scattering to probe CD's aggregation level and molecular dynamics simulation to determine the nature of interactions between SUPRA-DES components. We discover that CDs are homogeneously distributed in bulk and that HB interactions, together with the electrostatic ones, play a major role in determining mutual interaction between components. However, dispersive forces act in synergy with HB to accomplish a fundamental task in hindering hydrophobic interactions between neighbor CDs and maintaining the system homogeneity. The mechanism of mutual solvation of CD and levulinic acid is fully described, providing fundamental indications on how to extend the spectrum of SUPRA-DES combinations. Overall, this study provides the key to interpreting structural organization and solvation tunability in SUPRA-DES to extend the range of sustainable applications for these new, unique solvents.

Does variation in composition affect dynamics when approaching the eutectic composition?

Deep eutectic solvent is a mixture of two or more components, mixed in a certain molar ratio, such that the mixture melts at a temperature lower than individual substances. In this work, we have used a combination of ultrafast vibrational spectroscopy and molecular dynamics simulations to investigate the microscopic structure and dynamics of a deep eutectic solvent (1:2 choline chloride: ethylene glycol) at and around the eutectic composition. In particular, we have compared the spectral diffusion and orientational relaxation dynamics of these systems with varying compositions. Our results show that although the time-averaged solvent structures around a dissolved solute are comparable across compositions, both the solvent fluctuations and solute reorientation dynamics show distinct differences. We show that these subtle changes in solute and solvent dynamics with changing compositions arise from the variations in the fluctuations of the different intercomponent hydrogen bonds.

MD simulations explain the excess molar enthalpies in pseudo-binary mixtures of a choline chloride-based deep eutectic solvent with water or methanol

The addition of molecular liquid cosolvents to choline chloride (ChCl)-based deep eutectic solvents (DESs) is increasingly investigated for reducing the inherently high bulk viscosities of the latter, which represent a major obstacle for potential industrial applications. The molar enthalpy of mixing, often referred to as excess molar enthalpy HE—a property reflecting changes in intermolecular interactions upon mixing—of the well-known ChCl/ethylene glycol (1:2 molar ratio) DES mixed with either water or methanol was recently found to be of opposite sign at 308.15 K: Mixing of the DES with water is strongly exothermic, while methanol mixtures are endothermic over the entire mixture composition range. Knowledge of molecular-level liquid structural changes in the DES following cosolvent addition is expected to be important when selecting such “pseudo-binary” mixtures for specific applications, e.g., solvents. With the aim of understanding the reason for the different behavior of selected DES/water or methanol mixtures, we performed classical MD computer simulations to study the changes in intermolecular interactions thought to be responsible for the observed HE sign difference. Excess molar enthalpies computed from our simulations reproduce, for the first time, the experimental sign difference and composition dependence of the property. We performed a structural analysis of simulation configurations, revealing an intriguing difference in the interaction modes of the two cosolvents with the DES chloride anion: water molecules insert between neighboring chloride anions, forming ionic hydrogen-bonded bridges that draw the anions closer, whereas dilution of the DES with methanol results in increased interionic separation. Moreover, the simulated DES/water mixtures were found to contain extended hydrogen-bonded structures containing water-bridged chloride pair arrangements, the presence of which may have important implications for solvent applications.

Does Viscosity Drive the Dynamics in an Alcohol-Based Deep Eutectic Solvent?

Deep eutectic solvents, consisting of heterogeneous nanodomains of hydrogen-bonded networks, have evolved into a range of viscous fluids that find applications in several fields. As viscosity is known to influence the dynamics of other neoteric solvents like ionic liquids, understanding the effect of viscosity on deep eutectic solvents is crucial to realize their full potential. Herein, we combine polarization-selective pump-probe spectroscopy, two-dimensional infrared spectroscopy, and molecular dynamics simulations to elucidate the impact of viscosity on the dynamics of an alcohol-based deep eutectic solvent, ethaline. We compare the solvent fluctuation and solute reorientation time scales in ethaline with those in ethylene glycol, an ethaline constituent. Interestingly, we find that the solute's reorientation apparently scales the bulk viscosity of the solvent, but the same is not valid for the overall solvation dynamics. Using the variations in the estimated intercomponent hydrogen bond lifetimes, we show that a dissolved solute does not sense the bulk viscosity of the deep eutectic solvent; instead, it senses domain-specific local viscosity determined by the making and breaking of the hydrogen bond network. Our results indicate that understanding the domain-specific local environment experienced by the dissolved solute is of utmost importance in deep eutectic solvents.

Predicting the Surface Tension of Deep Eutectic Solvents Using Artificial Neural Networks

Studies on deep eutectic solvents (DESs), a new class of "green" solvents, are attracting increasing attention from researchers, as evidenced by the rapidly growing number of publications in the literature. One of the main advantages of DESs is that they are tailor-made solvents, and therefore, the number of potential DESs is extremely large. It is essential to have computational methods capable of predicting the physicochemical properties of DESs, which are needed in many industrial applications and research. Surface tension is one of the most important properties required in many applications. In this work, we report a relatively generalized artificial neural network (ANN) for predicting the surface tension of DESs. The database used can be considered comprehensive because it contains 1571 data points from 133 different DES mixtures in 520 compositions prepared from 18 ions and 63 hydrogen bond donors in a temperature range of 277-425 K. The ANN model uses molecular parameter inputs derived from the conductor-like screening model for real solvents (S _σ-profiles). The training and testing results show that the best performing ANN architecture consisted of two hidden layers with 15 neurons each (9-15-15-1). The proposed ANN was excellent in predicting the surface tension of DESs, as R ² values of 0.986 and 0.977 were obtained for training and testing, respectively, with an overall average absolute relative deviation of 2.20%. The proposed models represent an initiative to promote the development of robust models capable of predicting the properties of DESs based only on molecular parameters, leading to savings in investigation time and resources.

Molecular Dynamic Insights into the Distinct Solvation Structures of Aromatic and Aliphatic Compounds in Monoethanolamine-Based Deep Eutectic Solvents

Deep eutectic solvents (DESs) are developing as an alternate medium for aromatic extraction, especially benzene and thiophene from aliphatic hydrocarbon mixtures. In this work, molecular dynamics (MD) simulations were first used to investigate the solvation structure of benzene, thiophene, and n-hexane in monoethanolamine-based DESs. It reveals the liquid structures in the adjacent neighbor shells, which is a function of electron-withdrawing sulfur attached to thiophene and the π-electron cloud of benzene. The intermolecular forces between aromatic, aliphatic, and DES components are analyzed in van der Waals and hydrogen bond interactions. The chloride ions serve as a charge carrier bridge between choline and monoethanolamine precursors. The solvation of benzene, thiophene, and n-hexane in the DESs depends on volume expansion and minor solvent structural changes. Density functional theory results provided information on the mechanism of short-range interactions between organic solutes and studied DES. It aids in understanding the structural orientations of a DES with the addition of solutes, essential to the formation of DES. The solvation shell structure and characteristics were investigated in tandem with the possibility of benzene and thiophene clustering. The ¹H NMR and 2D ¹H-¹H-NOESY were used to investigate the intermolecular interactions between benzene, thiophene, and n-hexane with monoethanolamine-based solvents. It concludes that high-ordered DES1 is more inclined to higher solubility than lower-ordered ones with a higher molar ratio of monoethanolamine. The solvation was reduced because the entropy gain was not maximized in the lower ordered DESs.

Deuteron magnetic resonance study of glyceline deep eutectic solvents: Selective detection of choline and glycerol dynamics

Glyceline, a green solvent considered for various electrochemical applications, represents a multi-component glass former. Viewed from this perspective, the choline cation and the hydrogen bond donor glycerol, the two major constituents forming this deep eutectic solvent, were studied using nuclear magnetic resonance in a selective manner by means of suitably deuteron-labeled isotopologues. Carried out from far above to far below the glass transition temperature, measurements and analyses of the spin-lattice and spin-spin relaxation times reveal that the reorientational dynamics of the components, i.e., of glycerol as well as of chain deuterated choline chloride are slightly different. Possible implications of this finding regarding the hydrogen-bonding pattern in glyceline are discussed. Furthermore, the deuterated methyl groups in choline chloride are exploited as sensitive probes of glyceline's supercooled and glassy states. Apart from spin relaxometry, a detailed line shape analysis of the CD₃ spectra yields valuable insights into the broad intermolecular and intramolecular energy barrier distributions present in this binary mixture.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

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.