# Semiclassical Boltzmann theory of spin Hall effects in giant Rashba systems

###### Abstract

For the spin Hall effect arising from strong band-structure spin-orbit coupling, a semiclassical Boltzmann theory reasonably addressing the intriguing disorder effect called side-jump is still absent. In this paper we describe such a theory of which the key ingredient is the spin-current-counterpart of the semiclassical side-jump velocity (introduced in the context of the anomalous Hall effect). Applying this theory to spin Hall effects in a two-dimensional electron gas with giant Rashba spin-orbit coupling, we find largely enhanced spin Hall angle in the presence of magnetic impurities when only the lower Rashba band is partially occupied.

###### pacs:

72.10.-d, 72.10.Bg, 72.25.-b## I Introduction

It is now generally accepted that three mechanisms – intrinsic, side-jump and skew scattering – contribute to both the spin Hall effect (SHE) and anomalous Hall effect (AHE) Nagaosa2010 (); Sinova2015 (). Among the three mechanisms, the side-jump mechanism is of special interest because it originates from scattering but can, in some simple cases Nagaosa2010 (); Sinova2015 (); Kovalev2010 (); Culcer2010 (); Yang2011 (), be independent of both the disorder density and scattering strength. In particular, when the SHE or AHE arises from strong spin-orbit coupling in the band-structure, the side-jump belongs to the category of disorder-induced interband-coherence effect which has recently been an important topic in condensed matter physics Nagaosa2010 (); Sinova2015 (); Kovalev2010 (); Sinitsyn (); Hou2015 (); Culcer2017 (); Xiao2017SOT ().

In investigating transport phenomena in solids, the semiclassical Boltzmann approach is appealing due to its conceptual intuition Ziman1960 (). In the study of SHE-AHE, how to incorporate side-jump effects into the semiclassical formalism is an attracting theoretical issue Nagaosa2010 (). In the study of AHE, the renewed semiclassical theory addressing this issue has proven useful in obtaining physical pictures Sinitsyn (). In such a theory, the quantum mechanical information on side-jump is coded in the expressions of gauge-invariant classical concepts such as the coordinate-shift and side-jump velocity Sinitsyn (). On the other hand, in the field of SHE when the spin is not conserved due to strong spin-orbit coupling in the band structure, such as in a Rashba two-dimensional electron gas (2DEG), a semiclassical description to side-jump SHE is still absent Zhang2005 (); noteZhang2005 (). Although the modified Boltzmann equation Sinitsyn () developed in studying the AHE can be directly applied to SHE, the spin–current-counterpart of the side-jump velocity in this case has not been addressed before.

In the present paper we formulate a semiclassical Boltzmann framework of SHE when the spin is not conserved due to strong band-structure spin-orbit coupling. This semiclassical theory takes into account interband-coherence effects induced by both the dc uniform electric field and weak static disorder. We work out the spin-current-counterpart of the side-jump velocity based on scattering-induced modifications of conduction-electron states. When the electric field turns on, this quantity contributes one part of the side-jump SHE.

As applications we consider the SHE in a 2DEG with giant Rashba spin-orbit coupling and short-range impurities. We focus on the enhancement of spin Hall angle when the Fermi energy is tuned down towards and below the band crossing point in giant Rashba 2DEGs with magnetic disorder. The spin Hall angle which measures the generation efficiency of the transverse spin current from the longitudinal electric current is the figure of merit of the SHE. Giant Rashba spin-orbit coupling energy comparable to or even larger than the Fermi energy is possible in the polar semiconductor BiTeX (X=Cl, Br and I) family and related surfaces and interfaces Eremeev2012 (); Sakano2013 (); Wu2014 (). Thus these systems are promising to realize efficient conversion of charge current into spin current.

The paper is organized as follows. In Sec. II we outline the semiclassical formulation of SHE. Section III introduces the Rashba model and calculates the SHE. Section IV concludes the paper.

## Ii Semiclassical formulation

