Hi Barış, you can for instance make the lead_onsite different from onsite and then print syst[lat(i,j)] for different regions to check. Regards, Ousmane
Le sam. 6 juil. 2019 à 22:43, Barış Özgüler firstname.lastname@example.org a écrit :
How can I make sure that there is no overlap between the lead and normal region? Lead_shape() returns -L <= x <= L and we want the lead to satisfy x <= 0. If we don't use color_site(site), then we would get the attached plot. The black region not only represents the normal region but also represents the lead. However, the lead may also cover the region of x > 0 as the function Lead_shape() suggests. Is there another way to check what part of the shape is scattering region, what part is lead?
-- A. Barış Özgüler
PhD Candidate in Physics
University of Wisconsin–Madison Madison, Wisconsin 53706, USA