Skip to content

Hubbard Bands and Spectral-Weight Transfer

Hubbard bands are many-body electron-addition and electron-removal sectors, not ordinary one-electron bands separated by a larger gap. Correlations can move weight between those sectors while the exact total remains fixed. Three quantities must therefore be kept separate: an integrated spectral weight, the residue of a coherent quasiparticle pole, and the charge gap or compressibility that decides whether the state is insulating.

Required background. Use the Hubbard atomic and strong-coupling limits and spectral moments and sum rules. Helpful background. Dyson equations define self-energy poles and quasiparticle residues.

Begin with one Hubbard site in grand-canonical form,

Kat=HatμN=Unnμ(n+n).K_{\mathrm{at}}=H_{\mathrm{at}}-\mu N =U n_\uparrow n_\downarrow -\mu(n_\uparrow+n_\downarrow).

Here zz is a complex frequency, GR(ω)=G(ω+i0+)G^R(\omega)=G(\omega+i0^+), and ω=0\omega=0 is the declared chemical potential. At zero temperature it is useful to retain the Lehmann labels “addition” and “removal” explicitly. With nonnegative excitation costs Ωm±\Omega_m^\pm measured using K=HμNK=H-\mu N,

Aiσ+(ω)=2πmm,N+1ciσ0,N2δ(ωΩm+),A^+_{i\sigma}(\omega) =2\pi\sum_m \left|\langle m,N+1|c_{i\sigma}^\dagger|0,N\rangle\right|^2 \delta(\omega-\Omega_m^+), Aiσ(ω)=2πmm,N1ciσ0,N2δ(ω+Ωm),A=A++A.A^-_{i\sigma}(\omega) =2\pi\sum_m \left|\langle m,N-1|c_{i\sigma}|0,N\rangle\right|^2 \delta(\omega+\Omega_m^-), \qquad A=A^++A^-.

Usually A+A^+ lies above and AA^- below the chemical potential. The sector labels are nevertheless more fundamental than the sign of ω\omega: at an exact level crossing, an addition pole and a removal pole can both occur at ω=0\omega=0.

For spin σ\sigma, adding or removing an electron costs one energy when the opposite spin is absent and another when it is present. Conditioning on those two sectors gives the exact atomic Green function

Gat,σ(z)=1nˉσˉz+μ+nˉσˉz+μU,nˉσˉ=nσˉ.G_{\mathrm{at},\sigma}(z) =\frac{1-\bar n_{\bar\sigma}}{z+\mu} +\frac{\bar n_{\bar\sigma}}{z+\mu-U}, \qquad \bar n_{\bar\sigma}=\langle n_{\bar\sigma}\rangle.

Equivalently,

Gat,σ(z)=z+μU(1nˉσˉ)(z+μ)(z+μU).G_{\mathrm{at},\sigma}(z) =\frac{z+\mu-U(1-\bar n_{\bar\sigma})} {(z+\mu)(z+\mu-U)}.

With Aσ(ω)=2ImGσR(ω)A_\sigma(\omega)=-2\operatorname{Im}G^R_\sigma(\omega), the spectrum per spin is

Aat,σ(ω)=2π(1nˉσˉ)δ(ω+μ)+2πnˉσˉδ(ω+μU).A_{\mathrm{at},\sigma}(\omega) =2\pi(1-\bar n_{\bar\sigma})\delta(\omega+\mu) +2\pi\bar n_{\bar\sigma}\delta(\omega+\mu-U).

The residues are occupation probabilities, not adjustable peak heights. They obey

dω2πAat,σ(ω)=1.\int_{-\infty}^{\infty}\frac{\mathrm d\omega}{2\pi} A_{\mathrm{at},\sigma}(\omega)=1.

At paramagnetic half filling, nˉσˉ=1/2\bar n_{\bar\sigma}=1/2 and particle–hole symmetry fixes μ=U/2\mu=U/2. The poles then sit at ω=±U/2\omega=\pm U/2 with equal per-spin weight:

