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 accurately describe EPR spectra. The most common parameters are the second-order zero-field splitting parameters D and E.

In addition to the terms of D and E, a commonly used complete set of high-order operators are the extended Stevens operators (ESO) [eqn] for [eqn] and [eqn].

The high-order part of the spin Hamiltonian for one electron spin written in terms of ESOs is

[eqn]

where each of the operators [eqn] is hermitian, and the associated coefficients [eqn] are always real-valued. Orders with k>6 are practically irrelevant.

In contrast to spherical tensor operators (analogous to p+, p0, and p- orbitals), these operators are tesseral (analogous to px, py and pz). Positive k corresponds to cosine tesseral components, and negative k corresponds to sine tesseral components.

Be careful when using literature values for the coefficients, since conventions for the normalization and the phase of these operators vary widely. Also, many papers contain a significant amount of typos.

In EasySpin, matrix representations of the extended Stevens operators are provided by the function stev and can be used in most simulation functions.

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 via

[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 all extended Stevens operators for k=2, 4 and 6 in terms of polynomials in S+, S- and Sz. A more general list of extended Stevens operators (including odd k) can be found in Altshuler/Kozyrev. Abragam/Bleaney contains only a partial, but compatible list. Ryabov (1999) gives a formula for computing Stevens operator polynomials for arbitrary k and q.

[eqn]
References

Here are the some useful references for extended Stevens operators:

Note that there are typos in Altshuler/Kozyrev and in Abragam/Bleaney, as discussed by Rudowicz (2004).