Considering the linear response of the spin current polarized in one particular direction (z direction is chosen in the following) to a weak dc uniform electric field in non-degenerate multiband electron systems in the weak disorder regime, one has the semiclassical formula

(1) |

where is the amount of spin current carried by the conduction-electron state denoted by , is the semiclassical distribution function.

The conduction-electron state may be modified by the electric field and static impurity scattering, can thus deviate from the customary pure-band value Ziman1960 (); Zhang2005 (): with the spin-current operator. Here is the Bloch state, is the band index, is the crystal momentum, and are the plane-wave and periodic parts of , respectively. Following the recipe based on the quantum-mechanical perturbation theory for the electric-field modified Bloch state and the Lippmann-Schwinger equation for the scattering modified conduction-electron state in Ref. Xiao2017SOT, , in the weak disorder regime nontrivial corrections caused by interband-coherence effects to read:

(2) |

The intrinsic correction arises from the interband-virtual-transition (electron charge e) induced by the electric field Zhang2005 (), with the energy of Bloch state and the velocity operator. Thus

(3) |

is an electric-field-induced interband-coherence effect.

The extrinsic correction originates from the interband-coherence during the elastic electron-impurity scattering process. The scattering-induced modification to conduction-electron states is captured by the Lippmann-Schwinger equation describing the scattering state with the T-matrix related to the disorder potential and disorder-free Hamiltonian . denotes the scattering-induced modification to the Bloch state. Thus is related to the values of and in the lowest nonzero order in the disorder potential. Here denotes the average over disorder configurations and we assume that the statistical average of the disorder potential is zero (nonzero value only shifts the origin of total energy) . Only the terms containing interband matrix elements of represent the disorder-induced interband-coherence effects, therefore note-SO scattering ()

(4) |

It has been shown Xiao2017SOT () that the side-jump velocity which is an important ingredient in the semiclassical theory of AHE Sinitsyn () can also be obtained in this way () and thus shares the same origin. can therefore be deemed as the spin–current-counterpart of the side-jump velocity in the case of band-structure spin-orbit coupling.

The properly modified steady-state linearized Boltzmann equation in the presence of weak static disorder has been proposed as Sinitsyn ():

(5) |

where is the band velocity, is the Fermi distribution function, is the coordinate-shift in the scattering process () Sinitsyn () and the semiclassical scattering rate (). Up to the linear order of the electric field one has the decomposition Sinitsyn (); Xiao2017AHE ()

(6) |

with the normal part of the out-of-equilibrium distribution function satisfying the Boltzmann equation in the absence of and the anomalous distribution function related to . It is now clear Sinitsyn () that is a disorder-induced interband-coherence effect and so is .

Given that the semiclassical formulation is relevant in the weak disorder regime, Eq. (1) reduces to Xiao2017SOT ()

(7) |

up to the zeroth order of total impurity density and scattering strength. represents the value of in the lowest Born order Sinitsyn (). In higher Born orders, some additional contributions to appear and are responsible for the transverse transport due to the breakdown of the principle of microscopic detailed balance. The analysis of these higher-Born-order contributions under the non-crossing approximation has been detailed in Ref. Xiao2017AHE, . Here we only mention that there is an interband-coherence scattering effect called “intrinsic-skew-scattering induced side-jump” appearing in the third Born order under the Gaussian disorder. Below we set which is just the intrinsic contribution to the spin current independent of the disorder Zhang2005 (), and because it is related to the spin–current-counterpart of the side-jump velocity. In general case of SHE induced by strong band-structure spin-orbit coupling, is just one part of the side-jump SHE arising from disorder-induced interband-coherence effects. Other two semiclassical contributions to the side-jump SHE (from the anomalous distribution function and the intrinsic-skew-scattering induced side-jump) note-sj () and the skew scattering SHE arising from non-Gaussian disorder are all included in the first term of Eq. (7) Sinitsyn (); Xiao2017AHE ().

