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:2202.00226 · The European Physical Journal B

The plain and simple parquet approximation: single-and multi-boson exchange in the two-dimensional Hubbard model

Friedrich Krien, Anna Kauch

Published Apr 1, 2022• 22 citations
arXivHTMLPDF

Abstract

Abstract The parquet approach to vertex corrections is unbiased but computationally demanding. Most applications are therefore restricted to small cluster sizes or rely on various simplifying approximations. We have recently shown that the bosonization of the parquet diagrams provides interpretative and algorithmic advantages over the original purely fermionic formulation. Here, we present first results of the numerical implementation of this method by applying it to the half-filled Hubbard model on the square lattice at weak coupling. The improved algorithmic performance allows us to evaluate the parquet approximation for a $$16\times 16$$ 16 × 16 lattice, retaining the full momentum and frequency structure of the various vertex functions. We discuss their symmetries and consider parametrizations of their momentum dependence using the truncated-unity approximation. Graphical abstract

DOI
10.1140/epjb/s10051-022-00329-6

Related papers

Truncated-unity parquet equations: Application to the repulsive Hubbard model

Christian Eckhardt, G. A. H. Schober, Jannis Ehrlich +1 more

The parquet equations are a self-consistent set of equations for the effective two-particle vertex of an interacting many-fermion system. The application of these equations to bulk models is, however, demanding due to the complex emergent momentum and frequency structure of the vertex. Here, we show how a channel decomposition by means of truncated unities, which was developed in the context of the functional renormalization group to efficiently treat the momentum dependence, can be transferred to the parquet equations. This leads to a significantly reduced complexity and memory consumption scaling only linearly with the number of discrete momenta. We apply this technique to the half-filled repulsive Hubbard model on the square lattice and present approximate solutions for the channel-projected vertices and the full reducible vertex.

arXiv:1802.09797•Physical review. B./Physical review. B•Aug 24, 2018•26 citations

Boson-exchange parquet solver for dual fermions

Friedrich Krien, Angelo Valli, Patrick Chalupa-Gantner +3 more

We present and implement a parquet approximation within the dual-fermion formalism based on a partial bosonization of the dual vertex function which substantially reduces the computational cost of the calculation. The method relies on splitting the vertex exactly into single-boson exchange contributions and a residual four-fermion vertex, which physically embody, respectively, long- and short-range spatial correlations. After recasting the parquet equations in terms of the residual vertex, these are solved using the truncated-unity method of Eckhardt et al. [Phys. Rev. B 101, 155104 (2020)], which allows for a rapid convergence with the number of form factors in different regimes. While our numerical treatment of the parquet equations can be restricted to only a few Matsubara frequencies, reminiscent of Astretsov et al. [Phys. Rev. B 101, 075109 (2020)], the one- and two-particle spectral information is fully retained. In applications to the two-dimensional Hubbard model the method agrees quantitatively with a stochastic summation of diagrams over a wide range of parameters.

arXiv:2008.04184•Physical review. B./Physical review. B•Nov 18, 2020•51 citations

Truncated unity parquet solver

Christian Eckhardt, Carsten Honerkamp A, Karsten Held +1 more

We present an implementation of a truncated unity parquet solver (TUPS) which solves the parquet equations using a truncated form-factor basis for the fermionic momenta. This way fluctuations from different scattering channels are treated on an equal footing. The essentially linear scaling of computational costs in the number of untruncated bosonic momenta allows us to treat system sizes of up to $76\ifmmode\times\else\texttimes\fi{}76$ discrete lattice momenta, unprecedented by previous unbiased methods that include the frequency dependence of the vertex. With TUPS, we provide the first numerical evidence that the parquet approximation might indeed respect the Mermin-Wagner theorem and further systematically analyze the convergence with respect to the number of form factors. Using a single form factor seems to qualitatively describe the physics of the half-filled Hubbard model correctly, including the pseudogap behavior. Quantitatively, using a single or a few form factors only is not sufficient at lower temperatures or stronger coupling.

arXiv:1912.07469•Physical review. B./Physical review. B•Apr 6, 2020•64 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

Dynamical vertex approximation in its parquet implementation: Application to Hubbard nanorings

Angelo Valli, Thomas Schäfer, Patrik Thunström +5 more

