Skip to main content

Full text of "On the way to a Gutzwiller density functional theory"

See other formats


On the Way to a Gutzwiller Density Functional 


f^*) ' Werner Weber 1 , Jorg Biinemann 2 , and Florian Gebhard 2 


1 Institut fur Physik, Universitat Dortmund, D-44221 Dortmund, Germany 
^ ' 2 Fachbereich Physik, Philipps-Universitat Marburg, D-35032 Marburg, Germany 

H- •> ' 


Abstract. Multi-band Gutzwiller-correlated wave functions reconcile the contrasting 
concepts of itinerant band electrons versus electrons localized in partially filled atomic 
q_) . shells. The exact evaluation of these variational ground states in the limit of large 

coordination number allows the identification of quasi-particle band structures, and the 
calculation of a variational spinwave dispersion. The study of a generic two-band model 
elucidates the co-operation of the Coulomb repulsion and the Hund's-rule exchange for 

"ti ' itinerant ferromagnetism. We present results of calculations for ferromagnetic nickel, 

using a realistic 18 spin-orbital basis of 4s, Ap and 3d valence electrons. The quasi- 
particle energy bands agree much better with the photo-emission and Fermi surface 

i-rt ' data than the band structure obtained from spin-density functional theory (SDFT). 


O , 1 Exchange versus Correlations 

More than 50 years ago two basically different scenarios had emerged from 
^ ■ early quantum-mechanical considerations on electrons in metals with partly filled 

<"*"> i d bands. 

en . 

Scenario I: As proposed by Slater jj] and Stoner M, band theory alone was 
argued to account for itinerant ferromagnetism. Due to the Pauli princi- 
ple, electrons with parallel spins cannot come arbitrarily close to each other 
("Pauli" or "exchange hole"), and, thus, a ferromagnetic alignment of the 
electron spins reduces the total Coulomb energy with respect to the param- 
agnetic situation ( "exchange field energy" ) . 
Scenario II: As emphasized by van Vleck p], electronic correlations are im- 
portant in narrow-band materials. Due to the strong electron-electron in- 
teraction, charge fluctuations in the atomic d shells are strongly suppressed 
Q | ("minimum polarity model"). The atomic magnetic moments arise due to 

O ■ the local Coulomb interactions (in particular, Hund's-rule couplings) and 

the electrons' motion through the crystal may eventually align them at low 
enough temperatures. 


In principle, such a dispute can be resolved in natural sciences. The correspond- 
ing theories have to be worked out in detail, and their results and predictions 
have to be compared to experiments. 

This was indeed done for scenario I y,g|. The (spin-) density functional theory 
is a refined band theory which describes some iron group metals with consid- 
erable success. Unfortunately, progress for scenario II was much slower. It calls 



2 Werner Weber, Jorg Biinemann, and Florian Gebhard 

for a theory of correlated electrons, i.e., a genuine many-body problem has to 
be solved. It was only recently that reliable theoretical tools became available 
which allow to elucidate scenario II in more detail P,|7|,p|,p PlO| , pd[ | . 

A first step in this direction was the formulation of appropriate model Hamil- 
tonians which allowed to discuss matters concisely, e.g., the Hubbard model p2 13 ljJla . 

This model covers both aspects of d electrons on a lattice: they can move through 
the crystal, and they strongly interact when they sit on the same lattice site. 
The model is discussed in more detail in Sec. |[ 

Even nowadays, it is impossible to calculate exact ground-state properties 
of such a model in three dimensions. In 1963/1964 Gutzwiller introduced a 
trial state to examine variationally the possibility of ferromagnetism in such 
a model |12|,[t3|. His wave function covers both limits of weak and strong cor- 
relations and should, therefore, be suitable to provide qualitative insights into 
the magnetic phase diagram of the Hubbard model. Gutzwiller-correlated wave 
functions for multi-band Hubbard models are defined and analyzed in Sec. |3|. 

The evaluation of multi-band Gutzwiller wave functions itself poses a most 
difficult many-body problem. Perturbative treatments pq , |l7[ | are constrained to 
small to moderate interaction strengths. The region of strong correlations could 

only be addressed within the so-called "Gutzwiller approximation" |12,13J,|l8[ 
and its various extensions p9J , pO| . Some ten years ago, the Gutzwiller approxi- 
mation was found to become exact for the one-band Gutzwiller wave function 
in the limit of infinite spatial dimensions, d — * oo pUE2LE3[, and Gebhard p4| 
developed a compact formalism which allows the straightforward calculation of 
the variational ground-state energy in infinite dimensions. Recently, Gebhard's 
approach was generalized by us to the case of multi-band Gutzwiller wave func- 
tions pi]. Thereby, earlier results by Biinemann and Weber p5|, based on a 
generic extension of the Gutzwiller approximation |26j ], were found to become 
exact in infinite dimensions p7J]. 

As shown in Sect. for a two-band toy model, the Gutzwiller variational 
scheme approach also allows the calculation of spinwave spectra |p8| . In this 
way, the dispersion relation of the fundamental low-energy excitations can be 
derived consistently. Albeit the description is based on itinerant electrons, the 
results for strong ferromagnets resemble those of a Heisenberg model for localized 
spins whereby a unified description of localized and itinerant aspects of electrons 
in transition metals is achieved. 

In Sect. H we discuss results from a full-scale calculation for nickel. The ad- 
ditional local correlations introduced in the Gutzwiller scheme lead to a much 
better description of the quasi-particle properties of nickel than in previous cal- 
culations based on spin-density functional theory. 

2 Hamilton Operator 

Our multi-band Hubbard model [Q is defined by the Hamiltonian 

On the Way to a Gutzwiller Density Functional Theory 3 

Here, c£ a creates an electron with combined spin-orbit index a = 1,...,2JV 
(N = 5 for 3d electrons) at the lattice site i of a solid. 

The most general case is treated in Rcf. |1Q]. In this work we assume for 
simplicity that different types of orbitals belong to different representations of 
the point group of the respective atomic state (e.g., s, p, d{e g ), d(i,2 g ))- In this 
case, different types of orbitals do not mix locally, and, thus, the local crystal field 
is of the from tlf = e a 5 a ^ a ' . Consequently, we may later work with normalized 
single-particle product states |<?o) which respect the symmetry of the lattice, 

(Na+<A CT ,|#o) = <W<* • ( 2 ) 

We further assume that the local interaction is site-independent 

&;at= E ^'^^i^W** • ( 3 ) 

^1 ) (7 '2 :C3 ,CT4 

This term represents all possible local Coulomb interactions. 

As our basis for the atomic problem we choose the configuration states 

10 = 1 01, 02, . . .) = c+ CTl c+ CT2 • • • |vacuum) (en < <r 2 < • • •) , (4) 

which are the "Slater determinants" in atomic physics. The diagonalization of 
the Hamiltonian H^t is a standard exercise [gj|. The eigenstates \T) obey 

\r) =J2T I ,r\I) , (5) 

