In prior work, AIMD has successfully been used to model the melting of the ZIF-4 crys-tal, and to produce atomistic configurations of melt-quenched glasses. All these models have been generated and equilibrated in the (N, V, T ) ensemble, and thus are all at the same density (taken to be the crystallographic density 1.22 g cm -3 ). Although this difficulty in capturing changes in density has been identified as an important limitation of the ab initio methodology, such glass models still represent the most chemically accurate atomistic description of ZIF-4 glasses published to date in the literature. Here, we have computed the bulk modulus K for four of these systems: the crystal and three melt-quenched glasses, chosen to preserve the diversity found in the coordination environments. As only three different discrete values were found for the Zn-N coordination numbers of the 10 original glass models, we chose one model per value. This restriction was guided by the significant computational cost of ab initio MD.