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.

Quantitative Simulations of Siloxane Adsorption in Metal-Organic Frameworks

We present a transferable force field (FF) for simulating the bulk properties of linear and cyclic siloxanes and the adsorption of these species in metal-organic frameworks (MOFs). Unlike previous FFs for siloxanes, our FF accurately reproduces the vapor-liquid equilibria of each species in the bulk phase. The quality of our FF combined with the Universal Force Field using standard Lorentz-Berthelot combining rules for MOF atoms was assessed in a wide range of MOFs without open metal sites, showing good agreement with dispersion-corrected density functional theory calculations. Predictions with this FF show good agreement with the limited experimental data for siloxane adsorption in MOFs that is available. As an example of using the FF to predict adsorption properties in MOFs, we present simulations examining entropy effects in binary linear and cyclic siloxane mixture coadsorption in the large-pore MOF with structure code FOTNIN.

Temperature Extrapolation of Henry's Law Constants and the Isosteric Heat of Adsorption

Computational screening of adsorbent materials often uses the Henry's law constant (K_H) (at a particular temperature) as a first discriminator metric due to its relative ease of calculation. The isosteric heat of adsorption in the limit of zero pressure (q_st^∞) is often calculated along with the Henry's law constant, and both properties are informative metrics of adsorbent material performance at low-pressure conditions. In this article, we introduce a method for extrapolating K_H as a function of temperature, using series-expansion coefficients that are easily computed at the same time as K_H itself; the extrapolation function also yields q_st^∞. The extrapolation is highly accurate over a wide range of temperatures when the basis temperature is sufficiently high, for a wide range of adsorbent materials and adsorbate gases. Various results suggest that the extrapolation is accurate when the extrapolation range in inverse-temperature space is limited to |β - β₀ | < 0.5 mol/kJ. Application of the extrapolation to a large set of materials is shown to be successful provided that K_H is not extremely large and/or the extrapolation coefficients converge satisfactorily. The extrapolation is also able to predict q_st^∞ for a system that shows an unusually large temperature dependence. The work provides a robust method for predicting K_H and q_st^∞ over a wide range of industrially relevant temperatures with minimal effort beyond that necessary to compute those properties at a single temperature, which facilitates the addition of practical operating (or processing) conditions to computational screening exercises.

Performance-Based Screening of Porous Materials for Carbon Capture

Computational screening methods have changed the way new materials and processes are discovered and designed. For adsorption-based gas separations and carbon capture, recent efforts have been directed toward the development of multiscale and performance-based screening workflows where we can go from the atomistic structure of an adsorbent to its equilibrium and transport properties at different scales, and eventually to its separation performance at the process level. The objective of this work is to review the current status of this new approach, discuss its potential and impact on the field of materials screening, and highlight the challenges that limit its application. We compile and introduce all the elements required for the development, implementation, and operation of multiscale workflows, hence providing a useful practical guide and a comprehensive source of reference to the scientific communities who work in this area. Our review includes information about available materials databases, state-of-the-art molecular simulation and process modeling tools, and a complete catalogue of data and parameters that are required at each stage of the multiscale screening. We thoroughly discuss the challenges associated with data availability, consistency of the models, and reproducibility of the data and, finally, propose new directions for the future of the field.

Assessment of Performance of Density Functionals for Predicting Potential Energy Curves in Hydrogen Storage Applications

The availability of accurate computational tools for modeling and simulation is vital to accelerate the discovery of materials capable of storing hydrogen (H₂) under given parameters of pressure swing and temperature. Previously, we compiled the H2Bind275 data set consisting of equilibrium geometries and assessed the performance of 55 density functionals over this data set (Veccham, S. P.; Head-Gordon, M. J. Chem. Theory Comput. 2020, 16, 4963-4982). As it is crucial for computational tools to accurately model the entire potential energy curve (PEC), in addition to the equilibrium geometry, we extended this data set with 389 new data points to include two compressed and three elongated geometries along 78 PECs for H₂ binding, forming the H2Bind78 × 7 data set. By assessing the performance of 55 density functionals on this significantly larger and more comprehensive H2Bind78 × 7 data set, we identified the best performing density functionals for H₂ binding applications: PBE0-DH, ωB97X-V, ωB97M-V, and DSD-PBEPBE-D3(BJ). The addition of Hartree-Fock exchange improves the performance of density functionals, albeit not uniformly throughout the PEC. We recommend the usage of ωB97X-V and ωB97M-V density functionals as they offer good performance for both geometries and energies. In addition, we also identified B97M-V and B97M-rV as the best semilocal density functionals for predicting H₂ binding energy at its equilibrium geometry.

