Gutzwiller wave function - davidar/scholarpedia GitHub Wiki

The Gutzwiller wave function (GWF) is a variational many-particle wave function which describes interacting particles on a rigid lattice. It provides an illustrative example for Landau's Fermi-liquid theory. Based on the GWF, the dispersion of quasi-particle excitations in the Fermi liquid state of transition metals can be calculated and compared to angle-resolved photo-emission (ARPES) experiments.

Table of Contents Definition Toy model description of the hydrogen molecule Gutzwiller correlator Expectation values and variational ground-state energy Evaluation Exact results in one spatial dimension Average double occupancy Momentum distribution Variational energy for the single-band Hubbard model Spin correlations and Haldane-Shastry model Exact results in the limit of infinite spatial dimensions Simplifications Variational ground-state energy Landau-Gutzwiller quasi-particles Relation to other methods Gutzwiller Approximation (GA) and Brinkman-Rice transition Kotliar-Ruckenstein slave-boson mean-field theory Applications Liquid Helium-3 Band structure of nickel Further applications Correlated superconductors Atoms in optical lattices Appendix: historical note References Further reading See also

In order to illustrate the basic idea behind the GWF, consider the simplest interacting electron model, the hydrogen molecule, <math>H_2\ ,</math> where the two sites represent the protons. In a toy model description, the electrons can occupy only the <math>1s</math> orbital of either proton with spin <math>\sigma=\uparrow,\downarrow\ .</math> The four possible configurations are shown in <figref>Gutzwiller_Wave_Function_GWF-hydrogen.jpg</figref>. Naturally, the basis set is far too small for a quantitative description of the hydrogen molecule.

The electrons can tunnel between the protons with amplitude <math>-t</math> (<math>t>0</math>), and the ground state <math>|\Psi\rangle</math> is a linear combination of the four configurations. In the Hund-Mulliken molecular-orbital (MO) description, all four configurations have the same probability. The MO ground state <math>|\Psi_0\rangle</math> is incorrect because the configurations <math>|{\rm iii}\rangle</math> and <math>|{\rm iv}\rangle</math> correspond to a negatively charged <math>H^-</math> atom next to a bare proton (<math>H^+</math>). Such configurations where the <math>1s</math>-orbital is doubly occupied are energetically unfavorable because the We added electrons on the same site repel each other due to their mutual Coulomb interaction (generally called the Hubbard <math>U\ ,</math> <math>U>0</math>). Therefore, the weight of doubly occupied orbitals is reduced in the ground state <math>|\Psi\rangle\ .</math>

If one starts with some simple single-particle product wave function <math>|\Psi\rangle_0</math> like the MO function, a better wave function <math>|\Psi\rangle</math> can be obtained by reducing the weight of the configurations containing double occupancies, such as <math>|{\rm iii}\rangle</math> and <math>|{\rm iv}\rangle</math> in <figref>Gutzwiller_Wave_Function_GWF-hydrogen.jpg</figref>. The operator that counts the number of doubly occupied sites is <math>\hat{D}</math> and the ground state can be written in the form of a GWF (M.C. Gutzwiller, 1963),

<math>\label{eq:def-Psi-hydro}

|\Psi\rangle = g^{\hat{D}}|\Psi_0\rangle \; , </math>

where <math>g=g(U)</math> is a function of the electrons' intra-orbital Coulomb interaction. The toy model Hamiltonian is the two-site Hubbard model,

<math>\label{eq:two-siteH}

\hat{H} = \hat{T}+U\hat{D} = (-t)\sum_{\sigma=\uparrow,\downarrow} \left( \hat{c}_{1,\sigma}^+\hat{c}_{2,\sigma} + \hat{c}_{2,\sigma}^+\hat{c}_{1,\sigma} \right) + U \left( \hat{c}_{1,\uparrow}^+\hat{c}_{1,\uparrow}\hat{c}_{1,\downarrow}^+\hat{c}_{1,\downarrow}+ \hat{c}_{2,\uparrow}^+\hat{c}_{2,\uparrow}\hat{c}_{2,\downarrow}^+\hat{c}_{2,\downarrow} \right) \; . </math>

The MO ground state is <math>|\Psi_0\rangle=(|{\rm i}\rangle -|{\rm ii}\rangle +|{\rm iii}\rangle +|{\rm iv}\rangle)/2</math> and the GWF reads

