Anatomy of Microscopic Structure of Ethaline Deep Eutectic Solvent Decoded through Molecular Dynamics Simulations

Deep eutectic solvents on a tightrope: balancing the entropy and enthalpy of mixing

The large melting point depressions characterising deep eutectic solvents (DESs) are related to negative deviations from ideal mixing behaviour characterised by the excess Gibbs energy. Favourable excess Gibbs energies result from a balance between the excess entropy and enthalpy of mixing, which was experimentally determined for three choline chloride (ChCl) based mixtures using calorimetry. While the excess Gibbs energy of H2O + ChCl is enthalpy dominated, those of ethylene glycol (EG) + ChCl and 1,3-propanediol + ChCl are entropy dominated. Molecular dynamics simulations using polarisable force-fields show intermolecular hydrogen bonds between DES constituents for H2O + ChCl and EG + ChCl. Hence, inter-species hydrogen bonding does not guarantee enthalpy-dominated melting point depressions. We suggest future research to focus on tuning the entropy-enthalpy balance via the chemical nature of the DES constituents.

Exploring Unusual Confined Chemistry: Excited-State Proton Transfer Reaction in Supported Liquid Membrane with Deep Eutectic Solvents

Supported liquid membrane (SLM) incorporating ionic liquids (ILs) or deep eutectic solvents (DESs) offers a promising method for ion and (bio)chemical separations and CO2 capture. However, a molecular understanding of whether chemical reactions occur in these confined media is crucial. We report excited-state proton transfer (ESPT) reaction of a photoacid, HPTS, in various DES-based SLMs (pore size ∼280 nm) using steady-state and time-resolved fluorescence spectroscopy. Our findings reveal that, while the ESPT is unfavorable in bulk DESs, it occurs substantially in SLM-containing DESs. Time-resolved area normalized emission spectra (TRANES) show that ESPT time ranges from 2.6 to 7.5 ns and is greatly affected by changes in DES constituents. The results suggest that HPTS interacts with ordered DES structures formed by long-range interfacial effects in membrane pores, making it suitable for ESPT. Furthermore, it is found that the dynamics of solvent relaxation in confined DESs are significantly slower than in bulk liquids. This observation, together with a large red-edge excitation shift, supports the impact of long-range interfacial effects on the DES structure inside membrane pores. Given the task-specific properties of DESs, the incorporation of these solvents into SLM pores can be a useful strategy for investigating new chemical processes.

Structure, Microheterogeneity, and Transport Properties of Ethaline Decoded by X-ray/Neutron Scattering and MD Simulation

Ethaline, a deep eutectic solvent (DES) composed of choline chloride (ChCl)-ethylene glycol (EG) in a 1:2 molar ratio, is garnering significant interest for its wide potential applications. The nature of liquid formation and the structure of H-bonds within ethaline were investigated by X-ray scattering (XRS), neutron scattering (NS), and MD simulations. The sum of the dissociation energy barriers of Ch-EG (3.31 kJ·mol-1) and EG-Cl (4.28 kJ·mol-1) exceeds that of Ch-Cl (5.97 kJ·mol-1). This results in a more pronounced solvation of ChCl by EG compared to ChCl association, facilitating the solubilization of ChCl crystals by EG to form a DES. A partial radial distribution function (PDF) reveals that Cl- solvation is dominated by the hydroxyl group of EG, while the methylene group dominates Ch+ solvation. The spatial distribution function (SDF) shows that the distribution of EG and Cl- around Ch+ partially overlaps with that of the quaternary ammonium group. However, the center of mass distance of Ch-Cl (4.95 Å) is significantly lower than that of Ch-EG (5.65 Å), suggesting a favorable advantage for Cl- in this competition. Chain and ring structure distributions provide direct evidence of the microheterogeneity of ethaline. Hydroxyl groups on the EG promote the formation of a chain structure in ethaline, while methylene groups favor a ring structure. H-bond, carbon H-bond, and Cl- bridge bond restrict Cl- diffusion. This new understanding is crucial for a deeper comprehension of the microstructure of ethaline and for elucidating its mechanisms in applications.

Endowing the Operability of Supercapacitors at High Temperatures by Regulating the Solvation Structure in Dilute Hybrid Electrolyte with Trimethyl Phosphate Cosolvent

The redox stabilities of different oxygen donor solvents (C═O, P═O and S═O) and lithium salt anions for supercapacitors (SCs) electrolytes have been compared by calculating the frontier molecular orbital energy. Among six lithium difluoro(oxalate)borate (LiDFOB)-based mono-solvent electrolytes, the dilute LiDFOB-1,4-butyrolactone (GBL) electrolyte exhibits the highest operating voltage but suffers from electrolyte breakdown at elevated temperatures. Trimethyl phosphate (TMP) exhibits the highest redox stability and a strongly negative electrostatic potential (ESP), making it suitable for promoting the dissolution of LiDFOB as expected. Therefore, TMP is selected as a co-solvent into LiDFOB-GBL electrolyte to regulate Li+ solvation structure and improve the operability of electrolytes at high temperatures. The electrochemical stable potential window (ESPW) of 0.5 m LiDFOB-G/T(5/5) hybrid electrolyte can reach 5.230 V. The activated carbon (AC)-based symmetric SC using 0.5 m LiDFOB-G/T(5/5) hybrid electrolyte achieves a high energy density of 54.2 Wh kg-1 at 1.35 kW kg-1 and the capacitance retention reaches 89.2% after 10 000 cycles. The operating voltage of SC can be maintained above 2 V when the temperature rises to 60 °C.

