Constrained reaction coordinate (Bluemoon) dynamics were used to characterize the free energy profile of benzene in silicalite-l at 300K along the mean reaction path for diffusion. The reaction path was found empirically by fitting a parametric curve through the mean positions of the benzene centre of mass. Both rigid and flexible zeolite lattices were investigated. In both cases the primary adsorption site was located at the intersection of the straight and sinusoidal channels. Two other distinct adsorption sites exist in the straight channel and four in the sinusoidal channel. Lattice flexibility was found to have a very strong influence on the relative free energies of the adsorption sites and on the free energies of the transition states connecting them. The free energies, combined with estimates of the transmission coefficient, were used to obtain rate constants for diffusion between the main adsorption sites in the flexible lattice. Subsequent diffusive Monte Carlo simulations, using these parameters, provided the self-diffusion coefficient, and its components, at 300K. A diffusivity of 3.36 x 10-14m2/s was obtained; in excellent accord with the best experimental value of 2.2 x 10-14m2/s. Diffusion was fastest in the Y direction with Dyy/Dxx = 1.49 and Dyy/Dzz = 5.65 being deduced.

Bluemoon Simulations of Benzene in Silicalite-l Prediction of Free Energies and Diffusion Coefficients. T.R.Forester, W.Smith: Journal of the Chemical Society - Faraday Transactions, 1997, 93[17], 3249-57