Active learning of alchemical adsorption simulations; towards a universal adsorption model

Adsorption is a fundamental process studied in materials science and engineering because it plays a critical role in various applications, including gas storage and separation. Understanding and predicting gas adsorption within porous materials demands comprehensive computational simulations that are often resource intensive, limiting the identification of promising materials. Active learning (AL) methods offer an effective strategy to reduce the computational burden by selectively acquiring critical data for model training. Metal-organic frameworks (MOFs) exhibit immense potential across various adsorption applications due to their porous structure and their modular nature, leading to diverse pore sizes and chemistry that serve as an ideal platform to develop adsorption models. Here, we demonstrate the efficacy of AL in predicting gas adsorption within MOFs using "alchemical" molecules and their interactions as surrogates for real molecules. We first applied AL separately to each MOF, reducing the training dataset size by 57.5% while retaining predictive accuracy. Subsequently, we amalgamated the refined datasets across 1800 MOFs to train a multilayer perceptron (MLP) model, successfully predicting adsorption of real molecules. Furthermore, by integrating MOF features into the AL framework using principal component analysis (PCA), we navigated MOF space effectively, achieving high predictive accuracy with only a subset of MOFs. Our results highlight AL's efficiency in reducing dataset size, enhancing model performance, and offering insights into adsorption phenomenon in large datasets of MOFs. This study underscores AL's crucial role in advancing computational material science and developing more accurate and less data intensive models for gas adsorption in porous materials.

High throughput screening of pure silica zeolites for CF4 capture from electronics industry gas

The capture and separation of CF4 from CF4/N2 mixture gas is a crucial issue in the electronics industry, as CF4 is a commonly used etching gas and the ratio of CF4 to N2 directly affects process efficiency. Utilizing high-throughput computational screening techniques and grand canonical Monte Carlo (GCMC) simulations, we comprehensively screened and assessed 247 types of pure silicon zeolite materials to determine their adsorption and separation performance for CF4/N2 mixtures. Based on screening, the relationships between the structural parameters and adsorption and separation properties were meticulously investigated. Four indicators including adsorption selectivity, working capacity, adsorbent performance score (APS), and regenerability (R%) were used to evaluate the performance of adsorbents. Based on the evaluation, we selected the top three best-performing zeolite structures for vacuum swing adsorption (LEV, AWW and ESV) and pressure swing adsorption (AVL, ZON, and ERI) processes respectively. Also, we studied the preferable adsorption sites of CF4 and N2 in the selected zeolite structures through centroid density distributions at the molecule level. We expect the study may provide some valuable guidance for subsequent experimental investigations on adsorption and separation of CF4/N2.

Conversion of waste cork to N-doped porous carbons by urea-assisted hydrothermal method for enhanced VOC capture

Converting waste resources into porous carbon for pollutants capture is an effective strategy to achieve the environmental goal of "treating waste with waste". Cork is an ideal precursor of porous carbons due to its ordered honeycomb-like cell structure and layered composition distribution. Herein, N-doped porous carbons (PCs) were prepared via two steps of urea-assisted hydrothermal carbonization and chemical activation to mitigate volatile organic compounds (VOCs) pollution. Results indicated that the obtained PC4-800 exhibited remarkable features for adsorption including high total pore volume (0.97 cm3/g) and specific surface area (1864.89 m2/g), as well as abundant N-containing functional groups. The excellent pore structure was primarily owing to the corrosion of the carbon matrix by the gas produced from the reaction of K2CO3 and N-containing functional groups. The adsorption results showed that the PC4-800 have an outstanding toluene adsorption capacity (867.03 mg/g) that outperforming majority of adsorbents previously reported. There are substantial pores in N-doped PCs with a pore width of 1.71-2.28 nm, which is 3 to 4 times the molecular dynamic diameter of toluene, and plays a crucial role in the absorption process. Moreover, the promotional influence of N-functional groups on the toluene adsorption process was verified through DFT calculation by Gaussian imitating, where N-6 generated π-electron enrichment sites on the surface of N-doped PCs, facilitating π-π dispersion with the benzene ring in toluene. This study provides a new strategy to convert waste cork into high-performance adsorbents for VOCs removal.

Atomistic Insights into the Effect of Functional Groups on the Adsorption of Water by Activated Carbon for Heat Energy Storage

