Hi Daniel, Khani, If you want to use multiple degrees of freedom in the leads, specify a conservation law with the matrix that has different eigenvalues for each combination of the quantities of interest. For example, if your sites have 4 orbitals, like Khani describes, you can specify syst = kwant.Builder(conservation_law=np.identity(np.arange(4)), ...).
Hope this helps, Anton On Tue, 22 Jun 2021 at 15:45, Daniel Hauck <[email protected]> wrote: > > Hi, > > I know this is a very old discussion, but I'm in a similar position and > want to look at a superconducting system with spin. > > Have you been able to resolve this? > > Best Regards, > Daniel
