Electronic double refraction due to the Rashba effect: Analytical and numerical results

Peng-rui SHAO , Wen-ji DENG

Front. Phys. ›› 2007, Vol. 2 ›› Issue (1) : 44 -50.

PDF (531KB)
Front. Phys. ›› 2007, Vol. 2 ›› Issue (1) :44 -50. DOI: 10.1007/s11467-007-0015-6
Condensed Matter and Materials Physics
Electronic double refraction due to the Rashba effect: Analytical and numerical results
Author information +
History +
PDF (531KB)

Abstract

By analogy with the classic effect of the double refraction of light, we investigate the relevant effect of an electron entering from the Non-Rashba region to the Rashba region in two-dimensional systems. It is shown that the effect of electronic double refraction is determined by a combined parameter γ=mλ F α/2π2, rather than both the Rashba coefficient α and wavelength λ F  of a Fermi electron, separately. For the case of normal incidence, the analytical expressions for the wavefunction of the electron are presented; it is predicted that the Rashba spin-orbit coupling can induce a current perpendicular to the normal incident direction of the electron. Moreover, the general case of incident electron with any given momentum and spin state are studied numerically in detail, including the abrupt changes of spin direction and the two-step characters for reflection.

Graphical abstract

Keywords

spin-orbit coupling interaction / spin polarized current / electronic double refraction

Cite this article

Download citation ▾
Peng-rui SHAO, Wen-ji DENG. Electronic double refraction due to the Rashba effect: Analytical and numerical results. Front. Phys., 2007, 2 (1) : 44-50 DOI:10.1007/s11467-007-0015-6

登录浏览全文

4963

注册一个新账户 忘记密码

1 Introduction

Since Datta and Das proposed, as a spintronic analogy of electro-optical modulators, the original idea of spin-field-effect transistor (s-FET) in 1990 [1], it has been widely believed that the Rashba effect [2, 3] may play an important role in producing and controlling the spin-polarized currents in nonmagnetic semiconductor devices. This belief had inspired not only the revival of interest about the physics of spin-orbit interaction in solids [414], but also the designs of device setups based upon the Rashba spin-precession [1526], although none of them has been realized in the laboratory because of the difficulty of spin-polarized injection from ferromagnetic contacts to semiconductors [2729].

Recently, an alternative way to generate spin-polarized currents in heterostructures by means of the spin-double refraction effect is explored [3035]. It is suggested that a lateral interface connecting two regions with different strength of Rashba spin-orbit coupling interaction can be used as a spin polarizer for electrons in two-dimensional semiconductor heterostructures, and the incident wave of the electron can be split as two beams of electronic wave, each of which is completed spin polarized. The transmission and reflection waves strongly depend on the initial states of electrons, and can be manipulated by the strength of Rashba coupling.

In this work, we carry out, compared with the well-known effect of the double refraction of light [36], some investigations on the phenomena of electronic double refraction occurring at an interface between the Non-Rashba region and Rashba region of two dimensional electron gas. Section 2 is devoted to describe our model system and present some basic formulae. In Section 3, we present some results about the characters of electronic double refraction by analytical investigation and numerical simulations. Some further discussions and conclusions are presented in Section 4. For completeness, an appendix is added to discuss the problem of cross-term currents in the Rashba region.

2 Model and formalism

We consider a two-dimensional system that contains an interface between the Non-Rashba region and Rashba region with spin-orbit coupling interaction, as depicted in Fig. 1. The Rashba effect describes the influence of a large local electric field arising from the structure inversion asymmetry (SIA), which can also be modulated by the applied voltage across the heterostructure [4, 5, 7]. Under the approximation of effective mass, the Hamiltonian of an electron moving in such a system is always taken as a simple form [9], i.e.

(1)H=12m*(p^x2+p^y2)+α(σxp^yσyp^x)

here p^x,p^y and m* represent the momentum operator and effective mass of the electron, respectively; =h/2π denotes the Planck’s constant. The symbol α stands for the Rashba coefficient measuring the strength of Rashba effect, which is assumed to be zero in the Non-Rashba region of x < 0 and a positive constant in the Rashba region of x > 0; σx,σy and σz are the well-known Pauli matrices.

The Schrödinger equation for an electron described by the Hamiltonian of Eq. (1) can be written in a scaleless form as:

(2)iΨt=[k^k^2ασ(n×k)]Ψ

