CO2solvation free energy using quasi-chemical theory

Transport and Energetics of Carbon Dioxide in Ionic Liquids at Aqueous Interfaces

A major hurdle in utilizing carbon dioxide (CO2) lies in separating it from industrial flue gas mixtures and finding suitable storage methods that enable its application in various industries. To address this issue, we utilized a combination of molecular dynamics simulations and experiments to investigate the behavior of CO2 in common room-temperature ionic liquids (RTIL) when in contact with aqueous interfaces. Our investigation of RTILs, [EMIM][TFSI] and [OMIM][TFSI], and their interaction with a pure water layer mimics the environment of a previously developed ultrathin enzymatic liquid membrane for CO2 separation. We analyzed diffusion constants and viscosity, which reveals that CO2 molecules exhibit faster mobility within the selected ILs compared to what would be predicted solely based on the viscosity of the liquids using the standard Einstein-Stokes relation. Moreover, we calculated the free energy of translocation for various species across the aqueous-IL interface, including CO2 and HCO3-. Free energy profiles demonstrate that CO2 exhibits a more favorable partitioning behavior in the RTILs compared to that in pure water, while a significant barrier hinders the movement of HCO3- from the aqueous layer. Experimental measurement of the CO2 transport in the RTILs corroborates the model. These findings strongly suggest that hydrophobic RTILs could serve as a promising option for selectively transporting CO2 from aqueous media and concentrating it as a preliminary step toward storage.

Carbonic Anhydrase Robustness for Use in Nanoscale CO2 Capture Technologies

Continued dependence on crude oil and natural gas resources for fossil fuels has caused global atmospheric carbon dioxide (CO2) emissions to increase to record-setting proportions. There is an urgent need for efficient and inexpensive carbon sequestration systems to mitigate large-scale emissions of CO2 from industrial flue gas. Carbonic anhydrase (CA) has shown high potential for enhanced CO2 capture applications compared to conventional absorption-based methods currently utilized in various industrial settings. This study aims to understand structural aspects that contribute to the stability of CA enzymes critical for their applications in industrial processes, which require the ability to withstand conditions different from those in their native environments. Here, we evaluated the thermostability and enzyme activity of mesophilic and thermophilic CA variants at different temperature conditions and in the presence of atmospheric gas pollutants like nitrogen oxides and sulfur oxides. Based on our enzyme activity assays and molecular dynamics simulations, we see increased conformational stability and CA activity levels in thermostable CA variants incubated week-long at different temperature conditions. The thermostable CA variants also retained high levels of CA activity despite changes in solution pH due to increasing NO and SO2 concentrations. A loss of CA activity was observed only at high concentrations of NO/SO2 that possibly can be minimized with the appropriate buffered solutions.

Control of the Structural Charge Distribution and Hydration State upon Intercalation of CO2 into Expansive Clay Interlayers

Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO₂ preferentially into their interlayers and therefore potentially act as a material for CO₂ separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO₂ into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO₂ intercalation. Our free energy calculations indicate that CO₂ favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO₂/H₂O ratio in clay interlayers as compared to that in bulk water. CO₂ molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO₂ intercalation. The structural charge distribution significantly affects the amount of CO₂ intercalated in the interlayers. These results provide a mechanistic understanding of CO₂ intercalation in clays for CO₂ separation, capture, and storage.

Probing dissolved CO2(aq) in aqueous solutions for CO2 electroreduction and storage

CO₂ dissolved in aqueous solutions CO₂(aq) is important to CO₂ capture, storage, photo-/electroreduction in the fight against global warming and to CO₂ analysis in drinks. Here, we developed microscale infrared (IR) spectroscopy for in situ dynamic quantitating CO₂(aq). The quantized CO₂(g) rotational state transitions were observed to quench for CO₂(aq), accompanied by increased H₂O IR absorption. An accurate CO₂ molar extinction coefficient ε was derived for in situ CO₂(aq) quantification up to 58 atm. We directly measured CO₂(aq) concentrations in electrolytes under CO₂(g) bubbling and high-pressure conditions with high spectral and time resolutions. In KHCO₃ electrolytes with CO₂(aq) > ~1 M, CO₂ electroreduction (CO₂RR) to formate reached >98% Faradaic efficiencies on copper (Cu₂O/Cu)-based electrocatalyst. Furthermore, CO₂ dissolution/desolvation kinetics showed large hysteresis and ultraslow reversal of CO₂(aq) supersaturation in aqueous systems, with implications to CO₂ capture, storage, and supersaturation phenomena in natural water bodies.

The role of hydrophobic hydration in the free energy of chemical reactions at the gold/water interface: Size and position effects

