The geometry was locally optimised using scipy, applying a harmonic potential to the bonds and no angular potential in the interests of computational speed . A single bond switching move introduces a defect known as a Stone-Wales defect into a hexagonal network, and is shown schematically in Fig. 20 . Stone-Wales defects have been studied in great detail as a major type of defect in 2D silica glasses .