where T/,r are the elements of the unitary matrix which diagonalizes the atomic 
Hamiltonian matrix with entries {I\H^ ai \I') . Then, 


The atomic properties, i.e., eigenenergies Ep, eigenstates \r), and matrix ele- 
ments Ti t p, are essential ingredients of our solid-state theory. 

3 Multi-band Gutzwiller Wave Functions 

3.1 Variational Ground-State Energy 

Gutzwiller-correlated wave functions are written as a many-particle correla- 
tor Pq acting on a normalized single-particle product state |#o), 

\& G ) = P G \$ ). (7) 

The single-particle wave function |<Pq} which obeys (|2|) contains many configura- 
tions which are energetically unfavorable with respect to the atomic interactions. 
Hence, the correlator Pq is chosen to suppress the weight of these configurations 

4 Werner Weber, Jorg Bunemann, and Florian Gebhard 

to minimize the total ground-state energy of O). In the limit of strong cor- 
relations the Gutzwiller correlator Pq should project onto atomic eigenstates. 
Therefore, the proper multi-band Gutzwiller wave function with atomic correla- 
tions reads 


p i& = II >^r r = II I 1 + ( A ^ - !) **;?] = 1 + Y1 ( A ^ - X ) **r ■ (8) 

r r r 

The 2 2N variational parameters Xi.r per site are real, positive numbers. For 
^i-,r ¥" an d au other Xi-r — all atomic configurations at site i but |i~b) are 
removed from |^o)- Therefore, by construction, \&g) covers both limits of weak 
and strong coupling. In this way it incorporates both itinerant and local aspects 
of correlated electrons in narrow-band systems. 

The class of Gutzwiller-correlated wave functions as specified in (0) was eval- 
uated exactly in the limit of infinite dimensions in Rcf. |lQj. The expectation 
value of the Hamiltonian (|l|) reads |3(| 

{H) m* G ) (9) 

Here, n° a — ($o\ni ;<T \<I>o) is the local particle density in \$o). The local q factors 
are given by |Hj 

m r m r > ^ fi fi' I o o 

x r rVu,)W)^r + ',/^ , (10) 

where mf.j (jn®. r ) is the probability to find the configuration \I) (the atomic 
eigenstate \T)) on site i in the single-particle product state |^o)- The fermionic 
sign function fl = (I\Ja\c£\I) gives a minus (plus) sign if it takes an odd (even) 
number of anticommutations to shift the operator c+ to its proper place in the 
sequence of electron creation operators in \I U a). 

Eqs. (0) and (10) show that we may replace the original variational parame- 

ters Xi : r by their physical counterparts, the atomic occupancies rrii-r- They are 
related by the simple equation |10[ 

m %r = A-. r m? r . (11) 

The probability for an empty site (\I\ = 0) is obtained from the completeness 

m i; = 1 - ^ m i\r • (12) 


On the Way to a Gutzwiller Density Functional Theory 
The probabilities for a singly occupied site (|/| = 1) are given by BQ] 

rUi-a = n?. 


Y^ m *;i > ( 13a ) 


^ r V l;r 

m° . (13b) 

The parameters ra i; and TO i;0 . must not be varied independently. All quantities 
in (9) are now expressed in terms of the atomic multi-particle occupancies m^r 
(\r > 2), the local densities n^. a , and further variational parameters in \<Po). 

It is seen that the variational ground-state energy can be cast into the form 
of the expectation value of an effective single-particle Hamiltonian with renor- 
malized electron transfer amplitudes t\ 


tif = y/qi^y/qj^tif . (14) 


Therefore, |^o) is the ground state of H c g whose parameters have to be de- 
termined self-consistently from the minimization of (<Po\H c g\<Po) with respect to 
m^r and n° CT . For the optimum set of parameters, H°^ defines a band structure 
for correlated electrons. Similar to density-functional theory, this interpretation 
of our ground-state results opens the way to detailed comparisons with experi- 
mental results; see Sect, Ol. 

3.2 Spinwaves 

The variational principle can also be used to calculate excited states [[HJ. If \<P) 
is the ferromagnetic, exact ground state with energy E$, the trial states 

W{q))=S-\<P) (15) 

are necessarily orthogonal to \<P), and provide an exact upper bound to the first 
excited state with momentum q and energy e(q) 

e{q) - Es{q) = Wfimd) ~ ■ ( } 

Here, S~ — (<§i~) + = X)z 6 CX P( — i^O^n, i^z 6 t flips a spin from up to down 
in the system whereby it changes the total momentum of the system by q. In 
this way, the famous Bijl-Feynman formula for the phonon-roton dispersion in 
superfluid Helium was derived |3^| . In the case of ferromagnetism the excitation 
energies E s (q) can be identified with the spinwave dispersion if a well-defined 
spinwave exists at all J28[. Experimentally this criterion is fulfilled for small 
momenta q and energies E s (q). 


Werner Weber, Jorg Biinemann, and Florian Gebhard 

Unfortunately, we do not know the exact ground state or its energy in gen- 
eral. However, we may hope that the Gutzwiller wave function |^g) is a good 
approximation to the true ground state. Then, the states 

\* G (q)) = §7\*c 


will provide a reliable estimate for E s (q), 

E s (q) » ET(q) = 

