Pith. Sign in

REVIEW

Evolution of the spectrum of the Hubbard model with filling

Strongly Correlated Electrons · Submitted 23 Feb 2006 · 29 references

Pith's verdict: The claimed narrow band near the Fermi level is built on an uncontrolled interpolation at divergent frequencies, but the paper is transparent and the underlying physics is likely real.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDcond-mat/0602537 (version 1)
Pith IDF2GRJVET, Pith's permanent record of this paper and its review
CitationA. Sherman. "Evolution of the spectrum of the Hubbard model with filling." arXiv:cond-mat/0602537 (2006). https://pith.science/paper/F2GRJVET
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

Doping the Hubbard model below half-filling creates a narrow spin-polaron band near the Fermi level.

The reading

This paper tries to establish how the one-band Hubbard model's electron spectrum changes as the system is doped away from half-filling. The author develops a diagram expansion in the hopping constant whose building blocks are site cumulants of electron operators, derives a Larkin-type equation for the Green's function, and solves it in one loop on a 2D square lattice. The result is concrete: when the chemical potential falls low enough, the four bands already known at half-filling are joined by a fifth, narrow band at the Fermi level for momenta near the magnetic Brillouin-zone boundary. The new band's width is set by the superexchange constant $J=4t^2/U$, it becomes narrower as electrons are removed, and its dispersion follows the spin-polaron band of the $t$-$J$ model. If correct, this gives a direct route from the Hubbard model to the low-energy quasiparticle that governs the leading photoemission edge, without adding exchange terms by hand.

What carries the argument

The machinery is the site-cumulant diagram technique: every on-site average of electron creation and annihilation operators that appears in the expansion in the hopping constant $t$ is replaced by its connected part, a cumulant, so the expansion is organized by local correlations rather than by bare Wick contractions. The full Green's function satisfies the Larkin equation $G(k,z)=K(k,z)/[1-t_k K(k,z)]$, where $K$ is the sum of irreducible diagrams; the author keeps $K$ in one-loop order, with a $K_2$-type loop, and uses the Hubbard-I Green's function inside the loop. Divergencies in higher-order cumulants at $\omega=-\mu$ and $U-\mu$ make $\operatorname{Im} K$ positive there in this truncation; the paper repairs this by interpolating $\operatorname{Im} K$ between the regular regions, reconstructing the real part by Kramers-Kronig, and scaling the tails. The sharp features near those two frequencies are precisely what grows when the Fermi level enters the broad minimum of $\operatorname{Im} K$, producing the new band.

Extended reading notes

Core claim

The paper's central claim is that a narrow quasiparticle band appears near the Fermi level in the moderately to strongly repulsive one-band Hubbard model once the electron concentration is lowered past a threshold ($\mu \approx 0.17U$ for $t=-U/8$; $\mu\approx0.27U$ for $t=-U/4$). This band is distinct from the four bands present at half-filling; it forms for momenta near the boundary of the magnetic Brillouin zone, touches the Fermi level at that boundary, and has a bandwidth comparable to $J=4t^2/U$ that shrinks as the electron concentration drops. The author identifies this band with the spin-polaron band of the $t$-$J$ model: its dispersion, width, and filling dependence all match. The paper also argues that the low-energy feature seen only weakly in Monte-Carlo and cluster calculations is the same object, appearing here at somewhat larger doping because the one-loop approximation underestimates the electron concentration.

Load-bearing premise

The calculation assumes that the divergent pieces of higher-order cumulants at $\omega=-\mu$ and $U-\mu$ cancel in the full series, so that a hand interpolation of the damped function across these two points is an acceptable substitute for the exact answer — and the new band sits precisely at those points.

Editorial extensions

If this is right

  • The Hubbard model on its own produces a spin-polaron-like quasiparticle near the Fermi level at moderate hole doping; no explicit $t$-$J$ or Heisenberg exchange term is needed for this feature to exist.
  • The fifth band appears only after a doping threshold, and its bandwidth decreases with further electron removal, giving a testable doping dependence for the near-Fermi spectral weight.
  • Because the band's maximum energies sit at the magnetic Brillouin-zone boundary and touch the Fermi level there, the model predicts that the leading photoemission edge at moderate doping should be dominated by momenta near that boundary.
  • The threshold scales with $t/U$ ($\mu\approx0.17U$ for $t=-U/8$ versus $\mu\approx0.27U$ for $t=-U/4$), so the onset of the narrow band moves to higher chemical potential as hopping strengthens.
  • The cumulant and Larkin construction extends to multi-band Hubbard-type models such as a Cu-O plane, where similar narrow near-Fermi bands could be sought for each orbital subset.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the fifth band is truly a spin-polaron, its quasiparticle weight should be carried by local configurations with an antiparallel spin in the surrounding exchange cloud; a real-space spin-correlation calculation would test this beyond the width comparison.
  • The interpolation across $\omega=-\mu$ and $U-\mu$ is a free input; comparing two or more interpolation schemes would quantify how much of the band's intensity and dispersion is an artifact of that choice.
  • The paper's own concentration estimate shows the one-loop approximation underestimates electron number; applying a sum-rule-preserving chemical-potential shift would move the predicted band to lower hole concentrations, possibly reconciling it with cluster methods that see the feature at much smaller doping.
  • If the band touches the Fermi level at the magnetic Brillouin-zone boundary, then in a doped antiferromagnet the lowest hole excitations should live at those momenta; angle-resolved photoemission with sufficient resolution at low doping could look for that dispersion directly.