<math>\label{eq:GWF-two-site}

|\Psi\rangle = \frac{1}{2} \left[|{\rm] \;. </math>

For <math>U\to\infty\ ,</math> <math>g\to 0\ ,</math> and <math>|\Psi\rangle</math> reduces to the Heitler-London valence-bond (VB) state, which is a linear combination of the configurations <math>|{\rm i}\rangle</math> and <math>|{\rm ii}\rangle</math> in <figref>Gutzwiller_Wave_Function_GWF-hydrogen.jpg</figref>. The GWF \eqref{eq:def-Psi-hydro} smoothly interpolates between the Hund-Mulliken and the Heitler-London wave functions.

For the two-site Hubbard model, the variational energy is given by

<math>\label{eq:2site-Evar}

E_{\rm var}(g) = -\frac{4tg}{1+g^2} + U \frac{g^2}{1+g^2} \; . </math>

The variational minimum of the energy is obtained from the condition <math>\partial E_{\rm var}(g)/\partial g=0\ ,</math> which is fulfilled for

<math>\label{eq:gopt-twosite}

g_{\rm opt}(U) = \frac{1}{4t}\left[-U] \; . </math>

For the two-site Hubbard model, the variational upper bound becomes the exact ground-state energy,

<math>\label{eq:groundstateenergy-twosite}

E_{\rm var}(g_{\rm opt}(U))=E_0(U)=\frac{1}{2} \left(U - \sqrt{U^2+(4t)^2\,}\right)\; . </math>

For large interaction strengths, <math>U\gg 4t\ ,</math> the Gutzwiller parameter and the ground-state energy become small, <math>g\to 2t/U\ ,</math> <math> E_0\to -4t^2/U\ ,</math> but they remain finite, of order <math>1/U\ .</math>

Crystal lattices contain a large number <math>L</math> of atoms and a similarly large number of electrons. To a first approximation, each electron can be treated as moving in an average potential. The resulting single-electron product state <math>|\Psi_0\rangle</math> is the lattice generalization of the Hund-Mulliken molecular-orbital (MO) wave function (also known as Hartree-Fock wave function). For example, <math>|\Psi_0\rangle=|\hbox{FS}\rangle</math> denotes the Fermi sea of noninteracting electrons.

The spectrum and eigenstates of the individual atoms on the lattice sites are <math>E_{\Gamma}</math> and <math>|\Gamma\rangle\ .</math> In the single-band example, there are four atomic states, <math>|1\rangle=|H^+\rangle\ ,</math> <math>|2\rangle=|H_{\uparrow}\rangle\ ,</math> <math>|3\rangle=|H_{\downarrow}\rangle\ ,</math> and <math>|4\rangle=|H^-\rangle\ ,</math> with the atomic energies <math>E_1=E_2=E_3=0</math> and <math>E_4=U\ .</math>

The Gutzwiller correlator (GC) reduces the weight of those configurations in <math>|\Psi_0\rangle</math> which are energetically unfavorable, i.e., which have a large atomic energy. The general GWF is thus defined by (J. Bünemann, F. Gebhard, and W. Weber, 1998)

<math>\label{eq:def}

|\Psi_{\rm G}\rangle = \hat{P}_{\rm G}|\Psi_0\rangle \; . </math>

The symbol <math>\hat{P}_{\rm G}</math> is an operator that modifies the single-particle product wave function <math>|\Psi_0\rangle\ .</math> The operator <math>\hat{m}_{\Gamma,\vec{R}}= |\Gamma\rangle_{\vec{R}} {}_{\vec{R}}\langle \Gamma|</math> picks out a particular atom on site <math>\vec{R}</math> and then projects the atom onto the atomic configuration <math>\Gamma\ .</math> The GC assigns a weight factor <math>\lambda_{\Gamma,\vec{R}}</math> to each atomic configuration,

<math>\label{eq:def-PG1}

\hat{P}_{\rm G} = \prod_{\vec{R}} \hat{P}_ \quad ,\quad \hat{P}_ = \prod_{\Gamma} \lambda_{\Gamma,\vec{R}}^{\hat{m}_{\Gamma,\vec{R}}} = \sum_{\Gamma} \lambda_{\Gamma,\vec{R}} \hat{m}_{\Gamma,\vec{R}} \; . </math>

The GWF \eqref{eq:def-Psi-hydro} for the one-band model is recovered for <math>\lambda_{1}=\lambda_{2}=\lambda_{3}=1</math> and <math>\lambda_{4}=g</math> for the four hydrogen atomic states <math>\bigl\{ H^+,H_{\uparrow},H_{\downarrow}, H^-\bigr\}\ .</math>

The Gutzwiller correlator suppresses charge fluctuations which are too large in <math>|\Psi_0\rangle\ .</math> For example, a Fermi-gas description of the <math>3d</math>-electrons of nickel employs the Fermi-sea ground state. It predicts a certain probability to detect five-fold ionized nickel ions <math>\hbox{Ni}^{5+}</math> in a nickel crystal. The energy <math>E_</math> of the corresponding atomic configuration is more than one hundred eV. Since the ground state has the lowest possible total energy, the probability to find five-fold ionized nickel ions <math>\hbox{Ni}^{5+}</math> in the ground state of a nickel crystal must be exponentially small. This is ensured by the Gutzwiller correlator. In the GWF, only atomic configurations with approximately 10 electrons in the <math>3d</math>-shell and the <math>4s</math>-shell have non-negligible probability.

The Gutzwiller correlator belongs to the class of Jastrow-Feenberg correlators which are used, e.g., for the investigation of superfluid Helium-4.

Physical quantities such as the magnetization <math>M_z</math> are determined by expectation values of corresponding quantum-mechanical operators <math>\hat{M_z}\ .</math> Quite generally, the Gutzwiller variational expressions for expectation values of observables are

<math>\label{eq:expectationvalues}

A_{\rm var}=\langle \hat{A}\rangle_{\rm G} = \frac{\langle \Psi_0 | \hat{P}_{\rm G}^+ \hat{A} \hat{P}_{\rm G}|\Psi_0 \rangle} {\langle \Psi_0 | \hat{P}_{\rm G}^+ \hat{P}_{\rm G}|\Psi_0 \rangle} \; . </math>

The variational parameters in the GWF are obtained as follows: Starting from <math>|\Psi_0\rangle</math> and <math>\hat{P}_{\rm G}\ ,</math> the variational energy <math>E_{\rm var}=\langle\hat{H}\rangle</math> for a given model Hamiltonian, e.g., for the Hubbard model (read a historical note on its invention), must be calculated. Then, <math>E_{\rm var}</math> must be minimized with respect to all variational parameters to find the optimal variational energy, <math>E_{\rm var}^{\rm opt}\ .</math>

The evaluation of expectation values within the GWF poses a many-particle problem which is unsolvable in general. The numerical evaluation of expectation values is possible on finite lattice with the help of the Variational Monte Carlo (VMC) method as introduced by C. Gros et al., 1987, and by H. Yokoyama and H. Shiba, 1987; for a review, see Edegger et al., 2007.

As in standard Feynman-Dyson perturbation theory, the analytical evaluation of expectation values starts from a series expansion around the noninteracting limit, <math>\lambda_{\vec{R}}=1\ ,</math> and the individual orders are expressed in terms of diagrams. In contrast to the standard calculation of Green functions for interacting electron systems, the GWF permits various choices for the bare vertex. This flexibility of the GWF allows its exact evaluation in one spatial dimension and in the limit of infinite dimensions. In his three articles, Gutzwiller provided the first realistic picture of the ferromagnetism of nickel and its alloys with copper; further developments have been due to other authors.

For electrons on a ring, expectation values for the Gutzwiller-correlated Fermi sea can be evaluated for all electron densities <math>n=n_{\uparrow}+n_{\downarrow}\ ,</math> magnetizations <math>m=(n_{\uparrow}-n_{\downarrow})\ ,</math> and interaction parameters <math>g</math> without making further approximations. When expectation values are expanded in terms of <math>(g^2-1)</math> and particle-hole symmetry is used, all coefficients of the series expansion can be determined (W. Metzner and D. Vollhardt, 1988) (F. Gebhard and D. Vollhardt, 1988).

The average number of doubly occupied sites is defined by <math>\overline{d} = \langle \hat{D}/L\rangle_{\rm G} = (1/L) \sum_{\ell=1}^L \langle \hat{n}_{\ell,\uparrow}\hat{n}_{\ell,\downarrow} \rangle_{\rm G}\ .</math> The result for the Gutzwiller-correlated Fermi sea in one dimension for all interaction parameters <math>g\geq 0\ ,</math> electron densities <math>0\leq n\leq 1\ ,</math> and magnetizations <math>0\leq m\leq n</math> is (M. Kollar and D. Vollhardt, 2002)

<math>\label{eq:overlined}

\overline{d}(g,n,m)= \frac{g^2}{2(1-g^2)^2}\left[-(1-g^2)(n-m) +\ln\left(\frac{1-(1-g^2)m}{1-(1-g^2)n}\right)\right] \; . </math>

For <math>g=1</math> (noninteracting limit), the result reproduces the Hartree-Fock value, <math>\overline{d}(g=1,n,m) = (n^2-m^2)/4\ .</math> For <math>g\to 0</math> (strong-coupling limit), the double occupancy vanishes, <math>\overline{d}(g=0,n,m) = 0\ .</math>

The momentum distribution <math>n_{k,\sigma}= \langle \hat{n}_{k,\sigma}\rangle_{\rm G}</math> is <math>2\pi</math>-periodic and inversion symmetric for a symmetric dispersion relation <math>\varepsilon(-k)= \varepsilon(k)</math> of the underlying Hubbard model.

The exact formulae for the Gutzwiller-correlated Fermi sea are rather involved. For the paramagnetic half-filled Fermi sea, <math>n=1</math> and <math>m=0\ ,</math> and with <math>0\leq \tilde{k}=2|k|/\pi\leq 2\ ,</math> <math>G=1-g^2\ ,</math> they can be written as

<math>\label{eq:momdistribution}

n_<table k\leq \pi/2,\sigma}(g)=\frac{g^2+4g+1}{2(1+g)^2}></table>+ \frac{g^2}{(1+g)^2} \frac{4}{\pi\sqrt{(2-G)^2-(\tilde{k}G)^2}} K\left[\frac{G\sqrt{1-\tilde{k}^2}}{\sqrt{(2-G)^2-(\tilde{k}G)^2}}\right] \; , </math>

