Bold diagrammatic Monte Carlo technique for frustrated spin systems

Using fermionic representation of spin degrees of freedom within the Popov-Fedotov approach, we develop an algorithm for Monte Carlo sampling of skeleton Feynman diagrams for Heisenberg-type models. Our scheme works without modifications for any dimension of space, lattice geometry, and interaction range, i.e., it is suitable for dealing with frustrated magnetic systems at finite temperature. As a practical application, we compute uniform magnetic susceptibility of the antiferromagnetic Heisenberg model on the triangular lattice and compare our results with the best available high-temperature expansions. We also report results for the momentum dependence of the static magnetic susceptibility throughout the Brillouin zone.