Dear Kwant Users,
I am trying to calculate the bulk density of states for an infinite(in x) homogenous 2D superconductor with spin orbit and Zeeman field with 3 sites in the y direction. I have defined the leads and the system according to the superconducting hamiltonian. For the same, I calculate the dos for the first layer(in y) using the following piece of code: def ldos(syst,energy): ldos1=kwant.ldos(syst, energy) #print(ldos) X1 = np.sum(ldos1[0:4]) return X1 Now, the result I get is physically correct, although I would want the DOS to be a little smoother. In the traditional way of doing things, a retarded greens function would have a small imaginary part, increasing which would give me a smoother curve. Is there any way in kwant, that I can control such a broadening parameter?
Thanking you in advance, Wils