High-order operators

In systems with electron spin S>1/2, higher-order terms have to be included in the spin Hamiltonian in order to properly describe EPR spectra. The most common parameters are the second-order zero-field splitting parameters D and E and the fourth-order parameters a and F.

In addition to the terms of D, E, a and F, a commonly used complete set of high-order operators are the extended Stevens operators (ESO) [eqn] for [eqn]. The associated coefficients [eqn] are always real. The high-order part of the spin Hamiltonian written in terms of ESOs is

[eqn]

In Easyspin, matrix representations of the extended Stevens operators are provided by the function stev and are used by the simulation functions eigfields, resfields and pepper.

Relation to conventional parameters

The commonly used second-order zero-field splitting parameters D and E are related to the [eqn] coefficients by

[eqn]

The associated second-order Hamiltonian terms have the forms

[eqn]

The traditional cubic and axial fourth-order parameters a and F are related to the [eqn] coefficients.

[eqn]

The associated fourth-order Hamiltonian terms have the forms

[eqn]
[eqn]

(see Abragam/Bleaney, p. 142, 437; Pilbrow, p. 125, 129). The axes x, y, and z are the fourfold symmetry axes of the cubic part of the crystal field.

Polynomial expressions

The following table lists the most common Stevens operators in terms of polynomials in S+, S- and Sz. A more general list of extended Stevens operators (including odd k) can be found in Al'tshuler/Kozyrev. Abragram/Bleaney contains only a partial, but compatible list. Ryabov (1999) has devised a formula for computing Stevens operator polynomials for arbitrary 0<=k and -k<=q<=k.

[eqn]
References

Here are the most important references for extended Stevens operators

Note that there are typos and mistakes both in Al'tshuler/Kozyrev and in Abragam/Bleaney, as discussed by Rudowicz (2004).