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

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

W31,22 = S

16m e4

 

 

%

 

c

%

%

 

h3l2εL2

 

I (n1

,n2

,λ) ,

(2.16)

 

%

%

%

I

(d )

%

 

%

%

 

 

 

 

I (n1

,n3 ,λ) 2β

 

(n1,n3

,λ) =

 

 

2πn1

 

2πn3

 

π

 

2

 

2

(θ )

 

 

 

 

 

s

(θ )+ p

 

 

%

 

%

 

 

 

 

 

 

 

 

 

 

= 2β k%1dk%1

k%3dk%3 dθ

 

 

dξ ξ ×

 

 

0

 

0

π

 

|s2 (θ )p2 (θ )|

 

×

 

 

(M% (d ) (ξ))2

 

 

 

 

.

 

{ξ2 [s2 (θ) p2 (θ)]2}{[s2 (θ) + p2 (θ)]2 ξ2}

(2.17)

Коэффициент β описывает слабую интерференцию прямого и обменного вкладов в вероятность перехода 31 22. Как показывают более детальные вычисления, этот коэффициент слабо отличается от единицы. Очевидно,

что

I

(d )

% %

%

 

(d )

 

% %

 

%

 

 

 

 

 

анализ

соотношений

(2.16,

 

(n1,n3

,λ) = I

 

(n3 ,n1

,λ) . Детальный

2.17) показывает, что при λ%

d 1 имеет место следующее приближенное со-

отношение:

 

 

 

 

%

%

%

 

 

rn1n3

 

 

 

 

 

 

 

 

 

 

 

 

 

(d )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

(n1,n3 ,λ)

 

 

%

%

,

 

 

 

 

(2.18)

 

 

 

 

 

 

 

 

%

 

%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

+ n3 + b

 

 

 

 

 

 

 

где r 0.04437, b 4.88. При λ% t 30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d )

 

 

 

 

 

 

 

 

 

%

%

 

 

 

 

 

 

 

 

 

 

 

I

% %

%

 

 

 

 

r n1n3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

(n1,n3 ,

λ)

%

 

%

 

 

 

(2.19)

 

 

 

 

 

 

 

 

 

 

 

%

 

 

