> Since the site.pos of honeycomb(electron)-sublattice-a is the same as the > site.pos of honeycomb(hole)-sublattice-c, should this method not also work > for this lattice?
The problem with this is that lat(a, b) isn't a site with coordinates a & b, but rather a site that belongs to the unit cell (a, b). That you could fix in your code by using site.tag instead of site.pos. But even then the orbital approach is better. Best, Anton On Thu, 14 Jan 2021 at 10:36, <henry....@physik.rwth-aachen.de> wrote: > > Hey Anton, > > Thank you for the reply. > > It is indeed a honeycomb lattice I am working with, and I attempting to > connecting the two sublattices the honeycomb lattices to create a > superconducting lead. > > Since the site.pos of honeycomb(electron)-sublattice-a is the same as the > site.pos of honeycomb(hole)-sublattice-c, should this method not also work > for this lattice? > > Maybe I will just rewrite it in the orbital method. > > Regards, > Henry