Structure–Property Relationship of Piezoelectric Properties in Zeolitic Imidazolate Frameworks: A Computational Study

Metal–organic frameworks (MOFs) are a class of nanoporous crystalline materials with very high structural tunability. They possess a very low dielectric permittivity ε_r due to their porosity and hence are favorable for piezoelectric energy harvesting. Even though they have huge potential as piezoelectric materials, a detailed analysis and structure–property relationship of the piezoelectric properties in MOFs are lacking so far. This work focuses on a class of cubic non-centrosymmetric MOFs, namely, zeolitic imidazolate frameworks (ZIFs) to rationalize how the variation of different building blocks of the structure, that is, metal node and linker substituents affect the piezoelectric constants. The piezoelectric tensor for the ZIFs is computed from ab initio theoretical methods. From the calculations, we analyze the different contributions to the final piezoelectric constant d₁₄, namely, the clamped ion (e₁₄⁰) and the internal strain (e₁₄^int) contributions and the mechanical properties. For the studied ZIFs, even though e₁₄ (e₁₄⁰ + e₁₄^int) is similar for all ZIFs, the resultant piezoelectric coefficient d₁₄ calculated from piezoelectric constant e₁₄ and elastic compliance constant s₄₄ varies significantly among the different structures. It is the largest for CdIF-1 (Cd²⁺ and −CH₃ linker substituent). This is mainly due to the higher elasticity or flexibility of the framework. Interestingly, the magnitude of d₁₄ for CdIF-1 is higher than II–VI inorganic piezoelectrics and of a similar magnitude as the quintessential piezoelectric polymer polyvinylidene fluoride.