where k^i, α=m*lα/2, σ is the vector of Pauli matrices and n denotes the unit vector of the positive direction of z axis. The parameter l, which is defined as the unit of length, can be chosen as any characteristic length of the system. Because only the electrons at the Fermi surface are taken into account in this work, we define the wavelength of Fermi electrons as length unit, i.e., l=λF. The units of other relevant quantities, such as energy, momentum, wave vector and time should be E˜=2/m*λF2, p˜=/λF, k˜=1/λF, and T˜=/E˜, respectively. In any uniform region with constant Rashba coefficient, an explicit expression for the density of probability current can be derived from Eq. (2), which is

(3)J=i2(Ψ)Ψi2Ψ(Ψ)+α Ψ(n×σ)Ψ

For an incident electron with any given momentum k0= (k0x,k0y) and spin state |e0, the energy ω=k02 and momentum component k0y of the electron are conserved, the wavefunction of eigenstates can be assumed a general form as Ψ=ei(k0yyωt)ψ. The wavefunction of the electron in the Non-Rashba region can always take the form of

(4)ψI=|e0eik0xx+r|ereik0xx,    for  x<0

here |er describes the spin-state of reflected wave, and r denotes the reflection coefficient.

In a homogeneous region with constant Rashba coefficient, the eigenvalues of momentum k=(kx, k0y) for state ψ=eikxx|e can be drawn from the time-independent Schödinger equation [k22ασ(n×k)]|e=k02|e, which yields two eigenvalues of k, i.e., k±2k02=± 2α k±, and k±=(k02+α2)1/2±α. In the presentation of σz, the corresponding spinors |e± can be generally written out as:

(5)|e±=1N±((k0y+ik±x)         k±)

here k±x=(k±2k0y2)1/2 and N± denote the x components of wave-vectors of the two refracted waves and normalized constants, respectively. The other two modes with k±x= (k±2k0y2)1/2 are ignored because no reflected wave exists in the right region for the obvious reasons of physics. Then the function for diffracted waves can be assumed as:

(6)ψII=t+|e+eik+xx+t|eeikxx,    for  x>0

Noting that k+x=(k+2k0y2)1/2 stays real, the expression k+ in Eq. (5) can be simplified as:

(7)|e+=(ei(ϕ++π/2)cosπ4       sinπ4)

here ϕ+ denotes the refraction angle of k+ wave, and cos ϕ+=k+x/k+. In contrast, there is a critical angle ϕc= arcsin (k/k0) for the k mode. If kk0yk, the incident angle is not greater than the critical value, then the relevant spin state is similar to that of k+ mode, i.e.,

(8)|e=(ei(ϕπ/2)cosπ4       sinπ4)

and ϕ denotes the refraction angle of k wave, and cosϕ =kx/k.Otherwise kx=ik˜x and k˜x(k0y2k2)1/2, which means that the k mode becomes an evanescent wave, and the corresponding refraction angle is defined as ϕ=π/2. The spin state defined in Eq. (8) should be modified as:

(9)|e=1N(k0yk˜x      k)

with normalization constant N=[2k0y(k0yk˜x)]1/2. The spin directions e+ and e of double refracted waves are always perpendicular to the direction of the momentum k±=(k±x,k0y), i.e., φ+=ϕ++π/2 and φ=ϕ+π/2, as shown in Fig.1.

According to the formula of Eq. (3), the density of probability current of the incident, reflected, and refracted waves defined in Eqs. (4), (6) can be found as J0 = k0, Jr = |r| 2 kr, and J±=|t±|2(1α/k±)k± (if kx becomes a pure imaginary number, then Jx = 0). Therefore, the reflecting coefficient RJrx/J0x and the two refracting coefficients T±J±x/J0x can be expressed as

(10)R=|r|2,    T±=(1+γ2)1/2|t±|2cosϕ±cosϕ0

here we introduce a new parameter γa/k0. In appendix A, we will deal with the problem of cross-term currents between the incident wave and reflected wave, and between the two refracted waves. It will be shown exactly that R+T++T=1.

To solve the Schrödinger equation completely, we still need another two boundary conditions:

(11)Ψ|0+=Ψ|0,    Ψx|0+= Ψx|0+iα σyΨ|0

The latter can guarantee the conservation of probability current [9]. By inserting the wavefunctions of Eqs. (4), (6) into Eq. (11), we can derive a group of equations to find the amplitudes t+, tand the reflected wave as below:

(12)  (n+k+xk++cosϕ0γσy)t+|e++(nkxk+cosϕ0γσy)t|e=2cosϕ0|e0

(13)r|er=t+|e++t|e|e0