Gat,σ(z)=12(1z+U/2+1zU/2)=zz2(U/2)2,G_{\mathrm{at},\sigma}(z) =\frac12\left(\frac1{z+U/2}+\frac1{z-U/2}\right) =\frac{z}{z^2-(U/2)^2}, Aat,σ(ω)=πδ(ω+U/2)+πδ(ωU/2).A_{\mathrm{at},\sigma}(\omega) =\pi\delta(\omega+U/2)+\pi\delta(\omega-U/2).

Hopping broadens these exact transitions into lower and upper Hubbard bands. It also changes their internal moments and, away from the filled-lower-band half-filled limit, can change their integrated weights. Calling the effect “broadening” alone misses that redistribution.

The apparent form of the atomic self-energy depends on what is included in the reference propagator. For the unshifted Hamiltonian above, take

G0,σ1(z)=z+μ,Gσ1=G0,σ1Σfull,σ.G_{0,\sigma}^{-1}(z)=z+\mu, \qquad G_\sigma^{-1}=G_{0,\sigma}^{-1}-\Sigma_{\mathrm{full},\sigma}.

Inverting the exact Green function gives

Σfull,σ(z)=Unˉσˉ+U2nˉσˉ(1nˉσˉ)z+μU(1nˉσˉ).\Sigma_{\mathrm{full},\sigma}(z) =U\bar n_{\bar\sigma} +\frac{U^2\bar n_{\bar\sigma}(1-\bar n_{\bar\sigma})} {z+\mu-U(1-\bar n_{\bar\sigma})}.

The first term is the static Hartree contribution. If it is absorbed into a shifted chemical potential,

μH,σ=μUnˉσˉ,Σcorr,σ=Σfull,σUnˉσˉ,\mu_{\mathrm H,\sigma}=\mu-U\bar n_{\bar\sigma}, \qquad \Sigma_{\mathrm{corr},\sigma} =\Sigma_{\mathrm{full},\sigma}-U\bar n_{\bar\sigma},

then

Gσ1(z)=z+μH,σΣcorr,σ(z).G_\sigma^{-1}(z)=z+\mu_{\mathrm H,\sigma} -\Sigma_{\mathrm{corr},\sigma}(z).

At paramagnetic half filling, μH,σ=0\mu_{\mathrm H,\sigma}=0 and

Σfull(z)=U2+U24z,Σcorr(z)=U24z.\Sigma_{\mathrm{full}}(z)=\frac U2+\frac{U^2}{4z}, \qquad \Sigma_{\mathrm{corr}}(z)=\frac{U^2}{4z}.

Thus the often-quoted U2/(4z)U^2/(4z) is the Hartree-subtracted self-energy, not the full self-energy paired with G01=z+U/2G_0^{-1}=z+U/2. The particle–hole shift μU/2\mu-U/2 coincides with the Hartree shift only at paramagnetic half filling.

For a Hubbard dispersion εk\varepsilon_{\mathbf k} and frequencies measured relative to the chemical potential, define

Mm,σ(k)=dω2πωmAσ(k,ω).M_{m,\sigma}(\mathbf k) =\int_{-\infty}^{\infty}\frac{\mathrm d\omega}{2\pi}\, \omega^m A_\sigma(\mathbf k,\omega).

The exact zeroth and first moments are

M0,σ(k)=1,M1,σ(k)=εkμ+Unσˉ.M_{0,\sigma}(\mathbf k)=1, \qquad M_{1,\sigma}(\mathbf k) =\varepsilon_{\mathbf k}-\mu+U\langle n_{\bar\sigma}\rangle.

