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