These categories of defects can be combined, performing multiple bond switching moves and edge removals to create more complex defect structures, for example the combined bond switch and edge removal shown in Fig. . Furthermore, the combined defects can be correlated or uncorrelated with one another to lead to more complex scenarios; in this work we choose uncorrelated defects unless otherwise stated. However, the introduction of multiple defects in the graph representation required optimisation to a local energy minimum. At nodes, the energy of an edge was given by a simple harmonic potential