Dear kwanters, I am using kwant for handle a k.p Hamiltonian, for the conductance calculations. I start from a k.p model, and transform it to a tight-binding model using kwant.continuum tools [procedures listed here https://kwant-project.org/doc/1/tutorial/discretize]. Now I want to study the conductance of the system under magnetic field. I try to do Peierls substitution for hoppings between site i and site j (i-j being parallel to the lead). However, I find that using the kwant.continuum scheme, there are no information regarding the lattice information such as "lat = kwant.lattice.honeycomb() ". This prevents the application of Peierls substitution by the following way:
_____________________________ def hopx(site1, site2, B): y = site1.pos[1] return - np.exp(-1j * B * y) syst[kwant.builder.HoppingKind((1, 0), lat.a, lat.b)] = hopx ____________________________ My question is how to add magnetic field during the calculations for the systems generated from a continuum model by kwant.continuum? Thank you very much!