Dear Adel, thank you for the reply. Actually, as you can see below, I have defined the bilayer lattice in 2D shape:
blg = kwant.lattice.general([(a0*cos_30, a0*sin_30), (a0*cos_30, -a0*sin_30)],
[(0, 0), (acc, 0), (acc, 0), (2 * acc, 0)],
norbs=1)
Just for your consideration I will send my code to you by your personal email.
Because, here I couldn't find any option to attach a file.
Thank you very much for your consideration!
Best,
Mohamad
