How sensitive are physical properties of choline chloride-urea mixtures to composition changes: Molecular dynamics simulations and Kirkwood-Buff theory

Deep eutectic solvents for the extraction of polyphenols from food plants

Deep Eutectic Solvents (DESs) offer a promising, sustainable alternative for extracting polyphenols from food plants, known for their health benefits. Traditional extraction methods are often costly and involve toxic solvents. This review discusses the basic concepts, preparation techniques, and factors influencing the effective and safe use of DESs in polyphenol extraction. DESs' adaptability allows integration with other green extraction technologies, such as microwave- and ultrasound-assisted extractions, enhancing their efficiency. This adaptability demonstrates the potential of DESs in the sustainable extraction of bioactive compounds. Current research indicates that DESs could play a significant role in the sustainable procurement of these compounds, marking an important advancement in food science research and development. The review underscores DESs as a realistic, eco-friendly alternative in the realm of natural extraction technologies, offering a significant contribution to sustainable practices in food science.

Mutual Diffusivities of Mixtures of Carbon Dioxide and Hydrogen and Their Solubilities in Brine: Insight from Molecular Simulations

H2-CO2 mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H2-CO2 mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H2 due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO2, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO2 content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO2-H2-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO2 and H2 in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO2 in NaCl brine decreased in the ternary system compared to the binary CO2-NaCl brine system, the solubility of H2 in NaCl brine increased less in the ternary system compared to the binary H2-NaCl brine system. The cooperative effect of H2-CO2 enhances the H2 solubility while suppressing the CO2 solubility. The water content in the gas phase was found to be intermediate between H2-NaCl brine and CO2-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO2-H2 mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H2-CO2 mixtures.

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.

Feasibility of replacing proton exchange membranes with pressure-driven membranes in membrane electrochemical reactors for high salinity organic wastewater treatment

Membrane electrochemical reactor (MER) shows superiority to electrochemical oxidation (EO) in high salinity organic wastewater (HSOW) treatment, but requirement of proton exchange membranes (PEM) increases investment and maintenance cost. In this work, the feasibility of using low-cost pressure-driven membranes as the separation membrane in MER system was systematically investigated. Commonly used pressure-driven membranes, including loose membranes such as microfiltration (MF) and ultrafiltration (UF), as well as dense membranes like nanofiltration (NF) and reverse osmosis (RO), were employed in the study. When tested in a contamination-free solution, MF and UF exhibited superior electrochemical performance compared to PEM, with comparable pH regulation capabilities in the short term. When foulant (humic acid, Ca2+ and Mg2+) presented in the feed, UF saved the most energy (43 %) compared to PEM with similar removal rate of UV254 (∼85 %). In practical applications of MER for treating nanofiltration concentrate (NC) of landfill leachate, UF saved 27 % energy compared to PEM per cycle with the least Ca2+ and Mg2+ retention in membrane and none obvious organics permeation. For fouled RO and PEM with ion transport impediment, water splitting was exacerbated, which decreased the percentage of oxidation for organics. Overall, replacing of PEM with UF significantly reduce the costs associated with both the investment and operation of MER, which is expected to broaden the practical application for treating HSOW.

Molecular Modelling of Ionic Liquids: Situations When Charge Scaling Seems Insufficient

Charge scaling as an effective solution to the experiment-computation disagreement in molecular modelling of ionic liquids (ILs) could bring the computational results close to the experimental reference for various thermodynamic properties. According to the large-scale benchmark calculations of mass density, solvation, and water-ILs transfer-free energies in our series of papers, the charge-scaling factor of 0.8 serves as a near-optimal option generally applicable to most ILs, although a system-dependent parameter adjustment could be attempted for further improved performance. However, there are situations in which such a charge-scaling treatment would fail. Namely, charge scaling cannot really affect the simulation outcome, or minimally perturbs the results that are still far from the experimental value. In such situations, the vdW radius as an additional adjustable parameter is commonly tuned to minimize the experiment-calculation deviation. In the current work, considering two ILs from the quinuclidinium family, we investigate the impacts of this vdW-scaling treatment on the mass density and the solvation/partition thermodynamics in a fashion similar to our previous charge-scaling works, i.e., scanning the vdW-scaling factor and computing physical properties under these parameter sets. It is observed that the mass density exhibits a linear response to the vdW-scaling factor with slopes close to -1.8 g/mL. By further investigating a set of physiochemically relevant temperatures between 288 K and 348 K, we confirm the robustness of the vdW-scaling treatment in the estimation of bulk properties. The best vdW-scaling parameter for mass density would worsen the computation of solvation/partition thermodynamics, and a marginal decrease in the vdW-scaling factor is considered as an intermediate option balancing the reproductions of bulk properties and solvation thermodynamics. These observations could be understood in a way similar to the charge-scaling situation. i.e., overfitting some properties (e.g., mass density) would degrade the accuracy of the other properties (e.g., solvation free energies). Following this principle, the general guideline for applying this vdW-tuning protocol is by using values between the density-derived choice and the solvation/partition-derived solution. The charge and current vdW scaling treatments cover commonly encountered ILs, completing the protocol for accurate modelling of ILs with fixed-charge force fields.

Neutron Diffraction Study of Indole Solvation in Deep Eutectic Systems of Choline Chloride, Malic Acid, and Water

Deep eutectic systems are currently under intense investigation to replace traditional organic solvents in a range of syntheses. Here, indole in choline chloride‐malic acid deep eutectic solvent (DES) was studied as a function of water content, to identify solute interactions with the DES which affect heterocycle reactivity and selectivity, and as a proxy for biomolecule solvation. Empirical Potential Structure Refinement models of neutron diffraction data showed [Cholinium]⁺ cations associate strongly with the indole π‐system due to electrostatics, whereas malic acid is only weakly associated. Trace water is sequestered into the DES and does not interact strongly with indole. When water is added to the DES, it does not interact with the indole π‐system but is exclusively in‐plane with the heterocyclic rings, forming strong H‐bonds with the ‐NH group, and also weak H‐bonds and thus prominent hydrophobic hydration of the indole aromatic region, which could direct selectivity in reactions.

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.

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.

Deep eutectic solvents-The vital link between ionic liquids and ionic solutions

When selecting a solvent for a given solute, the strongly held idiom "like dissolves like", meaning that polar solvents are used for polar solutes, is often used. This idea has resulted from the concept that most molecular solvents are homogeneous. In a deep eutectic solvent (DES), however, both components can be ionic or non-ionic, polar or non-polar. By tuning the components, DESs can solubilize a wide variety of solutes, often mixing hydrophobic and hydrophilic components, and the mixture can be designed to control phase behavior. The liquids often contain significant short-length order, and preferential solvation of one component often occurs. The addition of small polar molecules such as water or alcohols results in non-homogeneous liquids, which have significantly decreased viscosity and increased ionic conductivity. Accordingly, the areas covered in this special issue focus on structure and dynamics, solvation, the mobility of charged species, and the ability to obtain controllable phase behavior by adding polar diluents or using hydrophobic 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.