Efficient computation of Kubo conductivity for incommensurate 2D heterostructures

We introduce a numerical method for computing conductivity via the Kubo formula for incommensurate 2D bilayer heterostructures using a tight-binding framework. We begin by 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 computational cost in terms of model 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 the momentum space algorithm for a wide class of materials. We then showcase our method by computing conductivity for twisted bilayer graphene (tBLG) for small twist angles.