Hi,

I am trying to calculate the conductance through coupled quantum nanowires. I have basically considered a square lattice with hopping tx in x direction and modulating the coupling between the wires through ty i.e. the hopping in y direction. My system requires  periodic boundary condition along x direction and open boundary condition in y direction. Can you please help me in  how this can be done using kwant to meet the above boundary conditions. I have mentioned my code below. Since I am  new to Kwant, it would be of much help to me if you help me solving this.

Regards
Deepti Rana

def make_sys1(a=1,W=4,barrierpos=1):
   
   
    def onsite_normal(site,p):
        return( (2 * (p.tx+p.ty) - p.mu) * pauli.s0sz+p.Ez*pauli.szs0))

    def onsite_sc(site,p):
        return (2 * (p.tx+p.ty ) * pauli.s0sz+p.Ez*pauli.szs0)+p.delta*pauli.s0sx)

    def onsite_barrier(site,p):
        return((2*(p.tx+p.ty)+ p.Vbarrier-p.mu)*pauli.s0sz+p.Ez*pauli.szs0))

    def hopx(site1, site2, p):
        return -p.tx * pauli.s0sz +1j* p.alphax * pauli.sysz

    def hopy(site1, site2, p):
        return -p.ty * pauli.s0sz -1j * p.alphay * pauli.sxsz
  
   
    lat = kwant.lattice.square(a,norbs=4)
    syst = kwant.Builder()
   
    syst[(lat(x, y) for x in range(barrierpos)
         for y in range(W))]=onsite_barrier
    syst[kwant.builder.HoppingKind((1, 0), lat, lat)] = hopx
    syst[kwant.builder.HoppingKind((0, 1), lat, lat)] = hopy
   
    sym_left = kwant.TranslationalSymmetry((-a, 0))
    lead0 = kwant.Builder(sym_left, conservation_law=-pauli.s0sz)
    lead0[(lat(0, j) for j in range(W))]=onsite_normal
    lead0[kwant.builder.HoppingKind((1, 0), lat, lat)] = hopx
    lead0[kwant.builder.HoppingKind((0, 1), lat, lat)] = hopy
   
    sym_right = kwant.TranslationalSymmetry((a, 0))
    lead1 = kwant.Builder(sym_right)
    lead1[(lat(0, j) for j in range(W))]=onsite_sc
    lead1[kwant.builder.HoppingKind((1, 0), lat, lat)] = hopx
    lead1[kwant.builder.HoppingKind((0, 1), lat, lat)] = hopy
    syst.attach_lead(lead0)
    syst.attach_lead(lead1)

    syst = syst.finalized()
    return syst