Dear all, I have written a code in Kwant for a bilayer graphene quantum dot (in 2D manner) to calculate some electronic properties. I can get the results correctly for, e.g, conductance, spectrum, and LDOS. But there is a problem in plotting CURRENT density and I get the following error.
~~ need at least one array to concatenate Has anyone had a similar experience ? The current density part of code is: ________________________________________________________________ ### Current plot params = dict(B=Bfield, r0=[0,0], h0=0) wf = kwant.wave_function(syst, energy=0.01, params=params) psi = wf(0)[0] J_0 = kwant.operator.Current(syst) current = J_0(psi, params=params) kwant.plotter.current(syst, current) _________________________________________________________________
