Efficient computation of Kubo conductivity for incommensurate 2D heterostructures

Abstract

Here we introduce a numerical method for computing conductivity via the Kubo Formula for incommensurate 2D bilayer heterostructures using a tight-binding framework. We begin with deriving the momentum space formulation and Kubo Formula from the real space tight-binding model using the appropriate Bloch transformation operator. We further discuss the resulting algorithm along with its convergence rate and computation cost in terms of parameters such as relaxation time and temperature. In particular, we show that for low frequencies, low temperature, and long relaxation times conductivity can be computed very efficiently using momentum space for a wide class of materials. We then demonstrate our method by computing conductivity for twisted bilayer graphene (tBLG) for small twist angles.