Hi,

On 04/10/2018 01:50 PM, elch...@auth.gr wrote:
> Hello everyone,
>
> I'm not sure what I am doing wrong, but here is the output when using
> interpolate_current():
>
> --------------------------------------------------------------------
>>>> kwant_sys = wraparound.wraparound(kwant_model).finalized()
>>>> ham_mat = kwant_sys.hamiltonian_submatrix(sparse=False)
>>>> evecs = sla.eigsh(ham_mat, k=48, which='SM')[1]
>>>> J = kwant.operator.Current(kwant_sys)
>>>> current = J(evecs[:, 9])
>>>> IJ = kwant.plotter.interpolate_current(kwant_sys, J)

You're passing "J" to "interpolate current", when you should pass "current"

>
> and I think I will try also using ipyvolume or mayavi. But I would
> like to know if there is a problem with this being a closed system. 

No, no problem. Of course if you don't break time reversal symmetry
(e.g. with magnetic field) then you're eigenstates will have 0 current
everywhere.


Happy Kwanting,


Joe

Reply via email to