Point defects in twisted bilayer graphene: A density functional theory study

We have used ab initio density functional theory, incorporating van der Waals corrections, to study twisted bilayer graphene (TBLG) where Stone-Wales defects or monovacancies are introduced in one of the layers. We compare these results to those for defects in single-layer graphene or Bernal stacked graphene. The energetics of defect formation is not very sensitive to the stacking of the layers or the specific site at which the defect is created, suggesting a weak interlayer coupling. However, signatures of the interlayer coupling are manifested clearly in the electronic band structure. For the gamma gamma Stone-Wales defect in TBLG, we observe two Dirac cones that are shifted in both momentum space and energy. This up/down shift in energy results from the combined effect of a charge transfer between the two graphene layers and a chemical interaction between the layers, which mimics the effects of a transverse electric field. Charge density plots show that states near the Dirac points have significant admixture between the two layers. For Stone-Wales defects at other sites in TBLG, this basic structure is modified by the creation of minigaps at energy crossings. For a monovacancy, the Dirac cone of the pristine layer is shifted up in energy by similar to 0.25 eV due to a combination of the requirement of the equilibration of Fermi energy between the two layers with different numbers of electrons, charge transfer, and chemical interactions. Both kinds of defects increase the density of states at the Fermi level. The monovacancy also results in spin polarization, with magnetic moments on the defect of 1.2-1.8 mu(B).