A molecular dynamics investigation of the unusual concentration dependencies of Fick diffusivities in silica mesopores

Molecular Dynamics (MD) simulations were carried out to determine the self-diffusivitiy, D(i,self), the Maxwell-Stefan diffusivity, D(i), and the Fick diffusivity, D(i), for methane (Cl), ethane (C2), propane (C3), n-butane (nC4), n-pentane (nC5), n-hexane (nC6), n-heptane (nC7), and cyclohexane (cC6) in cylindrical silica mesopores for a range of pore concentrations. The MD simulations show that zero-loading diffusivity D(i)(0) is consistently lower, by up to a factor of 20, than the values anticipated by the classical Knudsen formula. The concentration dependence of the Fick diffusivity, D(i) is found to be unusually complex, and displays a strong minimum in some cases; this characteristic can be traced to molecular clustering. (C) 2010 Elsevier Inc. All rights reserved.