ham_so

Spin-orbit interaction Hamiltonian.

Syntax

H = ham_so(Sys)
H = ham_so(Sys,eSpins)
H = ham_so(Sys,eSpins,'sparse')

Description

This function computes the Hamiltonian of all spin-orbit coupling interactions (SOCs) in the spin system described by the structure Sys.

If the vector eSpins is given, the SOC of only the specified electron spins is returned (1 is the first electron spin, 2 the second, etc). Otherwise, all electrons are included.

If 'sparse' is given as third argument, the Hamiltonian matrix is returned in sparse format instead of full format.

Examples

The spin-orbit coupling Hamiltonian between an electron spin-1/2 with orbital angular momentum L=2 is

Sys.S = 1/2;
Sys.L = 1;
Sys.soc = 13e3;
ham_so(Sys)

ans =
   1.0e+03 *
    6.5000         0         0         0         0         0
         0         0         0    9.1924         0         0
         0         0   -6.5000         0    9.1924         0
         0    9.1924         0   -6.5000         0         0
         0         0    9.1924         0         0         0
         0         0         0         0         0    6.5000