Influence of adsorption thermodynamics on guest diffusivities in nanoporous crystalline materials

Exploring the adsorption behavior of molecular hydrogen on CHA-zeolite by comparing the performance of various force field methods

Molecular hydrogen (H2) adsorption plays a crucial role in numerous applications, including hydrogen storage and purification processes. Understanding the interaction of H2 with porous materials is essential for designing efficient adsorption systems. In this study, we investigate H2 adsorption on CHA-zeolite using a combination of density functional theory (DFT) and force field-based molecular dynamics (MD) simulations. Firstly, we employ DFT calculations to explore the energetic properties and adsorption sites of H2 on the CHA-zeolite framework. The electronic structure and binding energies of H2 in various adsorption configurations are analyzed, providing valuable insights into the nature of the adsorption process. Subsequently, force field methods are employed to perform extensive MD simulations, allowing us to study the dynamic behavior of H2 molecules adsorbed on the CHA-zeolite surface. The trajectory analysis provides information on the diffusion mechanisms and mobility of H2 within the porous structure, shedding light on the transport properties of the adsorbed gas. Furthermore, the combination of DFT and MD results enables us to validate and refine the force field parameters used in simulations, improving the accuracy of the model, and enhancing our understanding of the H2-CHA interactions. Our comprehensive investigation into molecular hydrogen adsorption on CHA-zeolite using density functional theory and molecular dynamics simulations yields valuable insights into the fundamental aspects of the adsorption process. These findings contribute to the development of advanced hydrogen storage and separation technologies, paving the way for efficient and sustainable energy applications.

Highlighting the Anti-Synergy between Adsorption and Diffusion in Cation-Exchanged Faujasite Zeolites

Using configurational-bias Monte Carlo simulations of adsorption equilibrium and molecular dynamics simulations of guest diffusivities of CO₂, CH₄, N₂, and O₂ in FAU zeolites with varying amounts of extra-framework cations (Na⁺ or Li⁺), we demonstrate that adsorption and diffusion do not, in general, proceed hand-in-hand. Stronger adsorption often implies reduced mobility. The anti-synergy between adsorption and diffusion has consequences for the design and development of pressure-swing adsorption and membrane separation technologies for CO₂ capture and N₂/O₂ separations.

Structure and dynamics of a water/methanol mixture confined in zeolitic imidazolate framework ZIF-8 from atomistic simulations

A classical atomistic simulation study is reported for the microscopic structure and dynamics of a water/methanol mixture confined in flexible nanoporous zeolitic imidazolate framework ZIF-8. Both the radial density distribution and vivid two-dimensional density profile demonstrate that methanol molecules can roughly be viewed as "embedded" between two layers of water molecules to form a "sandwich" structure. The reason for the formation of such a specific structure is explained based on the hydrogen-bonding state and the strength of various hydrogen bonds. The investigation of guest molecular diffusion shows that the self-diffusion coefficient of confined water is generally one to two orders of magnitude smaller than that of bulk water. In addition, the dependence of the self-diffusion coefficient on loading is non-monotonic: the self-diffusion coefficient firstly shows a significant increase and then decreases at higher loading. Moreover, both the structure and dynamics of the hydrogen bond (HB) network of confined water molecules are investigated in a spatially resolved manner. The results indicate that both the HB structure and dynamics of water molecules near the ZIF-8 surface deviate significantly from those of bulk water. However, while water molecules located at the pore center are relatively similar to bulk water molecules with respect to the HB structure, they exhibit strong slowdown in HB dynamics when compared with bulk water. This simulation study elucidates in detail the structural and dynamical properties of a water/methanol mixture in nanoscopic ZIF-8 confinement, which is expected to provide a deep insight into the role of porous fillers, such as ZIF-8, in improving the performance of the dehydration of alcohols via pervaporation and other related processes.

Using Molecular Simulations to Unravel the Benefits of Characterizing Mixture Permeation in Microporous Membranes in Terms of the Spreading Pressure