Adsorption heat storage holds great promise for solar energy applications. The development of new adsorbent materials is currently the research focus in this area. The present work designs several activated carbon models with different functional groups, including -OH, -NH2, -COOH, and -SO3H, and explores the influence of functional groups' categories and numbers on the water adsorption capacity of the activated carbon using the GCMC method. The adsorption mechanism between functional groups and water molecules is analyzed using density functional theory. The results show that the functional groups could significantly improve the water adsorption capacity of activated carbon due to the hydrogen bond between functional groups and water molecules. In the scope of this paper, under low pressure, the activated carbon with -SO3H exhibits the best adsorption capacity, followed by the activated carbon with -COOH. Under low and medium pressure, increasing the number of -SO3H functional groups could increase the water adsorption capacity; however, when the pressure is high, increasing the functional group numbers might decrease the water adsorption capacity. As the temperature increases, the water adsorption capacity of activated carbons decreases, and the activated carbon with -SO3H is proven to have excellent application prospects in heat energy storage.

Judicious design functionalized 3D-COF to enhance CO2 adsorption and separation

The effects of functional groups (including OH, OCH₃ , NH₂ , CH₂ NH₂ , COOH, SO₃ H, OCO(CH₂ )₂ COOH(E-COOH), and (CH₂ )₄ COOH(c-COOH)) in 3D covalent organic frameworks (3D-COFs) on CO₂ adsorption and separation are investigated by grand canonical Monte Carlo (GCMC) simulations and density functional theory calculations. The results indicate that interaction between CO₂ and the framework is the main factor for determining CO₂ uptakes at low pressure, while pore size becomes the decisive factor at high pressure. The binding energy of CO₂ with functionalized linker is correlated to CO₂ uptake at 0.3 bar and 298 K on 3D-COF-1, suggesting functional groups play a key role in CO₂ capture in microporous 3D-COFs. Moreover, CO₂ selectivity over CH₄ , N₂ , and H₂ can be significantly enhanced by functionalization, where CH₂ NH₂ , COOH, SO₃ H, and E-COOH enhance CO₂ adsorption more effectively at 1 bar. Among them, SO₃ H is the most promising functional group in 3D-COFs for CO₂ separation.

H2O-prompted CO2 capture on metal silicates in situ generated from SBA-15

A series of metal silicates, NaMSi₁₀Ox (M = Cu, Mn and Ni), were prepared by in situ doping of metals into mesoporous SBA-15 under a hydrothermal process, displaying a continuous framework of SiO₄ structure with a narrow pore size distribution. These metal silicate materials were tested for CO₂ adsorption behavior in the absence and presence of water. The results exhibited that the effect of H₂O on the CO₂ capture capability of metal silicates depends on the types of metal inserted into SBA-15. Compared to the dry condition, H₂O addition enhances CO₂ uptake dramatically for NaCuSi₁₀Ox by 25%, and slightly for NaNiSi₁₀Ox (∼10%), whereas little effect is shown on NaMnSi₁₀Ox. The metal silicate materials are stable after adsorption of CO₂ under wet conditions, which is benefited from their synthesis method, hydrothermal conditions. The improvement of CO₂ uptake on metal silicates by H₂O is attributed to the competitive and synergistic adsorption mechanism on the basis of IR investigations, where initially adsorbed H₂O acts as a promoter for further CO₂ capture through a hydration reaction, i.e., formation of bicarbonate and carbonates on the surface of the samples. These observations provide new possibilities for the design and synthesis of porous metal silicate materials for CO₂ capture under practical conditions where moisture is present.

Porous metal silicates prepared by an in situ doping strategy of metals into SBA-15 under hydrothermal conditions display efficient CO₂ capture performances in the absence and presence of moisture.

Separation of CF4/N2, C2F6/N2, and SF6/N2 Mixtures in Amorphous Activated Carbons Using Molecular Simulations

The capture and separation of CF₄, C₂F₆, and SF₆ and their mixtures containing nitrogen is a challenging process. To solve this, we propose the use of saccharose coke-based carbons as membranes for the adsorption and separation of these gases. By means of advanced techniques of Monte Carlo and molecular dynamics simulations, we have studied the adsorption and diffusion of CF₄, C₂F₆, and SF₆ as well as their mixtures with nitrogen in three HRMC carbon models, namely, CS400, CS1000, and CS1000a. We have computed the adsorption isotherms of the single components and the heat of adsorption as a function of the adsorbed concentration. We have also calculated the competitive adsorption of fluoride molecules and nitrogen at two different molar fractions, 0.1 and 0.9. We have computed the transport properties of the adsorbed gases in terms of the self-diffusivities and corrected diffusivities. The performance of the membranes for the targeted separations has been characterized by the calculation of the permselectivity. Our results indicate that the activated amorphous carbon CS1000a is an efficient adsorbent for the capture of the fluoride adsorbates as well as their purification from nitrogen-based mixtures.