where n± k±/k0 = (1+γ2)1/2 ± γ, φ0 and θ0 denote the spin direction of the incident electron. If kx is real, then k±x/k±=cosϕ±, and n±sinϕ±=sinϕ0, which means that the Snell’s law of refraction holds for double refraction, and that the parameters n+ and n play the role of refraction index for the electron.

It is worthwhile to emphasize, based upon the discussion above, that the effect of electronic double refraction is determined by the parameter γ=m*λFα/2π2, a combination of the Rashba coefficient α and wavelength λF of the incident electron, rather than both of them separately.

3 Reflection and double refraction

A beam of light passing through some uniaxial crystals such as calcite or quartz, is split into two rays. The two waves are each plane polarized, and travel with different wavelengths, or have different indices of refraction [36].

Similarly, the matter wave of an electron entering the Rashba region from the Non-Rashba region can also be doubly refracted; one is the k+ mode and the other is k mode. They are each spin polarized and described by different refraction indices n+ and n.

The quantity of double refraction is depicted by the difference of two refraction indices. For a calcite crystal, the two indices of ordinary light and extraordinary light are no = 1.66 and ne = 1.49, respectively; and then Δn=0.17.

When taken into account, a typical semiconductor system, such as InAs system [10], the relevant parameters are chosen as typical values [26] α = 0.6 eV·Å, EF = 99.54 meV, m*=0.05 m0, then λF=174 Å, kF = 3.62×10–2 Å, α= 0.683 and γ0.1.The difference of refraction indices Δnn+n0.20, which should be large enough for the occurrence of double refraction.

In principle, all information about the double refraction of an electron entering the Rashba region, as depicted in Fig.1, can be dug out from Eqs. (12) and (13). The next two subsections are devoted to study some typical cases with different incident angles.

3.1 Case of normal incidence

For the case of normal incidence, i.e., k0y = 0, and e+ = j, e = –j, the problem can be solved analytically. The reflection wave

(14)r|er=(1+γ2)1/21(1+γ2)1/2+1|e0

undergoes a phase shift of π radian, but its spin state doesn’t change. The transmission wave

(15)ψII=2(1+γ2)1/2+1|χII

with |χII=eik+xj|e0|j+eikxj|e0|j behaves as a spin precession around the –y axes with a space period Lx=lπ/α=2mα/2.

According to Eq. (10), the reflection and refraction coefficients can be calculated analytically, i.e.,

(16)R=[(1+γ2)1/21(1+γ2)1/2+1]2

(17)T±=2(1+γ2)1/2[(1+γ2)1/2+1]2(1±sinθ0sinφ0)

the reflection is independent of the spin state of the incident electron, but both of the two refraction coefficients T+ and T depend on it.

In addition, the spin-orbit interaction can induce a current perpendicular to the normal incident direction of the electron. The y component of current

(18)Jy=Jy0sin (2α x+δ)

oscillates periodically in the Rashba region with amplitude

(19)Jy0=α(1sinθ02sinφ02)1/2[21+(1+γ2)1/2]2

and phase shift

(20)δ=arcsin[cosφ0sinθ0(1sinθ02sinφ02)1/2]

We predict that this kind of persistent current of charge can be observed in experiment.

3.2 Case of oblique incidence

The case of oblique incidence can only be studied numerically. For a practical system of InAs semiconductor [12], we find that γ0.1 and R104 following Eq. (10), which means a very weak reflection. To show more explicitly the effect, we should assume larger values of the Rashba constant α and/or lower energy values of EF to reach large enough values of parameter γ, for example γ=0.2, 0.5, 0.7, 1.0.

For the transportation of light wave in uniaxial crystals, a plane wave can always be treated as a superposition of a couple of orthogonal polarized plane waves [36], i.e., the ordinary ray (o-ray) and the extraordinary ray (e-ray). However, the spin eigenstates of an electron transporting in the Rashba region are not orthogonal, the two modes of k+ and k are mixed together by the spin-orbit coupling interaction, and we cannot deal with an incident wave of electron as two independent polarized waves. As a prominent consequence, the spin state of the reflected wave undergoes an abrupt change at the interface x = 0 in a complex way on the initial spin states.

