Quantitatively
modeling adsorbate diffusion through zeolitic imidazolate
frameworks (ZIFs) must account for the inherent flexibility of these
materials. The lack of a transferable intramolecular ZIF force field
(FF) for use in classical simulations has previously made an accurate
simulation of adsorbate diffusion in many ZIFs impossible. We resolve
this problem by introducing a density functional theory parameterized
force field (FF) for ZIFs named the intraZIF-FF, which includes perturbations
to the class I force fields previously used to model ZIFs. This FF
outperforms ad hoc force fields at predicting ab initio relative energies
and atomic forces taken from fully periodic ab initio molecular dynamics
simulations of SALEM-2, ZIF-7, ZIF-8, and ZIF-90. We use the intraZIF-FF
to predict the infinite dilution self-diffusion coefficients of 30
adsorbates with molecular diameters ranging from 2.66 to 7.0 Å
in these 4 ZIFs. These results greatly increase the number of adsorbates
for which accurate information about molecular diffusion in ZIFs is
available.