The simulations were performed in the NVT ensemble where the temperature was slowly cooled from k B T = 1.12125 over the transition point k B T ≈ 1 with temperature steps of T i+1 = T i × 0.975. The simulations were run for 8 hours using 64 cores on Imperial College London's Research Computing Service facilities which resulted in 9 temperatures being sampled. For the first 8 temperatures, 10 8 timesteps were sampled, where a random number of particles are moved in the cell for each timestep. As the time for each timestep to be completed depends on the number of interactions being calculated, once the polyhedra form a cluster, the time taken to complete each timestep increases. Due to the finite wall time of the simulations, the simulations were terminated before finishing all time steps. The cluster outputted on termination is considered to be representative of the phase behaviour at all lower temperatures.