A variety of force field based simulations have been used to study the location and diffusion of benzene adsorbed in a model zeolite HY (Si/Al ) 2.43), namely: molecular docking; equilibrium and nonequilibrium molecular dynamics; and Monte Carlo umbrella sampling. Multiple adsorption sites are found, with benzene facially coordinated to one or two H(1) or H(2) protons in the supercage. Some slight adsorption onto the 12-membered ring windows is also observed, in accordance with infrared measurements. The minimum energy path at low temperature proceeds via a creeping of the molecule along the zeolite wall between stable sites, with an activation energy varying between 10 and 20 kJ mol -1 . This type of creeping motion is observed both for intracage and intercage diffusion. Cartwheel jumps between sites are seen to proceed with higher activation energies of approximatively 30 kJ mol -1 . Multiple paths from site to site open as the temperature increases. This results in a strong temperature dependence of the potential of mean force in the zeolite cage, as calculated by umbrella sampling. Nonequilibrium molecular dynamics simulations initialized at the transition state between two states show that the molecules do not relax in a single final state but in a multiplicity of states; only cage-to-cage jumps keep a sense, as a majority of molecules relax in the final cage. Due to the multiplicity of possible cage-to-cage paths, the temperature dependence of the cage-to-cage rate constants is deeply non-Arrhenius.

No abstract

A scheme to obtain approximate analytical functions for the atomic distributed multipole moments of the crystallographically different atoms within aluminosilicate and aluminophosphate sieves is discussed. Respective atomic multipole moments are derived within the CRYSTAL95 ab initio periodic Hartree-Fock code with different basis sets, from minimal STO-3G to 6-21G * . In order to illustrate the possible applications, distributed analyses are carried out for various structural models from all-siliceous zeolites and aluminophosphates with ratio Al/P = 1 to hydrogen forms of aluminosilicates. Simple approximate forms based on charge and geometry coordinates are proposed for the high-order moments of each atom, which are further required for the calculation of the electrostatic field within the structures. The possibility to use this analytical approach to evaluate the electrostatic field within embedded cluster models is also shortly discussed.

Separation of gas mixtures using a range of zeolite membranes: A molecular-dynamics study