Dear Adel, Thank you for the link. I have another question regarding the translation symmetry. symlead=kwant.TranslationalSymmetry((-n*a,0)) lead=kwant.Builder(symlead) sys.attach(lead) sys.attach(lead.reversed())
If n>1, then, there is some extra addition of the layers of the lead to the scattering region (different numbers on right and left side of the scattering region) when I use kwant.plot(sys) and sys.finalized().sites(). If I use params for defining the onsite and hopping term for cases with n>1, will the kwant be able to recognize this layers or do I need to define by myself which layer is the main scattering regions and which one is the additional lead layers added due to the translation symmetry? Thank you Hari
