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

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

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

c энергиями меньше E , в то время как основной вклад в ток вно-

сят носители, занимающие намного более мелкие состояния с энергиями вблизи Etrans . Можно привести также примеры (прямоугольное энергетическое распределение ЛС), когда времена teq _ μ и teq _ D

практически совпадают, так что квазидисперсионный режим отсутствует, см. раздел 5.4.

Следует заметить, что уравнение (3.39) описывает все указанные выше режимы переноса, т. е., в порядке убывания степени неравновесности – дисперсионный, квазидисперсионный и квазиравновесный режимы (уравнение (3.39) может описать и переход к

равновесному режиму, если заменить DF (t ) на коэффициент полной диффузии, D(t )= D0θ1 (t )+ DF (t )). Оно включает сильно не-

равновесный (дисперсионный) и квазиравновесный режимы как предельные случаи при малых и больших временах, соответственно. Можно назвать его уравнением неравновесного транспорта. Очевидно, проявления неравновесности носителей заряда значительно шире, чем эффекты дисперсионного транспорта. То обстоятельство, что квазидисперсионный режим переноса имеет признаки, характерные для обоих предельных случаев, в некоторой степени разрешает противоречия в интерпретации экспериментальных результатов ВПМ для некоторых модельных полимеров (например, поливинилкарбазол), когда транспорт характеризуется то как сильно неравновесный [1:13], то как квазиравновесный [1:14].

4.4. Диффузия, стимулированная переменным полем

Полученные выше уравнения, описывающие неравновесный транспорт (см. (3.39), (4.18)), справедливы не только в случае постоянного электрического поля F0, но и в том случае, если его на-

пряжённость F (t ) достаточно медленно изменяется со временем:

характерное время его изменения ω1 больше, чем характерное время установления квазиравновесия teq _ D . В противном случае

стимулированная полем диффузия имеет ряд особенностей, которые исследованы в данном разделе для двух предельных случаев –

181

квазиравновесного ( t > teq _ D ) и сильно неравновесного ( t < teq _ μ )

транспорта. В последнем случае о стимулированной полем диффузии можно говорить потому, что значительная дисперсия, превышающая дисперсию от обычной диффузии, может возникать без сдвига «центра тяжести» распределения носителей. Следует заметить, что условие

ωthop

1,

(4.26)

где thop v01 = ω0 exp(2γa) – характерное время прыжка между проводящими состояниями, предполагается выполненным в данном параграфе, что не исключает из рассмотрения случая ωeq _ D >1 ,

поскольку thop teq _ D .

4.4.1. Режим квазиравновесного транспорта

В этом случае ( t > teq _ D ) в уравнении (3.39б) можно считать

θ1 (t )≈ θ1 ()≡ θ , θ2 (t ) ≈ θ2 () ≡β и λ = 0 . Энергетическое распределение ловушек считаем достаточно глубоким для реализации условия N0 (x,t ) N (x,t ) . Вместе с тем, быстропеременность по-

ля приводит к тому, что N0 (x,t )t нельзя, вообще говоря, считать малой величиной в сравнении с N (x,t )t . Из уравнений

(3.31), (3.33) получено следующее уравнение для полной плотности носителей [4:11]:

 

 

 

 

 

 

 

 

 

 

N (x,t)

eq F (t )N

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

−θμ0

F (t )

 

dt F (t

)exp −ωt (t t

)

 

x2

N (x,t)= 0,

(4.27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

где

μ

eq

= μ

θ ≡ μ

τ

0

t

rel

,

ω ≡ τ 1θ2β. Характерное время измене-

 

 

 

0

0

 

 

 

t

0

 

 

 

 

 

 

ния функции N (x,t )

не может быть меньше времени освобожде-

ния

trel = τ0

θ , поскольку подавляющее

большинство

носителей

182

локализовано на ловушках. Следовательно, N (x,t )N (x,t ), где

угловые скобки означают усреднение на временном интервале от t до t +trel . Выполняя такое усреднение в уравнении (4.27), получим

 

N (x,t )

 

F (t)

N

D

(t )

2

N (x,t)= 0,

 

(4.28)

 

t

 

x

x2

 

 

 

eq

 

F

 

 

 

 

 

DF (t )= (τ0

 

2

 

 

t

′ ′

 

) .

(4.29)

 

 

 

 

 

 

 

trel )μ0