To be more clear we can consider the case of randomly distributed scalar pointlike Gaussian disorder with density and average strength . Then , , , and the third-Born-order contribution to behaves as (thus is called the intrinsic-skew-scattering Sinitsyn (); Xiao2017AHE ()). In this case the side-jump SHE may consist of three semiclassical contributions in the zeroth order of both the impurity density and scattering strength: , and the intrinsic-skew-scattering induced side-jump.

## Iii Model calculation

### iii.1 Model

The model Hamiltonian of a Rashba 2DEG is , where is the 2D wavevector, is the effective mass, is the vector of Pauli matrices, the Rashba coefficient. The internal eigenstates read , where label the two bands , and .

For the corresponding wave number in band is given as . Here measures the momentum splitting of two Rashba bands, whereas . The density of states of band takes the form , with .

For , the iso-energy surface slices the spectrum two rings of radii , where denote the two monotonic segments (Fig. 1), . The density of states of branch reads .

The conventional definition of the spin current as an anti-commutator of velocity and spin is employed: . It is purely off-diagonal in band-index space in this model: , thus the SHE in Eq. (7) is determined only by

(8) |

The Boltzmann equation can be conveniently solved by using variables for and for . Correspondingly, for . If , in the lowest Born order the energy-integrated elastic scattering rate is . Whereas if there exists elastic scattering between the two branches and , and one has . Assuming isotropic disorder potential, transport-time type solutions to exist Xiao2016PRB (). For we have

(9) |

where the transport time is determined by

(10) |

For we have

(11) |

with the transport time decided by

(12) |

### iii.2 Calculations

We consider the impurity potential is produced by randomly distributed short-range scatters at , i.e., with and the unity matrix in spin space Lu2013 (). Here the short-range potential is approximated by the delta-potential. We assume Gaussian disorder approximation and isotropic magnetic scattering Lu2013 (); Inoue2006 (). and are the concentrations of nonmagnetic and magnetic impurities, respectively. and are the average strengths for the nonmagnetic and magnetic scattering, respectively. The external electric field is applied in x direction.

#### iii.2.1 Nonmagnetic impurities

When , straightforward calculation leads to the spin-current-counterpart of the side-jump velocity

(13) |

with . The transport time reads Xiao2016PRB () , and then the side-jump spin Hall current is

(14) |

which completely cancels out the intrinsic spin Hall current . This just reproduces the well-known Sinova2015 () vanishing spin Hall current in the semiclassical Boltzmann theory for the first time.

When , the intrinsic spin Hall current is . Meanwhile the spin-current-counterpart of the side-jump velocity reads

(15) |

and thus the side-jump spin Hall current

(16) |

again cancels out the intrinsic one. This also coincides with the zero SHE obtained by the Kubo formula Grimaldi2006 ().

#### iii.2.2 Magnetic impurities

For isotropic delta-like magnetic impurity potential, since the contributions from and cancel out in Eq. (4), the spin–current-counterpart of the side-jump velocity is given by

(17) |

The transport time is given by

(18) |

for , and

(19) |

for , with .

When both Rashba bands are partially occupied, the side-jump spin Hall current

(20) |

enhances the total spin Hall current to . This is the same as the weak-disorder-limit value of that obtained by Kubo diagrammatic calculations Inoue2006 (); Wang2007 (). The longitudinal electric current is , the spin Hall angle is therefore

(21) |

When only the lower Rashba band is partially occupied, the side-jump and the total spin Hall currents are

(22) |

and , respectively. The longitudinal electric current is , thus

(23) |

Although is a small quantity in giant Rashba systems, the factor can be very small leading to a large spin Hall angle when is located close to the band bottom of the lower Rashba band. For instance, if , leads to , which is quite large Ebert2015 (); Fert2011 (). Smaller and smaller may lead to larger . However, the quantitative analysis of this possibility is beyond the scope of the semiclassical theory which is valid only in the weak disorder regime. From the above equation, goes to infinity as goes to the band bottom of the lower Rashba band. But this low carrier density limit is actually beyond the Boltzmann regime, and more rigorous microscopic treatments are called for.

