Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Никитенко Нестационарные процессы переноса и 2011

.pdf
Скачиваний:
73
Добавлен:
12.11.2022
Размер:
4.79 Mб
Скачать

I (t )≈ −(pi η0

κr03 eμ0 )cd ( f )

d

 

1

,

(6.26)

dt τ(t )

 

 

 

 

где множитель cd ( f ) зависит от размерности и поля. Нетрудно

показать, используя (6.12) и (6.25), что в отсутствие внешнего поля вероятность выживания пары и интенсивность люминесценции не

зависят от размерности: Ω(t )= exp w0 (t )3 ,

I

(

t

)

≈ −p η

piη0

 

dw0 (t )

exp w

(

t

)

3

,

f <<1 .

(6.27)

 

 

 

 

i 0

3

 

dt

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (6.27) следует

 

 

d 1

 

 

 

 

 

 

 

 

I (t )≈ −(piη0 κr03

 

3eμ0 )

,

 

 

 

f

<<1 , t >> tin .

(6.28)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt τ(t )

 

 

 

 

 

 

 

 

В случае экспоненциально распределённых ЛС получим степенной закон убывания: I (t) t(1) . В частности, при α = 0,5 I (t ) t32 ,

что обычно объясняют диффузионным разделением геминальных пар в режиме нормального транспорта [1:80]. Как видно, такая зависимость может возникать как проявление дисперсионного транспорта, что естественно для неупорядоченной среды.

6.4.2. Температурная зависимость кинетики люминесценции

Скорость реакции, то есть число рекомбинационных событий в единицу времени на единицу объёма, можно вычислить, очевидно, как производную по времени от плотности геминальных пар, сохранившихся к моменту t после импульсной генерации при t = 0 . Cкорость реакции пропорциональна интенсивности люминесценции I(t). Используя уравнения (6.12), (6.23), в предельном случае ξ→1, a 0 (модель Онзагера), получим (при F0 = 0 , но с учётом

диффузии)

I (t )= −2πη0 pi

{limx0 p0 (x,t ) w2 (t ) }.

(6.29)

t

 

 

 

В случае короткого импульса генерации эта зависимость имеет резкий максимум, соответствующий времени наиболее интенсивной БР. С ростом напряжённости внешнего поля возрастает Ω

221

доля разделяющихся пар, что приводит к полевому гашению люминесценции в случае d = 2,3 .

При больших временах ( t >> tin ), как было показано выше, см.

уравнения (6.27), (6.28), интенсивность люминесценции определяется быстротой опустошения глубоких ловушек, что даёт степен-

ной закон убывания t(1) в случае экспоненциально распределённых ЛС либо достаточно близкую к степенной временную зависимость в случае гауссовского распределения. Такая ситуация, повидимому, реализуется в сопряжённом полимере MeLPPP [1:55; 6:14]. Надо заметить, что образование геминальных пар в сопряжённых полимерах – не единственный механизм фотогенерации носителей заряда (большая часть экситонов за очень короткое ( < τ0 ) время распадается на свободные заряды, двигаясь по делока-

лизованным состояниям вдоль полимерной цепи) [6:15, 16].

103

T, K

 

 

298

102

203

 

83

I(t)101

t-3/2

100

10- 1

10-4

10-3

10-2

t, c

Рис. 6.4. Временные зависимости экспериментально измеренной задержанной флуоресценции в сопряжённом полимере MeLPPP и темпа геминальной рекомбинации [6:19], вычисленные согласно (6.29), в зависимости от темпера-

туры. ω0 =1012 с-1, остальные параметры показаны на рисунке

Однако геминальные пары образуются в результате перескока на другую (соседнюю) цепь. По-видимому, причиной задержанной флуоресценции в MeLPPP является именно близнецовая (геминальная) рекомбинация [1:55]. Оптическое излучение генерирует

222

