Computation of gas solubilities in choline chloride urea and choline chloride ethylene glycol deep eutectic solvents using Monte Carlo simulations

Unveiling CO2 capture in tailorable green neoteric solvents: An ensemble learning approach informed by quantum chemistry

In the relentless battle against the impending climate crisis, deep eutectic solvents (DESs) have emerged as beacons of hope in the realm of green chemistry, igniting a resurgence of scientific exploration. These versatile compounds hold the promise of revolutionizing carbon capture, effectively countering the rising tide of carbon dioxide (CO2) emissions responsible for global warming and climate instability. Their adaptability offers a tantalizing prospect, as they can be finely tailored for a multitude of applications, thereby encompassing the uncharted territory of potential DESs. Navigating this unexplored terrain underscores the vital need for predictive computational methods, which serve as our guiding compass in the expansive landscape of DESs. Thermodynamic modeling and solubility prognostications stand as our unwavering navigational aides on this treacherous odyssey. In this direction, the COSMO-RS model intertwined with the captivating Stochastic Gradient Boosting (SGB) algorithm. Together, they unveil the elusive truths pertaining to CO2 solubility in DESs, forging a path toward a sustainable future. Our quest is substantiated by two exhaustive datasets, a repository of knowledge encompassing 1973 and 2327 CO2 solubility data points spanning 132 and 150 distinct DESs respectively, encapsulating a spectrum of conditions. The SGB models, incorporating features derived from COSMO-RS, as well as accounting for pressure and temperature variables, furnishes predictions that harmonize seamlessly with experimental CO2 solubility values, boasting an impressive Average Absolute Relative Deviation (AARD) of a mere 0.85% and 2.30% respectively. When juxtaposed with literature-reported methodologies like different EoS, as well as Computational Solvation, and machine learning (ML) models, our SGB model emerges as the epitome of reliability, offering robust forecasts of CO2 solubility in DESs. It emerges as a potent tool for the design and selection of DESs for CO2 capture and utilization, heralding a sustainable and environmentally conscientious future in the battle against climate change.

Interfacial Tensions, Solubilities, and Transport Properties of the H2/H2O/NaCl System: A Molecular Simulation Study

Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H2), such as underground H2 storage and H2O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H2O, (b) self-diffusivities of infinitely diluted H2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O, and (c) solubilities of H2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O. The force fields used are the TIP4P/2005 for H2O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H2. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.

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

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

Novel hybrid QSPR-GPR approach for modeling of carbon dioxide capture using deep eutectic solvents

In recent years, deep eutectic solvents (DESs) have garnered considerable attention for their potential in carbon capture and utilization processes. Predicting the carbon dioxide (CO2) solubility in DES is crucial for optimizing these solvent systems and advancing their application in sustainable technologies. In this study, we presented an evolving hybrid Quantitative Structure-Property Relationship and Gaussian Process Regression (QSPR-GPR) model that enables accurate predictions of CO2 solubility in various DESs. The QSPR-GPR model combined the strengths of both approaches, leveraging molecular descriptors and structural features of DES components to establish a robust and adaptable predictive framework. Through a systematic evolution process, we iteratively refined the model, enhancing its performance and generalization capacity. By incorporating experimental CO2 solubility data in varied DES compositions and temperatures, we trained the model to capture the intricate solubility behaviour precisely. The analytical capability of the evolving hybrid model was validated against an extensive dataset of experimental CO2 solubility values, demonstrating its superiority over individual QSPR and GPR models. The model achieves high accuracy, capturing the complex interactions between CO2 and DES components under varying thermodynamic conditions. The versatility of the evolving hybrid model was highlighted by its ability to accommodate new experimental data and adapt to different DES compositions and temperatures. The proposed QSPR-GPR model presented a powerful tool for predicting CO2 solubility in DES, providing valuable insights for designing and optimizing solvent systems in carbon capture technologies. The model's remarkable performance enhances our understanding of CO2 solubility mechanisms and contributes to sustainable solutions for mitigating greenhouse gas emissions. As research in DESs progresses, the evolving hybrid QSPR-GPR model offers a versatile and accurate means for predicting CO2 solubility, supporting advancements in carbon capture and utilization processes towards a greener and more sustainable future.

Everything You Wanted to Know about Deep Eutectic Solvents but Were Afraid to Be Told

