Molecular Dynamics Studies of Diffusion in Model Cylindrical Pores at Very Low Densities

Abstract

Molecular dynamics simulation has been used to study diffusion of methane at ambient temperature in cylindrical pores at very low densities. The cylinders were modelled as a continuum solid which interacts with the methane in the radial direction only. At the lowest densities, the VACF method does not yield reliable values of the self diffusion coefficient, D_s , but a suitable choice of time step and run length enables values of D_s to be found from MSD plots that are below the classical Knudsen diffusion coefficients. When density is increased, D_s passes through a maximum although the adsorption isotherm remains inside the Henry law region. Maxima are found for two cylinder radii and for two adsorbent field strengths. The existence of a maximum is attributed to transient intermolecular interactions. Analysis of a molecular trajectory demonstrates that long diffusion paths can be triggered by the rare event of an intermolecular encounter which forces a molecule into the repulsive part of the wall potential. At sufficiently high density, subsequent collisions quench the tendency towards long paths, and D_s decreases again. The issue of simulation artefact as a source of these observations is discussed.