Spin polarization induced by an electric field in the presence of weak localization effects

Spin polarization induced by an electric field in the presence of weak localization effects

Daniele Guerci Juan Borge Roberto Raimondi roberto.raimondi@uniroma3.it Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy Nano-Bio Spectroscopy group, Departamento de Fìsica de Materiales, Universidad del Paìs Vasco UPV/EHU, E-20018 San Sebastiàn, Spain
Abstract

We evaluate the spin polarization (Edelstein or inverse spin galvanic effect) and the spin Hall current induced by an applied electric field by including the weak localization corrections for a two-dimensional electron gas. We show that the weak localization effects yield logarithmic corrections to both the spin polarization conductivity relating the spin polarization and the electric field and to the spin Hall angle relating the spin and charge currents. The renormalization of both the spin polarization conductivity and the spin Hall angle combine to produce a zero correction to the total spin Hall conductivity as required by an exact identity. Suggestions for the experimental observation of the effect are given.

keywords:
Spin-orbit coupling, weak localization, 2DEG, spintronics
journal: Journal of LaTeX Templates

1 Introduction

Weak localization (WL) is the result of quantum interference corrections to the semiclassical theory of transportbergman1984 (); lee1985 (). It manifests itself in good conductors as a negative or positive correction to the electrical conductivity depending on the symmetry properties of the system. The functional form varies with the effective dimensionality of the sample, behaving as a square root of temperature in three dimensions and logarithmically in two dimensionsgorkov1979 (). In the latter case, the resummation of the logarithmic correction via the renormalization group leads eventually to the Anderson localization transition in dimensionsabrahams1979 (); efetov1980 (). In the presence of spin-orbit coupling (SOC), the correction is positive and hence manifests as an antilocalizing behaviorhikami1980 (). SOC affects WL because it yields a finite spin relaxation time, which introduces a cutoff in the logarithmic singularity associated with the so-called triplet channel of the particle-particle ladder, known as the Cooperon. Since the singlet and the triplet channels contribute to WL with opposite signs, the elimination of the triplet leaves the singlet alone, which then produces the antilocalizing behavior. In metallic conductors and doped semiconductors SOC was traditionally attributed to the electric field of impurities, which do not affect the nature of the electron eigenstates. In the last two decades, however, the two-dimensional electron gas (2DEG) has become one of the most analyzed model systems for electrical transport, due to the possibility of realizing it in semiconducting systems, and more recently at metallicRojasSanchez2013 () and oxidesCaviglia10 () interfaces. The realization of the 2DEG leads to the breaking of inversion symmetry with respect to the axis, say the z axis, perpendicular to the 2DEG plane, say the x and y plane. In these circumstances, in the presence of a finite spin-orbit interaction, Bychkov and Rashba have proposed a model HamiltonianRashba84 (), which, besides the standard effective-mass kinetic energy term, contains a spin-orbit coupling term linear in momentum

(1)

where is the vector of the components of the momentum operator, is the effective mass and a SOC constant with and the standard Pauli matrices. The Rashba Hamiltonian Eq.(1) has been extensively studied over the last twenty years, especially aiming at the development of new spintronic functionalitiesAwschalom2007 (). In this respect the spin Hall effect (SHE)Dyakonov71 (); Hirsch99 (); Zhang00 (); Murak04 (); Sinova04 () and the current-induced spin polarization effectEdelstein1990 (); Lyanda-Geller89 () (known also as the Edelstein or inverse spin-galvanic effect) have been the focus of an intensive dedicated research. These effects, whose precise definition will be given later on, manifest due to the coupling of charge and spin degrees of freedom and hence introduce, besides the standard electrical conductivity, new transport parameters. These are defined as the linear coefficients relating the spin polarization and the spin current to the applied electric field

(2)

where the double brackets indicate the quantum and statistical average. and are referred to as the spin polarization or Edelstein and the spin Hall conductivities, respectively. As for the electrical conductivity, these transport parameters can be studied with the well-known impurity technique. One advantage of this techniqueAbrikosov:QFT (), based on standard diagrammatic perturbation theory, is the appearance of the semiclassical Drude-Boltzmann theory of transport at the leading approximation in an expansion of the small parameter , where and are the Fermi energy and the elastic scattering time, the only two parameters characterizing a disordered Fermi gas. In such an expansion, WL arises in the next-to-leading approximation in the expansion in .

