Dissipationless Quantum Spin Current at Room Temperature

See allHide authors and affiliations

Science  05 Sep 2003:
Vol. 301, Issue 5638, pp. 1348-1351
DOI: 10.1126/science.1087128


Although microscopic laws of physics are invariant under the reversal of the arrow of time, the transport of energy and information in most devices is an irreversible process. It is this irreversibility that leads to intrinsic dissipations in electronic devices and limits the possibility of quantum computation. We theoretically predict that the electric field can induce a substantial amount of dissipationless quantum spin current at room temperature, in hole-doped semiconductors such as Si, Ge, and GaAs. On the basis of a generalization of the quantum Hall effect, the predicted effect leads to efficient spin injection without the need for metallic ferromagnets. Principles found here could enable quantum spintronic devices with integrated information processing and storage units, operating with low power consumption and performing reversible quantum computation.

Our work is driven by the confluence of the important technological goals of quantum spintronics (1, 2) with the quest of generalizing the quantum Hall effect (QHE) to higher dimensions. The QHE is a manifestation of quantum mechanics observable at macroscopic scales. In contrast to most transport coefficients in solid-state systems, which are determined by the elastic and the inelastic scattering rates, the Hall conductance σH in QHE is quantized and completely independent of any scattering rates in the system, where the transport equation is given by jα = σHϵαβ Eβ [jα and Eβ (α,β = 1,2) are the charge current and the electric fields respectively, ϵαβ is the fully antisymmetric tensor in two dimensions]. Though dissipative transport coefficients are expressed in terms of states in the vicinity of the Fermi level, the nondissipative quantum Hall conductance is expressed in terms of equilibrium response of all states below the Fermi level. The topological origin of the QHE is revealed through the fact that the Hall conductance can also be expressed as the first Chern number of a U(1) gauge connection defined in momentum space (3). Recently, the QHE has been generalized to four spatial dimensions (4). In that case, an electric field Eν induces an SU(2) spin current Math through the nondissipative transport equation Math, where Math is the t'Hooft tensor, explicitly given by Math and σ is a dissipationless transport coefficient. The quantum Hall response in that system is physically realized through the spin-orbit coupling in a time-reversal symmetric system. At the boundary of this four-dimensional quantum liquid, when both the electric field and the spin current are restricted to the three-dimensional sub-space, the dissipationless response is given by Math(1) This fundamental response equation shows that it is possible to induce a purely topological and dissipationless spin current by an electric field in the physical, three-dimensional space.

We consider a realization of this electric field–induced topological spin current in conventional hole-doped semiconductors. In a large class of semiconductors, including Si, Ge, GaAs, and InSb, the valence bands are fourfold degenerate at the Γ point (Fig. 1). The effective Luttinger Hamiltonian (5) for holes is given by Math(2) where Si is the spin-3/2 matrix. We take the hole picture, and reverse the sign of the energy. Good quantum numbers for this Hamiltonian are the helicity Math and the total angular momentum Math. This kinetic Hamiltonian is diagonalized in the basis where the helicity operator λ is diagonal and the eigenvalue is given by Math For a given wave vector k, the Hamiltonian (Eq. 2) has two eigenvalues, Math and Math forming Kramers doublets. They are referred to as the light-hole (LH) and heavy-hole (HH) bands. In semiconductors with zincblende structure, such as GaAs, inversion symmetry breaking causes an additional tiny splitting in the LH and HH bands. We can neglect the inversion symmetry breaking when the temperature is much higher than this splitting. The band structure of semiconductors deviates from the spherical to the cubic symmetry. We also neglect this effect for simplicity, because physics described below are not so much affected by it.

Fig. 1.

Approximate band structure of GaAs. We neglect the small splitting due to inversion symmetry breaking. We also neglect the anisotropy of the bands. The conduction band (CB) is twofold degenerate. The valence bands consist of the HH, the LH, and the split-off (SO) bands, each of which is twofold degenerate. We consider the p-GaAs, and the Fermi momentum for each band is labeled as Embedded Image and Embedded Image, respectively. The Fermi energy shown corresponds to n = 1019 cm3.

