Nucleation of zeolitic imidazolate frameworks: from molecules to nanoparticles

We have studied the clusters involved in the initial stages of nucleation of Zeolitic Imidazolate Frameworks, employing a wide range of computational techniques. In the pre-nucleating solution, the prevalent cluster is the ZnIm₄ cluster (formed by a zinc cation, Zn²⁺, and four imidazolate anions, Im^-), although clusters such as ZnIm₃, Zn₂Im₇, Zn₂Im₇, Zn₃Im₉, Zn₃Im₁₀, or Zn₄Im₁₂ have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm₃, have a tetrahedrally coordinated Zn²⁺ cation. Small ZnIm clusters are less stable than the ZnIm₄ cluster. The first cluster that is found to be more stable than ZnIm₄ is the Zn₄₁Im₈₈ cluster, which is a disordered cluster with glassy structure. Bulk-like clusters do not begin to be more stable than glassy clusters until much larger sizes, since the larger cluster we have studied (Zn₁₄₄Im₂₈₈) is still less stable than the glassy Zn₄₁Im₈₈ cluster, suggesting that Ostwald's rule (the less stable polymorph crystallizes first) could be fulfilled, not for kinetic, but for thermodynamic reasons. Our results suggest that the first clusters formed in the nucleation process would be glassy clusters, which then undergo transformation to any of the various crystal structures possible, depending on the kinetic routes provided by the synthesis conditions. Our study helps elucidate the way in which the various species present in solution interact, leading to nucleation and crystal growth.