Hi all,

I want to plot the the Charge density using kwant.plot method, I want to know how to plot the Charge density in a system with multiple orbitals per site.

Here is some pesuo codes:

..... syst = kwant.Builder() ..... model= kwant.wraparound.wraparound(syst).finalized() ham=model.hamiltonian_submatrix

(eigen,evecs) = la.eigh(ham)

def site_size(i): return 2 * wf[i] / wf.max()

kwant.plot(model, site_size=site_size, site_color=(0, 0, 1, 0.3),fig_size=(10,10)); ********* In my codes, I just found len(model.sites)=6, that is the number of sites in the unit cell, however, due to there are multiple orbitals per site, the len(evecs[:,i])=22 orbitals. So I want to know the real infortion when using kwant.plot? Like what I just posted, it only plot evecs[:6,i]**2 functions?

Thanks Xiangru