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

Determinant Monte Carlo for irreducible Feynman diagrams in the strongly correlated regime

Fedor Šimkovic, Evgeny Kozik

Published Sep 4, 2019• 57 citations
arXivHTMLPDF

Abstract

We develop a numerically exact method for the summation of irreducible Feynman diagrams for fermionic self-energy in the thermodynamic limit. The technique, based on the diagrammatic determinant Monte Carlo and its recent extension to connected diagrams, allows us to reach high $(\ensuremath{\sim}10)$ orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to high orders reveals a nontrivial analytic structure of the self-energy and enables its controlled reconstruction with arbitrary momentum resolution in the nonperturbative regime of essentially strong correlations, which has recently been reached with ultracold atoms in optical lattices.

DOI
10.1103/physrevb.100.121102

Related papers

Determinant Monte Carlo for irreducible Feynman diagrams

Fedor Šimkovic, Evgeny Kozik

We introduce a numerically exact method for efficient summation of irreducible Feynman diagrams for the fermionic self-energy in the thermodynamic limit. The technique, based on the Diagrammatic Determinant Monte Carlo, allows us to reach unprecedentedly high ($\sim 10$) orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to high orders reveals a non-trivial analytic structure of the self-energy and enables its controlled reconstruction with arbitrary momentum resolution in the non-perturbative regime of essentially strong correlations, which has recently been reached with ultracold atoms in optical lattices.

•Dec 28, 2017•0 citations

Determinant Monte Carlo algorithm for dynamical quantities in fermionic quantum systems

Alice Moutenet, Wéi Wú, Michel Ferrero

We introduce a Monte Carlo determinantal algorithm that allows to efficiently sample dynamical quantities, such as the self-energy, of fermionic quantum systems in their thermodynamic limit. The method expresses the sum of a factorial number of one-particle-irreducible diagrams as a recursive sum of determinants with exponential complexity. By comparing results of this algorithm for the lattice Hubbard model with those obtained from diagrammatic Monte Carlo and the recently introduced connected determinant approach, we show that we can reach higher perturbation orders and greater accuracy for the same computational effort.

arXiv:1712.10304?context=cond-mat•Dec 29, 2017•0 citations

Determinant Diagrammatic Monte Carlo Algorithm in the Thermodynamic Limit

Riccardo Rossi

We present a simple trick that allows us to consider the sum of all connected Feynman diagrams at fixed position of interaction vertices for general fermionic models, such that the thermodynamic limit can be taken analytically. With our approach one can achieve superior performance compared to conventional diagrammatic Monte Carlo algorithm, while rendering the algorithmic part dramatically simpler. By considering the sum of all connected diagrams at once, we allow for massive cancellations between different diagrams, greatly reducing the sign problem. In the end, the computational effort increases only exponentially with the order of the expansion, which should be contrasted with the factorial growth of the standard diagrammatic technique. We illustrate the efficiency of the technique for the two-dimensional Fermi-Hubbard model.

arXiv:1612.05184•Physical Review Letters•Jul 25, 2017•110 citations

Determinant Monte Carlo algorithms for dynamical quantities in fermionic systems

Alice Moutenet, Wéi Wú, Michel Ferrero

We introduce and compare three different Monte Carlo determinantal algorithms that allow one to compute dynamical quantities, such as the self-energy, of fermionic systems in their thermodynamic limit. We show that the most efficient approach expresses the sum of a factorial number of one-particle-irreducible diagrams as a recursive sum of determinants with exponential complexity. By comparing results for the two-dimensional Hubbard model with those obtained from state-of-the-art diagrammatic Monte Carlo, we show that we can reach higher perturbation orders and greater accuracy for the same computational effort.

arXiv:1712.10304•Physical review. B./Physical review. B•Feb 12, 2018•49 citations

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

Efficient one-loop-renormalized vertex expansions with connected determinant diagrammatic Monte Carlo

Fedor Šimkovic, Riccardo Rossi, Michel Ferrero

