Diffusion of alkane mixtures in zeolites: validating the Maxwell-Stefan formulation using MD simulations

Abstract

Molecular dynamics (MD) simulations have been carried out for pure components, binary, ternary, and quaternary mixtures containing methane, ethane, propane, and n-butane in FAU zeolite at 300 Κ for a range of molecular loadings ϑ, approaching saturation limits. The n-dimensional matrix of Maxwell-Stefan (M-S) diffusivities [Δ], defined by (N) = -ρ[Δ][Γ](∇ ^ϑ), was determined along with the self-diffusivities, D_i,self. Additionally, configurational-bias Monte Carlo (CBMC) simulations were carried out to obtain the pure component sorption isotherms and the saturation capacities ϑ_i,sat. From the information on Δ_ij, D_i,self, and ϑ_i,sat, the various M-S diffusivities were determined:(1) component Ð _i, reflecting the interactions of the species i with the zeolite, self-exchange Ð _ii, and (2) binary exchange Ð _ij. The obtained data underline the major advantage of the M-S formulation that at a given occupancy, θ = Σⁿ_i=1ϑ _i/ϑ _i,sat within the zeolite, the Ð_i has nearly the same value for species i whether this species is present on its own or in a mixture with other species. The same advantage holds, too, for the self-exchange Ð _ii; the value at a given occupancy, ϑ, is the same whether determined from pure component, binary, or ternary mixture data. For all binary and ternary mixtures studied, it was verified that the binary exchange coefficient Ð_ij can be interpolated from the corresponding values of the self-exchange parameters Ð_ii and Ð_jj using a generalization of the interpolation formula developed earlier (Skoulidas et al., Langmuir, 2003, 19, 7977). We also demonstrate that if the occupancy dependence of the pure component parameters Ð_i and Ð_ii are modeled properly, this information is sufficient to provide very good estimates of the matrix [Δ] for mixtures with 2, 3, or 4 components over the entire range of loadings. Simulations of mixture diffusion of alkanes in MFI and LTA confirm that the above-mentioned advantages of the M-S formulation also hold for these zeolite topologies.