Renormalization group methods are well-established tools for the (numerical) investigation of the low-energy properties of correlated quantum many-body systems, allowing us to capture their scale-dependent nature. The functional renormalization group (FRG) allows us to continuously evolve a microscopic model action to an effective low-energy action as a function of decreasing energy scales via an exact functional flow equation, which is then approximated by some truncation scheme to facilitate computation. Here, we report on our implementation of multiloop FRG, an extended truncation scheme recently developed for electronic FRG calculations, within the pseudofermion functional renormalization group (pf-FRG) framework for interacting quantum spin systems. We discuss in detail the conceptual intricacies of the flow equations generated by the multiloop truncation, as well as essential refinements to the integration scheme for the resulting integrodifferential equations. To benchmark our approach, we analyze antiferromagnetic Heisenberg models on the pyrochlore, simple cubic, and face-centered cubic lattice, discussing the convergence of physical observables for higher-loop calculations and comparing with existing results where available. Combined, these methodological refinements systematically improve the pf-FRG approach to one of the numerical tools of choice when exploring frustrated quantum magnetism in higher spatial dimensions.
We derive functional renormalization group (fRG) schemes for Fermi systems which are based on the two-particle irreducible approach to the quantum many-body problem. In a first step, the cutoff is introduced in the non-interacting propagator as it is commonly done in fRG based on one-particle irreducible vertex functions. The most natural truncation of the resulting infinite hierarchy of flow equations is shown to be fully equivalent to self-consistent perturbation theory. An earlier suggested alternative truncation strategy is considered as well. In a second step, the cutoff is introduced in the two-particle interaction. Again two truncation procedures are investigated, one of which was derived before. In the latter, the mean-field solution of the many-body problem is considered as the starting point of the renormalization group flow. We compare the performance and the required numerical resources for solving the coupled flow equations for all the approximate schemes by applying them to the problem of the quantum anharmonic oscillator. In a functional integral representation, this model has a formal similarity to the quantum many-body problem. The perspectives for applying the derived two-particle irreducible fRG approaches to zero- and one-dimensional systems of correlated fermions are discussed.
The spin-boson model has nontrivial quantum phase transitions in the sub-Ohmic regime. For the bath spectra exponent $0\ensuremath{\leqslant}s<1/2$, the bosonic numerical renormalization group (BNRG) study of the exponents $\ensuremath{\beta}$ and $\ensuremath{\delta}$ are hampered by the boson-state truncation, which leads to artificial interacting exponents instead of the correct Gaussian ones. In this paper, guided by a mean-field calculation, we study the order-parameter function $m(\ensuremath{\tau}=\ensuremath{\alpha}\ensuremath{-}{\ensuremath{\alpha}}_{c},\ensuremath{\epsilon},\ensuremath{\Delta})$ using BNRG. Scaling analysis with respect to the boson-state truncation
Nb, the logarithmic discretization parameter $\ensuremath{\Lambda}$, and the tunneling strength $\ensuremath{\Delta}$ are carried out. Truncation-induced multiple-power behaviors are observed close to the critical point, with artificial values of $\ensuremath{\beta}$ and $\ensuremath{\delta}$. They cross over to classical behaviors with exponents $\ensuremath{\beta}=1/2$ and $\ensuremath{\delta}=3$ on the intermediate scales of $\ensuremath{\tau}$ and $\ensuremath{\epsilon}$, respectively. We also find $\ensuremath{\tau}/{\ensuremath{\Delta}}^{1\ensuremath{-}s}$ and $\ensuremath{\epsilon}/\ensuremath{\Delta}$ scalings in the function $m(\ensuremath{\tau},\ensuremath{\epsilon},\ensuremath{\Delta})$. The role of boson-state truncation as a scaling variable in the BNRG result for $0\ensuremath{\leqslant}s<1/2$ is identified and its interplay with the logarithmic discretization revealed. Relevance to the validity of quantum-to-classical mapping in other impurity models is discussed.