Adsorption isotherm sensitivity to small changes in zeolite structure

Connecting theory and simulation with experiment for the study of diffusion in nanoporous solids

Nanoporous solids are ubiquitous in chemical, energy, and environmental processes, where controlled transport of molecules through the pores plays a crucial role. They are used as sorbents, chromatographic or membrane materials for separations, and as catalysts and catalyst supports. Defined as materials where confinement effects lead to substantial deviations from bulk diffusion, nanoporous materials include crystalline microporous zeotypes and metal–organic frameworks (MOFs), and a number of semi-crystalline and amorphous mesoporous solids, as well as hierarchically structured materials, containing both nanopores and wider meso- or macropores to facilitate transport over macroscopic distances. The ranges of pore sizes, shapes, and topologies spanned by these materials represent a considerable challenge for predicting molecular diffusivities, but fundamental understanding also provides an opportunity to guide the design of new nanoporous materials to increase the performance of transport limited processes. Remarkable progress in synthesis increasingly allows these designs to be put into practice. Molecular simulation techniques have been used in conjunction with experimental measurements to examine in detail the fundamental diffusion processes within nanoporous solids, to provide insight into the free energy landscape navigated by adsorbates, and to better understand nano-confinement effects. Pore network models, discrete particle models and synthesis-mimicking atomistic models allow to tackle diffusion in mesoporous and hierarchically structured porous materials, where multiscale approaches benefit from ever cheaper parallel computing and higher resolution imaging. Here, we discuss synergistic combinations of simulation and experiment to showcase theoretical progress and computational techniques that have been successful in predicting guest diffusion and providing insights. We also outline where new fundamental developments and experimental techniques are needed to enable more accurate predictions for complex systems.

Competitive Adsorption of Xylenes at Chemical Equilibrium in Zeolites

The separation of xylenes is one of the most important processes in the petrochemical industry. In this article, the competitive adsorption from a fluid-phase mixture of xylenes in zeolites is studied. Adsorption from both vapor and liquid phases is considered. Computations of adsorption of pure xylenes and a mixture of xylenes at chemical equilibrium in several zeolite types at 250 °C are performed by Monte Carlo simulations. It is observed that shape and size selectivity entropic effects are predominant for small one-dimensional systems. Entropic effects due to the efficient arrangement of xylenes become relevant for large one-dimensional systems. For zeolites with two intersecting channels, the selectivity is determined by a competition between enthalpic and entropic effects. Such effects are related to the orientation of the methyl groups of the xylenes. m-Xylene is preferentially adsorbed if xylenes fit tightly in the intersection of the channels. If the intersection is much larger than the adsorbed molecules, p-xylene is preferentially adsorbed. This study provides insight into how the zeolite topology can influence the competitive adsorption and selectivity of xylenes at reaction conditions. Different selectivities are observed when a vapor phase is adsorbed compared to the adsorption from a liquid phase. These insight have a direct impact on the design criteria for future applications of zeolites in the industry. MRE-type and AFI-type zeolites exclusively adsorb p-xylene and o-xylene from the mixture of xylenes in the liquid phase, respectively. These zeolite types show potential to be used as high-performing molecular sieves for xylene separation and catalysis.

Flexible Force Field Parameterization through Fitting on the Ab Initio-Derived Elastic Tensor

Constructing functional forms and their corresponding force field parameters for the metal–linker interface of metal–organic frameworks is challenging. We propose fitting these parameters on the elastic tensor, computed from ab initio density functional theory calculations. The advantage of this top-down approach is that it becomes evident if functional forms are missing when components of the elastic tensor are off. As a proof-of-concept, a new flexible force field for MIL-47(V) is derived. Negative thermal expansion is observed and framework flexibility has a negligible effect on adsorption and transport properties for small guest molecules. We believe that this force field parametrization approach can serve as a useful tool for developing accurate flexible force field models that capture the correct mechanical behavior of the full periodic structure.

Insight into the adsorption mechanism of benzene in HY zeolites: the effect of loading

An interesting two-stage adsorption mechanism, defined as “ideal adsorption” and “insertion adsorption”, was first proposed for the benzene/HY system by Metropolic Monte Carlo simulations at loadings below and above an “inflection point”.

Understanding and solving disorder in the substitution pattern of amino functionalized MIL-47(V)

Electronic energies and elastic constants of four amino functionalized MIL-47(V) supercells were computed using the plane wave density functional theory to determine the influence of the substituent positions on the organic linker.

Molecular screening of alcohol and polyol adsorption onto MFI-type zeolites

Configurational-bias grand canonical Monte Carlo (CB-GCMC) simulations and expanded ensemble (EE)-CB-GCMC simulations were performed to obtain adsorption isotherms of alcohols and polyols onto MFI-type zeolites from the gas phase and aqueous solution. In adsorption from both phases, Henry's constants and heats of adsorption at infinite dilution for straight-chain alcohols, diols, and triols in silicalite-1 are found to increase, and the saturation loadings decrease with increasing carbon number. Adsorption of straight-chain alcohols is more favorable than that of branched-chain alcohols. Henry's constants increase with increasing number of hydroxyl groups for gas-phase adsorption but decrease for adsorption from aqueous solution due to the strong hydrophilic solvent effect of water. The location of the hydroxyls does not affect significantly the adsorption from aqueous solution but does so in gas-phase adsorption. The saturation pressures for gas-phase adsorption decrease by orders of magnitude from the alcohols to the triols. Nonframework cations increase the adsorption of the small alcohols by an order magnitude at low concentrations (<1 mg/mL), but result in only a factor of 2 increase for larger alcohols like butanol at low concentrations (<0.03 mg/mL), and then decrease the adsorption at higher concentrations. Overall, the simulated results are in reasonable agreement with available experimental data.