The separation performance of microporous crystalline materials in membrane constructs is dictated by a combination of mixture adsorption and intracrystalline diffusion characteristics; the permeation selectivity S _perm is a product of the adsorption selectivity S _ads and the diffusion selectivity, S _diff. The primary objective of this article is to gain fundamental insights into S _ads and S _diff by use of molecular simulations. We performed configurational-bias Monte Carlo (CBMC) simulations of mixture adsorption equilibrium and molecular dynamics (MD) simulations of guest self-diffusivities of a number of binary mixtures of light gaseous molecules (CO₂, CH₄, N₂, H₂, and C₂H₆) in a variety of microporous hosts of different pore dimensions and topologies. Irrespective of the bulk gas compositions and bulk gas fugacities, the adsorption selectivity, S _ads, is found to be uniquely determined by the adsorption potential, Φ, a convenient and practical proxy for the spreading pressure π that is calculable using the ideal adsorbed solution theory for mixture adsorption equilibrium. The adsorption potential Φ is also a proxy for the pore occupancy and is the thermodynamically appropriate yardstick to determine the loading and composition dependences of intracrystalline diffusivities and diffusion selectivities, S _diff. When compared at the same Φ, the component permeabilities, Π _i for CO₂, CH₄, and N₂, determinable from CBMC/MD data, are found to be independent of the partners in the various mixtures investigated and have practically the same values as the values for the corresponding unary permeabilities. In all investigated systems, the H₂ permeability in a mixture is significantly lower than the corresponding unary value. These reported results have important practical consequences in process development and are also useful for screening of materials for use as membrane devices.

Diffusion in nanopores: correlating experimental findings with “first-principles” predictions

Measurement of molecular diffusion in nanoporous host materials, which are typically inhomogeneous and anisotropic, often involves an intricate web of factors and relations to be taken into account since the associated diffusivities are a function of the diffusion path of the guest molecules during a given observation time. Depending on the observation time, therefore, the result of the experimental measurement can point to completely different conclusions about the underlying diffusion phenomena. The risk of misinterpretation of the experimental data, by correlating them with irrelevant phenomena, may be reduced if there is an option to compare the data with the results of totally independent measurements. The present communication addresses this issue with reference to the particular potentials of pulsed field gradient NMR and microimaging by infrared microscopy as techniques of microscopic diffusion measurement.

Highlighting the Influence of Thermodynamic Coupling on Kinetic Separations with Microporous Crystalline Materials

The main focus of this article is on mixture separations that are driven by differences in intracrystalline diffusivities of guest molecules in microporous crystalline adsorbent materials. Such "kinetic" separations serve to over-ride, and reverse, the selectivities dictated by mixture adsorption equilibrium. The Maxwell-Stefan formulation for the description of intracrystalline fluxes shows that the flux of each species is coupled with that of the partner species. For n-component mixtures, the coupling is quantified by a n × n dimensional matrix of thermodynamic correction factors with elements Γ _ij ; these elements can be determined from the model used to describe the mixture adsorption equilibrium. If the thermodynamic coupling effects are essentially ignored, i.e., the Γ _ij is assumed to be equal to δ _ij , the Kronecker delta, the Maxwell-Stefan formulation degenerates to yield uncoupled flux relations. The significance of thermodynamic coupling is highlighted by detailed analysis of separations of five different mixtures: N₂/CH₄, CO₂/C₂H₆, O₂/N₂, C₃H₆/C₃H₈, and hexane isomers. In all cases, the productivity of the purified raffinate, containing the tardier species, is found to be significantly larger than that anticipated if the simplification Γ _ij = δ _ij is assumed. The reason for the strong influence of Γ _ij on transient breakthroughs is traceable to the phenomenon of uphill intracrystalline diffusion of more mobile species. The major conclusion to emerge from this study is that modeling of kinetic separations needs to properly account for the thermodynamic coupling effects.

Efficient configurational-bias Monte-Carlo simulations of chain molecules with "swarms" of trial configurations

