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:cond-mat/0612196 · Physical Review B

Dual fermion approach to nonlocal correlations in the Hubbard model

Alexey N. Rubtsov, Mikhail I. Katsnelson, Alexander I. Lichtenstein

Published Jan 11, 2008• 350 citations
arXivHTMLPDF

Abstract

A diagrammatic technique is developed to describe nonlocal effects (e.g., pseudogap formation) in the Hubbard-like models. In contrast to cluster approaches, this method utilizes an exact transition to the dual set of variables, and it therefore becomes possible to treat the irreducible vertices of an effective single-impurity problem as small parameters. This provides a very efficient interpolation between weak coupling (band) and atomic limits. The antiferromagnetic pseudogap formation in the Hubbard model is correctly reproduced by just the lowest-order diagrams.

DOI
10.1103/physrevb.77.033101

Related papers

Dual fermion approach to the two-dimensional Hubbard model: Antiferromagnetic fluctuations and Fermi arcs

Alexey N. Rubtsov, Mikhail I. Katsnelson, Alexander I. Lichtenstein +1 more

We present an efficient diagrammatic method to describe nonlocal correlation effects in lattice fermion Hubbard-type models, which is based on a change of variables in the Grassmann path integrals. The new fermions are dual to the original ones and correspond to weakly interacting quasiparticles in the case of strong local correlations in the Hubbard model. The method starts with dynamical mean-field theory as a zeroth-order approximation and includes nonlocal effects in a perturbative way. In contrast to cluster approaches, this method utilizes an exact transition to a dual set of variables. It therefore becomes possible to treat vertices of an effective single-impurity problem as small parameters. This provides a very efficient interpolation between bandlike weak-coupling and atomic limits. The method is illustrated on the two-dimensional Hubbard model. The antiferromagnetic pseudogap, Fermi-arc formations, and non-Fermi-liquid effects due to the Van Hove singularity are correctly reproduced by the lowest-order diagrams. Extremum properties of the dual fermion approach are discussed in terms of the Feynman variational principle.

arXiv:0810.3819•Physical Review B•Jan 30, 2009•144 citations

Non-Perturbative Many-Body Approach to the Hubbard Model and Single-Particle Pseudogap

Yu. M. Vil’k, A.–M. S. Tremblay

A new approach to the single-band Hubbard model is described in the general context of many-body theories. It is based on enforcing conservation laws, the Pauli principle and a number of crucial sum-rules. More specifically, spin and charge susceptibilities are expressed, in a conserving approximation, as a function of two constant irreducible vertices whose values are found self-consistently. The Mermin-Wagner theorem in two dimensions is automatically satisfied. The effect of collective modes on single-particle properties is then obtained by a paramagnon-like formula that is consistent with the two-particle properties in the sense that the potential energy obtained from $TrΣG$ is identical to that obtained using the fluctuation-dissipation theorem for susceptibilities. The vertex corrections are included through constant irreducible vertices. The theory is in quantitative agreement with Monte Carlo simulations for both single-particle and two-particle properties. In the two-dimensional renormalized classical regime, spin fluctuations lead to precursors of antiferromagnetic bands (shadow bands) and to the destruction of the Fermi-liquid quasiparticles in a wide temperature range above the zero-temperature phase transition. The analogous phenomenon of pairing pseudogap can occur in the attractive model in two dimensions when the pairing fluctuations become critical. Other many-body approaches are critically compared. It is argued that treating the spin fluctuations as if there was a Migdal's theorem can lead to wrong predictions, in particular with regard to the the single-particle pseudogap.

arXiv:cond-mat/9702188•Journal de Physique I•Nov 1, 1997•336 citations

Dual-space cluster-diagrammatic approach to nonlocal electronic correlations

Félix Fossati, E. A. Stepanov

