OPLS Force Field for Choline Chloride-Based Deep Eutectic Solvents

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Nature of a Tetrabutylammonium Chloride-Levulinic Acid Deep Eutectic Solvent

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

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

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

Lysozyme stability in various deep eutectic solvents using molecular dynamics simulations

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

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

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

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

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

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

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

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

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

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

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

Deep Eutectic Solvents: Properties and Applications in CO2 Separation

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

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

CONTEXT

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

METHODS

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Solubilities and Transport Properties of CO2, Oxalic Acid, and Formic Acid in Mixed Solvents Composed of Deep Eutectic Solvents, Methanol, and Propylene Carbonate

Recently, deep eutectic solvents (DES) have been considered as possible electrolytes for the electrochemical reduction of CO₂ to value-added products such as formic and oxalic acids. The applicability of pure DES as electrolytes is hindered by high viscosities. Mixtures of DES with organic solvents can be a promising way of designing superior electrolytes by exploiting the advantages of each solvent type. In this study, densities, viscosities, diffusivities, and ionic conductivities of mixed solvents comprising DES (i.e., reline and ethaline), methanol, and propylene carbonate were computed using molecular simulations. To provide a quantitative assessment of the affinity and mass transport of CO₂ and oxalic and formic acids in the mixed solvents, the solubilities and self-diffusivities of these solutes were also computed. Our results show that the addition of DES to the organic solvents enhances the solubilities of oxalic and formic acids, while the solubility of CO₂ in the ethaline-containing mixtures are in the same order of magnitude with the respective pure organic components. A monotonic increase in the densities and viscosities of the mixed solvents is observed as the mole fraction of DES in the mixture increases, with the exception of the density of ethaline-propylene carbonate which shows the opposite behavior due to the high viscosity of the pure organic component. The self-diffusivities of all species in the mixtures significantly decrease as the mole fraction of DES approaches unity. Similarly, the self-diffusivities of the dissolved CO₂ and the oxalic and formic acids also decrease by at least 1 order of magnitude as the composition of the mixture shifts from the pure organic component to pure DES. The computed ionic conductivities of all mixed solvents show a maximum value for mole fractions of DES in the range from 0.2 to 0.6 and decrease as more DES is added to the mixtures. Since for most mixtures studied here no prior experimental measurements exist, our findings can serve as a first data set based on which further investigation of DES-containing electrolyte solutions can be performed for the electrochemical reduction of CO₂ to useful chemicals.

Ion Correlation in Choline Chloride-Urea Deep Eutectic Solvent (Reline) from Polarizable Molecular Dynamics Simulations

In recent years, deep eutectic solvents (DESs) emerged as highly tunable and environmentally friendly alternatives to common ionic liquids and organic solvents. In this work, a polarizable model based on the CHARMM Drude polarizable force field is developed for a 1:2 ratio mixture of choline chloride/urea (reline) DES. To successfully reproduce the structure of the liquid as compared to first-principles molecular dynamics simulations, a damping factor was introduced to correct the observed over-binding between the chloride and the hydrogen bonding site of choline. Investigated radial distributions reveal the formation of hydrogen bonds between all the constituents of reline and similar interactions for chloride and urea's oxygen atoms, which could contribute to the melting point depression of the mixture. Predicted dynamic properties from our polarizable force field are in good agreement with experiments, showing significant improvements over nonpolarizable models. Similar to some ionic liquids, an oscillatory behavior in the velocity autocorrelation function of the anion is visible, which can be interpreted as a rattling motion of the lighter anion surrounded by the heavier cations. The obtained results for ionic conductivity of reline show some degree of correlated ion motion in this DES. However, a joint diffusion of ion pairs cannot be observed during the simulations.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

Evolution of microscopic heterogeneity and dynamics in choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are an emerging class of non-aqueous solvents that are potentially scalable, easy to prepare and functionalize for many applications ranging from biomass processing to energy storage technologies. Predictive understanding of the fundamental correlations between local structure and macroscopic properties is needed to exploit the large design space and tunability of DESs for specific applications. Here, we employ a range of computational and experimental techniques that span length-scales from molecular to macroscopic and timescales from picoseconds to seconds to study the evolution of structure and dynamics in model DESs, namely Glyceline and Ethaline, starting from the parent compounds. We show that systematic addition of choline chloride leads to microscopic heterogeneities that alter the primary structural relaxation in glycerol and ethylene glycol and result in new dynamic modes that are strongly correlated to the macroscopic properties of the DES formed.

Tailoring the macroscopic properties of deep eutectic solvents requires knowing how these depend on the local structure and microscopic dynamics. The authors, with computational and experimental tools spanning a wide range of space- and timescales, shed light into the relationship between micro and macroscopic properties in glyceline and ethaline.

Insights on choline chloride-based deep eutectic solvent (reline) + primary alcohol mixtures: a molecular dynamics simulation study

Deep eutectic solvents (DESs) emerged as green solvents for new generation technologies owing to their high chemical and thermal stability. Addition of restricted amount of organic solvents into the DESs plays a significant role in the improvement of thermodynamic and the transport properties to work as a potential solvent in process industries. In this paper, molecular dynamics (MD) simulations were performed to understand the thermophysical and transport properties of choline chloride-based DES (reline) and primary alcohol (methanol and ethanol) mixture in relation to microscopic structure. Density, radial distribution function, coordination number, average number of H-bond, diffusion coefficient and spatial distribution function was calculated in order to understand the structure and involvement of H-bond network at an atomic level. H-bond and spatial distribution function analyses revealed that the chloride ion prefers to be spatially distributed around hydroxyl group of alcohol and found to be more pronounced upon increase in alcohol concentration. As a consequence, it was observed that the H-bonds between Cl^- and urea decreases overall with the loading of alcohol and effect is more pronounced beyond a concentration of 0.4. Self-diffusion values for choline, Cl^- and urea were found to be increased significantly upon increase in concentration of alcohol from 0.6 to 0.8. Overall, our simulation points to the interplay and interactions between the chloride ions and the solvents in determining the structural and transport properties of choline chloride-based DES.

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.