Hi Derek,
The code to plot the wave function is: wf = kwant.wave_function(sys,energy=0)(0) rho = kwant.operator.Density(sys) density = rho(wf) kwant.plotter.density(sys, density)
When I do print(np.shape(wf)) I get (0, 300); where 300 is the number of sites in this example lattice. Thus it seems like the wave function information might be contained in the returned object, but I am not sure how to access it. Please let me know if you have any idea on what I am doing incorrectly here
The size of the first axis of 'wf' is zero because there are no propagating modes originating from lead #0 at energy 0 in your system. Try plotting the band structure of lead #0 with: kwant.plotter.bands(sys.leads[0]) to see at what energies the lead supports propagating modes. Happy Kwanting, Joe