Investigation of the Air Separation Properties of Zeolites Types A, X and Y by Monte Carlo Simulations

Capture, Storage, and Release of Oxygen by Metal–Organic Frameworks (MOFs)

Oxygen is a critical gas for medical and industrial settings. Much of today's global oxygen supply is via inefficient technologies such as cryogenic distillation, membranes or zeolites. Metal–organic frameworks (MOFs) promise a superior alternative for oxygen separation, as their fundamental chemistry can in principle be tailored for reversible and selective oxygen capture. We evaluate the characteristics for reversible and selective uptake of oxygen by MOFs, focussing on redox‐active sites. Key characteristics for separation can also be seen in MOFs for oxygen storage roles. Engineering solutions to release adsorbed oxygen from the MOFs are discussed including Temperature Swing Adsorption (TSA), Pressure Swing Adsorption (PSA) and the highly efficient Magnetic Induction Swing Adsorption (MISA). We conclude with the applications and outlooks for oxygen capture, storage and release, and the likely impacts the next generation of MOFs will have on industry and the broader community.

Cytochrome P450-The Wonderful Nanomachine Revealed through Dynamic Simulations of the Catalytic Cycle

This Account addresses the catalytic cycle of the enzyme cytochrome P450 (CYP450) as a prototypical biological machine with automatic features. CYP450 is a nanomachine that uses dioxygen and two reducing and two proton equivalents to oxidize a plethora of molecules (so-called substrates) as a means of supplying bio-organisms with essential molecules (e.g., brain neurotransmitters, sex hormones, etc.) and protecting biosystems against poisoning. An enticing property of CYP450s is that entrance of an oxidizable substrate into the active site initiates a series of events that constitute the catalytic cycle, which functions "automatically" in a regulated sequence of events culminating in the production of the oxidized substrates (e.g., hydroxylated, epoxidized, etc.), oftentimes with remarkable stereo- and regioselectivities. It is timely to demonstrate how theory uses molecular dynamics (MD) simulations and quantum-mechanical/molecular-mechanical (QM/MM) calculations to complement experiments and elucidate the choreography by which the protein regulates the catalytic cycle. CYP450 is a heme enzyme that contains a ferric ion (Fe^III) coordinated by a porphyrin ligand, a water molecule, and a cysteinate ligand that is provided by a strategic residue of the encapsulating protein. While many of the individual steps are sufficiently well-understood, we shall provide here an overview of the factors that cause all of the steps to be sequentially coordinated. To this end, we use examples from three different CYP450 enzymes: the bacterial ones CYP450_BM3 and CYP450_CAM and the mammalian enzyme CYP450_3A4. The treatment is limited to the catalytic cycle, as aspects of two-state reactivity were reviewed previously (e.g., Shaik , S. ; et al. Chem. Rev. 2005 , 105 , 2279 ). What are the principles that govern the seeming automatic feature? For example, how do substrate entrance and binding gate the enzyme? How does the reductase attachment to the enzyme affect the next steps? What triggers the attachment of the reductase? How does the electron transfer (ET) that converts Fe^III to Fe^II occur? Is the ET coordinated with the entrance of O₂ into the active site? What is the mechanism of the latter step? Since the entrance of the substrate expels the water molecules from the active site, how do water molecules re-enter to form a proton channel, which is necessary for creating the ultimate oxidant Compound I? How do mutations that disrupt the water channel nevertheless create a competent oxidant? By what means does the enzyme produce regio- and stereoselective oxidation products? What triggers the departure of the oxidized product, and how does the exit occur in a manner that generates the resting state ready for the next cycle? This Account shows that the entrance of the substrate triggers all of the ensuing events.

Mechanism of Kinetically Controlled Capillary Condensation in Nanopores: A Combined Experimental and Monte Carlo Approach

