Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Спин-орбитальные взаимодействия в полупроводниковых квантовых точках и кольцах

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
808 Кб
Скачать

состоящей из четырех блоков:

 

 

 

 

 

 

 

 

f1(m,ri)

 

g1(m,ri)

 

f21(m,ri)

 

g21(m,ri)

,

S˜1 =

f1(m,ri)

 

g1(m,ri)

 

f21(m,ri)

 

g21(m,ri)

 

 

(m + 1,ri) f1(m + 1,ri)

 

 

 

 

 

g1

 

g21(m + 1,ri)

 

f21(m + 1,ri)

 

 

g1(m + 1,ri) f1(m + 1,ri)

g21(m + 1,ri)

 

f21(m + 1,ri)

 

 

 

 

 

 

f22(m,ri)

 

g22(m,ri)

0

0

 

 

 

 

 

 

S˜2 =

f22

(m,ri)

g22(m,ri)

0

0 ,

 

 

 

 

 

 

 

 

(m + 1,ri)

 

 

 

 

 

 

 

 

 

 

g22

f22(m + 1,ri) 0 0

 

 

 

 

 

 

 

g22(m + 1,ri)

f22(m + 1,ri) 0 0

 

 

 

 

 

 

S˜3

 

0

0

f21(m,1)

g21(m,1)

 

,

 

 

 

 

 

 

= 0 0

f21(m,1)

g21(m,1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 g21(m + 1,1) f21(m + 1,1)

 

 

 

 

 

 

 

0 0 g21(m + 1,1) f21(m + 1,1)

.

 

 

 

 

f22(m,1)

 

 

g22(m,1)

 

f2(m,1)

 

g2(m,1)

 

S˜4

=

 

f22(m,1)

 

 

g22(m,1)

 

f2(m,1)

 

g2(m,1)

 

 

 

 

g22(m + 1,1) f22(m + 1,1)

 

 

 

 

 

 

 

g2(m + 1,1) f2(m + 1,1)

 

 

 

 

g22(m + 1,1) f22(m + 1,1)

g2(m + 1,1) f2(m + 1,1)

 

Значения энергии e являются результатом численного решения уравнения

˜

(2.35)

det S(m,e,v,a,ri) = 0,

после чего находится аналитическое решение системы (2.33) для коэффициентов, с помощью которых определяются радиальные волновые функции u(r) и w(r).

Рис. 2.3 демонстрирует функции u(r) (сплошная линия) и w(r) (штриховая линия) для m = 1, a = 10, ri = 0,8, v = 100 и e = 21,4511.

 

2

 

 

 

 

 

1

 

 

 

 

W(R)

 

 

 

 

 

),

0

 

 

 

 

