Because of their wide diversity and tunability, both structurally and chemically, metalorganic frameworks (MOFs) have been proposed over the past two decades for a large variety of applications such as gas storage, fluid separation, and heterogeneous catalysis. Out of this vast family of nano-and meso-porous materials, composed of inorganic nodes connected by organic linkers to form three-dimensional architectures, several classes demonstrated high chemical and thermal stability, 2 such as zeolitic imidazolate frameworks (ZIFs). However, while their mechanical stability is also essential for MOFs to fully achieve their potential in industrial-scale processes -in particular for processing (extrusion and pellet formation) -, the study of how these materials respond to mechanical stress is comparatively still emerging. Particularly promising in this regard, is the growing family of non-crystalline MOFs, coming as a large variety of MOF glasses and MOF gels (forming aerogels or monoliths upon drying) and which can be formed by various production routes. While conserving many of the intrinsic advantages of their crystalline counterparts, these states yield the potential of increased mechanical robustness and would allow for greater ease of processing, notably by circumventing the performance drop due to the necessary densification of the MOF powders. While a series of experimental techniques have successfully been employed to determine the mechanical properties of crystalline MOFs, there is a lack of studies on the amorphous phases. Some methods such as those relying on high-pressure X-ray diffraction are not straightforwardly applicable to disordered materials, and many others require large bulk glass samples which are challenging to prepare. Fortunately, the same mechanical properties can be studied by computational simulations, which have been used extensively as a complement to experiments on MOF crystals and can lead to systematic studies of the structure-property relationships. Despite the challenges in the determination of the framework structure of amorphous states at the microscopic scale, several methodologies are available to generate the atomistic models which are a prerequisite of any computational determination of the mechanical properties. One family of approaches to computing finite temperature mechanical properties consists in using Molecular Dynamics (MD), either to mimic in silico the high-pressure experiments or by using the strain-fluctuation method. Unlike most computational works which led to the determination of mechanical properties at 0 K, often using the stress-strain approach, the temperature dependence of these methods allows for a better prediction and understanding of the mechanical strength and stability of the studied materials, as the mechanical properties at operational temperatures may deviate significantly from those computed at 0 K. It is key for disordered MOFs which can display complex pressure and temperature behaviors at both low and high tem-perature, with the occurence of multiple phase transitions. In this work, we applied a range of methods and MD schemes -ab initio, reactive force fields, classical force fields -to extract the finite temperature mechanical properties of ZIF crystals and glasses. We first report how ab initio MD can be used to obtain reference values, which have never been determined for a MOF glass in the literature before. We then investigate two MD schemes with lower computational cost: classical and reactive force fields (FF). We conclude that reactive force fields are inapplicable to low-density models, while classical force fields reproduce with good agreement ab initio results. We find that ZIF-4 glasses have a larger bulk modulus than the crystal, and confirm previous studies that the density is larger for the glass phases. In addition, we confirm the persistent relationship between density and bulk modulus, visible across every model, MD scheme and method.