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