For the atom, the two delta functions give M1,σ=μ+UnˉσˉM_{1,\sigma}=-\mu+U\bar n_{\bar\sigma} directly. At symmetric half filling this vanishes because equal weights lie at opposite energies. Harris and Lange developed a band-resolved moment expansion and used its zeroth moments to show the hopping-induced transfer of integrated weight away from exactly one electron per site Harris and Lange 1967, §§ IV–V, printed pp. 301–308. At exactly one electron per site and zero temperature, the lower band is filled and the upper band empty, so their total weights retain the atomic counting even though their shapes and moments change.

These identities are stronger checks than a visually plausible two-peak spectrum. A numerical or experimental decomposition should reproduce the total zeroth moment and, after its energy convention is fixed, the first moment before lower- and upper-band weights are interpreted.

A controlled lattice bridge: paramagnetic DMFT

Section titled “A controlled lattice bridge: paramagnetic DMFT”

The atomic solution alone does not describe a metal or a transition. A controlled bridge is supplied by the homogeneous one-band Hubbard model in the large-coordination limit. With hopping scaled so that jtij2\sum_j|t_{ij}|^2 remains finite, dynamical mean-field theory (DMFT) makes the self-energy local while preserving its frequency dependence:

GR(k,ω)=1ω+μεkΣR(ω),G^R(\mathbf k,\omega) =\frac{1}{\omega+\mu-\varepsilon_{\mathbf k}-\Sigma^R(\omega)}, Gloc(z)=dερ0(ε)z+μεΣ(z).G_{\mathrm{loc}}(z) =\int\mathrm d\varepsilon\, \frac{\rho_0(\varepsilon)}{z+\mu-\varepsilon-\Sigma(z)}.

The self-consistent impurity bath vanishes in the zero-hopping limit, recovering the exact atom Georges et al. 1996, §§ II.C and III, printed pp. 20–27. At finite hopping it allows repeated motion into and out of the local environment, so the atomic poles evolve into incoherent Hubbard bands and—on a metallic branch—a low-energy coherent resonance Georges et al. 1996, § VII.C.1, printed pp. 62–64. This is a controlled large-coordination statement, not a claim that a local self-energy is exact in two dimensions Georges et al. 1996, § IX.A, printed pp. 107–108.

At zero temperature a Fermi-liquid solution has the low-frequency expansion

ΣR(ω,0)=ΣR(0,0)+(1Z1)ωiCω2+,C>0,\Sigma^R(\omega,0) =\Sigma^R(0,0)+(1-Z^{-1})\omega -iC\omega^2+\cdots, \qquad C>0,

and hence a coherent dispersion

EkZ[εk+ReΣR(0,0)μ].E_{\mathbf k}\simeq Z\left[\varepsilon_{\mathbf k} +\operatorname{Re}\Sigma^R(0,0)-\mu\right].

The coherent scale is therefore of order ZWZW, whereas the Hubbard bands remain at energies of order UU. In the symmetry-restricted, half-filled paramagnetic DMFT problem at T=0T=0, the metallic branch approaches its endpoint with Z0Z\to0: the coherent scale collapses while the sum rule forces the missing pole weight into incoherent sectors Zhang, Rozenberg, and Kotliar 1993, pp. 1666–1668. This bounded result reconstructs the advertised coherent-to-incoherent evolution.

It is not a universal transition law. For finite 0<T<Tcrit0<T<T_{\mathrm{crit}}, the paramagnetic DMFT metal and insulator coexist and the physical first-order transition is selected by thermodynamics, not by continuing a single Z0Z\to0 curve Georges et al. 1996, § VII.D.1, printed pp. 65–66. The zero-temperature endpoint instead occurs at Uc=Uc2U_c=U_{c2} Georges et al. 1996, § VII.E, printed pp. 70–71. On a bipartite lattice, antiferromagnetic order can preempt the paramagnetic transition Georges et al. 1996, § VII.D.3, printed pp. 69–70. In finite dimensions, nonlocal correlations can also create momentum-selective loss of coherence that a single local ZZ cannot represent. The DMFT mapping page develops the fixed-point construction and these control tests.

