Towards Predicting the Sequential Appearance of Zeolitic Imidazolate Frameworks Synthesized by Mechanochemistry

We use computational materials methods to study the sequential appearance of zinc-based zeolitic imidazolate frameworks (ZIFs) generated in the mechanochemical conversion process. We consider nine ZIF topologies, namely RHO, ANA, QTZ, SOD, KAT, DIA, NEB, CAG and GIS, combined with the two ligands 2-methylimidazolate and 2-ethylimidazolate. Of the 18 combinations obtained, only six (three for each ligand) were actually observed during the mechanosynthesis process. Energy and porosity calculations based on density functional theory, in combination with the Ostwald rule of stages, were found to be insufficient to distinguish the experimentally observed ZIFs. We then show, using classical molecular dynamics, that only ZIFs withstanding quasi-hydrostatic pressure P >= 0.3 GPa without being destroyed were observed in the laboratory. This finding, along with the requirement that successive ZIFs be generated with decreasing porosity and/or energy, provides heuristic rules for predicting the sequences of mechanically generated ZIFs for the two ligands considered.

Automated summary

This study investigates the sequential appearance of zinc-based zeolitic imidazolate frameworks (ZIFs) in the mechanochemical conversion process using computational materials methods. The results show that energy and porosity calculations alone are insufficient to distinguish the experimentally observed ZIFs, while only ZIFs able to withstand quasi-hydrostatic pressure and possess lower porosity and/or energy were observed in the laboratory. These findings provide heuristic rules for predicting the sequences of mechanically generated ZIFs.