U(R

 

 

 

 

 

 

-1

 

 

 

 

 

-2

 

 

 

 

 

0.0

0.5

1.0

1.5

2.0

 

 

 

R

 

 

Рисунок 2.3. Радиальные волновые функции для m = 1, ri = 0,8

21

Видно, что радиальные волновые функции быстро убывают под барьерами.

Табл. 2.2 показывает зависимость энергии e от отношения радиусов кольца ri и параметра Рашбы a для фиксированного значения глубины потенциальной ямы v = 100 при двух значениях углового момента m = 0, 1.

Таблица 2.2. Уровни энергии в квантовом кольце для v = 100

ri

a

 

 

e

 

 

 

 

 

 

 

 

 

 

 

m = 0

 

 

 

 

 

 

 

 

0,2

0

8,72

12,77 36,62

43,42 79,52

88,94

0,2

2

6,97

12,47 34,95

43,08 77,90

88,50

0,2

10

17,01

14,43

8,23

20,60 53,15

64,08

0,5

0

19,23

21,22 71,61

74,30

 

0,5

2

17,54

20,91 69,98

73,91

 

0,5

10

7,93

2,62 45,67

50,05

 

0,8

0

54,30

55,61

 

 

 

0,8

2

52,66

55,24

 

 

 

0,8

10

26,57

32,94

 

 

 

 

 

 

 

 

 

 

 

 

m = 1

 

 

 

 

 

 

 

 

0,2

0

12,77

21,74 43,42

57,91 88,94

 

0,2

2

9,56

22,71 40,43

58,84 86,19

 

0,2

10

21,43

12,83

9,39

36,34 62,39

 

0,5

0

21,22

27,03 74,30

81,32

 

0,5

2

18,15

28,05

71,39

82,15

 

0,5

10

12,64

2,34 45,78

59,81

 

0,8

0

55,61

59,49

 

 

 

0,8

2

52,69

60,38

 

 

 

0,8

10

21,45

40,12

 

 

 

 

 

 

 

 

 

 

2.2.Взаимодействие Рашбы в присутствии магнитного

поля

Уравнение Шредингера для электрона в двумерной аксиально-

симметричной потенциальной яме имеет вид

 

(

P2

+ Vc(ρ) + VR + VZ )Ψ = E Ψ,

(2.36)

2Me

где оператор P задан соотношением (1.6), а оператор взаимодействия Рашбы VR представлен выражением (1.7). Здесь имеем обычное вы-

22

ражение для взаимодействия Зеемана

 

VZ = 1 Bz,

(2.37)

2

 

где g – эффективный гиромагнитный фактор, µB – магнетон Бора и σz – матрица Паули.

Как было показано в [42], уравнение (2.36) также допускает разделение переменных (2.6) в полярных координатах ρ, φ. Кроме того, к перечню безразмерных величин (2.3) следует добавить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qeρo2

 

 

 

 

 

 

 

gMe

 

 

 

 

 

 

 

 

 

 

 

 

 

b =

 

 

B, ζ =

 

 

.

 

 

 

 

2~

4Me

Имеем следующие радиальные уравнения:

 

 

 

d2u

 

 

 

1 du

 

 

 

 

 

m

 

 

2

 

 

 

 

 

 

+

 

 

 

 

+(e − vc(r)

(

 

+ br)

4ζb)u =

 

 

d r2

r

dr

r

 

 

 

 

 

 

 

 

 

 

 

 

 

dw

 

m + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= a(

 

 

+

 

 

 

w + brw),

 

 

 

 

 

 

 

 

 

 

 

dr

 

r

 

 

 

d2w

 

1 dw

 

 

 

m + 1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

+(e − vc(r) (

 

 

 

 

+ br) + 4ζb)w =

 

d r2

r

dr

 

 

r

()

=a −dudr + mr u + bru ,

(2.38)

(2.39)

(2.40)

где обезразмеренный удерживающий потенциал vc(r) имеет вид (2.9) и (2.10).

Ищем радиальные волновые функции u(r) и w(r), регулярные в начале координат r → 0 и убывающие на бесконечности r → ∞.

Следуя [42], используем подстановки

 

 

 

 

u(r) = (br2)|m|/2 exp(

 

br2

)f(r),

 

 

2

(2.41)

w(r) = (br2)|m+1|/2 exp(

br2

)g(r),

 

 

 

2

 

 

которые приводят к вырожденным гипергеометрическим уравнениям для прямоугольных потенциалов (2.9) и (2.10).

23

Вслучае квантовой точки рассмотрим две области 0 < r < 1 (область 1) и 1 < r < ∞ (область 2) раздельно.

Вобласти 1, используя известные свойства

M(α,β,ξ)

 

dM(α,β,ξ)

=

β − α

M(α,β + 1),

 

 

 

 

 

β

 

 

 

dM(α,β,ξ)

(β − 1 − ξ)M(α,β,ξ) + ξ

 

 

 

= (β − 1)M(α − 1,β − 1)

 

 

вырожденных гипергеометрических функций M(α,β,ξ) [88], можно показать, что подходящие частные решения системы радиальных уравнений (2.39) и (2.40) строятся так:

 

 

 

 

f1(r) = c+f+(r) + cf

(r),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

1

1

 

 

 

 

 

 

 

 

 

 

g1(r) = (

a

 

 

 

)

(c1+g1+(r) + c1g1(r)),

 

 

(2.42)

 

 

2

 

 

 

 

 

 

b

 

 

 

где f1±, g1± определены равенствами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1±(r) = M(m + 1 − ki±,m + 1,br2),

2

(2.43)

 

 

 

(

ki±

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M m + 1 − ki±,m + 2,br

 

 

 

 

g1±(r) =

 

 

 

m + 1

 

 

 

 

 

( ki± + e/4b + ζ 1/2

 

)

 

для m = 0,1,2, . . . , тогда как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

±

 

 

 

 

 

 

 

 

 

±

 

 

 

 

 

,br2

,

 

 

 

 

 

f1

 

(r) = M(1 − ki ,

− m + 1

 

2

 

)

 

 

(2.44)

 

 

 

 

 

 

 

 

i

(

 

 

 

 

 

 

)

 

 

 

 

 

 

 

g±(r) =

 

 

m M −ki±, − m,br

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

k± + e/4b + ζ 1/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для m = 1, − 2, − 3, . . . Здесь используется обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

a2

 

 

 

 

 

 

 

 

a2

 

4b

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ki± =

4b e + 2 ± ae +

 

4 + ( a ) (ζ

1/2)2

.

(2.45)

Соответствующие волновые функции (2.41) являются регулярными в начале координат.

24

Аналогично в области 2 при применении известных свойств

U(α,β,ξ) dU(α,β,ξ) = U(α,β + 1), dξ

(β − 1 − ξ)U(α,β,ξ) + ξ dU(α,β,ξ) = −U(α − 1,β − 1)

вырожденных гипергеометрических функций U(α,β,ξ) [88] находятся решения

 

 

 

 

f2(r) = c+f+(r) + cf(r),

 

 

 

где

 

 

 

 

 

 

2

2

2

2

 

 

)

 

 

 

 

 

 

 

 

a

(

 

 

 

 

 

 

 

 

 

 

 

 

g2(r) = (

2

 

) c2+g2+(r) + c2g2(r) ,

 

 

 

b

 

 

 

f2±(r) = U(m + 1 − ko±,m + 1,br

2

 

 

 

 

 

2),

 

 

 

 

 

g±(r) =

 

U(m + 1 − ko±,m + 2,br )

 

 

 

 

 

2

 

 

−ko± + (e − v)/4b + ζ − 1/2

для m = 0,1,2, . . . , и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

±

 

 

 

 

 

 

±

 

 

 

 

,br2

,

 

 

 

 

 

f2 (r) = U(1

− ko

, − m + 12

)

)

 

 

 

 

 

g±(r) =

 

 

 

 

U(−ko±, − m,br

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

−ko± + (e − v)/4b + ζ − 1/2

 

 

 

 

 

 

для m = 1, − 2, − 3, . . . Здесь используется обозначение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ko± =

1

 

a2

± ae − v +

 

a2

4b 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4b e − v + 2

 

4 + ( a ) (ζ − 1/2)2 .

(2.46)

(2.47)

(2.48)

(2.49)

Соответствующие волновые функции u(r) и w(r) имеют правильное поведение на бесконечности.

Условия непрерывности для радиальных волновых функций u(r), w(r) и их производных u(r), w(r) в граничной точке r = 1 могут быть записаны следующим образом:

f1(1) = f2(1), f1(1) = f2(1), g1(1) = g2(1), g1(1) = g2(1),

где штрихами обозначены производные по r. Отсюда получаем систему алгебраических уравнений

M4(m,e,v,a,b) X = 0

(2.50)

25

для четырех коэффициентов X = {c+1 ,c1 ,c+2 ,c2 }. Здесь M4(m,e,v,a,b) – 4 × 4 матрица

 

 

 

f1+(1)

f1(1)

f2+(1)

f2(1)

 

 

 

 

 

+

f1−′(1)

+

 

 

(2.51)

M4

=

f1+(1)

−f2+(1)

−f2(1) .

 

 

 

1

1

2

2

 

 

 

 

 

g1 (1)

g1(1)

g2 (1)

g2(1)

 

 

 

 

g+(1)

g−′(1)

g+(1)

g(1)

 

 

Следовательно, точное уравнение для определения энергий e(m,v,a,b) имеет вид

det M4(m,e,v,a,b) = 0.

(2.52)

Это уравнение не может быть решено аналитически, но легко исследуется численно.

После того как значения e(m,v,a,b) найдены из (2.52), получаем значения коэффициентов c+1 ,c1 ,c+2 ,c2 из системы (2.50):

 

2

 

1 3

 

+

 

 

 

 

 

c

 

 

f

+(1)

 

 

 

 

1

 

= c+M1(m,e,v,a,b)

1

 

 

 

(2.53)

 

c+

g1+(1)

 

,

c2

 

g1

(1)

 

 

где

 

 

 

 

+

 

 

 

 

 

M3

=

f1

(1)

−f2+ (1)

−f2

(1)

.

(2.54)

g(1)

g2

(1)

g2(1)

 

 

1

 

+

 

 

 

 

 

 

 

g1(1)

−g2

(1)

−g2−′(1)+

 

 

Остающийся произвол в выборе коэффициента c1 устраняется стан-

дартной нормировкой (2.23). Таким образом, радиальные волновые функции u(r) и w(r) построены полностью.

В дополнение к аналитическим результатам представим результаты численного решения уравнения (2.52). На рис. 2.4–2.7 сплошные линии соответствуют первым уровням энергии, а штриховые — вторым уровням при заданном значении m. В соответствии с [42] выбирается относящееся к арсениду галлия GaAs значение параметра:

g = 0,44, (2.55)

тогда получаем ζ = 0,00737. При выборе ρo = 30 нм имеет место следующее соответствие между безразмерной и размерной величиной:

b = 1 → B = 1,45649 Тл.

(2.56)

26

На рис. 2.4 и 2.5 представлена зависимость уровней энергии от

величины магнитной индукции при фиксированном значении пара-

метра Рашбы.

 

 

 

 

30

 

 

 

 

 

 

 

m 1

 

 

 

25

 

 

 

 

 

20

 

m ‡ -2

 

 

 

 

 

m 0

 

 

 

 

 

 

 

e

 

 

 

 

 

15

 

m 1

 

 

 

 

 

 

 

 

10

m

-2

m ‡ -1

 

 

 

 

 

 

 

 

 

 

 

 

 

m ‡ -1

 

 

5

 

 

m 0

 

 

 

0

1

2

3

4

 

 

 

b

 

 

 

Рисунок 2.4. Зависимость e от b при a = 1 для v = 400

 

30

m 1

 

 

 

 

 

 

 

 

 

25

 

 

 

 

 

 

 

m ‡ -2

 

 

 

 

 

 

m 0

 

 

20

 

 

 

 

 

e

 

 

 

 

 

15

 

m 1

 

 

 

 

 

 

 

 

 

m

-2

m ‡ -1

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

m ‡ -1

 

 

5

 

 

m 0

 

 

 

 

 

 

 

 

0

1

2

3

4

 

 

 

b

 

 

 

Рисунок 2.5. Зависимость e от b при a = 1 для v =

 

27

На рис. 2.6 и 2.7 показана зависимость уровней энергии от пара-

метра Рашбы при фиксированном значении магнитной индукции.

40

m 1

 

 

 

 

 

 

 

30

 

 

 

 

 

 

m 0

 

 

e

m 1

 

 

 

20

 

 

 

 

 

m ‡ -2

 

 

 

10

 

 

m ‡ -1

 

 

 

 

 

 

 

m ‡ -1

m ‡ -2

 

 

 

m 0

 

 

0

 

 

 

 

0

1

2

3

4

 

 

a

 

 

 

Рисунок 2.6. Зависимость e от a при b = 3 для v = 400

 

 

m 1

 

 

 

40

 

 

 

 

30

 

 

 

 

 

 

m 0

 

 

e

m 1

 

 

 

 

 

 

 

20

m ‡ -2

 

 

 

 

 

 

 

10

 

 

m ‡ -1

 

 

 

m ‡ -1

m ‡ -2

 

 

 

 

 

 

 

m 0

 

 

0

 

 

 

 

0

1

2

3

4

 

 

a

 

 

 

Рисунок 2.7. Зависимость e от a при b = 3 для v =

 

Заметим, что характер указанных зависимостей сохраняется при варьировании фиксированных параметров в широком диапазоне. Вид-

28

но, что для всех состояний кривые, соответствующие v = (потенциальная яма с непроницаемой стенкой), лежат выше кривых, соответствующих v = 400.

Рис. 2.8 и 2.9 дают примеры нормированных условием (2.23) радиальных волновых функций u(r) и w(r) при фиксированном значении a = 1, b = 1 для v = 400, m = 0 и m = 1.

3

 

 

 

 

 

 

 

m 0

 

 

 

 

 

 

2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

uHrL

 

 

 

m 0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

m 1

 

 

 

-1

 

 

 

 

 

 

 

-2

 

 

m 1

 

 

 

 

0.0

0.2

0.4

0.6

0.8

1.0

1.2

1.4

 

 

 

 

r

 

 

 

 

Рисунок 2.8. Радиальные волновые функции u(r)

 

2

 

 

m 0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

m 1

 

 

 

 

0

 

 

 

 

 

 

 

wHrL

 

 

 

 

 

 

 

 

 

 

m 0

 

 

 

 

-1

 

 

 

 

 

 

 

-2

 

 

m 1

 

 

 

 

 

 

 

 

 

 

 

0.0

0.2

0.4

0.6

0.8

1.0

1.2

1.4

 

 

 

 

r

 

 

 

 

Рисунок 2.9. Радиальные волновые функции w(r)

 

29

Соответствие сплошных и штриховых линий уровням энергии то же, что и прежде. Следует отметить, что радиальные волновые функции быстро убывают вне потенциальной ямы.

В случае квантового кольца рассмотрим три области 0 < r < ri (область 1), ri < r < 1 (область 2) и 1 < r < ∞ (область 3) раздельно.

Легко показать, что решения радиальных уравнений с правильным поведением около начала координат в первой области строятся так:

 

f1(r) = c+f+(r) + cf(r),

 

 

 

 

1

 

1

 

1

1

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

g1(r) = (

2

 

) (c1+g1+(r) + c1g1(r)),

 

 

где

b

 

 

 

f1±(r) = M(m + 1 − ko±,m + 1,br2),

2

)

 

 

ko±

M(m + 1

 

ko±,m + 2,br

 

g1±(r) = (

 

)

 

 

ko± + (e

 

1/2

 

m + 1

 

 

 

v)/4b + ζ

 

 

 

 

 

 

 

 

 

 

(2.57)

(2.58)

для m = 0,1,2, . . . , и

 

 

 

 

 

 

 

±

 

 

±

,br2

)

,

f1

(r) = M

(1 − ko , − m + 1

2

)

(2.59)

g±(r) =

 

m M(−ko±, − m,br

 

 

 

 

 

 

1

 

−ko± + (e − v)/4b + ζ − 1/2

 

 

для m = 1, − 2, − 3, . . . Здесь величины ko± имеют вид (2.49).

В области 3 решения радиальных уравнений с правильным поведением на бесконечности строятся согласно соотношениям

f3(r) = c+f+(r) + cf(r),

 

3

2

3

2

 

a

 

 

(2.60)

g3(r) = (

2

 

)

(c3+g2+(r) + c3g2(r)),

b

где функции f2±(r) и g2±(r) имеют вид (2.47) при положительных

и(2.48) при отрицательных значениях m.

Вобласти 2 нет ограничений на выбор частных решений, и искомые функции строятся так:

f2(r) = c+ f+

(r) + cf

(r) + c+ f+

(r) + cf

(r),

 

 

21

21

 

21

21

 

22

22

 

22

22

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

(2.61)

g2(r) = (

2

 

) (c21+ g21+

(r) + c21g21

(r) + c22+ g22+

(r) + c22g22

(r)),

b

30

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]