(y G \S+HS-\* G ) (* G \H\V C 
(<P G \S+S^\V G ) " (^gI^g) 


Naturally, E^ ar (q) does not obey any strict upper-bound principles. 

The actual calculation of the variational spinwave dispersion is rather in- 
volved. However, explicit formulae are available p8| which can directly be applied 
once the variational parameters have been determined from the minimization of 
the variational ground-state energy. 




1 1 

i | i 

FM '_ 









- PM 



■ — ■■ 


A i 





\ FM 




\. - 







4 6 



Fig. 1. Phase diagram as a function of U and J for the Hartree-Fock-Stoner theory 
(HF) and the Gutzwiller wave function (GW) for (a) n = 1.17 and (b) n = 1.40; PM: 
paramagnet, FM: ferromagnet 

On the Way to a Gutzwiller Density Functional Theory 7 

4 Results for a Generic Two-Band Model 

4.1 Ground-State Properties 

The atomic Hamiltonian for a two-band model (b — 1, 2) can be cast into the 

+ J ^_^ c l, CT C 2,-cr C l,-cr C 2,CT + ■'C I C l , | C l , | C 2, X C 2, T "^ C 2,T C 2, J. C 1, J. C 1,T ) «•*■") 


For two d(e g ) orbitals, -ff a t exhausts all possible two-body interaction terms. 
Since we assume that the model describes two degenerate d(e g ) orbitals, the 
following restrictions are enforced by the cubic symmetry |E9j: (i) J = Jq, and 
(ii) U — U' = 2J. Therefore, there are two independent Coulomb parameters, 
the local Coulomb repulsion U (of the order of 10 eV) and the local exchange 
coupling J (of the order of leV, as typical for atomic Hund's rule couplings). 
For the one-particle part H\ we use an orthogonal tight-binding Hamiltonian 
with first and second nearest neighbor hopping matrix elements, resulting in a 
bandwidth W = 6.6 eV. 

In the following we concentrate on two band- fillings, (a), n = 1.17, where the 
non-interacting density of states (DOS) shows a pronounced peak at the Fermi 
energy, most favorably for ferromagnetism, and, (b), n = 1.40, a position near 
the DOS peak, where the DOS exhibits a positive curvature as a function of the 

In Fig. P we display the J-U phase diagram for both fillings. It shows that 
Hartree-Fock theory always predicts a ferromagnetic instability. In contrast, the 
correlated-electron approach strongly supports the ideas of van Vleck M and 
Gutzwiller [ fl3| : (i) a substantial on-site exchange J is required for the occurrence 
of ferromagnetism if, (ii), realistic Coulomb repulsions U are assumed. At the 
same time the comparison of Figs, [lk and Hb shows the importance of band- 
structure effects which are the basis of the Stoner theory. The ferromagnetic 
phase in the U-J phase diagram is much bigger when the density of states at 
the Fermi energy is large. Therefore, the Stoner mechanism for ferromagnetism 
is well taken into account in our correlated-electron approach. 

In Fig. y, we display the energy differences between the paramagnetic and 
ferromagnetic ground states ("condensation energy", -E C ond) as a function of 
the interaction strength for J = 0.2U. This quantity should be of the order of 
the Curie temperature which is in the range of 100 K-1000 K in real materials. 
The Hartree-Fock-Stoner theory yields such small condensation energies only 
in the range of U ~ 4eV; for larger U, E con a is of order U. In any case, the 
interaction parameter U has to be tuned very precisely to give condensation 
energies which concur with experimental Curie temperatures |lj. In contrast, for 
the Gutzwiller-correlated wave function, we find relatively small condensation 
energies E con a — 0.5 • 10 3 K even for interaction values as large as twice the 

Werner Weber, Jorg Biinemann, and Florian Gebhard 

1 i ' i 

i. | i | i | 

^ 2000 

| n^O.29 



" HF 

f n o =0.35 





CD 1000 



o 50 ° 




i . V \: 

, -^r^f " 

4 6 




Fig. 2. Condensation energy as a function of U for J = 0.2(7 for the Hartree-Fock 
theory (HF) and the Gutzwiller wave function (GW) for n — 1.17 (full lines) and n — 
1.40 (dashed lines) 

bandwidth (U « 12 eV). Moreover, the dependence of the condensation energy 
on U is rather weak such that uncertainties in U do not drastically influence the 
estimates for the Curie temperature. 

4.2 Spinwave Dispersions 

In Fig. we show E^ aT ((q, 0,0)), the variational spinwave dispersion (filf), in 
x direction for the model parameters n = 1.17, J = 0.2U, and the four different 
values U/eV — 7.8, 10, 12, 13.6 which correspond to a magnetization per band of 
m = 0.12,0.20,0.26,0.28. This quantity is defined as < m = (n bA - n bd )/2 < 
n/A. Note that our last case corresponds to an almost complete ferromagnetic 
polarization. The data fit very well the formula 

Er((q,0,0)) = Dq 2 (l+Pq 2 ) + O( q 6 ) 


in qualitative agreement with experiments on nickel |33j . The corresponding val- 
ues D = 1.4 eV A 2 and D = 1.2 eV A 2 for m = 0.26 and m = 0.28, respectively, 
are of the right order of magnitude for nickel where D — 0.43 eV A 2 . As lattice 
constant of our simple-cubic lattice we chose a — 2.5 A. 

As shown in the inset of Fig. pi the dispersion relation is almost isotropic for 
q values up to half the Brillouin zone boundary p8[ , in particular for large mag- 
netizations. This is in contrast to the strong dependence of the electron-transfer 
amplitudes tij on the lattice direction. This implies for strong ferromagnets that 
the collective motion of the local moments is similar to that of localized spins 
in an insulator |34] . Such ferromagnetic insulators are conveniently described by 
the Heisenberg model with exchange interactions between neighboring sites (i, j) 
on a cubic lattice, 

H s = -jJ2SiS j . (21) 

On the Way to a Gutzwiller Density Functional Theory 









> 0.4 


- w 0.2 


j( m=Q2.y* 


0.2 0.4 

LU 0.3 

q / A"' y 

/> -^^*m=0.28 


%/ jr 


n n 







Fig. 3. Variational spinwave dispersion in x direction, E™((q, 0,0)), for the two-band 
model defined in Sect. §; n = 1.17, J = 0.2(7, and the values U/eV = 7.8, 10, 12, 13.6 
correspond to m — 0.12,0.20,0.26,0.28. The lattice constant is a — 2.5 A. Inset: 
£ s var ((g,0,0)) and £ s var ((g/ \/2, g/\/2, 0)) for m = 0.2 and m = 0.28, respectively. The 
spinwave dispersion is almost isotropic 

For such a model one finds D — 2SJa 2 . The length of the effective local spins 
can be calculated from \& G ) as S(S + 1) « 0.95 (S = 0.6) for m > 0.20 @. 
Therefore, J w Z?/(1.2a 2 ), which gives the typical value J = 0.17 eV. For an 
estimate of the Curie temperature Tq we use the result from quantum Monte- 
Carlo calculations plj 

T c = 1.44 JS 2 (22) 

for spins S on a simple-cubic lattice. In this way we find Tq ~ 0.5 J = 0.09 eV = 
1-10 3 K. This is the same order of magnitude as the condensation energy for these 
values of the interaction, E^ond = 5 • 10 2 K , see Sect. [|. Given the arbitrariness 
in the relation between E con a and Tc, and the application of the Heisenberg 
model to our itinerant-electron system, we may certainly allow for difference of 
a factor two in these quantities. Nevertheless, the results of this section clearly 
show that, (i), E con( i gives the right order of magnitude for Tc, and that, (ii), 
the spinwave dispersion of strong itinerant ferromagnets resembles the physics 
of localized spins. 

5 Correlated Band-Structure of Nickel 

5.1 Discrepancies Between Experiment and SDFT 

Of all the iron group magnetic metals, nickel is the most celebrated case of 
discrepancies between the results from experiment and from spin-density func- 
tional theory (SDFT) [[36|. From very early on, the photo-emission data have 
indicated that the width of the occupied part of the d bands is approximately 
W* cc = 3.3 eV |3 whereas all SDFT results yield values of W * c ^p FT = 4.5 eV 
or larger P,p7[| . Similarly, the low temperature specific heat data p8[ give a much 

10 Werner Weber, Jorg Bunemann, and Florian Gebhard 

larger value of N*(Ep), the quasi-particle density of states at the Fermi energy 
(3.0 vs. 1.9 states/(eV atom)), which indicates a quasi-particle mass enhance- 
ment by a factor of approximately 1.6. Here, the Sommerfeld formula is used 
to convert the specific heat data; the theoretical value follows directly from the 
quasi-particle band structure. Furthermore, very detailed photo-emission stud- 
ies at symmetry points and along symmetry lines of the Brillouin zone show 
discrepancies to SDFT results for individual band-state energies which are of 
similar magnitude as seen in the overall d bandwidth. 

The studies revealed even bigger discrepancies in the exchange splittings of 
majority spin and minority spin bands. The SDFT results give a rather isotropic 
exchange splitting of about 600 meV P, p7y39[| . In contrast, the photo-emission 
data show small and highly anisotropic exchange splittings between 160 meV for 
pure d(e g ) states such as X 2 and 330 meV for pure dfog) states, the latter value 
estimated from the exchange splitting of A3 states along r to L jl(j,[Il|. The 
much larger and much too isotropic exchange splitting of the SDFT results has 
further consequences. 

1. The experimental magnetic moment of the strong ferromagnet Ni is fi — 
0.61/ie; yet of relevance is its spin-only part ^ sp j n _ on i y = 0.55/ib p2| - The 
SDFT result is /^pin-only = 0.59^b @1 , an overestimate related to the too 
large exchange splitting. 

2. the X2 state of the minority spin bands lies below Ep |43|| , whereas all SDFT 
results predict it to lie above the Fermi level j|,[4j|j4f|T As a consequence, 
the SDFT Fermi surface exhibits two hole ellipsoids around the X point of 
the Brillouin zone whereas in the dc-Haas-van-Alphen experiments only one 
ellipsoid has been found [ flipql . 

3. The strong t2 g -e g anisotropy is also reflected in the total d hole spin density, 
i.e., in the observation that the d- hole part of the Ni magnetic moment has 
81% d(t2g) and 19% dje g ) character J4j, whereas the SDFT results give a 
ratio of 74% to 26% J47J. 

In the late 70's and early 80's various authors have investigated in how far many- 
body effects improve the agreement between theory and experiment, see, e.g., 
Refs. |^,^|. For example, Cooke et al. Q introduced an anisotropic exchange 
splitting as a fit parameter. 

5.2 Present Status of the Gutzwiller-DFT 

Limitations: By construction, the Gutzwiller approach naturally combines 
with density- functional theory (DFT) which provides a basis of one-particle wave 
functions and a 'bare' band structure. The Gutzwiller-DFT introduces important 
local correlations and provides a variational ground-state energy, a quasi-particle 
band structure, and a spin-wave dispersion. 

Nevertheless, the Gutzwiller-DFT has its own limitations which we collect 
here for further reference. 

1. It starts from a model Hamiltonian whose parameters need to be determined 
from a DFT calculation; we shall comment on this procedure below. 

On the Way to a Gutzwiller Density Functional Theory 11 

The true ground state is approximated by a variational many-body wave 
function; however, our experience from the two-band model supports our 
hope that the variational freedom of our wave function is big enough to 
capture the essential features of itinerant ferromagnetism in real materials 
as well. 

The variational ground-state energy is evaluated exactly only in the limit of 
infinite dimensions; however, from the one-band case, we expect 1/d correc- 
tions to be small J24| . 

Similar in spirit to density-functional theory, we interpret the ground-state 
energy in terms of a quasi-particle band structure; it should be kept in mind, 
though, that this quantity is, in general, not identical to the quasi-particle 
dispersion in the sense of standard many-body theory |50(| . 
Most dynamic quantities, e.g., the spectral function, cannot be determined 

within our approach; the example of the spinwave dispersion in Sect. 4.2 
shows, however, that we can calculate low-order moments of spectral func- 
tions consistently. 

Despite all these restrictions, the method remedies many problems of the SDFT 

in the description of the quasi-particle band structure of nickel, see Sect. 5.3 

Parameterization of the One-Particle Hamiltonian: In the present study, 
we determine the hopping matrix elements t°'° in (Jlj) from a least squares' fit 
to the energy bands obtained from a density-functional-theory calculation for 
non-magnetic nickel. An orthogonal nine orbital basis is used, and the root- 
mean-square deviation of the d band energies is about 60 meV. 

A more complete description should include the flexibility of the wave func- 
tions to relax in the magnetic state. This could be achieved by enhancing the 
orbital basis by Ad states. Moreover, spin-orbit coupling is of significance in 
nickel, as it leads to a 10% enhancement of the total magnetic moment. In 
principle, the spin-orbit coupling, or, more generally, an arbitrarily large orbital 
basis can be treated within our formalism pQ|, yet it leads to complications 
such as local q factors which now depend on two spin-orbital indices instead of 
one as in (|l0|). These extensions not only enhance the numerical complexity of 
the problem but also require different methods for extracting the single-particle 
Hamiltonian from DFT, for example by a more direct evaluation of DFT results 
obtained from local basis methods. 

Since we start from a DFT basis, the 'bare' band structure incorporates al- 
ready some important exchange and correlation effects. In particular, we may 
expect that the non-local Coulomb terms are well taken into account because the 
electron-electron interaction is screened at a length scale of the order of the in- 
verse Fermi wave number. In this way, we can restrict all explicit Coulomb inter- 
action terms in H to local interactions. This assumption is supported by the fact 
that the Hartree-Fock approximation becomes exact in infinite dimensions for 

density-density interactions, V r<T ' (T (r ^ 0) = J2i ni l(T ni+r,a' — > V^p (r ^ 0) |plj| . 
Therefore, we expect that interaction terms beyond the purely local Hubbard 

12 Werner Weber, Jorg Biinemann, and Florian Gebhard 

interaction should be properly taken into account in the density-functional ap- 
proach in three dimensions. However, the proper treatment of the "double count- 
ing" problem for both local and non-local interactions remains a serious prob- 
lem for all methods which try to combine density-functional approaches with 
model-based many-particle theories; see, e.g., the contributions by Lichtenstein, 
Vollhardt, and Potthoff in this volume. 

Chemical Potentials: In the translationally invariant system under investiga- 
tion, the local occupation densities are the same as their system averages, 

(fk,*) = {Na)/L , (23) 

where N„ = ^. c^. G c i . a counts the number of electrons with spin-orbit index a. 
Therefore, we may equally work with chemical potentials fi a for each spin-orbit 
index in the Hamiltonian 

H gc =H-J2 M.A . (24) 


In this grand-canonical view, the chemical potentials rather than the particle 
densities act as variational parameters. Naturally, not all of these parameters 
may be varied independently. For example, as a consequence of the hybridization 
of the Asp and the 3d electrons, the 3d levels would be depleted for a strong d-d 
repulsion which needs to be compensated using one of the parameters. Presently 
we keep fixed the values of the 4s and 4p partial charges, and thus also the 3d 
total charge, to the values of the non-magnetic calculation. This is achieved by 
using two of the four chemical potentials for 4s and 4p electrons. 

As can be seen from (H), the chemical potentials act as a shift of the 'bare' 
(DFT) values of the fields e CT , 



H* • (25) 

A b = ef t - ef, . (26) 

In this way, the variational approach naturally contains the flexibility to ad- 
just the magnetic (or "exchange") splitting between majority bands (6, t) and 
minority bands (b, J,) 

_ fc M — t b,l ■ 

In particular, we may allow for an anisotropy in the exchange splittings of the 
d(e g ) and dfag) electrons. 

Interaction Parameters of the Atomic Hamiltonian: Presently we employ 
only the on-site Coulomb interaction within the 3d shell, i.e., all interactions 
within the 4s, 4p shell and between 4sp and 3d are neglected. In spherical atom 
approximation, which is found to be well justified, all matrix elements in (||) can 
either be expressed as a function of the Slater integrals F(k) (k = 0,2,4) or of 
the Racah parameters A, B, C [g9|. We use C/B ~ 4-5 |2J| and determine A and 
C in order to give an optimal agreement with experimental data (effective mass 

On the Way to a Gutzwiller Density Functional Theory 13 

and bandwidth, condensation energy, ti g j e g ratio of the d part of the magnetic 
moment, Fermi surface topology). 

Currently, there is a big debate on the magnitude of the interaction param- 
eters. In principle, the interaction parameters could also be deduced from DFT 
results. However, there is no consensus on how to calculate these parameters 
consistently. For example, they could be calculated from atomic or Wannier 
functions, or they could be found using constrained DFT methods (see, e.g., 
Ref. ||). 

Minimization: The number of multi-electron states \T) is 2 2N = 2 10 . Be- 
cause of the cubic site symmetry, the number of independent variational param- 
eters nip reduces to approximately 200 for the paramagnetic and to approxi- 
mately 400 for the ferromagnetic cases. These "internal" variational parameters 
obey 2 N + 1 sum rules (O) and ( |13a| ) ; in cubic symmetry there remain three 
for the paramagnetic and five for the ferromagnetic cases. There is freedom to 
choose those mr 1 which, through the sum rules, are dependent on the other mr- 
It is advisable to pick those mr 1 which can be expected to have large values. 
This avoids unphysical negative values of mr 1 to occur during the variational 

The chemical potentials of ( P5|) are the "external" variational parameters. 
In the present case these are eight, however three are fixed to yield the total 
4s, 4p, and 3d densities, such that the space of the external parameters is five- 
dimensional. Given a fixed set of external variational parameters, the procedure 
to determine the internal ones begins to put them equal to their uncorrelated 
values mr — mP r . Thus, q° a — 1. Note that q s . a — q p . a — 1 always holds, 
as there is no interaction for 4s, 4p orbitals. From this, the 'bare' (DFT) band 
structure and |^o) barc follow as an initial guess for the quasi-particle band struc- 
ture and one-particle product state. Then, the following self-consistent scheme 
is employed: 

1. Calculate the ground-state energy for |<£n)a ld where a labels the set of exter- 
nal variational parameters. This requires momentum-space integrations up 
to the respective Fermi surface. 

2. Minimize the ground-state energy (J9) with respect to the internal variational 

3. Calculate the q factors and derive |^o)'q OW as the ground state of the H c s ( |14J ) 
with the renormalized hopping matrix elements t^j] repeat steps 1-3 until 
convergence to |^o)a P * is reached. 

Self-consistency is usually reached rather quickly, i.e., |^o)a pt is found after three 
to five iterations. 

The global minimum, |^o)g P b a i * s f° un( i by a search through the space of the 
external variational parameters keeping the average d and sp occupations. This 
search can be sped up by first optimizing with respect to the most important 
external variational parameter which is the isotropic exchange splitting Ad = 
{A e + A t2 )/2, putting the difference to zero as a first approximation. 

14 Werner Weber, Jorg Bunemann, and Florian Gebhard 

In a second step, the anisotropy of the exchange splitting is investigated, 
i.e., we introduce A e and A t2 in the minimization procedure, keeping Ad at 
the value of A° d p obtained in the first optimization step. The searches for A d p , 
and for Z\° pt and A° p can be carried out starting with |<?o) barc - Only then the 
self-consistency procedure for \<Po) opt has to be launched. 

Typical energy gains are (in meV): 

^barc _ £bare^opt) _ 10 _1 0, (27a) 

E^(AT) - E^(A°?\A^) « 5-10, (27b) 

£ barc (Z\°f, A°f a ) - E° pt (Alf,A^) « 5-10 . (27c) 

The energy gains from the variations of A s and A p are of the order of 0.1 meV. 

5.3 Comparison to Experiments 

The results for nickel of our DFT-based Gutzwiller calculations agree best with 
experiment when we choose the following values of the interaction parameters: 
A « 10-12eV, C w 0.1-0.4eV with C/B « 4.5 ||. The width of the d bands 
is predominantly determined by A (essentially the Hubbard U) via the values of 
the hopping reduction factors qd,a- The exchange splittings and, consequently, 
the magnetic moment are mainly governed by C and to some extend also by A. 
The Racah parameter C causes the Hund's-rule splitting of the d 8 multiplets; in 
the hole picture, d 8 is the only many-particle configuration which is significantly 
occupied (by 1.90 electrons), while 5.94 electrons are in d 9 , 0.89 electrons are in 
d 10 , and 1.18 electrons have s or p character. 

In our present study, the parameter C is found to be rather small (0.1 eV) 
compared to A in order to reproduce the measured spin-only moment /^spin-only = 
0.55. Larger values of C move the minimum of the total energy curve l?tot 
vs. magnetization to to values of to « 0.60-0.65/iB- 

There are two points to discuss here. The first concerns the large value of A, 
which seems incompatible with the position of the satellite peak in the photo- 
emission data at about 6 eV below the Fermi energy E-p p& . Model calculations 
for this many-body excitation peak use values of U ~ 3-5 eV. However, these 
models use single of few d band models, excluding hybridization with the 4s, Ap 
bands, see, e.g., Ref. |49(| . When, in our calculation, the hybrization effects are 
switched off, and only the d band contribution to the total energy matters, we 
also find that values of A ta 3-5 eV agree best with experiment, and A w 10 eV 
would be way out of a reasonable range of parameter values. 

The second point concerns the shape of the total energy curve Etotirn) at 
large values of m in the limit of strong ferromagnetism. In this limit, the increase 
of the magnetic moment is fed from the d admixture in the majority 4s, Ap bands. 
Compared to analogous curves obtained from SDFT, the curvature at large to 
values is much smaller in our results. We presume that the larger SDFT curvature 
is related to the balance between 4s, Ap and 3d electrons. It is well known that 
this balance in a delicate manner determines the stability of transition metals 

On the Way to a Gutzwiller Density Functional Theory 15 

as well as of noble metals; see, e.g., Ref. [p4|, and the discussion of this problem 
in Ref. |55} ]. The balance between Asp and 3d electrons is the more influenced 
the larger the exchange splitting fields are because the minority band 3d level is 
shifted towards the 4s, 4p levels and the majority band 3d level is shifted away. 
Only in first order of the splitting energy, we can expect that no change in the 
overall 4s, 4p population happens, as is imposed by the choice of our 4s, Ap 
chemical potentials. Presently, the flow between 4s, Ap and 3d electrons cannot 
be described with our model Hamiltonian as the electron-electron interaction 
within the 4s, Ap shell and between Asp and 3d is not included. 

The exchange splittings not only determine the magnetic moment but also 
influence strongly the shape of the single-particle bands in the vicinity of Ep 
(not the overall bandwidth). For the detailed comparison with photo-emission 
data we have thus either chosen calculations with small C values (0.1 eV), where 
the minimum of _E tot (m) yields m — 0.55/ie, or, for larger C values, with a fixed 
moment constraint, using the experimental spin-only moment of /^spin-only = 
0.55. The resulting quasi-particle bands do not differ much from each other. 
There is however a tendency that values C w 0.4 eV and larger appear to agree 
somewhat better with the bulk of the photo-emission data. 

Generally, the Gutzwiller results agree much better with experiment than the 
SDFT results. For example, this is the case for, (i), the value for the quasi-particle 
density of states at the Fermi energy (^g-dft(-^f) = 2.6 vs. 3.0 states/(eV 
atom)), (ii), the positions of individual quasi-particle energies, (hi), the values 
of the exchange splittings, (iv), their t2 g -e g anisotropy, and, (v), the t2 g /e g ratio 
of the d part of the magnetic moment ((£2g/e s )G-DFT = 83/17 vs. 81/19). As a 
consequence of the small d(e g ) exchange splitting, the X21 state lies below Ep 
and, thus, the Fermi surface exhibits only one hole ellipsoid around X, in nice 
agreement with experiment. 

The large anisotropy of the exchange splittings is a result of our ground-state 
energy optimization, which allows A t2 and A e to be independent variational 
parameters. We find A t2 rs 3Z\ e rs 800 meV. Note that these values enter 
\<P ) haTC and are renormalized by factors g^j « 0.7, qd,i ~ 0.6, when |<?o)° pt is 
reached. This also implies that the width of the majority spin bands is about 10% 
bigger than that of the (higher lying) minority spin bands. It causes a further 
reduction of the exchange splittings of states near E-p , especially for those with 
strong t2 g character. Note that this band dispersion effect causes larger exchange 
splittings near the bottom of the d bands, e.g., 0.45 eV splitting of X\ and 0.74 eV 
splitting of X3. There, however, the quasi-particle linewidths have increased to 
1.25eV and 1.4eV, respectively p7|, so that an exchange splitting near the 
bottom of the d bands could, so far, not be observed experimentally. 

The large anisotropy may originate from peculiarities special to Ni with its 
almost completely filled d bands and its fee lattice structure. Near the top of 
the d bands, the ti g states dominate because they exhibit the biggest hopping 

integrals to nearest neighbors, t dda sw 0.5 eV. The e g states have t dd7T « — 0.3eV 

to nearest neighbors, and t dda fa 0.1 eV to next-nearest neighbors; the latter 

are small because of the large lattice distance to second neighbors. The e g states 

16 Werner Weber, Jorg Bunemann, and Florian Gebhard 

also mix with the nearest-neighbor ti g states with tjJ -type coupling. Therefore, 
the system can gain more band energy by avoiding occupation of anti-bonding 
t2 g states in the minority spin bands via large values of A t2 , at the expense 
of allowing occupation of less anti-bonding e g states via small A e values. This 
scenario should not apply to materials with a bcc lattice structure which have 
almost equal nearest and next-nearest neighbor separations. Since the bands in 
nickel are almost completely filled, the suppression of charge fluctuations actually 
reduces the number of atomic configurations where the Hund's-rule coupling is 
active. It is also in this respect that nickel does not quite reflect the generic 
situation of other transition metals with less completely filled d bands. 

The results for nickel presented here must be seen as preliminary inasmuch 
some important interaction terms were not yet included; see Sect. |5.2| . However, 
the present study already shows that the Gutzwiller-DFT is a working approach. 
It should allow us to resolve many of the open issues in itinerant ferromagnetism 
in nickel and other transition metals. 

6 Conclusions and Outlook 

Which scenario for itinerant ferromagnetism in transition metals is the correct 

Band theory along the lines of Slater and Stoner could be worked out in 
much detail whereas a correlated-elcctron description of narrow-band systems 
was lacking until recently. Our results for a two-band model and for nickel show 
that the van-Vleck scenario is valid. Band theory alone does not account for the 
strong electronic correlations present in the material which lead to the observed 
renormalization of the effective mass, exchange splittings, bandwidths, and Fermi 
surface topology. Moreover, charge fluctuations are indeed small, and large local 
moments are present both in the paramagnetic and the ferromagnetic phases. 

Roughly we may say that the electrons' motion through the crystal leads to 
a ferromagnetic coupling of pre-formed moments which eventually order at low 
enough temperatures. In this way, strong itinerant ferromagnets resemble ferro- 
magnetic insulators as far as their low-energy properties are concerned: spinwaves 
exist which destroy the magnetic long-range order at the Curie temperature. 

Our present scheme allows us a detailed comparison with data from refined 
photo-emission experiments on nickel which are currently carried out [|5q| . It 
should be clear that our approach is applicable not only to nickel but to all 
other itinerant electron systems. 

Despite all recent progress much work remains to be done. The present imple- 
mentation of the Gutzwiller-DFT needs to be improved by the inclusion of more 
orbits, their mutual Coulomb interaction terms, and the spin-orbit coupling. Ul- 
timately, some of the principle limitations of our variational approach will have 
to be overcome by a fully dynamic theory. Most probably, such a theory will re- 
quire enormous numerical resources such that a fully developed Gutzwiller-DFT 
will always remain a valuable tool to study ground-state properties of correlated 
electron systems. 

On the Way to a Gutzwiller Density Functional Theory 17 


We gratefully acknowledge helpful discussions with all participants of the Her- 
aeus seminar Ground-State and Finite-Temperature Bandferromagnetism. This 
project is supported in part by the Deutsche Forschungsgemeinschaft under 
WE 1412/8-1. 


1. J.C. Slater, Phys. Rev. 49, 537 (1936); ibid., 931 (1936). 

2. E.C. Stoner, Proc. Roy. Soc. A 165, 372 (1938); for early reviews, see J.C. Slater, 
Rev. Mod. Phys. 25, 199 (1953) and E.P. Wohlfarth, ibid., 211 (1953). 

3. J.H. van Vleck, Rev. Mod. Phys. 25, 220 (1953). 

4. V.L. Moruzzi, J.F. Janak, and A.R. Williams, Calculated Electronic Properties of 
Metals (Pergamon Press, New York, 1978). 

5. See also the contributions in this volume by O. Erikson; R. Wu; J. Kiibler and 
K.H. Bennemann; R. Brinzanik; P.J. Jensen. 

6. W. Nolting, W. Borgiel, V. Dose, and Th. Fauster, Phys. Rev. B 40, 5015 (1989); 
W. Borgiel and W. Nolting, Z. Phys. B 78, 241 (1990). 

7. H. Hasegawa, J. Phys. Soc. Jpn 66, 3522 (1997); Phys. Rev. B 56, 1196 (1997); 
R. Fresard and G. Kotliar, Phys. Rev. B 56, 12909 (1997). 

8. Th. Obermeier, Th. Pruschke, and J. Keller, Phys. Rev. B 56, 8479 (1997); 
Th. Maier, M.B. Zolfl, Th. Pruschke, and J. Keller, Euro. Phys. J. B 7, 377 
(1999); M.B. Zolfl, Th. Pruschke, J. Keller, A.I. Poteryaev, LA. Nekrasov, and 
V.I. Anisimov, Phys. Rev. B 61, 12810 (2000). 

9. D. Vollhardt, N. Bliimer, K. Held, M. Kollar, J. Schlipf, M. Ulmke, and J. Wahle, 
Adv. in Solid-State Phys. 38, 383 (1999); LA. Ne krasov, K. Held, N Bliimer, A.I. 

Poteryaev, V.I. Anisimov, D. Vollhardt, preprint |cond-mat/0005207j (2000). 

10. J. Biinemann, W. Weber, and F. Gebhard, Phys. Rev. B 57, 6896 (1998). 

11. See also the contributions in this volume by A.I. Lichtenstein; D.M. Edwards 
and A. CM. Green; D. Vollhardt; W. Nolting, M. Potthoff, T. Herrmann, and 
T. Wegner; A.M. Oles and L.L. Feiner. 

12. M.C. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963). 

13. M.C. Gutzwiller, Phys. Rev. 134, A923 (1964); ibid. 137, A1726 (1965). 

14. J. Hubbard, Proc. Roy. Soc. London Ser. A 276, 238 (1963); ibid. 277, 237 (1964). 

15. J. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 

16. G. Stollhoff and P. Fulde, J. Chem. Phys. 73, 4548 (1980); G. Stollhoff and P. Thal- 
meier, Z. Phys. B 43, 13 (1981); A.M. Oles and G. Stollhoff, Phys. Rev. B 29, 314 
(1984); for further details on the "local ansatz" technique, see P. Fulde, Electron 
Correlations in Molecules and Solids, Springer Series in Solid-State Sciences 100 
(Springer, Berlin, 1991). 

17. D. Baeriswyl and K. Maki, Phys. Rev. B 31, 6633 (1985); D. Baeriswyl, 
J. Carmelo, and K. Maki, Synth. Met. 21, 271 (1987). 

18. D. Vollhardt, Rev. Mod. Phys. 56, 99 (1984). 

19. K.A. Chao and M.C. Gutzwiller J. Appl. Phys. 42 1420 (1971); K.A. Chao, Phys. 
Rev. B 4 4034 (1971); ibid. 1088 (1973); J. Phys. C 7 127 (1974). 

20. P. Fazekas, Lecture Notes on Electron Correlation and Magnetism, Series in Mod. 
Cond. Matt. Phys. 5 (World Scientific, Singapore, 1999), gives an introduction to 
the theory of ferromagnetism, and a concise description and some applications of 
the Gutzwiller approximation. 

18 Werner Weber, Jorg Bunemann, and Florian Gebhard 

21. W. Metzner and D. Vollhardt, Phys. Rev. Lett. 59, 121 (1987); Phys. Rev. B 37, 
7382 (1988). 

22. W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 (1989). 

23. For a review, see F. Gebhard, The Mott Metal-Insulator Transition (Springer, 
Berlin, 1997). 

24. F. Gebhard, Phys. Rev. B 41, 9452 (1990). 

25. J. Bunemann and W. Weber, Phys. Rev. B 55, 4011 (1997). 

26. J. Bunemann, Eur. Phys. J. B 4, 29 (1998). 

27. J. Bunemann, F. Gebh ard, and W. Weber J. Phys. Cond. Matt. 8, 7343 (1997). 

28. J. Bunemann, preprint |cond-mat/0005154| (2000). 

29. S. Sugano, Y. Tanabe, and H. Kamimura, Multiplets of Transition-Metal Ions in 
Crystals, Pure and Applied Physics 33 (Academic Press, New York, 1970). 

30. This holds for our symmetry-restricted basis. 

31. A. Messiah, Quantum Mechanics, 3rd printing (North Holland, Amsterdam, 1965). 

32. R.P. Feynman, Statistical Mechanics, Frontiers in Physics 36 (Benjamin, Reading, 

33. R.D. Lowde and C.G. Windsor, Adv. Phys. 19, 813 (1970). 

34. See, e.g., S.V. Halilov, H. Eschrig, AY. Perlov, and P.M. Oppeneer, Phys. 
Rev. B 58, 293 (1998). 

35. K. Chen, A.M. Ferrenberg, and D.P. Landau, Phys. Rev. B 48, 3249 (1993). 

36. For a review, see S. Hufher, Photoelectron Spectroscopy (Springer, Berlin, 1995). 

37. W. Eberhardt and E.W. Plummer, Phys. Rev. B 21, 3245 (1980). 

38. M. Dixon, F.E. Hoare, T.M. Holden, and D.E. Moody, Proc. R. Soc. A 285, 561 

39. J. Callaway in Physics of Transition Metals, ed. by P. Rhodes (Conf. Ser. Notes 55, 
Inst, of Physics, Bristol, 1981), p. 1. 

40. M. Donath, Surface Science Reports 20, 251 (1994). 

41. K.-P. Kamper, W. Schmitt, and G. Giintherodt, Phys. Rev. B 42, 10696 (1990). 

42. H.A. Mook, Phys. Rev. 148, 495 (1966). 

43. R. Raue, H. Hopster, and R. Clanberg, Phys. Rev. Lett. 50, 1623 (1983). 

44. C.S. Wang and J. Callaway, Phys. Rev. B 15, 298 (1977). 

45. E.P. Wohlfahrt in Handbook of Magnetic Materials 1, ed. by E.P. Wohlfarth 
(North Holland, Amsterdam, 1980). 

46. D.C. Tsui, Phys. Rev. 164, 561 (1967). 

47. O. Jepsen, J. Madsen, and O.K. Andersen, Phys. Rev. B 26, 2790 (1982). 

48. J.F. Cooke, J.W. Lynn, and H.L. Davis, Phys. Rev. B 21, 4118 (1980). 

49. A. Liebsch, Phys. Rev. B 23, 5203 (1981). 

50. A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems 
(McGraw-Hill, New York, 1971). 

51. E. Miiller-Hartmann, Z. Phys. B 74, 507 (1989); ibid. 76, 211 (1989). 

52. G. Vielsack and W. Weber, Phys. Rev. B 54, 6614 (1996). 

53. Some preliminary results can be found in J. Bunemann, F. Gebhard, and W. We- 
ber, Found. Phys. 30 (Dec. 2000). 

54. D.G. Pettifor, J. Magn. Magn. Mat. 15-18,847 (1980). 

55. J. Hafner, From Hamiltonians to Phase Diagrams: The Electronic and Statistical 
Mechanical Theory of Sp-Bonded Metals and Alloys (Springer Series in Solid-State 
Sciences 70, 1987), pp. 72. 

56. R. Claessen, private communication (2000).