Hi,

> We consider the values of the  conservation law based on the
> off-diagonal term of Hamiltonian matrix, details are as follows:
>
>  
>
> H = tinyarray.array([[0, 0, 0, 1, 0, 0],
>
>                                 [0, 0, 0, 0, 1, 0],
>
>                                 [0, 0, 0, 0, 0, 1],
>
>                                 [1, 0, 0, 0, 0, 0],
>
>                                 [0, 1, 0, 0, 0, 0],
>
>                                 [0, 0, 1, 0, 0, 0]])
>

I don't really understand; is this the full Hamiltonian?

> spin_block = tinyarray.array([[0, 0, 0, 1, 0, 0],
>
>                                              [0, 0, 0, 0, 1, 0],
>
>                                              [0, 0, 0, 0, 0, 1],
>
>                                              [-1, 0, 0, 0, 0, 0],
>
>                                              [0, -1, 0, 0, 0, 0],
>
>                                              [0, 0, -1, 0, 0, 0]])
>
> And we add “conservation_law= - spin_block” in  “lead=kwant.Builder()”
>
>  
>
> The program can not work, and prompt error:
>
>  
>
> IndexError: index 2 is out of bounds for axis 0 with size 2
>

There is not enough information for me to be able to help you. Please
attach a complete Kwant script.


Happy Kwanting,


Joe

Attachment: signature.asc
Description: OpenPGP digital signature

Reply via email to