The Zhou potential [22] was used in this study. For each macroscopic grain boundary structure, the energies of 100 to 10000 initial configurations were minimized using the conjugate-gradient method in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code [23] at 0 K