OPLS Force Field for Choline Chloride-Based Deep Eutectic Solvents

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.

Two Faces of the Two-Phase Thermodynamic Model

The quantum harmonic model and the two-phase thermodynamic method (2PT) are widely used to obtain quantum-corrected properties such as isobaric heat capacities or molar entropies. 2PT heat capacities were calculated inconsistently in the literature. For water, the classical heat capacity was also considered, but for organic liquids, it was omitted. We reanalyzed the performance of different quantum corrections on the heat capacities of common organic solvents against experimental data. We have pointed out serious flaws in previous 2PT studies. The vibrational density of states was calculated incorrectly causing a 39% relative error in diffusion coefficients and 45% error in the 2PT heat capacities. The wrong conversion of isobaric and isochoric heat capacities also caused about 40% error but in the other direction. We have introduced the concept of anharmonic correction (AC), which is simply the deviation of the classical heat capacity from that of the harmonic oscillator model. This anharmonic contribution is around +30 to 40 J/(mol K) for water depending on the water model and −8 to −10 J/(mol K) for hydrocarbons and halocarbons. AC is unrealistically large, +40 J/(K mol) for alcohols and amines, indicating some deficiency of the OPLS force field. The accuracy of the computations was also assessed with the determination of the self-diffusion coefficients.

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.

Effects of water on the solvation and structure of lipase in deep eutectic solvents containing a protein destabilizer and stabilizer

Aqueous deep eutectic solvent (DES) solutions emerge as new media for biocatalysis. The large number of DESs provides a space for designing solutions with desired features. One challenge for this design is to understand the fundamental relationship between the water effect on biocatalysis and the DES compositions. We investigate the solvation and structure of a lipase protein in two DESs containing a protein destabilizer (choline : urea (1 : 2)) and stabilizer (choline : glycerol (1 : 2)) and their 1 : 1 aqueous solution using molecular dynamics simulations. The lipase protein in the pure aqueous solution is simulated as the reference. The lipase protein remains folded in both DESs and their aqueous solutions. In both DESs, water molecules weaken the solvation shell of the lipase protein by reducing the protein-DES hydrogen bond lifetimes. However, the water molecules change the surface area and conformation of the active site on the lipase protein differently in the two DESs. Our simulations indicate that the impact on active sites plays an important role in differentiating the effect of water on biocatalysis in aqueous DESs.

Vapor pressures and vapor phase compositions of choline chloride urea and choline chloride ethylene glycol deep eutectic solvents from molecular simulation

Despite the widespread acknowledgment that deep eutectic solvents (DESs) have negligible vapor pressures, very few studies in which the vapor pressures of these solvents are measured or computed are available. Similarly, the vapor phase composition is known for only a few DESs. In this study, for the first time, the vapor pressures and vapor phase compositions of choline chloride urea (ChClU) and choline chloride ethylene glycol (ChClEg) DESs are computed using Monte Carlo simulations. The partial pressures of the DES components were obtained from liquid and vapor phase excess Gibbs energies, computed using thermodynamic integration. The enthalpies of vaporization were computed from the obtained vapor pressures, and the results were in reasonable agreement with the few available experimental data in the literature. It was found that the vapor phases of both DESs were dominated by the most volatile component (hydrogen bond donor, HBD, i.e., urea or ethylene glycol), i.e., 100% HBD in ChClEg and 88%-93% HBD in ChClU. Higher vapor pressures were observed for ChClEg compared to ChClU due to the higher volatility of ethylene glycol compared to urea. The influence of the liquid composition of the DESs on the computed properties was studied by considering different mole fractions (i.e., 0.6, 0.67, and 0.75) of the HBD. Except for the partial pressure of ethylene glycol in ChClEg, all the computed partial pressures and enthalpies of vaporization showed insensitivity toward the liquid composition. The activity coefficient of ethylene glycol in ChClEg was computed at different liquid phase mole fractions, showing negative deviations from Raoult's law.

Manipulating meso-scale solvent structure from Pd nanoparticle deposits in deep eutectic solvents

