On the Performance of Confined Deep Eutectic Solvents and Ionic Liquids for Separations of Carbon Dioxide from Methane: Molecular Dynamics Simulations

Interfacial Nanostructure and Hydrogen Bond Networks of Choline Chloride and Glycerol Mixtures Probed with X-ray and Vibrational Spectroscopies

The molecular distribution at the liquid-vapor interface and evolution of the hydrogen bond interactions in mixtures of glycerol and choline chloride are investigated using X-ray photoelectron spectroscopy (XPS) and Raman spectroscopy. Nanoscale depth profiles of supersaturated deep eutectic solvent (DES) mixtures up to ∼2 nm measured by ambient-pressure XPS show the enhancement of choline cation (Ch+) concentration by a factor of 2 at the liquid-vapor interface compared to the bulk. In addition, Raman spectral analysis of a wide range of DES mixtures reveals the conversion of gauche-conformer Ch+ into the anti-conformer in relatively lower ChCl concentrations. Finally, the depletion of Ch+ from the interface (probing depth = 0.4 nm) is demonstrated by aerosol-based velocity map imaging XPS measurements of glyceline and water mixtures. The nanostructure of liquid-vapor interfaces and structural rearrangement by hydration can provide critical insight into the molecular origin of the deep eutectic behavior and gas-capturing application of DESs.

Quantum Chemical Simulations of CO2 and N2 Capture in Reline, a Prototypical Deep Eutectic Solvent

Deep eutectic solvents such as reline are an emerging class of low-cost, environmentally friendly solvents with tunable properties that are potentially applicable for the capture and separation of CO2. Experimental measurements showed that a reline-based membrane contactor can capture and separate CO2 via physisorption through a dissolution process with 96.7% purity from a mixed gas containing CO2 and N2 (50:50% molar ratio). We examine the nature of the interaction of CO2 and N2 with reline employing quantum chemical methods. We focus on explaining the mechanism by which CO2 and N2 bind to reline and the reason for the high selectivity for absorption of CO2 compared to N2. We analyze the dynamics, energetics, and binding motifs for CO2 and N2 in reline employing density functional theory, density functional tight binding, and ab initio molecular dynamics. We also investigate the effect of reline on the vibrational spectra of CO2 and reline. Our simulations indicate that the selective capture of CO2 from the mixture of CO2 and N2 is due to the interplay between attractive electrostatic and charge polarization forces with opposing entropic effects, which shift the energetic balance and make the N2 absorption unfavorable in reline.

The correlation between the micelle morphology of surface-active ionic liquids with self-assembly and thermodynamic characteristics: coarse-grained MD simulation and experiment

Surface-active ionic liquids (SAILs) show great promise as novel green solvents due to their low vapor pressure, high thermal stability, high electrical conductivity, and bio-friendly nature to replace traditional volatile organic solvents in industrial processes. In the present work, the combination of coarse-grained (CG) molecular dynamics (MD) simulations with conductivity measurements was employed to explain the correlation between the micelle morphology and physicochemical and thermodynamic properties of self-assembly. A homologous series of SAIL molecules, 1-n-alkyl-3-methylimidazolium bromide [Cnmim][Br] (n = 4, 6, 8, 10, and 12), were chosen at various concentrations to shed light on this issue. Simultaneously two factors of concentration and alkyl chain length affected the morphology to control the physical and thermodynamic features. Moreover, the nature of the headgroup for two SAILs with the longest alkyl chain was assessed by shifting from imidazolium into ammonium. First, the critical micelle concentration (CMC), the degree of counterion dissociation of micelles, and the standard Gibbs energy of micellization of SAILs were determined using conductivity data. The micelle morphology such as the aggregation number, micelle radius, and moment of inertia was computed before, around, and after the CMC by MD simulation. Simulated results in accordance with the experimental measurements provide a quantitative understanding of the micellar properties. Increasing the alkyl chain length was associated with a non-spherical bigger micelle while the ammonium-based surfactant with a lower repulsion between neighboring monomers in micelles induced bigger and more spherical aggregates. Raising the SAIL concentration did not considerably influence the sphericity of the micelle except for the SAIL with the longest tail. The umbrella sampling method calculated the potential of mean force (PMF) for pulling a monomer of SAIL from a pre-assembled micelle into the solution. The dissociation energy of a SAIL monomer from a micelle increased with the tail length or with shifting into the ammonium head group and was substantially influenced by micelle morphology. Comparison between a sphere micelle with an oval one demonstrated that the dissociation of a SAIL monomer from a non-spherical shape needed a higher amount of energy. An improved understanding of how the shape of the SAIL micelles controls the physicochemical properties and stability helps to extend their application to different chemical processes.

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

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

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

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

