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

Annual Conference on Magnetism and Magnetic Materials (MMM)

Projection Operator Formalism: Application to the Hubbard Model

A. J. Fedro, Robert S. Wilson

Published Jan 1, 1974• 0 citations
Publisher

Abstract

We present a method of doing many body theory based on the commutator projection operator discussed by Kim and Wilson. Such a method immediately yields the Dyson equation leaving one the task of evaluating the self energy term. Such a term was previously handled by ordinary perturbation theory. We find, on the other hand, an exact set of differential equations involving the self energy and its time derivatives in Zubarev‐like fashion which can be solved self‐consistently in every order of perturbation in contrast to chain breaking procedures which are not in general perturbationally correct. We apply our method to the Hubbard model and find the Esterling and Lange (EL) type result. That is, the form of EL is correct but their T − L(k⇒) expression is not complete. It is worthwhile to point out here that our calculation is done with finite hopping so we do not encounter the degeneracy problem of EL.

DOI
10.1063/1.3141798

Related papers

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

Degenerate Mass Operator Perturbation Theory in the Hubbard Model

D. M. Esterling, Robert Lange

The "atomic limit" (kinetic energy zero) of the Hubbard Hamiltonian has extreme degeneracy under very general circumstances. The kinetic energy or hopping term removes this degeneracy. It has been shown by Harris and Lange that there are important first-order manifestations of this fact. For instance, nonlocal correlations exist which invalidate usual Green's function truncation schemes. We present a mass operator perturbation scheme devised on the principle that, first of all, the role of kinetic energy as a breaker of degeneracy must be accounted for. A sort of Taylor's expansion of the mass operator in terms of the kinetic energy is found, the coefficients of which are functionals of the exact Green's functions at zero kinetic energy. The ambiguity due to the degeneracy is removed from these functions self consistently and thereby the degeneracy properly treated. We present a Green's function which accounts for all first-order effects of the kinetic energy.

Reviews of Modern Physics•Oct 1, 1968•35 citations

Nonperturbative solutions to the Hubbard Model

Lori Ann Thomas

Nonperturbative solutions to the Hubbard model are found by using the functional derivative method. A system of closed equations is obtained for the Green's functions of the Hubbard model. Exact expressions for the self-energy are derived which involve only first order functional derivatives. A scheme is proposed for systematically calculating self-energy corrections. We apply the solution to finite rings of two and four lattice sites and compare to the results of numerical calculations on these systems.

•Jan 1, 1987•0 citations

Functional-derivative study of the Hubbard model. I. Perturbation method and first-order approximation

Tadashi Arai, Morrel H. Cohen, Monica Tosi

In order to study the properties of the Hubbard model for narrow bands, a systematic treatment of the equations of motion of the Green's functions appropriate to that model has been developed. Higher-order Green's functions are reduced to functional derivatives of the basic Green's function $G$ and calculated iteratively in a perturbation scheme which takes the Hubbard I solution ${G}_{0}$ as the zeroth-order Green's function. A zeroth-order approximation to the self-energy correction obtained by inserting ${G}_{0}$ into the functional derivatives is compared with various existing solutions. The perturbation scheme is further extended to an infinite order and the self-energy is calculated exactly up to terms linear in the hopping motion $\ensuremath{\epsilon}$, a result which has not been obtained previously. The self-energy correction in this final result is drastically different from the zeroth-order solution, demonstrating the importance of the infinite-order iterative procedure. Finally, the electron correlations included in the final result are discussed in terms of diagrams.

Physical review. B, Solid state•Feb 15, 1977•8 citations

Functional-derivative study of the Hubbard model. II. Self-consistent equation and its complete solution

Tadashi Arai, Morrel H. Cohen