Metal/water interfaces catalyze a large variety of chemical reactions, which often involve small hydrophobic molecules. In the present theoretical study, we show that hydrophobic hydration at the Au(100)/water interface actively contributes to the reaction free energy by up to several hundreds of meV. This occurs either in adsorption/desorption reaction steps, where the vertical distance from the surface changes in going from reactants to products, or in addition and elimination reaction steps, where two small reactants merge into a larger product and vice versa. We find that size and position effects cannot be captured by treating them as independent variables. Instead, their simultaneous evaluation allows us to map the important contributions, and we provide examples of their combinations for which interfacial reactions can be either favored or disfavored. By taking a N₂ and a CO₂ reduction pathway as test cases, we show that explicitly considering hydrophobic effects is important for the selectivity and rate of these relevant interfacial processes.

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.

Understanding CO2 capture kinetics and energetics by ionic liquids with molecular dynamics simulation

The interactions between ionic liquids (ILs) and CO₂ were studied by molecular dynamics simulations. Several key characters, including the volumes of cations and anions, the length of the alkyl chain have been assessed on CO₂ capture capability.

Machine Learning-Guided Approach for Studying Solvation Environments

Molecular-level understanding and characterization of solvation environments are often needed across chemistry, biology, and engineering. Toward practical modeling of local solvation effects of any solute in any solvent, we report a static and all-quantum mechanics-based cluster-continuum approach for calculating single-ion solvation free energies. This approach uses a global optimization procedure to identify low-energy molecular clusters with different numbers of explicit solvent molecules and then employs the smooth overlap for atomic positions learning kernel to quantify the similarity between different low-energy solute environments. From these data, we use sketch maps, a nonlinear dimensionality reduction algorithm, to obtain a two-dimensional visual representation of the similarity between solute environments in differently sized microsolvated clusters. After testing this approach on different ions having charges 2+, 1+, 1-, and 2-, we find that the solvation environment around each ion can be seen to usually become more similar in hand with its calculated single-ion solvation free energy. Without needing either dynamics simulations or an a priori knowledge of local solvation structure of the ions, this approach can be used to calculate solvation free energies within 5% of experimental measurements for most cases, and it should be transferable for the study of other systems where dynamics simulations are not easily carried out.

Quasi-Chemical Theory with Cluster Sampling from Ab Initio Molecular Dynamics: Fluoride (F-) Anion Hydration

Accurate predictions of the hydration free energy for anions typically has been more challenging than that for cations. Hydrogen bond donation to the anion in hydrated clusters such as F(H₂O) _n ^- can lead to delicate structures. Consequently, the energy landscape contains many local minima, even for small clusters, and these minima present a challenge for computational optimization. Utilization of cluster experimental results for the free energies of gas-phase clusters shows that even though anharmonic effects are interesting they need not be of troublesome magnitudes for careful applications of quasi-chemical theory to ion hydration. Energy-optimized cluster structures for anions can leave the central ion highly exposed, and application of implicit solvation models to these structures can incur more serious errors than those for metal cations. Utilizing cluster structures sampled from ab initio molecular dynamics simulations substantially fixes those issues.

Structure, Dynamics, and Hydration Free Energy of Carbon Dioxide in Aqueous Solution: A Quantum Mechanical/Molecular Mechanics Molecular Dynamics Thermodynamic Integration (QM/MM MD TI) Simulation Study

The solvation of carbon dioxide in solution represents a key step for the capture and fixation CO₂ in nature, which may be further influenced by the formation of (bi)carbonate species and/or the formation of CO₂ clusters in solution. The latter processes are strongly dependent on the exact environment of the liquid state (e.g., pH value, solvated ions, etc.) and may interfere with the experimental determination of structural, dynamical, and thermodynamic properties. In this work a hybrid quantum mechanical/molecular mechanical (QM/MM) simulation approach at correlated ab initio level of theory resolution-of-identity second-order Møller-Plesset Perturbation Theory (RI-MP2) has been applied in the framework of thermodynamic integration (TI) to study structure, dynamics, and the hydration free energy of a single carbon dioxide molecule in aqueous solution. A detailed analysis of the individual QM/MM potential energy contributions demonstrate that the overall potential remains highly consistent over the entire sampling phase and that no artificial contributions are influencing the determination of the hydration free energy. The latter value of 0.01 ± 0.92 kcal/mol was found in very good agreement with the values of 0.06 and 0.24 kcal/mol obtained via quasi-chemical theory and experimental measurements, respectively. In order to obtain detailed information about the C- and O _C-water interaction, conically restricted regions with respect to the main axis of the CO₂ molecule have been employed in structural analysis. The presented data not only provide detailed information about the hydration properties of CO₂ but act as a critical validation of the simulation technique, which will be beneficial in the study of nonaqueous solvents such as pure and aqueous NH₃ solutions, which have been suggested as potential candidates to capture CO₂ from anthropogenic sources.

Density Fluctuation in Aqueous Solutions and Molecular Origin of Salting-Out Effect for CO2