Original abstract

The diagram technique for the one-band Hubbard model is formulated for the case of moderate to strong Hubbard repulsion. The expansion in powers of the hopping constant is expressed in terms of site cumulants of electron creation and annihilation operators. For Green's function an equation of the Larkin type is derived and solved in a one-loop approximation for the case of two dimensions and nearest-neighbor hopping. With decreasing the electron concentration in addition to the four bands observed at half-filling a narrow band arises near the Fermi level. The dispersion of the new band, its bandwidth and the variation with filling are close to those of the spin-polaron band in the t-J model.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,224 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.48 · human participants 0.01 · controlled comparison 0.06 · randomized assignment 0.01 · states quantitative result 0.22 · compares with prior results 0.84 · mentions code or data availability 0.01.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The diagram technique for the one-band Hubbard model is formulated for the case of moderate to strong Hubbard repulsion. The expansion in powers of the hopping constant is expressed in terms of site cumulants of electron creation and annihilation operators. For Green's function an equation of the Larkin type is derived and solved in a one-loop approximation for the case of two dimensions and nearest-neighbor hopping. With decreasing the electron concentration in addition to the four bands observed at half-filling a narrow band arises near the Fermi level. The dispersion of the new band, its bandwidth and the variation with filling are close to those of the spin-polaron band in the t-J model.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The diagram technique for the one-band Hubbard model is formulated for the case of moderate to strong Hubbard repulsion. method 0.67 0.82 0.04 asserted 1.00
s1 The expansion in powers of the hopping constant is expressed in terms of site cumulants of electron creation and annihilation operators. method 0.97 0.84 0.03 asserted 1.00
s2 For Green's function an equation of the Larkin type is derived and solved in a one-loop approximation for the case of two dimensions and nearest-neighbor hopping. method 0.97 0.93 0.06 asserted 1.00
s3 With decreasing the electron concentration in addition to the four bands observed at half-filling a narrow band arises near the Fermi level. result 1.00 0.86 0.47 asserted 1.00
s4 The dispersion of the new band, its bandwidth and the variation with filling are close to those of the spin-polaron band in the t-J model. interpretation 0.56 0.81 0.35 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 0.92 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Diagrams of the first four orders of expansion (4). Averages which appear in higher-order terms of equati￾on (4) can be transformed in the same manner. The average in the k-th order term which con￾tains k + 1 creation and annihi￾lation operators is represented by the sum of all possible prod￾ucts of cumulants with the sum of orders equal to k + 1. All possible distributions of opera￾tors between the cumulants in these products have to be taken into account. The sign of a term in this sum is determined by the number of permutations of fermion operators, which bring the sequence of…
Figure 2. The imaginary part of K(ω) calculated using equations (10), (18) and (19) for a 100×100 lattice, t = −U/8 and T = 0.001U (the dashed lines). (a) µ = 0.5U, η = 0.01U. (b) µ = 0.1U, η = 0.02U. The solid lines show the corrected ℑK(ω) (see text). of ℑK(ω) at the singular frequencies can be reconstructed using an interpolation and its values in the regular region. Examples of such interpolation are given in figure 2. As seen from figure 2a, at half-filling, µ = U 2 , ℑK(ω) has two broad minima. With decreasing the chemical potential from this value these minima shift with respect to t
Figure 3. Panels (a) and (c): the spectral function A(kω) calculated for momenta along the symmetry lines of the square Brillouin zone in a 100×100 lattice for t = −U/8, T = 0.001U, η = 0.02U, µ = 0.2U (a) and µ = 0.05U (c). The dispersions of maxima in the panels (a) and (c) are shown in the panels (b) and (d), respectively. Here darker areas correspond to larger intensities of maxima. 10

Figures reproduced from arXiv:cond-mat/0602537 by the authors.

References
How the literature uses it

Papers whose main result relates to this one's

How each label was checked · JSON

Found among the papers closest to this one by title and main claim. Pith's classifier compared each pair's main results and a second model family confirmed the relation; checked against two further families, shown conflicts agreed 118 times in 120 and shown same results 218 times in 239. The quoted sentence is the other paper's main claim.

States the same result 1 paper

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.