The inner production ere0 is adopted to represent the change of spin direction of reflected electron, and the unit vector e0 and er of the spin direction of incident and reflected electron can be calculated as e0=e0|σ|e0 and er=er|σ|er, respectively. Although the spin direction remains unchanged after reflection for the normal incident electron with any spin state, it does alter dramatically for the oblique incident electron. For example, we do numerical calculations about the inner production er·e0 of an electron with incident angle ϕ0=π/6, this function varies with the polar angle θ0 and azimuth angle φ0 of the spin direction of the incident electron as shown in Fig. 2. In contrast with er·e0≡1 for a normal incident electron, the function for an oblique incident electron can hold er·e0 = 1 only for two special initial spin states.

To demonstrate more clearly the laws of reflection and refraction, we also calculate the coefficients of reflection and refraction for the incident electron with different incident angles and spin states. The three plots in Fig. 3 show the typical way of the coefficients depending on the initial state of an incident electron for a given incident angle ϕ0= π/6 and system parameter γ=0.7. Compared with Eq. (16) and Eq. (17), the formulae for normal incidence, the reflection coefficient R is no longer a constant, while the two refraction coefficients T+ and T vary with θ0 and φ0 in a way similar to that described by Eq. (17). We have also confirmed the equality R+T++T=1 by calculating the quantity T++T, which compensates the variation of R exactly.

The refraction coefficients T+, T and reflection coefficient R averaged over all possible spin states with equal probability are some other interested quantities. As shown in Fig. 4, all of them are monotonic functions of the incident angle ϕ0. Because there exists a critical angle ϕc=arcsinn for the k mode, both of the reflection coefficient R and refraction coefficient T behave like two-step functions. According to the definition n=(1+γ2)1/2γ after Eq. (13), we obtain the corresponding critical angles ϕc/π=0.31, 0.21, 0.17 and 0.14 for different parameters γ=0.2, 0.5, 0.7 and 1.0. It can be seen obviously from Fig. 4 that, the total reflection for k mode, and the abrupt changes of the curves of T+ and R, consistently occurs at these incident angles; however, the coefficient T+ smoothly decreased with the incident angle ϕ0.

4 Discussions and conclusions

In this work, we systematically study the effect of double refraction of an electron occurring at an intersurface between two half infinite regions, from the Non-Rashba (normal metal) region to the Rashba (semiconductor) region.

Similar to the well-known double refraction effect of light, an incident wave of electron with given momentum and spin state partially enters the Rashba region, and is divided into a couple of completed spin polarized currents. Both of the effects for light and electron share a common origin that, for any given incident wave, the two modes of refraction waves have different refraction indices, i.e., neno for light and n+ > n for electron. Otherwise, the laws of double refraction for electron due to Rashba spin-orbit coupling interaction re more complicated than that of the light effect. This distinction is simply attributed to the fact that the spin states of k+ and k modes for electron in the Rashba region are not orthogonal, but the o-ray and e-ray in uniaxial crystals are always polarized perpendicularly.

We derive a set of basic equations, i.e., Eqs. (12), (13), to study the electronic double refraction at the single interface between the Non-Rashba region and Rashba region. It is shown explicitly that this effect is determined by a combined parameter γ=m*λFα/2π2 instead of two independent parameters α and αF. Therefore, we suggest a new way to enhance the effect of double refraction of electron due to Rashba spin-orbit interaction. Besides the conventional way by increasing the gate voltage [4], we can also control this effect by changing the Fermi energy of the incident electron.

The exact results for the case of normal incidence are derived analytically. It is found that the influence of the interface at x = 0 does not alter the spin direction of reflected or refracted electron for any initial spin state. It is predicted that the spin-orbit interaction can induce a current perpendicular to the normal incident direction of the electron.

The numerical calculations in the case of oblique incidence reveal more details about the law of electronic double reflection. We find that for an electron with nonzero incident angle, the spin state of the reflected wave often undergoes an abrupt change at the interface x = 0, see Fig. 2. As shown in Fig. 3, the coefficients of reflection and refraction depend on the incident angle and the initial spin state of the electron in a rather complex way. The origin of the two-step character of reflection and one refraction mode, as shown in Fig. 4, is simply attributed to the total reflection of this refraction mode.

The problem of probability current is re-studied in Section 2 and the appendix A. At first, we find that the identity R+T++T=1 requires a subtle understanding about the relevant coefficients T+, T and R; secondly, we prove exactly that the cross term of the current remains zero even taking into account the spin-orbit coupling interaction.

Finally, we expect that the electronic double refraction effect can also be used to design a new kind of spin-filter device similar to the most commonly used Nicol prism for the production of linearly polarized light [36].

References

[1]

Datta S. and Das B., Appl. Phys. Lett., 1990, 56: 665

[2]

