Although classical MD simulations are routinely used for crystalline MOFs, they are unable by design to simulate bond breaking or reformation, a key mechanism in the formation of ZIF-4 glasses by melt-quenching. To develop a method able to screen the mechanical properties of a large number of ZIF glasses, it would be convenient and consistent to use the same MD scheme for the glass model generation and the subsequent determination of its