Hi again, I f you want translational symmetry, you need to use wraparound module and you will need to find your result only on one unit cell. You will have also to do integration on the Brillouin zone.
I hope this helps Adel On Thu, May 7, 2020 at 5:12 PM Naveen Yadav <naveengunwa...@gmail.com> wrote: > Dear sir, > I understand what you have said. But how can I maintain translation > symmetry ?because my system is 3D and for plotting current I access the > sites of 3D system using *list(sys.sites) *and plot the current for 2D > slice. Please suggest me. > > Best Regards > Naveen Yadav > Research Scholar > Department of Physics & Astrophysics > University of Delhi > New Delhi-110007 > > On Thu, May 7, 2020, 19:06 Abbout Adel <abbout.a...@gmail.com> wrote: > >> Dear Naveen, >> >> What you get is what is expected. You do not have translational symmetry. >> Rewrite your code by keeping only two dimensions and you will see why. >> >> I hope this helps, >> Adel >> >> >> On Thu, May 7, 2020 at 12:52 PM Naveen Yadav <naveengunwa...@gmail.com> >> wrote: >> >>> Dear KWANT Developers, >>> I am trying to plot the current density. The procedure is >>> straightforward. I have attached leads to the scattering region (leads have >>> same onsite and hopping as of the scattering region) as >>> >>> >>> >>> >>> >>> >>> >>> >>> >>> >>> >>> *sys[(lat(z,y,x) for z in range(H) for y in range(W)for x in >>> range(L))]=onsitesys[kwant.builder.HoppingKind((1, 0, 0), lat, lat)] = >>> hoppingzsys[kwant.builder.HoppingKind((0, 1, 0), lat, lat)] = >>> hoppingysys[kwant.builder.HoppingKind((0, 0, 1), lat, lat)] = >>> hoppingxlead=kwant.Builder(kwant.TranslationalSymmetry((1,0,0)))lead[(lat(z,y,x) >>> for z in range(H) for y in range(W)for x in >>> range(L))]=onsitelead[kwant.builder.HoppingKind((1, 0, 0), lat, lat)] = >>> hoppingzlead[kwant.builder.HoppingKind((0, 1, 0), lat, lat)] = >>> hoppingylead[kwant.builder.HoppingKind((0, 0, 1), lat, lat)] = >>> hoppingxsys.attach_lead(lead, add_cells=80)sys.attach_lead(lead.reversed())* >>> [image: image.png] >>> As you can see from the plot there is a discontinuity in the plot at >>> (W=20). Why it is so even if the scattering region and leads have same >>> onsite and hopping and there is translation symmetry throughout ? >>> >>> -- >>> Best Regards, >>> Naveen Yadav >>> Research Scholar >>> Department of Physics & Astrophysics >>> University Of Delhi >>> New Delhi-110007 >>> >> >> >> -- >> Abbout Adel >> > -- Abbout Adel
