Is it possible to obtain the edge or surface states of a TI from its LDOS on a semi-infinite slab? I'm attaching a single lead to a scattering region consistent of a single layer of the same H as the lead, then I ask for the LDOS at the scattering region, which would be the surface of the semi-infinite system.
The problem is that kwant uses the propagating modes of the lead to calculate the LDOS, and consequently I only get the surface states on energy ranges where the bulk bands exist, but I'm interested in the LDOS at the surface and within the bulk gap.