We present a technique that enables the evaluation of perturbative expansions based on one-loop-renormalized vertices up to large expansion orders. Specifically, we show how to compute large-order corrections to the random phase approximation in either the particle-hole or particle-particle channels. The algorithm's efficiency is achieved by the summation over contributions of all symmetrized Feynman diagram topologies using determinants, and by integrating out analytically the two-body long-range interactions in order to yield an effective zero-range interaction. Notably, the exponential scaling of the algorithm as a function of perturbation order leads to a polynomial scaling of the approximation error with computational time for a convergent series. To assess the performance of our approach, we apply it to the nonperturbative regime of the square-lattice fermionic Hubbard model away from half-filling and report, as compared to the bare interaction expansion algorithm, significant improvements of the Monte Carlo variance as well as the convergence properties of the resulting perturbative series.

arXiv:2008.07377•Physical review. B./Physical review. B•Nov 13, 2020•15 citations
We develop a numerically exact method for the summation of irreducible Feynman diagrams for fermionic self-energy in the thermodynamic limit. The technique, based on the diagrammatic determinant Monte Carlo and its recent extension to connected diagrams, allows us to reach high (\ensuremath∼10)(\ensuremath{\sim}10)(\ensuremath∼10) orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to high orders reveals a nontrivial analytic structure of the self-energy and enables its controlled reconstruction with arbitrary momentum resolution in the nonperturbative regime of essentially strong correlations, which has recently been reached with ultracold atoms in optical lattices.
We introduce a numerically exact method for efficient summation of irreducible Feynman diagrams for the fermionic self-energy in the thermodynamic limit. The technique, based on the Diagrammatic Determinant Monte Carlo, allows us to reach unprecedentedly high (∼10\sim 10∼10) orders of the weak-coupling expansion for the self-energy of the two-dimensional Hubbard model. Access to high orders reveals a non-trivial analytic structure of the self-energy and enables its controlled reconstruction with arbitrary momentum resolution in the non-perturbative regime of essentially strong correlations, which has recently been reached with ultracold atoms in optical lattices.
We introduce a Monte Carlo determinantal algorithm that allows to efficiently sample dynamical quantities, such as the self-energy, of fermionic quantum systems in their thermodynamic limit. The method expresses the sum of a factorial number of one-particle-irreducible diagrams as a recursive sum of determinants with exponential complexity. By comparing results of this algorithm for the lattice Hubbard model with those obtained from diagrammatic Monte Carlo and the recently introduced connected determinant approach, we show that we can reach higher perturbation orders and greater accuracy for the same computational effort.
We present a simple trick that allows us to consider the sum of all connected Feynman diagrams at fixed position of interaction vertices for general fermionic models, such that the thermodynamic limit can be taken analytically. With our approach one can achieve superior performance compared to conventional diagrammatic Monte Carlo algorithm, while rendering the algorithmic part dramatically simpler. By considering the sum of all connected diagrams at once, we allow for massive cancellations between different diagrams, greatly reducing the sign problem. In the end, the computational effort increases only exponentially with the order of the expansion, which should be contrasted with the factorial growth of the standard diagrammatic technique. We illustrate the efficiency of the technique for the two-dimensional Fermi-Hubbard model.
We introduce and compare three different Monte Carlo determinantal algorithms that allow one to compute dynamical quantities, such as the self-energy, of fermionic systems in their thermodynamic limit. We show that the most efficient approach expresses the sum of a factorial number of one-particle-irreducible diagrams as a recursive sum of determinants with exponential complexity. By comparing results for the two-dimensional Hubbard model with those obtained from state-of-the-art diagrammatic Monte Carlo, we show that we can reach higher perturbation orders and greater accuracy for the same computational effort.
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 present a technique that enables the evaluation of perturbative expansions based on one-loop-renormalized vertices up to large expansion orders. Specifically, we show how to compute large-order corrections to the random phase approximation in either the particle-hole or particle-particle channels. The algorithm's efficiency is achieved by the summation over contributions of all symmetrized Feynman diagram topologies using determinants, and by integrating out analytically the two-body long-range interactions in order to yield an effective zero-range interaction. Notably, the exponential scaling of the algorithm as a function of perturbation order leads to a polynomial scaling of the approximation error with computational time for a convergent series. To assess the performance of our approach, we apply it to the nonperturbative regime of the square-lattice fermionic Hubbard model away from half-filling and report, as compared to the bare interaction expansion algorithm, significant improvements of the Monte Carlo variance as well as the convergence properties of the resulting perturbative series.