alex.black9980@gmail.com wrote:
I have imaginary parts in my hopping due to the spin-orbit coupling and then tried to calculate using ldos according to your answer, but it gives an error as : ldos for non-Hermitian Hamiltonians is not implemented yet.
Right, I forgot about that. As far as I can reconstruct we have disabled LDOS computation for non-Hermitian Hamiltonians because we considered it ill-defined or at least dangerous. (But I may be wrong, it was a long time ago.) LDOS in Kwant is computed as the sum over the squared absolute values of the wave functions that are the response of all the incoming modes at a given energy. You should be able to compute the same thing by calling ’kwant.wave_function’ and summing over abs(...)**2 for all modes and all leads. (Do check that you get the same result as the built-in 'ldos’ function for Hermitian Hamiltonians.) Is this what you would like to compute? I am not sure whether computing LDOS in this above way makes much physical sense, but it should be possible. Christoph