The choice of the maximal temperature, necessarily above the melting temperature, is a trade-off between the need to gather statistics on relatively rare events during the relatively short times explored in MD, and the necessity to preserve the physical consistency of the model. A value of 1500 K was chosen for both the previous ab initio 10 and ReaxFF works. While ab initio MD simulations of the liquid ZIF have been performed up to 2000 K while preserving the integrity of the imidazolate linkers, no published work investigated the temperature stability of ReaxFF simulations of the ZIF liquid.