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

Combinatorial Mori-Zwanzig Theory

Yuanran Zhu

Published Oct 15, 2022• 4 citations
arXivHTMLPDF

Abstract

We introduce a combinatorial version Mori-Zwanzig theory and develop from it a family of self-consistent evolution equations for the correlation function or Green's function of interactive many-body systems. The core idea is to use an ansatz to rewrite the memory kernel (self-energy) of the regular Mori-Zwanzig equation as a function composition of the correlation (Green's) function. Then a series of algebraic combinatorial tools, especially the commutative and noncommutative Bell polynomials, are used to determine the exact Taylor series expansion of the composition function. The resulting combinatorial Mori-Zwanzig equation (CMZE) yields novel non-perturbative expansions of the equation of motion for the correlation (Green's) function. The structural equation for deriving such a combinatorial expansion resembles the combinatorial Dyson-Schwinger equation and may be viewed as its temporal-domain analogue. After introducing the abstract word and tree representation of the CMZE, we show its wide-range application in classical, stochastic, and quantum many-body systems. In all these examples, the new self-consistent expansions we obtained with the CMZE are similar to the diagrammatic skeleton expansions used in quantum many-body theory and lattice statistical field theory. We expect such a new framework can be used to calculate the correlation (Green's) function for strongly correlated/interactive many-body systems.

DOI
10.48550/arxiv.2210.08395

Related papers

Approximate but accurate quantum dynamics from the Mori formalism. II. Equilibrium time correlation functions

Andrés Montoya−Castillo, David R. Reichman

The ability to efficiently and accurately calculate equilibrium time correlation functions of many-body condensed phase quantum systems is one of the outstanding problems in theoretical chemistry. The Nakajima-Zwanzig-Mori formalism coupled to the self-consistent solution of the memory kernel has recently proven to be highly successful for the computation of nonequilibrium dynamical averages. Here, we extend this formalism to treat symmetrized equilibrium time correlation functions for the spin-boson model. Following the first paper in this series [A. Montoya-Castillo and D. R. Reichman, J. Chem. Phys. 144, 184104 (2016)], we use a Dyson-type expansion of the projected propagator to obtain a self-consistent solution for the memory kernel that requires only the calculation of normally evolved auxiliary kernels. We employ the approximate mean-field Ehrenfest method to demonstrate the feasibility of this approach. Via comparison with numerically exact results for the correlation function Czz(t)=Re⟨σz(0)σz(t)⟩, we show that the current scheme affords remarkable boosts in accuracy and efficiency over bare Ehrenfest dynamics. We further explore the sensitivity of the resulting dynamics to the choice of kernel closures and the accuracy of the initial canonical density operator.

arXiv:1610.04242•The Journal of Chemical Physics•Feb 24, 2017•40 citations

Self-consistent Dyson equation and self-energy functionals: An analysis and illustration on the example of the Hubbard atom

Walter Tarantino, Pina Romaniello, J. A. Berger +1 more

Perturbation theory using self-consistent Green's functions is one of the most widely used approaches to study many-body effects in condensed matter. On the basis of general considerations and by performing analytical calculations for the specific example of the Hubbard atom, we discuss some key features of this approach. We show that when the domain of the functionals that are used to realize the map between the noninteracting and the interacting Green's functions is properly defined, there exists a class of self-energy functionals for which the self-consistent Dyson equation has only one solution, which is the physical one. We also show that manipulation of the perturbative expansion of the interacting Green's function may lead to a wrong self-energy as a functional of the interacting Green's function, at least for some regions of the parameter space. These findings confirm and explain numerical results of Kozik et al. for the widely used skeleton series of Luttinger and Ward [Phys. Rev. Lett. 114, 156402 (2015)]. Our study shows that it is important to distinguish between the maps between sets of functions and the functionals that realize those maps. We demonstrate that the self-consistent Green's functions approach itself is not problematic, whereas the functionals that are widely used may have a limited range of validity.

arXiv:1703.05587•Physical review. B./Physical review. B•Jul 18, 2017•55 citations

Non-Hermitian Green’s function theory with N -body interactions: The coupled-cluster similarity transformation

Christopher J. N. Coveney, David P. Tew