Boundary-Monte Carlo Method for Neutral and Charged Confined Fluids

In this work, we describe a new Monte Carlo (MC) simulation method to investigate highly coupled fluids in confined geometries at a constant chemical potential. This method is based on so-called multi-scale Hamiltonian methods, wherein the chemical potential is determined using a more amenable Hamiltonian for a fluid in an "outer" region, which facilitates standard methods, such as grand canonical MC simulations or Widom's particle insertion method. The (inner region) fluid of interest is placed in diffusive contact with the simpler outer fluid via a boundary zone wherein the Hamiltonian is transformed. The current method utilizes an ideal fluid for the outer regions, which allows for implicit rather than explicit simulations. Only the boundary and inner region need explicit consideration; hence, the nomenclature used is boundary-Monte Carlo. We illustrate the utility of the method for simple neutral and charged fluids in cylindrical and planar pores. In the latter case, we use a dense room-temperature ionic liquid model and illustrate how the boundary zone establishes a proper Donnan equilibrium between inner and outer fluids in the presence of charged planar electrodes. Thus, the method allows direct calculation of properties such as the differential capacitance, without the need for additional difficult calculations of the requisite Donnan potential.

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.

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

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

New Functionalized Ionic Liquids Based on POSS for the Detection of Fe3+ Ion

This study reports a novel series of imidazolium ionic liquids (ILs) based on polyhedral oligomeric silsesquioxanes (POSS) towards effective detection of metal ions, especially Fe3+ ion. 1H, 13C, 29Si NMR, high resolution mass spectra (HRMS) and Fourier transform infrared spectra (FTIR) were applied to confirm the structures of the ILs based on POSS (ILs-POSS). The three ILs-POSS were synthesized via a green chemistry approach-a thiol-ene "click" reaction without metal ions as catalysts. Furthermore, the spherical vesicle structures of the ILs-POSS were observed and caused by self-assembly behaviors. Through comprehensive characterizations, these ILs-POSS have performed excellent thermal stabilities and low glass transition temperatures. In addition, we found it very meaningful that the limits of detection (LODs) of the three ILs-POSS for the detection of the Fe3+ ion were 7.91 × 10-8 M, 1.2 × 10-7 M, and 1.2 × 10-7 M, respectively. These data illustrate that these ILs-POSS have great potential for the detection of the Fe3+ ion. In conclusion, this study not only prepared novel ILs-POSS, but also provided new materials as fluorescent sensors in the detection of Fe3+.

Deep Eutectic Solvents: A Review of Fundamentals and Applications

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

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

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

Local Grand Canonical Monte Carlo Simulation Method for Confined Fluids

We describe a new local grand canonical Monte Carlo method to treat fluids in pores in chemical equilibrium with a reference bulk. The method is applied to Lennard-Jones particles in pores of different geometry and is shown to be much more accurate and efficient than other techniques such as traditional grand canonical simulations or Widom's particle insertion method. It utilizes a penalty potential to create a gas phase, which is in equilibrium with a more dense liquid component in the pore. Grand canonical Monte Carlo moves are employed in the gas phase, and the system then maintains chemical equilibrium by "diffusion" of particles. This creates an interface, which means that the confined fluid needs to occupy a large enough volume so that this is not an issue. We also applied the method to confined charged fluids and show how it can be used to determine local electrostatic potentials in the confined fluid, which are properly referenced to the bulk. This precludes the need to determine the Donnan potential (which controls electrochemical equilibrium) explicitly. Prior approaches have used explicit bulk simulations to measure this potential difference, which are significantly costly from a computational point of view. One outcome of our analysis is that pores of finite cross-section create a potential difference with the bulk via a small but nonzero linear charge density, which diminishes as ∼1/ln(L), where L is the pore length.