Спин-орбитальные взаимодействия в полупроводниковых квантовых точках и кольцах
.pdfРезультаты численного решения уравнения (3.26) представлены в табл. 3.2–3.4, которые демонстрируют зависимость уровней энергии e от отношения радиусов кольца ri, значения квантового числа m и глубины потенциальной ямы v при a = 0.
Таблица 3.2. Зависимость e от ri при v = 100
m |
ri |
|
e |
|
|
|
|
|
|
0 |
0,1 |
6,44373 |
29,3430 |
67,2040 |
0 |
0,5 |
19,2307 |
71,6109 |
|
0 |
0,9 |
80,7533 |
|
|
|
|
|
|
|
1 |
0,1 |
12,1534 |
40,4349 |
82,7073 |
1 |
0,5 |
21,2235 |
74,2974 |
|
1 |
0,9 |
82,0011 |
|
|
|
|
|
|
|
2 |
0,1 |
21,6651 |
57,2693 |
|
2 |
0,5 |
27,0296 |
81,3150 |
|
2 |
0,9 |
85,5640 |
|
|
|
|
|
|
|
Таблица 3.3. Зависимость e от ri при v = 400
m |
ri |
|
|
|
e |
|
|
0 |
0,1 |
8,40562 |
36,9454 |
85,0236 |
151,822 |
235,986 |
334,551 |
0 |
0,5 |
26,8322 |
107,288 |
235,075 |
384,438 |
|
|
0 |
0,9 |
218,218 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0,1 |
13,5274 |
45,7776 |
96,8514 |
166,519 |
253,899 |
355,996 |
1 |
0,5 |
28,7345 |
109,360 |
237,247 |
387,653 |
|
|
1 |
0,9 |
219,339 |
|
|
|
|
|
2 |
0,1 |
23,9099 |
64,1730 |
122,141 |
197,573 |
289,541 |
392,547 |
2 |
0,5 |
34,3903 |
115,577 |
243,741 |
394,771 |
|
|
2 |
0,9 |
222,701 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3.4. Зависимость e от ri при v = 700
m |
ri |
|
|
|
|
e |
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0,1 |
8,99511 |
39,2512 |
90,2916 |
161,661 |
252,768 |
362,667 |
489,615 |
629,427 |
0 |
0,5 |
29,2605 |
117,701 |
262,291 |
456,199 |
671,751 |
|
|
|
0 |
0,9 |
296,625 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
0,1 |
13,9367 |
47,3751 |
100,653 |
173,853 |
266,753 |
378,764 |
508,587 |
652,099 |
1 |
0,5 |
31,1499 |
119,735 |
264,374 |
458,342 |
674,608 |
|
|
|
1 |
0,9 |
297,740 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
0,1 |
24,4912 |
65,8332 |
125,629 |
204,061 |
301,176 |
416,680 |
549,380 |
692,161 |
2 |
0,5 |
36,7773 |
125,841 |
270,626 |
464,762 |
681,873 |
|
|
|
2 |
0,9 |
301,084 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
41
Из данных таблиц видно, что число допустимых связанных состояний увеличивается с ростом v и уменьшается с ростом ri.
3.2.Взаимодействия Рашбы и Дрессельхауса равных интенсивностей в присутствии магнитного поля
При рассмотрении состояний электрона с учетом внешнего однородного постоянного магнитного поля в уравнении Шредингера (3.2) следует заменить оператор импульса p на оператор P вида (1.6), а также учесть взаимодействие Зеемана VZ вида (2.37).
Полный гамильтониан рассматриваемой задачи представим в виде суммы невозмущенного гамильтониана
Px2 + Py2 |
+ α(σx − σy)(Px + Py)/~ + Vc(ρ), |
(3.29) |
H0 = 2Me |
где Vc(ρ) – удерживающий аксиально-симметричный прямоугольный потенциал, и возмущения VZ .
Так как не удается найти точное решение полного уравнения
(H0 + VZ )Ψ = E Ψ, |
(3.30) |
сначала получим точное решение невозмущенного уравнения
H0Ψ0 = E0Ψ0, |
(3.31) |
а потом учтем влияние взаимодействия Зеемана в рамках теории возмущений.
Заметим, что и при учете внешнего магнитного поля в случае невозмущенного уравнения Шредингера (3.31) с гамильтонианом (3.29) сохраняются оба интеграла движения (3.3) и (3.4), и искомые решения Ψ±0 (x,y) уравнения (3.31) имеют вид (3.5). Специфика этих решений определяется конкретными выражениями для радиальных волновых функций u(ρ).
Перейдя к безразмерным величинам (2.3) и (2.38), запишем радиальное уравнение
d2u |
|
1 du |
|
m |
2 |
|
|||||
|
+ |
|
|
|
|
+(e0 + a2 |
− vc(r) − ( |
|
|
+ br) )u = 0 |
(3.32) |
d r2 |
r |
dr |
r |
||||||||
для прямоугольных потенциальных ям (2.9) и (2.10), моделирующих квантовую точку и квантовое кольцо. Здесь введено обозначение
e0 |
= |
2Me ρo2 |
E0. |
(3.33) |
|
||||
|
~2 |
|
|
|
42
В каждой из областей изменения r радиальное уравнение (3.32) совпадает с уравнением для радиальной волновой функции электрона в однородном магнитном поле без учета спин-орбитального взаимодействия [89], если заменить e0 на e0 + a2 −vc(r). Поэтому, следуя [89], представим искомую функцию u(r) в виде
u(r) = (br2)|m|/2 exp( |
br2 |
)f(r), |
|
−2 |
(3.34) |
где f(r) выражается через вырожденные гипергеометрические функции. Из уравнения (3.32) видно, что выполняется следующее соотношение:
e0(m,v,a,b) − e0(−m,v,a,b) = 4bm. |
(3.35) |
Выражение для взаимодействия Зеемана в безразмерных величинах принимает вид
vz = |
2Me ρo2 |
VZ = 4ζbσz. |
(3.36) |
|
|||
~2 |
|
|
|
Напоминаем, что и при наличии магнитного поля каждый уровень энергии невозмущенной системы является двукратно вырожденным с собственными функциями (3.5).
Вследствие того, что σz n± = n , в базисе собственных векторов |Ψ+0 и |Ψ−0 невозмущенного гамильтониана для диагональных матричных элементов выполняются равенства Ψ±0 |σz|Ψ±0 = 0. Недиагональные матричные элементы имеют вид
Ψ0+|σz|Ψ0− = Ψ0−|σz|Ψ0+ = ∆(m,v,a,b), |
(3.37) |
|||
где |
|
|
|
|
∆(m,v,a,b) = |
∫0 ∞ J0(2ar)u2(r)r dr |
. |
(3.38) |
|
|
||||
|
∫0 |
∞ u2(r)r dr |
||
Тогда для возмущенных уровней энергии получаем |
|
|||
e± = e0 ± e′, |
(3.39) |
|||
где |
|
|
(3.40) |
|
e′ ≡ e′(m,v,a,b) = 4ζb∆(m,v,a,b). |
||||
43
Для поправки e′ выполняется соотношение
e′(−m,v,a,b) = e′(m,v,a,b). |
(3.41) |
Нормированные собственные функции нулевого приближения, которые соответствуют собственным значениям e±, описываются формулами
|
Ψ± = Ψ0+ ± Ψ0− . |
|
|
|
(3.42) |
||||||
|
|
|
√ |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
||||
Заметим, что в предельном случае α = 0 выражения (3.39) и |
(3.42) |
||||||||||
становятся точными. |
|
|
|
|
|
|
|
|
|
|
|
В случае круговой квантовой точки функция f(r) имеет вид |
|||||||||||
c U(γo,βm,b)M(γi,βm,br2), |
0 < r < 1, |
|
|||||||||
f(r) = {c M(γi,βm,b)U(γo,βm,br2), |
1 < r < ∞, |
(3.43) |
|||||||||
где |
|
|
|
|
e0 + a2 |
|
|
||||
γ = |
m + |m| + 1 |
|
, |
|
|||||||
|
|
(3.44) |
|||||||||
i |
2 |
|
− |
4b |
|
|
|||||
γo = γi + v/4b, βm = |m| + 1. |
|
||||||||||
Множитель c выбирается |
так, |
чтобы |
обеспечить нормировку |
||||||||
Ψ±0 |Ψ±0 = 1.
Условие непрерывности для первой производной волновой функ-
ции в граничной точке r = 1 приводит к уравнению |
|
|
|||||||
γoU(γo + 1,βm + 1,b) |
+ |
|
γi |
|
M(γi + 1,βm + 1,b) |
= 0 |
(3.45) |
||
|
|
|
|
|
|
|
|||
|
U(γo,βm,b) |
βm |
|
M(γi,βm,b) |
|||||
|
|
|
|
|
|||||
для определения энергии e0(m,v,a,b).
Теперь представим результаты вычислений для квантовой точки при ζ = −0,00737 и v = 400.
На рис. 3.3 и 3.4 в безразмерных единицах показана зависимость невозмущенной энергии e0 и относительной поправки −e′/e0 к энергии от магнитной индукции b при фиксированном значении параметра a = 1.
44
140 |
|
|
|
|
120 |
|
|
|
|
100 |
|
|
|
|
80 |
|
|
m ‡ 2 |
|
0 |
|
|
|
|
e |
|
|
|
|
60 |
|
|
m ‡ 1 |
|
|
|
|
|
|
40 |
|
|
m ‡ 2 |
|
|
|
|
m ‡ 0 |
|
20 |
|
|
m ‡ 1 |
|
|
|
|
|
|
|
|
|
m ‡ 0 |
|
0 |
1 |
2 |
3 |
4 |
0 |
||||
|
|
b |
|
|
Рисунок 3.3. Зависимость e0 |
от b при a = 1 для квантовой точки |
|
||
|
0.015 |
|
|
|
|
|
|
|
|
m ‡ 0 |
|
0 |
0.010 |
|
|
|
|
e |
|
|
|
|
|
-e¢ • |
|
|
|
|
|
|
|
|
|
|
|
|
0.005 |
|
|
m ‡ 2m ‡ 1m ‡ 0 |
|
|
|
|
|
m ‡ 1 |
|
|
|
|
|
m ‡ 2 |
|
|
0.000 |
1 |
2 |
3 |
4 |
|
0 |
||||
|
|
|
b |
|
|
Рисунок 3.4. Зависимость −e′/e0 |
от b при a = 1 для квантовой точки |
||||
45
На рис. 3.5 и 3.6 представлена зависимость e0 и |
−e′/e0 от интен- |
|||||||
сивности спин-орбитального взаимодействия a при фиксированном |
||||||||
значении магнитной индукции b = 1. |
|
|
|
|||||
|
140 |
|
|
|
|
|
|
|
|
120 |
|
|
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
80 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
e |
|
|
|
m ‡ 2 |
|
|
|
|
60 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
40 |
|
|
|
m ‡ 1 |
|
|
|
|
|
|
|
|
m ‡ 2 |
|
|
|
|
20 |
|
|
|
m ‡ 0m ‡ 1 |
|
|
|
|
0 |
|
|
|
m ‡ 0 |
|
|
|
|
|
|
0.5 |
1.0 |
|
1.5 |
2.0 |
|
|
0.0 |
|
|
|
||||
|
|
|
|
|
a |
|
|
|
|
Рисунок 3.5. Зависимость e0 от a при b = 1 для квантовой точки |
|||||||
|
0.010 |
|
|
|
|
|
|
|
|
0.008 |
|
|
|
|
|
|
|
0 |
0.006 |
m ‡ 0 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
¢ •e |
0.004 |
|
|
|
|
|
|
|
e |
|
|
|
|
|
m ‡ 0 |
|
|
- |
|
|
|
|
|
m ‡ 1 |
|
|
|
|
|
|
|
|
m ‡ 2 |
||
|
0.002 |
m |
‡ |
1 |
|
|
||
|
|
|
|
|
|
|||
|
|
m ‡ 2 |
|
|
|
|
||
|
0.000 |
|
|
|
|
|
|
|
|
-0.002 |
|
|
|
|
|
|
|
|
0.0 |
|
|
0.5 |
1.0 |
|
1.5 |
2.0 |
|
|
|
|
|
a |
|
|
|
Рисунок 3.6. Зависимость −e′/e0 от a при b = 1 для квантовой точки |
||||||||
46 |
|
|
|
|
|
|
|
|
В случае кругового квантового кольца функция f(r) имеет вид
|
c1f1(r), |
0 < r < ri, |
|
|
f(r) = |
c21f21(r) + c22f22(r), |
ri < r < 1, |
(3.46) |
|
|
|
1 < r < |
∞ |
. |
|
c3f3(r), |
|
||
Здесь c1, c21, c22, c3 – произвольные коэффициенты, а функции f1(r), f21(r), f22(r), f3(r) выражаются через вырожденные гипергеометрические функции следующим образом:
f1(r) = M(γo,βm,br2),
f21(r) = M(γi,βm,br2), f22(r) = U(γi,βm,br2), (3.47)
f3(r) = U(γo,βm,br2).
Условия непрерывности для функции u(r) и ее первой производной в граничных точках r = ri и r = 1 приводят к системе алгебраических уравнений
˜ |
(m,e0,v,a,b,ri) X = 0 |
(3.48) |
T4 |
˜
для четырех коэффициентов X = {c1,c21,c22,c3}; T4(m,e0,v,a,b,ri) – 4×4 матрица вида
T˜4 |
|
f1(ri) |
f21(ri) |
f22(ri) |
0 |
. |
(3.49) |
|
= |
f1′ |
(ri) |
−−f21′ (ri) |
−−f22′ (ri) |
0 |
|||
|
|
|
0 |
f21(1) |
f22(1) |
|
|
|
|
|
|
−f3(1) |
|
||||
|
|
|
0 |
f21′ (1) |
f22′ (1) |
−f3′(1) |
|
|
Следовательно, значения энергии e0(m,v,a,b,ri) определяются уравнением
˜ |
(m,e0,v,a,b,ri) = 0. |
(3.50) |
det T4 |
Приведем графические иллюстрации для кругового квантового
кольца при ζ = −0,00737, ri = 0,5 и v = 400.
На рис. 3.7 и 3.8 представлена зависимость невозмущенной энергии e0 и относительной поправки −e′/e0 к энергии от магнитной индукции b при фиксированном значении a = 1. На рис. 3.9 и 3.10 представлена зависимость e0 и −e′/e0 от интенсивности спин-орбитального взаимодействия a при фиксированном значении b = 1.
47
140 |
|
|
m ‡ 2 |
|
|
|
|
|
|
120 |
|
|
m ‡ 1 |
|
|
|
|
m ‡ 0 |
|
100 |
|
|
|
|
80 |
|
|
|
|
0 |
|
|
|
|
e |
|
|
|
|
60 |
|
|
|
|
|
|
|
m ‡ 2 |
|
40 |
|
|
m ‡ 1 |
|
|
|
|
m ‡ 0 |
|
20 |
|
|
|
|
0 |
1 |
2 |
3 |
4 |
0 |
||||
|
|
b |
|
|
Рисунок 3.7. Зависимость e0 |
от b при a = 1 для квантового кольца |
|||
|
0.0015 |
|
|
m ‡ 0 |
|
|
|
|
|
m ‡ 1 |
|
0 |
0.0010 |
|
|
|
|
•e |
|
|
m ‡ 2 |
|
|
-e¢ |
|
|
|
|
|
|
|
|
m ‡ 1m ‡ 0 |
|
|
|
|
|
|
m ‡ 2 |
|
|
0.0005 |
|
|
|
|
|
0.0000 |
|
|
|
|
|
0 |
1 |
2 |
3 |
4 |
|
|
|
b |
|
|
Рисунок 3.8. Зависимость −e′/e0 |
от b при a = 1 для квантового кольца |
||||
Табл. 3.5 показывает существенную зависимость невозмущенной энергии e0 и относительной поправки −e′/e0 к энергии от отношения радиусов кольца ri при фиксированных значениях a = 1 и b = 1. В скобках приведены значения величины −e′/e0. Значению ri = 0 соответствует квантовая точка. В табл. 3.5 приведены только два первых уровня энергии. Заметим, что при ri = 0,9 уже вторые уровни энергии лежат выше потенциальной ямы, безразмерная глубина v
48
которой в рассматриваемом случае равна 400. Если магнитное поле |
|||||
отсутствует, то связанные состояния имеют место только при e0 < v. |
|||||
Наличие магнитного поля приводит к появлению дискретных уровней |
|||||
энергии, превосходящих по величине глубину потенциальной ямы. |
|||||
|
140 |
|
|
|
|
|
120 |
|
m ‡ 2 |
|
|
|
|
m ‡ 1 |
|
|
|
|
|
|
|
|
|
|
100 |
|
m ‡ 0 |
|
|
|
|
|
|
|
|
|
80 |
|
|
|
|
|
0 |
|
|
|
|
|
e |
|
|
|
|
|
60 |
|
|
|
|
|
40 |
|
m ‡ 2 |
|
|
|
|
|
m ‡ 1 |
|
|
|
20 |
|
m ‡ 0 |
|
|
|
0 |
0.5 |
1.0 |
1.5 |
2.0 |
|
0.0 |
||||
|
|
|
a |
|
|
Рисунок 3.9. Зависимость e0 от a при b = 1 для квантового кольца |
|||||
|
0.0010 |
m ‡ 0 |
|
|
|
|
m ‡ 1 |
|
|
|
|
|
|
|
|
|
|
|
0.0008 |
m ‡ 2 |
|
|
|
|
0.0006 |
m ‡ 0 |
|
|
|
0 |
|
|
|
|
|
e |
|
|
|
|
|
|
m ‡ 1 |
|
|
|
|
¢• |
0.0004 |
|
|
|
|
m ‡ 2 |
|
|
|
||
-e |
|
|
|
|
|
|
0.0002 |
|
|
|
|
|
0.0000 |
|
|
|
|
|
-0.0002 |
|
|
|
|
|
0.0 |
0.5 |
1.0 |
1.5 |
2.0 |
|
|
|
a |
|
|
Рисунок 3.10. Зависимость −e′/e0 |
от a при b = 1 для квантового кольца |
||||
|
|
|
|
|
49 |
Таблица 3.5. Зависимость e0 и −e′/e0 от ri при a = 1 и b = 1
ri |
m = 0 |
m = 1 |
||
|
|
|
|
|
|
|
e0, |
(−e′/e0) |
|
0 |
4,48334 |
26,952 |
14,677 |
45,9308 |
|
(5,162 × 10−3) |
(7,679 × 10−4) |
(1,357 × 10−3) |
(4,378 × 10−4) |
0,1 |
7,72294 |
36,3199 |
14,9000 |
47,1566 |
|
(2,737 × 10−3) |
(5,504 × 10−4) |
(1,327 × 10−3) |
(4,216 × 10−4) |
0,5 |
26,4059 |
106,878 |
30,3106 |
110,949 |
|
(5,687 × 10−4) |
(1,391 × 10−4) |
(4,937 × 10−4) |
(1,340 × 10−4) |
0,9 |
218,124 |
404,788 |
221,245 |
407,21 |
|
(3,812 × 10−5) |
(−1,216 × 10−5) |
(3,757 × 10−5) |
(−2,038 × 10−5) |
В качестве развития изложенного подхода можно указать на возможность учета поправок более высокого порядка в рамках теории возмущений.
