Hierarchical Crystal Structure Prediction of Zeolitic Imidazolate Frameworks Using DFT and Machine-Learned Interatomic Potentials

TL;DR

This study employs crystal structure prediction combined with machine-learned interatomic potentials to explore and identify known and novel structures of zinc imidazolate frameworks, aiding experimental design.

Contribution

It introduces a high-throughput CSP approach using MLIPs to discover new polymorphs and match experimental structures in ZIFs, expanding structural understanding.

Findings

Identified 9609 energy minima with 1484 network topologies.

Matched all but one known ZnIm2 structures within the search.

Discovered 855 previously unreported topologies.

Abstract

Crystal structure prediction (CSP) is emerging as a powerful method for the computational design of metal-organic frameworks (MOFs). In this article we employ CSP to perform high-throughput exploration of the crystal energy landscape of zinc imidazolate (ZnIm2). As the most polymorphic member of the zeolitic imidazolate framework (ZIF) family, ZnIm2 has at least 24 reported structural and topological forms, and new polymorphs still being regularly discovered. With the aid of custom-trained machine-learned interatomic potentials (MLIPs) we have performed a high-throughput sampling of over 3 million randomly-generated crystal packing arrangements and identified 9609 energy minima characterized by 1484 network topologies, including 855 topologies that have not been reported before. All but one experimentally-reported structures of ZnIm2, falling within the search boundaries, were ultimately matched with the predicted structures, demonstrating the power of the CSP method in sampling experimentally-relevant ZIF structures. Finally, through a combination of topological analysis, density and porosity considerations, we have identified a set of structures representing promising targets for future experimental screening. as well as demonstrated how structures of mechanochemically-synthesized MOFs could be identified via matching experimental powder diffraction patterns with simulated patterns from the predicted structures.