Echap1 3
Echap1 3
Echap1 3
Fig. 1.3. The hcp and dhcp crystal structures. In the latter, the
B and C sites have hexagonal symmetry, while the A sites have local
cubic symmetry, for an ideal c/a ratio.
Fig. 1.4. The Brillouin zones for the hexagonal and fcc lattices.
effective procedure for the rare earths is to adopt one of the linear meth-
ods of Andersen (1975). In the following, we will use the Atomic Sphere
Approximation (ASA) which will allow us to illustrate the construction
and characteristics of the energy bands in a transparent way. This ap-
proximation, and the closely-related Linear Muffin-Tin Orbitals Method
(LMTO), which allows computationally very efficient calculations of ar-
bitrarily precise energy bands, for a given potential, have been concisely
described by Mackintosh and Andersen (1980) and, in much more detail,
by Skriver (1984).
In a close-packed solid, the electrons may to a very good approxima-
tion be assumed to move in a muffin-tin potential, which is spherically
symmetric in a sphere surrounding each atomic site, and constant in the
interstitial regions. We recall that the atomic polyhedron, or Wigner–
Seitz cell, is bounded by the planes which perpendicularly bisect the
vectors joining an atom at the origin with its neighbours, and has the
same volume as the atomic sphere, whose radius S is chosen accordingly.
If we surround each site in the crystal with an atomic sphere, the po-
tential within each of these overlapping regions will, to a high degree of
accuracy, be spherically symmetric. Neglecting the spin, we may there-
fore write the solutions of the Schrödinger equation for a single atomic
sphere situated at the origin in the form
where the radial function Rl (ε, r) satisfies eqn (1.2.12) and is a function
of the continuous energy variable ε. Examples of such radial functions
are shown in Fig. 1.5.
20 1. ELEMENTS OF RARE EARTH MAGNETISM
Dl (ε) + l + 1
pl (ε) = , (1.3.4)
Dl (ε) − l
Rl (ε, S)
Dl (ε) = S . (1.3.5)
Rl (ε, S)
1.3 THE METALLIC STATE 21
with √ ∗
Sl m ,lm (R) = gl m ,lm 4π(−i)λ Yλµ (R̂)(R/S)−λ−1 ,
where
and
λ ≡ l + l ; µ ≡ m − m .
From (1.3.3) and (1.3.7), the required tail-cancellation occurs if
[Pl (ε)δl l δm m − Slk m ,lm ]ajk
lm = 0, (1.3.9)
lm
22 1. ELEMENTS OF RARE EARTH MAGNETISM
Fig. 1.6. The canonical bands for the fcc structure. The band structure
in the metal may be obtained by placing, scaling, and distorting the
canonical bands according to the values of the corresponding potential
parameters Cnl , µnl , and γnl , and finally hybridizing them. The extent of
the bands, according to the Wigner–Seitz rule, is indicated on the right.
1.3 THE METALLIC STATE 23
In this determinantal equation for the band structure εj (k), the infor-
mation about the structure is separated from that on the potential. The
structure constants Slk m ,lm are canonical in the sense that they de-
pend only on the crystal structure and not, for example, on the lattice
constant, as may be seen from the definition (1.3.8), and the potential
function Pl (ε) is determined entirely by the potential within the atomic
sphere. We shall consider these two terms in turn.
If we include values of l up to 3, i.e. s, p, d, and f partial waves, the
structure constants form a square matrix with 16 rows and columns. The
terms with l = l fall into 4 blocks, and these submatrices may be diag-
onalized by a unitary transformation from the lm to an lj representation.
k
The (2l+1) diagonal elements Slj of each sub-block are the unhybridized
canonical l bands. The canonical bands for the fcc structure are shown
in Fig. 1.6. If hybridization is neglected, which corresponds to setting
to zero the elements of Slk m ,lm with l = l , eqn (1.3.11) takes the simple
form
k
Pl (ε) = Slj . (1.3.12)
Dl (Vnl ) = l (1.3.13)
and, within a particular region, the energy Cnl of the centre of the band
is fixed by the condition that Pl (Cnl ) = 0, or
The allowed k-values corresponding to this energy are just those for
24 1. ELEMENTS OF RARE EARTH MAGNETISM
Fig. 1.7. The nth period of the logarithmic derivative function Dl (ε),
and the corresponding potential function Pl (ε). The bottom, centre, and
top of the nl band are defined respectively by Pl (Bnl ) = −2(2l+1)(l+1)/l
(Dl (Bnl ) = 0), Pl (Cnl ) = 0, and Pl (Anl ) = l (Dl (Anl ) = −∞).
k
which Slj = 0 and, since the average over the Brillouin zone may be
shown to vanish, i.e.
2l+1
k
Slj dk = 0, (1.3.15)
j=1 BZ
are then known respectively as the top and bottom of the nl band, even
though this designation is not precisely accurate.
Over the energy range Anl –Bnl , the potential function may be
parametrized with reasonable accuracy as
1 ε − Cnl
Pl (ε) . (1.3.18)
γnl ε − Vnl
The pure l bands are thus obtained from the corresponding canonical
bands by fixing the position with Cnl , scaling the bandwidth by µnl S 2 ,
and distorting them with γnl .
Hybridization between bands of different l is taken into account
by including the structure constants with l = l in (1.3.11), causing
a repulsion between energy levels with the same k and symmetry, as
specified by the labels in Fig. 1.6. Bands of the same symmetry are
thus not allowed to cross, and strong hybridization instead creates an
energy gap. In addition, weak hybridization gives rise to a mixing and
repulsion between bands which do not cross in the absence of hybridiza-
tion. In order to complete the calculation of the band structure, the
inaccuracies due to approximating the atomic polyhedron by a sphere,
and to neglecting higher partial waves, may be conveniently treated to-
gether by perturbation theory. In practice, the energy bands are not
of course calculated step-wise as described above, but all the steps are
performed simultaneously on a computer. Nevertheless, the conceptual
description of the procedure as a placing, scaling and distortion of the
26 1. ELEMENTS OF RARE EARTH MAGNETISM
6s 6p 5d 4f
constant, was therefore consistent with the standard model. However the
positron-annihilation experiments of Gustafson and Mackintosh (1964)
showed that the change in f occupancy, when the transition was induced
by a change in temperature, was much less than one, and indeed that
the results in both phases were consistent with about one f electron per
ion. Similar results were obtained by Gustafson et al. (1969) when the
transition was driven by pressure, and they concluded that it involves
not primarily a change in the f occupancy but rather a change in the f
state, from being localized in the γ-phase to being an itinerant band elec-
tron in the α-phase. This idea was taken up by Johansson (1974) who,
from a consideration of spectroscopic, cohesive and thermodynamic evi-
dence, proposed that the γ–α transition should be considered as a Mott
localized–delocalized transition among the f electrons. Glötzel (1978)
used density-functional theory to calculate the ground state properties
and showed that the equation of state in the α-phase can be accounted
for rather satisfactorily by including the f electrons in the band struc-
ture, and furthermore that a transition to a spin-polarized state should
occur at a lattice constant close to that of γ–Ce, though at a (nega-
tive) pressure considerably lower than that deduced from experiment.
Eriksson et al. (1990) have recently shown that this discrepancy may be
substantially reduced by including the l–l coupling, which is responsible
for the second of Hund’s rules, in the calculation of the 4f bands. This
leads to a ground state in γ–Ce in which the 4f electrons are almost
fully polarized, thus occupying the Hund’s-rule ground state on each
site. Despite the fact that they are described in the band picture, they
may thus be considered as localized, making very little contribution to
the cohesive properties. The calculated atomic volumes in both phases
are in good agreement with experiment. Podloucky and Glötzel (1983)
found a cohesive energy for α–Ce in accord with the measured value,
while that of a ‘promotional’ state with no f electrons is far too small.
They were also able to account for the Compton-scattering experiments
of Kornstädt et al. (1980), who had verified that the change in f occu-
pancy at the transition is small. Skriver (1985) calculated the crystal
structure and equation of state of α–Ce up to high pressures, finding
very good agreement with experiment (Staun Olsen et al. 1985), pro-
vided that the f bands are included, but very poor agreement if the
f electrons are promoted to the d bands, or are assumed to be local-
ized, and therefore to make a negligible contribution to the electronic
pressure. The relative stability at high pressures of low-symmetry con-
figurations such as the α–U structure, which is observed experimentally,
is a strong indicator that there are f electrons in the conduction bands,
as in the light actinides, where they play a decisive role in determining
the structure (Skriver 1985).
1.3 THE METALLIC STATE 29
Fig. 1.9. The band structure of dhcp Pr, calculated by Skriver. The
energy bands in the vicinity of the Fermi level are predominantly d-like,
and the 4f states are assumed to be localized and therefore do not appear,
in contrast to Fig. 1.8.
Fig. 1.10. The intersections of the Fermi surfaces, for the two spin
states of the seventh band in dhcp Pr, with the faces of the Brillouin zone
of Fig. 1.4. The surfaces are generated by a rigid splitting of the energy
bands of Fig. 1.9 by 10 mRy. The unshaded majority-spin surface is a
single sheet, whereas the exchange splitting modifies the topology of the
shaded minority-spin surface, giving rise to a closed lens at M, a small
electron pocket, and an irregular tube along HK.
and clear evidence for this transition has been observed in the dHvA
effect. The changes of the Fermi surface in a magnetic field, and partic-
ularly the enhancement of the effective masses by the interaction with
the 4f moments (Wulff et al. 1988), which we will discuss further in Sec-
tion 7.3, give an average value of I of about 9 mRy, with a variation of
some 30% over different bands and orbits. The agreement between the
measured and calculated electron orbits is such that shifts in the energy
bands of only a few mRy are required to bring the two into concordance,
and this is comparable to typical values for transition metals (Mackin-
tosh and Andersen 1980). The experimental study of the dHvA effect in
Pr, which is the most elaborate which has yet been undertaken for a rare
earth metal, has thus led to the important conclusion that energy-band
32 1. ELEMENTS OF RARE EARTH MAGNETISM
The Fermi surface of paramagnetic Lu, in which the 4f states are all
filled, has also been studied by the dHvA effect (Johanson et al. 1982)
and found to be in semi-quantitative agreement with the calculations of
Tibbetts and Harmon (1982). Since the results of band structure calcu-
lations have been confirmed experimentally at the Fermi level in widely
separated elements in the rare earth series, it is reasonable to suppose
that they will also be successful in accounting for other ground-state
properties. Characteristic band energies for the trivalent lanthanides,
calculated by Skriver (1983) at a common atomic volume close to the
equilibrium value for Gd, are shown in Fig. 1.12. In this figure, the effect
of the change in potential is thus separated from that of the interatomic
spacing. The most notable feature is the fall in energy of the s band
relative to the d band with increasing atomic number, which results in a
decrease of the occupation of the latter, with consequences, as we shall
34 1. ELEMENTS OF RARE EARTH MAGNETISM
see, for the crystal structure. The reason for the fall in the band en-
ergies is the increase of the nuclear charge with atomic number, which
is incompletely screened by the additional f electrons. The potential
veff (r) in (1.2.12) is therefore on average increasingly negative, and in
order to maintain an unchanged boundary condition, as expressed by
the logarithmic derivative, the band energies must decrease accordingly.
1.3 THE METALLIC STATE 35
This effect is relatively modest for the d bands, but much greater for
the s and p bands. The relative shift of the s and d bands is reduced by
the adjustment of the lattice to its equilibrium configuration, but only
by a small amount. As may be seen from Fig. 1.12, (Bd − Bs ) increases
from 101 mRy for La to 373 mRy for Lu at constant S, whereas the cor-
responding values for the equilibrium atomic volumes are 136 mRy and
380 mRy. The band masses also change across the series; µd at constant
volume increases from 2.1 in La to 3.0 in Lu, so that the d bands narrow
as they fall, while µs increases slightly with atomic number, but remains
below 1 throughout (Skriver 1983).
The canonical-band theory may be used to calculate the electronic
pressure and its partitioning between the different angular momentum
components. According to the force theorem (see Mackintosh and Ander-
sen 1980) the change in the total energy, due to an infinitesimal change
in the lattice constant, may be determined as the difference in the band
energies, calculated while maintaining the potential unchanged. We may
thus write
εF
dU = δ εN (ε)dε, (1.3.24)
where N (ε) is the total electronic density of states, and δ indicates the
restricted variation with a frozen potential. The electronic pressure is
then given by
dU
P=− , (1.3.25)
dΩ
where Ω is the volume of the atomic polyhedron. The expression (1.3.21)
for the canonical-band energies then leads to the approximate result for
the l partial pressure:
δCl δ ln µl S 2
3Pl Ω = −nl + nl (εl − Cl ) , (1.3.26)
δ ln S δ ln S
the pressure while the d and f states provide the binding, through their
negative, attractive partial pressure. This difference is essentially due to
the fact that the s and p wavefunctions have a positive curvature at the
atomic sphere, over the energy range of the corresponding bands, as illus-
trated in Fig. 1.5, while the d and f functions have a negative curvature.
Consequently, a decrease in volume causes an increase in the logarithmic
derivative for the former and a decrease for the latter, and since Dl (E)
is a decreasing function of energy, the s and p bands must rise and the d
and f bands fall in order to maintain the boundary condition. Equation
(1.3.25) then immediately accounts for the signs of the corresponding
partial pressures. The attractive f pressure for α–Ce is substantial; if
it is removed, the lattice expands to a volume greater than that of γ–
Ce. The partial pressures at a constant atomic volume for the trivalent
rare earths are shown in Fig. 1.13. As may be seen, it is primarily the
decrease in the s and p pressures, which has its origin in the incompletely
Fig. 1.13. The partial 6s, 6p, and 5d pressures for the trivalent rare
earths, calculated for a common atomic volume close to the equilibrium
value for Gd, after Skriver (1983). It is the decrease in the s and p
pressures which gives rise to the lanthanide contraction.
1.3 THE METALLIC STATE 37
Fig. 1.14. The equilibrium atomic radii for the rare earth metals, after
Skriver (1983). The full circles indicate the experimental values. The full
line is a calculation including the s, p, and d partial pressures, while the
broken line indicates that the f contribution is also taken into account.
Fig. 1.15. The occupation numbers of the 5d states for the trivalent
lanthanides, at the observed equilibrium atomic volumes, after Skriver
(1983). For Ce, the 4f electrons are included in the energy bands. The
experimentally observed crystal structures are labelled by h, s, d, and
f, for hcp, Sm-structure, dhcp, and fcc, respectively. The empirical d-
occupation numbers which separate the different structures are indicated
by the lines on the right.
1.3 THE METALLIC STATE 39