Hole doping: exact counting and dynamical transfer

Section titled “Hole doping: exact counting and dynamical transfer”

Let xx be the hole density, so the electron density is n=1xn=1-x. In the zero-hopping, large-UU, zero-temperature counting limit, double occupancy is absent: a fraction xx of sites is empty and a fraction 1x1-x is singly occupied. For 0<x<10<x<1, this strict atomic statement is a canonical mixture—or an ensemble at the μ=0\mu=0 degeneracy between empty and singly occupied sites—not a nondegenerate grand-canonical atomic ground state at generic μ\mu.

The spin-summed weights per site are

ProcessCounting reasonAtomic weight
Electron removalOne removable electron on each singly occupied site1x1-x
Low-energy electron additionEither spin can fill each empty site2x2x
Upper-Hubbard-band additionOnly the opposite spin can join each singly occupied site1x1-x

The total removal weight is 1x1-x, and the full addition weight is 2(1x)=1+x2-(1-x)=1+x. Their sum is the exact spin-summed zeroth moment, 22. The factor 2x2x is state counting, not a quasiparticle residue Eskes, Meinders, and Sawatzky 1991, printed pp. 1035–1036.

To state the finite-hopping result unambiguously, choose a separated-band window

tΩLU,t\ll\Omega_L\ll U,

and define the zero-temperature low-energy addition weight using the Lehmann addition sector,

Waddlow(ΩL)=1Li,σ0ΩLdω2πAiσ+(ω).W_{\mathrm{add}}^{\mathrm{low}}(\Omega_L) =\frac1L\sum_{i,\sigma} \int_0^{\Omega_L}\frac{\mathrm d\omega}{2\pi}\, A^+_{i\sigma}(\omega).

The A+A^+ label matters at the μ=0\mu=0 atomic degeneracy because a sign-of-frequency split cannot distinguish addition from removal there. At nonzero temperature, photoemission and inverse-photoemission weights are separated by the equilibrium factors f(ω)A(ω)f(\omega)A(\omega) and [1f(ω)]A(ω)[1-f(\omega)]A(\omega), subject to probe matrix elements.

The strong-coupling transformation gives more than a definition of the transferred weight. Use the prerequisite convention Ht=ijσtijciσcjσH_t=-\sum_{ij\sigma}t_{ij}c_{i\sigma}^\dagger c_{j\sigma}, let PP project onto states without doublons, and transform the creation operator together with the Hamiltonian:

ciσ,eff=PeSciσeSP=P(ciσ+[S(1),ciσ]+)P.c_{i\sigma,\mathrm{eff}}^\dagger =Pe^S c_{i\sigma}^\dagger e^{-S}P =P\left(c_{i\sigma}^\dagger+[S^{(1)},c_{i\sigma}^\dagger]+\cdots\right)P.

Using completeness within the low-energy sector converts the integrated addition spectrum into an equal-time expectation value of these matched operators. Let ΨP|\Psi_P\rangle be the normalized projected representative of the transformed low-energy ground state, with PΨP=ΨPP|\Psi_P\rangle=|\Psi_P\rangle, and define

εkinP=1LΨPHtΨP<0.\varepsilon_{\mathrm{kin}}^P =\frac1L\langle\Psi_P|H_t|\Psi_P\rangle<0.

then, through first order in t/Ut/U,

Waddlow=2x2εkinPU+O ⁣(t2U2)2x+α,W_{\mathrm{add}}^{\mathrm{low}} =2x-\frac{2\varepsilon_{\mathrm{kin}}^P}{U} +O\!\left(\frac{t^2}{U^2}\right) \equiv2x+\alpha, α=2εkinPU+O ⁣(t2U2).\alpha=-\frac{2\varepsilon_{\mathrm{kin}}^P}{U} +O\!\left(\frac{t^2}{U^2}\right).