We shall consider the effect of a uniform electric field E. Our full Hamiltonian is thus given by H = H0 + V(x), where V(x) = eE·x, and –e is the charge of an electron. We assume that the split-off band is totally occupied. We first define a 4 by 4 unitary matrix U(k) which diagonalizes the kinetic Hamiltonian H0. U(k) is defined by U(k)(k·S) U(k) = k Sz. In the spherical coordinates where k = k(sinθcosϕ, sinθsinϕ, cosθ), U(k) can be expressed as U(k) = exp(iθSy) exp(iϕSz). Under this unitary transformation, the new Hamiltonian H̃ ≡ U(k) HU(k) becomes Math Math Eigenvalues of Sz physically describe the helicity Math in the original basis. The kinetic part H0 now becomes diagonal, in the representation where Sz is diagonal. Because x = iδk, the potential term becomes V(), where the covariant derivative is defined by = iδkà and à = –i U(kk U(k). Because à is a pure gauge potential, there is no curvature associated with it. Up to this point, the transformation is exact. We now consider adiabatic transport and make a corresponding approximation. As is usually assumed in the transport theory, we neglect the interband transitions, i.e., the off–block-diagonal matrix elements of à connecting the LH and HH bands. Then we arrive at a nontrivial adiabatic gauge connection A [supporting online material (SOM) text], which takes a block-diagonal form in the LH and HH subspace. Because each band is twofold degenerate, the gauge connection is, in general, non-Abelian. However, A has no matrix elements connecting the λ = 3/2 and λ = –3/2 states in the HH band, because the gauge field à only connects states with helicity difference Δλ = 0, ±1. Therefore, the non-Abelian structure is only present in the LH band. For simplicity of presentation, we shall first make an additional, Abelian approximation (AA), in which only the diagonal components in A are retained. Afterward, we shall give our final results, including fully the non-Abelian corrections.

Within the AA, A is a diagonal 4 by 4 matrix in the helicity basis. Because a band-touching point acts as a Dirac magnetic monopole in momentum space (6), each diagonal component of A(k) is given by that of a Dirac monopole at k = 0, with the monopole strength eg given by λ. The associated field strength is given by Math(4) The effective Hamiltonian takes the form Math(5) Henceforth, xi denotes a covariant derivative in momentum space: xi = Di = iδ/δkiAi(k). The definition of xi has changed by projecting the original Hamiltonian H onto the HH or LH band. Whereas Heff seems to be trivial, its nontrivial dynamics are revealed through the nontrivial commutation relations Math(6) Such a situation also happens in the Gutzwiller projection of the SO(5) model (7). It also resembles the nontrivial commutation relation between the position operators of a two-dimensional electron gas projected onto the lowest Landau level (8), where Fij = Bϵij, and B is the external magnetic field. This general algebraic structure, called “noncommutative geometry,” also underlies the four-dimensional QHE model (4). In our present context, the non-commutativity between the three-dimensional coordinates arises from the magnetic monopole in momentum space, and it is a natural generalization of the QHE to three dimensions.

The equation of motion for holes can be derived easily from Eqs. 5 and 6 as Math(7) The last term, proportional to Fij, is a topological term, describing the effect of the magnetic monopole on the orbital motion. It represents a “Lorentz force” in momentum space, making the hole velocity noncollinear with its momentum, in contrast to the usual situations. In fact, if we interchange the roles of x and k in this term, this term becomes the Lorentz force for a charged particle moving in the presence of a magnetic monopole in real space. This set of equations can be integrated analytically (SOM text), and the resulting trajectory is shown in Fig. 2. The hole motion in real space obtains a shift perpendicular to S. This shift is analogous to the deflection of a charged particle by a magnetic monopole in a direction perpendicular to the plane spanned by its position and velocity vectors (9). It causes a spin current perpendicular to both E and S. For example, for E parallel to the +z direction, the spin current for each band at zero temperature, with spin parallel to the x axis, flowing to the y direction is given by Math(8) which is obtained by summing contributions from all the filled states. Here, we assumed that the equilibrium momentum distribution is attained by the random impurity scattering that causes the charge relaxation. Eq. 7 describes only the ballistic motion, and scattering by random impurities would lead to additional contributions to the spin current. As one can see from the detailed discussions in the SOM text, these extrinsic effects are not only small, but also scale with a higher power of kFn1/3, where n is the hole density. Therefore, by plotting σs/n1/3 against n and extrapolating to the limit of n → 0, the constant intercept would uniquely determine our predicted dissipationless spin conductivity.

Fig. 2.

The real-space trajectory of the hole obtained by solving Eq. 7, The electric field E = e 1V is parallel to the +z direction. Due to the noncommutative relation between the components of the position operator xi, i.e., [xi, xj] = – iF ij, with Fij being the gauge curvature defined in the momentum space, the hole obtains a transverse velocity whose direction depends on the helicity Embedded Image, indicated by the thick arrows. The broken line is parallel to k.

It is worth noting that this AA becomes exact in zero-gap semiconductors, e.g., α-Sn. In this class of materials, the bottom of the conduction band and the top of the valence band correspond to the LH and HH bands in other semiconductors like GaAs. These two bands touch at k = 0. In this case, p-doping introduces holes only into the HH band, and the AA becomes exact.

The spin current induced by the electric field can also be understood in terms of the conservation of the total angular momentum Math. As remarked earlier, J commutes with H0. When E is parallel to the z direction, Jz also commutes with the potential. Therefore, substituting Math, we obtain Embedded Image(9) The second term, representing the torque, vanishes in our case because points along the z direction. The first term Math vanishes in usual problems; however, it does not in our case, due to the noncollinearity of the velocity and the momentum. Furthermore, the first term, describing the time derivative of the orbital angular momentum Math, is proportional to the spin current. The third term Embedded Image , describing the time derivative of the spin angular momentum S, can be easily evaluated from the acceleration equation in Eq. 7. Therefore, we see that the conservation of the total angular momentum in Eq. 9 directly implies the spin current in Eq. 8. The spin current flows in such a way that the change of L exactly cancels the change of S.

We now discuss the correction due to the non-Abelian nature of the gauge connection of the LH band. Remarkably, even though the gauge connection is non-Abelian, the associated field strength is Abelian and gives a correction factor of –3, compared with the AA (1012). The equation of motion modified accordingly agrees with that obtained by generalizing the wave packet formalism (13) to the non-Abelian case. This non-Abelian correction gives the following result for the spin current: Math(10) Here, we defined σs to have the same dimension as the electrical conductivity, to facilitate comparison. The spin current equation is rotationally invariant, with the covariant form given in Eq. 1, and is the central result of our paper. In contrast with similar effects (14, 15), this spin current has a topological character; the spin conductivity σs in Eq. 1 is independent of the mean free path and relaxational rates, and all states below the Fermi energy contribute to the spin current, where each contribution is determined purely by the gauge curvature in momentum space, similar to the QHE (3). Assuming the hole density n = 1019 cm3, the mobility of the holes at room temperature in GaAs is μ = 50 cm2/V·s (16), and the conductivity is σ = e n μ = 80 Ω1 cm1. On the contrary, the spin Hall conductivity σs in Eq. 10 is estimated as σs ∼ 80 Ω1 cm1, being of the same order with σ. For lower carrier concentration, σs becomes larger than σ; for n = 1016 cm3, we have σ = 0.6 Ω1 cm1 and σs = 7 Ω1 cm1. At finite temperature, Eq. 10 is modified only through the Fermi distribution function nλ(k). Because the typical energy difference between the LH and HH bands at the same wavenumber is about 0.1 eV, which largely exceeds the energy scale of the room temperature ∼0.025 eV, our predicted effect remains of the same order even at room temperature.

The nondissipative spin transport equation Eq. 1 does not violate the time-reversal symmetry T. Our microscopic Hamiltonian H, the electric field E, and the spin current are all T invariant. Therefore, the electric field and the spin current can be related by a T symmetric, dissipationless transport coefficient σs. This situation is to be contrasted with Ohm's law. As the charge current is odd under T, while the electric field is even, they can only be related by a T antisymmetric, dissipative transport coefficient, namely, the charge conductivity. One of the main objectives in quantum computing is to achieve reversible computation (17, 18). From the above analysis, we see that there is a fundamental difference between the ordinary irreversible electronics computation based on Ohm's law and the reversible spintronics computation based on Eq. 1. The time reversal symmetry property encoded in Eq. 1 could provide a fundamental principle for the reversible quantum computation.

This spin current is also useful for spin injection into semiconductors. Although effective spin injection is necessary for spintronic devices, it has been an elusive issue (2). Usage of ferromagnetic metals is not practical because most of the spin polarizations are lost at the interface due to conductivity mismatch between metal and semiconductor (19, 20). Spin injection from ferromagnetic semiconductors such as Ga1xMnxAs has been successful (2123). Nevertheless, Tc is at most 110 K for Ga1xMnxAs, still too low for practical use at room temperature. Thus, it is desirable to find an effective method for spin injection. The electric field–induced spin current serves as a spin injector, because it creates a spin current inside the semiconductor. One might worry that the short relaxation time τs ∼ 100 fs (24) of the hole spins. This shortness of τs is because the strong spin-orbit interaction in the valence band combines the relaxation of momenta and spins (25). Most of the efforts on the spintronics in GaAs have been focused on electrons in the conduction bands, which have much longer spin-relaxation time (∼100 ps) (26). Nevertheless, our spin current is free from such rapid relaxation of spins of holes, because it is a purely quantum mechamical effect with equilibrium spin-momentum distribution. Only when the spin-momentum distribution deviates from equilibrium, e.g., in spin accumulation at boundaries of the sample, does the rapid relaxation of hole spins become effective.

One can consider following experimental setups for detection of the constant spin supply from p-GaAs. When the electric field is applied along the z direction and the electric current Jz is induced, the sx spin current jyx will flow along the y direction. One possibility is to see the spin-dependent electric transport through a ferromagnetic electrode with the magnetization M along ±x direction attached to the positive y side of the sample. With a lead connecting this electrode and the other (negative y) side of p-GaAs as shown in Fig. 3A, one should see a change of the electric current I depending on the direction of M. The ratio of I when M is along ±x-direction, I(+x)/I(–x), is expected to be well larger than unity. For the ferromagnetic electrode, ferromagnetic metals are not efficient because of the conductance mismatch (20). Instead, ferromagnetic semiconductors will be suitable. Another possibility is to measure circular polarization of light emitted via recombination with electrons. This can be achieved by a similar experimental setup in (21), replacing (Ga,Mn)As by p-GaAs, where the quantum well structure of (In,Ga)As is sandwiched by p-GaAs and n-GaAs (Fig. 3B). The spin current injected along the y direction will be recombined with the electrons supplied from the n-GaAs in the (In,Ga)As quantum well.

Fig. 3.

Experimental setups for the detection of spin current induced by an electric field. (A) Detection by attaching a ferromagnetic electrode. The dependence of current I flowing into the electrode on the direction of the magnetization M is to be measured. (B) Detection by measuring the polarization of the emitted light in the quantum well of (In,Ga)As. The emitted light with circular polarization is shown as an open arrow. Right and left circular polarization will be switched when the direction of the external electric field is reversed.

When the system is not connected to the leads along the y direction, spins accumulate near the edges of the sample. This spin polarization can, in principle, be measured by the Kerr rotation. The spin distribution is determined by a balance between the spin current supply and the spin relaxation. At room temperature, because τs = 100 fs (24) is rather short, the area density of spin accumulation at the sample surface, jxyτs, is too small to be observed. However, there are several ways to make τs longer. One is to lower the temperature. Another is to inject the spin into n-GaAs through the p-n junction, as demonstrated recently (27, 28). The spin lifetime of electrons in GaAs is 100 ps, 103 times longer than that of holes. Thereby, the spin current can be detected by the Kerr rotation through surface reflection.

Supporting Online Material

SOM Text


References and Notes

View Abstract

Stay Connected to Science

Navigate This Article