Deep Eutectic Solvents (DESs) are complex solutions that present unique challenges compared to traditional solvents. Unlike most aqueous electrolytes and ionic liquids, DESs have delicate hydrogen bond networks that are responsible for their highly sensitive compositional dependence on the melting point. Prior work has demonstrated a unique nanoscale structure both experimentally and theoretically that brings both challenges and opportunities to their adoption in traditional electrochemical processes. In this study, we use in situ sample-rotated ultra-small angle x-ray scattering to resolve the near-interface solvent structure after electrodepositing Pd nanoparticles onto a glassy carbon electrode in choline chloride:urea and choline chloride:ethylene glycol DESs. Our results indicate that a hierarchical solvent structure can be observed on the meso-scale in the choline chloride:urea and choline chloride:ethylene glycol systems. Importantly, this extended solvent structure increases between -0.3 V and -0.5 V (vs Ag/AgCl) and remains high until -0.9 V (vs Ag/AgCl). Experimentally, the nature of this structure is more pronounced in the ethylene glycol system, as evidenced by both the x-ray scattering and the electrochemical impedance spectroscopy. Molecular dynamics simulations and dipolar orientation analysis reveal that chloride delocalization near the Pd interface and long-range interactions between the choline and each hydrogen bond donor (HBD) are very different and qualitatively consistent with the experimental data. These results show how the long-range solvent-deposit interactions can be tuned by changing the HBD in the DES and the applied potential.

Molecular dynamic study of alcohol-based deep eutectic solvents

The applicability of deep eutectic solvents is determined by their physicochemical properties. In turn, the properties of eutectic mixtures are the result of the components' molar ratio and chemical composition. Owing to the relatively low viscosities displayed by alcohol-based deep eutectic solvents (DESs), their application in industry is more appealing. Modeling the composition-property relationships established in polyalcohol-based mixtures is crucial for both understanding and predicting their behavior. In this work, a physicochemical property-structure comparison study is made between four choline chloride polyalcohol-based DESs, namely, ethaline, propeline, propaneline, and glyceline. Physicochemical properties obtained from molecular dynamic simulations are compared to experimental data, whenever possible. The simulations cover the temperature range from 298.15 to 348.15 K. The simulated and literature experimental data are generally in good agreement for all the studied DESs. Structural properties, such as radial and spatial distribution functions, coordination numbers, hydrogen bond donor (HBD)-HBD aggregate formation, and hydrogen bonding are analyzed in detail. The higher prevalence of HBD:HBD and HBD:anion hydrogen bonds is likely to be the major reason for the relatively high density and viscosity of glyceline as well as for lower DES self-diffusions.

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.

Significant Improvement in CO2 Absorption by Deep Eutectic Solvents as Immobilized Sorbents: Computational Analysis

To find an alternative way for improving the efficacy of deep eutectic solvents (DESs) to dissolve carbon dioxide, a computational study of DES systems comprising choline chloride and different hydrogen-bond donors (ethylene glycol and glycerol) immobilized on hydrophobic (graphite) and hydrophilic (titanium dioxide) solid surfaces was performed. This research provides quantitative molecular understanding of the role of the DES thickness and also the type of solid support in CO₂ sorption and diffusion using molecular dynamics simulations. In general, the proposed model based on supported DESs immobilized on different supports was developed to correlate the solubility of CO₂ in DESs based on choline chloride. The simulated systems illustrate that CO₂ molecules mainly accumulate at the gas/DES interface in short times, whereas diffusion of CO₂ to the bulk DESs is slower as the thickness of the immobilized DES increases. In addition, the CO₂ absorption capacity of both DESs coated on the TiO₂ surface is larger than that on the graphite surface. Structural and dynamic characteristics were determined using density profiles, distribution functions, orientational analysis, and mean-square displacements. We further demonstrate the effective interaction parameters associated with CO₂ capture by DESs via density functional theory.

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.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

Effect of temperature on concentrated electrolytes for advanced lithium ion batteries

