Density functional study of graphene overlayers on SiC

Despite the ongoing graphene boom of the last three years our understanding of epitaxial graphene grown on SiC substrate is only beginning to emerge. Along with experimental methods such as low energy electron diffraction (LEED), scanning tunneling microscopy (STM) and angle resolved photoemission spectroscopy (ARPES), ab initio calculations help to uncover the geometric and electronic structure of the graphene/SiC interface. In this chapter we describe the density-functional calculations we performed for single and double graphene layers on Si- and C-terminated 6H-SiC surfaces. Experimental data reveal a pronounced difference between the two surface terminations. On a Si-terminated surface the interface adopts a 6 root 3 x 6 root 3 unit cell whereas the C-face supports misoriented (turbostratic) graphene layers. It has been recently realized that, on the Si-face, the large commensurate cell is subdivided into patches of coherently matching to the substrate carbon atoms. In our calculations we assumed the coherent match geometry for the whole interface plane. This reduces the periodic unit to the root 3 x root 3R30 degrees cell but requires a substantial stretching of the graphene sheet. Although simplified, the model provides a qualitative picture of the bonding and of the interface electron energy spectrum. We find that the covalent bonding between the carbon layer and the substrate destroys the massless relativistic electron energy spectrum, the hallmark of a freestanding graphene. Hence the first carbon layer cannot be responsible for the graphene-type electron spectrum observed by ARPES and rather plays a role of a buffer between the substrate and the subsequent carbon sheets. The true graphene spectrum appears with the second carbon layer which exhibits a weak van der Waals bonding to the underlying structure. For Si-terminated substrate, we find that the Fermi level is pinned by the interface state at 0.45 eV above the graphene Dirac point, in agreement with experimental data. This renders the interface metallic. On the contrary, for a C-face the coherent match model predicts the Fermi level exactly at the Dirac point. However, this does not necessarily apply to the turbostratic graphene layers that normally grow on the C-terminated substrate. (C) 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.