As the identification of the building units (Zn atoms and imidazolates) of the glass models could not be performed by molsys, due to large fluctuations in the bond angles and lengths as well as the presence of defects, it was carried out with our python library aMOF , 40 interfaced with molsys. aMOF is available online on , where the version used for this work is v1.0.0.