Solvation Shell Anatomy of H2S and CO Dissolved in Reline and Ethaline Deep Eutectic Solvents

Rising atmospheric concentrations of anthropogenic hydrogen sulfide (H2S) and carbon monoxide (CO) as a result of industrialization have encouraged researchers to explore innovative technologies for capturing these gases. Deep eutectic solvents (DESs) are an alternative media for mitigating H2S and CO emissions. Herein, we have employed ab initio molecular dynamics simulations to investigate the structures of the nearest-neighbor solvation shells surrounding H2S and CO when they are dissolved in reline and ethaline DESs. We aim to delineate the structural arrangement responsible for favorable H2S and CO capture by analyzing the key interactions between H2S and CO solutes with various components of the DESs. We observe that in the reline-H2S system, chloride and carbonyl oxygen of urea are found to have the closest distance interaction with hydrogen atoms of the H2S solute. The sulfur atom of H2S is found to be predominantly solvated by hydrogen and oxygen atoms of urea molecules and the hydroxyl hydrogen of choline cations. The chloride ions and ethylene glycol molecules predominantly govern the solvation of H2S in the ethaline-H2S system. In both the DESs, H2S is solvated by the hydroxyl group of the choline cations rather than by their ammonium group. In the reline-CO system, all the atoms of urea along with chloride dominate the immediate solvation shell around CO. In the ethaline-CO system, hydroxyl oxygen and hydrogen atoms of ethylene glycol are found in the nearest solvation structure around CO. Both the DESs exhibit a stronger solvent-solute charge-transfer tendency toward the H2S solute compared to CO.

Bulk nanostructure of a deep eutectic solvent with an amphiphilic hydrogen bond donor

Neutron diffraction with empirical potential structure refinement (EPSR) show the deep eutectic solvent (DES) 1 : 4 choline chloride : butyric acid is amphiphilically nanostructured. Nanostructure results from solvophobic interactions between the alkyl chains of the butyric acid hydrogen bond donor (HBD) and is retained with addition of 10 wt% water. EPSR fits to the diffraction data is used to produce a three-dimensional model of the liquid which is interrogated to reveal the interactions leading to the solvophobic effect, and therefore nanostructure, in this DES at atomic resolution. The model shows electrostatic and hydrogen bond interactions cause the cation, anion and HBD acid group to cluster into a polar domain, from which the acid alkyl chains are solvophobically excluded into theapolar domain. The polar and apolar domains percolate through the liquid in a bicontinuous sponge-like structure. The effect of adding 10 wt% water is probed, revealing that water molecules are sequestered around the cation and anion within the polar domain, while the neat bulk structure is retained. Alkyl chain packing in the apolar domain becomes slightly better-defined indicating water marginally strengthens solvophobic segregation. These findings reveal bulk self-assembled nanostructure can be produced in DESs via an amphiphilic HBD.

Strong Replaces Weak: Design of H-Bond Interactions Enables Cryogenic Aqueous Zn Metal Batteries

Despite the numerous advantages of aqueous Zn batteries, their practical application under cryogenic conditions is hindered by the freezing of the electrolyte because the abundance of hydrogen bonds (H-bonds) between H2O molecules drives the aqueous system to transform to an orderly frozen structure. Here, a design of H-bond interactions based on the guiding ideology of "strong replaces weak" is proposed. The strong H-bonds formed between introduced eutectic components and water molecules break down the weak H-bonds in the original water molecule network, which contributes to an ultralow freezing point and a high ionic conductivity of 1.7 mS cm-1 at -40 °C. Based on multiperspective theoretical simulations and tailor-made in situ cooling Raman characterizations, it has been demonstrated that substituting weak H-bonds with strong H-bonds facilitates the structural reshaping of Zn2+ solvation and remodeling of the H-bond network in the electrolyte. Endowed with this advantage, reversible and stable Zn plating/stripping behaviors could be realized at -40 °C, and the full cells display a high discharge capacity (200 mA h g-1) at -40 °C with ∼75% capacity retention after 1000 cycles. This study will expand the design philosophy of antifreezing aqueous electrolytes and provide a perspective to promote the adoption of Zn metal batteries for cryogenic environment large-scale energy storage.

Solvation Shell Structures of Ammonia in Reline and Ethaline Deep Eutectic Solvents