The zeroth-order term counts the two ways to fill an empty site. The commutator term describes an added electron that makes a virtual doublon and then reaches a neighboring hole; it is typically O(xt/U)O(xt/U) and positive for the usual kinetic-energy convention. This derivation and its local form appear in Randeria et al. 2005, printed p. 137001-2, eqs. (5)–(6), especially eq. (6). It is the operator-level version of dynamical spectral-weight transfer developed on the superexchange page.

Because the full addition weight remains 1+x1+x, the separated upper-band addition weight is

WaddUHB=1xα+O ⁣(t2U2).W_{\mathrm{add}}^{\mathrm{UHB}} =1-x-\alpha +O\!\left(\frac{t^2}{U^2}\right).

Thus weight moves downward without violating the total sum rule. The leading correction vanishes in the no-hole projected state because PHtPPH_tP cannot move an electron at exactly one electron per site. At finite temperature or finite doublon density, the simple three-process atomic partition acquires thermally or virtually activated contributions even though the total addition, removal, and full zeroth-moment identities remain exact.

The diagram summarizes the exact counting and the controlled leading transfer. Inspect the arrow rather than the sketched line shapes: α\alpha moves addition weight between separated sectors, while ZkFZ_{\mathbf k_F} is defined in a different, momentum-resolved inset.

The spin-summed total spectral weight remains two: hole doping creates low-energy addition weight two x, finite hopping transfers alpha from upper-band to low-energy addition, and a separate quasiparticle-pole inset keeps Z at the Fermi momentum distinct from either integrated weight.

Spectral-weight accounting for the repulsive one-band model at zero temperature. Atomic hole counting gives removal 1x1-x, low-energy addition 2x2x, and upper-band addition 1x1-x; in a separated-band strong-coupling expansion, α\alpha transfers addition weight downward while the spin-summed total stays 22. The line shapes and transfer arrow are schematic, α\alpha depends on the cutoff and projected kinetic energy, and neither 2x+α2x+\alpha nor a symmetric low-energy integral equals the momentum-resolved pole residue ZkFZ_{\mathbf k_F}.

Open the diagram at full size or download its semantic data.

RegimeTotal removalLow-energy additionUpper-band additionStatus
Half-filled atom110011Exact at t=0t=0, T=0T=0
Hole-doped atom1x1-x2x2x1x1-xExact no-doublon counting at t=0t=0, T=0T=0
Separated bands, finite hopping1x1-x in total2x+α2x+\alpha1xα1-x-\alphaAddition weights through the declared t/Ut/U order and cutoff
Fermi-level quasiparticleNot an integrated sectorNot ZkFZ_{\mathbf k_F}Not applicablePole residue only when a quasiparticle exists

When Hubbard bands overlap, no unique ΩL\Omega_L separates them and α\alpha becomes a window-dependent descriptor rather than a sharply band-resolved quantity. The cutoff dependence follows from the definition itself. Exact-diagonalization results show growth beyond the static 2x2x count and sensitivity to doping and hybridization; temperature can enter indirectly when it changes that hybridization Meinders, Eskes, and Sawatzky 1993, printed pp. 3917–3921 and 3925–3926.

Integrated weight is not quasiparticle residue

Section titled “Integrated weight is not quasiparticle residue”

Let LL be the number of sites and

Aσ(ω)=1LkAσ(k,ω)A_\sigma(\omega) =\frac1L\sum_{\mathbf k}A_\sigma(\mathbf k,\omega)

the local spectrum per site and per spin. It obeys dωAσ/(2π)=1\int\mathrm d\omega\,A_\sigma/(2\pi)=1. A symmetric low-energy integral,

Wlow(Ω)=ΩΩdω2πσAσ(ω),W_{\mathrm{low}}(\Omega) =\int_{-\Omega}^{\Omega}\frac{\mathrm d\omega}{2\pi} \sum_\sigma A_\sigma(\omega),

