Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1807.02697 · SciPost Physics

Multiloop functional renormalization group for the two-dimensional Hubbard model: Loop convergence of the response functions

Agnese Tagliavini, Cornelia Hille, Fabian B. Kugler, Sabine Andergassen, A. Toschi, Carsten Honerkamp A

Published Jan 18, 2019• 74 citations
arXivHTMLPDF

Abstract

We present a functional renormalization group (fRG) study of the two dimensional Hubbard model, performed with an algorithmic implementation which lifts some of the common approximations made in fRG calculations. In particular, in our fRG flow; (i) we take explicitly into account the momentum and the frequency dependence of the vertex functions; (ii) we include the feedback effect of the self-energy; (iii) we implement the recently introduced multiloop extension which allows us to sum up all the diagrams of the parquet approximation with their exact weight. Due to its iterative structure based on successive one-loop computations, the loop convergence of the fRG results can be obtained with an affordable numerical effort. In particular, focusing on the analysis of the physical response functions, we show that the results become independent from the chosen cutoff scheme and from the way the fRG susceptibilities are computed, i.e., either through flowing couplings to external fields, or through a “post-processing” contraction of the interaction vertex at the end of the flow. The presented substantial refinement of fRG-based computation schemes paves a promising route towards future quantitative fRG analyses of more challenging systems and/or parameter regimes.

DOI
10.21468/scipostphys.6.1.009

Related papers

Quantitative functional renormalization group description of the two-dimensional Hubbard model

Cornelia Hille, Fabian B. Kugler, Christian J. Eckhardt +5 more

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 twodimensional 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 momentumdependent 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.

arXiv:2002.02733•Physical Review Research•Sep 8, 2020•61 citations

Multiloop functional renormalization group from single bosons

Fraboulet, Kilian, Al-Eryani, Aiman, Heinzelmann, Sarah +2 more

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.

arXiv:2512.11190•Dec 12, 2025•0 citations

Quantitative functional renormalization group description of the two-dimensional Hubbard model

Cornelia Hille, Fabian B. Kugler, Christian Eckhardt +5 more

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.

•Jan 1, 2020•55 citations

Single-boson exchange functional renormalization group application to the two-dimensional Hubbard model at weak coupling

Kilian Fraboulet, Sarah Heinzelmann, Pietro M. Bonetti +4 more

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.

arXiv:2206.14478•The European Physical Journal B•Dec 1, 2022•16 citations

Multiloop functional renormalization group for general models

Fabian B. Kugler, Jan von Delft

We present multiloop flow equations in the functional renormalization group (fRG) framework for the four-point vertex and self-energy, formulated for a general fermionic many-body problem. This generalizes the previously introduced vertex flow [F. B. Kugler and J. von Delft, Phys. Rev. Lett. 120, 057403 (2018)] and provides the necessary corrections to the self-energy flow in order to complete the derivative of all diagrams involved in the truncated fRG flow. Due to its iterative one-loop structure, the multiloop flow is well suited for numerical algorithms, enabling improvement of many fRG computations. We demonstrate its equivalence to a solution of the (first-order) parquet equations in conjunction with the Schwinger-Dyson equation for the self-energy.

arXiv:1707.04536•Physical review. B./Physical review. B•Jan 31, 2018•88 citations

Two-loop functional renormalization group approach to the one- and two-dimensional Hubbard model

A. A. Katanin

We consider the application of the two-loop functional renormalization group (fRG) approach to study the low-dimensional Hubbard model. This approach accounts for both the universal and nonuniversal contributions to the RG flow. While the universal contributions were studied previously within the field-theoretical RG for the one-dimensional Hubbard model with linearized electronic dispersion and the two-dimensional Hubbard model with flat Fermi surface, the nonuniversal contributions to the flow of the vertices and susceptibilities appear to be important at large momenta scales. The two-loop fRG approach is also applied to the two-dimensional Hubbard model with a curved Fermi surface and the van Hove singularities near the Fermi level. The vertices and susceptibilities in the end of the flow of the two-loop approach are suppressed in comparison with both the one-loop approach with vertex projection and the modified one-loop approach with corrected vertex projection errors. The results of the two-loop approach are closer to the results of one-loop approach with the projected vertices, the difference of the results of one- and two-loop fRG approaches decreases when going away from the van Hove band filling. The quasiparticle weight remains finite in two dimensions at not too low temperatures above the paramagnetic ground state.

arXiv:0806.0300•Physical Review B•Jun 11, 2009•36 citations
We present a functional renormalization group (fRG) study of the two dimensional Hubbard model, performed with an algorithmic implementation which lifts some of the common approximations made in fRG calculations. In particular, in our fRG flow; (i) we take explicitly into account the momentum and the frequency dependence of the vertex functions; (ii) we include the feedback effect of the self-energy; (iii) we implement the recently introduced multiloop extension which allows us to sum up all the diagrams of the parquet approximation with their exact weight. Due to its iterative structure based on successive one-loop computations, the loop convergence of the fRG results can be obtained with an affordable numerical effort. In particular, focusing on the analysis of the physical response functions, we show that the results become independent from the chosen cutoff scheme and from the way the fRG susceptibilities are computed, i.e., either through flowing couplings to external fields, or through a “post-processing” contraction of the interaction vertex at the end of the flow. The presented substantial refinement of fRG-based computation schemes paves a promising route towards future quantitative fRG analyses of more challenging systems and/or parameter regimes.
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 twodimensional 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 momentumdependent 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.
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.
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.
We present multiloop flow equations in the functional renormalization group (fRG) framework for the four-point vertex and self-energy, formulated for a general fermionic many-body problem. This generalizes the previously introduced vertex flow [F. B. Kugler and J. von Delft, Phys. Rev. Lett. 120, 057403 (2018)] and provides the necessary corrections to the self-energy flow in order to complete the derivative of all diagrams involved in the truncated fRG flow. Due to its iterative one-loop structure, the multiloop flow is well suited for numerical algorithms, enabling improvement of many fRG computations. We demonstrate its equivalence to a solution of the (first-order) parquet equations in conjunction with the Schwinger-Dyson equation for the self-energy.
We consider the application of the two-loop functional renormalization group (fRG) approach to study the low-dimensional Hubbard model. This approach accounts for both the universal and nonuniversal contributions to the RG flow. While the universal contributions were studied previously within the field-theoretical RG for the one-dimensional Hubbard model with linearized electronic dispersion and the two-dimensional Hubbard model with flat Fermi surface, the nonuniversal contributions to the flow of the vertices and susceptibilities appear to be important at large momenta scales. The two-loop fRG approach is also applied to the two-dimensional Hubbard model with a curved Fermi surface and the van Hove singularities near the Fermi level. The vertices and susceptibilities in the end of the flow of the two-loop approach are suppressed in comparison with both the one-loop approach with vertex projection and the modified one-loop approach with corrected vertex projection errors. The results of the two-loop approach are closer to the results of one-loop approach with the projected vertices, the difference of the results of one- and two-loop fRG approaches decreases when going away from the van Hove band filling. The quasiparticle weight remains finite in two dimensions at not too low temperatures above the paramagnetic ground state.