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

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

Johan Carlström

Published Sep 2, 2019• 1 citation
OpenAlex

Abstract

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.

Related papers

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

Johan Carlström

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.

arXiv:1909.00816•Physical review. B./Physical review. B•May 19, 2021•9 citations

Diagrammatics in the Dual Space, or There and Back Again

E. A. Stepanov

Accurately describing many-body effects in multi-orbital systems remains a major challenge in theoretical condensed matter physics. At present, there is a significant methodological gap between the numerical tools used in ab initio computational materials science and those developed to study strong electronic correlations. The former can treat realistic, large-scale systems but typically neglect many-body effects, while the latter focus on simplified models with only a few degrees of freedom, as only such models can be solved accurately in the presence of strong interactions. The purpose of this thesis is to bridge these two approaches and establish a systematic theoretical framework for realistic correlated electronic materials. This involves a full-cycle methodology that begins with constructing ab initio interacting models from density-functional theory, solving them using dynamical mean-field theory to capture local correlations, and extending beyond to incorporate non-local collective electronic fluctuations. To this end, we introduce the "dual" approach to strong correlations, which includes the dual fermion, dual boson, and dual triply irreducible local expansion methods. The central idea of the dual theories is to shift the reference point of the conventional Feynman diagrammatic expansion from a non-interacting electronic system to an interacting but exactly solvable one. Integrating out this reference system recasts the expansion in an effective dual space, where all diagrammatic building blocks are renormalized by the corresponding impurity quantities. This procedure transforms a non-perturbative expansion in the original variables into a perturbative expansion in terms of dual fermionic and bosonic fields, exact in both the weak- and strong-coupling limits. In this thesis, we collect and systematize the major developments of dual techniques achieved to date.

arXiv:2512.12389•Dec 13, 2025•0 citations

Strong-coupling approach for strongly correlated electron systems

Andrij M. Shvaika

A perturbation theory scheme in terms of electron hopping, which is based on the Wick theorem for Hubbard operators, is developed. Diagrammatic series contain single-site vertices connected by hopping lines and it is shown that for each vertex the problem splits into the subspaces with ``vacuum states'' determined by the diagonal Hubbard operators and only excitations around these vacuum states are allowed. The rules to construct diagrams are proposed. In the limit of infinite spatial dimensions the total auxiliary single-site problem exactly splits into subspaces that allows to build an analytical thermodynamically consistent approach for a Hubbard model. Some analytical results are given for the simple approximations when the two-pole (alloy-analogy solution) and four-pole (Hartree-Fock approximation) structure for Green's function is obtained. Two poles describe contribution from the Fermi-liquid component, which is dominant for small electron and hole concentrations (``overdoped case'' of high-${T}_{c}'\mathrm{s}),$ whereas other two describe contribution from the non-Fermi liquid and are dominant close to half-filling (``underdoped case'').

arXiv:cond-mat/9911255•Physical review. B, Condensed matter•Jul 15, 2000•14 citations

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

Grassmannization of classical models

Lode Pollet, Mikhail N. Kiselev, Nikolay V. Prokof'ev +1 more

Applying Feynman diagrammatics to non-fermionic strongly correlated models with local constraints might seem generically impossible for two separate reasons: (i) the necessity to have a Gaussian (non-interacting) limit on top of which the perturbative diagrammatic expansion is generated by Wick’s theorem, and (ii) Dyson’s collapse argument implying that the expansion in powers of coupling constant is divergent. We show that for arbitrary classical lattice models both problems can be solved/circumvented by reformulating the high-temperature expansion (more generally, any discrete representation of the model) in terms of Grassmann integrals. Discrete variables residing on either links, plaquettes, or sites of the lattice are associated with the Grassmann variables in such a way that the partition function (as well as all correlation functions) of the original system and its Grassmann-field counterpart are identical. The expansion of the latter around its Gaussian point generates Feynman diagrams. Our work paves the way for studying lattice gauge theories by treating bosonic and fermionic degrees of freedom on equal footing.

arXiv:1605.01953•New Journal of Physics•Nov 10, 2016•6 citations

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
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.
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.
Accurately describing many-body effects in multi-orbital systems remains a major challenge in theoretical condensed matter physics. At present, there is a significant methodological gap between the numerical tools used in ab initio computational materials science and those developed to study strong electronic correlations. The former can treat realistic, large-scale systems but typically neglect many-body effects, while the latter focus on simplified models with only a few degrees of freedom, as only such models can be solved accurately in the presence of strong interactions. The purpose of this thesis is to bridge these two approaches and establish a systematic theoretical framework for realistic correlated electronic materials. This involves a full-cycle methodology that begins with constructing ab initio interacting models from density-functional theory, solving them using dynamical mean-field theory to capture local correlations, and extending beyond to incorporate non-local collective electronic fluctuations. To this end, we introduce the "dual" approach to strong correlations, which includes the dual fermion, dual boson, and dual triply irreducible local expansion methods. The central idea of the dual theories is to shift the reference point of the conventional Feynman diagrammatic expansion from a non-interacting electronic system to an interacting but exactly solvable one. Integrating out this reference system recasts the expansion in an effective dual space, where all diagrammatic building blocks are renormalized by the corresponding impurity quantities. This procedure transforms a non-perturbative expansion in the original variables into a perturbative expansion in terms of dual fermionic and bosonic fields, exact in both the weak- and strong-coupling limits. In this thesis, we collect and systematize the major developments of dual techniques achieved to date.
A perturbation theory scheme in terms of electron hopping, which is based on the Wick theorem for Hubbard operators, is developed. Diagrammatic series contain single-site vertices connected by hopping lines and it is shown that for each vertex the problem splits into the subspaces with ``vacuum states'' determined by the diagonal Hubbard operators and only excitations around these vacuum states are allowed. The rules to construct diagrams are proposed. In the limit of infinite spatial dimensions the total auxiliary single-site problem exactly splits into subspaces that allows to build an analytical thermodynamically consistent approach for a Hubbard model. Some analytical results are given for the simple approximations when the two-pole (alloy-analogy solution) and four-pole (Hartree-Fock approximation) structure for Green's function is obtained. Two poles describe contribution from the Fermi-liquid component, which is dominant for small electron and hole concentrations (``overdoped case'' of high-Tc′s),{T}_{c}'\mathrm{s}),Tc′​s), whereas other two describe contribution from the non-Fermi liquid and are dominant close to half-filling (``underdoped case'').
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.
Applying Feynman diagrammatics to non-fermionic strongly correlated models with local constraints might seem generically impossible for two separate reasons: (i) the necessity to have a Gaussian (non-interacting) limit on top of which the perturbative diagrammatic expansion is generated by Wick’s theorem, and (ii) Dyson’s collapse argument implying that the expansion in powers of coupling constant is divergent. We show that for arbitrary classical lattice models both problems can be solved/circumvented by reformulating the high-temperature expansion (more generally, any discrete representation of the model) in terms of Grassmann integrals. Discrete variables residing on either links, plaquettes, or sites of the lattice are associated with the Grassmann variables in such a way that the partition function (as well as all correlation functions) of the original system and its Grassmann-field counterpart are identical. The expansion of the latter around its Gaussian point generates Feynman diagrams. Our work paves the way for studying lattice gauge theories by treating bosonic and fermionic degrees of freedom on equal footing.
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.