Next-Generation Accurate, Transferable, and Polarizable Potentials for Material Simulations

PHAHST (potentials with high accuracy, high speed, and transferability) intermolecular potential energy functions have been developed from first principles for H₂, N₂, the noble gases, and a metal-organic material, HKUST-1. The potentials are designed from the outset to be transferable to heterogeneous environments including porous materials, interfaces, and material simulations. This is accomplished by theoretically justified choices for all functional forms, parameters, and mixing rules, including explicit polarization in every environment and fitting to high quality electronic structure calculations using methods that are tractable for real systems. The models have been validated in neat systems by comparison to second virial coefficients and bulk pressure-density isotherms. For inhomogeneous applications, our main target, comparisons are presented to previously published experimental studies on the metal-organic material HKUST-1 including adsorption, isosteric heats of adsorption, binding site locations, and binding site energies. A systematic prescription is provided for developing compatible potentials for additional small molecules and materials. The resulting models are recommended for use in complex heterogeneous simulations where existing potentials may be inadequate.

CO2 capture on natural zeolite clinoptilolite: Effect of temperature and role of the adsorption sites

In this study, the adsorption capacity of the low-cost zeolite clinoptilolite was investigated for capturing carbon dioxide (CO₂) emitted from industrial processes at moderate temperature. The CO₂ adsorption capacity of clinoptilolite (a commercial natural zeolite) and ion-exchanged (with Na⁺ and Ca²⁺) clinoptilolite were tested under both dynamic (using a fixed-bed reactor operating with 10% vol. CO₂ in N₂) and equilibrium conditions (measuring single component adsorption isotherms). The dynamic CO₂ adsorption capacity of bare clinoptilolite and ion-exchanged clinoptilolite were evaluated in the temperature range from 293 K to 338 K and the obtained breakthrough curves were compared with those of the commercial zeolite 13X (Z13X). Although the adsorption capacity of Z13X exceeded those of bare clinoptilolite and ion-exchanged clinoptilolite at 293 K, the clinoptilolite exhibited the highest CO₂ uptake at a moderate temperature of 338 K (i.e. 25 % higher than Z13X). This feature appears in agreement with the lower isosteric heat of CO₂ adsorption on clinoptilolite compared to the other samples. The surface species affecting the q_iso and adsorption capacity were investigated through the FTIR spectroscopy using CO₂ as probe molecule. As a whole, it has been observed that CO₂ forms linear adducts onto K⁺ and Mg²⁺ cations of the bare clinoptilolite, and carbonate-like species onto its basic sites. With the Na-exchanged clinoptilolite, Na⁺ ions led to a decrease in surface basicity and to the formation of both single (Na⁺···OCO) and dual (Na⁺···OCO⋯Na⁺) cationic sites available for the formation of linear adducts. As a result of the remarkable adsorption capacity of clinoptilolite at 338 K, this material appears to be a promising adsorbent for the direct CO₂ removal from different flue gases sources operating at such temperatures.

High Throughput Methods in the Synthesis, Characterization, and Optimization of Porous Materials

Porous materials are widely employed in a large range of applications, in particular, for storage, separation, and catalysis of fine chemicals. Synthesis, characterization, and pre- and post-synthetic computer simulations are mostly carried out in a piecemeal and ad hoc manner. Whilst high throughput approaches have been used for more than 30 years in the porous material fields, routine integration of experimental and computational processes is only now becoming more established. Herein, important developments are highlighted and emerging challenges for the community identified, including the need to work toward more integrated workflows.

Molecular dynamics simulation of carbon dioxide diffusion in NaA zeolite: assessment of surface effects and evaluation of bulk-like properties

Molecular dynamics simulations were carried out for a finite sample of NaA zeolite in contact with bulk carbon dioxide in a wide range of temperatures and CO2 contents. Density and diffusion profiles were obtained to estimate the depth at which the external surfaces of the zeolite affect CO2 diffusion in porous space. The approximate depth of surface effects for NaA zeolite was estimated as ca. 2 nm, though this figure may vary depending on temperature and adsorbed gas density. Diffusion coefficients and diffusion activation energies were calculated for CO2 and Na+ in the bulk-like region of the zeolite. Diffusion activation energy for carbon dioxide demonstrated a non-monotonic dependence on the amount of adsorbed gas.

