Using a Green's-function equation of motion approach the two-dimensional Hubbard model is studied. The decoupling of the hierarchy is performed according to the conserving of the first four spectral moments. Beyond the Roth's two-pole approximation, the damping effect occurs naturally. The energy dispersions at various fillings, the spectral weight, and the density of states are calculated. Comparisons with the quantum Monte Carlo results and the cluster perturbation theory show that our approach can capture some essential features of the Hubbard model.Using a Green's-function equation of motion approach the two-dimensional Hubbard model is studied. The decoupling of the hierarchy is performed according to the conserving of the first four spectral moments. Beyond the Roth's two-pole approximation, the damping effect occurs naturally. The energy dispersions at various fillings, the spectral weight, and the density of states are calculated. Comparisons with the quantum Monte Carlo results and the cluster perturbation theory show that our approach can capture some essential features of the Hubbard model.

近代物理研究所和兰州重离子加速器实验室年报：英文版Jan 1, 20000 citations

Spectral function of the two-dimensional Hubbard model: A density matrix renormalization group plus cluster perturbation theory study

Chun Yang, Adrian E. Feiguin

We study the spectral function of the two-dimensional Hubbard model using cluster perturbation theory, and a density matrix renormalization group as a cluster solver. We reconstruct the two-dimensional dispersion at and away from half-filling using $2\ifmmode\times\else\texttimes\fi{}L$ ladders, with $L$ up to 80 sites, yielding results with unprecedented resolution in excellent agreement with quantum Monte Carlo. The main features of the spectrum can be described with a mean-field dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom.We study the spectral function of the two-dimensional Hubbard model using cluster perturbation theory, and a density matrix renormalization group as a cluster solver. We reconstruct the two-dimensional dispersion at and away from half-filling using $2\ifmmode\times\else\texttimes\fi{}L$ ladders, with $L$ up to 80 sites, yielding results with unprecedented resolution in excellent agreement with quantum Monte Carlo. The main features of the spectrum can be described with a mean-field dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom.

arXiv:1512.06102Physical review. B./Physical review. BFeb 24, 201630 citations

The sn-pole approximation in the Composite Operator Method

Ferdinando Mancini

A well-established method to deal with highly correlated systems is based on the expansion of the Green's function in terms of spectral moments. In the context of the Composite Operator Method one approximation is proposed: a set of n composite fields is assumed as fundamental basis and the dynamics is considered up to the order s. The resulting Green's function has a sn-pole structure. The truncation of the hierarchy of the equations of motion is made at the s-th order and the first s-1 equations are treated exactly. A theorem, which rules the conservation of the spectral moments, is presented. The procedure is applied to the Hubbard model and a recurrence relation for the calculation of its electronic spectral moments is derived.A well-established method to deal with highly correlated systems is based on the expansion of the Green's function in terms of spectral moments. In the context of the Composite Operator Method one approximation is proposed: a set of n composite fields is assumed as fundamental basis and the dynamics is considered up to the order s. The resulting Green's function has a sn-pole structure. The truncation of the hierarchy of the equations of motion is made at the s-th order and the first s-1 equations are treated exactly. A theorem, which rules the conservation of the spectral moments, is presented. The procedure is applied to the Hubbard model and a recurrence relation for the calculation of its electronic spectral moments is derived.

arXiv:cond-mat/0007341Jul 21, 20000 citations

The Hubbard Model in the Two-Pole Approximation

Adolfo Avella, Francesco Paolo Mancini, D. Villani +2 more

The two-dimensional Hubbard model is analyzed in the framework of the two-pole expansion. It is demonstrated that several theoretical approaches, when considered at their lowest level, are all equivalent and share the property of satisfying the conservation of the first four spectral momenta. It emerges that the various methods differ only in the way of fixing the internal parameters and that it exists a unique way to preserve simultaneously the Pauli principle and the particle–hole symmetry. A comprehensive comparison with respect to some general symmetry properties and the data from quantum Monte Carlo analysis shows the relevance of imposing the Pauli principle.The two-dimensional Hubbard model is analyzed in the framework of the two-pole expansion. It is demonstrated that several theoretical approaches, when considered at their lowest level, are all equivalent and share the property of satisfying the conservation of the first four spectral momenta. It emerges that the various methods differ only in the way of fixing the internal parameters and that it exists a unique way to preserve simultaneously the Pauli principle and the particle–hole symmetry. A comprehensive comparison with respect to some general symmetry properties and the data from quantum Monte Carlo analysis shows the relevance of imposing the Pauli principle.