достаточно близкие (r0<6 нм <<rc) геминальные пары. Распределение ЛС g (E ) описывается гауссовской функцией (1.84), причём значения ширины распределения σ = 0,055 эВ и параметра локализации M0 γ3 = 0,0005 близки к величинам, применявшимся ранее к

описанию данного материала [6:17, 18]. Предложенная модель [6:19] позволила объяснить экспериментальные данные по прибли-

зительно

степенному

спаду задержанной флуоресценции,

I (t) tn ,

n =1,1÷2,7

при T =80 ÷300 К, как результат БР, кон-

тролируемой освобождением «близнецов» с ловушек, распределённых по энергии, см. рис. 6.4. При низкой (T =80 К) температу-

ре получено I (t ) ln (ω0t )4 t1 и, таким образом, n 1,1 , что характерно для режима прыжков вниз по энергии [1:48; 6:20, 21].

6.5. Ток поляризации близнецовых пар

6.5.1. Бездиффузионное приближение

Используя уравнения (6.10) и (6.11) и безразмерные величины

(6.14), нетрудно показать, что в случае

d =1,2,3 ток поляризации

геминальных пар можно вычислить следующим образом [6:11]:

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

j (t )= (epi r0

κ)π2d 2

 

d γγ∫ dxxd ϕ(x, γ,t ), d

= 2,3 ;

(6.30а)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 1

 

0

 

 

 

 

 

 

 

 

 

(

 

)

 

(

i 0

 

)

 

 

 

(

 

)

 

 

(

 

)

d =1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

t

 

=

 

ep r

κ

 

t

0

dxx

ϕ

 

x,1,t

 

−ϕ

 

x, 1,t

,

(6.30б)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае сферически симметричной генерации при достаточно малых значениях параметров x0 и f существует временной ин-

тервал, на котором величина j (t ) отрицательна. Происхождение

данного эффекта нестационарной отрицательной фотопроводи-

мости (НОФП) легко понять, заметив, что в предельном случае низких температур и слабых полей все пары рекомбинируют и

P()= ∫ dt ' j (t ')= 0 [6:1, 1:83]. В этом случае

0

j (t ) dτ(t)dt t(1−α) при t << tin , где характерное время БР tin

223

(оно же время инверсии знака j (t )) определяется условием (6.7), и j (t ) dτ1 (t )dt t(1) при t >> tin . Рис. 6.5 показывает влияние величины F0 на характер временной зависимости j (t ) при d =1,2,3 .

Рис. 6.5. Полевая зависимость низкотемпературной кинетики тока поляризации геминальных пар. κr03 eμ0τ0 = 4 104 , kTE1 = 0,5, κr02 F0 e : 1 – <<1; 2 – 1,05;

3 – 4

Как видно из рисунка, при низких температурах и слабых полях понижение размерности не приводит к существенному изменению

формы кривых j (t ). В то же время с ростом напряжённости внеш-

него поля характер временной зависимости переходного тока существенно меняется (знакопеременность исчезает). Физическая причина в том, что с ростом напряжённости внешнего поля возрастает вероятность разделения пары, а также уменьшается характерное время разделения. Поскольку разделившиеся носители движутся в направлении внешнего поля, временной интервал существования НОФП сокращается, и в случае, если напряжённость прило-

женного поля превышает критическую величину F (r0 , d ), ток по-

