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:2509.08501

Numerical Verification of Perturbative Schwinger-Dyson Resummation on Lattice Models

T. Banks, Anish Suresh

Published Sep 10, 2025• 0 citations
arXivHTMLPDF

Abstract

We investigate an approximation to the Schwinger-Dyson (SD) equations of the collective Coulomb field of the large $N$ Homogeneous Electron Fluid. The large $N$ approximation transforms the infinite SD hierarchy into a set of closed, equations for 1 and 2-pt correlators. In this paper, the dynamics of a toy model -- a small, square Euclidean lattice with periodic boundary conditions -- are considered. The Markov Chain Monte Carlo numerical method evaluated the 1 and 2-pt correlation functions on a $2 \times 2$ and $3 \times 3$ lattice. The derived equations are checked with the correlator values, and an agreement at $N \sim 100$ to order $10^{-3}$ was found. The agreement can be further strengthened by increasing runs in the Markov Chain Monte Carlo method.

DOI
10.48550/arxiv.2509.08501

Related papers

Schwinger Dyson Summation of Perturbation Expansions

Tom Banks

We introduce a hierarchical system of approximations for summing both conventional perturbation theory and large N vector expansions of models in quantum field theory and condensed matter physics. Each stage of the hierarchy consists of a closed set of nonlinear equations for one particle irreducible correlation functions with no more than K points and captures the perturbation expansion of each of those functions up to some finite order. As K goes to infinity the full Schwinger Dyson hierarchy is approached. We present an argument that for ordinary finite dimensional integrals, the procedure converges to the exact answer in this limit, despite the fact that no finite order of the hierarchy sees effects from non-leading saddle points of the integral. Some potential applications to the calculation of critical exponents, to low dimensional condensed matter systems, and to the phase structure of the homogeneous electron gas are outlined, but no detailed calculations are presented.

arXiv:2410.19683•Oct 25, 2024•0 citations

Schwinger-Dyson approach to nonequilibrium classical field theory

Krastan B. Blagoev, Fred Cooper, John F. Dawson +1 more

In this paper we discuss a Schwinger-Dyson (SD) approach for determining the time evolution of the unequal time correlation functions of a nonequilibrium classical field theory, where the classical system is described by an initial density matrix at time $t=0.$ We focus on $\ensuremath{\lambda}{\ensuremath{\varphi}}^{4}$ field theory in $1+1$ space-time dimensions where we can perform exact numerical simulations by sampling an ensemble of initial conditions specified by the initial density matrix. We discuss two approaches. The first, the bare vertex approximation (BVA), is based on ignoring vertex corrections to the SD equations in the auxiliary field formalism relevant for $1/N$ expansions. The second approximation is a related approximation made to the SD equations of the original formulation in terms of $\ensuremath{\varphi}$ alone. We compare these SD approximations as well as a Hartree approximation with exact numerical simulations. We find that both approximations based on the SD equations yield good agreement with exact numerical simulations and cure the late time oscillation problem of the Hartree approximation. We also discuss the relationship between the quantum and classical SD equations.

arXiv:hep-ph/0106195•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Nov 26, 2001•67 citations

Systematic Resummation of the Large N expansion of Vector Models: Application to the Hubbard model and $2 + 1 $ dimensional QED

Tom Banks

We introduce a hierarchy of closed equations for charge density correlation functions in the Hubbard model and $2 + 1$ dimensional QED. Each step in the hierarchy can be considered a large $N$ truncation of an exact, but infinite set of equations relating all $k-$point charge correlators. $N$ is the number of fermion spin components. Each step in the hierarchy sums up an infinite number of large $N$ diagrams, including all diagrams up to some fixed order, for $k$ point functions with $k \leq K$. Higher point functions are replaced with their leading large $N$ behavior. The simplest truncation gives a closed nonlinear equation for the $2$ point function.

arXiv:2003.13136•Mar 29, 2020•0 citations

Green functions in the renormalized many-body perturbation theory

V. Janiš

