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

TL;DR

This study combines density-functional and Green's function methods to analyze the electronic and magnetic properties of a substitutional vacancy in graphene, revealing localized states, magnetic moments, and long-range wave functions.

Contribution

It provides a detailed theoretical analysis of vacancy-induced electronic states and magnetism in graphene using combined DFT and Green's function approaches, highlighting the long-range nature of localized states.

Findings

Localized mid-gap states (Vσ) split by crystal field and Jahn-Teller distortion

Sharp resonance state (Vπ) in the band structure

Magnetic moment of approximately 1.7 μB due to localized electrons

Abstract

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.