Effects of Intrinsic Flexibility on Adsorption Properties of Metal-Organic Frameworks at Dilute and Nondilute Loadings

Molecular simulation of adsorption in nanoporous materials has become a valuable complement to experimental studies of these materials. In almost all cases, these simulations treat the adsorbing material as rigid. We use molecular simulations to examine the validity of this approximation for the adsorption in metal-organic frameworks (MOFs) that have framework flexibility without change in their unit cells because of thermal vibrations. All nanoporous materials are subject to this kind of framework flexibility. We examine the adsorption of nine molecules (CO₂, CH₄, ethane, ethene, propane, propene, butane, Xe, and Kr) and four molecular mixtures (CO₂/CH₄, ethane/ethene, propane/propene/butane, and Xe/Kr) in 100 MOFs at dilute and nondilute adsorption conditions. Our results show that single-component adsorption uptakes at nondilute conditions are only weakly affected by framework flexibility, but adsorption selectivities at both dilute and nondilute conditions can be significantly affected by flexibility. The most dramatic impacts of framework flexibility occur for adsorption uptake in the limit of dilute adsorption. These results suggest that the importance of including framework flexibility when attempting to make quantitative predictions of adsorption selectivity in MOFs and similar materials may have been underestimated in the past.

Solid-Waste-Derived Carbon Dioxide-Capturing Materials

Solid sorbents are considered to be promising materials for carbon dioxide capture. In recent years, many studies have focused on the use of solid waste as carbon dioxide sorbents. The use of waste resources as carbon dioxide sorbents not only leads to the development of relatively low-cost materials, but also eliminates waste simultaneously. Different types of waste materials from biomass, industrial waste, household waste, and so forth were used as carbon dioxide sorbents with sufficient carbon dioxide capture capacities. Herein, progress on the development of carbon dioxide sorbents produced from waste materials is reviewed and covers key factors, such as the type of waste, preparation method, further modification method, carbon dioxide sorption performance, and kinetics studies. In addition, a new research direction for further study is proposed. It is hoped that this critical review will not merely sum up the major research directions in this field, but also provide significant suggestions for future work.

Experimental Investigation of the Adsorption Characteristics of Mixed Coal and Variations of Specific Surface Areas before and after CH4 Adsorption

Tectonic coal is a kind of soft coal that is generated during tectonic movement. Gas outbursts usually occur in seams containing both virgin coal and tectonic coal. To reveal the adsorption characteristics of this type of coal seam (containing both virgin coal and tectonic coal), both tectonic coal and virgin coal were collected from the same longwall face and a series of laboratory tests were conducted, including coal sorption tests and pore specific surface measurements. Both the tectonic coal and virgin coal were crushed into coal powder (0.18–0.25 mm) for the coal sorption tests. In these laboratory tests, different mass ratios between tectonic coal and virgin coal were tested. We found that with the increase of the percentage of tectonic coal, the adsorption volume showed a rising trend, reached its maximum value, and then decreased. The specific surface areas of the mixed coal samples had the same evolution trends as those of the adsorption volume. From the laboratory tests, we found that when the mass ratio of virgin coal to tectonic coal was 1:1, both the adsorption volume and the specific surface areas reached their maximum values. Due to the percentage variation of the tectonic coal in the panel with the advancement of the longwall face, when the tectonic coal accounted for 50% of the total coal, the gas content would rise. Thus, proper measures should be adopted for outburst hazards control. The mathematical model between the change of specific surface area and the stress and strain of pore expansion before and after gas adsorption was established, and the relationship between the change of pore structure and gas emission before and after gas adsorption was obtained. It provides a theoretical basis for further research on coal and gas outburst mechanisms.

Adsorption and Diffusion of Benzene in Mg-MOF-74 with Open Metal Sites

