In order to offer a meaningful comparison, we considered several glass models: the same ab initio glass models studied with AIMD, and MLP glasses obtained by melt-quenching in the previous section. For consistency between ab initio and MLP glasses, we used for both 10 single cell models and averaged results. These models are thoroughly presented in the supporting information. We checked that the use of a supercell did not seem to significantly impact the resulting bulk moduli K and density at zero pressure ρ 0 , as exemplified on Table in the case of the ZIF-4 crystal. The P -V plot with the standard deviations as well as the fitted EoS is shown on Figure in the case of the crystal. As shown on Table and Figure , we obtain a bulk modulus of K crystal = 1.69 GPa, in reasonable agreement with AIMD (1.39 GPa ) and previous experimental works -which found a value of 2.01 GPa with high-pressure crystallography 62 and 1.42 GPa with mercury intrusion. Further interpretation is limited, as the precise value determined depends on the computational methodology and -for experiments -also on subtle differences between MOFs from different batches or on the experimental setup. However, the comparison between two different systems within the same methodology is meaningful, and calculations with the NequIP MLP reproduced the trend previously identified with other MD schemes, where glasses have larger K than the crystalline phase.