properties. An option is the use of reactive force fields, which have been proposed to generate models of ZIF glasses by melt-quenching with the development of a ReaxFF parametrization for ZIFs. They have previously been used to compute mechanical properties and are a natural candidate for the task at hand. Able to explore larger spatial and time scales, they could in principle be used to obtain the bulk modulus with several methods while reducing the finite size effects inherent to AIMD. However, we showed in a previous work the atypical structural properties of the glasses obtained with this approach, and there is still no in-depth validation of the mechanical properties obtained with this force field in the existing literature. In order to offer a meaningful comparison to other MD schemes, we consider several glass models: a ReaxFF glass obtained in our previous work, the same three ab initio glass models studied in the previous section and a glass obtained by Reverse Monte Carlo (RMC) modeling. These models, as well as the crystalline system, went through a careful preparation process and are thoroughly presented in the supporting information.