epraya.JLorbit#

JLorbit(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 (jax.np.array) – Pauli matrix of the spin x component.

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

  • sz (jax.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 (jax.np.array) – Matrix of the spin orbit interaction with dimension dim.

Example

>>> import epraya as epr
>>> Ham, Exp, Vary=epr.Jstart()
>>> Ham.S=1/2
>>> Ham.L=1/2
>>> Ham.lc=20
>>> sx,sy,sz=epr.JPauli(Ham.S)
>>> dim=int((2*Ham.S+1)*(2*Ham.L+1))
>>> print(epr.JLorbit(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]]