Because of increasing atmospheric anthropogenic ammonia (NH₃) emission, researchers are devising new techniques to capture NH₃. Deep eutectic solvents (DESs) are found as potential media for NH₃ mitigation. In the present study, we have carried out ab initio molecular dynamics (AIMD) simulations to decipher the solvation shell structures of an ammonia solute in reline (1:2 mixture of choline chloride and urea) and ethaline (1:2 mixture of choline chloride and ethylene glycol) DESs. We aim to resolve the fundamental interactions which help stabilize NH₃ in these DESs, focusing on the structural arrangement of the DES species in the nearest solvation shell around NH₃ solute. In reline, the hydrogen atoms of NH₃ are preferentially solvated by chloride anions and the carbonyl oxygen atoms of urea. The nitrogen atom of NH₃ renders hydrogen bonding with hydroxyl hydrogen of the choline cation. The positively charged head groups of the choline cations prefer to stay away from NH₃ solute. In ethaline, strong hydrogen bonding interaction exists between the nitrogen atom of NH₃ and hydroxyl hydrogen atoms of ethylene glycol. The hydrogen atoms of NH₃ are found to be solvated by hydroxyl oxygen atoms of ethylene glycol and choline cation. While ethylene glycol molecules play a crucial role in solvating NH₃, the chloride anions remain passive in deciding the first solvation shell. In both the DESs, choline cations approach NH₃ from their hydroxyl group side. We observe slightly stronger solute-solvent charge transfer and hydrogen bonding interaction in ethaline than those in reline.

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.

The Role of Water Content of Deep Eutectic Solvent Ethaline in the Anodic Process of Gold Electrode

Traditional coupling of ligands for gold wet etching makes large-scale applications problematic. Deep eutectic solvents (DESs) are a new class of environment-friendly solvents, which could possibly overcome the shortcomings. In this work, the effect of water content on the Au anodic process in DES ethaline was investigated by combining linear sweep voltammetry (LSV) and electrochemical impedance spectroscopy (EIS). Meanwhile, we employed atomic force microscopy (AFM) to image the evolution of the surface morphology of the Au electrode during its dissolution and passivation process. The obtained AFM data help to explain the observations about the effect of water content on the Au anodic process from the microscopic perspective. High water contents make the occurrence of anodic dissolution of gold at higher potential, but enhances the rate of the electron transfer and gold dissolution. AFM results reveal the occurrence of massive exfoliation, which confirms that the gold dissolution reaction is more violent in ethaline with higher water contents. In addition, AFM results illustrate that the passive film and its average roughness could be tailored by changing the water content of ethaline.

Ethaline and related systems: may be not “deep” eutectics but clearly interesting ionic liquids

Ethaline, the 1:2 molar ratio mixture of ethylene glycol (EG) and choline chloride (ChCl), is generally regarded as a typical type III deep eutectic solvent (DES). However, careful differential scanning calorimetry (DSC) of EG + ChCl mixtures surprisingly revealed that the liquidus lines of the phase diagram apparently follow the predictions for an ideal binary non-electrolyte mixture. Applying broad-band dielectric relaxation spectroscopy to room-temperature solutions of ChCl, and of the related salts choline iodide and chlorocholine chloride, in EG up to saturation, we explored the possible reasons for this conundrum. It appears that in these solutions free ions are rather scarce. Instead, contact ion pairs and larger aggregates predominate.

Why Does Ethaline Apparently Behave as an Ideal Binary Mixture?

Phase diagram mixtures of the salt choline chloride (ChCl) with ethylene glycol (EG) surprisingly seem to behave as ideal binary nonelectrolyte mixtures [Agieienko, V.; Buchner, R. Phys. Chem. Chem. Phys.2022, 24, 5265]. To shed some light on this conundrum, results of broad-band dielectric relaxation spectroscopy (DRS) and quantum-chemical calculations are reported for solutions of ChCl, choline iodide (ChI), and chlorocholine chloride (ClChCl), in EG up to saturation at 298.15 K. The data revealed that all three solutes are only weakly solvated in the sense that on average per equivalent of solute only one EG OH-group is dynamically affected. While contact ion pairs are significant for solute concentrations of ≲1 M, free cation concentrations are rather low. Instead, over the entire concentration range a large fraction of the dipolar cations could not be detected by DRS. We argue that the latter are embedded in large solute aggregates, explaining thus the phase diagram of ChCl + EG and the very low ionicity of all systems.

Structure of Deep Eutectic Solvents (DESs): What We Know, What We Want to Know, and Why We Need to Know It

Deep eutectic solvents (DESs) are a tunable class of solvents with many advantageous properties including good thermal stability, facile synthesis, low vapor pressure, and low-to-negligible toxicity. DESs are composed of hydrogen bond donors and acceptors that, when combined, significantly decrease the freezing point of the resulting solvent. DESs have distinct interfacial and bulk structural heterogeneity compared to traditional solvents, in part due to various intramolecular and intermolecular interactions. Many of the physiochemical properties observed for DESs are influenced by structure. However, our understanding of the interfacial and bulk structure of DESs is incomplete. To fully exploit these solvents in a range of applications including catalysis, separations, and electrochemistry, a better understanding of DES structure must be obtained. In this Perspective, we provide an overview of the current knowledge of the interfacial and bulk structure of DESs and suggest future research directions to improve our understanding of this important information.

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.

Solvent Organization around Methane Dissolved in Archetypal Reline and Ethaline Deep Eutectic Solvents as Revealed by AIMD Investigation