Salt-concentrated electrolytes are emerging as promising electrolytes for advanced lithium ion batteries (LIBs) that can offer high energy density and improved cycle life. To further improve these electrolytes, it is essential to understand their inherent behavior at various operating conditions of LIBs. Molecular dynamics (MD) simulations are extensively used to study various properties of electrolytes and explain the associated molecular-level phenomena. In this study, we use classical MD simulations to probe the properties of the concentrated electrolyte solution of 3 mol/kg lithium hexafluorophosphate (LiPF₆) salt in the propylene carbonate solvent at various temperatures ranging from 298 to 378 K. Our results reveal that the properties such as ionic diffusivity and molar conductivity of a concentrated electrolyte are more sensitive to temperature compared to that of dilute electrolytes. The residence time analysis shows that temperature affects the Li⁺ ion solvation shell dynamics significantly. The effect of temperature on the transport and dynamic properties needs to be accounted carefully while designing better thermal management systems for batteries made with concentrated electrolytes to garner the advantages of these electrolytes.

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.

Thermodynamic, transport, and structural properties of hydrophobic deep eutectic solvents composed of tetraalkylammonium chloride and decanoic acid

With the emergence of hydrophobic deep eutectic solvents (DESs), the scope of applications of DESs has been expanded to include situations in which miscibility with water is undesirable. Whereas most studies have focused on the applications of hydrophobic DESs from a practical standpoint, few theoretical works exist that investigate the structural and thermodynamic properties at the nanoscale. In this study, Molecular Dynamics (MD) simulations have been performed to model DESs composed of tetraalkylammonium chloride hydrogen bond acceptor and decanoic acid hydrogen bond donor (HBD) at a molar ratio of 1:2, with three different cation chain lengths (4, 7, and 8). After fine-tuning force field parameters, densities, viscosities, self-diffusivities, and ionic conductivities of the DESs were computed over a wide temperature range. The liquid structure was examined using radial distribution functions (RDFs) and hydrogen bond analysis. The MD simulations reproduced the experimental density and viscosity data from the literature reasonably well and were used to predict diffusivities and ionic conductivities, for which experimental data are scarce or unavailable. It was found that although an increase in the cation chain length considerably affected the density and transport properties of the DESs (i.e., yielding smaller densities and slower dynamics), no significant influence was observed on the RDFs and the hydrogen bonds. The self-diffusivities showed the following order for the mobility of the various components: HBD > anion > cation. Strong hydrogen bonds between the hydroxyl and carbonyl groups of decanoic acid and between the hydroxyl group of decanoic acid and chloride were observed to dominate the intermolecular interactions.

Group contribution and atomic contribution models for the prediction of various physical properties of deep eutectic solvents

The urgency of advancing green chemistry from labs and computers into the industries is well-known. The Deep Eutectic Solvents (DESs) are a promising category of novel green solvents which simultaneously have the best advantages of liquids and solids. Furthermore, they can be designed or engineered to have the characteristics desired for a given application. However, since they are rather new, there are no general models available to predict the properties of DESs without requiring other properties as input. This is particularly a setback when screening is required for feasibility studies, since a vast number of DESs are envisioned. For the first time, this study presents five group contribution (GC) and five atomic contribution (AC) models for densities, refractive indices, heat capacities, speeds of sound, and surface tensions of DESs. The models, developed using the most up-to-date databank of various types of DESs, simply decompose the molecular structure into a number of predefined groups or atoms. The resulting AARD% of densities, refractive indices, heat capacities, speeds of sound and surface tensions were, respectively, 1.44, 0.37, 3.26, 1.62, and 7.59% for the GC models, and 2.49, 1.03, 9.93, 4.52 and 7.80% for the AC models. Perhaps, even more importantly for designer solvents, is the predictive capability of the models, which was also shown to be highly reliable. Accordingly, very simple, yet highly accurate models are provided that are global for DESs and needless of any physical property information, making them useful predictive tools for a category of green solvents, which is only starting to show its potentials in green technology.

Using coarse-grained molecular dynamics to understand the effect of ionic liquids on the aggregation of Pluronic copolymer solutions

