С. Д. Ханин
МОДЕЛИ ЯВЛЕНИЙ ПРЫЖКОВОГО ПЕРЕНОСА И МЕТОДЫ АНАЛИЗА ЭЛЕКТРОННЫХ СВОЙСТВ НЕУПОРЯДОЧЕННЫХ СИСТЕМ С СИЛЬНОЙ ЛОКАЛИЗАЦИЕЙ НОСИТЕЛЕЙ ЗАРЯДА
Развиты модельные представления о нестационарных и неравновесных электронных явлениях в неупорядоченных диэлектриках. Показаны возможности экспериментального определения механизмов локализации и малой подвижности носителей заряда.
Исследования электронных свойств неупорядоченных систем занимают в последние десятилетия одно из центральных мест в физике конденсированного состояния вещества. Это обусловлено как их практической важностью, так и проблемами в развитии теории конденсированного состояния применительно к системам, в строении которых отсутствует дальний порядок [1-3]. Здесь приходится отказаться от использования ряда основополагающих понятий и представлений зонной теории. Так, понятия квазиимпульса, зон Бриллюэна лишаются смысла; в некристаллических телах не существует строго запрещенных зон энергий. Единственным сохраняющим свою силу в теории понятием является плотность электронных состояний. Это позволяет, в частности, сохраняя для некристаллических тел зонную модель энергетического спектра электронов, определить квази-запрещенные области энергии, где плотность состояний меньше, чем в разрешенных зонах, и сами состояния локализованы.
Физика неупорядоченных систем указывает на то, что при нахождении уровня Ферми в области локализованных состояний система обладает неметаллическими свойствами, причем для широкого круга объектов характерен прыжковый механизм переноса, характеризующийся низкой подвижностью [4-6].
При всей важности знаний об энергетическом спектре, характере состояний, механизмах локализации носителей заряда и кинетических параметрах вещества возможности их экспериментального определения остаются весьма ограни-
ченными. Особенно сложной эта задача становится при анализе систем с сильной локализацией носителей заряда, когда одновременно присутствуют и взаимосвязаны различные механизмы локализации.
В настоящей работе развиваются модели явлений прыжкового переноса, открывающие возможности экспериментального определения сильной электрон-фононной связи (поляронного механизма локализации), плотности электронных локализованных состояний и подвижности носителей заряда в неупорядоченных системах. Возможности предлагаемых экспериментальных методов иллюстрируются результатами анализа кинетических свойств аморфных оксидов переходных металлов.
Наличие поляронов малого (сверхмалого) радиуса предлагается выявлять на основе анализа частотной зависимости проводимости в режиме мультиплетной прыжковой электропроводности при приложении к образцу сильного постоянного электрического поля. Неупорядоченным системам свойствен эффект уменьшения дифференциальной проводимости с увеличением напряженности электрического поля, предшествующий ее экспоненциальному росту (рис. 1) [7]. Этот эффект связан с перестройкой путей протекания тока и перераспределением электронной плотности.
На рис. 2 представлены соответствующие экспериментальные результаты по электронной проводимости аморфного оксида тантала на постоянном токе в сильных электрических полях.
Рис. 1. Схематическое изображение перестройки пути протекания тока
в неупорядоченной системе в сильном электрическом поле (а) и качественные зависимости проводимости от напряженности поля (б) при изотропном (1) и направленном (2) протекании. Сплошной кривой показана результирующая зависимость у(Е)
1 2 3 4 Е-10-6, В-см1
Рис.2. Зависимость проводимости оксида тантала от напряженности электрического поля при различных температурах: 298 К (1); 323 К (2); 343 К (3); 363 К (4)
Можно полагать, что влияние постоянного смещения Е на степенную частотную зависимость прыжковой проводимости контролируется безразмерным
барометрическим параметром ви/кт' где
и — падение напряжения на кластере состояний, участвующих в прыжковом переносе, размером Ь, равное вЬ. При этом сама частотная зависимость у(Е,а) описывается следующей формулой:
у(Е,о) = у(0,о)/ (вЕЬ/кТ ). (1)
Входящая в это выражение безразмерная функция /(х) ^ 1 при х = вЕукГ ^ 0 .
Ясно, что соотношение (1) может быть справедливо лишь в случае кластеров больших размеров — в режиме мульти-плетной прыжковой электропроводности,
_ 1
когда Ь >> N 3 (здесь N — концентрация локализованных состояний, участвующих в прыжковом переносе). Кроме того, напряженность электрического поля Е не должна быть слишком велика, так, чтобы _ 1
произведение вEN 3 было соизмеримо с кТ, так как в противном случае (на участке экспоненциального роста у(Е,0)) токопрохождение является направленным и постоянное поле не оказывает влияния
на пространственное распределение носителей заряда.
Таким образом, можно ожидать влияния постоянного электрического поля на частотную зависимость прыжковой электропроводности на достаточно низких частотах и существования граничной частоты, начиная с которой эффект постоянного поля перестает проявляться.
Граничная частота со0 может быть определена из следующего уравнения [8, 9]:
сеЕф0) = 1 кТ '
(2)
Ь(со) = N- ^
1п
- Я
(3)
где показатель I в приближении эффективной среды равен единице, а в скейлин-говой теории отождествляется с критическим индексом для среднеквадратичного смещения носителя при случайном блуждании по узлам перколяционной сетки в допороговой области и равен 2у или 2у - в (здесь V = 0,9 и в = 0,4 — крити-
ческие индексы корреляционной длины и мощности бесконечного кластера, соответственно).
Подставляя (3) в (2), приходим к искомому выражению для граничной частоты с:
ю0 = Жс ехр< aN
Уз
(
с
eEN
-1з
кТ
(4)
где С — численный коэффициент.
Действительно, при со > с0 размер
проводящего кластера Ь(со) настолько
мал, что неоднородность распределения электронов внутри него не существенна.
Поэтому при СО>Щ) постоянное смещающее напряжение слабо влияет на проводимость на переменном токе, и у(Е,ю) = ^(0,ю). И, наоборот, при со <с0 кластер достаточно велик, для того чтобы перераспределение носителей заряда в нем, под воздействием постоянного смещения, сильно влияло на частотную зависимость прыжковой проводимости.
В трехмерной модели мультиплетной прыжковой электропроводности характерный размер кластера локализованных состояний, участвующих в переносе, зависит от частоты следующим образом [6]:
Рассмотрим температурную зависимость граничной частоты со0. Как следует из соотношения (4), она может быть различной, в зависимости от характера изменения критической вероятности прыжка Жс (минимальной вероятности перескоков в бесконечном кластере состояний) с температурой, т. е. типа прыжков носителей заряда между локализованными состояниями. В случае отсутствия значимого барьера при межузельном перескоке (слабой электрон-фононной связи), когда Жс слабо зависит от О, как следует из (4), со0 падает с ростом Т. В случае же надбарьерных перескоков (при сильной электрон-фононной свя-
зи), когда Жс ~ ехр| -
Е,.
'кТ
, где
Е, —
высота потенциального барьера (энергия активации межузельного поляронного перескока), температурная зависимость С определяется двумя конкурирующими факторами: экспоненциальным ростом Wc с увеличением Т и падением показателя экспоненциального множителя в формуле (4). При достаточно большой высоте барьера активационный рост Жс (Т) становится доминирующим фактором, и с ростом Т граничная частота с 0 увеличивается.
На рис. 3 представлена экспериментальная зависимость проводимости аморфного оксида тантала от частоты без приложения и с приложением постоянного смещения. Как видно из рисунка, она отвечает предсказаниям модели: экспериментальные данные обнаруживают на
низких частотах эффект падения проводимости при приложении постоянного поля напряженности, соответствующей области падения статической прыжковой проводимости, и наличия граничной частоты, начиная с которой данный эффект уже не проявляется.
Рис.3. Частотная зависимость проводимости оксидной пленки Та2О5 толщиной 170 нм без постоянного смещающего напряжения (1, 3) и с постоянным смещающим напряжением величиной 50 В (2, 4) при различных температурах: 298 К (1, 2) и 358 К (3, 4)
Таким образом, характер температурной зависимости граничной частоты а>0 может сигнализировать о сильной элек-трон-фононной связи и, в частности, о наличии в системе поляронных состояний.
Как показывают представленные данные, граничная частота а>0 для аморфного оксида тантала возрастает с увеличением температуры. Это указывает на сильную электрон-фононную связь в аморфном металлооксиде. В совокупности с другими результатами это свидетельствует о наличии поляронных состояний, что является, наряду со структурной разупоря-доченностью, причиной сильной локали-
зации носителей заряда в исследуемом материале.
Плотность локализованных состояний g(E) в неупорядоченных материалах с по-ляронным и андерсоновским механизмами локализации может быть восстановлена из температурных и концентрационных зависимостей статической прыжковой проводимости. Предлагаемая в данной работе методика экспериментального определения g(E) основывается на модельных представлениях о том, что в широком температурном диапазоне наблюдается как Л-протекание в области сравнительно высоких температур, так и Л-£-протекание при низких температурах (рис. 4), причем температурная зависимость прыжковой проводимости с переменной длиной прыжка
_ 1
подчиняется закону Мотта «Т 4» [10, 11]. Соответствующие этим представлениям экспериментальные данные для аморфных оксидов тантала и ниобия представлены на рис. 5.
Исходя из этой модели на основании температурной зависимости проводимости можно по величине параметра в законе Мотта Т0 = р/kg(^)а3 (здесь а — радиус локализации электронных состояний, р — определенный численный коэффициент) [4] восстановить плотность состояний g(p) на уровне Ферми, находящемся в области локализованных состояний. С целью определения g(E) положение уровня Ферми следует изменять, сканируя им распределение плотности состояний. Это можно осуществлять посредством направленного изменения состава вещества.
Сравнивания величины энергий активации прыжковой электропроводности в режиме Л-протекания, в которую вносит вклад энергия связи поляронов малого (сверхмалого) радиуса и разность энергий уровня максимума плотности состояний и уровня Ферми, можно определить последнюю. Полученные таким образом данные позволяют рассчитать искомую плотность состояний g(E).
Рис. 4. Схематическое изображение плотности электронных состояний и электронные переходы между локализованными состояниями при сравнительно высоких (1) и низких (2) температурах
а
у, Ом 1 -м 1
10-
10-
10-
J_I_I_I_
0,26
0,28
б
0,30 г1/4, к-1/4
Рис. 5. а — температурная зависимость статической проводимости оксидных пленок тантала (1) и ниобия (2); б — температурная зависимость проводимости оксидных пленок тантала (1) и ниобия (2) в области сравнительно низких температур
Предлагаемая методика была использована применительно к аморфному оксиду ниобия, где для изменения положения уровня Ферми использовалась восстановительная термообработка. Последняя, согласно данным рентгеновской фотоэлектронной спектроскопии, повышает концентрацию заполненных электронами локализованных состояний, отвечающих катионам основы в непредельной степени окисления [12]. Полученные результаты указывают на наличие в квазизапрещенной зоне аморфного оксида ниобия максимума плотности электронных локализованных состояний шириной приблизительно 300 мэВ. На основании этих данных можно объяснить наблюдаемые на эксперименте различия в проводимости и энергии активации электропроводности аморфного оксида ниобия и родственного ему оксида тантала.
Для определения подвижности носителей заряда в неупорядоченных системах предлагается методика, состоящая в измерении частотной зависимости диэлектрических потерь образца после неоднородного по его объему оптического возбуждения в области фундаментального поглощения [13, 14].
В основе развиваемой методики лежит следующая модель. Плоский образец диэлектрика толщиной й подвергается воздействию облучения, создающего избыточные носители заряда вблизи одной из поверхностей х = 0; толщина обогащенного фотоиндуцированного носителями заряда слоя й1 < й. Далее облучение прекращается, и образовавшийся электронный пакет, после захвата электронов на локализованные состояния, за счет диффузии или дрейфа начинает распространяться по направлению к противоположной поверхности образца х = й. Контакт диэлектрика с электродом на облучаемой поверхности считается для электронов блокирующим, так что носители заряда не могут покинуть образец.
Определяются активная и реактивная составляющие проводимости как функ-
ции частоты переменного поля = у' и Imу(<я) = у" в различные моменты времени после облучения и производная от них величина тангенса угла диэлектрических потерь
tgS(co, t) = Y ( , ) . Учитывая предпосыл-
Y' Кt)
ки модели, можно полагать, что эволюция частотной зависимости tgS(c$) со временем зависит от диффузионного расплыва-ния пакета и несет информацию о подвижности носителей заряда.
Дальнейший расчет определяется конкретной зависимостью у' и у' от координаты х. Рассмотрим простейшую двухслойную модель, полагая, что при 0 < х < d1, у' = у' и у" = у' , а при d1 < х < d у' = у'2 и у" = у''. Параметры di, у'2 и у'2 зависят,
разумеется, о времени наблюдения t. Вопрос о зависимости результатов наблюдения от формы пакета обсуждается ниже.
Выражение для tgS в рассматриваемой модели двухслойного диэлектрика имеет вид:
tgS =
dj [ )2 +(у'')' ]+ г d 2у 2 [ )2 ')21 , (5)
dlУl"[(у2)' ')2 ]+ d у [ 22
где й 2 = й — й1 — толщина свободного от избыточных носителей заряда слоя диэлектрика.
Будем полагать, что проводимость обогащенной неравновесными носителями заряда области от частоты не зависит и равна статической проводимости. Зависимость же у2 (а) должна быть существенной, но поскольку в свободной от фо-тоиндуцированных носителей заряда у'2 <<у[, величину у[ в конце расчета можно устремить к нулю, и для простоты мы временно положим ее также независящей от частоты. Реактивные составляющие проводимостей слоев примем равными их поляризационным составляющим.
В этом приближении частотная зависимость tgS(cд) принимает вид:
(6)
а2 + Ь/2 /а + ь/2'
где а = 2п[лЛ (г')2 + л2£2 (г!)2 ]; Ь = 8П [й^^2 + й ];
а2 = ^)2 + ¿2/2 Ш; Ь2 =П + 4)-
Здесь / = , е\ и £2 — абсолютные диэлектрические проницаемости слоев диэлектрика, которые в силу фотодиэлектрического эффекта могут быть различны.
На рис. 6 приведен график зависимости tgS(/), соответствующий (6). Такая зависимость с минимумом и максимумом имеет место в области параметров задачи,
когда — > 9— . Если это неравенство не
Ь1 Ь2
выполняется, tgS с ростом частоты монотонно падает.
Характерные точки на кривой имеют следующие координаты:
Ь,
/ =
^ тах
/т =
л тт
Х
tgS . =
о тт
а1 + Ь1 а2Ь2
tgS(cд)
(7)
(8)
Наиболее важным в плане объяснения экспериментальной ситуации является случай, когда проводимость обогащенной неравновесными носителями заряда области намного выше, чем проводимость области, свободной от них, а толщина области обогащения мала по сравнению с полной толщиной образца, т. е.
у[ >>/'; Л <<й; ¿у >>йу2.
В этом случае (7) и (8) приобретают вид:
/ =
^ тах
2пЕ,
й1£2 2йв1
/т =
^ тт
2пб2
гуЛ
I
= 2
Г'Л 2
(9)
(10)
а
тах
1
г . г
'тт Атах
Рис. 6. Теоретическая частотная зависимость тангенса угла диэлектрических потерь образца, подвергнутого неоднородному оптическому возбуждению
Таким образом, значение /тах несет информацию о проводимости обогащенной фотоиндуцированными носителями заряда области, а величина tgSmax — о ширине этой области, т. е. о кинетических параметрах диэлектрика. Что же касается /тт и tg8ш_п, то эти величины содержат информацию о проводимости свободного от неравновесных носителей заряда слоя.
Важным обстоятельством является тот факт, что величина tgSmax не зависит от у'. Это позволяет извлекать из временной зависимости tgSmax информацию об изменении ширины электронного пакета в образце диэлектрика без искажающего влияния временной зависимости проводимости обогащенной области (при условии, что у[ >> у').
Обсудим теперь влияние на результаты формы пакета. Для простоты предположим, что реактивная составляющая проводимости у" не зависит от х, а активная составляющая у'(х) ^ 0 при х ^ да , что, в терминах двухслойной модели, соответствует е1 = / и у'г = 0 .
Выражение для tgS при этом имеет форму:
.1Г_Г
" ИГ [ ' (х
' (х)у "(Их
(11)
^ Ш1+И2'
Примем, что функция у'(х) имеет
вид у'(х2 = у0/^У( ^, где функция
/ (г) (г = х/а^) изменяется от 1 при г = 0 до 0 при г = 1.
Тогда (11) преобразуется к виду:
И,
= —а>т
/ (г
И [/(г)] +М2'
(12)
где т = /
Анализ частотной зависимости обнаруживает наличие максимума, причем
&^тах = (А =
а,
где
да
А = с [
/2 + с2
(13)
(14)
а с — численный коэффициент, определяемый из уравнения
'■/(/2 - с 2(
= 0.
(15)
0 (2 + с2 )2
Таким образом, если известна форма пакета /^^^, используя (14) и (15),
можно найти коэффициенты А и с. В рассматриваемой в настоящем разделе модели двухслойного диэлектрика, когда /(г) = 1 при 0 < г < 1 и /(г) = 0 при г > 1, коэффициент с равен 1, а коэффициент А = >/2.
На основании представленной модели может быть предложена следующая методика определения подвижности носителей заряда. Определяется зависимость величины tg 2£тах от времени, прошедшего после возникновения электронного пакета. На графике этой зависимости выделяется линейный участок (на больших временах), по наклону которого определяется подвижность носителей заряда.
Эта методика была использована для определения подвижности носителей заряда в аморфном оксиде тантала.
Как видно из рис. 7, частотная зависимость tgS оксида тантала, подвергнутого воздействию облучения, обнаруживает характерный для неоднородного диэлектрика максимум, положение и величина которого зависят от времени, прошедшего после прекращения облучения. На временах, превышающих 30 минут после воздействия излучения, временная зависимость tg2c>max выходит на линейный участок (рис. 8), что говорит о достижении стационарного значения подвижности носителей заряда. Оценка подвижности носи-
телей заряда в аморфном оксиде тантала дает величину и = 10— см /В с.
tg5102
Рис. 7. Частотная зависимость тангенса угла
диэлектрических потерь оксида тантала в исходном состоянии (1) и подвергнутого воздействию ультрафиолетового излучения в различные моменты наблюдения после прекращения облучения: 30 мин (2); 17 ч (3); 22 ч (4)
Рис. 8. Зависимость квадрата величины
максимума тангенса угла диэлектрических потерь оксида тантала, подвергнутого неоднородному оптическому возбуждению, от времени, прошедшего после его воздействия
Эксперимент показал, что температурная зависимость подвижности имеет активационный характер, свидетельствующий в пользу наличия в аморфных металлооксидах сильной электрон-фононной связи и существенного вклада в энергию активации электропроводности энергии межузельных поляронных перескоков.
Предлагаемая экспериментальная методика определения малой подвижности носителей заряда может оказаться полезной для анализа кинетических свойств широкого класса материалов, возможности определения подвижности которых ограничены. Так, в материалах с ¿и < 1 ш 2 / В резко отличаются дрейфовая и холловская подвижность [6], в результате чего традиционные измерения проводимости и коэффициента Холла не дают возможности получить информацию о и . Единственным методом измерения и в материалах с малой подвижностью является дисперсионный перенос. Однако его особенности в случае прыжкового механизма электропереноса затрудняют однозначную интерпретацию экспериментальных данных и оценку времени носителями заряда образца. Кроме того, время пролета для материалов с крайне низкой подвижно -стью может оказаться весьма большим, в то время как предлагаемая методика позволяет определять и на временах, когда электронный пакет еще не покидает образца.
БИБЛИОГРАФИЧЕСКИЕ ССЫЛКИ
1. Электронная теория неупорядоченных полупроводников / В. Л. Бонч-Бруевич, И. П. Звягин, Р. Кайпер, А. Г. Миронов, Р. Эндерлайн, Б. Эссер. — М., 1981.
2. Лившиц И. М., Гредескул С. А., Пастур Л. А. Введение в теорию неупорядоченных систем. — М., 1982.
3. ЗайманДж. Модели беспорядка. — М., 1982.
4. Шкловский Б. И., Эфрос А. Л. Электронные свойства легированных полупроводников. — М., 1979.
5. Звягин И. П. Кинетические явления в неупорядоченных полупроводниках. — М., 1984.
6. HotterH., Bryksin V. V. Hopping Conduction in Solids. — VCH, 1985.
7. Nonlinear current-voltage characteristics of Ta2O5 and Nb2O5 amorphous oxides / V. V. Bryksin,
A. V. Goltsev, S. D. Khanin et al. // Phys. Stat. Sol. (b). 1990. Vol. 161. № 2. P. 777-781.
8. Брыксин В. В., Карпухина Л. Г., Ханин С. Д. Частотная зависимость проводимости аморфных окислов тантала при наличии постоянного смещающего напряжения // Физика тв. тела. 1990. Т. 32. № 12. С. 3564-3570.
9. Брыксин В. В. Проводимость на переменном сигнале неупорядоченных систем в токовом состоянии // Физика тв. тела. 1990. Т. 32. № 9. С. 2570-2578.
10. Брыксин В. В., Дьяконов М. Н., Ханин С. Д. Прыжковый перенос в некристаллическом окисле тантала // Физика тв. тела. 1980. Т. 22. № 5. C. 1403-1410.
11. Электронный перенос в аморфных окислах переходных металлов / В. В. Брыксин, М. Н. Дьяконов, В. М. Муждаба, С. Д. Ханин // Физические явления в некристаллических полупроводниках: Материалы Международной конф. «Аморфные полупроводники—80». — Кишинев, 1980. С. 34-38.
12. Исследование электронных состояний окислов тантала и ниобия методом электронной спектроскопии / М. Н. Дьяконов, Ю. П. Коститков, Д. В. Костров, С. Д. Ханин // Изв. АН СССР. Неорг. матер. 1980. Т. 16. № 5. С. 881-884.
13. Эффект фотопамяти в анодных оксидных пленках / М. Н. Дьяконов, Л. Г. Карпухина,
B. М. Муждаба, С. Д. Ханин // Письма в ЖТФ. 1982. Т. 8. № 19. С. 1186-1189.
14. Bryksin V. V., Goltsev A. V., Khanin S. D. Relation between the tangent of the angle of dielectric losses and low mobility in dielectrics // Phil. Mag. (b). Vol. 64. № 1. P. 91-100.
S. Khanin
MODELS OF TRANSFER PHENOMENA AND METHODS OF THE ANALYSIS OF DISORDERED SYSTEMS PROPERTIES WITH CHARGE CARRIERS LOCALIZATION
In this paper the models of non-stationary and non-equilibrium electron phenomena in disordered dielectrics are developed. The possibilities of experimental identification of localization mechanisms and low mobility of charge carriers are presented.