Because of the rising concentration of harmful greenhouse gases like methane in the atmosphere, researchers are striving for developing novel techniques for capturing these gases. Recently, neoteric liquids such as deep eutectic solvents (DESs) have emerged as an efficient means of sequestration of methane. Herein, we have performed ab initio molecular dynamics (AIMD) simulations to elucidate the solvation structure around a methane molecule dissolved in reline and ethaline DESs. We aim to elicit the structural organization of different constituents of the DESs in the vicinity of methane, particularly highlighting the key interactions that stabilize such gases in DESs. We observe quite different solvation structures of methane in the two DESs. In ethaline, chloride ions play an active role in solvating methane. Instead, in reline, chloride ions do not interact much with the methane molecule in the first solvation shell. In reline, choline cations approach the methane molecule from their hydroxyl group side, whereas urea molecules approach methane from their carbonyl oxygen as well as amide group sides. In ethaline, ethylene glycol and Cl^- dominate the nearest neighbor solvation structure around the methane molecule. In both the DESs, we do not observe any significant methane-DES charge transfer interactions, apart from what is present between choline cation and Cl^- anion.

Studying the Formation of Choline Chloride- and Glucose-Based Natural Deep Eutectic Solvent at the Molecular Level

The liquid waste is the major source of waste, which usually generated from academic laboratories and industry during the extraction, separation, chemical synthesis, and pretreatment processes. These chemical and engineering processes require more solvents. In this regard, there is a need to develop more environmentally friendly, cheaper, non-toxic solvents that are harmless to humans and the environment. In this regard, deep eutectic solvents (DES) and their derivatives so-called natural deep eutectic solvents (NADES) are a new field in the search for green alternative solvents. In our work, the formation of choline chloride-based NADESs using density functional theory (DFT) calculations, and classical all-atom molecular dynamics (MD) simulation was studied in detail using Gaussian09 and Gromacs software's. Next, the ground state geometry optimizations were performed in the gas phase using DFT B3LYP 6-31 + G(d) level of theory. Moreover, classical all-atom MD simulations were implemented using Gromos force field. After the modeling and simulations, the DFT calculation results revealed the formation of NADESs via formation (creation) of binding between chlorine and choline, and chlorine and glucose. At the same time, the results of classical all-atom MD simulations, based on the time average of the equilibrated production run of MD simulations, stated that the nitrogen atom of choline ion and chloride ion has greater interactions, while chloride ion has also greater interaction with glucose during formation of NADES. The outcomes of both DFT and classical all-atom MD simulations are in good agreements.

Structure, Organization, and Heterogeneity of Water-Containing Deep Eutectic Solvents

The spectroscopy and structural dynamics of a deep eutectic mixture (KSCN/acetamide) with varying water content is investigated from 2D IR (with the C-N stretch vibration of the SCN^- anions as the reporter) and THz spectroscopy. Molecular dynamics simulations correctly describe the nontrivial dependence of both spectroscopic signatures depending on water content. For the 2D IR spectra, the MD simulations relate the steep increase in the cross-relaxation rate at high water content to the parallel alignment of packed SCN^- anions. Conversely, the nonlinear increase of the THz absorption with increasing water content is mainly attributed to the formation of larger water clusters. The results demonstrate that a combination of structure-sensitive spectroscopies and molecular dynamics simulations provides molecular-level insights into the emergence of heterogeneity of such mixtures by modulating their composition.

Molecular Dynamics Simulation Study of the Far-Infrared Spectrum of a Deep Eutectic Solvent

Deep eutectic solvents (DESs) are similar to ionic liquids (IL) in terms of physicochemical properties and technical uses. In ILs, far-infrared (FIR) spectroscopy has been utilized to reveal ionic interactions and even to produce a signature of the strengthening of the cation-anion hydrogen bond. However, for the situation of the DES, where the mixing of a salt and a molecular species makes the interplay between multiple intermolecular interactions even more complex, a full investigation of FIR spectra is still absent. In this work, the FIR spectrum of the DES, often referred to as ethaline, which is a 1:2 mixture of choline chloride and ethylene glycol, is calculated using classical molecular dynamics (MD) simulations and compared to experimental data. To explore the induced dipole effect on the computed FIR spectrum, MD simulations were run with both nonpolarizable and polarizable models. The calculation satisfactorily reproduces the position of the peak at ∼110 cm^-1 and the bandwidth seen in the experimental FIR spectrum of ethaline. The MD simulations show that the charge current is the most important contributor to the FIR spectrum, but the cross-correlation between the charge current and dipole reorientation also plays a role in the polarizable model. The dynamics of the chloride-ethylene glycol correlation span a wide frequency range, with a maximum at ∼150 cm^-1, but it participates as a direct mechanism only in the charge current-dipole reorientation cross-term. Anion correlations, whose dynamics are regulated via correlation with both ethylene glycol and choline, make the most significant contribution to the charge current mechanism. The MD simulations were also utilized to investigate the effect on the FIR spectrum of adding water to the DES and switching to a 1:1 composition.

An Overview of Structure and Dynamics Associated with Hydrophobic Deep Eutectic Solvents and Their Applications in Extraction Processes

