Search

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance. [ABSTRACT FROM AUTHOR]

Recent experiments show that the diffusivities of C8 and C10 n‐alkanes in ZIF‐8 are higher than those of C7 and C9, respectively. We investigated this unusual 'odd‐even' effect by simulating the adsorption of C7‐C14 n‐alkanes in ZIF‐8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest‐host energies, isosteric heats, and chain length distributions are analyzed for trends among the n‐alkanes. ZIF‐8 cages filled with n‐alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. Results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments. [ABSTRACT FROM AUTHOR]

While multiple external stimuli (e.g., temperature, light, pressure) have been reported to regulate gas adsorption, limited studies have been conducted on controlling molecular admission in nanopores through the application of electric fields (E-field). Here we show gas adsorption capacity and selectivity in zeolite molecular sieves can be regulated by an external E-field. Through E-field pre-activation during degassing, several zeolites exhibited enhanced CO2 adsorption and decreased CH4 and N2 adsorptions, improving the CO2/CH4 and CO2/N2 separation selectivity by at least 25%. The enhanced separation performance of the zeolites pre-activated by E-field was maintained in multiple adsorption/desorption cycles. Powder X-ray diffraction analysis and ab initio computational studies revealed that the cation relocation and framework expansion induced by the E-field accounted for the changes in gas adsorption capacities. These findings demonstrate a regulation approach to sharpen the molecular sieving capability by E-fields and open new avenues for carbon capture and molecular separations. The gas adsorption capacity and selectivity in zeolite molecular sieves can be regulated by an external electric field through the electric field-induced cation relocation and framework expansion. [ABSTRACT FROM AUTHOR]

A machine learning (ML) methodology that uses a histogram of interaction energies has been applied to predict gas adsorption in metal–organic frameworks (MOFs) using results from atomistic grand canonical Monte Carlo (GCMC) simulations as training and test data. In this work, the method is first extended to binary mixtures of spherical species, in particular, Xe and Kr. In addition, it is shown that single-component adsorption of ethane and propane can be predicted in good agreement with GCMC simulation using a histogram of the adsorption energies felt by a methyl probe in conjunction with the random forest ML method. The results for propane can be improved by including a small number of MOF textural properties as descriptors. We also discuss the most significant features, which provides physical insight into the most beneficial adsorption energy sites for a given application. [ABSTRACT FROM AUTHOR]

We present a fast and accurate, semi-analytical method for predicting hydrogen adsorption in nanoporous materials. For any temperature and pressure, the adsorbed amount is calculated as an integral over the energy density of adsorption sites (guest-host interactions) plus an average guest-guest term. The guest-host interaction energy is calculated using a classical force field with hydrogen modelled as a single-site probe. The guest-guest interaction energy is approximated using an average coordination number, which is regressed using Gaussian Process Regression (GPR). Local adsorption at each site is then modelled using a Langmuir isotherm, which when weighted with its probability density gives an accurate description of hydrogen adsorption. The method is tested on 933 metal-organic frameworks (MOFs) from the Computation-Ready Experimental (CoRE) MOF database at 77 K from 10 − 5 to 100 bar, and the results are compared against GCMC predictions. To demonstrate the utility of the method, we calculated hydrogen adsorption isotherms for 12,914 existing MOF structures, at two different temperatures at a speed about 100 times that of GCMC simulations and analyzed the results. We found 13 MOFs with predicted deliverable capacities exceeding the DOE target of 50 g/L for adsorption at 100 bar, 77 K and desorption at 5 bar, 160 K. [ABSTRACT FROM AUTHOR]

