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@example.com wrote:
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