where

<math>\label{eq:defK}

K(k)=\int\limits_{0}^{\pi/2}{\rm d}\phi [1-k^2\sin^2(\phi)]^{-1/2} </math>

is the complete elliptic integral of the first kind. Moreover, <math>n_{\pi/2<|k|\leq \pi,\sigma}(g) = 1-n_{\pi-|k|,\sigma}(g)\ .</math>

In general, the momentum distribution is discontinuous at the Fermi wave number <math>k_{\rm F}=\pi/2\ ,</math> as is characteristic for a metal. For the half-filled paramagnetic Fermi sea (<math>n=1\ ,</math> <math>m=0</math>), the jump is given by

<math>\label{eq:jump}

q(g)=n_{k=\pi/2^-,\sigma}(g)- n_{k=\pi/2^+,\sigma}(g)= \frac{4g}{(1+g)^2} \; . </math>

The discontinuity vanishes for strong coupling, <math>g=0\ ,</math> when all electrons are localized (Brinkman-Rice insulator).

From the exact solution in one dimension via Bethe Ansatz it is known that the Hubbard model away from half band-filling provides an example of a Luttinger Liquid whose momentum distribution is continuous with a divergent slope at the Fermi wave number. The GWF does not reproduce this generic behavior of one-dimensional metals.

