Dimensionality reduction of free energy profiles of benzene in silicalite-1: calculation of diffusion coefficients using transition state theory

This study describes the development and application of a new computational methodology for calculating the self-diffusivity of sorbate molecules strongly confined within shape-selective nanoporous materials. An umbrella sampling strategy, employing repulsive walls to confine the sorbate within specific regions of the pore space, is invoked to extract free energy profiles with respect to the sorbate degrees of freedom. Based on these profiles, it is shown how the multidimensional problem of translational diffusion of benzene in flexible silicalite-1 can be reduced first to a six-dimensional problem, then to a three-dimensional (3D) problem and finally, to a 1D problem. A 3D free energy distribution is accumulated as a function of the benzene centre of mass position and ultimately reduced to a set of 1D profiles for the benzene centre of mass along the pore axes. From these profiles, the rate constants for jumps executed by benzene between sorption sites are calculated using transition state theory; from the latter rate constants, the low-occupancy self-diffusivity is obtained using the MESoRReD method [Kolokathis PD, Theodorou DN. On solving the master equation in spatially periodic systems. J Chem Phys. 2012;137:034112]. The activation energy for diffusion and preferred orientations in the various sorption states in silicalite are in very favourable agreement with available experimental measurements.