Recent development of novel water-immiscible green solvents known as hydrophobic deep eutectic solvents (HDESs) has opened the gates for applications requiring media where presence of water is undesirable. Ever since they were prepared, researchers have used HDESs in diverse fields such as extraction processes, CO 2 sequestration, membrane formation, and catalysis. The microstructure and dynamics associated with the species comprising HDESs guide their suitability for specific applications. For example, varying the alkyl tail length of HDES components significantly affects the dynamics of the components and thus helps in tuning the efficiency of extraction processes. The development of HDESs is still in infancy and very few theoretical studies are available in the literature that help in understanding the structure and dynamics of HDESs. This review highlights the recent work focused on the microscopic structure and dynamics of HDESs and their potential applications, particularly in extraction processes. We have also provided a glimpse of how the integration of experiments and computational techniques can help understand the mechanism of extraction processes.

Structure and Dynamic Properties of a Glycerol-Betaine Deep Eutectic Solvent: When Does a DES Become an Aqueous Solution?

Deep eutectic solvents (DESs) are an emerging class of green solvents with a wide spectrum of potential applications whose properties may be further tailored through the addition of water. Here, we study, through molecular dynamics, the influence of water on the properties of a betaine-glycerol-water (B:G:W) DES (1:2:ζ; ζ = 0 to 100), aiming at getting insight into the structural and dynamic crossover between a DES and an aqueous solution. The density, shear viscosity, and diffusion coefficients are found to exhibit a non-linear dependence of ζ, similar to that observed for the solvation layers' composition. Each Gly and Bet are replaced, respectively, by ∼3 and ∼5 water molecules, with the highest rates of depletion being found for Gly around Bet and Gly around Gly. Above ζ = 7 (70 mol %; 29.5 wt %), a major structural transformation occurs, with the complete disruption of the second Bet-Gly solvation layer and the formation of a new second layer at a shorter distance, accompanied by a sudden change in the rate of increase of the components' diffusion. Nonetheless, opposite to other DES, our results indicate a smooth crossover between a DES and an aqueous solution.

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.

Understanding of Bulk and Interfacial Structures Ternary and Binary Deep Eutectic Solvents with a Constant Potential Method: A Molecular Dynamics Study

In the last decade, deep eutectic solvents (DESs) emerge as promising electrolytes in supercapacitors and rechargeable batteries due to their unique properties, wide electrochemical window, low viscosity, and high ionic...

Effect of hydrogen bond donor molecules ethylene glycerol and lactic acid on electrochemical interfaces in choline chloride based-deep eutectic solvents

Choline chloride (ChCl)-based-deep eutectic solvents (DESs) are widely used in electrochemical fields. In this work, the effect of two types of hydrogen bond donor (HBD) molecules, ethylene glycerol and lactic acid (LA), on electrochemical interfaces between the Au electrode and DESs has been investigated by employing voltammetry and electrochemical impedance spectroscopy. The anodic dissolution and passivation behaviors of the Au electrode are revealed in both ethaline and ChCl:LA. In ChCl:LA, the anodic dissolution of Au is slowed down, and the passivation film is relatively dense and stable due to the existence of the carboxyl group in HBD molecule LA. In the double layer region, the lifting and formation of Au(111) surface reconstruction and a disorder-order phase transition of the chloride ion adlayer were observed in the two DESs. Moreover, compared with ethaline, an extra pair of current peaks appears in ChCl-LA possibly due to the adsorption and desorption of LA on the Au(111) surface, which might imply the stronger interaction of LA with the Au electrode in ChCl:LA. HBD LA could even have marked an impact on the disorder-order phase transition of the chloride ion adlayer. The above results provide new insight into the significant effect of HBD molecules on the anodic dissolution and the passivation of the Au electrode and the electrochemical behaviors in the double layer region.

Fine-Tuning the Polarizable CL&Pol Force Field for the Deep Eutectic Solvent Ethaline

Polarizable force fields are gradually becoming a common choice for ionic soft matter, in particular, for molecular dynamics (MD) simulations of ionic liquids (ILs) and deep eutectic solvents (DESs). The CL&Pol force field introduced in 2019 is the first general, transferable, and polarizable force field for MD simulations of different types of DESs. The original formulation contains, however, some problems that appear in simulations of ethaline and may also have a broader impact. First, the originally proposed atomic diameter parameters are unbalanced, resulting in too weak interactions between the chlorides and the hydroxyl groups of the ethylene glycol molecules. This, in turn, causes an artificial phase separation in long simulations. Second, there is an overpolarization of chlorides due to strong induced dipoles that give rise to the presence of peaks and antipeaks at very low q-vector values (2.4 nm^-1) in the partial components of the structure factors. In physical terms, this is equivalent to overestimated spatial nanoscale heterogeneity. To correct these problems, we adjusted the chloride-hydroxyl radial distribution functions against ab initio data and then extended the use of the Tang-Toennis damping function for the chlorides' induced dipoles. These adjustments correct the problems without losing the robustness of the CL&Pol force field. The results were also compared with the nonpolarizable version, the CL&P force field. We expect that the corrections will facilitate reliable use of the CL&Pol force field for other types of DESs.

In Silico Elucidation of Molecular Picture of Water-Choline Chloride Mixture

