Semi-Empirical Pseudopotential Method for Graphene and Graphene Nanoribbons

We implemented a semi-empirical pseudopotential (SEP) method for calculating the band structures of graphene and graphene nanoribbons. The basis functions adopted are two-dimensional plane waves multiplied by several B-spline functions along the perpendicular direction. The SEP includes both local and non-local terms, which were parametrized to fit relevant quantities obtained from the first-principles calculations based on the density-functional theory (DFT). With only a handful of parameters, we were able to reproduce the full band structure of graphene obtained by DFT with a negligible difference. Our method is simple to use and much more efficient than the DFT calculation. We then applied this SEP method to calculate the band structures of graphene nanoribbons. By adding a simple correction term to the local pseudopotentials on the edges of the nanoribbon (which mimics the effect caused by edge creation), we again obtained band structures of the armchair nanoribbon fairly close to the results obtained by DFT. Our approach allows the simulation of optical and transport properties of realistic nanodevices made of graphene nanoribbons with very little computation effort.