#### iii.2.3 Both nonmagnetic and magnetic impurities

The coexistence of nonmagnetic and magnetic impurities may be the more realistic case Inoue2006 (); Lu2013 (). Only main results will be given in this case. Since there is no mixing between the nonmagnetic and magnetic scattering as pointed out by Inoue et al. Inoue2006 (), the spin-current-counterpart of the side-jump velocity is . The total spin Hall current reads , depending on the relative weight of different types of scattering Yang2011 (); Lu2013 (). The side-jump effect vanishes when , the same condition as that for the vanishing of the ladder vertex correction in the Kubo diagrammatic calculation in the spin- basis for the case Inoue2006 (); note-vertex ().

For Fermi energies above and below the band crossing point, the spin Hall angles are

and

respectively. Here we define . Tuning the ratio one can find that changes monotonically and continuously from the scalar-disorder-dominated case to the magnetic-disorder-dominated regime.

## Iv Discussion and Summary

Before concluding this paper, we comment on some important issues not mentioned in above sections.

First, the simple form of the semiclassical Boltzmann equation (5) is exactly valid only for isotropic bands and isotropic scattering Nagaosa2010 (); Xiao2017AHE (). In our opinion, in the presence of anisotropy a more generic and complicated form of the Boltzmann equation may be necessary, we refer the readers to Ref. Luttinger1957, for detailed discussions.

Second, the recently highlighted “coherent skew scattering” under the Gaussian disorder beyond the non-crossing approximation Ado2015 () is also included in the first term of Eq. (7). This additional contribution is also in the zeroth order of both the impurity density and scattering strength in the weak disorder limit in the presence of only one type of disorder, like the side-jump contribution, but is not an interband-coherence scattering effect note-sj (). Thus how to place this contribution into the classification of AHE-SHE mechanisms suggested in Refs. Nagaosa2010, ; Sinova2015, is still an open question. Therefore, in presenting our theory we avoid this issue. Fortunately, in the Rashba model considered in Sec. III the first term of Eq. (7) vanishes. Besides, we should remind the interested readers that this so-called “coherent skew scattering” has actually already been proposed sixty years ago by Kohn and Luttinger Luttinger1957 (); Luttinger1958 (). We will provide a comprehensive description of a semiclassical Boltzmann theory going beyond the non-crossing approximation in a future publication.

Finally, in the presence of spin-orbit coupling the electron spin is not conserved thus the spin current is not uniquely defined. The conventionally defined spin current adopted in this study is not a conserved transport current. A physically attracting definition of the conserved spin current has been suggested by Shi et al. Shi2006 () by introducing the torque dipole moment. However, disorder effects on the torque dipole spin current Nagaosa2006 () in the Bloch representation are hard to deal with under the uniform external electric field in the Boltzmann theory. We reserve these for future studies.

In summary, we have formulated a semiclassical Boltzmann framework of spin Hall effects induced by strong band-structure spin-orbit coupling in non-degenerate multiband electron systems in the weak disorder regime. We worked out the absent ingredient in previous semiclassical theories, i.e., the spin–current-counterpart of the semiclassical side-jump velocity. This gauge-invariant quantity arises from the interband-coherence during the elastic electron-impurity scattering, and contributes one part of the side-jump spin Hall effect.

Applying this theory to a 2DEG with giant Rashba spin-orbit coupling, we showed an enhanced spin Hall angle when only the lower Rashba band is partially occupied in the presence of magnetic impurities. We note that this energy regime below the band crossing point in Rashba systems and similar systems is of intense theoretical interest also from the standpoint of enhanced efficiency of spin-orbit torque and of Edelstein effect Murakami2012 (); Xiao2016FOP (); Zhang2016 (), as well as enhanced thermoelectric conversion efficiency Wu2014 (); Xiao2016PRB ().

###### Acknowledgements.

The author is indebted to Qian Niu, Dingping Li and Zhongshui Ma for insightful discussions.## References

