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

Strong-coupling diagrammatic Monte Carlo technique for correlated fermions and frustrated spins

Johan Carlström

Published May 19, 2021• 9 citations
arXivHTMLPDF

Abstract

We describe a controllable and unbiased strong-coupling diagrammatic Monte Carlo technique that is applicable to a wide range of fermionic systems and spin models. Unlike previous strong coupling methods that generally rely on the Grassmannian Hubbard-Stratonovich transformation, our construction is based on Wick's theorem and a recursive procedure to group contractions into effective connected vertices that are nonperturbative in all local physics and can be calculated exactly. The resulting expansion is described by simple diagrammatic rules that make it suitable for systematic treatment via stochastic sampling. Benchmarks against numerical linked cluster expansion display excellent agreement.

DOI
10.1103/physrevb.103.195147

Related papers

Strong-coupling diagrammatics for lattice fermions and spins based on irreducible vertices

Johan Carlström

We describe a controllable and unbiased strong coupling technique that is applicable to a wide range of fermionic systems and spin models. Unlike previous works that generally rely on the Grassmannian Hubbard-Stratonovich transformation, our construction is based on Wicks theorem and a recursive procedure to group contractions into effective irreducible vertices that are non-perturbative in all local physics and can be calculated exactly. The resulting expansion is described by simple diagrammatic rules that make it suitable for systematic treatment, and is independent of model parameters or doping. Comparison to approaches based on fermionization techniques indicates a dramatic advantage in terms of computational complexity, which will result in access to completely new parameter regimes and phenomena.

•Sep 2, 2019•1 citation

Diagrammatic Monte Carlo procedure for the spin-charge transformed Hubbard model

Johan Carlström

Using a dual representation of lattice fermion models that is based on spin-charge transformation and fermionization of the original description, I derive an algorithm for diagrammatic Monte Carlo simulation of strongly correlated systems. This scheme allows eliminating large expansion parameters, as well as large corrections to the density matrix that generally prevent diagrammatic methods from being efficient in this regime. As an example, I compute the filling factor for the Hubbard model at infinite on-site repulsion and compare the results to controllable data obtained from numerical linked-cluster expansion. I find excellent agreement between the two methods, as well as rapid convergence of the diagrammatic series. I also report results for the momentum distribution and kinetic energy of the electrons.

arXiv:1709.06415•Physical review. B./Physical review. B•Feb 12, 2018•18 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

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

Fedor Šimkovic, Evgeny Kozik

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.

arXiv:1712.10001•Physical review. B./Physical review. B•Sep 4, 2019•57 citations

Fluctuating field series: Towards calculations of correlated systems with high accuracy

Yana S. Lyakhova, S. D. Semenov, Alexander I. Lichtenstein +1 more

We introduce regular series expansion for weakly and moderately correlated fermionic systems based on the fluctuating local field approach. The method relies on explicitly accounting for the leading fluctuating mode(s) and is therefore suitable for medium-sized lattices. It can also be used as a solver for the cluster approximations for infinite-size lattices. We introduce a classical fluctuating field coupled to a fermionic collective mode(s). In this way, fluctuations in selected modes are treated in a nonperturbative way. Other degrees of freedom are accounted for in the diagram expansion performed for each value of the fluctuating field. The method is benchmarked for the $U/t=1$ and $U/t=2$ Hubbard lattices at half filling. Results for susceptibility in the antiferromagnetic channel along with the spectral function are compared with the numerically exact quantum Monte Carlo data. Calculations up to the third order of the series expansion are performed and show a uniform convergence to the reference result for the susceptibility. This convergence is observed well below the dynamical mean-field theory (DMFT) Ne\'el temperature and makes our practically simple method applicable in a much wider temperature range than DMFT-based diagrammatic schemes.

arXiv:2406.12590•Physical review. B./Physical review. B•Dec 20, 2024•3 citations

Diagrammatic Methods for Strongly-Correlated Electrons in Realistic Systems

Jia Li

This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.

•Jan 1, 2021•0 citations
We describe a controllable and unbiased strong-coupling diagrammatic Monte Carlo technique that is applicable to a wide range of fermionic systems and spin models. Unlike previous strong coupling methods that generally rely on the Grassmannian Hubbard-Stratonovich transformation, our construction is based on Wick's theorem and a recursive procedure to group contractions into effective connected vertices that are nonperturbative in all local physics and can be calculated exactly. The resulting expansion is described by simple diagrammatic rules that make it suitable for systematic treatment via stochastic sampling. Benchmarks against numerical linked cluster expansion display excellent agreement.
We describe a controllable and unbiased strong coupling technique that is applicable to a wide range of fermionic systems and spin models. Unlike previous works that generally rely on the Grassmannian Hubbard-Stratonovich transformation, our construction is based on Wicks theorem and a recursive procedure to group contractions into effective irreducible vertices that are non-perturbative in all local physics and can be calculated exactly. The resulting expansion is described by simple diagrammatic rules that make it suitable for systematic treatment, and is independent of model parameters or doping. Comparison to approaches based on fermionization techniques indicates a dramatic advantage in terms of computational complexity, which will result in access to completely new parameter regimes and phenomena.
Using a dual representation of lattice fermion models that is based on spin-charge transformation and fermionization of the original description, I derive an algorithm for diagrammatic Monte Carlo simulation of strongly correlated systems. This scheme allows eliminating large expansion parameters, as well as large corrections to the density matrix that generally prevent diagrammatic methods from being efficient in this regime. As an example, I compute the filling factor for the Hubbard model at infinite on-site repulsion and compare the results to controllable data obtained from numerical linked-cluster expansion. I find excellent agreement between the two methods, as well as rapid convergence of the diagrammatic series. I also report results for the momentum distribution and kinetic energy of the electrons.
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 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.
We introduce regular series expansion for weakly and moderately correlated fermionic systems based on the fluctuating local field approach. The method relies on explicitly accounting for the leading fluctuating mode(s) and is therefore suitable for medium-sized lattices. It can also be used as a solver for the cluster approximations for infinite-size lattices. We introduce a classical fluctuating field coupled to a fermionic collective mode(s). In this way, fluctuations in selected modes are treated in a nonperturbative way. Other degrees of freedom are accounted for in the diagram expansion performed for each value of the fluctuating field. The method is benchmarked for the U/t=1U/t=1U/t=1 and U/t=2U/t=2U/t=2 Hubbard lattices at half filling. Results for susceptibility in the antiferromagnetic channel along with the spectral function are compared with the numerically exact quantum Monte Carlo data. Calculations up to the third order of the series expansion are performed and show a uniform convergence to the reference result for the susceptibility. This convergence is observed well below the dynamical mean-field theory (DMFT) Ne\'el temperature and makes our practically simple method applicable in a much wider temperature range than DMFT-based diagrammatic schemes.
This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.