epraya.Mjhval#

class Mjhval#

Bases: object

Jax container class for the hamiltonial parameters of two systems. Must be initialized with a variable like Ham. For the nth-system, it can be change using Ham.SN=1/2.

Note: The Hpp is the same for the two systems.

Warning: The Stevens operators constants are only defined for the first system.

Parameters:

  • S (float) – Spin value ex. (1/2,0,3/2).

  • g (array_like or float) – g value of system, can be float, for isotropic case or array for anisotropic.

  • I (float) – Nuclear spin value

  • L (float) – Angular momentum

  • A (array_like or float) – Hyperfine constant, float for isotropic and array for anisotropic

  • Q (array_like or float) – Quadrupole nuclear interaction constant, float for isotropic and array for anisotropic

  • D (array_like) – Zero field interaction constants D and E, two value array [0,0]

  • Bk2 (array_like) – Stevens k=-/+2 constants

  • Bk4 (array_like) – Stevens k=-/+4 constants

  • Bk6 (array_like) – Stevens k=-/+6 constants

  • lc (float) – Spin-orbit interaction constant

  • Hpp (array_like) – Peak to peak distance for the voigtian function using [Hg,Hl], for gaussian and lorentzian distance

  • eta (float) – weight of the gaussian contribution to the voigtian function, from 0 to 1. If eta is 0, the function is lorentzian and if eta is 1, the function is gaussian.

  • weight (float) – Dummy variable by the moment

  • Nucl (str) – Isotope of the sample. Can be the quantum number and the element or only the element (‘55Mn’ or ‘Mn’)

  • A1_2 (array_like) – Hipefine interaction between the spin 1 and nuclear spin 2.

  • X1_2 (array_like) – Electron-Electron interaction between spins 1 and 2.

Example

>>> import epraya as epr
>>> Ham=Mjhval()
>>> Ham.S1=1
Ham.I1=1/2
>>> Ham.g1=[2.003,1.8,1.5]
>>> Ham.S2=2
>>> Ham.I2=1
>>> Ham.g2=[2.003,1.8,1.5]
>>> Ham.A1_2=[200,300,200]
>>> Ham.X1_2=[500,1000,500]
>>> print(Ham)
Mjhval(S1=1, S2=2, g1=[2.003, 1.8, 1.5], g2=[2.003, 1.8, 1.5],
I1=0.5, I2=1, L1=0.0, L2=0.0, A1=0.0, A2=0.0,
Q1=Array([0, 0, 0], dtype=int32), Q2=Array([0, 0, 0], dtype=int32),
D1=Array([0, 0], dtype=int32), D2=Array([0, 0], dtype=int32),
Bk2=[0, 0, 0, 0, 0], Bk4=[0, 0, 0, 0, 0, 0, 0, 0, 0],
Bk6=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], lc1=0.0, lc2=0.0,
A1_2=[200, 300, 200], A2_1=0.0, X1_2=[500, 1000, 500],
Hpp=Array([0, 1], dtype=int32), eta=0.5, weight=0.0)