The pruned-enriched Rosenbluth method (PERM) is a popular and powerful Monte-Carlo technique for sampling flexible chain polymers of substantial length. In its original form, however, the method cannot be applied in Markov-chain Monte-Carlo schemes, which has rendered PERM unsuited for systems that consist of many chains. The current work builds on the configurational-bias Monte-Carlo (CBMC) method. The growth of a large set of trial configurations in each move is governed by simultaneous pruning and enrichment events, which tend to replace configurations with a low statistical weight by clones of stronger configurations. In simulations of dense brushes of flexible chains, a gain in efficiency of at least three orders of magnitude is observed with respect to CBMC and one order of magnitude with respect to recoil-growth approaches. Moreover, meaningful statistics can be collected from all trial configurations through the so-called "waste-recycling" Monte Carlo scheme.

Molecular Modelling for Reactor Design

Chemical reactor modelling based on insights and data on a molecular level has become reality over the last few years. Multiscale models describing elementary reaction steps and full microkinetic schemes, pore structures, multicomponent adsorption and diffusion inside pores, and entire reactors have been presented. Quantum mechanical (QM) approaches, molecular simulations (Monte Carlo and molecular dynamics), and continuum equations have been employed for this purpose. Some recent developments in these approaches are presented, in particular time-dependent QM methods, calculation of van der Waals forces, new approaches for force field generation, automatic setup of reaction schemes, and pore modelling. Multiscale simulations are discussed. Applications of these approaches to heterogeneous catalysis are demonstrated for examples that have found growing interest over the last few years, such as metal-support interactions, influence of pore geometry on reactions, noncovalent bonding, reaction dynamics, dynamic changes in catalyst nanoparticle structure, electrocatalysis, solvent effects in catalysis, and multiscale modelling.

The effect of gme topology on multicomponent adsorption in zeolitic imidazolate frameworks

We employ a simulation approach to study the adsorption of single, binary and ternary mixtures on eight gme Zeolitic Imidazolate Frameworks (ZIFs) at 298 K. Four adsorbate fluids were considered; carbon dioxide, methane, nitrogen and water. We compute the high pressure adsorption density profiles inside the micropore channels of each crystal. The profiles are compared directly for the different structures and adsorbate components and used to highlight the influence of the imidazolate ligands on pure and competitive adsorption. ZIFs with long ligands reveal an additional, accessible to the fluid space detected for the first time. This is a wedged volume on one direction of the pore walls, shaped thus because the long ligands tilt in order to be connected. We estimate the pressure required for water to become the dominating competing adsorbate within different crystal cavities. The simulated data for CO₂ adsorption in ZIF69 strongly correlate with a set of Raman spectroscopy intensity values that correspond to the same adsorbate-adsorbent system.

Connection Between Thermodynamics and Dynamics of Simple Fluids in Pores: Impact of Fluid-Fluid Interaction Range and Fluid-Solid Interaction Strength

Using molecular simulations, we investigate how the range of fluid-fluid (adsorbate-adsorbate) interactions and the strength of fluid-solid (adsorbate-adsorbent) interactions impact the strong connection between distinct adsorptive regimes and distinct self-diffusivity regimes reported in [Krekelberg, W. P.; Siderius, D. W.; Shen, V. K.; Truskett, T. M.; Errington, J. R. Langmuir2013, 29, 14527-14535]. Although increasing the fluid-fluid interaction range changes both the thermodynamics and the dynamic properties of adsorbed fluids, the previously reported connection between adsorptive filling regimes and self-diffusivity regimes remains. Increasing the fluid-fluid interaction range leads to enhanced layering and decreased self-diffusivity in the multilayer-formation regime but has little effect on the properties within film-formation and pore-filling regimes. We also find that weakly attractive adsorbents, which do not display distinct multilayer formation, are hard-sphere-like at super- and subcritical temperatures. In this case, the self-diffusivity of the confined and bulk fluid has a nearly identical scaling-relationship with effective density.

Insights into the Activity and Deactivation of the Methanol-to-Olefins Process over Different Small-Pore Zeolites As Studied with Operando UV-vis Spectroscopy