The variational energy provides an exact upper bound for the ground-state energy of the Hubbard model, which is exactly known for the Hubbard model from the Bethe Ansatz for the dispersion relation <math>\varepsilon(k)=-2t\cos(k)</math> (electron transfer between nearest neighbors only). For the paramagnetic case, <math>m=0\ ,</math> the comparison shows that the ground-state energy for half band filling, <math>n=1\ ,</math> deviates substantially from the exact result for large interactions, <math>U/t\to \infty</math> (W. Metzner and D. Vollhardt, 1988).

In this limit, double occupancies and empty sites are constrained to be adjacent to each other in the exact ground state, but this correlation is absent in the GWF (F. Gebhard and D. Vollhardt, 1988). To cure this problem, various extensions of the GWF have been proposed, e.g., the Baeriswyl-GWF wave function (M. Dzierzawa, D. Baeriswyl, and M. Di Stasio, 1995) and the Local-Ansatz wave functions (P. Fulde, 1995), which introduce correlations between neighboring sites. Typically, these wave functions can be evaluated analytically only in limiting cases. Otherwise, they must be treated numerically (Variational Monte Carlo) (B. Edegger, V.N. Muthukumar et al., 2007).

The <math>z</math>-component of the spin-spin correlation function is defined by <math>C^{\rm SS}(r) =

 (1/L)\sum_{\ell=1}^L \langle \hat{S}_{r+\ell}^{z}\hat{S}_{\ell}^{z}

