Dear Sir,
I was trying to plot 2D LDOS map for a 3D system. I am using
kwant.plotter.map and kwant.plotter.plot for plotting density and LDOS.
But there is an error message showing " Only 2D system can be plotted in this
way". The following part is included where I am getting this error. Thanks in
advance.
def density(sys, energy, lead_nr, params=params):
wf = kwant.wave_function(sys, energy, params=params)
return (abs(wf(lead_nr))**2).sum(axis=0)
d = density(sys,0.5, 0, params=params)
print (d)
kwant.plotter.map(sys, d, cmap='gist_heat_r', colorbar=True)#, site_lw=0)
ldos = kwant.ldos(sys, energy=0.5, params=params)
print (ldos)
kwant.plotter.plot(sys, ldos, num_lead_cells=2)
Regards,
priyanka
