Dear all,
First off, thanks for all the hard work in this. Very much appreciated!
I'm trying to write a system with a p-wave scattering region ( H =
(p^2+V-\mu)\sigma_z +1/2 {\alpha, p} \cdot \sigma) with normal metal leads
(\alpha = 0 in the leads). For simplicity, I'm going to talk about 1-D
systems (chains) here, although I'm working on 2D systems.
Because of the term linear in p, you can't take it directly to zero. For
example: Define the p-wave parameter to be zero in the hopping between
sites x=0 and x=1 (with lattice spacing=1). Then define the hopping between
sites x=1 and x=2 to to be =1. Then, this will force all of the hoppings to
undulate between 0 and 1. The way around this is to define it to be zero in
one hopping, then 0.5 in the next hopping, and then 1 in the hopping after
that.
My question is: Is this automatically implemented?
(Just in case, I'm putting this in by hand, not using HoppingKind, but
iterating over nearest neighbors and conditionally setting sys[lat(x,y),
lat(x+1,y) ] to values I want. )
Thank you very much,
Baris
