epraya.Lorbit#

Lorbit(sx, sy, sz, lamda, dim, l=0)#

Function for the spin orbit interaction, with an isotropic coupling constant. It’s recommended to use the Stevens Operators instead of this function

Parameters:

  • sx (np.array) – Pauli matrix of the spin x component.

  • sy (np.array) – Pauli matrix of the spin y component.

  • sz (np.array) – Pauli matrix of the spin z component.

  • lamda (float) – Isotropic coupling constant.

  • dim (int) – Dimension of the total hamiltonian.

  • l (float) – Angular momentum operator

Returns:

orbe (np.array) – Matrix of the spin orbit interaction with dimension dim.

Example

>>> import epraya as epr
>>> Ham,_,_=epr.Start()
>>> Ham.S=1/2
>>> Ham.L=1/2
>>> Ham.lc=20
>>> sx,sy,sz=epr.Pauli(Ham.S)
>>> dim=int((2*Ham.S+1)*(2*Ham.L+1))
>>> print(epr.Lorbit(sx,sy,sz,Ham.lc,dim,Ham.L))
([[ 5.+0.j,  0.+0.j,  0.+0.j,  0.+0.j],
[ 0.+0.j, -5.+0.j, 10.+0.j,  0.+0.j],
[ 0.+0.j, 10.+0.j, -5.+0.j,  0.+0.j],
[ 0.+0.j,  0.+0.j,  0.+0.j,  5.+0.j]])