Dear Pawel, you have to write: sym = kwant.TranslationalSymmetry((-dx,)) (and replace the lat(0, 0) by lat(0)) The extra comma is important! (-dx) is just a scalar, (-dx,) a tuple with one entry. TranslationalSymmetry needs a tuple/list/vector as input. Best, Michael On 21.05.2014 09:57, Pawel Wojcik wrote:
Dear all, I would like to simulate transport through a simple 1D structure, and i used chain() function to define the lattice but i dont know how to define leads. I tried to use transationalsymmetry() function but it does not work. Below i paste part of my simple code
def make_system(L=100,dx=1.0,m=0.0667):
lat=kwant.lattice.chain(dx) sys=kwant.Builder()
sys[(lat(i) for i in range(L))]=onsite sys[lat.neighbors()]=hopping
sym = kwant.TranslationalSymmetry((-dx,0)) #here i also tried to use kwant.TranslationalSymmetry((-dx)) l_lead=kwant.Builder(sym) l_lead[(lat(0,0))]=onsite l_lead[lat.neighbors()]=hopping sys.attach_lead(l_lead) sys.attach_lead(l_lead.reversed())
return sys.finalized()
I would be greatfull for advises. Thank you in advance Best Regards Pawel