Hi,

this is a follow up question to a previous (hijacked?) thread. For completeness I will provide all details, so that no cross referencing is necessary.

I want to compute the electrical conductance of a N-S interface embedded on a hexagonal lattice. The final result should look something like the picture I have attached. The BdG equation is specified in the .pdf.

At the moment I have both particle-hole and spin-up/down degrees of freedom, because later I want to add magnetic fields and such. To attack the problem I have tried to follow the steps outlined in tutorial 2.6 with the addition that I use 4 orbitals and a hexagonal lattice. The conductance should still be given as N - R_{ee,\uparrow} -  R_{ee,\downarrow} + R_{he,\uparrow} + R_{he,\downarrow}. Here e, h, \uparrow, and \downarrow refers to electron, hole, spin-up, and spin-down respectively.

In the documentation it says that smatrix.transmission((i, a), (j, b)) gives the transmission from block b of lead j to block a of lead i. I would therefore expect  smatrix.transmission((0, 2), (0, 0)) to give me the transmission of an incident electron of spin up in lead 0 to a reflected hole of spin up in lead 0. With this in mind I thought I could write

data.append(smatrix.submatrix((0, 0), (0, 0)).shape-

smatrix.transmission((0, 0), (0, 0)) -
smatrix.transmission((0, 1), (0, 0)) +
smatrix.transmission((0, 2), (0, 0)) +
smatrix.transmission((0,3), (0, 0)))
to compute the conductance. However, doing this gives me an out of bounds error.

As a last resort I tried to write
data.append(smatrix.submatrix((0, 0), (0, 0)).shape -
smatrix.transmission((0, 0), (0, 0)) +
smatrix.transmission((0, 1), (0, 0)))
as given in the tutorial 2.6. However, this gives me that the conductance is zero.

My question is then, how do I correctly interpret the smatrix.transmission function in my example so that I can compute the correct conductance?

I have attached my code in the .txt file.

Best,
Martin