We develop a self-consistent method to solve the basic equation for the self-energy correction of the Hubbard model obtained in the preceding paper. The term $\ensuremath{\pi}[\ensuremath{\Delta}]$ involving second functional derivatives is neglected and the quantities $〈{N}_{R\ensuremath{\sigma}}(t)〉$ and $〈{C}_{R\ensuremath{\sigma}}^{\ifmmode\dagger\else\textdagger\fi{}}(t){C}_{{R}^{\ensuremath{'}\ensuremath{'}}\ensuremath{\sigma}}(t)〉$ are initially assumed to be independent of the external fields $\ensuremath{\epsilon}(\ensuremath{\sigma})$ and $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$. Under these restrictions, the complete self-energy correction is shown to be expanded in powers of $\ensuremath{\epsilon}(\ensuremath{\sigma})$ and $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$ in the form $\ensuremath{\Sigma}{R{R}^{\ensuremath{'}}\ensuremath{\sigma}}^{}(t{t}^{\ensuremath{'}})={\ensuremath{\xi}}_{0}(R{R}^{\ensuremath{'}}\overline{\ensuremath{\sigma}}t){\ensuremath{\delta}}_{t{t}^{\ensuremath{'}}}+\ensuremath{\Sigma}{n=0}^{}\ensuremath{\Sigma}{{R}^{\ensuremath{'}\ensuremath{'}}}^{}\ensuremath{\Sigma}{{R}^{\ensuremath{'}\ensuremath{'}\ensuremath{'}}}^{}{\ensuremath{\xi}}^{(n)}(R{R}^{\ensuremath{'}};{R}^{\ensuremath{'}\ensuremath{'}}{R}^{\ensuremath{'}\ensuremath{'}\ensuremath{'}})\ensuremath{\epsilon}({R}^{\ensuremath{'}\ensuremath{'}}{R}^{\ensuremath{'}\ensuremath{'}\ensuremath{'}}\ensuremath{\sigma}t){\ensuremath{\delta}}_{t{t}^{\ensuremath{'}}}$, where ${\ensuremath{\xi}}_{0}$ consists of all possible terms linear in $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$, while ${\ensuremath{\xi}}^{(n)}$ is made up of all possible terms of the $n$th degree in $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$. Equations for ${\ensuremath{\xi}}_{0}$ and ${\ensuremath{\xi}}^{(n)}$ are solved exactly and the resulting series is summed analytically, yielding a compact and complete analytic solution for the restricted equation. The part which is linear in $\ensuremath{\epsilon}$ is shown to be equal to the perturbation result obtained in the preceding paper, confirming the claim that the perturbation result is exact through terms linear in $\ensuremath{\epsilon}$. The method is extended and the effect of $\frac{\ensuremath{\delta}〈N〉}{\ensuremath{\delta}\ensuremath{\epsilon}}$ and $\frac{\ensuremath{\delta}〈{C}^{\ifmmode\dagger\else\textdagger\fi{}}C〉}{\ensuremath{\delta}\ensuremath{\epsilon}}$ is included. The effect is found to eliminate the difficulty that the value of one of the terms in the self-energy correction is abnormally overestimated in the previous result in the split-band, half-filled limit.

Physical review. B, Solid state•Feb 15, 1977•5 citations

On a Self‐Consistent Approximation Procedure with Application to the Hubbard Model

K. Elk

Abstract Considering a very general derivation of Dyson's equation an approximation procedure is discussed, which exhibits a summary of different known methods. Concerning the Hubbard model a generalization of Roth's procedure follows. In particular it is possible to obtain a self‐consistent method.

physica status solidi (b)•Aug 1, 1974•9 citations
We present a method of doing many body theory based on the commutator projection operator discussed by Kim and Wilson. Such a method immediately yields the Dyson equation leaving one the task of evaluating the self energy term. Such a term was previously handled by ordinary perturbation theory. We find, on the other hand, an exact set of differential equations involving the self energy and its time derivatives in Zubarev‐like fashion which can be solved self‐consistently in every order of perturbation in contrast to chain breaking procedures which are not in general perturbationally correct. We apply our method to the Hubbard model and find the Esterling and Lange (EL) type result. That is, the form of EL is correct but their T − L(k⇒) expression is not complete. It is worthwhile to point out here that our calculation is done with finite hopping so we do not encounter the degeneracy problem of EL.
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.
The "atomic limit" (kinetic energy zero) of the Hubbard Hamiltonian has extreme degeneracy under very general circumstances. The kinetic energy or hopping term removes this degeneracy. It has been shown by Harris and Lange that there are important first-order manifestations of this fact. For instance, nonlocal correlations exist which invalidate usual Green's function truncation schemes. We present a mass operator perturbation scheme devised on the principle that, first of all, the role of kinetic energy as a breaker of degeneracy must be accounted for. A sort of Taylor's expansion of the mass operator in terms of the kinetic energy is found, the coefficients of which are functionals of the exact Green's functions at zero kinetic energy. The ambiguity due to the degeneracy is removed from these functions self consistently and thereby the degeneracy properly treated. We present a Green's function which accounts for all first-order effects of the kinetic energy.
Nonperturbative solutions to the Hubbard model are found by using the functional derivative method. A system of closed equations is obtained for the Green's functions of the Hubbard model. Exact expressions for the self-energy are derived which involve only first order functional derivatives. A scheme is proposed for systematically calculating self-energy corrections. We apply the solution to finite rings of two and four lattice sites and compare to the results of numerical calculations on these systems.
In order to study the properties of the Hubbard model for narrow bands, a systematic treatment of the equations of motion of the Green's functions appropriate to that model has been developed. Higher-order Green's functions are reduced to functional derivatives of the basic Green's function GGG and calculated iteratively in a perturbation scheme which takes the Hubbard I solution G0{G}_{0}G0​ as the zeroth-order Green's function. A zeroth-order approximation to the self-energy correction obtained by inserting G0{G}_{0}G0​ into the functional derivatives is compared with various existing solutions. The perturbation scheme is further extended to an infinite order and the self-energy is calculated exactly up to terms linear in the hopping motion $\ensuremath{\epsilon}$, a result which has not been obtained previously. The self-energy correction in this final result is drastically different from the zeroth-order solution, demonstrating the importance of the infinite-order iterative procedure. Finally, the electron correlations included in the final result are discussed in terms of diagrams.
We develop a self-consistent method to solve the basic equation for the self-energy correction of the Hubbard model obtained in the preceding paper. The term $\ensuremath{\pi}[\ensuremath{\Delta}]$ involving second functional derivatives is neglected and the quantities $〈{N}_{R\ensuremath{\sigma}}(t)〉$ and $〈{C}_{R\ensuremath{\sigma}}^{\ifmmode\dagger\else\textdagger\fi{}}(t){C}_{{R}^{\ensuremath{'}\ensuremath{'}}\ensuremath{\sigma}}(t)〉$ are initially assumed to be independent of the external fields $\ensuremath{\epsilon}(\ensuremath{\sigma})$ and $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$. Under these restrictions, the complete self-energy correction is shown to be expanded in powers of $\ensuremath{\epsilon}(\ensuremath{\sigma})$ and $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$ in the form $\ensuremath{\Sigma}{R{R}^{\ensuremath{'}}\ensuremath{\sigma}}^{}(t{t}^{\ensuremath{'}})={\ensuremath{\xi}}_{0}(R{R}^{\ensuremath{'}}\overline{\ensuremath{\sigma}}t){\ensuremath{\delta}}_{t{t}^{\ensuremath{'}}}+\ensuremath{\Sigma}{n=0}^{}\ensuremath{\Sigma}{{R}^{\ensuremath{'}\ensuremath{'}}}^{}\ensuremath{\Sigma}{{R}^{\ensuremath{'}\ensuremath{'}\ensuremath{'}}}^{}{\ensuremath{\xi}}^{(n)}(R{R}^{\ensuremath{'}};{R}^{\ensuremath{'}\ensuremath{'}}{R}^{\ensuremath{'}\ensuremath{'}\ensuremath{'}})\ensuremath{\epsilon}({R}^{\ensuremath{'}\ensuremath{'}}{R}^{\ensuremath{'}\ensuremath{'}\ensuremath{'}}\ensuremath{\sigma}t){\ensuremath{\delta}}_{t{t}^{\ensuremath{'}}}$, where ${\ensuremath{\xi}}_{0}$ consists of all possible terms linear in $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$, while ${\ensuremath{\xi}}^{(n)}$ is made up of all possible terms of the nnnth degree in $\ensuremath{\epsilon}(\overline{\ensuremath{\sigma}})$. Equations for ${\ensuremath{\xi}}_{0}$ and ${\ensuremath{\xi}}^{(n)}$ are solved exactly and the resulting series is summed analytically, yielding a compact and complete analytic solution for the restricted equation. The part which is linear in $\ensuremath{\epsilon}$ is shown to be equal to the perturbation result obtained in the preceding paper, confirming the claim that the perturbation result is exact through terms linear in $\ensuremath{\epsilon}$. The method is extended and the effect of $\frac{\ensuremath{\delta}〈N〉}{\ensuremath{\delta}\ensuremath{\epsilon}}$ and $\frac{\ensuremath{\delta}〈{C}^{\ifmmode\dagger\else\textdagger\fi{}}C〉}{\ensuremath{\delta}\ensuremath{\epsilon}}$ is included. The effect is found to eliminate the difficulty that the value of one of the terms in the self-energy correction is abnormally overestimated in the previous result in the split-band, half-filled limit.