Thermodynamically Consistent Methodology for Estimation of Diffusivities of Mixtures of Guest Molecules in...

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: Đ ₁, Đ ₂, and Đ ₁₂. The M-S diffusivities Đ ₁ and Đ ₂ characterize interactions of guest molecules with pore walls. The exchange coefficient Đ ₁₂ 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 Đ ₁, Đ ₂, and Đ ₁₂ using input data for the constituent unary systems. The dependence of the unary diffusivities Đ ₁ and Đ ₂ 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 Đ ₁ and Đ ₂ is assumed to hold; in this case, the occupancy, θ, is calculated using the ideal adsorbed solution theory. The exchange coefficient Đ ₁₂ 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 (H₂/CO₂, Ne/CO₂, CH₄/CO₂, CO₂/N₂, H₂/CH₄, H₂/Ar, CH₄/Ar, Ne/Ar, CH₄/C₂H₆, CH₄/C₃H₈, and C₂H₆/C₃H₈) in zeolites (MFI, BEA, ISV, FAU, NaY, NaX, LTA, CHA, and DDR) and MOFs (IRMOF-1, CuBTC, and MgMOF-74).