The nature and evolution of the hydrocarbon pool (HP) species during the Methanol-to-Olefins (MTO) process for three small-pore zeolite catalysts, with a different framework consisting of large cages interconnected by small eight-ring windows (CHA, DDR, and LEV) was studied at reaction temperatures between 350 and 450 °C using a combination of operando UV-vis spectroscopy and online gas chromatography. It was found that small differences in cage size, shape, and pore structure of the zeolite frameworks result in the generation of different hydrocarbon pool species. More specifically, it was found that the large cage of CHA results in the formation of a wide variety of hydrocarbon pool species, mostly alkylated benzenes and naphthalenes. In the DDR cage, 1-methylnaphthalene is preferentially formed, while the small LEV cage generally contains fewer hydrocarbon pool species. The nature and evolution of these hydrocarbon pool species was linked with the stage of the reaction using a multivariate analysis of the operando UV-vis spectra. In the 3-D pore network of CHA, the reaction temperature has only a minor effect on the performance of the MTO catalyst. However, for the 2-D pore networks of DDR and LEV, an increase in the applied reaction temperature resulted in a dramatic increase in catalytic activity. For all zeolites in this study, the role of the hydrocarbon species changes with reaction temperature. This effect is most clear in DDR, in which diamantane and 1-methylnaphthalene are deactivating species at a reaction temperature of 350 °C, whereas at higher temperatures diamantane formation is not observed and 1-methylnaphthalene is an active species. This results in a different amount and nature of coke species in the deactivated catalyst, depending on zeolite framework and reaction temperature.

Highlighting diffusional coupling effects in zeolite catalyzed reactions by combining the Maxwell–Stefan and Langmuir–Hinshelwood formulations

The combined phenomena of intra-crystalline adsorption, diffusion and reversible chemical reactions inside microporous crystalline zeolite catalyst particles are described by combining the Langmuir–Hinshelwood kinetics with the Maxwell–Stefan (M–S) diffusion formulation.

Interplay between hydrophilicity and surface barriers on water transport in zeolite membranes

A comprehensive understanding of molecular transport within nanoporous materials remains elusive in a broad variety of engineering and biomedical applications. Here, experiments and atomistic simulations are synergically used to elucidate the non-trivial interplay between nanopore hydrophilicity and surface barriers on the overall water transport through zeolite crystals. At these nanometre-length scales, these results highlight the dominating effect of surface imperfections with reduced permeability on the overall water transport. A simple diffusion resistance model is shown to be sufficient to capture the effects of both intracrystalline and surface diffusion resistances, thus properly linking simulation to experimental evidence. This work suggests that future experimental work should focus on eliminating/overcoming these surface imperfections, which promise an order of magnitude improvement in permeability.

Zeolite crystal with porous structure is predicted to be a good membrane material for water purification, but experiments show water uptake orders of magnitude smaller than the theory. Here, Fasano et al. attribute this disagreement to the additional diffusion resistance induced by surface defects.

Separation of benzene from mixtures with water, methanol, ethanol, and acetone: highlighting hydrogen bonding and molecular clustering influences in CuBTC

Configurational-bias Monte Carlo (CBMC) simulations are used to establish the potential of CuBTC for separation of water/benzene, methanol/benzene, ethanol/benzene, and acetone/benzene mixtures. For operations under pore saturation conditions, the separations are in favor of molecules that partner benzene; this is due to molecular packing effects that disfavor benzene. CBMC simulations for adsorption of quaternary water/methanol/ethanol/benzene mixtures show that water can be selectively adsorbed at pore saturation, making CuBTC effective in drying applications. Ideal Adsorbed Solution Theory (IAST) calculations anticipate the right hierarchy of component loadings but the quantitative agreement with CBMC mixture simulations is poor for all investigated mixtures. The failure of the IAST to provide reasonable quantitative predictions of mixture adsorption is attributable to molecular clustering effects that are induced by hydrogen bonding between water-water, methanol-methanol, and ethanol-ethanol molecule pairs. There is, however, no detectable hydrogen bonding between benzene and partner molecules in the investigated mixtures. As a consequence of molecular clustering, the activity coefficients of benzene in the mixtures is lowered below unity by one to three orders of magnitude at pore saturation; such drastic reductions cannot be adequately captured by the Wilson model, that does not explicitly account for molecular clustering. Molecular clustering effects are also shown to influence the loading dependence of the diffusivities of guest molecules.

