The crystal field ground state of the Dy ion in DyNiBC can be described by a quasi-quartet consisting of two doublets separated by a energy interval .
In order to calculate efficiently the single ion property (for small effective magnetic fields in comparison to the total CF splitting) the Hamiltonian is projected to the quasi quartet and may be written as
with the angular momentum operators given by the 4x4 matrices
The constants - can be computed from the crystal field parameters, if these are known. On the other hand, they are connected to the saturation magnetic moments by the following equations ( denotes the maximum of the argument values)
The module given in /examples/dyni2b2c/1ion_mod/quartett.c diagonalises the Hamiltonian (131) and calculates the thermal expectation value of the vector ,which is returned to the McPhas program.Martin Rotter 2017-01-10