WL effects in the presence of the Rashba SOC described by Eq.(1) have been analyzed by several authors, most of the attention having been focused on the electrical conductivity onlyIordanskii_1994 (); Edelstein95 (); Skvortsov_98 (); Lyanda98 (); Wenk_10 (); Araki_2014 (). It is the aim of the present work to extend this analysis to the other transport parameters mentioned above, whose experimental study has developed considerably in the last few yearsJungwirth2012 (); RojasSanchez2013 (). We find that and the spin Hall angle acquire logarithmic corrections which can be absorbed in terms of the renormalization of the scattering time appearing in the electrical conductivity . We emphasize that is not the full spin conductivity which would be measured in an experimentShen2013 (). As will be shown in the next Section, can be expressed in terms of and . The renormalizations of both and compensate in such a way that has no correction as expected on general argumentsDimitrova05 ().

The plan of the paper is as follows. In the next Section we introduce the disordered Rashba model and review the theory of and to the leading order in the parameter within the impurity technique. This is necessary to prepare the ground for the following Sections. Section III deals with the WL localization corrections in the presence of the Rashba SOC. The evaluation of the electrical conductivity is reviewed as an example. Section IV presents the calculation of the WL corrections to and . Section V provides a discussion of the results obtained, whereas technical points of the calculations are given in the appendices at the end of the paper. From now on, if not otherwise specified, we will work in natural units .

2 The disordered Rashba two-dimensional electron gas at leading order in

In the presence of scattering from impurities, the 2DEG Hamiltonian of Eq.(1) acquires an additional random potential term defined by the averages

(3)

where and are the coordinate operators, the two-dimensional density of states and the elastic scattering time. At leading order in the expansion parameter , the selfenergy is given by the selfconsistent Born approximation

(4)

where denotes the retarded and advanced Green functions. As discussed in Raimondi02 (); Raimondi05 (), in the presence of Rashba SOC the Green function has a nontrivial structure in spin space, whereas the selfenergy remains diagonal, , . Explicitly we have:

(5)

with

(6)

The Edelstein (EC) and spin Hall (SHC) conductivities are defined in terms of the spin polarization and spin Hall current induced by an applied electric field taken along the x axis for definiteness’s sake . The corresponding Kubo formulae are

(7)

and

(8)

where the bare vertices , and , denote the operators for spin polarization, spin current and charge current, respectively. The evaluation of the response functions (7-8) involves the standard bubble diagrams of the Green function lines obtained by the selfconsistent Born approximation (4) decorated by the insertion of the impurity ladder. This corresponds to the inclusion of the so-called vertex corrections, which lead to renormalized verticesRaimondi02 ().

The expression (7) for the EC becomes

(9)

where the vertex renormalization can be attributed either to the left spin vertex or to the right current vertex and we have dropped the dependence on the frequency argument of the Green function. In the former case, by using the renormalized spin vertex indicated by a capital letter , one obtainsEdelstein1990 ()

(10)

where the subscript in indicates the lowest order in the parameter . We have defined the parameter and introduced the D’yakonov-Perel relaxation time , being the Fermi momentum in the absence of the Rashba SOC. The model has two small parameters and with , in terms of which the above paramater can be expressed.

Similarly, for the SHC one has the expression

(11)

By using the renormalized spin current vertex Tse_IntrinSH_PRB06 (), one gets

(12)

where we indicated by the SHC in the absence of vertex corrections as first computed in Ref. Sinova04 (). It is then not difficult to see that the insertion of the result (10) into (12) gives a vanishing SHC. This latter result is actually expected following the argument derived by Dimitrova Dimitrova05 (). The commutation relation allows us to write

(13)

which, in stationary circumstances, implies .

