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.
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.
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.
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.
We provide a detailed exposition of our computational framework designed for the accurate calculation of real-frequency dynamical correlation functions of the single-impurity Anderson model in the regime of weak to intermediate coupling. Using quantum field theory within the Keldysh formalism to directly access the self-energy and dynamical susceptibilities in real frequencies, as detailed in our recent publication [Ge et al., Phys. Rev. B 109, 115128 (2024)], the primary computational challenge is the full three-dimensional real-frequency dependence of the four-point vertex. Our codebase provides a fully MPI+OpenMP parallelized implementation of the functional renormalization group (fRG) and the self-consistent parquet equations within the parquet approximation. It leverages vectorization to handle the additional complexity imposed by the Keldysh formalism, using optimized data structures and highly performant integration routines. Going beyond the results shown in the previous publication, the code includes functionality to perform fRG calculations in the multiloop framework, up to arbitrary loop order, including self-consistent self-energy iterations. Moreover, implementations of various regulators, such as hybridization, interaction, frequency, and temperature, are supplied.
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
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.
Using a leading algorithmic implementation of the functional renormalization group (fRG) for interacting fermions on two-dimensional lattices, we provide a detailed analysis of its quantitative reliability for the Hubbard model. In particular, we show that the recently introduced multiloop extension of the fRG flow equations for the self-energy and two-particle vertex allows for a precise match with the parquet approximation also for two-dimensional lattice problems. The refinement with respect to previous fRG-based computation schemes relies on an accurate treatment of the frequency and momentum dependences of the two-particle vertex, which combines a proper inclusion of the high-frequency asymptotics with the so-called truncated unity fRG for the momentum dependence. The adoption of the latter scheme requires, as an essential step, a consistent modification of the flow equation of the self-energy. We quantitatively compare our fRG results for the self-energy and momentum-dependent susceptibilities and the corresponding solution of the parquet approximation to determinant quantum Monte Carlo data, demonstrating that the fRG is remarkably accurate up to moderate interaction strengths. The presented methodological improvements illustrate how fRG flows can be brought to a quantitative level for two-dimensional problems, providing a solid basis for the application to more general systems.