We present the diagrammatic theory of the irreducible self-energy and Bethe-Salpeter kernel that naturally arises within the Green’s function formalism for a general N -body non-Hermitian interaction. In this work, we focus specifically on the coupled-cluster (CC) self-energy generated by the similarity transformation of the electronic structure Hamiltonian. The main result of this work is a Green’s function formalism for non-Hermitian many-body quantum systems. We develop the biorthogonal quantum theory to construct dynamical correlation functions where the time dependence of operators is governed by a non-Hermitian Hamiltonian. We extend the Gell-Mann and Low theorem to include non-Hermitian interactions and to generate perturbative expansions of many-body Green’s functions. We introduce the single-particle coupled-cluster Green’s function and derive the perturbative diagrammatic expansion for the non-Hermitian coupled-cluster self-energy in terms of the “noninteracting” reference Green’s function, Σ ̃ [ G 0 ] . From the exact equation of motion of the single-particle coupled-cluster Green’s function, we derive the self-consistent renormalized coupled-cluster self-energy, Σ ̃ [ G ̃ ] , and demonstrate its relationship to the perturbative expansion of the self-energy, Σ ̃ [ G 0 ] . Subsequently, we show that the usual electronic self-energy can be recovered from the coupled-cluster self-energy by neglecting the effects of the similarity transformation. We show how the coupled-cluster ground-state energy can be obtained from the coupled-cluster self-energy and provide an overview of the relationship between approximations for the coupled-cluster self-energy, ionization potential/electron affinity equation-of-motion coupled-cluster theory, and the G 0 W 0 approximation. As a result, we introduce the CC- G 0 W 0 self-energy by leveraging the connections between Green’s function and coupled-cluster theory. Finally, we derive the diagrammatic expansion of the coupled-cluster Bethe-Salpeter kernel.

arXiv:2503.06586•Physical Review Research•Feb 17, 2026•4 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

New self-consistent many-body perturbation theory: Application to the Hubbard model

A. J. Fedro, R. S. Wilson

A new many-body theory perturbationally consistent to all orders of the perturbation is presented and applied to the Hubbard problem as an example. Use is made of the recent Kim and Wilson commutator projection operator to derive the Dyson equation for the one-particle Green's function and an exact set of coupled differential equations for the damping or self-energy term which are solved self-consistently. The Hubbard $I$ result follows in our first approximation and the Esterling-Lange (EL)-type result in the second. Since we keep the hopping finite throughout the calculation, no degeneracy problem is encountered as in EL.

Physical review. B, Solid state•Mar 15, 1975•13 citations

Communication: The description of strong correlation within self-consistent Green's function second-order perturbation theory

Jordan J. Phillips, Dominika Zgid

We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.