Are deep eutectic solvents (DESs) a promising alternative to conventional solvents? Perhaps, but their development is hindered by a plethora of misconceptions. These are carefully analyzed here, beginning with the very meaning of DESs, which has strayed far beyond its original scope of eutectic mixtures of Lewis or Brønsted acids and bases. Instead, a definition that is grounded on thermodynamic principles and distinguishes between eutectic and deep eutectic is encouraged, and the types of precursors that can be used to prepare DESs are reviewed. Landmark works surrounding the sustainability, stability, toxicity, and biodegradability of these solvents are also discussed, revealing piling evidence that numerous DESs reported thus far, particularly those that are choline based, lack sufficient sustainability-related traits to be considered green solvents. Finally, emerging DES applications are reviewed, emphasizing their most remarkable feature: the ability to liquefy a solid compound with a target property, allowing its use as a liquid solvent.

Solubility of CO2 in Aqueous Formic Acid Solutions and the Effect of NaCl Addition: A Molecular Simulation Study

There is a growing interest in the development of routes to produce formic acid from CO₂, such as the electrochemical reduction of CO₂ to formic acid. The solubility of CO₂ in the electrolyte influences the production rate of formic acid. Here, the dependence of the CO₂ solubility in aqueous HCOOH solutions with electrolytes on the composition and the NaCl concentration was studied by Continuous Fractional Component Monte Carlo simulations at 298.15 K and 1 bar. The chemical potentials of CO₂, H₂O, and HCOOH were obtained directly from single simulations, enabling the calculation of Henry coefficients and subsequently considering salting in or salting out effects. As the force fields for HCOOH and H₂O may not be compatible due to the presence of strong hydrogen bonds, the Gibbs-Duhem integration test was used to test this compatibility. The combination of the OPLS/AA force field with a new set of parameters, in combination with the SPC/E force field for water, was selected. It was found that the solubility of CO₂ decreases with increasing NaCl concentration in the solution and increases with the increase of HCOOH concentration. This continues up to a certain concentration of HCOOH in the solution, after which the CO₂ solubility is high and the NaCl concentration has no significant effect.

A New Force Field for OH- for Computing Thermodynamic and Transport Properties of H2 and O2 in Aqueous NaOH and KOH Solutions

The thermophysical properties of aqueous electrolyte solutions are of interest for applications such as water electrolyzers and fuel cells. Molecular dynamics (MD) and continuous fractional component Monte Carlo (CFCMC) simulations are used to calculate densities, transport properties (i.e., self-diffusivities and dynamic viscosities), and solubilities of H2 and O2 in aqueous sodium and potassium hydroxide (NaOH and KOH) solutions for a wide electrolyte concentration range (0-8 mol/kg). Simulations are carried out for a temperature and pressure range of 298-353 K and 1-100 bar, respectively. The TIP4P/2005 water model is used in combination with a newly parametrized OH- force field for NaOH and KOH. The computed dynamic viscosities at 298 K for NaOH and KOH solutions are within 5% from the reported experimental data up to an electrolyte concentration of 6 mol/kg. For most of the thermodynamic conditions (especially at high concentrations, pressures, and temperatures) experimental data are largely lacking. We present an extensive collection of new data and engineering equations for H2 and O2 self-diffusivities and solubilities in NaOH and KOH solutions, which can be used for process design and optimization of efficient alkaline electrolyzers and fuel cells.

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

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

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.

Interfacial Properties of Hydrophobic Deep Eutectic Solvents with Water

Hydrophobic deep eutectic solvents (DESs) have recently gained much attention as water-immiscible solvents for a wide range of applications. However, very few studies exist in which the hydrophobicity of these DESs is quantified. In this work, the interfacial properties of hydrophobic DESs with water were computed at various temperatures using molecular dynamics simulations. The considered DESs were tetrabutylammonium chloride–decanoic acid (TBAC–dec) with a molar ratio of 1:2, thymol–decanoic acid (Thy–dec) with a molar ratio of 1:2, and dl-menthol–decanoic acid (Men–dec) with a molar ratio of 2:1. The following properties were investigated in detail: interfacial tensions, water-in-DES solubilities (and salt-in-water solubilities for TBAC–dec/water), density profiles, and the number densities of hydrogen bonds. Different ionic charge scaling factors were used for TBAC–dec. Thy–dec and Men–dec showed a high level of hydrophobicity with negligible computed water-in-DES solubilities. For charge scaling factors of 0.7 and 1 for the thymol and decanoic acid components of Thy–dec, the computed interfacial tensions of the DESs are in the following order: TBAC–dec (ca. 4 mN m^–1) < Thy–dec (20 mN m^–1) < Men–dec (26 mN m^–1). The two sets of charge scaling factors for Thy–dec did not lead to different density profiles but resulted in considerable differences in the DES/water interfacial tensions due to different numbers of decanoic acid–water hydrogen bonds at the interfaces. Large peaks were observed for the density profiles of (the hydroxyl oxygen of) decanoic acid at the interfaces of all DES/water mixtures, indicating a preferential alignment of the oxygen atoms of decanoic acid toward the aqueous phase.

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.

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.

