We study theoretically the effects of short-range electron-electron interactions on the electronic structure of graphene, in the presence of single substitutional impurities. Our computational approach is based on the $pi$ orbital tight-binding appro
ximation for graphene, with the electron-electron interactions treated self-consistently at the level of the mean-field Hubbard model. We compare explicitly non-interacting and interacting cases with varying interaction strength and impurity potential strength. We focus in particular on the interaction-induced modifications in the local density of states around the impurity, which is a quantity that can be directly probed by scanning tunneling spectroscopy of doped graphene. We find that the resonant character of the impurity states near the Fermi level is enhanced by the interactions. Furthermore, the size of the energy gap, which opens at high-symmetry points of the Brillouin zone of the supercell upon doping, is significantly affected by the interactions. The details of this effect depend subtly on the supercell geometry. We use a perturbative model to explain these features and find quantitative agreement with numerical results.
The response of thin films of Bi$_2$Se$_3$ to a strong perpendicular magnetic field is investigated by performing magnetic bandstructure calculations for a realistic multi-band tight-binding model. Several crucial features of Landau quantization in a
realistic three-dimensional topological insulator are revealed. The $n=0$ Landau level is absent in ultra-thin films, in agreement with experiment. In films with a crossover thickness of five quintuple layers, there is a signature of the $n=0$ level, whose overall trend as a function of magnetic field matches the established low-energy effective-model result. Importantly, we find a field-dependent splitting and a strong spin-polarization of the $n=0$ level which can be measured experimentally at reasonable field strengths. Our calculations show mixing between the surface and bulk Landau levels which causes the character of levels to evolve with magnetic field.
We report on microscopic tight-binding modeling of surface states in Bi$_2$Se$_3$ three-dimensional topological insulator, based on a sp$^3$ Slater-Koster Hamiltonian, with parameters calculated from density functional theory. The effect of spin-orbi
t interaction on the electronic structure of the bulk and of a slab with finite thickness is investigated. In particular, a phenomenological criterion of band inversion is formulated for both bulk and slab, based on the calculated atomic- and orbital-projections of the wavefunctions, associated with valence and conduction band extrema at the center of the Brillouin zone. We carry out a thorough analysis of the calculated bandstructures of slabs with varying thickness, where surface states are identified using a quantitative criterion according to their spatial distribution. The thickness-dependent energy gap, attributed to inter-surface interaction, and the emergence of gapless surface states for slabs above a critical thickness are investigated. We map out the transition to the infinite-thickness limit by calculating explicitly the modifications in the spatial distribution and spin-character of the surface states wavefunction with increasing the slab thickness. Our numerical analysis shows that the system must be approximately forty quintuple-layers thick to exhibit completely decoupled surface states, localized on the opposite surfaces. These results have implications on the effect of external perturbations on the surface states near the Dirac point.