arXiv:1405.3273•The Journal of Chemical Physics•Jun 25, 2014•123 citations
We introduce a combinatorial version Mori-Zwanzig theory and develop from it a family of self-consistent evolution equations for the correlation function or Green's function of interactive many-body systems. The core idea is to use an ansatz to rewrite the memory kernel (self-energy) of the regular Mori-Zwanzig equation as a function composition of the correlation (Green's) function. Then a series of algebraic combinatorial tools, especially the commutative and noncommutative Bell polynomials, are used to determine the exact Taylor series expansion of the composition function. The resulting combinatorial Mori-Zwanzig equation (CMZE) yields novel non-perturbative expansions of the equation of motion for the correlation (Green's) function. The structural equation for deriving such a combinatorial expansion resembles the combinatorial Dyson-Schwinger equation and may be viewed as its temporal-domain analogue. After introducing the abstract word and tree representation of the CMZE, we show its wide-range application in classical, stochastic, and quantum many-body systems. In all these examples, the new self-consistent expansions we obtained with the CMZE are similar to the diagrammatic skeleton expansions used in quantum many-body theory and lattice statistical field theory. We expect such a new framework can be used to calculate the correlation (Green's) function for strongly correlated/interactive many-body systems.
The ability to efficiently and accurately calculate equilibrium time correlation functions of many-body condensed phase quantum systems is one of the outstanding problems in theoretical chemistry. The Nakajima-Zwanzig-Mori formalism coupled to the self-consistent solution of the memory kernel has recently proven to be highly successful for the computation of nonequilibrium dynamical averages. Here, we extend this formalism to treat symmetrized equilibrium time correlation functions for the spin-boson model. Following the first paper in this series [A. Montoya-Castillo and D. R. Reichman, J. Chem. Phys. 144, 184104 (2016)], we use a Dyson-type expansion of the projected propagator to obtain a self-consistent solution for the memory kernel that requires only the calculation of normally evolved auxiliary kernels. We employ the approximate mean-field Ehrenfest method to demonstrate the feasibility of this approach. Via comparison with numerically exact results for the correlation function Czz(t)=Re⟨σz(0)σz(t)⟩, we show that the current scheme affords remarkable boosts in accuracy and efficiency over bare Ehrenfest dynamics. We further explore the sensitivity of the resulting dynamics to the choice of kernel closures and the accuracy of the initial canonical density operator.
Perturbation theory using self-consistent Green's functions is one of the most widely used approaches to study many-body effects in condensed matter. On the basis of general considerations and by performing analytical calculations for the specific example of the Hubbard atom, we discuss some key features of this approach. We show that when the domain of the functionals that are used to realize the map between the noninteracting and the interacting Green's functions is properly defined, there exists a class of self-energy functionals for which the self-consistent Dyson equation has only one solution, which is the physical one. We also show that manipulation of the perturbative expansion of the interacting Green's function may lead to a wrong self-energy as a functional of the interacting Green's function, at least for some regions of the parameter space. These findings confirm and explain numerical results of Kozik et al. for the widely used skeleton series of Luttinger and Ward [Phys. Rev. Lett. 114, 156402 (2015)]. Our study shows that it is important to distinguish between the maps between sets of functions and the functionals that realize those maps. We demonstrate that the self-consistent Green's functions approach itself is not problematic, whereas the functionals that are widely used may have a limited range of validity.
We present the diagrammatic theory of the irreducible self-energy and Bethe-Salpeter kernel that naturally arises within the Green’s function formalism for a general N -body non-Hermitian interaction. In this work, we focus specifically on the coupled-cluster (CC) self-energy generated by the similarity transformation of the electronic structure Hamiltonian. The main result of this work is a Green’s function formalism for non-Hermitian many-body quantum systems. We develop the biorthogonal quantum theory to construct dynamical correlation functions where the time dependence of operators is governed by a non-Hermitian Hamiltonian. We extend the Gell-Mann and Low theorem to include non-Hermitian interactions and to generate perturbative expansions of many-body Green’s functions. We introduce the single-particle coupled-cluster Green’s function and derive the perturbative diagrammatic expansion for the non-Hermitian coupled-cluster self-energy in terms of the “noninteracting” reference Green’s function, Σ ̃ [ G 0 ] . From the exact equation of motion of the single-particle coupled-cluster Green’s function, we derive the self-consistent renormalized coupled-cluster self-energy, Σ ̃ [ G ̃ ] , and demonstrate its relationship to the perturbative expansion of the self-energy, Σ ̃ [ G 0 ] . Subsequently, we show that the usual electronic self-energy can be recovered from the coupled-cluster self-energy by neglecting the effects of the similarity transformation. We show how the coupled-cluster ground-state energy can be obtained from the coupled-cluster self-energy and provide an overview of the relationship between approximations for the coupled-cluster self-energy, ionization potential/electron affinity equation-of-motion coupled-cluster theory, and the G 0 W 0 approximation. As a result, we introduce the CC- G 0 W 0 self-energy by leveraging the connections between Green’s function and coupled-cluster theory. Finally, we derive the diagrammatic expansion of the coupled-cluster Bethe-Salpeter kernel.
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.
A new many-body theory perturbationally consistent to all orders of the perturbation is presented and applied to the Hubbard problem as an example. Use is made of the recent Kim and Wilson commutator projection operator to derive the Dyson equation for the one-particle Green's function and an exact set of coupled differential equations for the damping or self-energy term which are solved self-consistently. The Hubbard III result follows in our first approximation and the Esterling-Lange (EL)-type result in the second. Since we keep the hopping finite throughout the calculation, no degeneracy problem is encountered as in EL.
We report an implementation of self-consistent Green's function many-body theory within a second-order approximation (GF2) for application with molecular systems. This is done by iterative solution of the Dyson equation expressed in matrix form in an atomic orbital basis, where the Green's function and self-energy are built on the imaginary frequency and imaginary time domain, respectively, and fast Fourier transform is used to efficiently transform these quantities as needed. We apply this method to several archetypical examples of strong correlation, such as a H32 finite lattice that displays a highly multireference electronic ground state even at equilibrium lattice spacing. In all cases, GF2 gives a physically meaningful description of the metal to insulator transition in these systems, without resorting to spin-symmetry breaking. Our results show that self-consistent Green's function many-body theory offers a viable route to describing strong correlations while remaining within a computationally tractable single-particle formalism.