ляризации всё время остаётся положительным. В бездиффузионном приближении (которое применимо при r0 << rc , для величин кри-

тической напряжённости F получено [6:11]:

224

κr02 F e =1,07 , d = 3 ; κr02 F e =1,04 , d = 2 ; κr02 F e =1 , d =1 . (6.31)

6.5.2 Изотропная проводимость (слабое поле)

Рост температуры, начального разделения пары либо напряжённости внешнего поля повышают вероятность разделения пары, поэтому при достаточно больших значениях указанных параметров переходный ток монотонно убывает, см. рис. 6.6; при этом дальняя асимптотика тока, как и начальная, определяется лишь энергетиче-

ской релаксацией носителей заряда: j (t) Ωdτ (t)dt Ωt(1α) , t >> tdis . Из уравнений (6.10), (6.11), (16.19), (16.20) следует (при a = 0 )

j (t )= pi

2π (erc )2

 

d

3

χ(x,t)

 

3

kT

F0

 

dxx

(6.32)

 

 

 

dt 0

 

 

 

Вычисления показывают [6:2, 6:12], что в случае слабого внешнего поля НОФП может возникать в системах с пониженной размерностью при значительно более высоких температурах, чем в системах

с

 

 

изотропной

проводимостью.

Так,

должно

быть

x

=κkTr

e2 < 0,16 при d =3 , и x

< 0,29 при d =1.

 

0

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

-4

 

Lg ( j /

j0 )

 

 

 

 

 

 

 

 

 

 

Рис. 6.6. Переходный ток, кон-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тролируемый близнецовой реком-

-8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

бинацией, в случае изотропной

 

 

 

2

 

 

3

 

 

 

 

 

 

 

 

 

проводимости. Параметры:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-12

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

(e2 E1 )3 κ2eμ0τ0 =103 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

κr E e2

= 0,4, j

= ep μ F .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

9

11

13

 

0 1

0

i 0

0

-13

 

 

 

 

kT E1 : 1 – 0,25; 2 – 0,35; 3 – 0,45

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

 

 

L g(ν0t)

Надо заметить, что в случае

-9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

достаточно больших значе-

 

 

Lg (j

j0 )

 

 

 

 

 

 

 

 

 

 

ний r

0,2 ,

когда ток по-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

ляризации монотонно убывает, его временную зависимость с хорошей точностью можно описать приближённой формулой [6:22]:

j (t )epiμ0 F0Ω(t )dτ(t ) dt .

(6.33)

225

t <tin

Для выявления экспериментальных условий, при которых эффект НОФП должен проявляться наиболее ярко, проведено также исследование зависимости максимального отрицательного тока j и

времени его достижения t от размерности, температуры, и параметров спектра ЛС. Анализ полученных результатов показывает, что абсолютная величина j будет наибольшей в материале с квазиодномерной проводимостью при температуре, меньшей критического значения: T Tkp (F0 , r0 ) [6:12]. Однако при малых значениях

начального разделения пары возникают экспериментальные трудности, связанные с малостью времени t .

В случае непрерывной генерации пар ( g0 штук в единице объёма в единицу времени) при t 0 нетрудно получить из уравнения (6.32):

t

jнепр(t )= (g0 pi )dt '

0

Рис. 6.7. Временные зависимости тока поляризации геминальных пар

врежиме непрерывного облучения

синтенсивностью g0 , нормиро-

ванные на j0 = eg0μ0τ0 F0 . Значения начального разделения r0 : 1 – 3 нм, 2 – 6 нм, 3 – 9 нм. Остальные пара-

метры:

α = 0,5 ,

ν0 =1010 с-1,

μ0τ0 = 7 1015 см2/В. Кривая 4 отвечает случаю гауссовского распределения значений r0 :

n(r0 ) = (4πr02 π32b3 )exp(r02 b2 ).

Кривая 5 построена для предельного случая низких температур, т.е. отсутствия диффузии ( r0 = 3 нм).

Стрелками показаны времена tin и tdis для r0 = 6 нм

j (t t ')=

2π

g0

(erc )2

3

χ(x,t ).

(6.34)

3

kT

F0 dxx

 

 

0

 

 

 

В случае непрерывного облучения материала ток возрастает при и t >tdis по асимптотическому закону jas τ(t) tα . В ши-

226

роком промежуточном интервале времён рост тока замедляется вследствие БР, а эффекту НОФП соответствует появление на кри-

вой j(t) локального максимума (кривая 1 на рис. 6.7) [6:23]. Сле-

дует заметить, что в реальных условиях влиятелен разброс значений r0 , который значительно нивелирует влияние БР на зависи-

мость j (t ) (см. кривую 4 на рис. 6.7). Кривая 4 соответствует случаю распределения начальных разделений согласно формуле n(r0 ) = (4πr02 π32b3 )exp(r02 b2 ), где b = 6 нм [6:10, 24].

6.5.3. Одномерная проводимость

Температурная зависимость тока поляризации геминальных пар в одномерном случае имеет те же качественные закономерности, что и в случае d = 3 . Параметрическая область существования

НОФП в случае d =1 показана на рис. 6.8. в осях x0 = κkTr0 e2 - f = κr02 F0 e сплошной линией. Для сравнения, пунктирные линии

на том же рисунке обозначают соответствующие результаты для трёхмерной (изотропной) проводимости.

Рис. 6.8. Параметрическая область

 

κkTr0 e2

 

 

 

 

 

существования эффекта неравновес-

0,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ной отрицательной фотопроводимо-

 

 

 

 

 

 

 

 

 

 

сти. Сплошная линия – d =1 (квази-

0, 2

 

 

 

 

 

 

 

 

 

одномерная проводимость), пунктир

 

 

 

 

 

 

 

 

 

d = 3 (изотропная проводимость).

0,1

 

 

 

 

 

 

 

 

 

κ = 4πεε0

 

 

 

 

 

 

κr

2F e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следует отметить, что в режи-

 

 

 

 

 

 

 

0

0

 

0,5

1

 

 

 

ме слабого ( f <<1 ) электри-

 

 

 

 

 

 

 

 

 

 

 

 

 

ческого поля эффект НОФП в случае d =1

может существовать

при значительно более высоких значениях температуры, чем в случае d = 3 (критические значения x0 равны соответственно 0,29 и

0,16) [6:2, 12]. Физическая причина такого различия состоит в том, что диффузионная составляющая тока в случае d =1 играет менее заметную роль, чем в случае d = 3 .

227

6.6. Вероятность разделения пары (квантовый выход)

6.6.1. Анизотропная проводимость

Кинетика БР в анизотропной среде рассмотрена выше в предельных случаях квазиодномерной либо квазидвумерной проводимости. В общем случае анизотропной электропроводности, даже для вероятности разделения пары в литературе известны только численные решения [1:74; 6:25]. C использованием приближения предписанной диффузии [6:26], в работе [6:27] получено следующее аналитическое выражение:

 

e2

3

 

D

 

 

 

 

3

x02 j

 

1 2

 

 

 

 

 

 

 

 

 

i

 

D D D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3kT i=1 κi

 

1

2

3

=1 Dj

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

,

(6.35)

Ω= exp

 

 

 

 

 

 

3

 

 

 

3

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

1 2

 

 

exp

 

e

 

 

x F

+

 

x0i

 

D F 2

)

 

 

 

 

2kT

 

 

 

 

 

 

 

 

i=1

0i i

 

i=1

D

(i=1 i i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где κi = 4πε0εi ,

Di ,

εi

 

и

Fi

значения коэффициента диффузии

носителей в проводящих состояниях, диэлектрической проницаемости и напряжённости поля вдоль координатных осей, x0i – началь-

ные координаты подвижного заряда геминальной пары. Направления координатных осей xi выбраны вдоль (совпадающих, по пред-

положению) главных осей тензоров диэлектрической проницаемости и коэффициента диффузии. В предельном случае отсутствия поля получено выражение

 

 

 

e2 3

D

 

3

x02 j

1 2

 

Ω

(0)

= exp

 

 

i

D D D

 

 

 

 

, (6.36)

 

 

κ

 

D

 

 

 

3kT i=1

 

 

1 2 3

j=1

 

 

 

 

 

 

 

 

 

 

i

 

 

 

j

 

 

 

которое в изотропном пределе переходит в известный результат

(0)

= exp(e

2

κkTr0 ),

r0

 

3

2 1 2

Онзагера, Ω

 

=

x0 j .

 

 

 

 

 

j=1

 

В реальных условиях начальное распределение термализованных носителей представляет собой результат усреднения по ансамблю генерированных близнецовых пар. Рассмотрено два конкретных

228

начальных распределения, представляющих собой два предельных случая: 1) сферически-симметричное распределение с начальным разделением r0 (6.3) (при d = 3 ) и 2) эллиптическое начальное рас-