We have implemented the dynamical vertex approximation $(\mathrm{D}\ensuremath{\Gamma}\mathrm{A})$ in its full parquet-based version to include spatial correlations on all length scales and in all scattering channels. The algorithm is applied to study the electronic self-energies and the spectral properties of finite-size one-dimensional Hubbard models with periodic boundary conditions (nanoscopic Hubbard rings). From a methodological point of view, our calculations and their comparison to the results obtained within dynamical mean-field theory, plain parquet approximation, and the exact numerical solution allow us to evaluate the performance of the $\mathrm{D}\ensuremath{\Gamma}\mathrm{A}$ algorithm in the most challenging situation of low dimensions. From a physical perspective, our results unveil how nonlocal correlations affect the spectral properties of nanoscopic systems of various sizes in different regimes of interaction strength.

arXiv:1410.4733•Physical Review B•Mar 6, 2015•99 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
Abstract The parquet approach to vertex corrections is unbiased but computationally demanding. Most applications are therefore restricted to small cluster sizes or rely on various simplifying approximations. We have recently shown that the bosonization of the parquet diagrams provides interpretative and algorithmic advantages over the original purely fermionic formulation. Here, we present first results of the numerical implementation of this method by applying it to the half-filled Hubbard model on the square lattice at weak coupling. The improved algorithmic performance allows us to evaluate the parquet approximation for a 16×1616\times 1616×16 16 × 16 lattice, retaining the full momentum and frequency structure of the various vertex functions. We discuss their symmetries and consider parametrizations of their momentum dependence using the truncated-unity approximation. Graphical abstract
The parquet equations are a self-consistent set of equations for the effective two-particle vertex of an interacting many-fermion system. The application of these equations to bulk models is, however, demanding due to the complex emergent momentum and frequency structure of the vertex. Here, we show how a channel decomposition by means of truncated unities, which was developed in the context of the functional renormalization group to efficiently treat the momentum dependence, can be transferred to the parquet equations. This leads to a significantly reduced complexity and memory consumption scaling only linearly with the number of discrete momenta. We apply this technique to the half-filled repulsive Hubbard model on the square lattice and present approximate solutions for the channel-projected vertices and the full reducible vertex.
We present and implement a parquet approximation within the dual-fermion formalism based on a partial bosonization of the dual vertex function which substantially reduces the computational cost of the calculation. The method relies on splitting the vertex exactly into single-boson exchange contributions and a residual four-fermion vertex, which physically embody, respectively, long- and short-range spatial correlations. After recasting the parquet equations in terms of the residual vertex, these are solved using the truncated-unity method of Eckhardt et al. [Phys. Rev. B 101, 155104 (2020)], which allows for a rapid convergence with the number of form factors in different regimes. While our numerical treatment of the parquet equations can be restricted to only a few Matsubara frequencies, reminiscent of Astretsov et al. [Phys. Rev. B 101, 075109 (2020)], the one- and two-particle spectral information is fully retained. In applications to the two-dimensional Hubbard model the method agrees quantitatively with a stochastic summation of diagrams over a wide range of parameters.
We present an implementation of a truncated unity parquet solver (TUPS) which solves the parquet equations using a truncated form-factor basis for the fermionic momenta. This way fluctuations from different scattering channels are treated on an equal footing. The essentially linear scaling of computational costs in the number of untruncated bosonic momenta allows us to treat system sizes of up to $76\ifmmode\times\else\texttimes\fi{}76$ discrete lattice momenta, unprecedented by previous unbiased methods that include the frequency dependence of the vertex. With TUPS, we provide the first numerical evidence that the parquet approximation might indeed respect the Mermin-Wagner theorem and further systematically analyze the convergence with respect to the number of form factors. Using a single form factor seems to qualitatively describe the physics of the half-filled Hubbard model correctly, including the pseudogap behavior. Quantitatively, using a single or a few form factors only is not sufficient at lower temperatures or stronger coupling.
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 have implemented the dynamical vertex approximation $(\mathrm{D}\ensuremath{\Gamma}\mathrm{A})$ in its full parquet-based version to include spatial correlations on all length scales and in all scattering channels. The algorithm is applied to study the electronic self-energies and the spectral properties of finite-size one-dimensional Hubbard models with periodic boundary conditions (nanoscopic Hubbard rings). From a methodological point of view, our calculations and their comparison to the results obtained within dynamical mean-field theory, plain parquet approximation, and the exact numerical solution allow us to evaluate the performance of the $\mathrm{D}\ensuremath{\Gamma}\mathrm{A}$ algorithm in the most challenging situation of low dimensions. From a physical perspective, our results unveil how nonlocal correlations affect the spectral properties of nanoscopic systems of various sizes in different regimes of interaction strength.
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.