New all-atom force field for molecular dynamics simulation of an AlPO(4)-34 molecular sieve

A force field of the triclinic framework of AlPO(4)-34, important in methanol-hydrocarbon conversion reactions, was developed using an empirical potential function. Molecular dynamics simulation of an AlPO(4)-34 triclinic framework segment of 1216 atoms, containing the template molecules isopropylamine and water, was performed with explicit consideration of atomic charges. The average RMS difference between instantaneous positions of the framework atoms during 1 ns simulation and their positions in the structure determined from single crystal X-ray diffraction was calculated, and the average structure of the flexible framework was determined. The computed Debye-Waller factors and simulated FTIR spectra are in good agreement with the experimental data. The new force field permits detailed molecular dynamics simulations of flexible, charged aluminophosphate molecular sieves which should lead to a better understanding of the catalytic processes and the crucial role played by templating molecules.