Choline chloride (ChCl) is a component of several deep eutectic solvents (DESs) having numerous applications. Recent studies have reported manifold promising use of aqueous choline chloride solution as an alternative to DES, where water plays the role of the hydrogen-bond donor. The characteristic physical properties of the DESs and aqueous DES originate from the "inter-" and intraspecies hydrogen-bond network formed by the constituents. However, a detailed molecular-level picture of choline chloride and water mixture is largely lacking in the literature. This motivates us to carry out extensive all-atom molecular dynamics simulations of the ChCl-water mixture of varying compositions. Our analyses clearly show an overall increase in the interspecies association with an increase in ChCl concentration. At higher concentrations, the trimethylammonium groups of choline are stabilized by a nonpolar interaction, whereas the hydroxyl groups preferentially interact with water. Chloride ions are found to be involved in two types of interactions: one where chloride ions intercalate two or more choline cations, and the other one where they are surrounded by five to six water molecules forming solvated chloride ions. However, the relative fractions of these two types of associations depend on the concentration of ChCl in the mixture. Another important structural aspect is the disruption of the hydrogen-bonded water network due to the presence of both choline cations and chloride ions. However, chloride ions participate to partially restore the tetrahedral arrangement of partners around water molecules.

Bulk and interfacial nanostructure and properties in deep eutectic solvents: Current perspectives and future directions

Deep eutectic solvents (DESs) are a tailorable class of solvents that are rapidly gaining scientific and industrial interest. This is because they are distinct from conventional molecular solvents, inherently tuneable via careful selection of constituents, and possess many attractive properties for applications, including catalysis, chemical extraction, reaction media, novel lubricants, materials chemistry, and electrochemistry. DESs are a class of solvents composed solely of hydrogen bond donors and acceptors with a melting point lower than the individual components and are often fluidic at room temperature. A unique feature of DESs is that they possess distinct bulk liquid and interfacial nanostructure, which results from intra- and inter-molecular interactions, including coulomb forces, hydrogen bonding, van der Waals interactions, electrostatics, dispersion forces, and apolar-polar segregation. This nanostructure manifests as preferential spatial arrangements of the different species, and exists over several length scales, from molecular- to nano- and meso-scales. The physicochemical properties of DESs are dictated by structure-property relationships; however, there is a significant gap in our understanding of the underlying factors which govern their solvent properties. This is a major limitation of DES-based technologies, as nanostructure can significantly influence physical properties and thus potential applications. This perspective provides an overview of the current state of knowledge of DES nanostructure, both in the bulk liquid and at solid interfaces. We provide definitions which clearly distinguish DESs as a unique solvent class, rather than a subset of ILs. An appraisal of recent work provides hints towards trends in structure-property relationships, while also highlighting inconsistencies within the literature suggesting new research directions for the field. It is hoped that this review will provide insight into DES nanostructure, their potential applications, and development of a robust framework for systematic investigation moving forward.

Effective and Selective Extraction of Quercetin from Onion (Allium cepa L.) Skin Waste Using Water Dilutions of Acid-Based Deep Eutectic Solvents

Deep Eutectic Solvents (DESs) are experiencing growing interest as substitutes of polluting organic solvents for their low or absent toxicity and volatility. Moreover, they can be formed with natural bioavailable and biodegradable molecules; they are synthesized in absence of hazardous solvents. DESs are, inter alia, successfully used for the extraction/preconcentration of biofunctional molecules from complex vegetal matrices. Onion skin is a highly abundant waste material which represents a reservoir of molecules endowed with valuable biological properties such as quercetin and its glycosylated forms. An efficient extraction of these molecules from dry onion skin from "Dorata di Parma" cultivar was obtained with water dilution of acid-based DESs. Glycolic acid (with betaine 2/1 molar ratio and L-Proline 3/1 molar ratio as counterparts) and of p-toluensulphonic acid (with benzyltrimethylammonium methanesulfonate 1/1 molar ratio)-based DESs exhibited more than 3-fold higher extraction efficiency than methanol (14.79 µg/mL, 18.56 µg/mL, 14.83 µg/mL vs. 5.84 µg/mL, respectively). The extracted quercetin was also recovered efficaciously (81% of recovery) from the original extraction mixture. The proposed extraction protocol revealed to be green, efficacious and selective for the extraction of quercetin from onion skin and it could be useful for the development of other extraction procedures from other biological matrixes.

Structural and dynamical changes observed when transitioning from an ionic liquid to a deep eutectic solvent

The microscopic molecular structure and dynamics of a new deep eutectic solvent (DES) composed of an ionic liquid (1-hexyl-3-methylimidazolium chloride) and an amide (trifluoroacetamide) at various molar ratios were investigated using linear and non-linear infrared spectroscopy with a vibrational probe. The use of the ionic liquid allows us to investigate the changes that the system undergoes with the addition of the amide or, equivalently, the changes from an ionic liquid to a DES. Our studies revealed that the vibrational probe in the DES senses a very similar local environment irrespective of the cation chemical structure. In addition, the amide also appears to perceive the same molecular environment. The concentration dependence studies also showed that the amide changes from being isolated from other amides in the ionic liquid environment to an environment where the amide-amide interactions are favored. In the case of the vibrational probe, the addition of the amide produced significant changes in the slow dynamics associated with the making and breaking of the ionic cages but did not affect the rattling-in-cage motions perceived by it. Furthermore, the concentration dependence of slow dynamics showed two regimes which are linked to the changes in the overall structure of the solution. These observations are interpreted in the context of a nanoscopic heterogeneous environment in the DES which, according to the observed dynamical regimes, appears at very large concentrations of the amide (molar ratio of greater than 1:1) since for lower amide molar ratios, the amide appears to be not segregated from the ionic liquid. This proposed molecular picture is supported by small angle x-ray scattering experiments.

