The pseudofermion functional renormalization group (PFFRG) method has proven to be a powerful numerical approach to treat frustrated quantum spin systems. In its usual implementation, however, the complex fermionic representation of spin operators introduces unphysical Hilbert-space sectors which render an application at finite temperatures inaccurate. In this work we formulate a general functional renormalization group approach based on Majorana fermions to overcome these difficulties. We, particularly, implement spin operators via an
SO(3) symmetric Majorana representation which does not introduce any unphysical states and, hence, remains applicable to quantum spin models at finite temperatures. We apply this scheme, dubbed pseudo-Majorana functional renormalization group (PMFRG) method, to frustrated Heisenberg models on small spin clusters as well as square and triangular lattices. Computing the finite-temperature behavior of spin correlations and thermodynamic quantities such as free energy and heat capacity, we find good agreement with exact diagonalization and the high-temperature series expansion down to moderate temperatures. We observe a significantly enhanced accuracy of the PMFRG compared to the PFFRG at finite temperatures. More generally, we conclude that the development of functional renormalization group approaches with Majorana fermions considerably extends the scope of applicability of such methods.
We implement an extension of the pseudofermion functional renormalization group method for quantum spin systems that takes into account two-loop diagrammatic contributions. An efficient numerical treatment of the additional terms is achieved within a nested graph construction which recombines different one-loop interaction channels. In order to be fully self-consistent with respect to self-energy corrections, we also include certain three-loop terms of Katanin type. We first apply this formalism to the antiferromagnetic ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice and benchmark our results against the previous one-loop plus Katanin approach. Even though the renormalization group (RG) equations undergo significant modifications when including the two-loop terms, the magnetic phase diagram, comprising N\'eel ordered and collinear ordered phases separated by a magnetically disordered regime, remains remarkably unchanged. Only the boundary position between the disordered and the collinear phases is found to be moderately affected by two-loop terms. On the other hand, critical RG scales, which we associate with critical temperatures
Tc, are reduced by a factor of $\ensuremath{\sim}2$ indicating that the two-loop diagrams play a significant role in enforcing the Mermin-Wagner theorem. Improved estimates for critical temperatures are also obtained for the Heisenberg ferromagnet on the three-dimensional simple cubic lattice where errors in
Tc are reduced by $\ensuremath{\sim}34%$. These findings have important implications for the quantum phase diagrams calculated within the previous one-loop plus Katanin approach which turn out to be already well converged.
The pseudo-fermion functional renormalization group is generalized to treat spin Hamiltonians with finite magnetic fields, enabling its application to arbitrary spin lattice models with linear and bilinear terms in the spin operators. We discuss in detail an efficient numerical implementation of this approach making use of the system's symmetries. Particularly, we demonstrate that the inclusion of small symmetry breaking magnetic seed fields regularizes divergences of the susceptibility at magnetic phase transitions. This allows the investigation of spin models within magnetically ordered phases at
T=0 in the physical limit of vanishing renormalization group parameter
Λ. We explore the capabilities and limitations of this method extension for Heisenberg models on the square, honeycomb and triangular lattices. While the zero-field magnetizations of these systems are systematically overestimated, the types of magnetic orders are correctly captured, even if the local orientations of the seed field are chosen differently than the spin orientations of the realized magnetic order. Furthermore, the magnetization curve of the square lattice Heisenberg antiferromagnet shows good agreement with error controlled methods. In the future, the inclusion of magnetic fields in the pseudo-fermion functional renormalization group, which is also possible in three dimensional spin systems, will enable a variety of additional interesting applications such as the investigation of magnetization plateaus.