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