\rangle_{\rm G}\ ,</math> where <math>\hat{S}_{\ell}^{z} =(\hat{n}_{\ell,\uparrow}- \hat{n}_{\ell,\downarrow})/2</math> is the operator for the z-component of the electron spin on site <math>\ell\ .</math>

While the correlations between double occupancies and holes are poorly described by the GWF in the strong-coupling limit, the spin-spin correlations at half band filling correctly show the characteristic behavior of Heisenberg-type spin models in one dimension. The spin-spin correlations of the Gutzwiller-projected paramagnetic Fermi sea (<math>g=0\ ,</math> <math>n=1\ ,</math> <math>m=0</math>) are given by (F. Gebhard and D. Vollhardt, 1988)

<math></math>\label{eq:spinCF}

C^{\rm SS}(r>0;g=0) = (-1)^r \fracdef} can be evaluated without further approximations. It is possible to set up a diagrammatic series expansion around the uncorrelated limit, <math>\lambda_{\Gamma,\vec{R}}=1\ .</math> In this case, not a single diagram must be calculated in the <math>Z\to \infty</math> limit. The theory remains nontrivial because the single-particle density matrix <math>P_{\vec{R},\sigma;\vec{R}',\sigma'}=\langle \hat{c}_{\vec{R},\sigma}^{+} \hat{c}_{\vec{R}',\sigma'}^{}\rangle_{\rm G}</math> and the average atomic occupancy <math>m_{\Gamma,\vec{R}} = \langle \hat{m}_{\Gamma,\vec{R}}\rangle_{\rm G}</math> are renormalized in the procedure (J. Bünemann, F. Gebhard, and W. Weber, 1998).