includes both addition and removal weight, and every coherent and incoherent contribution inside the chosen window. It is therefore not the same object as the addition-only WaddlowW_{\mathrm{add}}^{\mathrm{low}} above.

By contrast,

ZkF=[1ωReΣR(kF,ω)ω=0]1Z_{\mathbf k_F} =\left[ 1-\left. \partial_\omega\operatorname{Re}\Sigma^R(\mathbf k_F,\omega) \right|_{\omega=0} \right]^{-1}

is the residue of a Fermi-level quasiparticle pole when the derivative exists and a Fermi-liquid description applies. It is momentum resolved in general. Consequently, neither 2x2x, 2x+α2x+\alpha, nor Wlow(Ω)W_{\mathrm{low}}(\Omega) is ZkFZ_{\mathbf k_F}. The atomic Mott insulator makes the distinction stark: its Hubbard bands have finite integrated weights, while its Fermi-level self-energy pole leaves no quasiparticle residue.

An experimental peak area is a third object. Photoemission measures a matrix-element-weighted occupied intensity, schematically IM2fAI\propto |M|^2 fA, convolved with energy and momentum resolution and mixed with backgrounds. Inverse photoemission or absorption accesses addition weight with different matrix elements. A pole residue can be inferred only after those effects and the incoherent background are controlled.

The chapter’s reduction map keeps the full Hubbard spectrum distinct from a projected low-energy Hamiltonian. Its validity table gives the stop conditions for spectral and mechanism claims.

Proposed conclusionMinimum checksWhat survives if a check fails
“The sum rule is satisfied”Positivity, declared A=2ImGRA=-2\operatorname{Im}G^R, M0M_0, M1M_1, chemical-potential convention, and sufficient frequency rangeOnly a windowed intensity, not a normalized spectral function
“Weight moved between Hubbard sectors”A stable separating window, cutoff variation, temperature control, and matched orbital or probe matrix elementsA redistribution within the measured window
“A coherent quasiparticle exists”Momentum-resolved pole or Fermi-liquid self-energy, resolution below the coherent scale, and a controlled backgroundA narrow feature, not a certified ZZ
“The state is a Mott insulator”Charge gap or zero-temperature incompressibility, symmetry-restored persistence, atomic continuity or self-energy evidence, and exclusion of leading alternativesAn insulating or two-feature spectrum with unresolved mechanism

Analytic continuation is an inverse problem: broad Hubbard bands are usually more stable than a narrow low-energy resonance, while peak splitting and small gaps can depend on the prior, regularization, and covariance. A finite cluster adds discrete levels and limited momentum sectors; size, shape, boundary conditions, broadening, periodization, and interpolation must be varied before a thermodynamic or momentum-selective conclusion is drawn. Magnetic order can fold bands and open a Slater gap, and disorder, hybridization, phonon satellites, multiplets, and charge-transfer physics can all produce more than one spectral feature.

The strongest conclusion from two broad peaks alone is therefore modest: the measured or computed electron spectrum contains two resolved energy regions in the stated channel and resolution. Mottness requires the additional thermodynamic, symmetry, operator, and mechanism tests developed on the insulator-comparison page.

Separating addition and removal only by the sign of frequency. This works away from zero-cost degeneracies. At the μ=0\mu=0 atomic crossing, retain the Lehmann-sector label or take a controlled finite-hopping limit.

Calling 2x+α2x+\alpha a quasiparticle weight. It is an addition-sector integral over all coherent and incoherent states below a cutoff. ZkFZ_{\mathbf k_F} is the residue of one momentum-resolved pole.

Treating quasiparticle collapse as universal. The continuous Z0Z\to0 statement above belongs to the zero-temperature paramagnetic large-coordination DMFT branch. A finite-temperature first-order transition, magnetic preemption, or momentum-selective finite-dimensional transition can look different.