This study is aimed to enhance the understanding of the interaction between ionic liquids (ILs) and non-ionic Pluronic triblock copolymers in aqueous two-phase micellar systems (ATPMS) used for the selective separation/purification of hydrophobic biomolecules. The ILs allow a precise control of the cloud point phase separation temperature (CPT), particularly important when the stability of the molecule is highly dependent on temperature. The effect of choline-based ILs, with two different counter-anions, chloride and hexanoate, was evaluated using molecular dynamics simulations (MD) for F-68 and L-35 Pluronic aqueous solutions. The simulations revealed the role played by the anions during the Pluronic self-assembly, with choline chloride hindering Pluronic aggregation and the choline hexanoate favouring micelle formation and coalescence, in agreement with the experimental data. A detailed study of the accessible surface area of Pluronic showed a progressive dehydration of the Pluronic hydrophilic micelle corona in choline hexanoate mixtures promoting inter-micelle interactions and, consequently, micelle coalescence. With the addition of choline hexanoate, it was observed that the hydrophilic segments, which form the micelle corona, twisted towards the Pluronic micelle core. The electrostatic interaction is also shown to play a key role in this IL-Pluronic aqueous solution, as the hexanoate anions are accommodated in the Pluronic micelle core, while the choline cations are hosted by the Pluronic micelle corona, with the ions interacting with each other during the self-assembly process. In addition, a comparison study of F-68 and L-35 aqueous solutions shows that the IL impact depends on the length of the Pluronic hydrophilic segment. This work provides a realistic microscopic scenario of the complex interactions between Pluronic copolymers and ILs.

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.

Water in biodegradable glucose–water–urea deep eutectic solvent: modifications of structure and dynamics in a crowded environment

Molecular dynamics simulations have been performed on a highly viscous (η ∼ 255 cP) naturally abundant deep eutectic solvent (NADES) composed of glucose, urea and water in a weight ratio of 6 : 4 : 1 at 328 K.

Deep Eutectic Solvents: A Review of Fundamentals and Applications

Deep eutectic solvents (DESs) are an emerging class of mixtures characterized by significant depressions in melting points compared to those of the neat constituent components. These materials are promising for applications as inexpensive "designer" solvents exhibiting a host of tunable physicochemical properties. A detailed review of the current literature reveals the lack of predictive understanding of the microscopic mechanisms that govern the structure-property relationships in this class of solvents. Complex hydrogen bonding is postulated as the root cause of their melting point depressions and physicochemical properties; to understand these hydrogen bonded networks, it is imperative to study these systems as dynamic entities using both simulations and experiments. This review emphasizes recent research efforts in order to elucidate the next steps needed to develop a fundamental framework needed for a deeper understanding of DESs. It covers recent developments in DES research, frames outstanding scientific questions, and identifies promising research thrusts aligned with the advancement of the field toward predictive models and fundamental understanding of these solvents.

Modeling the Physicochemical Properties of Natural Deep Eutectic Solvents

Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point owing to specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents); hence, they are often considered to be a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts have been made to this end. This Review aims at structuring the present knowledge as an outline for future research. First, the key properties of NADES are reviewed and related to their structure on the basis of the available experimental data. Second, available modeling methods applicable to NADES are reviewed. At the molecular level, DFT and molecular dynamics allow density differences and vibrational spectra to be interpreted, and interaction energies to be computed. Additionally, properties at the level of the bulk medium can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed: models based on group-contribution methods and machine learning. A combination of bulk-medium and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension, and refractive index on the other. Multiscale modeling, combining molecular and macroscale methods, is expected to strongly enhance the predictability of NADES properties and their interaction with solutes, and thus yield truly tailorable solvents to accommodate (bio)chemical reactions.

Deep Eutectic Solvents Mixed with Fluorinated Refrigerants for Absorption Refrigeration: A Molecular Simulation Study

