Dear all,

I find in the discussion list that we can use kwant.plot to plot the hopping strength,

def hopping_strength(site1, site2):
    return sys.hamiltonian(site1, site2)
kwant.plot(sys, hop_lw=hopping_strength)

Can we plot the onsite energy of the system similarly?

Regards

Khani Hosein