A Strong Test of Atomically Detailed Models of Molecular Adsorption in Zeolites Using Multilaboratory Experimental Data for CO2 Adsorption in Ammonium ZSM-5

A recent international interlaboratory study led by the U.S. National Institute of Standards (NIST) reported CO₂ adsorption isotherms measured independently by 11 groups on reference material RM 8852, an ammonium ZSM-5 zeolite. Good reproducibility and high reliability of this experimental data provide a strong test for the ability of atomically detailed models to predict adsorption of CO₂ in zeolites. We developed force fields for CO₂ in ammonium zeolites based on first-principles calculations and also independently performed experiments with RM 8852 by microcalorimetry. At low pressures good agreement was obtained between predictions and experiments. At high pressures, however, deviations were observed. We show that the charge-balancing cations in the experimental material are the predominant source of the discrepancy between simulation and experiment at high pressures; the experimental sample treatment causes deammoniation. In addition, accounting for a small amount of noncrystalline mesoporosity in the zeolite brings predictions into much better agreement with experiments.

Multi-factor study of the effects of a trace amount of water vapor on low concentration CO2 capture by 5A zeolite particles

CO₂ adsorption amount is enhanced with below 0.1 ppm humidity, and water molecule partial charge is a dominant factor in adsorption.

Efficient and accurate description of adsorption in zeolites

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.

Polarizable Force Field for CO2 in M-MOF-74 Derived from Quantum Mechanics

On the short term, carbon capture is a viable solution to reduce human-induced CO₂ emissions, which requires an energy efficient separation of CO₂. Metal-organic frameworks (MOFs) may offer opportunities for carbon capture and other industrially relevant separations. Especially, MOFs with embedded open metal sites have been shown to be promising. Molecular simulation is a useful tool to predict the performance of MOFs even before the synthesis of the material. This reduces the experimental effort, and the selection process of the most suitable MOF for a particular application can be accelerated. To describe the interactions between open metal sites and guest molecules in molecular simulation is challenging. Polarizable force fields have potential to improve the description of such specific interactions. Previously, we tested the applicability of polarizable force fields for CO₂ in M-MOF-74 by verifying the ability to reproduce experimental measurements. Here, we develop a predictive polarizable force field for CO₂ in M-MOF-74 (M = Co, Fe, Mg, Mn, Ni, Zn) without the requirement of experimental data. The force field is derived from energies predicted from quantum mechanics. The procedure is easily transferable to other MOFs. To incorporate explicit polarization, the induced dipole method is applied between the framework and the guest molecule. Atomic polarizabilities are assigned according to the literature. Only the Lennard-Jones parameters of the open metal sites are parameterized to reproduce energies from quantum mechanics. The created polarizable force field for CO₂ in M-MOF-74 can describe the adsorption well and even better than that in our previous work.

Fitting electron density as a physically sound basis for the development of interatomic potentials of complex alloys

A new method to obtain physically sound EAM parameters using the density functional theory electron density as the starting point.

Temperature-regulated guest admission and release in microporous materials

While it has long been known that some highly adsorbing microporous materials suddenly become inaccessible to guest molecules below certain temperatures, previous attempts to explain this phenomenon have failed. Here we show that this anomalous sorption behaviour is a temperature-regulated guest admission process, where the pore-keeping group's thermal fluctuations are influenced by interactions with guest molecules. A physical model is presented to explain the atomic-level chemistry and structure of these thermally regulated micropores, which is crucial to systematic engineering of new functional materials such as tunable molecular sieves, gated membranes and controlled-release nanocontainers. The model was validated experimentally with H₂, N₂, Ar and CH₄ on three classes of microporous materials: trapdoor zeolites, supramolecular host calixarenes and metal-organic frameworks. We demonstrate how temperature can be exploited to achieve appreciable hydrogen and methane storage in such materials without sustained pressure. These findings also open new avenues for gas sensing and isotope separation.

CO2 Capture and Separations Using MOFs: Computational and Experimental Studies

This Review focuses on research oriented toward elucidation of the various aspects that determine adsorption of CO₂ in metal-organic frameworks and its separation from gas mixtures found in industrial processes. It includes theoretical, experimental, and combined approaches able to characterize the materials, investigate the adsorption/desorption/reaction properties of the adsorbates inside such environments, screen and design new materials, and analyze additional factors such as material regenerability, stability, effects of impurities, and cost among several factors that influence the effectiveness of the separations. CO₂ adsorption, separations, and membranes are reviewed followed by an analysis of the effects of stability, impurities, and process operation conditions on practical applications.

