Dear Kwant users, I write for advices concening apparent problems that I encountered programming with Kwant. I would like to calculate quantities related to transport (e.g. conductivity, Fano factor,.....) in 2D and 3D lattice systems hosting Dirac / Weyl points. I have tried to implement a couple of examples of semimetals with Kwant, where the lattices (scattering regions) are taken square or cubic, and the leads attached on entire faces of the lattices. I assumed both open and periodic (in the directions orthogonal to the leads) boundaries conditions for the lattices and for the leads. Generally, I have found strange things, like sudden falls to zero of the conductivities within the energies of the scattering areas, or suspicious conductivity asymmetries around the energies of the cones, even if the lattice spectra are symmetric. I am reasonably sure that the Hamiltonians are written correctly in Kwant, at least because I checked separately the spectra of the scattering zones and of the leads, comparing them with the analytics. Discussing at a conference, I got aware vaguely about possible problems that may arise when one wants to discretize in real space an Hamiltonian originally written in the space of momenta (also continuous). This is exaclty what I have done always; notice that I performed the discretization by hand, without exploiting the routine in Kwant. Do you know about similar issues with semimetals ? Do you have some advices ? Thank you very much to everybody and best regards, L. L.