Comparison and Validation of Force Fields for Deep Eutectic Solvents in Combination with Water and Alcohol Dehydrogenase

Deep eutectic solvents (DESs) have become popular as environmental-friendly solvents for biocatalysis. Molecular dynamics (MD) simulations offer an in-depth analysis of enzymes in DESs, but their performance depends on the force field chosen. Here, we present a comprehensive validation of three biomolecular force fields (CHARMM, Amber, and OPLS) for simulations of alcohol dehydrogenase (ADH) in DESs composed of choline chloride and glycerol/ethylene glycol with varying water contents. Different properties (e.g., protein structure and flexibility, solvation layer, and H-bonds) were used for validation. For two properties (viscosity and water activity) also experiments were performed. The viscosity was calculated with the periodic perturbation method, whereby its parameter dependency is disclosed. A modification of Amber was identified as the best-performing model for low water contents, whereas CHARMM outperforms the other models at larger water concentrations. An analysis of ADH's structure and interactions with the DESs revealed similar predictions for Amber and CHARMM.

Multiple evidences of dynamic heterogeneity in hydrophobic deep eutectic solvents

Hydrophobic deep eutectic solvents (HDESs) have gained immense popularity because of their promising applications in extraction processes. Herein, we employ atomistic molecular dynamics simulations to unveil the dynamics of DL-menthol (DLM) based HDESs with hexanoic (C6), octanoic (C8), and decanoic (C10) acids as hydrogen bond donors. The particular focus is on understanding the nature of dynamics with changing acid tail length. For all three HDESs, two modes of hydrogen bond relaxations are observed. We observe longer hydrogen bond lifetimes of the inter-molecular hydrogen bonding interactions between the carbonyl oxygen of the acid and hydroxyl oxygen of menthol with hydroxyl hydrogen of both acids and menthol. We infer strong hydrogen bonding between them compared to that between hydroxyl oxygen of acids and hydroxyl hydrogens of menthol and acids, marked by a faster decay rate and shorter hydrogen bond lifetime. The translational dynamics of the species in the HDES becomes slower with increasing tail length of the organic acid. Slightly enhanced caging is also observed for the HDES with a longer tail length of the acids. The evidence of dynamic heterogeneity in the displacements of the component molecules is observed in all the HDESs. From the values of the α-relaxation time scale, we observe that the molecular displacements become random in a shorter time scale for DLM-C6. The analysis of the self-van Hove function reveals that the overall distance covered by DLM and acid molecules in the respective HDES is more than what is expected from ideal diffusion. As marked by the shorter time scale associated with hole filling, the diffusion of the oxygen atom of menthol and the carbonyl oxygen of acid from one site to the other is fastest for hexanoic acid containing HDES.

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.

Liquid structure and dynamics in the choline acetate:urea 1:2 deep eutectic solvent

We report on the thermodynamic, structural, and dynamic properties of a recently proposed deep eutectic solvent, formed by choline acetate (ChAc) and urea (U) at the stoichiometric ratio 1:2, hereinafter indicated as ChAc:U. Although the crystalline phase melts at 36-38 °C depending on the heating rate, ChAc:U can be easily supercooled at sub-ambient conditions, thus maintaining at the liquid state, with a glass-liquid transition at about -50 °C. Synchrotron high energy x-ray scattering experiments provide the experimental data for supporting a reverse Monte Carlo analysis to extract structural information at the atomistic level. This exploration of the liquid structure of ChAc:U reveals the major role played by hydrogen bonding in determining interspecies correlations: both acetate and urea are strong hydrogen bond acceptor sites, while both choline hydroxyl and urea act as HB donors. All ChAc:U moieties are involved in mutual interactions, with acetate and urea strongly interacting through hydrogen bonding, while choline being mostly involved in van der Waals mediated interactions. Such a structural situation is mirrored by the dynamic evidences obtained by means of ¹H nuclear magnetic resonance techniques, which show how urea and acetate species experience higher translational activation energy than choline, fingerprinting their stronger commitments into the extended hydrogen bonding network established in ChAc:U.

Deep Eutectic Solvents: Molecular Simulations with a First-Principles Polarizable Force Field