arXiv:cond-mat/9708009International Journal of Modern Physics BJan 10, 199847 citations

A 4-pole approach to the Hubbard model within the Composite Operator Method

Adolfo Avella, Ferdinando Mancini, Gerardo Sica

The Hubbard model is studied within the Composite Operator Method framework in a 4-pole approximation. The operatorial basis is chosen according to both the hierarchy of equations of motion and the exact solution of the model reduced to the minimal cluster (2 sites) where all Hamiltonian terms are still active. Such a recipe amounts to include into the basis not only the two Hubbard operators, which usually take care of the scale of energy related to the on-site Coulomb repulsion U , but also two other non-local operators describing the Hubbard operators dressed by charge, spin and pair fluctuations on the nearest-neighboring sites and addressing the dynamically generated scale of energy related to the exchange J = 4 t 2 / U . The calculation framework is outlined and a very first comparison to numerical simulations, showing very good qualitative and quantitative agreement, is reported. As reference, the results in the 2-pole approximation are also reported and the differences discussed.The Hubbard model is studied within the Composite Operator Method framework in a 4-pole approximation. The operatorial basis is chosen according to both the hierarchy of equations of motion and the exact solution of the model reduced to the minimal cluster (2 sites) where all Hamiltonian terms are still active. Such a recipe amounts to include into the basis not only the two Hubbard operators, which usually take care of the scale of energy related to the on-site Coulomb repulsion U , but also two other non-local operators describing the Hubbard operators dressed by charge, spin and pair fluctuations on the nearest-neighboring sites and addressing the dynamically generated scale of energy related to the exchange J = 4 t 2 / U . The calculation framework is outlined and a very first comparison to numerical simulations, showing very good qualitative and quantitative agreement, is reported. As reference, the results in the 2-pole approximation are also reported and the differences discussed.

Journal of Physics Conference SeriesDec 14, 20121 citation

Spectral weight function for the two-dimensional Hubbard model

Steven R. White

We present results for the single-particle spectral weight function for the two-dimensional Hubbard model. Results were obtained using analytic expressions for the first two moments of the spectral weight function, along with numerical calculations using a combination of quantum Monte Carlo and the maximum-entropy method of analytic continuation. The expressions for the first two moments are surprisingly simple, involving only the parameters of the model, the chemical potential, and the filling. At half-filling, the squared width (variance) of the spectral weight is exactly ${\mathit{U}}^{2}$ /4. Numerical results are shown for 8\ifmmode\times\else\texttimes\fi{}8 lattices at a variety of fillings and momenta. The insulating antiferromagnetic gap in the spectral weights, readily apparent in the results at half-filling, is absent at a filling of 0.95.We present results for the single-particle spectral weight function for the two-dimensional Hubbard model. Results were obtained using analytic expressions for the first two moments of the spectral weight function, along with numerical calculations using a combination of quantum Monte Carlo and the maximum-entropy method of analytic continuation. The expressions for the first two moments are surprisingly simple, involving only the parameters of the model, the chemical potential, and the filling. At half-filling, the squared width (variance) of the spectral weight is exactly ${\mathit{U}}^{2}$/4. Numerical results are shown for 8\ifmmode\times\else\texttimes\fi{}8 lattices at a variety of fillings and momenta. The insulating antiferromagnetic gap in the spectral weights, readily apparent in the results at half-filling, is absent at a filling of 0.95.

Physical review. B, Condensed matterSep 1, 199185 citations

L

up to 80 sites, yielding results with unprecedented resolution in excellent agreement with quantum Monte Carlo. The main features of the spectrum can be described with a mean-field dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom.

{\mathit{U}}^{2}

/4. Numerical results are shown for 8\ifmmode\times\else\texttimes\fi{}8 lattices at a variety of fillings and momenta. The insulating antiferromagnetic gap in the spectral weights, readily apparent in the results at half-filling, is absent at a filling of 0.95.