In this work, we extend the dual triply irreducible local expansion (D-TRILEX) approach for correlated electronic systems by introducing a cluster reference system for the diagrammatic expansion. This framework allows us to consistently combine the exact treatment of short-range correlation effects within the cluster, with an efficient diagrammatic description of the long-range charge and spin collective fluctuations beyond the cluster. We demonstrate the effectiveness of our approach by applying it to the one-dimensional nano-ring Hubbard model, where the low dimensionality enhances non-local correlations. Our results show that the cluster extension of D-TRILEX accurately reproduces the electronic self-energy at momenta corresponding to the Fermi energy, in good agreement with the reliable Hirsch-Fye quantum Monte Carlo solution of the problem. We further compare this method with the more computationally demanding parquet dynamical vertex approximation and find that, our method yields substantially more accurate results at momenta associated with the Fermi surface. We show that the D-TRILEX diagrammatic extension drastically reduces the periodization ambiguity of cluster quantities when mapping back to the original lattice, compared to cluster dynamical mean-field theory (CDMFT). Furthermore, we identify the CDMFT impurity problem as the main source of the translational-symmetry breaking and propose a computational scheme for improving the starting point for the cluster-diagrammatic expansion.

arXiv:2507.06015•Physical review. B./Physical review. B•Feb 2, 2026•1 citation

Diagrammatic Monte Carlo approach for diagrammatic extensions of dynamical mean-field theory: Convergence analysis of the dual fermion technique

Jan Gukelberger, Evgeny Kozik, Hartmut Hafermann

The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.

arXiv:1611.07523•Physical review. B./Physical review. B•Jul 28, 2017•55 citations

Diagrammatic Monte Carlo for dual fermions

Sergei Iskakov, Andrey E. Antipov, Emanuel Gull

We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong coupling regime of the half-filled Hubbard model in two dimensions, illustrating that the method converges quickly where dynamical mean field theory is a good approximation, and show that corrections are large in the strong correlation regime at intermediate interaction. The fast convergence of dual corrections to dynamical mean field results illustrates the power of the approach and opens a practical avenue towards the systematic inclusion of nonlocal correlations in correlated materials simulations. An analysis of the frequency scale shows that only low-frequency propagators contribute substantially to the diagrams, putting the inclusion of higher order vertices within reach.

arXiv:1607.01401•Physical review. B./Physical review. B•Jul 5, 2016•50 citations

Cluster dual fermion approach to nonlocal correlations

Hartmut Hafermann, S. Brener, A. N. Rubtsov +2 more

A general cluster dual fermion approach to nonlocal correlations in crystals is formulated. The scheme allows the treatment of long-range correlations beyond the cluster dynamical mean field theory and nonlocal effects in realistic calculations of multiorbital systems. The simplest approximation is shown to exactly correspond to the free-cluster dynamical mean field theory. This approach is applied to the one-dimensional Hubbard model. Already the first dual-fermion correction to the free cluster leads to a drastic improvement of the calculated Green’s function.

