Tailoring ammonia capture in MOFs and COFs: A multi-scale and machine learning comprehensive investigation of functional group modification

Our study aims to examine the impact of ligand functionalization on the ammonia adsorption properties of MOFs and COFs, by combining multi-scale calculations with machine learning techniques. Density Functional Theory calculations were performed to investigate the interactions between ammonia (NH3) and a comprehensive set of 48 strategically chosen functional groups. In all of the cases, it is observed that functionalized rings exhibit a stronger interaction with ammonia molecule compared to unfunctionalized benzene, while -O2Mg demonstrates the highest interaction energy with ammonia (15 times stronger than the bare benzene). The trend obtained from the thorough DFT screening is verified via Grand Canonical Monte-Carlo calculations by employing interatomic potentials derived from quantum chemical calculations. Isosteric heat of adsorption plots provide a comprehensive elucidation of the adsorption process, and important insights can be taken for studies in fine-tuning materials for ammonia adsorption. Furthermore, a proof of concept machine learning (ML) analysis is conducted, which demonstrates that ML can accurately predict NH3 binding energies despite the limited amount of data. The findings derived from our multi-scale methodology indicate that the functionalization strategy can be utilized to guide synthesis towards MOFs, COFs, or other porous materials for enhanced NH3 adsorption capacity.

Effect of CNT over structural properties of SAPO-34 in MTO process: Experimental and molecular simulation studies

The hierarchical silicoaluminophosphate (SAPO-34) catalyst was synthesized using the mixtures of diethylamine (D) and butylamine (B) as a structure-directing agent (SDA), and carbon nanotube (CNT) as a secondary template in the hydrothermal method. The catalysts were characterized by Fourier-transform infrared spectroscopy (FT-IR), X-ray diffraction (XRD), field emission scanning electron microscopy (FESEM), transmission electron microscopy (TEM), energy-dispersive X-ray spectroscopy (EDX), N2 physisorption, and temperature-programmed desorption of ammonia (NH3-TPD) techniques and evaluated for the catalytic activity in the Methanol to Olefins (MTO) process. The results showed that the use of CNT as the secondary template improved the hierarchical structure of SAPO-34 due to increasing the external surface area and mesoporosity and decreasing the particle size and as a result, made better the performance of the prepared SAPO-34 zeolite in the MTO process. Among all the prepared samples, the CNT-B-D catalyst synthesized by mixing three templates displayed the highest ethylene and propylene selectivity of 49% and 34%, respectively. Also, using CNT in the synthesis of samples increased the catalytic stability. In addition, pure, binary, and ternary adsorption isotherms and diffusivities of the main products and reactants over the SAPO-34 were investigated by theoretical measurements, because sorption and diffusion affect the catalyst stability and C2-C3 selectivity in the MTO reaction. The higher diffusion rate of ethylene leads to following the aromatic-based cycle in the MTO process.

Experimental, Grand canonical Monte Carlo, and Density Functional Theory studies for comparing halloysite and kaolinite on heavy metals fate during pyrolysis of solid waste

In this study, we utilized low-cost halloysite (Hal) for the first time to enhance the solid-phase enrichment and stability of heavy metals (HMs) during solid waste pyrolysis through experimental and theoretical methods, and compared with kaolinite (Kao). Experimental results demonstrated that Hal was superior to Kao in improving the solid-phase enrichment of HMs. Specifically, the solid-phase enrichment of Cd increased by 32.6 % (500 °C) and 25.94 % (600 °C), while that of Pb and Zn increased by 17.37 %/16.83 % and 19.82 %/22.37 % (700/800 °C), respectively. Adding Hal reduced the proportion of HMs in the unstable fraction (F1 + F2), consequently lowering the environmental risk of biochar and the extractable state of HMs. Through Grand canonical Monte Carlo and Density Functional Theory (DFT) simulations, we analyzed the adsorption amounts, adsorption sites, and adsorption mechanisms of Cd/Pb compound on Hal/Kao surfaces, revealing that the primary factor influencing the adsorption performance of Hal and Kao was the difference in specific surface area. The adsorption amounts of HMs by Hal were significantly higher than Kao and decreased with increasing temperature, while the difference in adsorption performance caused by structural bending was negligible. The DFT results indicated that Cd and Pb monomers were stabilized by establishing covalent bonds with OH or reactive O atoms on the Al-(001) surface, whereas the covalent bonds with ionic bonding properties formed between Cl atoms and unsaturated Al atoms played a crucial role in stabilizing HM chlorides. Furthermore, the adsorption energy of Hal on HMs increased with the removal rate of OH. Our study highlights the potential of Hal in stabilizing HMs during pyrolysis without requiring any modifications, thereby avoiding the generation of modified waste solutions and unnecessary cost loss.

Experimental and theoretical calculations insight into acetone adsorption by porous carbon at different pressures: Effects of pore structure and oxygen groups

The influence of different pore size and oxygen groups for porous carbons on acetone adsorption at different pressure was studied by using experimental data and theoretical calculation, and the results were applied to prepare carbon-based adsorbents with superior adsorption capacity. First, we successfully prepared five types of porous carbons with different gradient pore structure but similar oxygen contents (4.9 ± 0.25 at.%). We found that the acetone uptake at different pressure depends on the different pore sizes. Besides, we demonstrate how to accurately decompose the acetone adsorption isotherm into multiple sub-isotherms based on different pore sizes. Based on the isotherm decomposition method, the acetone adsorption at 18 kPa is mainly in the form of pore-filling adsorption in the pore size range of 0.6-2.0 nm. When the pore size is greater than 2 nm, the acetone uptake mainly depends on the surface area. Second, porous carbons with different oxygen content, similar surface area and pore structure were prepared to study the influence of oxygen groups on acetone adsorption. The results show that the acetone adsorption capacity is determined by the pore structure at relatively high pressure, and the oxygen groups only slightly increase the adsorption capacity. However, the oxygen groups can provide more active sites, thereby enhancing acetone adsorption at low pressure.