Identifying Mottness from two peaks. Two structures do not establish incompressibility, a self-energy pole, or continuity to the interaction-blocked atomic limit. Test symmetry breaking, disorder, orbital character, and probe matrix elements.

For the atomic Hamiltonian, derive the general Green function and self-energy, verify M0M_0 and M1M_1, and specialize to paramagnetic half filling in both self-energy conventions.

Solution

When nσˉ=0n_{\bar\sigma}=0, the propagator is (z+μ)1(z+\mu)^{-1}; when nσˉ=1n_{\bar\sigma}=1, it is (z+μU)1(z+\mu-U)^{-1}. Averaging the two sectors gives

Gσ(z)=1nˉσˉz+μ+nˉσˉz+μU.G_\sigma(z) =\frac{1-\bar n_{\bar\sigma}}{z+\mu} +\frac{\bar n_{\bar\sigma}}{z+\mu-U}.

Using G0,σ1=z+μG_{0,\sigma}^{-1}=z+\mu in Dyson’s equation yields

Σfull,σ(z)=Unˉσˉ+U2nˉσˉ(1nˉσˉ)z+μU(1nˉσˉ).\Sigma_{\mathrm{full},\sigma}(z) =U\bar n_{\bar\sigma} +\frac{U^2\bar n_{\bar\sigma}(1-\bar n_{\bar\sigma})} {z+\mu-U(1-\bar n_{\bar\sigma})}.

The residues sum to one, so M0,σ=1M_{0,\sigma}=1. Weighting the poles by their energies gives

M1,σ=(1nˉσˉ)(μ)+nˉσˉ(Uμ)=μ+Unˉσˉ.M_{1,\sigma} =(1-\bar n_{\bar\sigma})(-\mu) +\bar n_{\bar\sigma}(U-\mu) =-\mu+U\bar n_{\bar\sigma}.

At half filling, μ=U/2\mu=U/2 and nˉσˉ=1/2\bar n_{\bar\sigma}=1/2, so M1=0M_1=0 and

Σfull(z)=U2+U24z.\Sigma_{\mathrm{full}}(z)=\frac U2+\frac{U^2}{4z}.

Subtracting the Hartree term and shifting μH=0\mu_{\mathrm H}=0 gives Σcorr(z)=U2/(4z)\Sigma_{\mathrm{corr}}(z)=U^2/(4z). Both conventions produce the same GG.

In the atomic hole-doped counting limit, derive the three spin-summed weights and the spin-summed first moment. Which part of the answer requires spin balance?

Solution

There are 1x1-x singly occupied sites, so removal has weight 1x1-x. There are xx empty sites and either spin can be added, giving low-energy addition weight 2x2x. Each singly occupied site accepts only the opposite spin at energy UU, giving upper-band addition weight 1x1-x. The sum is

(1x)+2x+(1x)=2.(1-x)+2x+(1-x)=2.

At the μ=0\mu=0 atomic crossing, the lower pole combines removal and low-energy addition and has spin-summed weight 1+x1+x; the upper pole has weight 1x1-x. Hence

M1,+M1,=U(1x).M_{1,\uparrow}+M_{1,\downarrow}=U(1-x).

With a generic algebraic energy origin the result is U(1x)2μU(1-x)-2\mu. This spin-summed identity does not require spin balance. Only the further statement M1,=M1,=U(1x)/2μM_{1,\uparrow}=M_{1,\downarrow}=U(1-x)/2-\mu assumes a paramagnetic or spin-balanced mixture.

Use the matched operator ciσ,effc_{i\sigma,\mathrm{eff}}^\dagger to explain why the low-energy addition weight is 2x2x at zeroth order and why its first correction is controlled by the projected kinetic energy. What consistency bounds must α\alpha obey?

Solution

Insert a complete set of no-doublon final states into the integrated Lehmann sum. This replaces the sum over final states by the equal-time product of matched projected operators. The leading operator PciσPPc_{i\sigma}^\dagger P can act only on an empty site, of which there are xx per site, and it has two spin choices. Its spin-summed contribution is therefore 2x2x.

