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 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

integer and the order determines the order of blocks. Is the exact value of the eigenvalues of any importance?

The only meaning of each eigenvalue is the number of the block of the conserved quantity. If memory serves me right, we adopted a convention where they have to be consecutive integers. Cheers, Anton

Best, Daniel

Am 22.06.21 um 16:52 schrieb Anton Akhmerov:

...

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 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