Stevens operator class of exactly solvable one-dimensional spin models

The new class of exactly solvable models of one-dimensional spin chains is found. Interactions of spins with crystalline electric fields (CEF) and spin-spin coupling between neighboring sites are taken into account for those models. Spin-spin coupling includes higher-order spin operators (the Stevens operators) and describes the interaction between zth components of the spins. We show that essential parts of the Hamiltonians of these models can be written in terms of exactly solvable transverse field Ising model (TFIM), in which certain combinations of the parameters of spin-spin interactions play the role of the TFIM exchange constant, and components of the CEFs play the role of external transverse magnetic field of the TFIM. The obtained results are valid for any values of spin (or total moments for rare-earth ions) and for the wide class of CEFs. The requirement for the exact solvability is that the space of single spin states should be split by the CEF into the subspaces with dimensions 2 (i.e., the only singlets and doublets are allowed). Cases of integer and half-integer spins are considered. For finite-size chains with integer spins and open ends we show numerically that energy spectra always contain the states whose energies do not depend on the CEF (the magnetic field for the TFIM). The energies (per spin) of these states decay as L-1, where L is the chain length. This behavior qualitatively distinguishes these states from known Majorana edge states and strong zero modes, where the energy decreases exponentially in L.