To achieve high curing rates without introducing artificially large distance cutoffs, we implemented a mixing procedure during which no cross-linking reaction was allowed. This involved an annealing run (Figure ) followed by an equilibration run (see Figure for an example iterative run). Specifically, for curing degrees below 90%, after every curing cycle the temperature was increased to 800 K over 30 ps in the NVT ensemble, kept constant in the NPT ensemble for 500 ps, and then reduced again to 503 K. When exceeding 90% curing, the annealing time in the NPT ensemble was progressively increased up to 2.5 ns, in 500-ps increments, until 95% conversion was achieved. Each annealing round was followed by an equilibration step involving 100 ps of NPT equilibration at 503 K, to let the density stabilize before the next stage of cross-linking (Figure ).