I was wondering if there is some link which explains how to put hoppings for a 3D system. The general procedure: sys[lat(i, j, k), lat(i - 1, j, k)] = -t doesn't seem to work and gives an error: Polyatomic object not callable.

You seem to be using a lattice with more than one sublattice, check out the graphene tutorial (kwant-project.org/doc/1.0/tutorial/tutorial4.html). Otherwise 3D isn't different from 2D.

Best, Anton

Thank you and regards, Girish

On Sun, Sep 27, 2015 at 8:53 AM, Anton Akhmerov anton.akhmerov@gmail.com wrote:

...

Dear Girish,

See the example of plotting an 1D wave function over here:

http://nbviewer.ipython.org/github/topocm/topocm_content/blob/master/w1_topo...

Best, Anton

On Fri, Sep 25, 2015 at 9:16 PM, Girish Sharma girishs@g.clemson.edu wrote:

...

Hi,

I have attached a small code, trying to visualize Majorana wavefunctions at the two end points of my 1D wire. The kwant.plotter doesn't really help, so I'm trying basic pyplot. However, due to random ordering of sites, I do not see the wavefunction peaking at the ends (see the last plot). Is there a way out to do this ordering correctly in PYTHON?

Thank you for your time !

Best regards, Girish