Hi Henry, You can access the sublattices of the polyatomic lattice via `.sublattices` attribute. What kind of lattices are you dealing with? You seem to be using the site position (x, y) as a site tag, which will only work with a square lattice with lattice constant one.
Assuming you are indeed doing that, you can indeed use HoppingKind: lead2[HoppingKind((0, 0), c, a)] = lead2[HoppingKind((0, 0), d, b)] = Delta. However, I rather suggest identifying all degrees of freedom at the same position into a single site, and use small matrices for the onsite values. This is demonstrated for example in the tutorial: https://kwant-project.org/doc/1/tutorial/superconductors Best, Anton On Thu, 14 Jan 2021 at 09:51, <henry....@physik.rwth-aachen.de> wrote: > > Hello, > > I have set up a system in which I want to define hoppings in a lead between > sites that originate from two sublattices: LatticeA has sublatices a and b, > LatticeB has sublattices c and d. Currently, I believe I am achieving this > via: > > for site in lead2.sites(): > (x,y) = site.pos > if str(site.family) == "<Monatomic lattice e0>": > lead2[c(x,y),a(x,y)] = Delta > if str(site.family) == "<Monatomic lattice e1>": > lead2[d(x,y),b(x,y)] = Delta > > and thereby connecting sublattices c and a, and d and b. > > Is this accurate? Can this be done with HoppingKind? And is there a way for > me to list all the hoppings created to check if this was successful?