- (1) J. Sinova, S. O. Valenzuela, J. Wunderlich, C. H. Back, and T. Jungwirth, Spin Hall effects, Rev. Mod. Phys. 87, 1213 (2015).
- (2) N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and N. P. Ong, Anomalous Hall effect, Rev. Mod. Phys. 82, 1539 (2010).
- (3) A. A. Kovalev, J. Sinova, and Y. Tserkovnyak, Anomalous Hall Effect in Disordered Multiband Metals, Phys. Rev. Lett. 105, 036601 (2010).
- (4) D. Culcer, E. M. Hankiewicz, G. Vignale, and R. Winkler, Side jumps in the spin Hall effect: Construction of the Boltzmann collision integral, Phys. Rev. B 81, 125332 (2010).
- (5) S. A. Yang, H. Pan, Y. Yao, and Q. Niu, Scattering universality classes of side jump in the anomalous Hall effect, Phys. Rev. B 83, 125122 (2011).
- (6) N. A. Sinitsyn, Semiclassical theories of the anomalous Hall effect, J. Phys.: Condens. Matter 20, 023201 (2008).
- (7) D. Hou, G. Su, Y. Tian, X. Jin, S. A. Yang, and Q. Niu, Multivariable Scaling for the Anomalous Hall Effect, Phys. Rev. Lett. 114, 217203 (2015).
- (8) D. Culcer, A. Sekine, and A. H. MacDonald, Interband coherence response to electric fields in crystals: Berry-phase contributions and disorder effects, Phys. Rev. B 96, 035106 (2017).
- (9) C. Xiao and Q. Niu, Semiclassical theory of spin-orbit torques in disordered multiband electron systems, preprint, arXiv: 1705.07947, to be published in Phys. Rev. B.
- (10) J. M. Ziman, Electrons and Phonons (Clarendon, Oxford, 1960).
- (11) S. Zhang and Z. Yang, Intrinsic Spin and Orbital Angular Momentum Hall Effect, Phys. Rev. Lett. 94, 066602 (2005).
- (12) For instance, one can verify that the semiclassical framework described in Ref. Zhang2005, cannot produce any extrinsic contribution to the SHE of the conventional spin current polarized in z direction in a disordered Rashba 2DEG. (Equations (2) – (7) in that paper).
- (13) S. V. Eremeev, I. A. Nechaev, Yu. M. Koroteev, P. M. Echenique, and E. V. Chulkov, Ideal two-dimensional electron systems with a giant Rashba-type spin splitting in real materials: Surfaces of Bismuth tellurohalides, Phys. Rev. Lett. 108, 246802 (2012).
- (14) M. Sakano, M. S. Bahramy, A. Katayama, T. Shimojima, H. Murakawa, Y. Kaneko, W. Malaeb, S. Shin, K. Ono, H. Kumigashira, R. Arita, N. Nagaosa, H. Y. Hwang, Y. Tokura, and K. Ishizaka, Strongly spin-orbit coupled two dimensional electron gas emerging near the surface of polar semiconductors, Phys. Rev. Lett. 110, 107204 (2013).
- (15) L. Wu, J. Yang, S. Wang, P. Wei, J. Yang, W. Zhang, and L. Chen, Thermopower enhancement in quantum wells with the Rashba effect, Appl. Phys. Lett. 105, 202115 (2014).
- (16) In the present paper we do not consider side-jump induced by spin-orbit scattering. A semiclassical treatment of this case can be found in P. M. Levy, H. Yang, M. Chshiev, and A. Fert, Spin Hall effect induced by Bi impurities in Cu: Skew scattering and side-jump, Phys. Rev. B 88, 214432 (2013).
- (17) C. Xiao, D. Li, and Z. Ma, Role of band-index-dependent transport relaxation times in anomalous Hall effect, Phys. Rev. B 95, 035426 (2017).
- (18) The reason why the side-jump AHE-SHE induced by band-structure spin-orbit coupling is defined as the sum of these three semiclassical contributions was detailed in Ref. Nagaosa2010, . Simply speaking, there are at least two motivations: one is the equivalence described in Ref. note-vertex, , the other is that all these three contributions belong to the disorder-induced interband-coherence effect (see Refs. Nagaosa2010, ; Kovalev2010, ; Hou2015, ; Xiao2017SOT, ).
- (19) C. Xiao, D. Li, and Z. Ma, Unconventional thermoelectric behaviors and enhancement of figure of merit in Rashba spintronic systems, Phys. Rev. B 93, 075150 (2016).
- (20) H.-Z. Lu and S.-Q. Shen, Extrinsic anomalous Hall conductivity of a topologically nontrivial conduction band, Phys. Rev. B 88, 081304(R) (2013).
- (21) J.-I. Inoue, T. Kato, Y. Ishikawa, H. Itoh, G. E.W. Bauer, and L.W. Molenkamp, Vertex Corrections to the Anomalous Hall Effect in Spin-Polarized Two-Dimensional Electron Gases with a Rashba Spin-Orbit Interaction, Phys. Rev. Lett. 97, 046604 (2006).
- (22) C. Grimaldi, E. Cappelluti, and F. Marsiglio, Off-Fermi surface cancellation effects in spin-Hall conductivity of a two-dimensional Rashba electron gas, Phys. Rev. B 73, 081303(R) (2006).
- (23) P. Wang, Y.-Q. Li, and X. Zhao, Nonvanishing spin Hall currents in the presence of magnetic impurities, Phys. Rev. B 75, 075326 (2007).
- (24) K Chadova, S. Wimmer, H. Ebert, and D. Kodderitzsch, Tailoring of the extrinsic spin Hall effect in disordered metal alloys, Phys. Rev. B 92, 235142 (2015).
- (25) A. Fert and P. M. Levy, Spin Hall Effect Induced by Resonant Scattering on Impurities in Metals, Phys. Rev. Lett. 106, 157208 (2011).
- (26) In the context of AHE induced by band-structure spin-orbit coupling, it has been established that the disorder-induced interband-coherence contribution (side-jump) calculated in the semiclassical Boltzmann theory is equivalent to the ladder vertex correction in the weak disorder limit to the bubble of the anomalous Hall conductivity in the non-chiral basis ( basis for the Rashba model). The present calculations suggest that this equivalence is also valid for the spin Hall conductivity of the conventionally defined spin current. In fact, this equivalence has been employed in the statement of Ref. Sinova2015, .
- (27) W. Kohn and J. M. Luttinger, Quantum Theory of Electrical Transport Phenomena, Phys. Rev. 108, 590 (1957).
- (28) I. A. Ado, I. A. Dmitriev, P. M. Ostrovsky, and M. Titov, Anomalous Hall effect with massive Dirac fermions, Europhys. Lett. 111, 37004 (2015).
- (29) J. M. Luttinger, Theory of the Hall Effect in Ferromagnetic Substances, Phys. Rev. 112, 739 (1958).
- (30) J. Shi, P. Zhang, D. Xiao, and Q. Niu, Proper definition of spin current in spin–orbit coupled systems, Phys. Rev Lett. 96. 076604 (2006).
- (31) N. Sugimoto, S. Onoda, S. Murakami, and N. Nagaosa, Spin Hall effect of a conserved current: Conditions for a nonzero spin Hall current, Phys. Rev. B 73, 113305 (2006).
- (32) K. Tsutsui and S. Murakami, Spin-torque efficiency enhanced by Rashba spin splitting in three dimensions, Phys. Rev. B 86, 115201 (2012).
- (33) C. Xiao, D. Li, and Z. Ma, Thermoelectric response of spin polarization in Rashba spintronic systems, Front. Phys. 11, 117201 (2016).
- (34) N. Zhang, Y. Wang, J. Berakdar and C. Jia, Giant spin-orbit torque and spin current generation in carriers at oxide interfaces, New J. Phys. 18, 093034 (2016).