The expression (12) shows that the vertex corrections for the SHC are associated to the EC. This connection between the two effects acquires a more transparent meaning by adopting the SU(2) gauge-field description of the SOCRaimondi_AnnPhys12 (). In such a picture, the SOC is expressed in terms of a non-Abelian gauge field , with and Tokatly_Color_PRL08 (); Gorini10 (). The first consequence of resorting to this language is the appearance of an SU(2) magnetic field , which arises from the non-commuting components of the Bychkov-Rashba vector potential. As in the normal Hall effect we have a spin Hall drift component of the spin current which can be described as (assuming )

(14)

with . In addition to the drift current, there is also a “diffusion current” due to spin precession around the effective spin-orbit field. Within the SU(2) formalism this current arises from the replacement of the ordinary derivative with the SU(2) covariant derivative in the expression for the diffusion current. The SU(2) covariant derivative, due to the gauge field, is

(15)

with a given quantity being acted upon. The normal derivative, , along a given axis is shifted by the commutator with the gauge field component along that same axis. As a result of the replacement diffusion-like terms, normally proportional to spin density gradients, arise even in uniform conditions and the diffusion contribution to the spin current turns out to be

(16)

where is the diffusion coefficient. In the diffusive regime the full spin current can thus be expressed as

(17)

with the spin Hall angle given by and . Relation Eq.(17), which implies the diffusive limit, coincides with the expression (12) when , i.e. when the spin-orbit induced spin splitting is much smaller than the disorder broadening of the levels.

3 The weak-localization correction to the electrical conductivity in the presence of Rashba spin-orbit coupling

In this Section we consider the WL corrections, which arise to first order in . In next Subsection we review the evaluation of the Cooperon propagator, whereas in the following one we apply it to the case of the electrical conductivity. The spin polarization and spin current response are considered in the next Section.

3.1 The structure of the Cooperon in spin space

Figure 1: Crossed ladder, also known as Cooperon. The indices and are the incoming spin indices, while and are the outgoing ones. Furthermore, the momentum is the sum of the initial and final external momenta in each Green function line.

In the absence of SOC, the Cooperon obeys the diffusion equation, which in momentum space reads

(18)

where is the sum of the initial and final momenta in each Green function line, as shown in Fig. 1. We remind that the particle-hole series of the maximally crossed diagrams can be transformed in the particle-particle ladder by reversing one of the two Green function lines and the associated momentum. The momentum corresponds then to the total momentum of the particle-particle pair of the Cooperon propagator. The SU(2) gauge-field point of view suggests that one can obtain the equation for the Cooperon by a minimal substitution procedure on the momenta of both Green function lines, i.e. and , where is the SU(2) spin-dependent vector potentials introduced in the previous Section. The equation for the Cooperon becomes then

(19)

where is the total spin of the two particles ( refers to the upper/lower Green function line). The Cooperon acquires a matrix structure acting in the four-dimensional Hilbert space resulting from the combination of the two spin one-half spaces. It is convenient to use the triplet-singlet basis , in terms of which the Cooperon becomes block-diagonal. Since for the singlet the total spin of the pair is zero , the one-dimensional corresponding block has the same form as in the absence of SOC. The triplet three-dimensional block can be obtained from Eq.(19) by expanding and inverting the square in the left hand side. One sees that, besides the standard diffusive terms going as , which are diagonal in spin space, linear-in- terms arise when considering the double product. Finally, the terms arising from the square of the spin-dependent vector potential give rise to -independent terms. The latter are also diagonal in spin space and yield a finite spin relaxation rate. The full expression for the Cooperon is obtained by inverting Iordanskii_1994 (); Edelstein95 (); Skvortsov_98 (); Lyanda98 (); Wenk_10 ()

(20)

where is the angle between and the -axis. According to the matrix form (20), the singlet state decouples from the triplet states and remains gapless, whereas the triplet sector acquires a gap proportional to . Furthermore, the linear-in- terms provide the mixing between the different triplet channels. These results are due to the D’yakonov-Perel spin relaxation mechanism, which mixes the triplet states and kills their singular contribution.

We notice that the previous expression (20) is valid in the regime . In A we derive the exact form of the Cooperon equation (65), whose solution tends to expression (20) as . However, the approximate solution (20) possesses all the physical features due to the presence of Rashba SOC and simplifies the calculations of the quantum corrections. Because of that, in the rest of this work we will use the Cooperon expression (20) to calculate the weak localization corrections.