We performed molecular simulations to investigate the adsorption and diffusion of benzene in metal-organic framework Mg-MOF-74. At 300 K and 20 Pa, the saturated loading of benzene reaches 8.2 mmol/g, almost twice of (12,12) single-walled carbon nanotube with a similar pore size, and 93% of the benzene molecules in Mg-MOF-74 can desorb at 390 K. The energy analysis indicates that the van der Waals contribution still dominates 70-80% of the total fluid-wall interaction energy compared with the Coulombic contribution. We further analyzed the structure of benzene confined in Mg-MOF-74 by the molecular snapshots, pair correlation functions, orientational order parameters, and local density profiles. It is found that low temperature and high pressure make the structure of adsorbed benzene more similar to that of the liquid benzene. Moreover, the benzene molecules in the contact adsorption layer lie flat on the surface of adsorbent, whereas those molecules near the pore center have no particular orientations. Due to the existence of open metal sites, the structures of adsorbed benzene are more compact and ordered than those of bulk liquid benzene. Consequently, the self-diffusion coefficient of saturated benzene in Mg-MOF-74 at 300 K is significantly lower than that of bulk liquid benzene and confined liquid benzene in slit pores and disordered carbons by 4-5 orders of magnitude. We investigated the separation and diffusion of benzene/cyclohexane in the mixture in Mg-MOF-74 and found that the pores almost completely adsorbed benzene, although its self-diffusion coefficient was slightly lower than that of cyclohexane.

Chemo-mechanical coupling in kerogen gas adsorption/desorption

Kerogen plays a central role in hydrocarbon generation in an oil/gas reservoir. In a subsurface environment, kerogen is constantly subjected to stress confinement or relaxation. The interplay between mechanical deformation and gas adsorption of the materials could be an important process for shale gas production but unfortunately is poorly understood. Using a hybrid Monte Carlo/molecular dynamics simulation, we show here that a strong chemo-mechanical coupling may exist between gas adsorption and mechanical strain of a kerogen matrix. The results indicate that the kerogen volume can expand by up to 5.4% and 11% upon CH4 and CO2 adsorption at 192 atm, respectively. The kerogen volume increases with gas pressure and eventually approaches a plateau as the kerogen becomes saturated. The volume expansion appears to quadratically increase with the amount of gas adsorbed, indicating a critical role of the surface layer of gas adsorbed in the bulk strain of the material. Furthermore, gas uptake is greatly enhanced by kerogen swelling. Swelling also increases the surface area, porosity, and pore size of kerogen. Our results illustrate the dynamic nature of kerogen, thus questioning the validity of the current assumption of a rigid kerogen molecular structure in the estimation of gas-in-place for a shale gas reservoir or gas storage capacity for subsurface carbon sequestration. The coupling between gas adsorption and kerogen matrix deformation should be taken into consideration.

Ca-Embedded C2N: an efficient adsorbent for CO2 capture

Carbon dioxide as a greenhouse gas causes severe impacts on the environment, whereas it is also a necessary chemical feedstock that can be converted into carbon-based fuels via electrochemical reduction. To efficiently and reversibly capture CO₂, it is important to find novel materials for a good balance between adsorption and desorption. In this study, we performed first-principles calculations and grand canonical Monte Carlo (GCMC) simulations, to systematically study metal-embedded carbon nitride (C₂N) nanosheets for CO₂ capture. Our first-principles results indicated that Ca atoms can be uniformly trapped in the cavity center of C₂N structure, while the transition metals (Sc, Ti, V, Cr, Mn, Fe, Co) are favorably embedded in the sites off the center of the cavity. The determined maximum number of CO₂ molecules with strong physisorption showed that Ca-embedded C₂N monolayer is the most promising CO₂ adsorbent among all considered metal-embedded materials. Moreover, GCMC simulations revealed that at room temperature the gravimetric density for CO₂ adsorbed on Ca-embedded C₂N reached 50 wt% at 30 bar and 23 wt% at 1 bar, higher than other layered materials, thus providing a satisfactory system for the CO₂ capture and utilization.

Carbon flakes based metal organic frameworks for H2, CH4 and CO2 gas storage: a GCMC simulation study

In this study, new metal organic frameworks have been designed by modifying the linker of IRMOF-1 with different carbon flakes (circular, rectangular and rhombus).

Porous Carbon Materials Based on Graphdiyne Basis Units by the Incorporation of the Functional Groups and Li Atoms for Superior CO2 Capture and Sequestration