I review the way the many-body Green functions are used to renormalize the perturbation theory of correlated fermions. The Green functions are introduced to implement systematically dynamical corrections to the static mean-field theory. The renormalizations enter the perturbation theory via self-consistent evaluations of one-particle and two-particle Green functions. The one-particle self-consistency is discussed within the Baym-Kadanoff construction. We point to an inherent dichotomy of the Baym-Kadanoff approach in the relation between one and two-particle Green functions. They are the Schwinger-Dyson equation, directly derived from the generating Luttinger-Ward functional, and the Ward identity. The latter imposes a severe restriction on the two-particle irreducible vertex and is not guaranteed in the Baym-Kadanoff approach. No approximate theory is capable to obey both relations. I further discuss the two-particle approach with two-particle renormalizations that offers a way to reconcile the Schwinger-Dyson equation and the Ward identity. The two-particle irreducible vertex from the critical Bethe-Salpeter equation becomes a generating functional in the two-particle renormalization of the perturbation theory. The irreducible vertex is used in the Ward identity to determine the order parameter beyond the critical point and in the Bethe-Salpeter equation representing the vertex of the Schwinger-Dyson equation to determine the spectral self-energy in such a way that the thermodynamic critical behavior from both equations is qualitatively the same. Finally, I discuss a static mean-field-like approximation with a renormalized effective interaction determined self-consistently that leads to a qualitatively correct description of strong electron correlations in impurity and lattice models.

arXiv:2109.04231•Sep 9, 2021•0 citations

Green’s-function Monte Carlo for lattice fermions: Application to the t−J model

C. Stephen Hellberg, Efstratios Manousakis

We develop a general numerical method to study the zero-temperature properties of strongly correlated electron models on large lattices. The technique, which resembles Green's-function Monte Carlo, projects the ground-state component from a trial wave function with no approximations. We use this method to determine the phase diagram of the two-dimensional $t\ensuremath{-}J$ model, using the Maxwell construction to investigate electronic phase separation. The shell effects of fermions on finite-sized periodic lattices are minimized by keeping the number of electrons fixed at a closed-shell configuration and varying the size of the lattice. Results obtained for various electron numbers corresponding to different closed shells indicate that the finite-size effects in our calculation are small. For any value of interaction strength, we find that there is always a value of the electron density above which the system can lower its energy by forming a two-component phase separated state. Our results are compared with other calculations on the $t\ensuremath{-}J$ model. We find that the most accurate results are consistent with phase separation at all interaction strengths.

arXiv:cond-mat/9911071•Physical review. B, Condensed matter•May 1, 2000•46 citations

Schwinger-Dyson Equations-of-Motion Method for Quantum Spin-1/2

Chuanru Dai, Wenxin Ding

Abstract The essential difficulty of a large portion of the strongly interacting systems can be traced down to the “noncanonicality” of the underlying operators which are used to formulate the theory, like t-J model, Hubbard-type models in the large-U limit, quantum spin, etc. We proposed a new Schwinger-Dyson equations-of-motion method to compute the dynamical correlation functions of quantum spins at zero temperature. This method can also compute all the static observables from the dynamical correlation functions. We propose that this method can formulate dynamical perturbation calculations for quantum spin-1/2 models.

