Dear all I was trying to plot the wavefunction on each site of a 3-D lattice and each site is described by a 4x4 matrix. However, it turned out the that plotted wave function is overlapped with the surrounding site and I could not identify the magnitude of wave function on each site. The code I had implemented for plotting was the same as the codes in " http://kwant-project.org/doc/1.0/tutorial/tutorial6" with a copy of specific plot command as below. I have varied the value "3" before wf[i] to tune the size however it turned out this value cannot tune the plot size of the wavefunction. I think I'm not doing it correctly. def site_size(i): return 3 * wf[i] / wf.max() kwant.plot(sys, site_size=site_size, site_color=(0, 0, 1, 0.3), hop_lw=0.1) This method works perfectly for 2-D lattice plot. However, it cannot display the wavefunction properly in 3-d lattice. My question is, how can I plot the wavefunction in 3-d lattice clearly? (Or maybe I can force the kwant to plot it in 2-d by configuring the kwant.plot?) Thank you for advance. Sincerely, Liming