F (t )dt F (t

)exp −ωt (t t

 

 

 

 

 

 

 

0

 

 

 

 

 

 

В случае постоянного поля из (4.29) следует

 

 

 

DF

= DFeq = θμ0

2 F0

2

ωt ,

F (t ) = F0 = const.

 

(4.30)

Этот предельный случай известен как полевая диффузия, или СПД. В обратном случае достаточно высоких частот, ωtrel 1 ,

F (t ) = 0 и расплывание пакета носителей происходит, согласно

(4.26), без заметного сдвига его «центра тяжести». На освободившийся носитель действует поле со случайной фазой, поэтому прыжки носителей на соседние ЛС в обе стороны оказываются равновероятными. В случае осциллирующего поля с амплитудой

F0 , F (t )= F0 cos(ωt ), из (4.27) получена не зависящая от времени

величина

= (θμ0

 

 

2 2ωt ) 1 +(ω ωt )2

 

 

 

 

DF

2 F0

.

 

 

(4.31)

При условии

ω ω

из (4.31) получаем

D

= θμ

2 F 2

2ω , т.е.

 

 

t

 

 

F

0

0

t

результат (4.28) с заменой F0 на «действующее значение» F0 2 . Физический механизм СПД особенно нагляден в случае моноэнер-

гетических ловушек

с

энергией

 

Et .

Из

(4.30)

следует

DF = (1 2)vrellF

2 , где

vrel = v0 exp(Et

 

kT )

характерная

частота

освобождения, lF = (μ0 F0 τ0 )

1 + ω2 τ0

2

– дрейфовый сдвиг до сле-

дующего захвата. В случае

ωτ0

1 величина lF

ограничена вре-

менем захвата на ловушки

τ0 ,

lF ≈μ0 F0 τ0 ,

а в обратном случае

длина lF близка к амплитуде колебаний носителя в переменном

183

поле, lF ≈ μ0 F0 ω (предполагается, что амплитуда превышает сред-

нее расстояние между ловушками, μ0 F0 ω >

Etrans

 

1 3

 

dEg (E )

).

 

 

−∞

 

 

 

 

 

 

 

В предельном случае низких частот, ωtrel

 

1 , средняя величина

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

F (t )F (t ), а коэффициент СПД, строго говоря, зависит от времени, см. (4.29). Центр тяжести пакета x(t ) колеблется с часто-

той ω, а его дисперсия растёт вследствие стохастического разброса времён освобождения носителей с ЛС, т.е. работает механизм полевой диффузии. Однако при выполнении условий

μ0 F0 ω x ωωt (ω2 + ωt 2 ) μ0 F0t и t (θωt )(ω2 + ωt 2 )ω2 ос-

цилляциями x(t ) и DF (t ) можно пренебречь, и коэффициент DF также определяется формулой (4.30).

Результаты, полученные для случая осциллирующего поля, показывают, что СПД является доминирующим по отношению к обычной диффузии механизмом расплывания пакета носителей при условии достаточно высокой напряжённости поля и низких темпе-

ратур: F02 > 2kT (ωt eμ0 )(1 + ω2 ωt 2 ).

В рассмотренном случае координатное распределение носителей заряда представляет собой гауссовский пакет с дисперсией

d

(

t

)

=

 

(

t

)

2

1 2 =

F

 

 

x x

 

 

2D t . Частотная зависимость СПД мо-

жет быть исследована специально поставленными времяпролётными экспериментами с использованием зависящего от времени поля

F (t )= F0 + FV cos(ωt ).

(4.32)

В этом случае следует ожидать появления зависящей от частоты

дисперсии дрейфующего пакета

 

 

 

 

d (t,ω)= μ0 {2(θt ωt ) F0

2 +(FV

2 2) (1 + ω2 ωt

2 )1

}1 2

, (4.33)

Экспериментальное исследование диффузии, стимулированной переменным полем, может быть полезно для выяснения механизма

184

электропроводности в конкретном материале и уточнения транспортных параметров.

4.4.2. Режим сильно неравновесного транспорта

По причине быстропеременности величины N0 (x,t ) (см. начало

раздела 1.5.1) уравнение N0 (x,t )≈ ∂t τ(t )N (x,t ) , применимое

в режиме СНТ в случае постоянного поля, следует уточнить следующим образом:

N0 (x,t ) dτ(t )dt N (x,t )+ τ(t ) N (x,t )t −∂N0 (x,t )t . (4.34)

Решив уравнение (4.34) относительно N0 (x,t ) и используя

(3.31) (без диффузии), в случае ωt

1 можно получить следующее

уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N (x,t )

 

dτ(t)

F

(

t

)

N

D

disp

(

t

)

 

3

 

τ

(

t

)

N = 0, (4.35а)

 

 

0

 

 

 

 

 

 

2

 

 

 

t

 

 

 

dt

 

 

x

F

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

t

 

 

′′

 

 

′′

 

 

 

 

 

 

disp

 

 

2

 

 

 

 

 

)exp dt