пределение

ρ(r0 ,0)= δ(u u0 ) 4πu0

2 ,

(6.37)

где u = (u12 +u22 + u32 )12 , ui = (xi kT 4πε0e2 )(ε1ε2ε3

u0 = (kT 4πε0e

2

 

ε3

3

2

3

1 2

 

) (ε1ε2

3)

(xi0

Di )Dj

εj

 

 

 

 

i=3

 

j=1

 

 

3

3Di )Dj εj , j=1

параметр, характеризующий размер эллипсоида начальной генерации.

Рис. 6.9. Температурная зависимость квантового выхода БР в случае анизотропной проводимости при F0 = 0 . Сплошные линии – сфериче-

ская генерация (a=r0), штриховые – эллипсоидальная (a=r0eff) . На вставке: зависимость энергии активации Ea от степени асимметрии D 7/ D ^:

1 – 1; 2 – 16; 3 – 36

Распределение (6.37) возникает, если для «горячих» носителей реализуется та же анизотропия подвижности, что и для термализованных. На рис. 6.9 представлены графики температурной зависимости вероятности разделения в отсутствие поля, см. уравнение (6.36), усреднённое по начальным координатам x0i согласно рас-

пределениям (6.3) и (6.37).

 

 

ε1 = ε2 = ε3 = ε,

