Hello Anton,
Thank you very much for your reply.
I am having a hard time with the "where" argument in
kwant.operator.Current. My system is not finalized. How can I find the
sites associated with the hopping? By their tag? Their position? Do
you have a piece of code for that?
Note that I am using TBModels and Wannier90 and the lattice has one
sublattice for each WF (or site, or orbital). norb is not defined, but
I see that
for x in lattice.sublattices:
x.norbs=1
works and looks correct.
Regards,
Eleni
Quoting Anton Akhmerov
Dear Eleni,
I think this should be doable with not too much of low level work. * Firstly you probably will want to select only the hoppings that are close to the 2D cut. You can do this by using "where" parameter to Current operator (see https://kwant-project.org/doc/1/reference/generated/kwant.operator.Current#k...) * After computing the current across those hoppings, you would need to call kwant.plotter.interpolate_current. This function only does the interpolation, but no plotting. Having that interpolation in 3D you need to slice it and probably take a projection onto 2D plane. * Finally, the resulting 2D array you can feed to kwant.plotter.streamplot
kwant.plotter.current is essentially a straightforward combination of kwant.plotter.interpolate_current and kwant.plotter.streamplot
Also: your problem sounds like a useful thing to do. Please let me know if you succeed, and what you ended up doing (or if you face any further problems as well).
Best, Anton
On Fri, Mar 23, 2018 at 5:22 PM,
wrote: Hello everyone,
Is there currently way to visualize the current through a 2D cut for a 3D system that does not involve messing with the kwant code?
I have a 3D system and would like to see the differences in the current when I switch some specific hoppings on and off.
Regards,
-- Dr. Eleni Chatzikyriakou Computational Physics lab Aristotle University of Thessaloniki elchatz@auth.gr - tel:+30 2310 998109
-- Dr. Eleni Chatzikyriakou Computational Physics lab Aristotle University of Thessaloniki elchatz@auth.gr - tel:+30 2310 998109