The cross term between PciσPPc_{i\sigma}^\dagger P and P[S(1),ciσ]PP[S^{(1)},c_{i\sigma}^\dagger]P contains one virtual doublon denominator 1/U1/U and one projected hop. Summing sites and spins gives

Waddlow=2x2εkinPU+O(t2/U2).W_{\mathrm{add}}^{\mathrm{low}} =2x-\frac{2\varepsilon_{\mathrm{kin}}^P}{U} +O(t^2/U^2).

For the stated hopping convention, εkinP<0\varepsilon_{\mathrm{kin}}^P<0, so α=2εkinP/U+\alpha=-2\varepsilon_{\mathrm{kin}}^P/U+\cdots is positive. A separated-band decomposition also requires 2x+α02x+\alpha\ge0 and 1xα01-x-\alpha\ge0. Failure of those bounds signals a sign error, an uncontrolled truncation, or a cutoff that no longer isolates the two sectors.

A finite cluster produces two broad peaks after analytic continuation. The low-energy integral changes by 20 percent when the cutoff is moved within the apparent valley; antiferromagnetic folding is present, and the narrow central feature is comparable to the continuation resolution. State what can and cannot be concluded.

Solution

One may report two resolved spectral regions for that cluster, continuation, temperature, and broadening. The cutoff sensitivity prevents a precise band-resolved α\alpha; magnetic folding leaves a Slater contribution viable; and the unresolved central feature cannot establish a quasiparticle pole or ZZ. A Mott attribution additionally requires size and cluster-shape convergence, moment and normalization checks, symmetry-restored behavior, a charge-gap or compressibility test, and comparison with hybridization, disorder, orbital, and matrix-element alternatives.

The superexchange and t–J page derives the transformed operators that generate the leading low-energy transfer. The DMFT mapping supplies the controlled local lattice construction, and the insulator comparison tests whether the resulting spectrum supports a Mott, Slater, band, charge-transfer, or Anderson mechanism.

  • H. Eskes, M. B. J. Meinders, and G. A. Sawatzky, “Anomalous Transfer of Spectral Weight in Doped Strongly Correlated Systems,” Physical Review Letters 67 (1991) 1035–1038, doi:10.1103/PhysRevLett.67.1035.
  • Antoine Georges, Gabriel Kotliar, Werner Krauth, and Marcelo J. Rozenberg, “Dynamical Mean-Field Theory of Strongly Correlated Fermion Systems and the Limit of Infinite Dimensions,” Reviews of Modern Physics 68 (1996) 13–125, doi:10.1103/RevModPhys.68.13; Open PDF.
  • A. Brooks Harris and Robert V. Lange, “Single-Particle Excitations in Narrow Energy Bands,” Physical Review 157 (1967) 295–314, doi:10.1103/PhysRev.157.295.
  • M. B. J. Meinders, H. Eskes, and G. A. Sawatzky, “Spectral-Weight Transfer: Breakdown of Low-Energy-Scale Sum Rules in Correlated Systems,” Physical Review B 48 (1993) 3916–3926, doi:10.1103/PhysRevB.48.3916.
  • Mohit Randeria, Rajdeep Sensarma, Nandini Trivedi, and Fu-Chun Zhang, “Particle-Hole Asymmetry in Doped Mott Insulators: Implications for Tunneling and Photoemission Spectroscopies,” Physical Review Letters 95 (2005) 137001, doi:10.1103/PhysRevLett.95.137001; Official PDF; Open preprint PDF.
  • X. Y. Zhang, Marcelo J. Rozenberg, and Gabriel Kotliar, “Mott Transition in the d=d=\infty Hubbard Model at Zero Temperature,” Physical Review Letters 70 (1993) 1666–1669, doi:10.1103/PhysRevLett.70.1666.