Estimation of Partial Charges in Small Zeolite Imidazolate Frameworks from Density Functional Theory Calculations

Zeolitic Imidazolate Frameworks (ZIFs) are the new frontier in the field of metal-organic materials. They incorporate the confining properties of the more traditional aluminosilicate zeolites together with the catalytic activity provided by transition metal ions and organic links. Computation of atomic point charges for these hybrid materials is important in the field of molecular simulations for substantial prediction of experimental results. However, due to the structural complexity of advanced materials in general, studies involving derivation of point charges for these materials are truly scarce. In this article, we have derived the atomic point charges of ZIF-8 through fitting of the quantum electrostatic potential obtained systematically from density functional theory (DFT) calculations both on finite clusters of increasing size and on the periodic system. For the periodic system, fluctuations on the atomic charges have been studied through ab initio molecular dynamics simulations. Using the latter approach, we have extended the study to ZIF-2 and ZIF-3, where it has been found that charge fluctuations are, as well as for ZIF-8, very narrow, therefore justifying the use of the point charge approximation for these materials.