c)

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1 + n3 + b (λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%

I

(d )

% %

 

%

%4

. Зависимости I

(d)

где c 7.37. При очень больших λ

 

(n1,n3,λ) λ

 

%

 

%

 

%

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

от n3

иλ

при n1 = 0.5 приведены на рис. 2.2 (а, б).

 

 

 

 

§ 2.3. Скорости межподзонной релаксации электронов в квантовых ямах

Межподзонная релаксация в квантовых ямах обусловлена, в первую очередь, взаимодействием электронов с оптическими фононами. Проблемы электрон-фононного взаимодействия в системах с квантовыми ямами рассматривались во многих публикациях (см., например, [42-52]). Сложность задачи обусловлена тем, что наличие конфайнмента существенным образом модифицирует не только электронный спектр системы, но также и колебательный спектр и электрон-колебательное взаимодействие. Необходим, в частности, корректный учет взаимодействия электронов как с конфайнментными, так и с интерфейсными модами. Межподзонное рассеяние рассматривалось в работах [45, 52]. В актуальном для нас случае, когда энергия продольного оптического фонона hωL мала по сравнению с энер-

23

гией E01 нижней подзоны в яме, выражение для скорости переходов между первой и второй подзонами, полученное в работе [45], принимает вид:

(R )

 

e2mcωLa

1

 

1

 

 

 

W21

2π

2

 

2

 

 

 

 

,

(2.20)

h

ε

ε0

 

 

 

 

 

 

 

 

 

где εи ε0 – высокочастотная и статическая диэлектрические проницаемости. В работе [52] вместо (2.20) получено выражение

W (G) 4π W (R )

ε0

Λ,

(2.21)

 

21

21

ε

 

где Λ – функция кинетической энергии электрона, принимающая в актуальной для нас области энергий значения ~ 1÷2×101. Обе формулы дают оценку W21 ~ 1011с-1 .

§ 2.4. Уравнения баланса для заселенностей электронов в подзонах размерного квантования

Кинетику фотопереходов в квантовых ямах мы будем описывать с помощью уравнений баланса для населенностей электронов в подзонах. Ясно, что такой подход к решению задачи является весьма грубым, т.к. явным образом не учитывается кинетика электронов внутри каждой из подзон. Вместе с тем, имеется ряд факторов, которые все же позволяют считать данное приближение удовлетворительным. Прежде всего, это малость вероятностей запрещенных непрямых внутризонных фотопереходов по сравнению с вероятностями разрешенных прямых межподзонных переходов. Внутриподзонные и непрямые межподзонные фотопереходы приводят к появлению высокоэнергетических хвостов функций распределения электронов в подзонах. Однако число электронов в этих хвостах незначительно, как из-за слабого внутриподзонного поглощения, так и благодаря коротким временам внутриподзонной релаксации. Роль этих хвостов в принятой модели фотонной лавины мала, и можно считать, что при низких температурах практически все электроны, находящиеся в какой-либо подзоне заполняют состояния на ее дне вплоть до уровня Ферми. При составлении уравнений баланса следует иметь в виду то обстоятельство, что при больших интенсивностях накачки j концентрации электронов во всех трех подзонах размерного квантования в квантовых ямах становятся сопоставимыми. При этом следует учитывать межподзонные переходы как с поглощением, так и с испусканием фотона. Сказанное не относится к переходам 2 1, так как электроны, попавшие в подзону 2 из подзоны 1 за счет слабого поглощения света и из подзоны 3 благодаря межподзонной релаксации и оже-процессам 31 22, сначала оказываются в состояниях, далеких от дна подзоны 2. Затем за времена ~ 10-13 с они релаксируют в область

24

вблизи дна подзоны, откуда вынужденные переходы в подзону 1 с испусканием фотона невозможны из-за отсутствия резонанса.

Система уравнений для концентраций электронов в трех подзонах зоны проводимости имеет вид:

n3 = −(W31 +W32 )n3 +σ23 j(n2 n3 ) W31,22 (n1,n3 ) ,

 

&

23 j(n2 n3 ) +W32n3 + 2W31,22 (n1,n3 ) +σ12 jn1 ,

 

n2 = −W21n2 σ

(2.22)

&

 

 

n&1 = −σ12 jn1 +W21n2 +W31n3 W31,22 (n1,n3 )

c начальными условиями ni(0) = n0δ1i. В формулах (2.22) Wij (i > j) – скорости релаксационных переходов из i-й подзоны в j-ю, σji – сечения оптических переходов между j-ой и i-ой подзонами. При численных расчетах удобно пользоваться уравнениями баланса именно в форме (2.22), хотя благодаря соотношению n1 + n2 + n3 = n0 можно свести (2.22) к системе из

двух нелинейных уравнений. Оценка сечений σji производится с учетом того, что матричные элементы переходов в случае поля, направленного вдоль оси роста наноструктуры Z, для квантовой ямы с бесконечными стенками равны:

ci ,ki Z cj ,k j =

8a ij(1δij )

 

 

i+ j

δ(ki k j ),

 

π 2 (i2 j2 )

1

(1)

 

(2.23)

а типичные ширины полос межподзонного поглощения составляют 10÷30 мэВ (см., например обзор [53]). Тогда для резонансных переходов

между соседними подзонами получим σi,i+1 4 см2/(пc МВт) (a/10 нм)2. Значение σ12 для нерезонансного непрямого перехода уменьшается с уве-

личением расстройки резонанса η [см. формулу (2.6)]. При η ~ 100 мэВ величина σ12 оказывается на 4÷5 порядков меньше, чем σ23.

Вообще говоря, положение пиков и интенсивности межподзонных переходов зависят от уровня легирования квантовой ямы, в том числе, благодаря многочастичным эффектам. Теоретическое рассмотрение соответствующего круга вопросов дано в обзоре [54]. Как показывает сравнение теории с экспериментом (см., например, [55]), вклад в коротковолновый сдвиг пиков межподзонного поглощения инфракрасного света вносят эффекты деполяризации (плазмонный сдвиг), эффекты экситонного типа (взаимодействие возбужденного электрона с дыркой в основном состоянии), а также эффекты прямого и обменного кулоновского взаимодействия. В рассматриваемой нами ситуации уровень легирования, т.е. величина n0, не меняется. В принципе, перераспределение электронов между подзонами 1, 2, 3 может оказывать некоторое влияние на величины σi,i+1. Повидимому, этим эффектом можно пренебречь в случае, когда плазменная частота электронов в квантовой яме ωp мала по сравнению с зазорами ωij между подзонами размерного квантования. В случае глубоких и узких

25

квантовых ям (см. ниже обсуждение в § 2.5) это условие выполняется при всех актуальных концентрациях n0. Впрочем, сколько-нибудь корректный учет этого эффекта весьма сложен и выходит за рамки нашего анализа.

Вернемся к формулам (2.16–2.19). Если безразмерный параметр

расстройки резонанса λ% d 1, то скорость W31,22 переходов оже-типа 31 22 слабо зависит от λ% и может быть представлена в виде

 

2

a

2

 

W31.22

4mc

e

 

n1n3 γ n1n3 .

(2.24)

 

 

 

hεL

 

 

В этом случае система уравнений баланса (2.22) приводится к более простой форме:

n3 = −(W31 +W32 )n3 +σ23 j(n2 n3 ) γ n3 (n0 n2 n3 ),

 

&

 

n&2 = −W21n2 σ23 j(n2 n3 ) +W32n3 + 2γ n3 (n0 n2 n3 ) +

(2.25)

+σ12 j(n0 n2 n3 ).

 

Относительная простота системы нелинейных уравнений (2.25) позволяет провести ее качественное исследование с помощью стандартных методов нелинейной динамики. Как нетрудно убедиться, для системы (2.25) имеются две стационарных точки, а именно, стабильный узел и нестабильный узел. При σ12 = 0 (отсутствует поглощение света на нерезо-

нансных переходах между первой и второй подзонами размерного квантования) и n0γ >Wpq эти две стационарные точки вырождаются в одну при

бифуркационном значении интенсивности j = jth , где

jth

= W21(n0γ +W31 +W32 ) .

(2.26)

 

σ

23

(n γ W W )

 

 

 

0

21

32

 

При σ12 0 бифуркационная интенсивность

jth оказывается ком-

%

 

′′

. Однако для актуальных значений σpq,

плексной величиной: jth = jth

+ i jth

Wqp и n0 имеет место неравенство:

′′

jth

<< jth . Поэтому при j, близких к jth ,

особенности решений нелинейной системы (2.25) сохраняются, становясь с увеличением σ12 все более сглаженными. Эти особенности проявляются в виде резкого возрастания квазиравновесных значений n2 и n3 и убывания n1 при увеличении j вблизи jth. При интенсивностях света, близких к jth происходит также резкое возрастание времени τeq установления квазиравновесных населенностей. Для населенностей np, близких к квазиравновесным

значениям np, выполняется соотношение npnp npexp(t τc ) , где τc

можно назвать временем эквилибрации. Для систем, подобных рассматриваемым нами, обычно τeq ~ 3 5 τc . При j = jth , и при выполнении неравен-

ства n0γ >Wpq получим:

26

τc

~

1

 

n0γσ23

.

(2.27)

2W21

σ12 (W21 +

2W31 +W32 )

 

 

 

 

В отличие от пороговой интенсивности jth время эквилибрации τc зависит от сечения σ12 оптического поглощения на переходах между подзонами 1 и 2, уменьшаясь при увеличении σ12.

§ 2.5. Результаты численного решения уравнений баланса и их обсуждение

Результаты численного решения системы нелинейных уравнений (2.22) даны на рис. 2.32.6. На рис. 2.3 представлены типичные зависимости концентраций ni от времени t, прошедшего с начала импульса накачки с интенсивностью j = 0.7 МВт/см2.

Рис. 2.3. Зависимости концентраций носителей в подзонах от времени с мо-

мента включения импульса накачки (пояснения в тексте)

В расчете использовались следующие значения параметров, фигури-

рующих в правых частях уравнений (2.22): n0 = 0.25·a2, a = 0.5·106

см,

η = 0.01 эВ,

W31 = 2·1011 c1, W32 = 0.8·1012 c1,

W21 = 0.9·1012 c1, β = 1.1,

σ12 = 0.003

см2/(пс·МВт), σ23 = 2 см2/(пс·МВт).

Видно, что при t > τeq

данном случае τeq 12.5 пс) устанавливается квазиравновесное распределение электронов в подзонах 1, 2, 3, причем концентрация электронов в подзонах 2 и 3 оказывается при выбранных значениях параметров выше, чем в нижней подзоне 1.

27