During melting with a 2.5 K/ps heating rate, we report that when the maximal temperature exceeds 1300 K for more than a few tens of ps, the system vaporizes (as shown on Fig- ). This effect is kinetically hidden when using the faster heating rates of 24 -96 K/ps that do not provide sufficient time for the system to equilibrate, explaining why it was not observed in the previous ReaxFF works. It suggests a tendency of this reactive force field to be overly favorable to zinc-imidazolate bond breaking events, which is exacerbated at such high temperatures. Changing the maximal temperature to 1100 K to avoid these undesired effects has significant implications for the computed structural properties. As shown on Table , the decoordination of the Zn atoms compared to the crystal is now in better agreement with ab initio data 10 than in previous ReaxFF works. Similarly, the Zn-N potential of mean force (PMF)