Ab initio parametrized MM3 force field for the metal-organic framework MOF-5

A new valence force field has been developed and validated for a particular class of coordination polymers known as nanoporous metal-organic frameworks (MOFs), introduced recently by the group of Yaghi. The experimental, structural, and spectroscopic data in combination with density functional theory calculations on several model systems were used to parametrize the bonded terms of the force field, which explicitly treats the metal-oxygen interactions as partially covalent as well as distinguishes different types of oxygens in the framework. Both the experimental crystal structure of MOF-5 and vibrational infrared spectrum are reproduced reasonably well. The proposed force field is believed to be useful in atomistic simulations of adsorption/diffusion of guest molecules inside the flexible pores of this important class of MOF materials.

Molecular simulation study of adsorption and diffusion on silicalite for a benzene/CO2 mixture

The adsorption and diffusion of a binary mixture of supercritical CO2 and benzene on silicalite (MFI-type) have been studied through the grand canonical Monte Carlo and molecular dynamics (MD) simulations. The adsorption behavior of pure CO2 on silicalite was discussed in detail from the adsorption isotherms, adsorption sites, interaction energies, and isosteric heats of adsorption. For the mixture, the influences of temperature, pressure and composition on the adsorption isotherms have been examined. The adsorption site behavior of the mixture has been analyzed, and benzene molecules get adsorbed preferentially in the more spacious channel intersection positions. These simulation results suggest that SC-CO2 fluid can be used as an efficient desorbent of larger aromatics in the zeolite material. The diffusion characteristic for the benzene/CO2 mixture was studied on the basis of MD simulation. It was found that the large coadsorbed benzene molecule has a pronounced effect on the CO2 diffusion in the mixture, while the mobility of benzene molecules is very small due to geometrical restrictions.

Adsorption behavior of the (+/-)-Tröger's base enantiomers in the phase system of a silica-based packing coated with amylose tri(3,5-dimethyl carbamate) and 2-propanol and molecular modeling interpretation

The binary adsorption isotherms of the enantiomers of Tröger's base in the phase system made of Chiral Technologies ChiralPak AD [a silica-based packing coated with amylose tri(3,5-dimethyl carbamate)] as the chiral stationary phase (CSP) and 2-propanol as the mobile phase were measured by the perturbation method. The more retained enantiomer exhibits a S-shaped adsorption isotherm with a clear inflection point, the concentration of the less retained enantiomer having practically no competitive influence on this isotherm: In the entire range of concentrations studied, dq2/dC1 approximately 0. By contrast, the less retained enantiomer has a Langmuir adsorption isotherm when pure. At constant mobile phase concentrations, however, its equilibrium concentration in the adsorbed phase increases with increasing concentration of the more retained enantiomer and dq1/dC2 > 0. This cooperative adsorption behavior, opposed to the classical competitive behavior, is exceedingly rare but was clearly demonstrated in this case. Two adsorption isotherm equations that account for these physical observations were derived. They are based on the formation of an adsorbed multi-layer, as suggested by the isotherm data. The excellent agreement between the experimental overloaded elution profiles of binary mixtures and the profiles calculated with the equilibrium-dispersive model validates this binary isotherm model. The adsorption energies calculated by molecular mechanics (MM) and by molecular dynamics (MD) indicate that the chiral recognition arising from the different interactions between the functional groups of the CSP and the molecules of the Tröger's base enantiomers are mainly driven by their Van der Waals interactions. The MD data suggest that the interactions of the (-)-Tröger's base with the CSP are more favored by 8+/-(5) kJ/mol than those of (+)-Tröger's base. This difference seems to be a contributing factor to the increased retention of the - enantiomer on this chromatographic system. The modeling of the data also indicates that both enantiomers can form high stoichiometry complexes while binding onto the stationary phase, in agreement with the results of the equilibrium isotherm studies.

Computational elucidation of the transition state shape selectivity phenomenon

The most commonly cited example of a transition state shape selective reaction, m-xylene disproportionation in zeolites, is examined to determine if the local spatial environment of a reaction can significantly alter selectivity. In the studied reaction, ZPE-corrected rate limiting energy barriers are 136 kJ/mol for the methoxide-mediated pathway and 109 to 145 kJ/mol for the diphenylmethane-mediated pathway. Both pathways are likely to contribute to selectivity and disfavor one product isomer (1,3,5-trimethylbenzene), but relative selectivity to the other two isomers varies with pore geometry, mechanistic pathway, and inclusion of entropic effects. Most importantly, study of one pathway in three different common zeolite framework types (FAU, MFI, and MOR) allows explicit and practically oriented consideration of pore shape. Variation of the environment shape at the critical transition states is thus shown to affect the course of reaction. Barrier height shifts on the order of 10-20 kJ/mol are achievable. Observed selectivities do not agree with the transition state characteristics calculated here and, hence, are most likely due to product shape selectivity. Further examination of the pathways highlights the importance of mechanistic steps that do not result in isomer-defining bonds and leads to a more robust definition of transition state shape selectivity.