Biomolecular Force Field Parameterization via Atoms-in-Molecule Electron Density Partitioning

Molecular mechanics force fields, which are commonly used in biomolecular modeling and computer-aided drug design, typically treat nonbonded interactions using a limited library of empirical parameters that are developed for small molecules. This approach does not account for polarization in larger molecules or proteins, and the parametrization process is labor-intensive. Using linear-scaling density functional theory and atoms-in-molecule electron density partitioning, environment-specific charges and Lennard-Jones parameters are derived directly from quantum mechanical calculations for use in biomolecular modeling of organic and biomolecular systems. The proposed methods significantly reduce the number of empirical parameters needed to construct molecular mechanics force fields, naturally include polarization effects in charge and Lennard-Jones parameters, and scale well to systems comprised of thousands of atoms, including proteins. The feasibility and benefits of this approach are demonstrated by computing free energies of hydration, properties of pure liquids, and the relative binding free energies of indole and benzofuran to the L99A mutant of T4 lysozyme.

DFT-Derived Force Fields for Modeling Hydrocarbon Adsorption in MIL-47(V)

Generic force fields such as UFF and DREIDING are widely used for predicting molecular adsorption and diffusion in metal-organic frameworks (MOFs), but the accuracy of these force fields is unclear. We describe a general framework for developing transferable force fields for modeling the adsorption of alkanes in a nonflexible MIL-47(V) MOF using periodic density functional theory (DFT) calculations. By calculating the interaction energies for a large number of energetically favorable adsorbate configurations using DFT, we obtain a force field that gives good predictions of adsorption isotherms, heats of adsorption, and diffusion properties for a wide range of alkanes and alkenes in MIL-47(V). The force field is shown to be transferable to related materials such as MIL-53(Cr) and is used to calculate the free-energy differences for the experimentally observed phases of MIL-53(Fe).

Advances in theory and their application within the field of zeolite chemistry

Zeolites are versatile and fascinating materials which are vital for a wide range of industries, due to their unique structural and chemical properties, which are the basis of applications in gas separation, ion exchange and catalysis. Given their economic impact, there is a powerful incentive for smart design of new materials with enhanced functionalities to obtain the best material for a given application. Over the last decades, theoretical modeling has matured to a level that model guided design has become within reach. Major hurdles have been overcome to reach this point and almost all contemporary methods in computational materials chemistry are actively used in the field of modeling zeolite chemistry and applications. Integration of complementary modeling approaches is necessary to obtain reliable predictions and rationalizations from theory. A close synergy between experimentalists and theoreticians has led to a deep understanding of the complexity of the system at hand, but also allowed the identification of shortcomings in current theoretical approaches. Inspired by the importance of zeolite characterization which can now be performed at the single atom and single molecule level from experiment, computational spectroscopy has grown in importance in the last decade. In this review most of the currently available modeling tools are introduced and illustrated on the most challenging problems in zeolite science. Directions for future model developments will be given.

Fine-tuning the theoretically predicted structure of MIL-47(V) with the aid of powder X-ray diffraction

The structural characterization of complex crystalline materials can be simplified by closely comparing theoretical and experimental diffraction patterns.

Quantitative prediction of physical properties of imidazolium based room temperature ionic liquids through determination of condensed phase site charges: a refined force field

Quantitative prediction of physical properties of room temperature ionic liquids through nonpolarizable force field based molecular dynamics simulations is a challenging task. The challenge lies in the fact that mean ion charges in the condensed phase can be less than unity due to polarization and charge transfer effects whose magnitude cannot be fully captured through quantum chemical calculations conducted in the gas phase. The present work employed the density-derived electrostatic and chemical (DDEC/c3) charge partitioning method to calculate site charges of ions using electronic charge densities obtained from periodic density functional theory (DFT) calculations of their crystalline phases. The total ion charges obtained thus range between -0.6e for chloride and -0.8e for the PF6 ion. The mean value of the ion charges obtained from DFT calculations of an ionic liquid closely matches that obtained from the corresponding crystal thus confirming the suitability of using crystal site charges in simulations of liquids. These partial charges were deployed within the well-established force field developed by Lopes et al., and consequently, parameters of its nonbonded and torsional interactions were refined to ensure that they reproduced quantum potential energy scans for ion pairs in the gas phase. The refined force field was employed in simulations of seven ionic liquids with six different anions. Nearly quantitative agreement with experimental measurements was obtained for the density, surface tension, enthalpy of vaporization, and ion diffusion coefficients.