Using molecular dynamics simulation, we studied the density fluctuations and cavity formation probabilities in aqueous solutions and their effect on the hydration of CO₂. With increasing salt concentration, we report an increased probability of observing a larger than the average number of species in the probe volume. Our energetic analyses indicate that the van der Waals and electrostatic interactions between CO₂ and aqueous solutions become more favorable with increasing salt concentration, favoring the solubility of CO₂ (salting in). However, due to the decreasing number of cavities forming when salt concentration is increased, the solubility of CO₂ decreases. The formation of cavities was found to be the primary control on the dissolution of gas, and is responsible for the observed CO₂ salting-out effect. Our results provide the fundamental understanding of the density fluctuation in aqueous solutions and the molecular origin of the salting-out effect for real gas.

Statistical Analyses of Hydrophobic Interactions: A Mini-Review

This review focuses on the striking recent progress in solving for hydrophobic interactions between small inert molecules. We discuss several new understandings. First, the inverse temperature phenomenology of hydrophobic interactions, i.e., strengthening of hydrophobic bonds with increasing temperature, is decisively exhibited by hydrophobic interactions between atomic-scale hard sphere solutes in water. Second, inclusion of attractive interactions associated with atomic-size hydrophobic reference cases leads to substantial, nontrivial corrections to reference results for purely repulsive solutes. Hydrophobic bonds are weakened by adding solute dispersion forces to treatment of reference cases. The classic statistical mechanical theory for those corrections is not accurate in this application, but molecular quasi-chemical theory shows promise. Finally, because of the masking roles of excluded volume and attractive interactions, comparisons that do not discriminate the different possibilities face an interpretive danger.

Dispersion- and Exchange-Corrected Density Functional Theory for Sodium Ion Hydration

A challenge in density functional theory is developing functionals that simultaneously describe intermolecular electron correlation and electron delocalization. Recent exchange-correlation functionals address those two issues by adding corrections important at long ranges: an atom-centered pairwise dispersion term to account for correlation and a modified long-range component of the electron exchange term to correct for delocalization. Here we investigate how those corrections influence the accuracy of binding free energy predictions for sodium-water clusters. We find that the dual-corrected ωB97X-D functional gives cluster binding energies closest to high-level ab initio methods (CCSD(T)). Binding energy decomposition shows that the ωB97X-D functional predicts the smallest ion-water (pairwise) interaction energy and larger multibody contributions for a four-water cluster than most other functionals - a trend consistent with CCSD(T) results. Also, ωB97X-D produces the smallest amounts of charge transfer and the least polarizable waters of the density functionals studied, which mimics the lower polarizability of CCSD. When compared with experimental binding free energies, however, the exchange-corrected CAM-B3LYP functional performs best (error <1 kcal/mol), possibly because of its parametrization to experimental formation enthalpies. For clusters containing more than four waters, "split-shell" coordination must be considered to obtain accurate free energies in comparison with experiment.

Octa-Coordination and the Aqueous Ba(2+) Ion

The hydration structure of Ba(2+) ion is important for understanding blocking mechanisms in potassium ion channels. Here, we combine statistical mechanical theory, ab initio molecular dynamics simulations, and electronic structure methods to calculate the hydration free energy and local hydration structure of Ba(2+)(aq). The predicted hydration free energy (-304 ± 1 kcal/mol) agrees with the experimental value (-303 kcal/mol) when a maximally occupied, unimodal inner solvation shell is treated. In the local environment defined by the first shell of hydrating waters, Ba(2+) is directly and stably coordinated by eight (8) waters. Octa-coordination resembles the crystal structure of Ba(2+) and K(+) bound in potassium ion channels, but differs from the local hydration structure of K(+)(aq) determined earlier.

Aqueous microsolvation of CdCl₂: density functional theory and Born-Oppenheimer molecular dynamics studies

We report a systematic study of aqueous microsolvation of CdCl2. The optimized structures and binding energies of the CdCl2-(H2O)n clusters with n = 1-24 have been computed at the B3PW91/6-31G** level. The solvation patterns obtained at the DFT level are verified at the MP2/AVTZ level for n < 6. Unlike HgCl2-(H2O)n case, where there are at most three Hg-O(w) orbital interactions, Cd also establishes four equatorial orbital interactions with water for n > 6 leading to a planar square bipyramid hexacoordination around Cd. The first solvation shell is fully attained with 12 water molecules. At the same level of theory the water binding energies are much larger than those previously found for HgCl2 due to the stronger Cd-O(w) interactions arising from the smaller core of Cd. For the largest system studied, CdCl2-(H2O)24, both penta- and hexa-coordination stable patterns around Cd are found. However, Born-Opphenheimer molecular dynamics simulations starting from these optimized geometries at 700 K reveal the greater stability of the Cd-pentacoordinated species, where a CdCl2-(H2O)3 trigonal bipyramid effective solute appears. The Cd-O(water) radial distribution function shows a bimodal distribution with two maxima at 2.4 Å and 4.2 Å, revealing the different coordination spheres, even with such a small number of solvating water molecules.