Adel Belayadi wrote:
Jerry wrote :
Then I would like to attach leads along the new X-axis (just normally connect to the rectangular sample that doesn't need adding extra sites). I put the translation symmetry rotated to X and confine the lead shape in Y. I expected the lead would just connect to the sample with the interface in Y direction. However, many extra sites are added and the interface is rendered in the old y direction, although the lead does extend in X.
I guess it might work if you add site family to your symmetry. your_sym.add_site_family(your_latice.sublattices[0], other_vectors=[v]), where v is a tuple and in your case it might be like v=(-1, 2)
Thanks, Adel. The following thread offers background information: https://mail.python.org/archives/list/kwant-discuss@python.org/thread/M4PUQD... Also note that the unwanted sites do not change the physics of the system since the Hamiltonian of the complete system (central region + leads) is the same no matter what the inclination of the lead unit cells is. Still, I agree that it would be better if attach_lead would avoid adding sites to the central region whenever this is possible. If someone feels motivated to work on this, I’d be happy to lend a hand. Here is an old and unfinished (but if I recall correctly mostly working) attempt to solve the problem: https://gitlab.kwant-project.org/cwg/kwant/-/commits/adrien_no_add Christoph