τ(t

) .

 

(4.35б)

 

DF

(t )= μ0

F (t )dt F (t

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

В случае постоянного поля

 

DF disp (t )= (μ0 F0 )2 τ(t ), и уравнение

(4.35а) является эквивалентной формой записи уравнения в приближении СНТ (1.40) (имеет то же решение (1.41)). В случае поля

F (t ), определяемого формулой (4.32), из (4.35б) при условии dτ(t )dt 1 получаем

DF

(t )= μ0

 

τ(t ){F0

 

+(F0

 

2) 1

+ ω τ

 

(t )

1

}.

(4.36)

disp

 

2

 

2

 

2

 

2

2

 

 

 

В случае F0 = 0 (осциллирующее поле) пространственная дисперсия носителей зависит от времени по закону

d (t )= (1 2)(μ0 FV

ω)

2

ln 1

+ ω τ

 

(t) ,

(4.37)

 

 

 

2

2

 

 

На начальном (после импульса генерации) интервале времени решение N (x,t ) уравнения (4.35) в условиях времяпролётного экс-

185

перимента имеет негауссовский вид [4:11],

 

 

 

 

 

N (x,t )= 2 πμ0 FV τ(t) K0

(

2

 

x

 

πμ0 FV τ(t )),

 

 

 

 

ωτ(t ) 1,

 

x

 

> μ0 FV τ(1 ω),

 

(4.38)

 

 

 

 

 

 

 

 

где K0 (x) – функция Макдональда.

Дисперсия распределения

(4.37), d (t )(1 2 )μ0 FV τ(t ), превышает (при условии ωτ(t ) 1 ) аналогичную величину, создаваемую обычной диффузией в дисперсионном режиме, dT (t )= 2(kTμ0 e)τ(t ) (в случае экспонен-

циально распределённых ЛС, τ(t ) (v0t )α ). Однако, как следует из (4.37), при ωτ(t ) 1 СПД практически не даёт вклада в расплывание пакета носителей: d (t )(μ0 FV ω){ln ωτ(t) }12 . Следует заме-

тить, что численная величина дисперсии может быть довольно значительной. Так, при FV =104 Всм, ω =106 с1 , μ0 = 0,1 см2 В с и

ωτ(t ) 1 получаем d (t )10 мкм.

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

бой степенной зависимости μ0 ωS , s 1, которая обусловлена прыжковым характером транспорта в неупорядоченной среде [1:5].

186

ГЛАВА 5. ПЕРЕХОДНЫЙ ТОК В УСЛОВИЯХ ВРЕМЯПРОЛЁТНОГО ЭКСПЕРИМЕНТА

5.1. Переходный ток при неравновесном транспорте

При условии t θ1 (t)θ2 (t) уравнения (2.13) и (4.18) приводят к следующей формуле для переходного тока:

 

 

 

 

 

 

 

 

 

t

 

 

 

)J

)eσ0λ(t ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.1а)

 

 

 

 