Journal of Physics Conference Series•Mar 1, 2023•0 citations
We investigate an approximation to the Schwinger-Dyson (SD) equations of the collective Coulomb field of the large NNN Homogeneous Electron Fluid. The large NNN approximation transforms the infinite SD hierarchy into a set of closed, equations for 1 and 2-pt correlators. In this paper, the dynamics of a toy model -- a small, square Euclidean lattice with periodic boundary conditions -- are considered. The Markov Chain Monte Carlo numerical method evaluated the 1 and 2-pt correlation functions on a 2×22 \times 22×2 and 3×33 \times 33×3 lattice. The derived equations are checked with the correlator values, and an agreement at N∼100N \sim 100N∼100 to order 10−310^{-3}10−3 was found. The agreement can be further strengthened by increasing runs in the Markov Chain Monte Carlo method.
We introduce a hierarchical system of approximations for summing both conventional perturbation theory and large N vector expansions of models in quantum field theory and condensed matter physics. Each stage of the hierarchy consists of a closed set of nonlinear equations for one particle irreducible correlation functions with no more than K points and captures the perturbation expansion of each of those functions up to some finite order. As K goes to infinity the full Schwinger Dyson hierarchy is approached. We present an argument that for ordinary finite dimensional integrals, the procedure converges to the exact answer in this limit, despite the fact that no finite order of the hierarchy sees effects from non-leading saddle points of the integral. Some potential applications to the calculation of critical exponents, to low dimensional condensed matter systems, and to the phase structure of the homogeneous electron gas are outlined, but no detailed calculations are presented.
In this paper we discuss a Schwinger-Dyson (SD) approach for determining the time evolution of the unequal time correlation functions of a nonequilibrium classical field theory, where the classical system is described by an initial density matrix at time t=0.t=0.t=0. We focus on $\ensuremath{\lambda}{\ensuremath{\varphi}}^{4}$ field theory in 1+11+11+1 space-time dimensions where we can perform exact numerical simulations by sampling an ensemble of initial conditions specified by the initial density matrix. We discuss two approaches. The first, the bare vertex approximation (BVA), is based on ignoring vertex corrections to the SD equations in the auxiliary field formalism relevant for 1/N1/N1/N expansions. The second approximation is a related approximation made to the SD equations of the original formulation in terms of $\ensuremath{\varphi}$ alone. We compare these SD approximations as well as a Hartree approximation with exact numerical simulations. We find that both approximations based on the SD equations yield good agreement with exact numerical simulations and cure the late time oscillation problem of the Hartree approximation. We also discuss the relationship between the quantum and classical SD equations.
We introduce a hierarchy of closed equations for charge density correlation functions in the Hubbard model and 2+12 + 12+1 dimensional QED. Each step in the hierarchy can be considered a large NNN truncation of an exact, but infinite set of equations relating all k−k-k−point charge correlators. NNN is the number of fermion spin components. Each step in the hierarchy sums up an infinite number of large NNN diagrams, including all diagrams up to some fixed order, for kkk point functions with k≤Kk \leq Kk≤K. Higher point functions are replaced with their leading large NNN behavior. The simplest truncation gives a closed nonlinear equation for the 222 point function.
I review the way the many-body Green functions are used to renormalize the perturbation theory of correlated fermions. The Green functions are introduced to implement systematically dynamical corrections to the static mean-field theory. The renormalizations enter the perturbation theory via self-consistent evaluations of one-particle and two-particle Green functions. The one-particle self-consistency is discussed within the Baym-Kadanoff construction. We point to an inherent dichotomy of the Baym-Kadanoff approach in the relation between one and two-particle Green functions. They are the Schwinger-Dyson equation, directly derived from the generating Luttinger-Ward functional, and the Ward identity. The latter imposes a severe restriction on the two-particle irreducible vertex and is not guaranteed in the Baym-Kadanoff approach. No approximate theory is capable to obey both relations. I further discuss the two-particle approach with two-particle renormalizations that offers a way to reconcile the Schwinger-Dyson equation and the Ward identity. The two-particle irreducible vertex from the critical Bethe-Salpeter equation becomes a generating functional in the two-particle renormalization of the perturbation theory. The irreducible vertex is used in the Ward identity to determine the order parameter beyond the critical point and in the Bethe-Salpeter equation representing the vertex of the Schwinger-Dyson equation to determine the spectral self-energy in such a way that the thermodynamic critical behavior from both equations is qualitatively the same. Finally, I discuss a static mean-field-like approximation with a renormalized effective interaction determined self-consistently that leads to a qualitatively correct description of strong electron correlations in impurity and lattice models.
We develop a general numerical method to study the zero-temperature properties of strongly correlated electron models on large lattices. The technique, which resembles Green's-function Monte Carlo, projects the ground-state component from a trial wave function with no approximations. We use this method to determine the phase diagram of the two-dimensional $t\ensuremath{-}J$ model, using the Maxwell construction to investigate electronic phase separation. The shell effects of fermions on finite-sized periodic lattices are minimized by keeping the number of electrons fixed at a closed-shell configuration and varying the size of the lattice. Results obtained for various electron numbers corresponding to different closed shells indicate that the finite-size effects in our calculation are small. For any value of interaction strength, we find that there is always a value of the electron density above which the system can lower its energy by forming a two-component phase separated state. Our results are compared with other calculations on the $t\ensuremath{-}J$ model. We find that the most accurate results are consistent with phase separation at all interaction strengths.
Abstract The essential difficulty of a large portion of the strongly interacting systems can be traced down to the “noncanonicality” of the underlying operators which are used to formulate the theory, like t-J model, Hubbard-type models in the large-U limit, quantum spin, etc. We proposed a new Schwinger-Dyson equations-of-motion method to compute the dynamical correlation functions of quantum spins at zero temperature. This method can also compute all the static observables from the dynamical correlation functions. We propose that this method can formulate dynamical perturbation calculations for quantum spin-1/2 models.