Rashba E. I., Fiz Tverd. Tela (Leningrad) , 1960, 2: 1224 [Sov. Phys. Solid State , 1960, 2: 1109]

[3]

Bychkov Y. L. and Rashba E. I., J. Phys. C, 1984, 17: 6039

[4]

Nitta J., Akazaki T., Takayanagi H., and Enoki T., Phys. Rev. Lett., 1997, 78: 1335

[5]

Schäpers Th., Engels G., Lange J., Klocke Th., Hollfelder M., and H. Lüth, J. Appl. Phys., 1998, 83: 4324

[6]

de E. A. Andrada e Silva and La Rocca G. C., Phys. Rev. B, 1999, 59: 15583(R)

[7]

Grundler D., Phys. Rev. Lett., 2000, 84: 6074

[8]

Grundler, D., ibid, 2001, 86: 1058

[9]

Zülicke U. and Schroll C., Phys. Rev. Lett., 2002, 88: 29701

[10]

Feve G., Oliver W. D., Aranzana M., and Yamamoto Y., Phys. Rev. B, 2002, 66: 155328

[11]

Benjamin R. and Benjamin C., Phys. Rev. B, 2004, 69: 085318

[12]

Papp G. and Peeters F. M., Appl. Phys. Lett., 2001, 78: 2184

[13]

Radantsev V. F., Ivankiv I. M., and Yafyasov A. M., Semiconductors , 2003, 37: 200

[14]

Yang W. and Chang K., Phys. Rev. B, 2006, 73: 113303

[15]

Kiselev A. A. and Kim K. W., Appl. Phys. Lett., 2001, 78: 775

[16]

Koga T., Nitta J., Takayanagi H., and Datta S., Phys. Rev. Lett., 2002, 88: 126601

[17]

Governale M., Boese D., U. Zülicke, and Schroll C., Phys. Rev. B, 2002, 65: 140403(R)

[18]

Ionicioiu R. and D’Amico I., Phys. Rev. B, 2003, 67: 041307(R)

[19]

David Z. -Y. Ting and X. Cartoixá, Phys. Rev. B, 2003, 68: 235320

[20]

Cartoixá X., David Z. -Y. Ting., and Chang Y. -C., Appl. Phys. Lett., 2003, 83: 1462

[21]

Schliemann J., Egues J. C., and Loss D., Phys. Rev. Lett., 2003, 90: 146801

[22]

Egues J. C., Burkard G., and Loss D., Appl. Phys. Lett., 2003, 82: 2658

[23]

Hall K. C., Lau W. H., Gündogdu K., Flatté M. E., and Boggess T. F., Appl. Phys. Lett., 2003, 83: 2937

[24]

Shen S. Q., Li Z. J., and Ma Z., Appl. Phys. Lett., 2004, 84: 996

[25]

Ma Z., Li P., and Shen S. Q., Phys. Rev. B, 2004, 70: 125318

[26]

Jiang Y. and Jalil M. B. A., J. Phys. C, 2003, 15: L31

[27]

Schmidt G., Ferrand D., Molenkamp L. W., Filip A. T., and van Wees B. J., Phys. Rev. B, 2000, 62: R4790

[28]

Bandyopadhyay S. and Cahay M., Appl. Phys. Lett., 2004, 85: 1433

[29]

Bandyopadhyay S. and Cahay M., Appl. Phys. Lett., 2004, 85: 1814

[30]

Khodas M., Shekhter A., and Finkel’stein A. M., Phys. Rev. Lett., 2004, 92: 086602

[31]

Shekhter A., Khodas M., and Finkel’stein A. M., Phys. Rev. B, 2005, 71: 125114

[32]

Marigliano Ramaglia V., Bercioux D., Cataudella V., De Filippis G., Perroni C. A., and Ventriglia F., Eur. Phys. J. B, 2003, 36: 365

[33]

Marigliano Ramaglia V., Bercioux D., Cataudella V., De Filippis G., and Perroni C. A., J. Phys. C, 2004, 16: 9143

[34]

Bercioux D. and Marigliano Ramaglia V., Superlattices and Micro- Structures , 2005, 37: 337

[35]

Lee M. and Choi M. S., Phys. Rev. B, 2005, 71: 153306

[36]

See, for example, Born M., and Wolf E., Principle of Optics (7th Ed.), Sec.1.5 and Chapt.15, Cambridge University Press, 1999

Rights & permissions

Higher Education Press and Springer-Verlag 2007

PDF (531KB)

1389

Accesses

0

Citation

Detail

Sections
Recommended

/