Hi Ousmane et al,

I actually want the block of unit cells to be perpendicular to the scattering region as in the Wanted_System.png, attached. Note that it is a representative plot, I didn't produce it via Kwant but I want to attach the red lead to the black region in that direction, not like in Unwanted_System.png, which is also attached and produced via the code:

##########################################

import kwant
import numpy as np
import matplotlib.pyplot as plt
from types import SimpleNamespace

def make_system(a, W, L):

    def shape(pos):
          (x, y) = pos

          return (0 <= y <= W/2 and -x/2 <= y <= (-x + L)/2) # or (-W <= y <= 0 and 0 <= x <= L)

    def onsite(site, par):
        return 4 * par.t - par.mu

    def hopx(par):
        return -par.t

    def hopy(par):
        return -par.t


    # lead
    lat = kwant.lattice.square(a, norbs=1)
    syst = kwant.Builder()
    syst[lat.shape(shape, (0, 0))] = onsite
    syst[kwant.builder.HoppingKind((1, 0), lat, lat)] = hopx
    syst[kwant.builder.HoppingKind((0, 1), lat, lat)] = hopy
    lead = kwant.Builder(kwant.TranslationalSymmetry((-a, -2*a)))

    def lead_shape(pos):
        (x, y) = pos
        return -W <= x <= 0

    def lead_hopx(site1, site2, par):
        return -par.t

    def lead_onsite(site, par):
        return 4 * par.t - par.mu

    lead[lat.shape(lead_shape, (-2*a, -a))] = lead_onsite
    lead[kwant.builder.HoppingKind((1, 0), lat, lat)] = lead_hopx
    lead[kwant.builder.HoppingKind((0, 1), lat, lat)] = hopy
    syst.attach_lead(lead)



    syst = syst.finalized()
    return syst


syst = make_system(a=5, W=100, L=100)
kwant.plot(syst)

##########################################

I couldn't make the lead rotated in the perpendicular direction to the surface. Is there a way to do it?

Thanks,

Barış