j (t )= J (t,0)+ dt λ(t

(t,t

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J (t,t)=

 

eσ0

exp −Λ(t,t) ×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

F0 M (t,t)L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 Erfc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

)

 

F0μ(t )+ λ(t)(L F0 M (t,t

))

+

 

 

 

 

 

 

D(t,t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, (5.1б)

× + F0μ(t )+ λ(t)(L F0 M (t,t)) +

 

 

 

 

 

 

 

 

 

L F M (t,t)

2

DF (t )2λ(t )(D(t,t))

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4D(t,t)

 

 

 

 

 

 

4πD(t,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Результат вычислений по формулам (2.13), (4.18) значительно упрощается, если в интеграле (2.13) заменить нижний предел на

−∞ , что справедливо в случае F0 M (t,0) 4πD (t,t) 1 (это не

исключает t ttr ), то есть когда при x = 0 концентрация близка к 0.

Надо заметить, что именно при этом условии и справедливы приближения, лежащие в основе уравнения (3.39), поскольку диффузионноподобный член в этом уравнении есть лишь способ описания дисперсии носителей заряда, возникающей в процессе их дрейфа в неравновесных условиях.

Практически важным является также случай генерации носителей однородно по объёму слоя ( N0 – начальная концентрация), ко-

торый реализуется при облучении быстрыми электронами в экспериментах по нестационарной радиационной электропроводности [1:13]. В том же приближении получено [3:28]

187

junif (t )= −(eN0 2L)(∂ ∂t ) L

dx(L x)2 N (x,t ),

 

(5.2)

−∞

 

 

 

 

 

В предельном случае DF 0 , т.е. пренебрегая СПД, из (5.1) не-

трудно получить

 

 

 

 

 

t

′ ′

) , t < t*

,

 

 

(5.3)

j (t ) (eσ0 F0 L)= μ(t )−λ(t )dt μ(t

)exp −Λ(t,t

0

 

 

 

 

 

j (t )

(eσ

F

 

L)= μ(t )

1 exp

−Λ t,t

0

(t )

 

 

 

0 0

 

 

 

 

 

 

{

 

 

 

 

 

}

 

(

)

t

 

 

(

 

 

)

 

 

(

 

)

 

*

 

 

 

 

 

 

 

 

 

 

−λ t

 

 

 

 

t

 

exp −Λ t,t

,

t t

 

,

 

 

 

 

dt μ

 

 

 

 

 

 

 

 

 

t0 (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где времена t0 (t), t*

определяются из уравнений

 

 

 

 

 

t

 

 

 

 

 

 

)= L,

 

 

 

 

 

 

 

 

 

 

F0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt μ(t

 

 

 

 

 

 

 

 

 

 

 

t0 (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.4)

(5.5а)

t*

)= L,

 

′ ′

(5.5б)

F0 dt μ(t

0

 

 

Уравнения (5.3)–(5.5) были получены ранее в модели МЗ без диффузии [1:17, 30], см. уравнения (1.34) и (2.13). Далее в этой главе уравнения (5.1)–(5.5) применяются для анализа кривых переходного тока в случае нескольких энергетических распределений ловушек.

5.2. Гауссовское распределение ловушек. Особенности квазидисперсионного режима транспорта носителей заряда

Чтобы получить простое приближение к достаточно громоздкой зависимости j(t), определяемой уравнениями (5.1), рассмотрим

наиболее важный практически случай умеренной неравновесности, t teq _ μ . Оценка интегралов в уравнениях (4.18), (5.1) в этом случае

при гауссовском распределении ЛС приводит к выражениям

p(x,t)G(x,t,0),

(5.6а)

188

 

eσ F

 

 

F M (t,0)L

 

j (t )

0 0

μ(t )exp tλ(t ) Erfc

0

 

.

(5.6б)

 

 

 

 

2L

 

 

 

2

DF (t,0)

 

 

 

 

 

 

 

 

 

На рис. 5.1 показано хорошее количественное согласие результатов счёта переходного тока по приближённой формуле (5.6) (пунктирная кривая) и согласно уравнениям (5.1) (сплошная линия), а также качественное согласие с результатом эксперимента (линия с кружками) и численного моделирования (чёрные кружки).

8

teq_μ

 

 

 

 

 

t0

 

 

 

6

 

 

 

 

 

 

 

 

 

4

 

 

t1/2

 

 

j(t)

 

 

 

 

 

2

 

 

 

 

 

 

 

ttr

 

 

 

0

 

 

 

 

 

0

50

100

150

200

250

t, MKC

Рис. 5.1. Временные зависимости переходного тока в условиях времяпролётного эксперимента. Сплошная кривая – результат вычислений из уравнения (5.1), пунктирная кривая (почти сливается со сплошной) – приближённая формула (5.6), линия с кружками и точки представляют результаты соответственно эксперимента и численного моделирования [1:14]. Штрихпунктирные линии показывают способ определения времён t0 и t1 2 , σkT = 3,5 , T = 312 К

Данные времяпролётного эксперимента и результаты численного моделирования были приведены в работе [1:14] Надо заметить, что применимость уравнения (5.6) ограничена по времени сверху

величиной (2 ÷3)ttr . Послепролётная асимптотика j(t) описывается уравнением (2.15б), в котором τ(t) определяется уравнениями (1.23), (1.26), (1.29). Величиной, непосредственно измеряемой в

189

эксперименте, является параметр дисперсии сигнала переходного тока W =(t12 t0 )t12 , см. раздел 2.1.2.

Используя уравнения (5.6), нетрудно получить следующую оценку:

W L1 π

1 2

 

ttr

 

 

 

dt DF (t

) .

(5.7)

0

 

 

 

 

 

 

 

Полагая, что ttr Lμeq F0 и DF (t) DF (ttr ), можно упростить уравнения. (5.7):

W ≈ πDF (ttr ) μeq F0 L.

(5.8)

Уравнение (5.8) позволяет экспериментально определить зависимость коэффициента полевой диффузии от времени. В предельном случае ttr > teq _ D это уравнение приводится к известному вы-

ражению [1:14; 2:13]

W ≈ πDFeq μeq F0 L ,

(5.9)

которое предсказывает зависимости от толщины и напряжённости поля W L0.5 и W F00.5 , соответственно, поскольку DF F02 .

Подобные зависимости действительно были отмечены как результат моделирования методом Монте-Карло в модели гауссовского беспорядка Бэсслера [1:14; 4:3] и получены экспериментально [3:3; 2:16; 5:1] при достаточно большой толщине образца и достаточно малой напряжённости поля. Это означает, что время пролёта достаточно велико для того, чтобы не только подвижность, но и коэффициент стимулированной полем диффузии успели достичь своих установившихся значений. Напротив, при уменьшении величин L

и F01 величина W перестаёт зависеть от данных параметров [2:13, 15]. Причина – в том, что коэффициент СПД ещё не достиг квази-

равновесного значения. Аппроксимируя зависимость

DF (t ) сте-

пенной

 

функцией, DF (ttr ) ttr S , и учитывая,

что

ttr L μF0 и

DF F0

2 , из уравнения (5.7) получаем W (L F0 )(S 1)

2 . Очевидно,

в случае

S =1 параметр W не зависит ни от

L , ни от F0 , что

 

 

190