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
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
-- Abbout Adel