Finally, to obtain additional mechanical properties, we applied the strain-fluctuation but found that it fails to highlight the difference in K between the crystal and the glasses as shown on Figure , which may be explained by the unreasonably large density of the crystal Despite this challenge, all finite difference methods are overall consistent, as illustrated on Figure , albeit with a large variability across glass configurations. Both FF adaptations yield similar results, suggesting they are not ill adapted and that this procedure is not very sensitive to such a small perturbation in the MOF-FF parametrization. We thus showed that the finite difference methods with classical FFs are ready to be deployed, and that the development of a specific MOF-FF parametrization for amorphous ZIFs could prove useful.