The unique properties of deep eutectic solvents make them useful in a variety of applications. In this work we develop a first-principles force field for reline, which is composed of choline chloride and urea in the molar ratio 1:2. We start with the symmetry adapted perturbation theory (SAPT) protocol and then make adjustments to better reproduce the structure and dynamics of the liquid when compared to first-principles molecular dynamics (FPMD) simulations. The resulting force field is in good agreement with experiments in addition to being consistent with the FPMD simulations. The simulations show that primitive molecular clusters are preferentially formed with choline-chloride ionic pairs bound with a hydrogen bond in the hydroxyl group and that urea molecules coordinate the chloride mainly via the trans-H chelating hydrogen bonds. Incorporating polarizability qualitatively influences the radial distributions and lifetimes of hydrogen bonds and affects long-range structural order and dynamics. The polarizable force field predicts a diffusion constant about an order of magnitude larger than the nonpolarizable force field and is therefore less computationally intensive. We hope this study paves the way for studying complex hydrogen-bonding liquids from a first-principles approach.

Partial Charges Optimized by Genetic Algorithms for Deep Eutectic Solvent Simulations

Deep eutectic solvents (DESs) are a class of solvents often composed of ammonium-based chloride salts and a neutral hydrogen bond donor (HBD) at specific ratios. These cost-effective and environmentally friendly solvents have seen significant growth in multiple fields, including organic synthesis, and in materials and extractions because of their desirable properties. In the present work, a new software called genetic algorithm machine learning (GAML) was developed that utilizes a genetic algorithm (GA) approach to facilitate the development of optimized potentials for liquid simulation (OPLS)-based force field (FF) parameters for eight unique DESs based on three ammonium-based salts and five HBDs at multiple salt:HBD ratios. As an initial test of GAML, partial charges were created for 86 conventional solvents based on neutral organic molecules that yielded excellent overall mean absolute deviations (MADs) of 0.021 g/cm³, 0.63 kcal/mol, and 0.20 kcal/mol compared to experimental densities, heats of vaporization (ΔH_vap), and free energies of hydration (ΔG_hyd), respectively. FFs for DESs constructed from ethylammonium, N,N-diethylethanolammonium, and N-ethyl-N,N-dimethylethanolammonium chloride salts were then parameterized using GAML with exceptional agreement achieved at multiple temperatures for experimental densities, surface tensions, and viscosities with MADs of 0.024 g/cm³, 4.2 mN/m, and 5.3 cP, respectively.

On the coupling between ionic conduction and dipolar relaxation in deep eutectic solvents: Influence of hydration and glassy dynamics

We have studied the ionic conductivity and the dipolar reorientational dynamics of aqueous solutions of a prototypical deep eutectic solvent (DES), ethaline, by dielectric spectroscopy in a broad range of frequencies (MHz-Hz) and for temperatures ranging from 128 to 283 K. The fraction of water in the DES was varied systematically to cover different regimes, starting from the pure DES and its water-in-DES mixtures to the diluted electrolyte solutions. Depending on these parameters, different physical states were examined, including low viscosity liquid, supercooled viscous liquid, amorphous solid, and freeze-concentrated solution. Both the ionic conductivity and the reorientational relaxation exhibited characteristic features of glassy dynamics that could be quantified from the deviation from the Arrhenius temperature dependence and non-exponential decay of the relaxation function. A transition occurred between the water-in-DES regime (<40 wt. %), where the dipolar relaxation and ionic conductivity remained inversely proportional to each other, and the DES-in-water regime (>40 wt. %), where a clear rotation-translation decoupling was observed. This suggests that for a low water content, on the timescale covered by this study (∼10^-6 to 1 s), the rotational and transport properties of ethaline aqueous solutions obey classical hydrodynamic scaling despite these systems being presumably spatially microheterogeneous. A fractional scaling is observed in the DES-in-water regime due to the formation of a maximally freeze-concentrated DES aqueous solution coexisting with frozen water domains at sub-ambient temperature.

Microscopic structural features of water in aqueous-reline mixtures of varying compositions

Reline, a mixture of urea and choline chloride in a 2 : 1 molar ratio, is one of the most frequently used deep eutectic solvents. Pure reline and its aqueous solution have large scale industrial use. Owing to the presence of active hydrogen bond formation sites, urea and choline cations can disrupt the hydrogen-bonded network in water. However, a quantitative understanding of the microscopic structural features of water in the presence of reline is still lacking. We carry out extensive all-atom molecular dynamics simulations to elucidate the effect of the gradual addition of co-solvents on the microscopic arrangements of water molecules. We consider four aqueous solutions of reline, between 26.3 and 91.4 wt%. A disruption of the local hydrogen-bonded structure in water is observed upon inclusion of urea and choline chloride. The extent of deviation of the water structure from tetrahedrality is quantified using the tetrahedral order parameter (q_tet). Our analyses show a monotonic increase in the structural disorder as the co-solvents are added. Increase in the q_tet values are observed when highly electro-negative hetero-atoms like nitrogen, oxygen of urea and choline cations are counted as partners of the central water molecules. Further insights are drawn from the characterization of the hydrogen-bonded network in water and we observe the gradual rupturing of water-water hydrogen bonds and their subsequent replacement by the water-urea hydrogen bonds. A negligible contribution from the hydrogen bonds between water and bulky choline cations has also been found. Considering all the constituents as the hydrogen bond partners we calculate the possibility of a successful hydrogen bond formation with a central water molecule. This gives a clear picture of the underlying mechanism of water replacement by urea.