We find the rule of capillary condensation from the metastable state in nanoscale pores based on the transition state theory. The conventional thermodynamic theories cannot achieve it because the metastable capillary condensation inherently includes an activated process. We thus compute argon adsorption isotherms on cylindrical pore models and atomistic silica pore models mimicking the MCM-41 materials by the grand canonical Monte Carlo and the gauge cell Monte Carlo methods and evaluate the rate constant for the capillary condensation by the transition state theory. The results reveal that the rate drastically increases with a small increase in the chemical potential of the system, and the metastable capillary condensation occurs for any mesopores when the rate constant reaches a universal critical value. Furthermore, a careful comparison between experimental adsorption isotherms and the simulated ones on the atomistic silica pore models reveals that the rate constant of the real system also has a universal value. With this finding, we can successfully estimate the experimental capillary condensation pressure over a wide range of temperatures and pore sizes by simply applying the critical rate constant.

Critical energy barrier for capillary condensation in mesopores: Hysteresis and reversibility

Capillary condensation in the regime of developing hysteresis occurs at a vapor pressure, Pcond, that is less than that of the vapor-like spinodal. This is because the energy barrier for the vapor-liquid transition from a metastable state at Pcond becomes equal to the energy fluctuation of the system; however, a detailed mechanism of the spontaneous transition has not been acquired even through extensive experimental and simulation studies. We therefore construct accurate atomistic silica mesopore models for MCM-41 and perform molecular simulations (gauge cell Monte Carlo and grand canonical Monte Carlo) for argon adsorption on the models at subcritical temperatures. A careful comparison between the simulation and experiment reveals that the energy barrier for the capillary condensation has a critical dimensionless value, Wc (*) = 0.175, which corresponds to the thermal fluctuation of the system and depends neither on the mesopore size nor on the temperature. We show that the critical energy barrier Wc (*) controls the capillary condensation pressure Pcond and also determines a boundary between the reversible condensation/evaporation regime and the developing hysteresis regime.

Molecular Dynamics Study of Hydrogen on Alkali-Earth Metal Cations Exchanged X Zeolites

The self-diffusion of hydrogen in Ca2+-, Mg2+- and Ba2+-exchanged X zeolites (Mg46X, Ca46X, and Ba46X) has been studied by molecular dynamics (MD) simulations for various temperatures and loadings. The results indicate that in the temperature range of 77–298 K and the loading range of 1–80 molecules/cell, the self-diffusion coefficients are found to range from1.2×10-9 m2·s−1to2.3×10-7 m2·s−1which are in good agreement with the experimental values from the quasielastic neutron scattering (QENS) and pulse field gradients nuclear magnetic resonance (PFG NMR) measurements. The self-diffusion coefficients decrease with loading due to packing of sorbate-sorbate molecules which causes frequent collusion among hydrogen molecules in pores and increases with increasing temperature because increasing the kinetic energy of the gas molecules enlarges the mean free path of gas molecule. The mechanism of diffusion of hydrogen molecules in these zeolites is transition diffusion. Knudsen diffusion occurs at low loading and the molecular bulk diffusion occurs at higher loading. For given temperature and loading, the self-diffusion coefficients decrease in the orderBa46X<Mg46X<Ca46X, due to the different sizes and locations of the divalent cations. Moreover, the effect of concentration of molecular hydrogen on self-diffusion coefficient also is analyzed using radial distribution function (RDF).

Molecular-Dynamics Simulation of Self-Diffusion of Molecular Hydrogen in X-Type Zeolite

The self-diffusion of hydrogen in NaX zeolite has been studied by molecular-dynamics simulations for various temperatures and pressures. The results indicate that in the temperature range of 77–293 K and the pressure range of 10–2700 kPa, the self-diffusion coefficients are found to range from 1.61 × 10−9 m2·s−1to 3.66 × 10−8 m2·s−1which are in good agreement with the experimental values from the quasielastic neutron scattering (QENS) and pulse field gradients nuclear magnetic resonance (PFG NMR) measurements. The self-diffusion coefficients decrease with increasing pressure due to packing of sorbate-sorbate molecules which causes frequent collusion among hydrogen molecules in pores and increase with increasing temperature because increasing the kinetic energy of the gas molecules enlarges the mean free path of gas molecule. The activated energy for hydrogen diffusion determined from the simulation is pressure-dependent.

