Abstract
The formation of extended defects in graphene from the coalescence of individual mobile vacancies can significantly alter its mechanical, electrical and chemical properties. We present the results of ab initio simulations which demonstrate that the strain created by multi-vacancy complexes in graphene determine their overall growth morphology when formed from the coalescence of individual mobile lattice vacancies. Using density functional theory, we map out the potential energy surface for the motion of mono-vacancies in the vicinity of multi-vacancy defects. The inhomogeneous bond strain created by the multi-vacancy complexes strongly biases the activation energy barriers for single vacancy motion over a wide area. Kinetic Monte Carlo simulations based on rates from ab initio derived activation energies are performed to investigate the dynamical evolution of single vacancies in these strain fields. The resultant coalescence processes reveal that the dominant morphology of multi-vacancy complexes will consist of vacancy lines running in the two primary crystallographic directions, and that more thermodynamically stable structures, such as holes, are kinetically inaccessible from mono-vacancy aggregation alone.