arXiv:0707.4022•Journal of Experimental and Theoretical Physics Letters•Jan 1, 2008•37 citations
A diagrammatic technique is developed to describe nonlocal effects (e.g., pseudogap formation) in the Hubbard-like models. In contrast to cluster approaches, this method utilizes an exact transition to the dual set of variables, and it therefore becomes possible to treat the irreducible vertices of an effective single-impurity problem as small parameters. This provides a very efficient interpolation between weak coupling (band) and atomic limits. The antiferromagnetic pseudogap formation in the Hubbard model is correctly reproduced by just the lowest-order diagrams.
We present an efficient diagrammatic method to describe nonlocal correlation effects in lattice fermion Hubbard-type models, which is based on a change of variables in the Grassmann path integrals. The new fermions are dual to the original ones and correspond to weakly interacting quasiparticles in the case of strong local correlations in the Hubbard model. The method starts with dynamical mean-field theory as a zeroth-order approximation and includes nonlocal effects in a perturbative way. In contrast to cluster approaches, this method utilizes an exact transition to a dual set of variables. It therefore becomes possible to treat vertices of an effective single-impurity problem as small parameters. This provides a very efficient interpolation between bandlike weak-coupling and atomic limits. The method is illustrated on the two-dimensional Hubbard model. The antiferromagnetic pseudogap, Fermi-arc formations, and non-Fermi-liquid effects due to the Van Hove singularity are correctly reproduced by the lowest-order diagrams. Extremum properties of the dual fermion approach are discussed in terms of the Feynman variational principle.
A new approach to the single-band Hubbard model is described in the general context of many-body theories. It is based on enforcing conservation laws, the Pauli principle and a number of crucial sum-rules. More specifically, spin and charge susceptibilities are expressed, in a conserving approximation, as a function of two constant irreducible vertices whose values are found self-consistently. The Mermin-Wagner theorem in two dimensions is automatically satisfied. The effect of collective modes on single-particle properties is then obtained by a paramagnon-like formula that is consistent with the two-particle properties in the sense that the potential energy obtained from TrΣGTrΣGTrΣG is identical to that obtained using the fluctuation-dissipation theorem for susceptibilities. The vertex corrections are included through constant irreducible vertices. The theory is in quantitative agreement with Monte Carlo simulations for both single-particle and two-particle properties. In the two-dimensional renormalized classical regime, spin fluctuations lead to precursors of antiferromagnetic bands (shadow bands) and to the destruction of the Fermi-liquid quasiparticles in a wide temperature range above the zero-temperature phase transition. The analogous phenomenon of pairing pseudogap can occur in the attractive model in two dimensions when the pairing fluctuations become critical. Other many-body approaches are critically compared. It is argued that treating the spin fluctuations as if there was a Migdal's theorem can lead to wrong predictions, in particular with regard to the the single-particle pseudogap.
In this work, we extend the dual triply irreducible local expansion (D-TRILEX) approach for correlated electronic systems by introducing a cluster reference system for the diagrammatic expansion. This framework allows us to consistently combine the exact treatment of short-range correlation effects within the cluster, with an efficient diagrammatic description of the long-range charge and spin collective fluctuations beyond the cluster. We demonstrate the effectiveness of our approach by applying it to the one-dimensional nano-ring Hubbard model, where the low dimensionality enhances non-local correlations. Our results show that the cluster extension of D-TRILEX accurately reproduces the electronic self-energy at momenta corresponding to the Fermi energy, in good agreement with the reliable Hirsch-Fye quantum Monte Carlo solution of the problem. We further compare this method with the more computationally demanding parquet dynamical vertex approximation and find that, our method yields substantially more accurate results at momenta associated with the Fermi surface. We show that the D-TRILEX diagrammatic extension drastically reduces the periodization ambiguity of cluster quantities when mapping back to the original lattice, compared to cluster dynamical mean-field theory (CDMFT). Furthermore, we identify the CDMFT impurity problem as the main source of the translational-symmetry breaking and propose a computational scheme for improving the starting point for the cluster-diagrammatic expansion.
The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.
We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong coupling regime of the half-filled Hubbard model in two dimensions, illustrating that the method converges quickly where dynamical mean field theory is a good approximation, and show that corrections are large in the strong correlation regime at intermediate interaction. The fast convergence of dual corrections to dynamical mean field results illustrates the power of the approach and opens a practical avenue towards the systematic inclusion of nonlocal correlations in correlated materials simulations. An analysis of the frequency scale shows that only low-frequency propagators contribute substantially to the diagrams, putting the inclusion of higher order vertices within reach.
A general cluster dual fermion approach to nonlocal correlations in crystals is formulated. The scheme allows the treatment of long-range correlations beyond the cluster dynamical mean field theory and nonlocal effects in realistic calculations of multiorbital systems. The simplest approximation is shown to exactly correspond to the free-cluster dynamical mean field theory. This approach is applied to the one-dimensional Hubbard model. Already the first dual-fermion correction to the free cluster leads to a drastic improvement of the calculated Green’s function.