Correlation of sorption behavior of nitrogen, oxygen, and argon with Ca2+ locations in zeolite a: a grand canonical Monte Carlo simulation study

The adsorption isotherms for nitrogen, oxygen, and argon in various NaCaA zeolite samples were calculated theoretically using the grand canonical Monte Carlo simulation method. The adsorption capacity, selectivity, and heat of adsorption of nitrogen increase with an increasing number of calcium cations in zeolite A. The heat of adsorption of nitrogen showed a sudden increase when the calcium ion exchange to zeolite A was around 60%. Adsorption isotherms, determined experimentally by the volumetric adsorption method, support theoretically predicted isotherms. These observations have been explained in terms of the interaction of the nitrogen molecule with Ca2+ ions and their locations in zeolite A.

Influence of extra-framework cations on the adsorption properties of X-faujasite systems: microcalorimetry and molecular simulations

Isotherms and differential enthalpies of adsorption are obtained for nitrogen at ambient temperature on monovalent (Li(+), Na(+), K(+)) and divalent (Ca(2+), Ba(2+), Sr(2+), Mn(2+)) substituted X-faujasite systems by microcalorimetry measurements. These experimental data are compared with those obtained by combining grand canonical Monte Carlo simulations and newly derived force fields for describing the interactions between the extra-framework cations and the adsorbates obtained from a simple model based only on the intrinsic properties of the cations. It is the first time that such good qualitative agreement is reported between experiment and simulation for a series of both monovalent and divalent cations.

Molecular dynamics-based approach to study the anisotropic self-diffusion of molecules in porous materials with multiple cage types: application to H2 in losod

The anisotropic self-diffusion of molecular hydrogen in the multiple cage clathrasil losod (LOS) is modeled by means of molecular dynamics (MD) simulations of up to 1 micros for the temperature range 900-1200 K while treating the framework as fully flexible. The LOS diffusion tensor is calculated employing an analytical method based on hopping rates. The diffusion in the c-direction of the unit cell is found to be approximately two times more rapid than in the a- and the b-directions, a characteristic of importance for the application of LOS as a membrane. The overall diffusion is based on five different hop types for which the individual hopping rates and diffusion barriers are calculated separately. We show explicitly that the shape and volume of the cages have a significant effect on the hopping rates and further that even small deformations of the circular Si6O6 apertures have a large influence on the energetic barrier for hydrogen diffusion. Compared to the single cage clathrasils dodecasil 3C (MTN) and sodalite (SOD), LOS has a lower diffusion rate. However, from a technical point of view this rate (at 573 K) is still fast enough for LOS to be interesting as a size-selective membrane or as a hydrogen-adsorption medium.

Self-diffusion of molecular hydrogen in clathrasils compared: Dodecasil 3C versus sodalite

The self-diffusion coefficient of molecular hydrogen through the all-silica microporous dodecasil 3C structure is calculated by means of molecular-dynamics (MD) calculations, allowing for full framework flexibility, in order to assess the material's feasibility as a hydrogen storage medium. The hydrogen uptake rate into dodecasil 3C is compared to that previously calculated for sodalite and it is found that the latter performs significantly better. The reason for this variation in performance is found to lie in intrinsic topological differences between each framework type. This is explicitly demonstrated by means of a simplified version of transition state theory helping to succinctly rationalize the MD data.

Effect of cation distribution on self-diffusion of molecular hydrogen in Na3Al3Si3O12 sodalite: A molecular dynamics study

