All ReaxFF calculations were performed on an eight-layer 10x10 (100) slab, where the two lowermost layers have been kept in their bulk structure. The entire outgoing row of atoms was occupied for each occupied surrounding site when generating the surface structure. In total, 2 10 = 1024 hopping and 2 14 = 16384 exchange processes are included in our kMC simulation.