при энергиях, по крайней мере, на порядок превышающих те, что актуальны для данной работы.
§ 3.3. Вероятности переходов оже-типа 3v → 11
Рассмотрим оже-переходы следующего типа: электрон с двумерным волновым вектором k3 , находящийся в подзоне 3 квантовой ямы, взаимодействует с электроном с волновым вектором kv = (kv , kv||).в валентной зоне области B В результате оба электрона переходят в подзону 1 квантовой ямы в состояния с волновыми векторами k11 и k12 . Переданный при этом импульс в плоскости, перпендикулярной оси роста наноструктуры, равен
q0 = k3 −k11 = k12 −kv . |
(3.14) |
Как и при выводе формул (2.16, 2.17) для вероятности переходов, 31 → 22, представим кулоновское взаимодействие в виде разложения в трехмерный ряд Фурье и выполним ряд стандартных преобразований. В результате получим следующее выражение для прямого матричного элемента, построенного на функциях (3.1-3.4):
(d) |
|
4πe2 |
|
∫dq|| ∫∫dz1dz2 |
eiq|| ( z1−z2 ) |
|
||||||
M k 3 , kv , k3 −q0 ,kv +q0 = |
|
|
|
|
|
|
× |
|
||||
aε |
L |
|
q2 |
+ q2 |
|
|||||||
|
|
|
|
|
|
0 |
|| |
|
|
|||
×β* (z )β (z )β* (z |
2 |
)β (z |
2 |
)Θ(cc)Θ(vc) |
, |
|
(3.15) |
|||||
3 1 |
1 1 |
v |
|
1 |
−q q |
|
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
(mn) |
|
* |
|
|
|
|
|
|
|
|
|
(3.16) |
Θq0 |
= ∫dr um K (r) un K+q0 (r) |
|
|
|
||||||||
– интегралы перекрывания блоховских амплитуд для m-й и n-й зон, которые для простоты считаются зависящими только от q0, но не от K. В свою очередь,
Θq(cc) 1, Θq(vc) χ q0 , χ ~ |
h |
πcv , |
(3.17) |
|||
|
||||||
|
|
|
m Eg |
|
||
где |
|
|
|
|
|
|
π = p + |
h |
[σ, VL ] , |
(3.18) |
|||
4mc2 |
||||||
|
|
|
|
|
||
p – оператор импульса, VL – периодический потенциал решетки, σ – матрица Паули. Выполняя в (3.15) интегрирование по dq||, получим
M (d ) |
|
(k |
v|| |
) |
= |
4π 2e2 |
χ |
|| |
(iI% |
+ 2 cosθ |
I% |
) |
||
k |
3 |
, kv , k3 −q0 |
, kv +q0 |
|
|
εLa |
2 |
|
1 |
|
q0π 2 |
|
||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
≡ 4π 2e2 χ M% (d ) |
|
|
|
|
(kv|| ), |
|
(3.19) |
||||
|
|
|
SεL |
L |
|
k 3 , kv , k3 −q0 , kv +q0 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
38
где
I%1 |
= ∫∫dz1dz2 e |
0 1 |
2 |
| |
sign(z1 − z2 )β3 |
|
(z1)β1(z1)βvkv|| (z2 )β1(z2 ) , |
(3.20) |
|||||
|
|
−q | z |
−z |
|
|
|
|
* |
|
(+)* |
|
||
|
I%2 = ∫∫dz1dz2 e |
|
0 1 |
2 |
| |
β3 |
(z1)β1 |
(z1)βvkv|| |
(z2 )β1(z2 ) . |
(3.21) |
|||
|
|
|
|
−q | z |
−z |
* |
|
|
(−)* |
|
|
||
Выражение для вероятности перехода запишем в виде:
где
f (k3 )
|
W |
|
= |
|
|
e4 |
|
|
d |
3 |
k |
|
|
f (k |
|
)d |
2 |
k |
|
|
|
|
d |
2 |
q |
× |
|
|
|
|
|
(2π)2 hεL2a7 ∫ |
|
v |
∫ |
3 |
|
3 |
∫ |
|
|
|
|||||||||||||||||
|
3v,11 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
(d ) |
|
|
|
|
|
|
|
|
(e) |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
||
×| Mk3 ,kv ,k3 −q0 ,kv +q0 (kv|| ) + Mk3 ,kv ,k3 −q0 ,kv +q0 (kv|| ) | |
|
|
|
|
|||||||||||||||||||||||||
|
va2 − |
h2kv2|| |
|
h2 |
(k3 −kv ) q0 − |
h2q2 |
|
|
|
h2k2 |
|
|
|
(3.22) |
|||||||||||||||
×δ |
|
|
|
+ |
|
|
|
0 |
− |
|
|
|
v |
|
, |
||||||||||||||
|
2mv |
mc |
mc |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2mr |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
= (E03 − 2E01 − E%g ), |
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.23) |
|||||||||
− функция распределения электронов в подзоне 3. Верхними ин-
дексами d и e отмечены прямой и обменный матричные элементы перехода. Вычисление семикратного интеграла в правой части (3.22) существенно упрощается, если учесть, что релаксация электронов внутри подзоны 3 происходит значительно быстрее других релаксационных процессов. В этом случае можно, не теряя точности вычислений, предположить, что все электроны в подзоне 3 находятся вблизи ее дна. Тогда в подынтегральной функции в правой части (3.22) можно положить k3 = 0. При этом
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ f (k3 )d 2k3 = (2π)2 n3a2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.24) |
||||||||||||||||||
|
|
Интегрирование по углам между k3 и q0 снимается с помощью |
|||||||||||||||||||||||||||||||||||||||||||||
δ -функции. После некоторых вычислений получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
1 |
|
|
W |
|
|
e4m |
χ2 kv( m) |
|
kv(||m) |
|
|
ϑ+ dq02 |
|
Mq(d ) (kv|| ) + Mq(e) (kv|| ) |
|
2 |
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
= |
|
3v,11 |
|
= |
c |
|| |
|
∫ |
dkv |
∫ |
dkv|| ∫ |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
, |
(3.25) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
τ |
aug |
|
S n |
|
2h3a5ε2 |
|
|
|
|
|
|
(q |
2 |
−ϑ ) (ϑ |
−q |
2 |
) |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
3 |
|
|
|
|
L |
0 |
|
|
0 |
|
|
ϑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
0 |
|
− |
|
|
+ |
0 |
|
|
|
|
|
|
|
|
||||||
|
|
|
ϑ |
= |
1 |
(k |
|
± ρ), ρ = |
k |
2 |
+ |
4m |
|
|
a |
2 |
− |
|
|
h2kv||2 |
− |
h2k2 |
|
|
|
, |
|
|
|
|
|
||||||||||||||||
|
|
|
v |
v |
|
c |
v |
|
|
|
|
|
|
|
|
v |
|
|
|
|
|
|
|
(3.26) |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2m |
2m |
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
± |
2 |
|
|
|
|
|
|
|
|
|
h2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
|
r |
|
|
|
|
|
|
|
|
|||||
|
(m) |
|
|
|
m |
4m a2 |
|
2 |
|
|
m |
|
|
|
|
|
|
(m) |
|
a 4m m |
|
|
|
|
|
|
|
||||||||||||||||||||
kv|| |
= |
|
|
v |
|
|
c2 |
v |
+ kv 1 |
+ |
2 |
c |
|
, |
|
kv |
|
|
|
= |
|
|
|
|
v |
|
|
c |
|
|
|
v . |
(3.27) |
||||||||||||||
|
2m |
|
|
|
|
h |
|
2m + m |
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
h |
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
v |
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
v |
|
|
|
|
|||||||||
Характерная зависимость τaug-1 ( v ) изображена на рис. 3.4.
Обсудим роль разогрева носителей за счет внутриподзонного поглощения света. В случае обычной межзонной оже-рекомбинации учет разо-
39
грева (точнее, квазиравновесной нефермиевской части функции распределения, см., например, [68]) принципиально важен, т.к. лишь носители с кинетической энергией, превышающей определенное пороговое значение могут участвовать в оже-процессе. Для рассматриваемых в настоящей работе процессов оже-типа это не так – в них могут участвовать электроны в подзоне 3, даже если их кинетическая энергия равна нулю. Роль разогрева
Рис. 3.4. Зависимость скорости оже-переходов 3v → 11 от
v = E03 − 2E01 − Ev
здесь в некотором смысле играет фотовозбуждение носителей на переходах между подзонами 2 и 3. Этот процесс учитывается в нашей модели явным образом. Нужно лишь, чтобы величина v была положительной, что и предполагается в работе. В рассматриваемом нами случае глубоких кван-
товых ям и достаточно больших hω из-за малости коэффициента внутриподзонного поглощения света (на 2-3 порядка меньше, чем для межподзонного поглощения) разогрев электронов светом оказывается слабым. Так в наиболее актуальной области интенсивностей j ~ jth изменение электронной температуры, как показывают грубые оценки, едва ли может превысить 10 К. Это, конечно, не может привести к каким-либо ощутимым последствиям для рассматриваемых в работе эффектов. В области же высоких интенсивностей разогрев (в широком смысле) электронов должен приниматься во внимание наряду с другими факторами, осложняющими анализ в этом случае (см. заключительную часть § 3.4). Однако действительно существенную роль внутриподзонное фотовозбуждение носителей могло бы, в принципе, играть лишь при малых значениях v d 30÷50 мэВ. Заметим, что возможное увеличение вероятностей процессов оже-типа могло бы лишь увеличить эффективность рассматриваемого механизма каскадно-
40
лавинной генерации пар, понижая пороговое значение jth интенсивности света.
§ 3.4. Уравнения баланса для заселенностей
При составлении уравнений баланса следует иметь в виду, что при больших интенсивностях накачки j концентрации электронов в подзонах 2 и 3 становятся сопоставимыми, и следует учитывать фотоиндуцированные переходы как с поглощением, так и с испусканием фотона. Сказанное не относится к переходам v → 1, 1 → 2 и 3 → с, так как конечные состояния для этих переходов практически не заполнены из-за отдаленности от минимумов соответствующих подзон в квантовой яме или в зоне проводимости. Обладая относительно большой кинетической энергией, электроны или дырки в этих конечных состояниях быстро теряют ее за счет внутризонной или внутриподзонной релаксации и покидают область резонанса.
Система уравнений для концентраций неравновесных дырок p в ва-
лентной зоне, электронов n1,2,3 в трех подзонах зоны проводимости и электронов в непрерывном спектре зоны проводимости nc имеет вид:
|
|
nc =σ3c jn3 −Wc3nc − dcnc ( p0 + p) , |
|
|
|
|
|
|
& |
|
|
|
|
|
& |
-1 |
|
|
|
|
|
n3 |
=Wc3nc −σ3c jn3 − (W31 +W32 +τaug )n3 |
|
|
|
|
|
|
+σ23 j(n2 − n3 ) −W31,22 (n1,n3 ) , |
|
|
|
|
n2 = −W21n2 −σ23 j(n2 − n3 ) +W32n3 + 2W31,22 (n1,n3 ) +σ12 jn1 , |
|
|||||
& |
|
|
|
|
|
(3.28) |
|
|
|
|
|
|
|
& |
= −σ12 jn1 |
−1 |
(q) |
j |
q |
|
n1 |
+W21n2 +W31n3 −W31,22 (n1,n3 ) + 2τaug n3 +σv1 |
|
|
|||
|
|
− d1n1( p0 + p) , |
|
|
|
|
|
|
p = n1 + n2 + n3 + nc , |
|
|
|
|
с начальными условиями: n2 , n3, nc , p = 0, n1 = n10 |
при t = 0. В формулах |
|||||
(3.28), Wc3 — скорость захвата электронов из непрерывного спектра зоны проводимости в подзону 3 в квантовой яме. Индекс q принимает значения 1 или 2 для случаев одно- и двухфотонного перехода между валентной зоной в области B и подзоной 1 в квантовой яме. В уравнениях (3.28) фигурируют «двумерные» концентрации p, p0, nc частиц непрерывного зонного спектра, которые отличаются от обычных концентраций p%, p%0 , n%c множи-
телем nw [см. § |
3.2, текст после формулы (3.11)]: nc = ncnw и т.п. Члены |
|
% |
−dcnc ( p0 + p) и |
−d1n1( p0 + p) описывают обычную бимолекулярную ре- |
комбинацию электронов в зоне проводимости и в подзоне 1 с дырками в валентной зоны области B. Величина p0 обозначает равновесную концентрацию дырок. Равновесную концентрацию электронов в зоне проводимости для простоты считаем равной нулю.
41
Состояния двух электронов в подзоне 2, куда они попадают в результате оже-перехода 31 → 22, в силу законов сохранения энергии и импульса характеризуются при не очень малых значениях расстройки резонанса ω −ω21 относительно большими значениями волновых векторов k21 и k22 .
Тем не менее, в рассматриваемой здесь модели, как показывают численные расчеты, концентрация n2 может при актуальных интенсивностях света j оказаться достаточной для заполнения этих состояний. При этом имеет место насыщение вероятностей оже-переходов, которое мы приближенно
описываем множителем [1+ (n2
nf ) r ] −1 . Приведенные ниже результаты
численных расчетов получены при использовании значений nf ≈ 3·1012 см-2, r = 4. Варьирование в разумных пределах значений nf и r, равно как и вероятности Wc3 захвата свободного электрона квантовой ямой, не приводит к качественному изменению результатов. Это обстоятельство, а также исследование зависимости решений (3.28) от других параметров (см. обсуждение в § 3.5), позволяет рассчитывать на то, что качественный характер результатов не изменится за счет многочастичных эффектов (см. обсуждение в § 2.4). В принципе, изменение заселенностей электронов в подзонах 1, 2, 3 может оказывать некоторое влияние на величины σi,i+1. Тем не менее, этот эффект едва ли может быть существенным в рассматриваемом нами случае глубоких квантовых ям, когда при всех актуальных значениях концентраций ni плазменные частоты электронов в яме ωp малы по сравнению с зазорами ωij между подзонами размерного квантования.
На рис. 3.5 показаны зависимости заселенностей n1,2,3 и nc от времени, прошедшего с момента включения накачки. Видно, что изменение за-
селенностей всех актуальных состояний происходит в некоторый момент времени (для выбранных значений параметров это 3.1 нс) скачкообразно. На рис. 3.6 изображены зависимости времени установления квазиравновесного распределения электронов τeq от интенсивности накачки j для одно- и двухфотонного варианта эффекта фотонной лавины.
На рис. 3.7 (а, б) приведены зависимости квазиравновесных заселенностей n1,2,3 и nc от интенсивности накачки j. В расчете были использованы следующие значения параметров, фигурирующих в правых частях уравне-
ний (3.28): a = 3·107 см, |
W31 = 0.02 пс−1, Wс3 = 0.01 пс−1, W32 = 0.07 пс−1, |
|
W21 = 0.1 пс−1, τaug = 2 пс, |
n10 = 0, p0 = 3·1010 см−2, |
σ12 = 0.003 см2/(пс·МВт), |
σ23 = 2 см2/(пс·МВт), σ3c = 0.025 см2/(пс·МВт), |
σv(1)1 =1.5 1010 МВт-1пс-1 , |
|
σv(2)1 =104 см2 МВт-2 пс-1 , d1 = 0.01 см2с−1, dс = 0.003 см2с−1, см. также § 3.2.
42