The diffusion of hydrogen in sodium aluminum sodalite (NaAlSi-SOD) is modeled using classical molecular dynamics, allowing for full flexibility of the host framework, in the temperature range 800-1200 K. From these simulations, the self-diffusion coefficient is determined as a function of temperature and the hydrogen uptake at low equilibrium hydrogen concentration is estimated at 573 K. The influence of the cation distribution over the framework on the hydrogen self-diffusion is investigated by comparing results employing a low energy fully ordered cation distribution with those obtained using a less ordered distribution. The cation distribution is found to have a surprisingly large influence on the diffusion, which appears to be due to the difference in framework flexibility for different cation distributions, the occurrence of correlated hopping in case of the ordered distribution, and the different nature of the diffusion processes in both systems. Compared to our previously reported calculations on all silica sodalite (all-Si-SOD), the hydrogen diffusion coefficient of sodium aluminum sodalite is higher in the case of the ordered distribution and lower in case of the disordered distribution. The hydrogen uptake rates of all-Si-SOD and NaSiAl-SOD are comparable at high temperatures (approximately 1000 K) and lower for all-Si-SOD at lower temperatures (approximately 400 K).

Computer simulations of adsorption and diffusion for binary mixtures of methane and hydrogen in titanosilicates

Equilibrium molecular dynamics (MD) simulations of equimolar mixtures of hydrogen and methane were performed in three different titanosilicates: naturally occurring zorite and two synthetic titanosilicates, ETS-4 and ETS-10. In addition, single-component MD simulations and adsorption isotherms generated using grand canonical Monte Carlo simulations were performed to support the mixture simulations. The goal of this study was to determine the best membrane material to carry out hydrogen/methane separations. ETS-10 has a three-dimensional pore network. ETS-4 and zorite have two-dimensional pore networks. The simulations carried out in this study show that the increased porosity of ETS-10 results in self-diffusion coefficients for both hydrogen and methane that are higher in ETS-10 than in either ETS-4 or zorite. Methane only showed appreciable displacement in ETS-10. The ability of the methane molecules to move in all three directions in ETS-10 was demonstrated by the high degree of isotropy shown in the values of the x, y, and z components of the self-diffusion coefficient for methane in ETS-10. From our simulations we conclude that ETS-10 would be better suited for fast industrial separations of hydrogen and methane. However, the separation would not result in a pure hydrogen stream. In contrast, ETS-4 and zorite would act as true molecular sieves for separations of hydrogen and methane, as the methane would not move through membranes made of these materials. This was indicated by the near-zero self-diffusion coefficient of methane in ETS-4 and zorite.

Molecular-dynamics analysis of the diffusion of molecular hydrogen in all-silica sodalite

In order to investigate the technical feasibility of crystalline porous silicates as hydrogen storage materials, the self-diffusion of molecular hydrogen in all-silica sodalite is modeled using large-scale classical molecular-dynamics simulations employing full lattice flexibility. In the temperature range of 700-1200 K, the diffusion coefficient is found to range from 1.610(-10) to 1.810(-9) m(2)/s. The energy barrier for hydrogen diffusion is determined from the simulations allowing the application of transition state theory, which, together with the finding that the pre-exponential factor in the Arrhenius-type equation for the hopping rate is temperature-independent, enables extrapolation of our results to lower temperatures. Estimates based on mass penetration theory calculations indicate a promising hydrogen uptake rate at 573 K.

The 'zero charge' partitioning behaviour of noble gases during mantle melting

Noble-gas geochemistry is an important tool for understanding planetary processes from accretion to mantle dynamics and atmospheric formation. Central to much of the modelling of such processes is the crystal-melt partitioning of noble gases during mantle melting, magma ascent and near-surface degassing. Geochemists have traditionally considered the 'inert' noble gases to be extremely incompatible elements, with almost 100 per cent extraction efficiency from the solid phase during melting processes. Previously published experimental data on partitioning between crystalline silicates and melts has, however, suggested that noble gases approach compatible behaviour, and a significant proportion should therefore remain in the mantle during melt extraction. Here we present experimental data to show that noble gases are more incompatible than previously demonstrated, but not necessarily to the extent assumed or required by geochemical models. Independent atomistic computer simulations indicate that noble gases can be considered as species of 'zero charge' incorporated at crystal lattice sites. Together with the lattice strain model, this provides a theoretical framework with which to model noble-gas geochemistry as a function of residual mantle mineralogy.