The graphdiyne family has attracted a high degree of concern because of its intriguing and promising properties. However, graphdiyne materials reported to date represent only a tiny fraction of the possible combinations. In this work, we demonstrate a computational approach to generate a series of conceivable graphdiyne-based frameworks (GDY-Rs and Li@GDY-Rs) by introducing a variety of functional groups (R = -NH₂, -OH, -COOH, and -F) and doping metal (Li) in the molecular building blocks of graphdiyne without restriction of experimental conditions and rapidly screen the best candidates for the application of CO₂ capture and sequestration (CCS). The pore topology and morphology and CO₂ adsorption and separation properties of these frameworks are systematically investigated by combining density functional theory (DFT) and grand canonical Monte Carlo (GCMC) simulations. On the basis of our computer simulations, combining Li-doping and hydroxyl groups strategies offer an unexpected synergistic effect for efficient CO₂ capture with an extremely CO₂ uptake of 4.83 mmol/g at 298 K and 1 bar. Combined with its superior selectivity (13 at 298 K and 1 bar) for CO₂ over CH₄, Li@GDY-OH is verified to be one of the most promising materials for CO₂ capture and separation.

Nanostructural control of methane release in kerogen and its implications to wellbore production decline

Despite massive success of shale gas production in the US in the last few decades there are still major concerns with the steep decline in wellbore production and the large uncertainty in a long-term projection of decline curves. A reliable projection must rely on a mechanistic understanding of methane release in shale matrix-a limiting step in shale gas extraction. Using molecular simulations, we here show that methane release in nanoporous kerogen matrix is characterized by fast release of pressurized free gas (accounting for ~30-47% recovery) followed by slow release of adsorbed gas as the gas pressure decreases. The first stage is driven by the gas pressure gradient while the second stage is controlled by gas desorption and diffusion. We further show that diffusion of all methane in nanoporous kerogen behaves differently from the bulk phase, with much smaller diffusion coefficients. The MD simulations also indicate that a significant fraction (3-35%) of methane deposited in kerogen can potentially become trapped in isolated nanopores and thus not recoverable. Our results shed a new light on mechanistic understanding gas release and production decline in unconventional reservoirs. The long-term production decline appears controlled by the second stage of gas release.

Extraction of Lanthanide and Actinide Ions from Aqueous Mixtures Using a Carboxylic Acid-Functionalized Porous Aromatic Framework

Porous aromatic frameworks (PAFs) incorporating a high concentration of acid functional groups possess characteristics that are promising for use in separating lanthanide and actinide metal ions, as required in the treatment of radioactive waste. These materials have been shown to be indefinitely stable to concentrated acids and bases, potentially allowing for multiple adsorption/stripping cycles. Additionally, the PAFs combine exceptional features from MOFs and inorganic/activated carbons giving rise to tunable pore surfaces and maximum chemical stability. Herein, we present a study of the adsorption of selected metal ions, Sr(2+), Fe(3+), Nd(3+), and Am(3+), from aqueous solutions employing a carbon-based porous aromatic framework, BPP-7 (Berkeley Porous Polymer-7). This material displays high metal loading capacities together with excellent adsorption selectivity for neodymium over strontium based on Langmuir adsorption isotherms and ideal adsorbed solution theory (IAST) calculations. Based in part upon X-ray absorption spectroscopy studies, the stronger adsorption of neodymium is attributed to multiple metal ion and binding site interactions resulting from the densely functionalized and highly interpenetrated structure of BPP-7. Recyclability and combustibility experiments demonstrate that multiple adsorption/stripping cycles can be completed with minimal degradation of the polymer adsorption capacity.

Adsorption and separation of binary and ternary mixtures of SO2, CO2 and N2 by ordered carbon nanotube arrays: grand-canonical Monte Carlo simulations

The selectivity value is found here to range from 4 to 16, indicating that the optimization of pore size tuning can increase it by 4 times.

Luminescent zinc metal–organic framework (ZIF-90) for sensing metal ions, anions and small molecules

We synthesize ZIF-90, which exhibits an intense blue luminescence excited under visible light. Luminescent studies indicate that ZIF-90 could be an efficient multifunctional fluorescence material for sensing metal cations, anions and organic small molecules.

Interaction between CO2 and ionic liquids confined in the nanopores of SAPO-11

TPD-CO₂ profiles of P_4,4,4,6Triz confined in the nanopores of SAPO-11 (30 wt%).

Synthesis optimization of the ultra-microporous [Ni3(HCOO)6] framework to improve its CH4/N2separation selectivity

The adsorptive separation selectivities of CH₄/N₂were successfully improved from 4.0–4.8 to 7.0–7.5viasynthesis optimization of ultra-microporous [Ni₃(HCOO)₆] frameworks.