Density functional theory was used to study interactions between phosphine and 17 functional groups, revealing that only the functional groups containing metal ions bind phosphine strongly. This led to a detailed examination of the adsorption of the pnictogen hydride gases ammonia, phosphine, and arsine on 33 metal catecholates in order to examine periodic trends in the binding strength for these adsorbates. Phosphine and arsine adsorption is primarily driven by donation of electron density from the adsorbate molecule to the metal atom, while ammonia binding involves both electron donation and a significant Coulomb attraction between the negatively charged N atom and the positive metal atom. This Coulomb effect results in notably different binding behavior for ammonia than for phosphine and arsine, which are quite similar to each other. Generally, for metals on the left side of the periodic table, the Coulomb effect dominates in ammonia adsorption, and metals on the right side of the periodic table have lower positive charge and accept greater amounts of electron density, making the electron-sharing contribution more important to the binding strength. We find that more electronegative metals, particularly Pd, Ir, Pt, and Au, which accept more electron density from the adsorbate, yield the best selectivity for the target molecules over water because water adsorption is mostly due to Coulomb attraction. We also analyze changes in the d orbital occupancy upon ammonia and phosphine adsorption. [ABSTRACT FROM AUTHOR]

A new software package, RASPA, for simulating adsorption and diffusion of molecules in flexible nanoporous materials is presented. The code implements the latest state-of-the-art algorithms for molecular dynamics and Monte Carlo (MC) in various ensembles including symplectic/measure-preserving integrators, Ewald summation, configurational-bias MC, continuous fractional component MC, reactive MC and Baker's minimisation. We show example applications of RASPA in computing coexistence properties, adsorption isotherms for single and multiple components, self- and collective diffusivities, reaction systems and visualisation. The software is released under the GNU General Public License. [ABSTRACT FROM PUBLISHER]

The continuous fractional component Monte Carlo method, which was designed to overcome difficulties with insertions and deletions of molecules, is modified to include configurational bias Monte Carlo methods and is further extended to binary systems. The modified method is shown to correctly predict adsorption of Ar in silicalite, Xe and Kr in HKUST-1, and enantiomers in a homochiral metal–organic framework. The modified method is also found to be approximately an order of magnitude more efficient in inserting and deleting molecules than traditional configurational bias grand canonical Monte Carlo simulations in dense systems. [ABSTRACT FROM AUTHOR]

We designed, synthesized, and characterized a new Zr-based metal-organic framework material, NU-1100, with a pore volume of 1.53 ccg−1 and Brunauer-Emmett-Teller (BET) surface area of 4020 m2g−1; to our knowledge, currently the highest published for Zr-based MOFs. CH4/CO2/H2 adsorption isotherms were obtained over a broad range of pressures and temperatures and are in excellent agreement with the computational predictions. The total hydrogen adsorption at 65 bar and 77 K is 0.092 g g−1, which corresponds to 43 g L−1. The volumetric and gravimetric methane-storage capacities at 65 bar and 298 K are approximately 180 vSTP/v and 0.27 g g−1, respectively. [ABSTRACT FROM AUTHOR]

Simulated adsorption isotherms for water in UiO-66 illustrate that defects in the form of missing linkers make this MOF more hydrophilic. Heats of adsorption and density plots further confirm the effect of defects on adsorption of water in UiO-66 at low loadings. In addition, water and CO2 isotherms indicate that not only the amount of defects but their locations within the material affect the loading of guest molecules. [ABSTRACT FROM AUTHOR]

Competitive coadsorption of water is a major problem in the deployment of adsorption-based CO2 capture. Water molecules may compete for adsorption sites, reducing the capacity of the material, and dehumidification prior to separating CO2 from N2 increases process complexity and cost. The development of adsorbent materials that can selectively adsorb CO2 in the presence of water would be a major step forward in the deployment of CO2 capture materials in practice. In this study, large-scale computational screening was carried out to search for metal-organic frameworks (MOFs) with high selectivity toward CO2 over H2O. Calculating framework charges for thousands of MOFs is a significant challenge, so initial screening used a fast, but approximate, charge calculation method. On the basis of the initial screening, 15 MOFs were selected, and Monte Carlo simulations were carried out to compute the adsorption isotherms for these MOFs using more accurate framework charges calculated by density functional theory. A detailed investigation was performed on the effect of using different methods for calculating partial charges, and it was found that electrostatic interactions contribute the majority of the adsorption energy of H2O in the selected MOFs. [ABSTRACT FROM AUTHOR]

