Dear all, I’m new in kwant and I want to use this package to calculate transport properties of some graphene-liked materials with multi orbital per atom. In order to define hopping in the form of matrix, we need to determine the relative coordinate of neighbor sites. If we build the system using latt = kwant.lattice.honeycomb Could anyone help me to show how can I make a hole on the center of main region and leadS (I would like to make a graphene nanoribbon which have holes)?like in attached fiqure. Please find attached file. Best, Nafise
Hi,
I’m new in kwant and I want to use this package to calculate transport properties of some graphene-liked materials with multi orbital per atom. In order to define hopping in the form of matrix, we need to determine the relative coordinate of neighbor sites. If we build the system using latt = kwant.lattice.honeycomb
Could anyone help me to show how can I make a hole on the center of main region and leadS (I would like to make a graphene nanoribbon which have holes)?like in attached fiqure. Please find attached file.
Making a hole in a system is shown in the graphene tutorial [1]. Happy Kwanting, Joe [1]: https://kwant-project.org/doc/1/tutorial/graphene
Dear Joseph Weston,
Thanks a lot for your quick response.
I wrote a code, it is in the following:
import kwant
import math
import matplotlib
d=1.32;
a1=2.62;
t=-0.033;
latt = kwant.lattice.general([(a1,0),(a1*0.5,a1*math.sqrt(3)/2)],
[(a1/2,-d/2),(a1/2,d/2)])
a,b = latt.sublattices
syst= kwant.Builder()
def rectangle(pos):
x, y = pos
z=x**2+y**2
return -3.2*a1
sym = kwant.TranslationalSymmetry(latt.vec((-1,0))) sym.add_site_family(latt.sublattices[0], other_vectors=[(-1, 2)]) sym.add_site_family(latt.sublattices[1], other_vectors=[(-1, 2)]) lead = kwant.Builder(sym)
Making a hole on the scatter reign is ok. I did the same procedure for lead to make hole on it but at the end there isn't any hole on the leads (it is like a rectangular). That is very kind of you if let me know what is the problem? how can I make it?
You haven't made the translational symmetry in the leads big enough. Think about how big the unit cell of the leads is, and how big the hole you are trying to make is. Joe
"sym = kwant.TranslationalSymmetry(latt.vec((-1,0))) sym.add_site_family(latt.sublattices[0], other_vectors=[(-1, 2)]) sym.add_site_family(latt.sublattices[1], other_vectors=[(-1, 2)]) lead = kwant.Builder(sym)" Would you please explain it a bit more. I am almost new in kwant, Thanks in advance for your reply.
participants (2)
-
Joseph Weston -
nafise.nour＠sci.ui.ac.ir