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