Adsorption and diffusion of the H2/CO2/CO/MeOH/EtOH mixture into the ZIF-7 using molecular simulation

The adsorption and diffusion of synthesis gas components (methanol, ethanol, H₂, CO_2, and CO molecules) in ZIF-7 by grand canonical Monte Carlo and molecular dynamics simulation were investigated. The initial diffusion coefficient at the beginning of the process depends on the kinetic diameter of the guest molecules. Also, the diffusion coefficient at equilibrium conditions probably depends on the interaction between the guest molecules with the ZIF-7 framework. The radial distribution function results indicate that the distribution of guest molecules in the framework is affected by the interaction between the guest molecules. These results indicate that the CO, CO₂, and H₂ guest molecules are adsorbed on both the Zn metal atom and the organic linker (especially the C₁ atom). In contrast, the organic linker is the most favorable adsorption site for methanol and ethanol guest molecules. In addition, the diffusion coefficient of guest molecules in binary mixtures is related to the attraction or repulsion between the guest molecules.

Enhancing sampling of water rehydration upon ligand binding using variants of grand canonical Monte Carlo

Water plays an important role in mediating protein-ligand interactions. Water rearrangement upon a ligand binding or modification can be very slow and beyond typical timescales used in molecular dynamics (MD) simulations. Thus, inadequate sampling of slow water motions in MD simulations often impairs the accuracy of the accuracy of ligand binding free energy calculations. Previous studies suggest grand canonical Monte Carlo (GCMC) outperforms normal MD simulations for water sampling, thus GCMC has been applied to help improve the accuracy of ligand binding free energy calculations. However, in prior work we observed protein and/or ligand motions impaired how well GCMC performs at water rehydration, suggesting more work is needed to improve this method to handle water sampling. In this work, we applied GCMC in 21 protein-ligand systems to assess the performance of GCMC for rehydrating buried water sites. While our results show that GCMC can rapidly rehydrate all selected water sites for most systems, it fails in five systems. In most failed systems, we observe protein/ligand motions, which occur in the absence of water, combine to close water sites and block instantaneous GCMC water insertion moves. For these five failed systems, we both extended our GCMC simulations and tested a new technique named grand canonical nonequilibrium candidate Monte Carlo (GCNCMC). GCNCMC combines GCMC with the nonequilibrium candidate Monte Carlo (NCMC) sampling technique to improve the probability of a successful water insertion/deletion. Our results show that GCNCMC and extended GCMC can rehydrate all target water sites for three of the five problematic systems and GCNCMC is more efficient than GCMC in two out of the three systems. In one system, only GCNCMC can rehydrate all target water sites, while GCMC fails. Both GCNCMC and GCMC fail in one system. This work suggests this new GCNCMC method is promising for water rehydration especially when protein/ligand motions may block water insertion/removal.

Insights into the Gas Adsorption Mechanisms in Metal-Organic Frameworks from Classical Molecular Simulations

Classical molecular simulations can provide significant insights into the gas adsorption mechanisms and binding sites in various metal-organic frameworks (MOFs). These simulations involve assessing the interactions between the MOF and an adsorbate molecule by calculating the potential energy of the MOF-adsorbate system using a functional form that generally includes nonbonded interaction terms, such as the repulsion/dispersion and permanent electrostatic energies. Grand canonical Monte Carlo (GCMC) is the most widely used classical method that is carried out to simulate gas adsorption and separation in MOFs and identify the favorable adsorbate binding sites. In this review, we provide an overview of the GCMC methods that are normally utilized to perform these simulations. We also describe how a typical force field is developed for the MOF, which is required to compute the classical potential energy of the system. Furthermore, we highlight some of the common analysis techniques that have been used to determine the locations of the preferential binding sites in these materials. We also review some of the early classical molecular simulation studies that have contributed to our working understanding of the gas adsorption mechanisms in MOFs. Finally, we show that the implementation of classical polarization for simulations in MOFs can be necessary for the accurate modeling of an adsorbate in these materials, particularly those that contain open-metal sites. In general, molecular simulations can provide a great complement to experimental studies by helping to rationalize the favorable MOF-adsorbate interactions and the mechanism of gas adsorption.

Gas adsorption in Mg-porphyrin-based porous organic frameworks: A computational simulation by first-principles derived force field

A novel type of porous organic frameworks, based on Mg-porphyrin, with diamond-like topology, named POF-Mgs is computationally designed, and the gas uptakes of CO₂ , H₂ , N₂ , and H₂ O in POF-Mgs are investigated by Grand canonical Monte Carlo simulations based on first-principles derived force fields (FF). The FF, which describes the interactions between POF-Mgs and gases, are fitted by dispersion corrected double-hybrid density functional theory, B2PLYP-D3. The good agreement between the obtained FF and the first-principle energies data confirms the reliability of the FF. Furthermore our simulation shows the presence of a small amount of H₂ O (≤ 0.01 kPa) does not much affect the adsorption quantity of CO₂ , but the presence of higher partial pressure of H₂ O (≥ 0.1 kPa) results in the CO₂ adsorption decrease significantly. The good performance of POF-Mgs in the simulation inspires us to design novel porous materials experimentally for gas adsorption and purification. © 2017 Wiley Periodicals, Inc.