Dear all,

I am using Kwant to discretize an effective mass approximation scheme.  Efforts to implement position dependent effective mass to model a metal-semiconductor hetrojunction have so far been futile. To avoid having to discretize the entire system with the high accuacy needed in the metal parts (the leads) i want to have different discretisation in different parts of the system.

So my question is if it is possible to have such an interface between discretizations and if so, how could it be implemented in Kwant?

Thank you in advance.

Emil Johansson

PS: Have tried the straight on approach with having a polyatomic lattice and in one part only having one of the monoatomic sublattices present. This however do no generate results in agreement with analytical results for a simple 2D wire test system. Attached a plot of the test system.