New Features of the Open Source Monte Carlo Software Brick-CFCMC: Thermodynamic Integration and Hybrid Trial Moves

We present several new major features added to the Monte Carlo (MC) simulation code Brick-CFCMC for phase- and reaction equilibria calculations (https://gitlab.com/ETh_TU_Delft/Brick-CFCMC). The first one is thermodynamic integration for the computation of excess chemical potentials (μ^ex). For this purpose, we implemented the computation of the ensemble average of the derivative of the potential energy with respect to the scaling factor for intermolecular interactions (). Efficient bookkeeping is implemented so that the quantity is updated after every MC trial move with negligible computational cost. We demonstrate the accuracy and reliability of the calculation of μ^ex for sodium chloride in water. Second, we implemented hybrid MC/MD translation and rotation trial moves to increase the efficiency of sampling of the configuration space. In these trial moves, short Molecular Dynamics (MD) trajectories are performed to collectively displace or rotate all molecules in the system. These trajectories are accepted or rejected based on the total energy drift. The efficiency of these trial moves can be tuned by changing the time step and the trajectory length. The new trial moves are demonstrated using MC simulations of a viscous fluid (deep eutectic solvent).

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

Deep eutectic solvents (DESs) have emerged as a cheaper and greener alternative to conventional organic solvents. Choline chloride (ChCl) mixed with urea at a molar ratio of 1:2 is one of the most common DESs for a wide range of applications such as electrochemistry, material science, and biochemistry. In this study, molecular dynamics simulations are performed to study the effect of urea content on the thermodynamic and transport properties of ChCl and urea mixtures. With increased mole fraction of urea, the number of hydrogen bonds (HBs) between cation-anion and ion-urea decreases, while the number of HBs between urea-urea increases. Radial distribution functions (RDFs) for ChCl-urea and ChCl-ChCl pairs shows a significant decrease as the mole fraction of urea increases. Using the computed RDFs, Kirkwood-Buff Integrals (KBIs) are computed. KBIs show that interactions of urea-urea become stronger, while interactions of urea-ChCl and ChCl-ChCl pairs become slightly weaker with increasing mole fraction of urea. All thermodynamic factors are found larger than one, indicating a non-ideal mixture. Our results also show that self- and collective diffusivities increase, while viscosities decrease with increasing urea content. This is mainly due to the weaker interactions between ions and urea, resulting in enhanced mobilities. Ionic conductivities exhibit a non-monotonic behavior. Up to a mole fraction of 0.5, the ionic conductivities increase with increasing urea content and then reach a plateau.

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.

Solubility of Carbon Dioxide, Hydrogen Sulfide, Methane, and Nitrogen in Monoethylene Glycol; Experiments and Molecular Simulation

Knowledge on the solubility of gases, especially carbon dioxide (CO₂), in monoethylene glycol (MEG) is relevant for a number of industrial applications such as separation processes and gas hydrate prevention. In this study, the solubility of CO₂ in MEG was measured experimentally at temperatures of 333.15, 353.15, and 373.15 K. Experimental data were used to validate Monte Carlo (MC) simulations. Continuous fractional component MC simulations in the osmotic ensemble were performed to compute the solubility of CO₂ in MEG at the same temperatures and at pressures up to 10 bar. MC simulations were also used to study the solubility of methane (CH₄), hydrogen sulfide (H₂S), and nitrogen (N₂) in MEG at 373.15 K. Solubilities from experiments and simulations are in good agreement at low pressures, but deviations were observed at high pressures. Henry coefficients were also computed using MC simulations and compared to experimental values. The order of solubilities of the gases in MEG at 373.15 K was computed as H₂S > CO₂ > CH₄ > N₂. Force field modifications may be required to improve the prediction of solubilities of gases in MEG at high pressures and low temperatures.