Diffusion of methane, ethane, propane and n-butane was studied within the micropores of several metal-organic frameworks (MOFs) of varying topologies, including the MOFs PCN-14, NU-125, NU-1100 and DUT-49. Diffusion coefficients of the pure components, as well as methane/ethane, methane/propane and methane/butane binary mixtures, were calculated using molecular dynamics simulations to understand the effect of the longer alkanes on uptake of natural gas in MOFs. The calculated self-diffusion coefficients of all four components are on the order of 10 −8 m 2 /s. The diffusion coefficients of the pure components decrease as a function of chain length in all of the MOFs studied and show different behaviour as a function of loading in different MOFs. The self-diffusivities follow the trend D PCN-14 NU-125≈NU-1100>DUT-49. By comparing the diffusion of pure methane and methane mixtures with the higher alkanes, it is observed that the diffusivity of methane is unaffected by the presence of the higher alkanes in the MOFs considered, indicating that the diffusion path of methane is not blocked by the higher alkanes present in natural gas. [ABSTRACT FROM AUTHOR]

Computational high-throughput screening was carried out to assess a large number of experimentally reported metal-organic frameworks (MOFs) and zeolites for their utility in hexane isomer separation. Through the work, we identified many MOFs and zeolites with high selectivity (SL+M > 10) for the group of n-hexane, 2-methylpentane, and 3-methylpentane (linear and monobranched isomers) versus 2,2-dimethylbutane and 2,3-dimethylbutane (dibranched isomers). This group of selective sorbents includes VICDOC (Fe2(BDP)3), a MOF with triangular pores that is known to exhibit high isomer selectivity and capacity. For three of these structures, the adsorption isotherms for a 10-component mixture of hexane and heptane isomers were calculated. Subsequent simulations of column breakthrough curves showed that the DEYVUA MOF exhibits a longer process cycle time than VICDOC MOF or MRE zeolite, which are previously reported, high-performing materials, illustrating the importance of capacity in designing MOFs for practical applications. Among the identified candidates, we synthesized and characterized a MOF in a new copper form with high predicted adsorbent capacity (qL+M > 1.2 mol/L) and moderately high selectivity (SL+M ≈ 10). Finally, we examined the role of pore shape in hexane isomer separations, especially of triangular-shaped pores. We show through the potential energy surface and three-dimensional siting analyses that linear alkanes do not populate the corners of narrow triangular channels and that structures with nontriangular pores can efficiently separate hexane isomers. Detailed thermodynamic analysis illustrates how differences in the free energy of adsorption contribute to shape-selective separation in nanoporous materials. [ABSTRACT FROM AUTHOR]

The methane deliverable capacity of adsorbent materials is a critical performance metric that will determine the viability of using adsorbed natural gas (ANG) technology in vehicular applications. ARPA-E recently set a target deliverable capacity of 315 cc(STP)/cc that a viable adsorbent material should achieve to yield a driving range competitive with incumbent fuels. However, recent computational screening of hundreds of thousands of materials suggests that the target is unattainable. In this work, we aim to determine whether the observed limits in deliverable capacity (~200 cc(STP)/cc) are fundamental limits arising from thermodynamic or material design constraints. Our efforts focus on simulating methane adsorption isotherms in a large number of systems, resulting in a broad exploration of different combinations of spatial distributions and energetics of adsorption sites. All systems were classified into five adsorption scenarios with varying degrees of realism in the manner that adsorption sites are created and endowed with energetics. The scenarios range from methane adsorption on discrete idealized lattice sites to adsorption in metal–organic frameworks with coordinatively unsaturated sites (CUS) provided by metalated catechol groups. Our findings strongly suggest that the ARPA-E target is unattainable, although not due to thermodynamic constraints but due to material design constraints. On the other hand, we also find that the currently observed deliverable capacity limits may be moderately surpassed. For instance, incorporation of CUS in IRMOF-10 is predicted to yield a 217 cc(STP)/cc deliverable capacity. The modified material has a ~0.85 void fraction and a heat of adsorption of ~15 kJ/mol. This suggests that similar, moderate improvements over existing materials could be achieved as long as CUS incorporation still maintains a relatively large void fraction. Nonetheless, we conclude that more significant improvements in deliverable capacity will require changes in the currently proposed operation conditions. [ABSTRACT FROM AUTHOR]