Thermodynamically Consistent Methodology for Estimation of Diffusivities of Mixtures of Guest Molecules in Microporous Materials

Abstract

The Maxwell-Stefan (M-S) formulation, that is grounded in the theory of irreversible thermodynamics, is widely used for describing mixture diffusion in microporous crystalline materials such as zeolites and metal-organic frameworks (MOFs). Binary mixture diffusion is characterized by a set of three M-S diffusivities: D1, D2, and D12. The M-S diffusivities D1 and D2 characterize interactions of guest molecules with pore walls. The exchange coefficient D12 quantifies correlation effects that result in slowing-down of the more mobile species due to correlated molecular jumps with tardier partners. The primary objective of this article is to develop a methodology for estimating D1, D2, and D12 using input data for the constituent unary systems. The dependence of the unary diffusivities D1 and D2 on the pore occupancy, θ, is quantified using the quasi-chemical theory that accounts for repulsive, or attractive, forces experienced by a guest molecule with the nearest neighbors. For binary mixtures, the same occupancy dependence of D1 and D2 is assumed to hold; in this case, the occupancy, θ, is calculated using the ideal adsorbed solution theory. The exchange coefficient D12 is estimated from the data on unary self-diffusivities. The developed estimation methodology is validated using a large data set of M-S diffusivities determined from molecular dynamics simulations for a wide variety of binary mixtures (H2/CO2, Ne/CO2, CH4/CO2, CO2/N2, H2/CH4, H2/Ar, CH4/Ar, Ne/Ar, CH4/C2H6, CH4/C3H8, and C2H6/C3H8) in zeolites (MFI, BEA, ISV, FAU, NaY, NaX, LTA, CHA, and DDR) and MOFs (IRMOF-1, CuBTC, and MgMOF-74).