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:2008.04184 · Physical review. B./Physical review. B

Boson-exchange parquet solver for dual fermions

Friedrich Krien, Angelo Valli, Patrick Chalupa-Gantner, Massimo Capone, Alexander I. Lichtenstein, A. Toschi

Published Nov 18, 2020• 51 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.102.195131

Related papers

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

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

Friedrich Krien, Anna Kauch

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

arXiv:2202.00226•The European Physical Journal B•Apr 1, 2022•22 citations

Parquetlike equations for the Hedin three-leg vertex

Friedrich Krien, Angelo Valli

Taking the competition and the mutual screening of various bosonic fluctuations in correlated electron systems into account requires an unbiased approach to the many-body problem. One such approach is the self-consistent solution of the parquet equations, whose numerical treatment in lattice systems is however prohibitively expensive. In a recent article it was shown that there exists an alternative to the parquet decomposition of the four-point vertex function, which classifies the vertex diagrams according to the principle of single-boson exchange (SBE) [Krien, Valli, and Capone, Phys. Rev. B 100, 155149 (2019)]. Here we show that the SBE decomposition leads to a closed set of equations for the Hedin three-leg vertex, the polarization, and the electronic self-energy, which sums self-consistently the diagrams of the Maki-Thompson type. This circumvents the calculation of four-point vertex functions and the inversion of the Bethe-Salpeter equations, which are the two major bottlenecks of the parquet equations. The convergence of the calculation scheme starting from a fully irreducible vertex is demonstrated for the Anderson impurity model.

arXiv:1909.02793•Physical review. B./Physical review. B•Dec 26, 2019•42 citations

Multiloop flow equations for single-boson exchange fRG

Marcel Gievers, Elias Walter, Anxiang Ge +2 more

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

arXiv:2201.04878•The European Physical Journal B•Jul 1, 2022•23 citations

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

Two-particle level diagrammatic approaches for strongly correlated systems

Herbert F. Fotso

This thesis presents the development of new numerical methods for the treatment of strongly correlated electron systems based on self-consistent approaches at both the one and the two-particle level such as the parquet formalism. The parquet formalism was solved for the first time on a two-dimensional cluster. When the fully irreducible vertex is approximated by the bare vertex, we obtain the parquet approximation. Its validity was investigated by comparing results that it produces to those of other conserving approximations such as the FLuctuation EXchange (FLEX) approximation or the Second Order Perturbation Theory (SOPT). We found that it provides a significant improvement of FLEX or SOPT and a satisfactory agreement with Quantum Monte Carlo results despite instabilities in the self-consistency at low temperatures and for strong Coulomb interaction. We use the parquet formalism to study the Quantum Critical Point at finite doping in the Hubbard model by decomposing the vertex into its contributions from different channels. We apply this decomposition to the pairing channel and we find that the dominant contribution to the vertex originates in the spin channel even at the quantum critical doping. Furthermore, we explore the divergence of the two parts of the pairing matrix at optimal doping and observe that the irreducible vertex decreases monotonically as the doping is increased while the bare susceptibility exhibits an algebraic divergence at the quantum critical doping supporting the Quantum Critical BCS scenario proposed by She and Zaanen. To circumvent the instabilities in the iteration of the parquet formalism, we explored the dual fermion approach introduced by Rubtsov et al. Here, we extended the formalism to the Dynamical Cluster Approximation, in the process introducing a small parameter in the dual fermions perturbation theory. We demonstrate the quality of the resulting Dual Fermion DCA through a systematic study of the cluster size dependence and of the different perturbative approximations. These efforts represent the initial steps in the development of the Multi-Scale Many Body approach that appropriately treats correlations at different length scales.

Civil War Book Review•Jul 1, 2011•1 citation
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.
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
Taking the competition and the mutual screening of various bosonic fluctuations in correlated electron systems into account requires an unbiased approach to the many-body problem. One such approach is the self-consistent solution of the parquet equations, whose numerical treatment in lattice systems is however prohibitively expensive. In a recent article it was shown that there exists an alternative to the parquet decomposition of the four-point vertex function, which classifies the vertex diagrams according to the principle of single-boson exchange (SBE) [Krien, Valli, and Capone, Phys. Rev. B 100, 155149 (2019)]. Here we show that the SBE decomposition leads to a closed set of equations for the Hedin three-leg vertex, the polarization, and the electronic self-energy, which sums self-consistently the diagrams of the Maki-Thompson type. This circumvents the calculation of four-point vertex functions and the inversion of the Bethe-Salpeter equations, which are the two major bottlenecks of the parquet equations. The convergence of the calculation scheme starting from a fully irreducible vertex is demonstrated for the Anderson impurity model.
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 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.
This thesis presents the development of new numerical methods for the treatment of strongly correlated electron systems based on self-consistent approaches at both the one and the two-particle level such as the parquet formalism. The parquet formalism was solved for the first time on a two-dimensional cluster. When the fully irreducible vertex is approximated by the bare vertex, we obtain the parquet approximation. Its validity was investigated by comparing results that it produces to those of other conserving approximations such as the FLuctuation EXchange (FLEX) approximation or the Second Order Perturbation Theory (SOPT). We found that it provides a significant improvement of FLEX or SOPT and a satisfactory agreement with Quantum Monte Carlo results despite instabilities in the self-consistency at low temperatures and for strong Coulomb interaction. We use the parquet formalism to study the Quantum Critical Point at finite doping in the Hubbard model by decomposing the vertex into its contributions from different channels. We apply this decomposition to the pairing channel and we find that the dominant contribution to the vertex originates in the spin channel even at the quantum critical doping. Furthermore, we explore the divergence of the two parts of the pairing matrix at optimal doping and observe that the irreducible vertex decreases monotonically as the doping is increased while the bare susceptibility exhibits an algebraic divergence at the quantum critical doping supporting the Quantum Critical BCS scenario proposed by She and Zaanen. To circumvent the instabilities in the iteration of the parquet formalism, we explored the dual fermion approach introduced by Rubtsov et al. Here, we extended the formalism to the Dynamical Cluster Approximation, in the process introducing a small parameter in the dual fermions perturbation theory. We demonstrate the quality of the resulting Dual Fermion DCA through a systematic study of the cluster size dependence and of the different perturbative approximations. These efforts represent the initial steps in the development of the Multi-Scale Many Body approach that appropriately treats correlations at different length scales.