Action-derived molecular dynamics simulations for the migration and coalescence of vacancies in graphene and carbon nanotubes

Abstract

We report the results of action-derived molecular dynamics simulations for the migration and coalescence processes of monovacancies in graphene and carbon nanotubes with different chiralities. In carbon nanotubes, the migration pathways and barriers of a monovacancy depend on the tube chirality, while there is no preferential pathway in graphene due to the lattice symmetry and the absence of the curvature effect. The probable pathway changes from the axial to circumferential direction as the chirality varies from armchair to zigzag. The chirality dependence is attributed to the preferential orientation of the reconstructed bond formed around each vacancy site. It is energetically more favourable for two monovacancies to coalesce into a divacancy via alternative movements rather than simultaneous movements. The energy barriers for coalescence are generally determined by the migration barrier for the monovacancy, although there are some variations due to interactions between two diffusing vacancies. In graphene and armchair nanotubes, two monovacancies prefer to migrate along different zigzag atomic chains rather than a single atomic chain connecting these vacancies. On the other hand, in zigzag tubes, the energy barrier for coalescence increases significantly unless monovacancies lie on the same circumference.