Dear Zhou,

The band structure can be divided into subgroups: the curves due to the translation in the y direction and the curves due to the number of sites in the unit cell.

When you do wraparound, you are fixing ky. So, the modes that remain are those due to the number of sites in the unit cell.

transmission.append(Sm.transmission(0,1)) sums over those modes since now ky is a parameter.

That is why it is a resolved transmission. I you want to get the total transmission, you need to sum over the all possible ky (integrating over the Brillouin zone).

I hope this helps,

Adel

On Mon, May 25, 2020 at 9:20 PM Zhou Jiaqi <jiaqi.zhou@uclouvain.be> wrote:

Dear Adel,

Thank you very much for your kind help. As you point out, wraparound can be the key to this problem. I learned a lot from your codes, and I am trying to convert it to my system.

BTW, could you please explain why you used 'trapz(transmission,Ky)' ? In my mind, 'transmission.append(Sm.transmission(0,1))' is enough for the transmission spectrum. So, what's the meaning of the integral using 'trapz' ? (In fact there is a code error for this instruction in my laptop even though I noticed the indentation)

Thanks again for your help!

Sincerely,
Jiaqi