Deriving quantum spin model for a zigzag-chain ytterbium magnet with anisotropic exchange interactions

Abstract: We derive a quantum spin Hamiltonian of the spin-1/2 zigzag chain realized in a rare earth ytterbium-based magnetic insulator, YbCuS2. This material undergoes a transition at 0.95K to an incommensurate magnetic phase with small moments, which does not conform to the nonmagnetic singlet ground state of the spin-1/2 Heisenberg model. We take account of octahedral crystal field effect, atomic spin-orbit coupling, and strong Coulomb interactions on Yb ions, and perform a four-order perturbation theory to evaluate the superexchange coupling constants. A small but finite anisotropic exchange coupling called {\Gamma}-term appears similarly to the case reported previously in other triangular-based magnets. By varying several material parameters, we figure out two important factors to enhance {\Gamma}-term. One is the splitting of excited f-states with two holes, which efficiently selects the perturbation terms associated with the lowest excited state having large total angular momentum, and accordingly with high spatial anisotropy. The other is the tilting of octahedra or distortion of S-Yb-S bond angle, which break the symmetry of the Slater-Koster overlap. These effects are systematically analyzed by the exchange anisotropy in units of pairs of octahedra within the local spin frame, which significantly reduces the complexity of directly referring to the Hamiltonian discussed in the global axis.