Electronic structure of the substitutional vacancy in graphene: Density-functional and Green's function studies

We study the electronic structure of graphene with a single substitutional vacancy using a combination of the density-functional, tight-binding, and impurity Green's function approaches. Density functional studies are performed with the all-electron spin-polarized linear augmented plane wave (LAPW) method. The three $sp^2 \sigma$ dangling bonds adjacent to the vacancy introduce localized states (V$\sigma$) in the mid-gap region, which split due to the crystal field and a Jahn-Teller distortion, while the $p_z \pi$ states introduce a sharp resonance state (V$\pi$) in the band structure. For a planar structure, symmetry strictly forbids hybridization between the $\sigma$ and the $\pi$ states, so that these bands are clearly identifiable in the calculated band structure. As for the magnetic moment of the vacancy, the Hund's-rule coupling aligns the spins of the four localized V$\sigma_1 \uparrow \downarrow$, V$\sigma_2 \uparrow$, and the V$\pi \uparrow$ electrons resulting in a S=1 state, with a magnetic moment of $2 \mu_B$, which is reduced by about $0.3 \mu_B$ due to the anti-ferromagnetic spin-polarization of the $\pi$ band itinerant states in the vicinity of the vacancy. This results in the net magnetic moment of $1.7 \mu_B$. Using the Lippmann-Schwinger equation, we reproduce the well-known $\sim 1/r$ decay of the localized V$\pi$ wave function with distance and in addition find an interference term coming from the two Dirac points, previously unnoticed in the literature. The long-range nature of the V$\pi$ wave function is a unique feature of the graphene vacancy and we suggest that this may be one of the reasons for the widely varying relaxed structures and magnetic moments reported from the supercell band calculations in the literature.