Вычисления

выполнены

для

случая

D1 = D2 D , D3 D , при различных значениях степени асиммет-

229

рии коэффициента диффузии (и подвижности) D D . Следует за-

метить, что в данном случае u

0

= r eff

r ,

 

 

 

 

 

 

0

c

 

 

 

3

 

 

 

1 2

 

 

 

 

 

 

 

2 } .

 

 

 

r0eff

={i=1 (D1 + D2 + D3 ) 3Di x0i

Как и в изотропном случае,

температурная

зависимость

квантового

выхода,

Ω

(0) exp(Ea kT ),

характеризуется определённой энергией ак-

тивации Ea , которая,

однако, сильно зависит от степени асиммет-

рии D D , особенно в случае эллипсоидальной начальной генера-

ции. Действительно,

при

D

 

D → ∞ из (6.36)

и (6.37)

следует

Ea

(1 3 3 )(e2 κr0eff )(D

D ). Как и следует ожидать, при пере-

ходе к квазиодномерному случаю Ω(0) 0 , D

D → ∞.

 

В случае слабого внешнего поля выражение (6.35) можно раз-

ложить в ряд по степеням Fi . В приближении, линейном по полю, получено

 

 

 

 

 

 

 

Ω

≈ Ω

 

(0) 1 + S

(

F

,

 

 

(6.38)

 

 

 

 

 

 

 

 

 

 

i )

 

 

 

 

 

 

 

 

e

3

3

Di

 

 

 

 

3

2

1 2

3

3

2

 

1 2

S (Fi

)=

 

 

(D1D2 D3 ) 1 2

 

Di Fi

)

+ ∑ x0i Fi

x0i

 

 

.

 

(kT )2

 

 

Di

 

 

6

i=1

κi

 

 

 

(i=1

 

i=1

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Даже в случае слабого внешнего поля величина Ω

 

 

(6.39)

зависит не

только от угла между направлением поля F0 и вектором начального разделения r0 , но и от углов между F0 и главными осями тензора D . В изотропном случае из (6.39) следует точный результат модели Онзагера: S = e3 F0 2κ(kT )2 .

Следует заметить, что не только для изотропного и эллипсоидального начальных распределений, но и для любого неполярного

3

начального распределения, xi0 Fi = 0 , величина S (Fi ) опреде-

i=1

ляется первым слагаемым в (6.39) и, таким образом, зависит только

230