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 existed in the straight channel, and four in the sinusoidal channel. Lattice flexibility was found to have a very marked 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 value of 3.36 x 10-14m2/s was obtained. This was 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 observed.

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