Dear all,

I want to plot the local density of states (LDOS) for a honeycomb system of a rectangular shape modeled by Hamiltonian H.

Because the Bravais lattice vectors are non-orthogonal, the honeycomb leads (but modeled by different Hamiltonian H') are attached in such a way, that sites have to be added to the scattering region. This added part, however, is modeled by H' whereas the rest (original scattering region) by H.

I know that I can just say that the added part is actually part of the leads (matter of definition).

But, if I want to plot LDOS using "kwant.ldos(sys, energy)", it plots LDOS also on this added part modeled by H', which obscures the LDOS of my interest on the original scattering region modeled by H.

How to make KWANT not to plot the LDOS on the artificially added part to the scattering region?

Thanks for an advice,
Tibor