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

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

фермиевскими функциями распределения, соответствующими мгновенным концентрациям nc (t), nc1 (t) и p(t) электронов и дырок. Так в форму-

лах (5.26) фигурирует функция распределения fcc электронов в зоне c, со-

 

1

 

ответствующая энергии εc (kcc(l ) ) в точке l-фотонного резонанса kcc(l )

между

1

1

 

зонами c и c1:

 

 

fcc1 ={exp (εc (kcc(l )1 ) μc (nc ))

kBT +1}1 .

(5.29)

Аналогичным образом вводятся функции распределения fc и fp электронов в зоне c и дырок в зоне v, соответствующие энергиям εc (k(vcn) ) и εv (k(vcn) ) в точке n-фотонного резонанса k(vcn) между зонами v и c. В случае

относительно больших квантов возбуждающего света (hω t 1 эВ) и не очень малых продолжительностей световых импульсов (τI t 1 пс) такое приближение оказывается не очень плохим с учетом малых сечений внутризонного поглощения света и коротких времен внутризонной релаксации

электронов и дырок. Разумеется, в области бóльших длин волн (hω d 0.1 эВ) пользоваться этим приближением нельзя (см. [67]).

Для вычисления химических потенциалов μc и μp ( μp > 0) для элек-

тронов и дырок, фигурирующих в формулах для функций распределения, воспользуемся трансцендентными уравнениями:

 

 

 

nc (t)

−Φ

 

(μ* )

= 0,

 

p(t)

−Φ

(μ* )

= 0

,

(5.30)

 

 

 

 

 

 

 

 

 

 

Nc (T )

 

1/ 2

 

 

c

 

 

 

 

N p (T )

1/ 2

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где, μ

= μ

c, p

k T ,

N

c

и N

p

обозначают эффективные плотности состоя-

c, p

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ний в зоне проводимости и валентной зоне:

 

 

 

 

 

 

 

 

Nc

 

2πmckBT

3/ 2

 

 

 

 

2πmvkBT

 

3/ 2

 

(5.31)

 

 

= 2

(2πh)

2

 

 

,

 

N p = 2

(2πh)

2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф1/2(x) – интеграл Ферми-Дирака:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ1/ 2 (x) =

2

 

 

z1 2dz

.

 

 

 

 

(5.32)

 

 

 

 

 

 

π 0 1+ exp(z x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, требуется получить самосогласованное решение системы нелинейных дифференциальных уравнений (5.26) и трансцендентных уравнений (5.30). Сложность этой задачи делает невозможным качественный анализ динамики системы типа того, что был проведён в § 2.4, и получение простых формул для пороговых интенсивностей света и времен установления квазиравновесных населенностей в электронной системе. Поэтому будем использовать результаты численных расчётов, которые приводятся в следующем параграфе.

83

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

На рис. 5.6-5.7 приведены результаты численного решения уравнений (5.26-5.30) для случая l = 2. На рис. 5.6 (а, б, в) сплошными линиями даются зависимости квазиравновесных концентраций электронов в обеих зонах проводимости (nc, nc1 ) и дырок в валентной зоне (p) от интенсивности све-

та. На этих рисунках штриховыми линиями показаны те же зависимости, рассчитанные с теми же параметрами, но без учета трехфотонных процессов оже-типа. На рис. 5.6 (г) дается зависимость времени установления квазиравновесных населенностей τeq от интенсивностей. Отметим следующие качественные особенности этих зависимостей. Имеется четко выраженное скачкообразное увеличение концентраций неравновесных носителей с ростом j вблизи порогового значения (при выбранных значениях параметров jth 2.7 1010 Вт/см2). При этом значении интенсивности происходит резкое возрастание времени τeq. Указанные особенности являются типичными для эффекта фотонной лавины, причем в рассматриваемом случае реализуется многофотонная лавина. Видно также, что концентрации неравновесных носителей при j > jth оказываются на 1-2 порядка выше, чем при обычных многофотонных межзонных переходах с вероятностями, рассчитанными при тех же значениях параметров.

На рис. 5.6 (г) видно, что при j ~ 1012 Вт/см2 время τeq ~ 1 пс. При уменьшении интенсивности j времена τeq быстро возрастают и для j ~ 1011 Вт/см2 оказываются уже в наносекундном диапазоне, тогда как максимальный практический интерес представляет взаимодействие с твердыми телами более коротких импульсов. Поэтому следует также рассмотреть генерацию неравновесных электрон-дырочных пар при продолжительностях лазерных импульсов τp, малых по сравнению с τeq.

На рис. 5.7 (а, б, в) показаны зависимости от интенсивности концентраций nc, nc1 и p, возникших в результате действия импульсов с продол-

жительностями 1, 5, 10 и 50 пс. Приведенные зависимости построены для импульсов, имеющих гауссову форму. Видно, что и здесь имеется область чрезвычайно резкого роста концентраций с увеличением интенсивности (при росте j на 30-40% концентрации возрастают на 4-5 порядков). При этом возможна ситуация, когда совсем небольшого увеличения интенсивности света оказывается достаточным для перехода от обратимого предпробойного возбуждения материала к его деструкции, возникающей при концентрациях неравновесных носителей, превышающей 1019 см3. Если обозначить через jthинтенсивность, вблизи которой скорость изменения,

например, p(j) максимальна, то имеем следующую приближенную зависи-

мость: jτ 13 .

th p

84

85

Рис. 5.6. Зависимости квазиравновесных концентраций электронов в зонах c (a), c1 (б), дырок в зоне v (в) и времени установления квазиравновесных концентраций (г) от интенсивности света для случая l = 2; штриховая линия (а, б, в) – те же зависимости, рассчитанные без учета процесса c1 + 3hω ccv ; в расче-

тах использовались следующие значения параметров: Eg = 5.5 эВ, Eg= 2.28 эВ,

Ñω = 1.17 эВ, mc =0.9 m, m = 0.03m , mv = 0.9 m,

d

c

= 3×10-4

см3/пс, dc=0.01

c

 

 

 

1

 

1

 

 

см3/пс, Wc1c =1.0 пc-1, T = 103 K, ε= 4.5, ε0 = 12.4, p0 = 1014 см-3, n0 = 1014 см-3

86

87