The functional renormalization group (fRG) is acknowledged as a powerful tool in quantum many-body physics and beyond. On the technical side, conventional implementations of the fRG rely on regulators for bare propagators only. Starting from Schwinger-Dyson and Bethe-Salpeter equations, we develop here an fRG formulation where both bare propagators and bare interactions can be dressed with regulators. The approach thus obtained is an extension of the multiloop fRG recently introduced for many-fermion systems. Using the single-boson exchange decomposition, we show that the underlying flow equations are simply interpreted as adding a regulator to the bosonic propagator and that such an extension scarcely changes the original structure of the flow equations. Overall, we provide a framework for implementing approaches that cannot be realized with conventional fRG methods, such as temperature flows for models with retarded interactions. For concrete applications, we analyze the loop convergence of our formulation against a conventional cutoff scheme for the Anderson impurity model. Finally, we devise a temperature-flow scheme that implements a cutoff in both the propagator and the bare interaction, and demonstrate its validity on a model of an Anderson impurity coupled to a phonon.
The functional renormalization group (fRG) is acknowledged as a powerful tool in quantum many-body physics and beyond. On the technical side, conventional implementations of the fRG rely on regulators for bare propagators only. Starting from Schwinger-Dyson and Bethe-Salpeter equations, we develop here an fRG formulation where both bare propagators and bare interactions can be dressed with regulators. The approach thus obtained is an extension of the multiloop fRG recently introduced for many-fermion systems. Using the single-boson exchange decomposition, we show that the underlying flow equations are simply interpreted as adding a regulator to the bosonic propagator and that such an extension scarcely changes the original structure of the flow equations. Overall, we provide a framework for implementing approaches that cannot be realized with conventional fRG methods, such as temperature flows for models with retarded interactions. For concrete applications, we analyze the loop convergence of our formulation against a conventional cutoff scheme for the Anderson impurity model. Finally, we devise a temperature-flow scheme that implements a cutoff in both the propagator and the bare interaction, and demonstrate its validity on a model of an Anderson impurity coupled to a phonon.
The functional renormalization group (fRG) is an established tool in the treatment of correlated electron systems, notably for the description of competing instabilities. In recent years, methodological advancements led to the multiloop extension of the fRG, which systematically includes loop corrections beyond the conventional one-loop truncation and yields a quantitatively accurate description of two-dimensional lattice systems. At the same time, the single-boson exchange (SBE) decomposition of the two-particle vertex has been shown to offer both computational and interpretative advantages paving the way to more affordable approximation schemes. We here apply their combination coined as multiloop SBE fRG to the two-dimensional Hubbard model at weak coupling. After providing a detailed account of the underlying formalism in physical channels, we analyze the results for the frequency- and momentum-dependent vertex functions. We find that the SBE approximation, i.e., neglecting the flow of the multi-boson exchange contributions, accurately reproduces the parquet approximation at loop convergence. The presented algorithmic improvement opens the route for the treatment of more challenging parameter regimes and more realistic models.
We present a reformulation of the functional renormalization group (fRG) for many-electron systems, which relies on the recently introduced single-boson exchange (SBE) representation of the parquet equations [F. Krien, A. Valli, and M. Capone, Phys. Rev. B 100, 155149 (2019)]. The latter exploits a diagrammatic decomposition, which classifies the contributions to the full scattering amplitude in terms of their reducibility with respect to cutting one interaction line, naturally distinguishing the processes mediated by the exchange of a single boson in the different channels. We apply this idea to the fRG by splitting the one-loop fRG flow equations for the vertex function into SBE contributions and a residual four-point fermionic vertex. Similarly as in the case of parquet solvers, recasting the fRG algorithm in the SBE representation offers both computational and interpretative advantages: The SBE decomposition not only significantly reduces the numerical effort of treating the high-frequency asymptotics of the flowing vertices, but it also allows for a clear physical identification of the collective degrees of freedom at play. We illustrate the advantages of an SBE formulation of fRG-based schemes, by computing through the merger of dynamical mean-field theory and fRG the susceptibilities and the Yukawa couplings of the two-dimensional Hubbard model from weak to strong coupling, for which we also present an intuitive physical explanation of the results. The SBE formulation of the one-loop flow equations paves a promising route for future multiboson and multiloop extensions of fRG-based algorithms.
Abstract: We illustrate the algorithmic advantages of the recently introduced single-boson exchange (SBE) formulation for the one-loop functional renormalization group (fRG), by applying it to the two-dimensional Hubbard model on a square lattice. We present a detailed analysis of the fermion-boson Yukawa couplings and of the corresponding physical susceptibilities by studying their evolution with temperature and interaction strength, both at half filling and finite doping. The comparison with the conventional fermionic fRG decomposition shows that the rest functions of the SBE algorithm, which describe correlation effects beyond the SBE processes, play a negligible role in the weak-coupling regime above the pseudo-critical temperature, in contrast to the rest functions of the conventional fRG. Remarkably, they remain finite also at the pseudo-critical transition, whereas the corresponding rest functions of the conventional fRG implementation diverge. As a result, the SBE formulation of the fRG flow allows for a substantial reduction of the numerical effort in the treatment of the two-particle vertex function, paving a promising route for future multiboson and multiloop extensions.
Abstract The recently introduced single-boson exchange (SBE) decomposition of the four-point vertex of interacting fermionic many-body systems is a conceptually and computationally appealing parametrization of the vertex. It relies on the notion of reducibility of vertex diagrams with respect to the bare interaction U, instead of a classification based on two-particle reducibility within the widely used parquet decomposition. Here, we re-derive the SBE decomposition in a generalized framework (suitable for extensions to, e.g., inhomogeneous systems or real-frequency treatments) following from the parquet equations. We then derive multiloop functional renormalization group (mfRG) flow equations for the ingredients of this SBE decomposition, both in the parquet approximation, where the fully two-particle irreducible vertex is treated as an input, and in the more restrictive SBE approximation, where this role is taken by the fully U-irreducible vertex. Moreover, we give mfRG flow equations for the popular parametrization of the vertex in terms of asymptotic classes of the two-particle reducible vertices. Since the parquet and SBE decompositions are closely related, their mfRG flow equations are very similar in structure. Graphic abstract
We present a unified C++ implementation of the functional renormalization group and the parquet equations within the single-boson exchange formalism for several paradigmatic tight-binding and impurity models at equilibrium. The implementation computes the full dynamic vertex and self-energies, with momentum dependence treated using a truncated unity framework. We implement multiple self-energy flow equations, cutoff schemes, and extensions ranging from the dynamical functional renormalization group to multiloop flow equations that incorporate cutoffs in both the propagator and the interaction. The codebase serves as a reference for recent developments in the fRG and parquet methods for correlated electron systems and provides a flexible foundation for developing new many-body approaches and extensions.
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.