Skip to content

DMFT Impurity Mapping and Self-Consistency

Dynamical mean-field theory (DMFT) replaces a lattice problem by a quantum impurity embedded in a self-consistent bath. The mapping becomes controlled when coordination tends to infinity with hopping scaled as tijt/zt_{ij}\sim t^*/\sqrt z, making the self-energy local while retaining its full frequency dependence.

Required background. Use the Hubbard model and its limits. Helpful background. Conserving self-energy constructions and quantum impurity models provide useful context.

In the large-coordination limit, nonlocal self-energy diagrams vanish after the hopping rescaling, so Σij(z)=δijΣ(z)\Sigma_{ij}(z)=\delta_{ij}\Sigma(z). The local lattice Green function is

Gloc(z)=dϵρ0(ϵ)z+μϵΣ(z).G_{\mathrm{loc}}(z)=\int\mathrm d\epsilon\, \frac{\rho_0(\epsilon)}{z+\mu-\epsilon-\Sigma(z)}.

Define a Weiss field G0\mathcal G_0 through

G01(z)=Gloc1(z)+Σ(z).\mathcal G_0^{-1}(z)=G_{\mathrm{loc}}^{-1}(z)+\Sigma(z).

The corresponding impurity action is

Simp=σ0βdτdτcσ(τ)G01(ττ)cσ(τ)+U0βdτnn.S_{\mathrm{imp}}=-\sum_\sigma\int_0^\beta\mathrm d\tau\,\mathrm d\tau' c_\sigma^\dagger(\tau)\mathcal G_0^{-1}(\tau-\tau')c_\sigma(\tau') +U\int_0^\beta\mathrm d\tau\,n_\uparrow n_\downarrow.

Solving the impurity produces GimpG_{\mathrm{imp}} and Σimp=G01Gimp1\Sigma_{\mathrm{imp}}=\mathcal G_0^{-1}-G_{\mathrm{imp}}^{-1}. Self-consistency requires Gimp=GlocG_{\mathrm{imp}}=G_{\mathrm{loc}} and Σ=Σimp\Sigma=\Sigma_{\mathrm{imp}}. This cavity construction and its infinite-dimensional proof are developed in Georges et al. 1996, §§ II–III.

For a semicircular noninteracting density of states of half-bandwidth 2t2t^*,

Δ(z)=t2Gloc(z),G01(z)=z+μΔ(z).\Delta(z)=t^{*2}G_{\mathrm{loc}}(z), \qquad \mathcal G_0^{-1}(z)=z+\mu-\Delta(z).

At U=0U=0, G=1/(zt2G)G=1/(z-t^{*2}G), so

G(z)=zz24t22t2,G(z)=\frac{z-\sqrt{z^2-4t^{*2}}}{2t^{*2}},

where the branch is fixed by G(z)1/zG(z)\sim1/z. This exactly reproduces the semicircular density of states and checks normalization. At finite UU, metallic and insulating fixed points can coexist over a temperature-dependent region, but locating them requires a converged impurity solution and a declared iteration protocol.

Track residuals in GG, Σ\Sigma, and observables, initialize both metallic and insulating branches, and verify causality and spectral moments. Bath discretization, Monte Carlo noise, impurity truncation, and analytic continuation are separate errors. DMFT does not become controlled in two dimensions merely because UU is large; short-range spatial correlations then motivate the cluster extensions.

Select the physical square-root branch in the noninteracting Bethe equation.

Solution

For large zz, z24t2=z2t2/z+\sqrt{z^2-4t^{*2}}=z-2t^{*2}/z+\cdots on the analytic branch outside the band. The minus sign gives G1/zG\sim1/z; the plus sign grows as z/t2z/t^{*2} and violates the spectral normalization.

  • 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, §§ II–III, doi:10.1103/RevModPhys.68.13.