Molecular dynamics simulations involving a flexible framework were used to simulate the diffusion of o- and p-xylene in purely siliceous CIT-1. The simulations were performed at 500K, investigating two loadings (corresponding to 1 or 0.25mol/unit-cell) for the ortho isomer and one loading (0.25mol/unit-cell) for the para isomer. In the former system, the diffusion coefficient decreased from 7.79 x 10-6cm2/s for the lower, to 3.56 x 10-6cm2/s for the higher loading (table 15). The diffusion coefficient for the para isomer (25.18 x 10-6cm2/s) was substantially greater. Graphical analyses revealed a jump diffusion mechanism which, in the case of the ortho isomer, took place in the 12 MR channels of the structure while, for the para isomer, incursions into the 10 MR channel were observed. The results of the molecular dynamics simulations were complemented by calculations of the activation energies for the diffusion of the two isomers in the two channels.
Selective Diffusion of C8 Aromatics in a 10 and 12 MR Zeolite. A Molecular Dynamics Study. G.Sastre, N.Raj, C.R.A.Catlow, R.Roque-Malherbe, A.Corma: Journal of Physical Chemistry B, 1998, 102[17], 3198–209
Table 14
Counter-Diffusion of 50:50 Xylene Mixtures in H-CIT-1
Diffusant | Temperature(K) | D(m2/s) |
p-xylene | 350 | 3.0 x 10-14 |
p-xylene | 375 | 3.8 x 10-13 |
p-xylene | 400 | 3.3 x 10-13 |
p-xylene | 425 | 1.2 x 10-12 |
o-xylene | 400 | 2.0 x 10-14 |