3.2 The electrical conductivity to order

(a)
(b)
(c)
Figure 2: Maximally crossed diagrams giving the weak localization corrections to the transport coefficients. In each diagram the ladder summation takes into account the renormalization of both vertices of the bubble. The crossed lines indicate the crossed ladder, also known as Cooperon. In diagram (a) the integration over the momentum of four Green functions yields the so-called Hikami box, whose full expression is obtained by adding the integrations over and of diagrams (b) and (c).

First we note that the resummation of the maximally crossed diagrams for the evaluation of the electrical conductivity gives rise to three diagrams, depicted in Fig. 2. The following discussion is quite general and applies also to the case of EC and SHC to be discussed in the next Section.

The weak localization correction to the electrical conductivity reads

(21)

where the subscript refers to each diagram in Fig. 2. As can be understood by looking at Fig. 1(a), the first diagram gives the following contribution

(22)

with

(23)

where is the direct product basis with , and

(24)

Notice that the renormalized charge current reads Raimondi02 (). In order to evaluate Eq.(24) we make three approximations: first, we cut off the sum over ; second, we expand for small ; third, we neglect the weak dependence of . The natural upper cut-off of the sum over is provided by the inverse of the mean free path, i.e. the length where the diffusive behaviour sets in. On the other hand, the lower cut-off of the sum is given by the inverse of the maximal size of a loop allowed to contribute to the coherent backscattering process. Then, for a system of linear size the lower cut-off is . For small the Cooperon is strongly depending on and gives a divergent contribution which is not modified by the weak dependence of the weight factors . Eventually, we arrive at the following expression for the conductivity correction given by the diagram in Fig. 1(a)

(25)

with the weight factors

(26)

and the -integrated Cooperon matrix elements

(27)

With the approximations introduced above, the contribution , given by the diagram of Fig. 1(b), is

(28)

with the weight factors

(29)

Finally for we find

(30)

where the weight factors can be expressed in terms of those of

(31)

Upon integration over the directions of , taking into account the Green function structure of Eq.(2), one discovers that several components of given in Eq.(26) vanish. A similar analysis can be done for the weight factors of Eq.(29) when considering integration over the directions of and . As a result, the only non vanishing weight factors are

(32)

and

(33)

where we have neglected terms beyond the leading order in and the first order in . By inserting the weight factors (32) in Eq. (25), the quantum correction, given by the diagram of Fig.1(a), is

(34)

The quantum correction, given by the diagrams of Figs. 1(b) and 1(c), is

(35)

To our knowledge the expression for diagrams of Figs. 1(b) and 1(c) has not been given before. In agreement with Edelstein95 () the contribution of the two diagrams 1(b) and 1(c) vanishes in the absence of Rashba SOC and does not exhibit scaling behavior for . So it can be ignored and the quantum correction to the electrical conductivity is given by Iordanskii_1994 (); Skvortsov_98 (); Lyanda98 (); Araki_2014 () .

Figure 3: The quantum correction to the electrical conductivity as a function of the variable for four different values of the ratio from bottom to top. In the axis we set , where .

The above correction (34) can be interpreted as a sum of a localizing contribution from the triplet sector and an antilocalizing contribution from the singlet channel. Therefore, for , the triplet modes acquire a gap that reduces their contribution and the total correction changes its sign and becomes antilocalizing at a certain . The critical value , which defines the crossover point between WL and WAL, depends on the ratio . In particular, as the ratio gets larger, the value decreasesSkvortsov_98 (); Araki_2014 (). This behavior is evident by looking at Fig. 3; for the quantum corrections are localizing, whereas for the quantum correction tends to

(36)

which implies the enhancement of the electrical conductivity and the diffusion coefficient.

4 The weak-localization corrections to the EC and SHC

The equation (21) for the quantum correction to the electrical conductivity can easily be generalized to the spin transport coefficients by substituting in the weight factors (26) and (29) the left charge current vertex with the renormalized spin density vertex or the renormalized spin current vertex. The correction to the Edelstein conductivity then reads

