ФИЗИКА КОНДЕНСИРОВАННОГО СОСТОЯНИЯ ВЕЩЕСТВА
Особенности частотной зависимости бесфононной прыжковой проводимости
М. А. Ормонт
Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра физики полупроводников. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
E-mail: [email protected]
Статья поступила 12.04.2012, подписана в печать 29.05.2012.
В рамках теории возмущений проведен расчет мнимой части бесфононной проводимости слаболегированного компенсированного полупроводника. Показано, что при использовании базиса локализованных функций атомного типа суперлинейной частотной зависимости вещественной части проводимости соответствует частотная зависимость мнимой части проводимости близкая к линейной. Найдено, что при частотах меньших частоты перехода (кроссовера) шсг от линейной к квадратичной частотной зависимости вещественной части проводимости, тангенс угла диэлектрических потерь слабо зависит от частоты и определяется отношением /гшсг к ширине примесной зоны. Показано, что из измерений тангенса угла диэлектрических потерь можно получить информацию о радиусе локализации примесных состояний.
Ключевые слова: бесфононная (резонансная) прыжковая проводимость, тангенс угла диэлектрических потерь.
УДК: 621.315.592. PACS: 72.20.Ее.
Введение
Известно, что из частотных зависимостей проводимости можно получить информацию об особенностях механизма переноса носителей заряда в среде. Приложенное к системе переменное электрическое поле напряженности Е = Е0 ехр(—вызывает смещенный по фазе отклик системы — например, плотность тока / = /оехр(—— у)), которая равна / = сгЕ; а = (71 + /(72 — комплексная проводимость системы, — угол сдвига фазы между током и напряжением, = (72/(71. Величину мощности диэлектрических потерь Р часто выражают через угол 7, дополняющий до 7г/2 угол сдвига фаз между током и напряжением Р = £¿(71/2 = (£|ег2)/2^7, где 7 = тг/2 - ц> — угол диэлектрических потерь, tg 7 = (71/(72 [1].
Исследования проводимости неупорядоченных систем, в частности слаболегированных компенсированных полупроводников, на изоляторной стороне перехода металл-диэлектрик дают универсальную степенную частотную зависимость вещественной части проводимости
СТ1 = Йе {о) = Аш*, (1)
где А, я — постоянные; как правило, 0 < я ^ 1 [2-6]. Степенная частотная зависимость проводимости (1) обычно свидетельствует о прыжковом механизме электронного переноса [6]; так, для случая прыжкового переноса по локализованным состояниям примесной зоны слаболегированных компенсированных полупроводников теория предсказывает частотную зависимость вещественной части проводимости типа (1) [2-7]. Степенная зависимость (1) с обычно связывается с прыжковой проводимостью по локализованным состояниям электронов на примесных центрах с участием фононов [2]. Аналогичная частотная зависимость
с получается при низких частотах и в случае
низкотемпературной бесфононной (резонансной) прыжковой проводимости при учете кулоновских корреляций локализованных носителей [4, 5, 7]. С ростом частоты теория бесфононной проводимости предсказывает переход (кроссовер) от линейной частотной зависимости проводимости (с к зависимости, близкой к квадратичной (я и 2) [3-5, 7]. Подобный переход наблюдался при возрастании частоты в области около 1 ТГц в легированном кремнии в окрестности перехода металл-изолятор [8, 9] и в металлических нанокомпо-зитах [10].
Анализ частотных зависимостей мнимой части проводимости неупорядоченных систем в изоляторном состоянии (проводимость на постоянном токе (7^ = 0) часто основывается на соотношениях Крамерса-Кро-нига [11], согласно которым при 0 < 5 < 1 степенной частотной зависимости вещественной части проводимости (1) соответствует такая же частотная зависимость ее мнимой части
<72 = 1т(<7) = -<711д(у). (2)
Из равенства (2) следует, что в этом случае угол сдвига фазы ц> между током и напряжением (угол диэлектрических потерь 7) не зависит от частоты
* V = -У = - = -(у) * (3)
приближенное равенство выполняется при значениях я, близких к единице [6]. Тогда \&2\/<у\ 1, т.е. угол сдвига фазы |<р|, мало отличается от 7г/2 (7«0).
Будем рассматривать случай низких температур, при которых проводимость слаболегированного компенсированного полупроводника по примесной зоне определяется бесфононной составляющей. Не останавливаясь
на деталях поведения вещественной части бесфононной проводимости, которые были весьма подробно обсуждены в работах [7, 12, 13], отметим наблюдавшуюся в широкой области частот суперлинейность (я > 1) частотной зависимости вещественной части проводимости ел (ш) [14, 15]. В этом случае соотношение (2) заведомо неприменимо и мнимую часть бесфононной проводимости <72 (о>) следует непосредственно рассчитывать по теории возмущений.
Измерения частотной зависимости низкотемпературной проводимости легированного компенсированного полупроводника на изоляторной стороне перехода металл-диэлектрик, проведенные в [16], показали, что при частотах меньших частоты кроссовера компоненты проводимости имеют близкие частотные зависимости |<72(о;)| ~ <71 (ш) ~ а/ с показателем степени я и 1. Однако величина их отношения \cr2\cr1 ~ 102 [16] превосходит значение, даваемое теорией релаксационной проводимости, более чем в 30 раз [17].
Принимая во внимание суперлинейность (я > 1) частотной зависимости вещественной части проводимости <71 (ш) и расхождение в значениях тангенса угла диэлектрических потерь ^gJ = o■l/a2, получаемых из теории и эксперимента, вычислим в рамках теории возмущений мнимую часть бесфононной проводимости <72 (ш) и отношение |о"21/о"1 ■
Расчет мнимой части низкотемпературной бесфононной проводимости
Мы рассматриваем случай компенсированного полупроводника п-типа, считая, что концентрация случайно расположенных в пространстве примесных центров мала, так что при рассмотрении проводимости на переменном токе можно воспользоваться парным приближением, когда полный ток выражается в виде суммы вкладов парциальных токов от отдельных пар центров локализации [18]. При решении задачи о бесфононной проводимости неупорядоченной системы можно использовать полный ортонормированный набор одноэлектронных функций {ф\}, соответствующих эффективному одночастичному гамильтониану Не, учитывающему случайное поле (А-представление) [18]. Полагая, что уровень Ферми находится ниже уровня протекания, т. е. попадает в область энергий, отвечающую локализованным состояниям, можно использовать усеченный базис из локализованных функций, взятых из набора {ф\} ■ В случае сильной локализации функции ф\ близки к волновым функциям основных состояний примесей; при этом число А отвечает номеру центра. Приближение слабого легирования и соответствующее ему неравенство а < (Л^)1/3 (Л^ — концентрация легирующей примеси, а — радиус локализации состояний) дает возможность пренебречь изменением потенциала, создаваемого всеми другими центрами, на радиусе локализации а волновой функции и позволяет учитывать смещения энергий электронов классическим образом, добавляя к соответствующей энергии электрона слагаемое -в(р(г\); <р(г\) — кулоновский потенциал, создаваемый заряженными центрами в точке расположения центра с номером А. При этом волновая функция примесного состояния остается неизменной и близкой к волновой функции, соответствующей изолированной
примеси. Разброс уровней, возникающий за счет беспорядка в расположении заряженных примесей, порядка (е2/к)Ыуз; при этом вклад в беспорядок в положении уровней за счет квантово-механической гибридизации состояний мал в силу экспоненциальной малости интеграла перекрытия волновых функций локализованных состояний примесей.
Согласно [5, 18], выражение для бесфононной проводимости неупорядоченной системы имеет вид
а(ш) =
1е2ш
~ЙГ
V Шп г)]/)]2 (^ЬМе/))
т, ¥1
где /, / — номера центров локализации; п — единичный вектор параллельный внешнему электрическому полю; пр(е) — средние числа заполнения состояний с энергией е;
е2
АЩ = Е] - е;--= <р/ - (5)
кг1{
— изменение энергии неупорядоченной системы при одноэлектронном переходе из начального состояния / в конечное /; Гц — расстояние между центрами / и /; к — диэлектрическая проницаемость среды; е — заряд электрона; ^ — энергии электрона в состояниях /,/; ег, е/ — самосогласованные энергии электрона, отвечающие состояниям /, / [19]. Выражение (4) получено с использованием адиабатической гипотезы; соответственно а — малая положительная величина, отвечающая адиабатически медленно возрастающему электрическому полю Е = Ео ехр(—/о;? + (а/й)0. Поскольку выражение (Нш - АНц + /а)^1 находится в (4) под знаком суммы (интеграла), то, выполнив согласно стандартной процедуре предельный переход Нта_»о> приходим к равенству
(Пш - АЩ + ¿а) = ~ Ы{ПШ ~ АН^> (6)
где Р — интеграл в смысле главного значения.
Соответственно компоненты проводимости равны
же2ш
0-1 м =-у— |<г|(л,г)|/)|2(Ме;)-Ме/))х
х 6(АЩ - Пш),
(7а)
е2ш ^ 2 (М^) - Ме/)) /7кЧ
<72М = —- |{г|(я,г)|/)| —7Т~.-Л и Ч ; (7Ь)
И)
(Нш-АЩ)
т, ¥1
согласно (6), при переходе к интегрированию в (7Ь) следует рассматривать интеграл в смысле главного значения.
Выражение для матричного элемента при Гц > а имеет приближенный вид [7]
(/|(п,г)|/) = М|(п,г)№/> » ехр СО5 0, (8)
где в — угол между векторами п и Гц. Матричный элемент имеет выраженный максимум при Гц = За; соответственно основной вклад в мнимую часть проводимости (7Ь) вносят слагаемые с Гц ~ гор\ > а и (7Ь) можно записать в виде
Гц
<72 М :
е2ш
~ЙГ
2М(«.г)1/>г-. Л ■ (9)
Перейдем в (9) от суммирования по парам к интегрированию по энергиям и пространственным координатам центров:
Í2
ст2 М;
2 2
dei de¡ йгц \(i\(n,r)\f}\2 х
(rtF(e¿) - nF(sf))
(e¿ - e¡ + huí)
(10)
de¡ def
(nF(£i) -nF{e¡))
(e¿ -£f + ñwi)
ñui/2 ¡í+huji-x
dy-
Í2
¡1+X
dx
dx
dyx =
íi
¡i+x
ñw¡/2 )i+ñuj i—x Í2
dx
{huj\ - 2x)
(11)
где ^ = (А-В+ЬшО/2, £,2 = (В^А + Пш0/2 — границы области интегрирования х е [£1, £2] ■ Будем рассматривать частоты, для которых выполняется неравенство Нш[ = е2/кгорх + Нш < 5 = В-А, т.е. ^<0, £,2>0, 161 <^2- Тогда интеграл (11) равен
dx
(fiuj 1 — 2х)
= hco¡
dx x
í 2
— \-2S, (12)
líil
где \£l\ = (B-A-ñwl)/2.
líil
где ñüoi = e2/nropt + пш. С учетом того что в принятой модели неупорядоченной среды ширина примесной зоны 5 существенно больше ширины куло-новской щели, плотность состояний можно считать постоянной: р(е) и ро при е е [Л,В], р(е) = 0 при eÍ[A,B], 5 = В-А [4, 19]. Напомним, что кулонов-ская щель, возникающая в одночастичной плотности состояний р(е), описывающей распределение самосогласованных энергий взаимодействующих локализованных носителей заряда в основном состоянии системы, является следствием дальнодействующего характера кулоновского взаимодействия [19]. Согласно [7], оптимальная длина прыжков, вносящих основной вклад в проводимость, соответствует переходам между состояниями с энергиями, лежащими вне кулоновской щели. Соответственно в случае столообразной модели плотности состояний кулоновские эффекты, приводящие к появлению кулоновской щели, играют малую роль.
В пределе низких температур в качестве начального состояния выберем основное состояние системы. В этом случае все состояния с энергиями e¡ < р заняты, а с энергиями e¡ > р — свободны (р — уровень Ферми) [19].
Перейдем к новым переменным х = (e¿ — e¡ + Ьш\)/2, у = (e¿ +e¡ + hüjl)/2; тогда e¿ = х + у - huj\, £f=y - х.
Разность чисел заполнения отлична от нуля и равна
nF{e{) -nF(ef) = 1
при условии £i< Р, £f > Р, Т.е. р + х < у < p + ñüüi - х, х < ñüüi/2;
nF(£i)^nF(£f) = ^l
при условии £i > Р, £f < Р, Т.е. p + ñüJ[ - X < у < р + X , Ншi/2 < х.
Тогда выражение для интегралов по энергии в (10) можно записать в виде
Интеграл (12) имеет смысл, если понимать под J ^f
ii
li.l
интеграл в смысле главного значения Р J ^f = 0; тогда
ii
(12) принимает вид
«2
dx-
= Нш 1 In
íi
5 — huj 1
25 «^25; (13)
приближенное равенство соответствует области частот <с 5.
Подставляя (8) в (10), получаем, что интеграл по пространственным координатам центров в (10) равен
drif\{i\(n,r)\f)\2 = j-ya5C3,
(14)
где С3 = Г(9).
Из вида подынтегральной функции в выражении (14) следует, что максимум подынтегрального выражения достигается при значении Гц = 4а и основной вклад в мнимую часть проводимости (7Ь) вносят слагаемые с Г1( ~ гор\ ~ 4а.
С учетом равенств (13) и (14) выражение для мнимой части проводимости (10) принимает вид
'5 + ПшЛ 25У (15)
(72М и е2шр20^^а5С3 {huji In
5 — ^
В области частот, для которых выполняется неравенство Нш\ = е2/кгор\ + Нш <с 5, выражение для мнимой части равно
а2(ш);
2 2 Va>pz0
3 • 27
а5С3 25.
(16)
Воспользуемся найденным в работе [7] (без учета гибридизации волновых функций удаленных центров) выражением для вещественной части проводимости
М = е2р2аьш (17)
ка
где С[ = С2 = Г(8), С3 = Г(9), | = 0.125 -численные коэффициенты. В этом случае оптимальная длина прыжка определяется параметрами системы и не зависит от частоты, поскольку с уменьшением расстояния между центрами в паре уменьшается и изменение дипольного момента системы при электронном переходе, а с увеличением расстояния между центрами происходит экспоненциальное уменьшение перекрытия волновых функций состояний, отвечающих центрам локализации [7]. Согласно (17), переход от линейной к квадратичной частотной зависимости вещественной части бесфононной проводимости происходит в окрестности частоты кроссовера шсг, определяемой равен-
С 2
ством энергий Ншсг = 2^ ^ [7]. Соответственно при
частотах, меньших частоты кроссовера шсг = ^ ^, имеем
a-i (ш):
о ка
(18)
Отметим, что при частотах ш < шсг неравенство + Ьш < 5 выполняется, поскольку кшсг <с 5. Соответственно отношение компонент проводимости |ст2|/с1 в области частот ш < шсг равно
(71 ТГ ПШСГ
т. е. определяется отношением ширины примесной зоны 6 к характерной кулоновской энергии йшсг = ^ ^.
Если для оценки учесть, что разброс уровней, возникающий за счет беспорядка в расположении заряженных примесей порядка (е2/к)Ы^3, то (19) принимает
вид 1 сх21 /га 102Ы^3а; т.е. из отношения компонент проводимости (тангенса угла диэлектрических потерь) можно получить информацию о радиусе локализации примесных состояний. Для приведенных в работе [16] значений концентрации примеси Л^ ~ 1018 см^3 и радиусов локализации а~ 100 А имеем \&2\/<у\ ~ Ю2; это согласуется с экспериментально полученными значениями \(Т2\/&1 [16]. Следует, однако, иметь в виду, что в этом случае при (Л^)^1/3 ~ а мы выходим за рамки применимости теории (а < (Л^)^1/3).
Заключение
Таким образом, в рамках теории возмущений показано, что при электронном прыжковом транспорте по локализованным состояниям примесной зоны мнимая и вещественная части проводимости степенным образом зависят от частоты; причем суперлинейной частотной зависимости вещественной части проводимости соответствует частотная зависимость мнимой части проводимости, близкая к линейной.
Из соотношений (16) и (18) следует, что в области ш<шсх компоненты проводимости а\(ш) и <72 (ш) имеют близкие частотные зависимости с показателем степени 5 га 1, т.е. |ст2(а;)| ~ <71 (ш) ~ и/ (я«1). Полученные выражения (16) и (18) для компонент комплексной
проводимости а(ш) согласуются с экспериментальными частотными зависимостями а\(ш) и ст2(ш) как по показателю степени s га 1, так и по величине отношения |<7г|/<71 [16].
Список литературы
1. Фрелих Г. Теория диэлектриков. Диэлектрическая проницаемость и диэлектрические потери. М., 1960.
2. Pollak AL, Geballe Т. Н. // Phys. Rev. 1961. 22. P. 1742.
3. MottN.F. 11 Phil. Mag. 1970. 22. P. 7.
4. Шкловский Б.И., Эфрос А.Л. // ЖЭТФ 1981. 81. С. 406.
5. Zvyagin LP. II Charge Transport in Disordered Solids with Applications in Electronics / Ed. by S. Baranovski. Chichester, 2006. P. 339.
6. Звягин И.П., Кинетические явления в неупорядоченных полупроводниках. М., 1984.
7. Ормонгп М.А. II Вестн. Моск. ун-та. Физ. Астрон. 2011. № 2. С. 57.
8. Lee AL, Stutzmann M.L. // Phys. Rev. Lett. 2001. 87. P. 056402.
9. Helgren E., Armitage N.P., Gruner G. // Phys. Rev. B. 2004. 69. P. 014201.
10. Reedijk J. A., Adriaanse L.J., Brom H.B. et al. 11 Phys. Rev. B. 1998. 57. P. R15116.
11. KronigR. de L. 11 J. Opt. Soc. Am. 1926. 12. P. 547.
12. Звягин И.П., Ормонгп М.А. 11 Вестн. Моск. ун-та. Физ. Астрон. 2008. № 4. С. 44.
13. Ормонгп М.А. II Вестн. Моск. ун-та. Физ. Астрон. 2010. № 2. С. 46.
14. Hering AL, Scheffler AL, Dressel M. et al. 11 Phys. Rev. B. 2007. 75. P. 205203.
15. Ritz E., Dressel M. 11 Phys. Stat. Sol. 2008. 5. P. 703.
16. Helgren E., Armitage N.P., Gruner G. // Phys. Rev. Lett. 2002. 89. P. 246601.
17. Эфрос А.Л. 11 ЖЭТФ. 1985. 89. C. 1834.
18. Бонч-Бруевич В.Л., Звягин И.П., Кайпер Р. и др. Электронная теория неупорядоченных полупроводников. М., 1981.
19. Efros A.L., Shklovskii B.I. // Electron-Electron Interactions in Disordered Systems / Ed. by A. L. Efros, M. Pollak. Amsterdam, 1985. P. 409.
Features of the frequency dependence of phononless hopping conductivity M. A. Ormont
Department of Semiconductor Physics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia. E-mail: [email protected].
Within the framework of perturbation theory, the imaginary part of the phononless conductivity of lightly doped compensated semiconductor has been calculated. It is shown that the use of localized atomic-like function basis leads to superlinear frequency dependence of the real part of conductivity that corresponds to approximately linear frequency dependence of the imaginary part of conductivity. It is found that at frequencies below the frequency of the transition (crossover) wcr from linear to quadratic frequency dependence of the real part of conductivity dielectric loss tangent is weakly dependent on the frequency and is determined by the ratio of hui„ to the width of the impurity band. It is shown that information on the radius of localization of impurity states can be obtained from the measurements of dielectric loss tangent.
Keywords: phononless (resonance) hopping conductivity, dielectric loss tangent. PACS: 72.20.Ее. Received 12 April 2012.
English version: Moscow University Physics Bulletin 5(2012).
Сведения об авторе
Ормонт Михаил Александрович — канд. физ.-мат. наук, доцент, доцент; тел.: (495) 939-41-18, e-mail: [email protected].