In the limit of infinite coordination number, the variational parameters <math>\lambda_{\Gamma,\vec{R}}</math> may be replaced by the physical expectation values <math>m_{\Gamma,\vec{R}}</math> for the occupation of an atomic configuration <math>\Gamma</math> on lattice site <math>\vec{R}\ .</math> Moreover, the local density matrix for noninteracting electrons <math>C^0_{\vec{R};\sigma,\sigma'}=P^0_{\vec{R},\sigma;\vec{R},\sigma'}</math> must obey certain constraints which can be included with the help of Lagrange parameters <math>\eta_{\vec{R};\sigma,\sigma'}\ .</math>

For a translationally invariant symmetric multi-band Hubbard model, the variational ground-state energy functional for a normalized single-particle product state <math>|\Psi_0\rangle</math> with fixed average particle density reads

<math>\label{eq:varfunctional}

E_{\rm var}\left(m_{\Gamma},\eta_{\sigma,\sigma'}, C_{\sigma,\sigma'},\left\{ |\Psi_0\rangle \right\}\right) =\langle \Psi_0 | \hat{T}^{\rm eff} | \Psi_0 \rangle + L \sum_{\Gamma} E_{\Gamma} m_{\Gamma} -L \sum_{\sigma,\sigma'}\eta_{\sigma,\sigma'} C^0_{\sigma,\sigma'}\quad , \quad \hat{T}_{\rm eff} =\sum_{\vec{k};\sigma,\sigma'} \varepsilon_{\sigma,\sigma'}^{\rm eff}(\vec{k}) \hat{c}_{\vec{k},\sigma}^+ \hat{c}_{\vec{k},\sigma'}^{} \; , </math>

where <math>\sigma=1,2,\ldots 2N</math> labels the atomic orbitals. In particular, <math>N=1,3,5</math> for atomic <math>s\ ,</math> <math>p\ ,</math> <math>d</math> shells. The minimization of the functional with respect to <math>|\Psi_0\rangle</math> shows that it is the ground state of the effective kinetic energy <math>\hat{T}^{\rm eff}</math> with the effective dispersion relation

<math>\label{eq:epseffective}

\varepsilon_{\sigma,\sigma'}^{\rm eff}(\vec{k}) = \sum_{\gamma,\gamma'} Q_{\gamma,\gamma'}^{\sigma,\sigma'} \varepsilon_{\gamma,\gamma'}^0(\vec{k}) +\eta_{\sigma,\sigma'} \; . </math>

The matrix <math>Q_{\gamma,\gamma'}^{\sigma,\sigma'}</math> is a known but, in general, complicated function of the variational parameters. The matrices <math>Q_{\gamma,\gamma'}^{\sigma,\sigma'}</math> and <math>\eta_{\sigma,\sigma'}</math> express the fact that the electron-electron interaction reduces the bandwidth of the bare bands with dispersion <math>\varepsilon_{\gamma,\gamma'}^0(\vec{k})</math> and changes their hybridization as well as their relative positions.

The minimization of the energy functional with respect to all variational parameters is a numerically demanding task for real materials because it requires the minimization of a functional with several thousands of parameters.

The Gutzwiller variational theory provides an explicit example of the Landau Fermi-liquid theory. The Gutzwiller correlator <math>\hat{P}_{\rm G}</math> continuously transforms the single-particle ground state <math>|\Psi_0\rangle</math> to the (variational) ground state for interacting particles <math>|\Psi\rangle_{\rm G}\ .</math> For example, the momentum distribution <math>n_{\vec{k},\sigma}</math> displays a jump discontinuity at the Fermi energy both in <math>|\Psi_0\rangle</math> and in <math>|\Psi\rangle_{\rm G}\ ,</math> see <figref>Gutzwiller_Wave_Function_GWF-momentum.jpg</figref> and eq. \eqref{eq:jump}. This jump represents the Fermi surface in the reciprocal (or Fourier) space where the individual single-electron states depend on the wave vector <math>\vec{k}\ .</math> The discontinuity is smaller than unity and depends on the strength of the coupling. Gutzwiller (1965) provided a relatively simple derivation of this generic result within the GWF. In contrast, Hubbard's original theory led to a band splitting into an upper and a lower Hubbard band, so that a system at half band filling would describe an insulator, irrespective of the strength of the Coulomb interaction. The Gutzwiller theory for the ground state of an interacting many-particle system provides the dominant Landau parameters <math>F_{0,1}^{s,a}\ ,</math> which determine the thermodynamics as well as the dispersion relation of the dominant hydrodynamic modes (zero and first sound). The excitation energies and the temperature must be small compared to the Fermi temperature.

In the spirit of the Landau Fermi-liquid theory, the Landau-Gutzwiller theory describes quasi-hole (quasi-particle) excitations as Gutzwiller-correlated holes (particles) in <math>|\Psi_0\rangle</math> (J. Bünemann, F. Gebhard, and R. Thul, 2003). Their dispersion relation is the same as that obtained for the effective kinetic energy <math>\hat{T}^{\rm eff}\ ,</math> eq. \eqref{eq:epseffective}. Therefore, <math>\varepsilon_{\sigma,\sigma'}(\vec{k})</math> defines the quasi-particle band structure, which can be compared to experimental data from angle-resolved photo-emission spectroscopy (ARPES).

A time-dependent version of the Gutzwiller wave function (G. Seibold et al., 2004) permits the direct calculation of the dynamic magnetic and charge correlation functions.

The general formulae in infinite dimensions considerably simplify for the Gutzwiller-correlated paramagnetic Fermi sea which is a variational ground state for the single-band Hubbard model. In a particle-hole symmetric system at half band filling, the average double occupancy <math>\overline{d}\ ,</math> the bandwidth reduction factor <math>q_{\sigma}\ ,</math> the dispersion of the quasi-particles <math>\varepsilon_{\sigma}^{\rm eff}(\vec{k})\ ,</math> and the momentum distribution <math>n_{\vec{k},\sigma}</math> are given by

<math>\label{overlinedinfty}

\overline{d}(g)= \frac{g}{2(1+g)} \quad , \quad q_{\sigma}(\overline{d}) = 8\overline{d}\left(1-2\overline{d}\right)=\frac{4g}{(1+g)^2}\quad , \quad \varepsilon_{\sigma}^{\rm eff}(\vec{k})= q_{\sigma}(\overline{d}) \varepsilon_{\sigma}^0(\vec{k})\quad ,\quad \hbox{and} \quad n_{\vec{k},\sigma}(\overline{d})= \left\{ \begin{array}{lcr} (1+q_{\sigma}(\overline{d}))/2 & \hbox{for} & \varepsilon(\vec{k}) \leq 0\\ (1-q_{\sigma}(\overline{d}))/2 & \hbox{for} & \varepsilon(\vec{k}) > 0 \end{array} \right. . </math>

These approximation-free results in infinite dimensions are identical to those obtained from the Gutzwiller Approximation (GA), which was based on a semi-classical counting of configurations (M.C. Gutzwiller, 1965) (J. Bünemann, 1998).

For the single-band Hubbard model, the GA describes a metal-to-insulator transition in which all electrons are localized above a finite critical strength of the Hubbard interaction (Brinkman-Rice transition) (W.F. Brinkman and T.M. Rice, 1970). For <math>U\geq U_{\rm BR}\ ,</math> the minimization of the variational ground-state energy functional leads to <math>\overline{d}^{\rm opt}(U>U_{\rm BR})=q_{\sigma}^{\rm opt}(U>U_{\rm BR})= E_{\rm var}^{\rm opt}(U>U_{\rm BR})=0\ ,</math> where <math>U_{\rm BR}</math> is of the order of the bandwidth of the noninteracting electrons.

The Brinkman-Rice transition is an artifact of the limit of infinite dimensions, i.e. the variational result <math> g_{\rm opt}>0</math> holds for the Hubbard model in finite dimensions for all finite interaction strengths <math>U>0\ .</math> The Brinkman-Rice transition cannot be removed by any finite-order expansion in <math>1/Z</math> (F. Gebhard, 1990). Nevertheless, the Brinkman-Rice transition provides an illustrative example of the breakdown of the Fermi liquid state at the Mott metal-to-insulator transition.

In the Kotliar-Ruckenstein slave-boson approach, each atomic configuration is represented by a boson, so that the Hubbard interaction becomes simple. In contrast, the kinetic energy becomes much more complicated in terms of the bosons because the motion of an electron from one site to another changes the atomic configuration and thus the boson number on both sites (G. Kotliar and A.E. Ruckenstein, 1986).

After an ingenious transformation of the boson transfer operators between sites, the replacement of the bosonic operators by average values (saddle-point approximation) leads to an effective Hamiltonian with the same dispersion relation as eq. \eqref{eq:epseffective}. Therefore, the Kotliar-Ruckenstein slave-boson mean-field theory is identical to the result of the Gutzwiller theory in the limit of infinite coordination number (J. Bünemann and F. Gebhard, 2007).

The application of the GWF to real materials (Gutzwiller theory) involves a three stage approximation. The starting point is the parameterization of the (multi-band) Hubbard model for which the bare dispersion relation <math>\varepsilon_{\sigma,\sigma'}^0(\vec{k})</math> and the energies of the atomic levels must be specified. The second approximation is the GWF itself, which is a variational ground state. Thirdly, the GWF is evaluated in infinite dimensions, but the corresponding expressions are applied to three-dimensional systems.

Like an electron, a Helium-3 atom is a spin-1/2 fermion. At low temperatures, Helium-3 is a normal liquid and also a Landau Fermi liquid. Its properties are well described by the Gutzwiller theory for the single-band Hubbard model because the Helium atoms repel each other strongly. The assumption that the fermions move on a lattice is an additional approximation in the liquid phase.

A reasonable agreement between the theoretical prediction and experimental data is obtained for the pressure dependence of the magnetic susceptibility <math>\chi_{\rm s}</math> and the compressibility <math>\kappa\ .</math> They are related to the Landau Fermi-liquid parameters <math>F_0^{\rm s,a}</math> (D. Vollhardt, 1984) via <math>\chi_{\rm s}/\chi_{\rm s}^0=(m^*/m) (1+F_0^{\rm a})^{-1}</math> and <math>\kappa/\kappa^0=(m^*/m) (1+F_0^{\rm s})^{-1}\ .</math> Here the upper index `0' refers to the noninteracting system. The two Landau parameters are predicted to be functions of the experimentally accessible mass enhancement <math>m^*(p)/m=1/q_{\sigma}(p)=1/(1-(U(p)/U_{\rm BR})^2)=1+F_1^{\rm s}/3\ ,</math> which increases with pressure <math>p\ .</math> The theoretical predictions <math>F_0^{\rm s}(p)=-1+1/(1-I(p))^2</math> and <math>F_0^{\rm a}=-1+1/(1+I(p))^2\ ,</math> with <math>I(p)=\sqrt{1-1/q_{\sigma}(p)}\ ,</math> reproduce the experimental values to an accuracy of 50% and 10%, respectively.

For nickel, the local spin-density approximation (LSDA) to density-functional theory (DFT) does not provide a good description of the quasi-particle bands as measured in ARPES experiments. Essentially all of the discrepancies are resolved by the Gutzwiller theory.

The Gutzwiller theory employs a multi-band Hubbard model with bare band structure <math>\epsilon_{\sigma,\sigma'}^0(\vec{k})\ ,</math> which is obtained from a paramagnetic LDA calculation. Only the bands close to the Fermi energy, i.e., the <math>3d\ ,</math> <math>4s\ ,</math> and <math>4p</math> bands, are taken into account in the Gutzwiller theory. In spherical approximation, the atomic spectrum of the <math>3d</math>-shell depends on three Racah parameters A, B, and C, whereby the parameters <math>B</math> and <math>C</math> are close to their values for free Ni<math>^{++}</math>-ions. The Racah-<math>A</math> corresponds to the Hubbard <math>U</math> and is a free parameter of the Gutzwiller theory. When only the correlations in the <math>3d</math> bands are considered, <math>A=8\, {\rm eV}{\rm -}10\, {\rm eV}</math> leads to a good agreement between the quasi-particle bands from the Gutzwiller theory and experimental ARPES data. This is shown in <figref>Gutzwiller_Wave_Function_GWF-ARPES.jpg</figref> (J. Bünemann, F. Gebhard et al., 2003).

The quality of the theoretical predictions improves when the spin-orbit coupling is taken into account. The theory explains the observed <111>-direction for the magnetic moment, the magnitude of the magnetic anisotropy energy, and the change in the Fermi-surface topology around the X-point in the presence of a strong external magnetic field along the <001>-direction (Gersdorf effect). A Fermi-surface cross section in the presence of the spin-orbit coupling is shown in<figref>Gutzwiller_Wave_Function_GWF-Fermisurfaces.jpg</figref> (J. Bünemann, F. Gebhard et al., 2008).

The Gutzwiller theory is not only applicable to ferromagnetism but to many other types of ground states with broken symmetry, e.g., to antiferromagnetism or to lattice disorder in which the ground state lacks translational symmetry. In the same way, Gutzwiller-correlated BCS (Bardeen-Cooper-Schrieffer) wave functions are candidates for superconductors with strong electronic correlations such as the high-temperature superconductors (P.W. Anderson, 1987). This idea was worked out in more detail by F.C. Zhang et al., 1988; for a review, see (B. Edegger, V.N. Muthukumar et al., 2007). The concept of gossamer superconductivity put forward by Bernevig, Laughlin, and Santiago, 2003, is also based on Gutzwiller-correlated BCS wave functions.

The GWF can equally be applied to the Bose-Hubbard model which is suitable for ultracold bosonic atoms in optical lattices. The Gutzwiller theory reproduces the mean-field result for the phase boundary between the superfluid phase and the Mott phase at zero temperature (D. Jaksch, C. Bruder et al., 1998). The mean-field approach becomes exact in the limit of infinite dimensions and provides a reasonable approximation to bosons in two-dimensional and three-dimensional confined geometries.

'''A historical note on the invention of the Hubbard model'''

Internal references

• Leon Cooper and Dmitri Feldman (2009) Bardeen-Cooper-Schrieffer theory. Scholarpedia, 4(1):6439.
• Giovanni Gallavotti (2008) Fluctuations. Scholarpedia, 3(6):5893.

Hubbard model

Category:Physics Category:Theoretical condensed matter Category:Multiple Curators