So, I know that to calculate surface band structure of a material, you have to take its 2D system eg. single or few-layered. But, how to find the surface band structure in the case of heterostructures. This paper does it "Proximity-induced topological state in graphene " but uses VASP instead. Is it possible to do the same using Quantum espresso?