(37)

with weight factors

(38)

and

(39)

Notice that in the right vertex it appears the renormalized charge current vertex , whereas in the leading order expression (9) the bare charge vertex appears. The non-vanishing weight factors are:

(40)

and if . In the above equation we have defined . Since the weight factors are diagonal in the indices and , only the diagonal elements contribute to the correction to the Edelstein conductivity Eq.(37). If we take into account only the singlet channel, which is the only mode with scaling behavior for , then all the diagonal elements are related: with and

(41)

Then, by performing the trace over the product of Pauli matrices in Eq.(37), one obtains

(42)

from which the quantum correction to the Edelstein conductivity reads

(43)

The above contribution increases the leading order value

(44)

which shows that the logarithmic correction to the Edelstein conductivity can be absorbed in terms of the renormalization of the elastic scattering time. Thus, Rashba SOC provides the spin relaxation mechanism, which causes the change of sign of the weak localization corrections and does not introduce a new scaling parameter.

For the evaluation of the spin Hall conductivity, we rely on the decomposition shown in the first line of Eq.(12), where here we replace the bare charge current verices with the renormalized one and insert a series of maximally crossed impurity lines between the two vertices. Hence, the resulting expression will have a Cooperon connected to a Hikami box with modified vertices. The resulting expression reads

(45)

where

(46)

Notice that the response function indicated by contains the renormalized charge current vertex in contrast with which has the bare charge current vertex. This is the reason why we only find the term and we don’t find in Eq.(45). The response function describes the influence of the WL corrections on the spin Hall angle.

The weight factors are

(47)

and

(48)

The only non-vanishing weight factors are:

(49)

and if . In the above equation we have defined , since the contribution of the first diagram of the Hikami boxes is equal to zero, . Finally, by replacing the weight factors with the in Eq.(42), the quantum correction to the drift component reads

(50)

where we have only taken the logarithimic contribution given by the singlet channel. It is then evident that the insertion of the results (50) and (43) into (45) gives the vanishing of the weak localization correction to the spin Hall conductivity

(51)

as required by the relation (13) in the static limit.

5 Discussion

In this Section we develop an interpretation of the structure of the weak localization results of the previous Section. By considering Eq.(17) connecting spin current and spin density, we assume for the WL corrections the structure

(52)

where the variation has been applied to both the coefficients connecting the averaged values of the observables and to the average of the observables themselves. In the limit , by noticing that , we identify the second term in the right hand side of the above equation with the second term in the right hand side of Eq.(51). In the same limit, i.e. keeping terms up to order , we see that the first and last term of Eq.(52) must necessarily cancel among themselves. Hence, one can identify the third term of Eq.(52) with the first term of Eq.(51), thus making clear the reason for introducing the response function describing the quantum correction to the spin Hall angle.

The direct measurement of the EC or its inverse is possibleRojasSanchez2013 (); ganichev2002 () and a possible test of the theory presented here would be the study of the EC as a function of an applied magnetic field perpendicular to the 2DEG. On the other hand, the SHC vanishes in the Rashba disordered 2DEG. It is then important to extend this analysis to a case where the SHC is finite, as for instance when also extrinsic SOC from impurities is presentRaimondi09 (). This however is beyond the scope of the present work.

Aknowledgements R.R. thanks Valentina Brosco and Lara Benfatto for discussions.

Appendix A The Cooperon in the presence of Rashba spin-orbit coupling

To describe the sum of maximally crossed diagrams, we have introduced the so-called Cooperon, corresponding to the diagram depicted in Fig. 1. In the main text we have derived an approximate solution to the Cooperon equation, by using an argument proposed by Kettemann and Wenk Wenk_10 (). In this appendix we want to obtain the exact derivation of the Cooperon.

As can be understood by looking at Fig. 1, the equation for the Cooperon can be written as

(53)

where and

(54)

where , , and . Since we expect the Cooperon to be large in vicinity of backscattering , we set and expand in powers of Q up to the quadratic order:

(55)

where the expansion coefficients are given by

(56)

Zeroth order:

By symmetry due to the integration over the angle of , there are only three integrals different from zero , and . In particular one has