Dear Dongwook, It is the easiest to achieve what you want while you are defining your scattering region, specifically by using the Builder.degree method. So if you have a builder syst with all the hoppings already added, your desired result would be achieved by syst[(site for site in syst.sites() if syst.degree(site))) != 4)] = boundary_potential
Best, Anton On Tue, Jan 3, 2017 at 4:08 AM, Dongwook Go <[email protected]> wrote: > Dear developers, > > I’m interested in investigating the influence of the edge potential of the > sample. > For example, in an example by tutorial “quantum_wire.py”, I would like to put > on-site potential only at the edge atoms. > > One could determine whether the site is edge or not by reading the number of > neighbors of the site by using “degree”. > For example, if the degree is 3 there is an on-site potentials , else 0. > Is there any example code from which I can learn to use this function? > > > Best, > Dongwook.
