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

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

Christian Eckhardt, G. A. H. Schober, Jannis Ehrlich, Carsten Honerkamp A

Published Aug 24, 2018• 26 citations
arXivHTMLPDF

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.

DOI
10.1103/physrevb.98.075143

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

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

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

Functional renormalization of spinless triangular-lattice fermions: N-patch vs. truncated-unity scheme

Nico Gneist, Dominik Kiese, Ravn Henkel +3 more

Abstract We study competing orders of spinless fermions in the triangular-lattice Hubbard model with nearest-neighbor interaction. We calculate the effective, momentum-resolved two-particle vertex in an unbiased way in terms of the functional renormalization group method and compare two different schemes for the momentum discretization, one based on dividing the Fermi surface into patches and one based on a channel decomposition. We study attractive and repulsive nearest-neighbor interaction and find a competition of pairing and charge instabilities. In the attractive case, a Pomeranchuk instability occurs at Van Hove filling and f-wave and p-wave pairing emerge when the filling is reduced. In the repulsive case, we obtain a charge density wave at Van Hove filling and extended p-wave pairing with reduced filling. The p-wave pairing solution is doubly degenerate and can realize chiral $$p+ip$$ p + i p superconductivity with different Chern numbers in the ground state. We discuss implications for strongly correlated spin-orbit coupled hexagonal electron systems such as moiré heterostructures. Graphic Abstract

arXiv:2205.12547•The European Physical Journal B•Sep 1, 2022•10 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
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 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
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.
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.
Abstract We study competing orders of spinless fermions in the triangular-lattice Hubbard model with nearest-neighbor interaction. We calculate the effective, momentum-resolved two-particle vertex in an unbiased way in terms of the functional renormalization group method and compare two different schemes for the momentum discretization, one based on dividing the Fermi surface into patches and one based on a channel decomposition. We study attractive and repulsive nearest-neighbor interaction and find a competition of pairing and charge instabilities. In the attractive case, a Pomeranchuk instability occurs at Van Hove filling and f-wave and p-wave pairing emerge when the filling is reduced. In the repulsive case, we obtain a charge density wave at Van Hove filling and extended p-wave pairing with reduced filling. The p-wave pairing solution is doubly degenerate and can realize chiral p+ipp+ipp+ip p + i p superconductivity with different Chern numbers in the ground state. We discuss implications for strongly correlated spin-orbit coupled hexagonal electron systems such as moiré heterostructures. Graphic 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.