Zeolite membranes – a review and comparison with MOFs

The latest developments in zeolite and MOF membranes are reviewed, with an emphasis on synthesis techniques. Industrial applications, hydrothermal stability, polymer-supported and mixed matrix membranes are some of the aspects discussed.

Uphill diffusion in multicomponent mixtures

Diffusional coupling effects may cause uphill transport of a species against its concentration gradient. Transient overshoots and serpentine equilibration trajectories in composition space are fingerprints of uphill diffusion. Serpentine diffusion trajectories may enter meta-stable zones, leading to spontaneous emulsification and the Ouzo effect.

Diffusion of interacting particles in discrete geometries: Equilibrium and dynamical properties

We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.110601]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self- and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self- and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Water transport control in carbon nanotube arrays

Based on a recent scaling law of the water mobility under nanoconfined conditions, we envision novel strategies for precise modulation of water diffusion within membranes made of carbon nanotube arrays (CNAs). In a first approach, the water diffusion coefficient D may be tuned by finely controlling the size distribution of the pore size. In the second approach, D can be varied at will by means of externally induced electrostatic fields. Starting from the latter strategy, switchable molecular sieves are proposed, where membranes are properly designed with sieving and permeation features that can be dynamically activated/deactivated. Areas where a precise control of water transport properties is beneficial range from energy and environmental engineering up to nanomedicine.

Microimaging of transient guest profiles to monitor mass transfer in nanoporous materials

The intense interactions of guest molecules with the pore walls of nanoporous materials is the subject of continued fundamental research. Stimulated by their thermal energy, the guest molecules in these materials are subject to a continuous, irregular motion, referred to as diffusion. Diffusion, which is omnipresent in nature, influences the efficacy of nanoporous materials in reaction and separation processes. The recently introduced techniques of microimaging by interference and infrared microscopy provide us with a wealth of information on diffusion, hitherto inaccessible from commonly used techniques. Examples include the determination of surface barriers and the sticking coefficient's analogue, namely the probability that, on colliding with the particle surface, a molecule may continue its diffusion path into the interior. Microimaging is further seen to open new vistas in multicomponent guest diffusion (including the detection of a reversal in the preferred diffusion pathways), in guest-induced phase transitions in nanoporous materials and in matching the results of diffusion studies under equilibrium and non-equilibrium conditions.

Fabrication of cuprous nanoparticles in MIL-101: an efficient adsorbent for the separation of olefin–paraffin mixtures

Cu⁺-doped MIL-101 was synthesized, and its application in the separation of olefin–paraffin was investigated for the first time.

Connection between thermodynamics and dynamics of simple fluids in highly attractive pores

Using molecular simulations, we investigate the structural and diffusive dynamics properties of a model fluid in highly absorptive cylindrical pores. At subcritical temperatures, self-diffusivity displays three distinct regimes as a function of average pore density ρ: (1) a decrease in self-diffusivity with increasing ρ at low ρ, (2) constant self-diffusivity with respect to varying ρ at moderate density, and (3) a decrease in self-diffusivity with increasing ρ at high density. These regimes are closely linked to the thermodynamic properties of the fluid in the pore, specifically, the adsorption isotherm, isosteric heat of adsorption, and the density profile. We show that these three diffusivity regimes qualitatively correspond to three distinct adsorption regimes: monolayer formation, multilayer adsorption, and pore filling, respectively. In addition, we find that the self-diffusivity is a universal function of the local film density in the monolayer formation regime at subcritical temperatures. The results of this work suggest a potential means to estimate the self-diffusivity over a broad pressure range using a limited number of experiments.