Материал: Эффект фотонной лавины в кристаллах и наноструктурах. Монография (Перлин), 2007, c.120

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

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

§ 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|| ( z1z2 )

 

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