PHYSICAL REVIEW LETTERS
PRL 105, 086401 (2010)
week ending
20 AUGUST 2010
Quantum Renormalization of the Spin Hall Effect
Bo Gu,1,2 Jing-Yu Gan,3 Nejat Bulut,4 Timothy Ziman,5 Guang-Yu Guo,6,7 Naoto Nagaosa,8,9 and Sadamichi Maekawa1,2
1
Advanced Science Research Center, Japan Atomic Energy Agency, Tokai 319-1195, Japan
2
JST, CREST, 3-Sanbancho, Chiyoda-ku, Tokyo 102-0075, Japan
3
Institute of Physics, Chinese Academy of Sciences, P.O. Box 603, Beijing 100190, China
4
Department of Physics, Izmir Institute of Technology, Urla 35430, Turkey
5
CNRS and Institut Laue Langevin, Boı̂te Postale 156, F-38042 Grenoble Cedex 9, France
6
Graduate Institute of Applied Physics, National Chengchi University, Taipei 116, Taiwan
7
Department of Physics, National Taiwan University, Taipei 106, Taiwan
8
Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan
9
Cross-Correlated Material Research Group and Correlated Electron Research Group, RIKEN-ASI, Wako 315-0198, Japan
(Received 30 November 2009; published 18 August 2010)
By quantum Monte Carlo simulation of a realistic multiorbital Anderson impurity model, we study the
spin-orbit interaction (SOI) of an Fe impurity in Au host metal. We show, for the first time, that the SOI is
strongly renormalized by the quantum spin fluctuation. Based on this mechanism, we can explain why the
gigantic spin Hall effect in Au with Fe impurities was observed in recent experiments, while it is not
visible in the anomalous Hall effect. In addition, we show that the SOI is strongly renormalized by the
Coulomb correlation U. Based on this picture, we can explain past discrepancies in the calculated orbital
angular momenta for an Fe impurity in an Au host.
DOI: 10.1103/PhysRevLett.105.086401
PACS numbers: 71.70.Ej, 75.30.Kz, 75.40.Mg
When a magnetic impurity with d orbitals is put into a
normal metal with s- or p-conduction bands, it forms a
virtual bound state hybridized with the conduction states,
leading to a broadening with width (hybridization energy) which is typically of the order of eV. This energy
scale competes with the Coulomb energy U between the
electrons in the d orbitals, which is also of the order of eV.
The latter tends to produce the spin moment, while the
former induces quantum fluctuation of that spin moment,
leading to the spin singlet. The competition between these
two interactions defines the time scale, or equivalently the
energy scale, of the quantum spin fluctuation, i.e., the
Kondo temperature TK , which can be much lower than
both and U [1]. In the language of the renormalization
group, this is described by the energy-dependent scaling of
the various operators. This means that even a weak interaction can be amplified in the low energy or temperature
scale and compete with the much larger energy scale due to
the correlation U.
In real systems, the d orbitals have fivefold degeneracy,
and these orbital degrees of freedom have often been
neglected in the analysis of the experimental results. The
reasoning is that either the crystal field splitting " is
much smaller than the hybridization energy or that it is
in the limit with " much larger than the Kondo temperature TK . However, the naive comparison between the bare
interaction strengths is dangerous since these are scaledependent running coupling constants. Similar nontrivial
behavior can also be expected for the relativistic spin-orbit
interaction.
Orbital degrees of freedom in the impurity scattering
lead to intriguing phenomena such as the anomalous Hall
0031-9007=10=105(8)=086401(4)
effect (AHE) and the spin Hall effect (SHE). A charge
current perpendicular to the applied electric field is produced in ferromagnetic metals (AHE), while a spin current
rather than a charge current is induced in semiconductors
and metals without magnetism (SHE). An extrinsic mechanism of these two effects arises from the skew scattering,
i.e., spin-dependent deflection of the scattered electrons
due to the spin-orbit interaction. It has been shown in the
case of the AHE that the resonant skew scattering due to
the virtual bound d states leads to a large Hall angle, of the
order of 0.01 [2], compared with the typical value of 103 .
In the AHE, the spin fluctuation is quenched, owing to the
ferromagnetism or the external magnetic field, and the ratio
of the spin-orbit interaction and basically determines
the Hall angle. On the other hand, in the SHE, the spin
fluctuation is active and the Kondo physics can be relevant
to the resonant skew scattering. Therefore, a crucial question is whether the Kondo effect and quantum spin fluctuation can produce an even larger spin Hall angle
compared with the AHE.
An important clue to this question comes from an experiment by Seki et al. [3] on Au=FePt, in which a spin
Hall angle of 0.114 was obtained. Motivated by this experiment, some of the present authors studied the Fe impurity in Au by a first-principles calculation [4]. The
Kondo effect of Fe impurities is a historic problem, with
a low TK around 0.4 K [5], high electric resistance at room
temperature [1], and an AHE with Hall angle of the order
of 0.01 [2]. However, a simple fivefold degeneracy of the
orbitals has been assumed to analyze the experiments. In
Ref. [4], on the other hand, the orbital-dependent Kondo
effect of Fe in Au was proposed to explain the nature of the
086401-1
Ó 2010 The American Physical Society
(ck ) is the creation (annihilation) operator of
the conduction electron with wave vector k and spin in
the band , dy (d ) is the creation (annihilation) operator of the localized electron at the impurity site with orbital
and spin , and n ¼ dy d . The host energy band
ðkÞ, the impurity energy levels , and the impurity-host
hybridization Vk , as the one-body problems, can be
properly obtained by the DFT calculations. U (U0 ) is the
on-site Coulomb repulsion within (between) the orbitals of
the impurity, and J is the Hund coupling between the
orbitals of the impurity. These many-body interactions
can be exactly treated by the QMC calculations, and thus
it becomes possible to accurately study the quantum fluctuations of the spins and orbitals of the impurity.
Considering the parameters used in the previous calculations for Fe in Au [4], and the relationship U ¼ U0 þ 2J
[12], in our following QMC calculations, we use the values
of U ¼ 5 eV, J ¼ 0:9 eV, and U0 ¼ 3:2 eV for most
(∑α|Vξkα|2)1/2 (eV)
(1)
(a)
2
0
L
1
z
where cyk
4
<(Mξ)2>
J X
n n 0 ;
2 Þ0 ;
0.5
0
-2
1
2
ξ=z
xz
yz
2 2
x -y
xy
Tχξ
k;;;
X
X
U0 X
þ n þ U n" n# þ
n n 0 0
2 Þ0 ;;0
;
cases, but we shall vary the values for a few cases in order
to clarify the role of correlations and reconcile with past
calculations.
Our DFT calculations are done by the code QUANTUMESPRESSO [13]. To calculate the impurity-host hybridization, we consider the supercell Au26 Fe, where the
exchange-correlation interactions are described by the
Perdew-Zunger local density approximation (LDA), and
the electron-ion interactions are represented by the RabeRappe-Kaxiras-Joannopoulos ultrasoft pseudopotentials.
(See Ref. [14] for details.)
Figure 1(a) shows the hybridization between orbitals
of an Fe impurity and an Au host, where the hybridization
matrix element has the form of Vk h’ jH0 j ðkÞi.
H0 is the one-particle part of Eq. (1), ’ is the state of the
Fe impurity, and ðkÞ is the Au host state with wave
vector k and band index . (See Ref. [14] for details.) It is
observed that, at the point (k ¼ 0), the hybridization
value of the ¼ eg ðz2 ; x2 y2 Þ orbitals of the Fe impurity
is smaller than that of the ¼ t2g (xz, yz, xy) orbitals.
Based on the above DFT/LDA calculation, we can determine approximate impurity levels ¼ 1:9 eV for
¼ eg (z2 , x2 y2 ) and ¼ 1:8 eV for ¼ t2g (xz,
yz, xy) with zero Fermi energy, where the Coulomb terms
included in the LDA calculations have been subtracted as
in Refs. [15,16]. The crystal field splitting " ¼ 0:1 eV is
in agreement with the previous LDA calculations [4,6].
To study the renormalization due to the correlation effect
for (i) the crystal field splitting and (ii) the relativistic spinorbit interaction, the magnetic behaviors of Eq. (1) are
calculated by the Hirsch-Fye QMC simulation for the
infinite host system [8]. Owing to calculation constraints,
we simplify to a three-orbital model with one eg and two
Γ
0.5
1
(b)
ξ = 1(z2)
2(p1)
3(p-1)
-1
log10T (eV)
0
(c)
0
-2
X
<nξ>
experimentally observed giant spin Hall signals [4]: It was
argued that the eg orbitals of Fe are in the Kondo limit and
t2g orbitals are in the mixed-valence region. The enhancement of the spin-orbit interaction by electron correlation in
the t2g orbitals leads to the giant spin Hall effect. However,
this proposal has been challenged theoretically by Ref. [6],
which suggests an effective 3-channel Kondo model, involving local and band electrons of t2g symmetry, and also
experimentally by the x-ray magnetic circular dichroism
(MCD) [7], which obtained a rather small value of the
orbital angular momentum in contrast to the large value
calculated in Ref. [4].
In order to resolve this confusing situation and quantify
the mechanism for the enhanced spin Hall effect, it is
essential to treat the quantum fluctuations of the spins
and orbitals systematically. This is impossible in the firstprinciples calculation, which assumes ordered spin and
orbital moments. In this Letter, we overcome this difficulty
by using the Hirsch-Fye quantum Monte Carlo (QMC)
simulation [8], combined with the density functional theory (DFT) [9,10], to study the renormalization due to
correlation effects. First, a single-impurity, multiorbital
Anderson model [11] is formulated within the DFT for
determining the host band structure, the impurity levels,
and the impurity-host hybridization. Second, the magnetic
behaviors of the Anderson impurity at finite temperatures
are calculated by QMC.
The single-impurity, multiorbital Anderson model is
defined as
X
X
H¼
ðVk dy ck þ H:c:Þ
ðkÞcyk ck þ
k;;
week ending
20 AUGUST 2010
PHYSICAL REVIEW LETTERS
PRL 105, 086401 (2010)
-1
log10T (eV)
0
-1
0
(d)
0.5
0
-2
log10T (eV)
FIG. 1 (color online). (a) Hybridization between the orbitals
of an Fe impurity and an Au host. We show the LDA result
including five d orbitals. (b) Temperature dependence of the
square of the magnetic moment hðMz Þ2 i, (c) the temperature
times susceptibility T , and (d) the occupation number hn i of
the orbitals of the Fe impurity. (b)–(d) are the QMC results
including three d orbitals. See text for details.
086401-2
t2g orbitals, which captures the essential physics. The eg
orbital is arbitrarily chosen as z2 , and the two t2g orbitals
are chosen as p1 and p1 , where the notation corresponds
to the transformational properties of t2g orbitals equivalent
to effective p orbitals [17]: p1 p1ffi2ffi ðxz iyzÞ, p0
ixy, and p1 p1ffi2ffi ðxz þ iyzÞ. The following results
are obtained with more than 105 Monte Carlo sweeps and
the Matsubara time step ¼ 0:25.
Figures 1(b)–1(d) show the temperature dependence of
the square of the magnetic moment hðMz Þ2 i, the temperature times susceptibility T , and the occupation number
hn i, which are defined as Mz ¼ n" n# , ¼
R
z
z
0 dhM ðÞM ð0Þi, n ¼ n" þ n# , respectively, and
the impurity levels ¼ 1:9 ( ¼ z2 ) and 1:8 eV ( ¼
p1 , p1 ) are used. It is found that the eg (t2g ) orbital has a
larger (smaller) magnetic moment, a much larger (smaller)
susceptibility with a much smaller (larger) Kondo temperature, and a larger (smaller) occupation number. In addition,
if we repeat the QMC calculation with the degenerate
impurity levels by hand, ¼ 1:85 eV for ¼ z2 , p1 ,
and p1 , nearly the same behaviors are observed. Thus it is
clear that the orbital-dependent Kondo effect comes
mainly from the renormalization of the impurity-host hybridization by correlations.
Next, we study problem (ii)—the renormalization of the
relativistic spin-orbit interaction due to the correlation. For
simplicity, we consider only the z component,
Hso ¼ ð=2Þð‘z z Þ;
week ending
20 AUGUST 2010
PHYSICAL REVIEW LETTERS
‘z z n2" n2# n3" þ n3# ;
(2)
where ‘z is the 3 3 z component of the angular moment
matrix of ‘ ¼ 1, and z is the 2 2 z component of the
Pauli matrix. ¼ 2 (3)denotes the p1 (p1 ) orbital. We
add Eq. (2) to Eq. (1), where the parameters in Eq. (1) are
taken as the same values of Figs. 1(b)–1(d). For an Fe atom,
the realistic value of the spin-orbit interaction is ¼
75 meV [18]. We also show the results with the smaller
value of ¼ 40 meV for comparison. Considering that the
value of is 2 orders of magnitude less than that of
impurity energies , it is not surprising that the temperature dependence of the square of the magnetic moment
hðMz Þ2 i, the susceptibility T , and the occupation number
hn i are nearly the same as those of Figs. 1(b)–1(d). In
contrast, as displayed in Fig. 2, a nonzero spin-orbit correlation function h‘z z i appears when the spin-orbit interaction within the t2g orbitals of the Fe impurity is included.
At temperature T ¼ 360 K, the lower limit of our present
calculations, we have h‘z z i ffi 0:44 with ¼ 75 meV,
and h‘z z i ffi 0:3 with ¼ 40 meV.
To compare the QMC results with the experiment, we
calculate the spin Hall angle s as follows. Since we
consider only two t2g orbitals with the z component of
the orbital angular moment ‘z ¼ 1, the spin-orbit interaction within the t2g orbitals gives rise to the difference in
the occupation numbers between the parallel (nP ) and
antiparallel (nAP ) states of the spin and angular momenta.
These occupation numbers are related to the phase shifts
P and
AP , respectively, as nPðAPÞ ¼ PðAPÞ = . These
quantities can be estimated as
h‘z z i ¼
P
AP ;
hn2 i þ hn3 i ¼
Pþ
AP ;
(3)
which are given in Figs. 2 and 1(d) , respectively. Putting
h‘z z i ¼ 0:44 for ¼ 75 meV, and hn2 i ¼ hn3 i ¼
0:65, we obtain P ¼ 1:35 and AP ¼ 2:73. Taking into account the estimate 1 ffi 0:1 of the phase shift for p-wave
scattering, and applying the equation [4]
s ¼
6 Im½ðe2i 1 1Þðe2i P e2i AP Þ=½25 15 cos2 P
10 cos2 AP , the spin Hall angle is thus obtained as s ffi
0:055, comparable to that observed in recent experiments.
To further understand how correlations renormalize the
effective relativistic spin-orbit interaction, we will now
vary the correlation energy and the spin polarization, while
keeping the temperature T ¼ 360 K and spin-orbit interaction ¼ 75 meV fixed.
Effect of correlation U.—For each U, we keep J=U ¼
0:9=5 and U ¼ U0 þ 2J. In Fig. 3 with hðM2z þ M3z Þ=2i ¼
0, we see that the spin-orbit correlation function jh‘z z ij
dramatically decreases as U decreases. This demonstrates
clearly how the correlation U renormalizes the relativistic
spin-orbit interaction. The results also explain past discrepancies in the calculated orbital angular momenta for an Fe
impurity in an Au host: A large value is calculated with
U ¼ 5 eV [4], while much smaller ones are obtained with
U ¼ 3 or 2 eV [6,19]. The virtual bound states of t2g have
majority and minority parts, whose contributions to h‘z z i
are opposite in sign. Increasing U will push up the minority
part away from the Fermi energy, decrease its contribution
to h‘z z i, and hence increase the total h‘z z i.
Effect of spin polarization.—To see clearly the effects of
quantum fluctuations, it is instructive to study the effects of
polarizing the Fe impurity spin, which can quench quantum fluctuations of the spins and orbitals. We apply a
0
-0.1
<lzσz>
PRL 105, 086401 (2010)
-0.2
-0.3
λ = 75 meV
40 meV
-0.4
-0.5
-2
-1
0
log10T (eV)
FIG. 2 (color online). Temperature dependence of spin-orbit
interactions for an Fe impurity in an Au host. The parameter in
Eq. (2) is taken as a realistic 75 meV [18], and as 40 meV for
comparison.
086401-3
PHYSICAL REVIEW LETTERS
PRL 105, 086401 (2010)
z z
<l σ >
-0.1
-0.2
-0.3
-0.4
z z
<l σ >
0
-0.1
-0.2
-0.3
-0.4
-0.5
AHE [Ref.2]
XMCD [Ref.7]
2
DFT [Ref.6]
3
U (eV)
4
5
DFT [Ref.4]
1
0.8
0.6 z /2 >
0.4 z M 3)
+
0.2
M2
<(
0
SHE [Ref.3]
FIG. 3 (color). For an Fe impurity in an Au host, we show the
spin-orbit interaction h‘z z i as a function of correlation energy
U and spin polarization hðM2z þ M3z Þ=2i. The QMC results of
h‘z z i are shown as squares with solid lines. In the phase
diagram plane of U to hðM2z þ M3z Þ=2i, the different colors
denote the different values of h‘z z i. Here, T ¼ 360 K and ¼
75 meV are fixed.
magnetic field h by hand on the orbitals 1ðz2 Þ, 2ðp1 Þ, and
3ðp1 Þ of the Fe impurity,
Hf ¼ hðn1" n1# þ n2" n2# þ n3" n3# Þ;
(4)
where h may be of a fraction of eV in order to polarize the
spins at the temperature of 360 K. We add Eq. (4) to the
previous Hamiltonian with the parameters in Eq. (1) taken
to be the same as in Fig. 2. It is found that, in Fig. 3 with
U ¼ 5 eV, when h is set up to 0.8 eV to make the spins
hMz i of ¼ 1ðz2 Þ, 2ðp1 Þ and 3ðp1 Þ gradually polarize
near the saturated value 1, the spin-orbit correlation function jh‘z z ij gradually decreases to 0. Thus, the quantum
fluctuation of the spin is crucially important in producing
an even larger spin Hall angle than the Hall angle of the
AHE. What happens is that the quantum fluctuation of the
spin leads to the motional narrowing of the spin-orbit
coupled virtual bound states, reducing the hybridization
energy with the conduction electrons and hence the energy width . So, the Hall angle will increase with decreasing [2].
Phase diagram.—Combining the above two parameters,
our QMC results for the spin-orbit correlation function
h‘z z i are noted as squares with solid lines in Fig. 3. In
the phase diagram plane of the correlation U to the spin
polarization hðM2z þ M3z Þ=2i, as also noted in Fig. 3, the
different colors note the different values of h‘z z i. It is
seen that the areas with large spin-orbit correlation should
possess strong correlation energy U and weak spin polarization. The measurement of the orbital angular momentum by x-ray MCD involves the application of a sum rule
for atomic wave functions: In order to eliminate interference between different impurities, it is necessary to polarize with an external magnetic field at low temperatures [7].
The spin Hall angle measurement [3] was done in a spinunpolarized state and at room temperature. From our phase
week ending
20 AUGUST 2010
diagram in Fig. 3, it is suggested that the orbital angular
momentum in a spin-polarized state could be much smaller
than that in a spin-unpolarized state. We emphasize that
this phase diagram is obtained for temperatures comparable to room temperature, as is relevant for spin Hall
applications; for quantitative comparison of very low temperature measurements, such as x-ray MCD and phase
decoherence, the scales can be renormalized.
To conclude, we have studied the SHE for an Fe impurity in an Au host by QMC simulation of a multiorbital
Anderson impurity model. Our phase diagram with respect
to correlation and spin polarization shows that large spinorbit correlations occur for strong correlation U and weak
spin polarization, which, in principle, may resolve the
current, somewhat confusing, situation of experiments
and calculations. More generally, we can, for the first
time, quantify an essential difference between the AHE
and the SHE. While we have focused on the specific case of
Fe in Au, the results suggest a more general mechanism,
and constraints on parameters, for obtaining a large SHE
by impurities.
This work was supported by Grant-in-Aids for Scientific
Research in the priority area ‘‘Spin current’’ with
No. 19048015, No. 19048009, No. 19048008, and
No. 21244053, the Next Generation Super Computing
Project, the Nanoscience Program from MEXT of Japan,
the Funding Program for World-Leading Innovative R&D
on Science and Technology (FIRST Program), and NSC of
Taiwan. We thank the Supercomputer Center at the ISSP,
University of Tokyo. T. Z. is grateful to the ICC-IMR.
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
[11]
[12]
[13]
[14]
[15]
[16]
[17]
[18]
[19]
086401-4
G. Gruener et al., Rep. Prog. Phys. 37, 1497 (1974).
A. Fert et al., J. Magn. Magn. Mater. 24, 231 (1981).
T. Seki et al., Nature Mater. 7, 125 (2008).
G. Y. Guo, S. Maekawa, and N. Nagaosa, Phys. Rev. Lett.
102, 036401 (2009).
J. W. Loram et al., Phys. Rev. B 2, 857 (1970).
T. A. Costi et al., Phys. Rev. Lett. 102, 056802 (2009).
W. D. Brewer et al., Phys. Rev. Lett. 93, 077205 (2004).
J. E. Hirsch and R. M. Fye, Phys. Rev. Lett. 56, 2521
(1986).
P. Hohenberg et al., Phys. Rev. 136, B864 (1964).
W. Kohn et al., Phys. Rev. 140, A1133 (1965).
P. W. Anderson, Phys. Rev. 124, 41 (1961).
S. Maekawa et al., Physics of Transition Metal Oxides,
Series in Solid State Sciences Vol. 144 (Springer,
New York, 2004), p. 16.
P. Giannozzi et al., http://www.quantum-espresso.org.
B. Gu et al., J. Phys. Conf. Ser. 200, 062007 (2010).
V. I. Anisimov et al., Phys. Rev. B 48, 16929 (1993).
G. K. H. Madsen et al., Europhys. Lett. 69, 777 (2005).
C. J. Ballhausen, Introduction to Ligand Field Theory
(McGraw-Hill, New York, 1962), p. 96.
A. T. Costa et al., Phys. Rev. B 78, 054439 (2008).
S. Shadov et al., Europhys. Lett. 82, 37001 (2008).