Molecular simulations were performed to evaluate mixtures of fluorinated refrigerants with deep eutectic solvents (DESs), for potential use in single-effect absorption refrigeration cycles that use low quality waste heat sources at temperatures of ∼80 °C. The refrigerants considered were the hydrofluorocarbon R245fa and the hydrofluoroolefins R1234zeE and HFO1336mzzE, whereas the DESs evaluated were 1:2 molar mixtures of choline chloride with either ethylene glycol (ethaline) or levulinic acid (levuline) as hydrogen bond donors (HBDs). Assuming the same cycle operating conditions, the waste heat cycle efficiency η was computed for all working fluid mixtures from molecular simulation results of the mixture densities and Henry's law constants of the refrigerants in the DESs, coupled with phase equilibrium calculations and the enthalpies of the pure refrigerants. The largest efficiency was obtained for the mixture R245fa-ethaline (η = 6.82), followed by R245fa-levuline (η = 4.64) and HFO1336mzzE-levuline (η = 2.10). These modest efficiency values could be further increased by tailoring the cycle operating conditions to each particular refrigerant-DES system, as well as optimizing our choice of working fluid mixtures, neither of which we attempted in this study. Strong interactions were observed between the chlorine anions and some of the hydrogen atoms of the refrigerants, but in general the cation-refrigerant and HBD-refrigerant interactions are weaker compared to the refrigerant-refrigerant interactions. Refrigerant molecules have the largest diffusivities and make the cations, anions and HBD to move faster compared to systems of DESs without refrigerant; in general, species in refrigerant-ethaline mixtures have larger diffusivities compared to those for refrigerant-levuline mixtures. We also computed waste heat cycle efficiencies for the same R134a-DES mixtures studied in our previous work, finding significant differences between the efficiencies determined from molecular simulation data and those determined before using the COSMO-RS approach using two standard parametrizations. This observation suggests that further work is needed to improve the accuracy of the COSMO-RS predictions for these systems.

Modeling Alcohol Dehydrogenase Catalysis in Deep Eutectic Solvent/Water Mixtures

The use of oxidoreductases (EC1) in non-conventional reaction media has been increasingly explored. In particular, deep eutectic solvents (DESs) have emerged as a novel class of solvents. Herein, an in-depth study of bioreduction with an alcohol dehydrogenase (ADH) in the DES glyceline is presented. The activity and stability of ADH in mixtures of glyceline/water with varying water contents were measured. Furthermore, the thermodynamic water activity and viscosity of mixtures of glyceline/water have been determined. For a better understanding of the observations, molecular dynamics simulations were performed to quantify the molecular flexibility, hydration layer, and intraprotein hydrogen bonds of ADH. The behavior of the enzyme in DESs follows the classic dependence of water activity (aW ) in non-conventional media. At low aW values (<0.2), ADH does not show any activity; at higher aW values, the activity was still lower than that in pure water due to the high viscosities of the DES. These findings could be further explained by increased enzyme flexibility with increasing water content.

Characterization of a trans-trans Carbonic Acid-Fluoride Complex by Infrared Action Spectroscopy in Helium Nanodroplets

The high Lewis basicity and small ionic radius of fluoride promote the formation of strong ionic hydrogen bonds in the complexation of fluoride with protic molecules. Herein, we report that carbonic acid, a thermodynamically disfavored species that is challenging to investigate experimentally, forms a complex with fluoride in the gas phase. Intriguingly, this complex is highly stable and is observed in abundance upon nanoelectrospray ionization of an aqueous sodium fluoride solution in the presence of gas-phase carbon dioxide. We characterize the structure and properties of the carbonic acid–fluoride complex, F^–(H₂CO₃), and its deuterated isotopologue, F^–(D₂CO₃), by helium nanodroplet infrared action spectroscopy in the photon energy range of 390–2800 cm^–1. The complex adopts a C_2v symmetry structure with the carbonic acid in a planar trans–trans conformation and both OH groups forming ionic hydrogen bonds with the fluoride. Substantial vibrational anharmonic effects are observed in the infrared spectra, most notably a strong blue shift of the symmetric hydrogen stretching fundamental relative to predictions from the harmonic approximation or vibrational second-order perturbation theory. Ab initio thermostated ring-polymer molecular dynamics simulations indicate that this blue shift originates from strong coupling between the hydrogen stretching and bending vibrations, resulting in an effective weakening of the OH···F^– ionic hydrogen bonds.

Thermophysical properties of glyceline–water mixtures investigated by molecular modelling

Water activity and shear viscosity of water–glyceline mixtures are important process parameters that can be effectively calculated using molecular modelling.

Investigation of glycerol hydrogen-bonding networks in choline chloride/glycerol eutectic-forming liquids using neutron diffraction

Neutron scattering reveals the persistent three-dimensional hydrogen-bonding network between glycerol molecules in the 1 : 2 choline chloride/glycerol eutectic.