Peierls phase from kwant.physics.magnetic_gauge
Dear KWANT community, Recently I have used KWANT to simulate electronic transport in T-junction with the perpendicular magnetic field applied to the nanostructure. As well known, the inclusion of the orbital effects in this case (by the peierls substitution) is not easy due to appropriate gauge which should be applied for the vertical leads. So I decided to use the newest functionality of KWANT included in kwant.physics.magnetic_gauge. I started testing this function from a simple nanowire (NW) with two leads applied on both sides of NW and the magnetic field perpendicular to NW. When calling ham=sys.hamiltonian_submatrix(params=params), similarly as in the documentation, everything works OK. The problem (error) appears when I am trying to calculate transport and call the function smatrix = kwant.smatrix(sys, energy=energy, params=params). I don't know how is the reason of this error because everything seems to be correct. Enclosed please find my notebook.NW-orbital-gauge.ipynb
I would be grateful for any help or hint what is wrong. Best regards, Pawel Wojcik