We therefore adapted our procedure to produce glass models, by running a simulation at a sufficiently high temperature to observe structural rearrangement while keeping the simulation duration short enough to avoid the unphysical events. We carried an MD melt-quenching simulation with a maximal temperature (T max ) of 1400 K, maintaining the simulation at this temperature for 100 ps: we observed melting while witnessing no imidazolate breaking during the process.