CC BY
37
constrained fields in the elements of constructions. Kyiv: 6. Pugachev, S. I. (Ed.) (1984). Pezokeramic
Scientific thinking, 280. transformers. Lviv: Shipbuilding, 256.
5. Ditkin, V. А., Prudnikov, А. О. (1965). Reference 7. Slepan, L. I. (1972). Non-stationary resilient waves.
Book on an operating calculation. Moscow: Higher school, 466. Lviv: Shipbuilding, 374.
Рекомендовано до публгкацИ доктортехтчних наук, професор Савгн В. Г.
Дата надходженнярукопису 27.01.2015
Бабаев Александр Арташесович, кандидат физико-математических наук, доцент кафедры теоретической механики, Национальный технический университет Украины "Киевский политехнический институт", пр. Победы, 37, г. Киев, Украина, 03056 E-mail: [email protected]
Кришталь Владимир Федорович, кандидат технических наук, доцент кафедры теоретической механики, Национальный технический университет Украины "Киевский политехнический институт", пр. Победы, 37, г. Киев, Украина, 03056 E-mail: [email protected]
Можаровская Тамара Николаевна, кандидат технических наук, доцент кафедры теоретической механики, Национальный технический университет Украины "Киевский политехнический институт", пр. Победы, 37, г. Киев, Украина, 03056 E-mail: [email protected]
УДК 537.32
DOI: 10.15587/2313-8416.2015.36332
ПЕРЕНОС ТЕПЛА ФОНОНАМИ В ТРАНСПОРТНОЙ МОДЕЛИ ЛАНДАУЭРА-ДАТТЫ-ЛУНДСТРОМА
© Ю. А. Кругляк
С позиций транспортной модели ЛДЛ строится обобщенная модель переноса тепла фононами. Аналогично фермиевскому окну электронной проводимости вводится понятие фермиевского окна фононной проводимости и через него выводится общее выражение для решеточной теплопроводности, в котором с самого начала фигурирует квант теплопроводности. Подчеркивается подобие и различия в построении теории электронной проводимости и теории теплопроводности
Ключевые слова: нанофизика, наноэлектроника, фононный транспорт, квант теплопроводности, коэффициент прохождения, фононные моды, дебаевская модель, рассеяние фононов
On the basis of Landauer-Datta-Lundstrom transport model the generalized model of heat transfer by phonons is formulated. Similarly to the Fermi window for electron conductivity the concept of the Fermi window for phonon conductivity is introduced and used to obtain the general expression for the lattice thermal conductivity with the quantum of thermal conductivity appearing at the very beginning. There are emphasized the similarity and differences in the construction of the theory of electron conductivity and the theory of heat conduction Keywords: nanophysics, nanoelectronics, phonon transport, thermal conductivity quantum, transmission coefficient, phonon modes, Debye model, phonon scattering
1. Введение
При транспорте электронов переносится как заряд, так и тепло. В металлах почти все тепло переносится электронами, а в полупроводниках только меньшая часть тепла переносится электронами, и почти все тепло переносится фононами.
Поток тепла, вызванный транспортом фононов вдоль проводника (ось х)
К , [Вт/м2] (1)
где коэффициентом пропорциональности потока тепла градиенту температуры является удельная решеточная теплопроводность к1. У алмаза, одного из лучших проводников тепла, кь « 2000 Вт/мК, у такого плохого проводника тепла как стекло кь «1
Вт/мК, а у воздуха кь «0.025 Вт/мК. Если электрическая проводимость твердых материалов меняется на 20 и больше порядков, то их решеточная теплопроводность изменяется всего на 3-4 порядка. Тем не менее, есть много общего в описании переноса тепла электронами и фононами в рамках транспортной модели Ландауэра - Датты -Лундстрома (ЛДЛ) [1-4]. Далее упомянем лишь необходимые нам сведения о фононах.
2. Обзор элементарной физики фононов
Физика фононов подробно изложена в [5-8]. Характерный ход дисперсии Е(к) для электронов качественно показан на рис. 1, а дисперсии Псо(д) для фононов - на рис. 2.
= vsq>
Vs =4Klas 1
m
(2)
f
Рис. 1. Упрощенная зонная структура. Ширина зоны проводимости помечена BW/Band Width. Часто используемая параболическая дисперсия показана пунктиром
Кристаллическая решетка периодична в реальном пространстве, а дисперсия - в £-про-странстве. Аналогично для фононов: только в ^-пространстве. В обоих случаях решения реализуются в зоне Бриллюэна и ширины обоих зон конечны.
Рис. 2. Упрощенный ход дисперсии для фононов. Пунктиром показаны простейшие модели Эйнштейна и Дебая
Подробнее ход дисперсии для фононов показан на рис. 3.
Для кристаллических 3D проводников при распространении колебаний по решетке возможны три состояния поляризации: одно продольное, в направлении распространения колебаний (longitu-dinal/L) и два поперечных (transverse/T). В области низких энергий фононные моды называют акустическими/A. Вблизи нулевых волновых чисел (q ^ 0) у акустических мод дисперсия линейная. Для продольной акустической моды LA
где скорость звука V определяется упругой постоянной кеЫ и массой те#. Характерные скорости
звука порядка 5 х 103 м/с, что примерно в 20 раз меньше скорости электронов. Из уравнений (2) следует, что у веществ с тяжелыми атомами скорость звука меньше и соответственно меньше ширина фононной зоны дисперсии по сравнению с веществами из легких атомов.
Рис. 3. Дисперсионные кривые для фононов
На рис. 3 показаны также три оптические моды. Дисперсия этих мод незначительна: частота оптических мод Ь и Т сравнительно мало зависит от д. Различие между акустическими и оптическими модами состоит в том, что в области малых д~0 соседние атомы смещаются в одном и том же направлении для акустических мод и в противоположных направлениях для оптических мод. В полярных веществах свет взаимодействует с оптическими модами; отсюда и их название. Скорость фононов определяется наклоном зависимости частоты от смещения а(д): акустические
моды характеризуются относительной большей скоростью фононов по сравнению с оптическими модами. Поэтому можно ожидать, что именно акустические моды переносят большую часть тепла.
В качестве реалистического примера приводим фононный спектр графена (рис. 4), построенный на основе экспериментальных данных, хорошо согласующихся с результатами расчетов методами теории функционала плотности [9-11].
В элементарной ячейке графена два атома углерода, что обусловливает наличие в спектре дисперсии графена шести фононных ветвей - трех оптических TO, ZO) и трех акустических (LA, TA, ZA). Последние связаны со смещением всей ячейки из положения равновесия, а оптическим модам отвечают сдвиги атомов в самой элементарной ячейке с сохранением положения центра масс. Фононные моды LA и TA соответствуют продольным и поперечным колебаниям атомов
углерода в плоскости графена, а мода ZA (flexural mode) - колебаниям атомов углерода в направлении, перпендикулярном направлению колебаний мод LA и TA, с выходом атомов из плоскости графена. Моды LA и TA характеризуются линейной дисперсией; соответствующие этим модам скорости звука составляют 21300 и 13600 м/c, соответственно. Все теплофизические свойства графена (тепло про-водность, теплоемкость и др.) определяются фононами.
Eph
см 1
/
Г/о
Ал /г
/ s1* /
' // yS lA
Г
Г
м к
Рис. 4. Фононный спектр графена, опирающийся на совокупные экспериментальные и расчетные данные
энергию электронов Е заменим на энергию фононов Ьсо. В случае электронов переносятся заряды с/. а в случае фононов - кванты энергии Ьсо : заменим д на Ьсо и внесем Ьсо под знак интеграла. Двойка в уравнении (3) для тока электронов соответствует вырождению по спину. Исключим ее из уравнения для фононов, а число состояний поляризации фононов включим в число фононных модМрЬ{Ьсо).
Окончательно для тепла, переносимого фононами, имеем:
Как и в случае электронов, нас интересует режим линейного отклика, или иначе квазиравновесного транспорта. В этом режиме Т ~ Т и П ~ п ~ П ~ П, что позволяет разложить п2 относительно п в ряд Тейлора и ограничиться первым порядком:
дпп
дТ
-AT,
(6)
так что нужная нам разность
3. Обобщенная модель Ландауэра-Датты-Лундстрома переноса тепла фононами
Электрический ток в обобщенной транспортной модели ЛДЛ [3, 4, 12]
I = J Td (E)Md (E) (f - f ) dE [А]
(3)
2д ~Н
пропорционален произведению коэффициента прохождения Те1 (Е) или иначе вероятности прохождения проводника по каналу проводимости энергии Е, числу мод проводимости Ме1 (Е) при этой энергии и разности фермиевских функций на контактах проводника.
Для описания транспорта фононов нам необходимо уравнение, аналогичное (3). Как и в случае электронов, контакты на концах проводника достаточно массивны настолько, чтобы поддерживать тепловое равновесие фононов при разных температурах контактов. Если в случае электронов состояния контактов заселяются в соответствии со статистикой Ферми-Дирака, то в случае фононов - в соответствии со статистикой Бозе-Эйнштейна
Тт
(4)
ё' — 1
Пусть температура на левом и правом контактах будет соответственно Т1 и Т2. Как и в случае электронов, предполагается, что контакты идеальные, другими словами, отражение фононов на контактах не происходит, и коэффициент прохождения действительно соответствует прохождению фононов через весь канал проводимости.
Теперь легко переписать выражение (3) применительно к фононам. Для потока фононов
дп0 дТ
-AT.
(7)
Производную дп0 / дТ находим дифференцированием бозевского распределения (4):
дп 0 д Т
Ьа Т
дпп
дфсо)
где
дпп
1
JuolkT
(8)
(9)
д(Гш) { кТ) (е Теперь уравнение (5) для малой разности температур на контактах можно переписать следующим образом:
Q = ~KL AT,
где решеточная теплопроводность
К, ~\Tph(tico)hipl,(Ьсо)
дп0 \
Ы1 дфсо) J
(10)
ä(na) [ВшП.(П)
Уравнение (10) есть закон теплопроводности Фурье, согласно которому поток тепла антинаправлен градиенту температуры. Полезно обратить внимание на следующую схожесть теплопроводности (11) и электронной проводимости
G =
2q
JTd (EMi (E) dE. (12)
Производная
W (E) = 1-
дE
(13)
известна как фермиевское окно проводимости, выделяющее те моды проводника, которые только и
п
могут внести вклад в электронную проводимость. Окно проводимости нормировано на единицу:
дЕ,
Я-
dE = 1.
(14)
Аналогичную роль окна проводимости фононов играет следующее выражение:
для
3 (ha
Wph(ha)=—\ —
дпп
„ ,- (15)
п' \кТ Д д(Ьа) в котором численный коэффициент необходим для нормировки на единицу выражения в квадратных скобках уравнения (11):
Ш-
- \d(ha)=-дфа)) 3
(16)
так что окончательно решеточная теплопроводность
тг2к2Т
KL =-
3h
-¡Tpli(ha)Mpli(ha)Wph(ha)d(ha) ■ (17)
В этом выражении численный коэффициент £0 = ж2к2Т / 3И « (9.456 х10~13Вт / К2)Т, (18) который мы обсудим позже, есть квант теплопроводности, впервые подтвержденный экспериментально в 2000 году [13].
Сравнивая выражения для решеточной теплопроводности (17) и электронной проводимости (12), опять убеждается в их схожести: и та и другая проводимости пропорциональны соответствующему кванту проводимости, умноженному на интеграл от произведения трех величин - коэффициента прохождения, числа мод и окна проводимости.
Профили окон проводимости для электронов и фононов показаны на рис. 5, 6.
Рис. 5. Окно проводимости для электронов (13) при двух разных температурах. Энергия Е может быть как больше фермиевской, так и меньше ее
Рис. 6. Окно проводимости для фононов (15) при двух разных температурах. Энергия фононов может быть только положительной
Окна проводимости для электронов и фононов не только схожи, но и имеют приблизительно одинаковую ширину в « +2кТ .
4. Теплопроводность массивных проводников
Теплопроводность диффузионных проводников является их ключевым свойством, определяющим работоспособность разнообразных электронных устройств. По аналогии с транспортом электронов (ф-ла (86) из [12]) коэффициент прохождения для фононов
Tph(ha) =
lph(ha) + L
(19)
где Ь - длина ЗЭ проводника, а ХрЬ(Ьа) есть длина
свободного пробега фонона частоты а. Кроме того, число фононных мод, очевидно, пропорционально площади поперечного сечения 3D проводника:
Мр11(Гш)кА. (20)
Умножая и деля уравнение (10) на Л/Ь, получаем
Q = -| klL 1 лЛГ
(21)
откуда следует постулированное в самом начале уравнение (1) для потока тепла:
л Qx L
dT dx
с удельной решеточной теплопроводностью
Kl = K L
L_
(22)
(23)
или после подстановки (19) в (17) в развернутом виде окончательно получаем
тгЬгТ г Л/ , (ha) к, =--lph(ba)Wph(tia)d(ba). (24)
ЗЛ •> А
Поскольку именно окно проводимости фононов Ш^(Ьсо) определяет какие именно
фононные моды могут принимать участие в проводимости тепла, естественно ввести среднее число активных фононных мод на единицу площади поперечного сечения проводника
, > rM Aha)
{MpiJл) = \— .-Wph(ha)d(ha). (25)
так что в компактном виде
ж2 k 2T 3h
{Mph /Л)((Äph)), (26)
где усредненное значение длины свободного пробега фононов
0
KL =
rM , (Ьсо)
-Xph(hm)Wph(hm)d(hm)
(Ы) =--• (27)
,M , Chco)
f \ Wph(hco)d(hco)
А
Итак, паре фононных транспортных уравнений (22) и (26) соответствует схожая по форме пара электронных транспортных уравнений [3, 4, 12]:
J = а d (Ер) xqdx
а =
2q h
M1А){ M».
(28)
(29)
Градиент электрохимического потенциала вызывает поток электронов, а градиент температуры - поток тепла. Термоэлектрические коэффициенты (26) и (29) имеют одинаковую структуру. Это всегда произведение соответствующего кванта прово-димости на число мод проводимости, участвующих в транспорте, и на усредненное значение средней длины свободного пробега соответствующего носи-теля (тока или тепла). К оценке этих трех величин для фононов мы вернемся позже.
5. Удельная теплоемкость
Хорошо известна пропорциональность между удельной теплопроводностью и удельной теплоемкостью при постоянном объеме Су [5-8]. Установим эту связь и покажем, что коэффициент пропорциональности между ними определяется усредненной должным образом средней длиной
свободного пробега фононов ((л^)) и средней скоростью фононов , а именно:
^ = |((Л 4 Ы с,.
(30)
Энергия фононов на единицу объема, очевидно, дается выражением:
00
Ерк = | (Гш)Бр11 {Тш) п0 {Тш) сЦПа), (31)
о
где ВрЬ{Ьа>) есть плотность фононных состояний. По определению,
д
Су = -£- = — f (tm)Dh(ha)n0(tm)d(tm) =
дТ oTJ
■\{tm)Dph{tia)
дп0{Ьсо) дТ
d(ilco) =
(32)
л2 k 2T 3
\Dph(ha)Wph(tm)d(tm),
где использовались уравнения (8) и (15).
Умножение и деление коэффициента теплопроводности (24) на теплоемкость (32) дает искомую пропорциональность:
1J M hJ
Mph(Tm)
А
Xph(ba)Wvh(tm)d(Tm)
ph
\Dph(Tm)Wph(Tm)d(Tm)
С, . (33)
Для получения окончательного результата в виде уравнения (30) и корректной интерпретации коэффициента пропорциональности между и нужно вернуться к связи между средней длиной свободного пробега и коэффициентом прохождения (19), имеющей одинаковый вид как для электронов, так и для фононов.
6. Средняя длина свободного пробега и коэффициент прохождения
Связь этих двух важнейших характеристик диффузионного транспорта проще всего проиллюстрировать на примере электронного транспорта по Ш проводнику в диффузионном режиме (рис. 7).
Рис. 7. К выводу связи между коэффициентом прохождения и средней длиной свободного пробега
на примере однородного 1D проводника
Рассмотрим однородный 1D проводник в диффузионном режиме. Левый контакт впрыскивает электроны с током I+ (x = 0). Доля электронов Т войдет в правый контакт с током I+ (x = L) = TI+ (x = 0). Оставшаяся доля вернется в левый контакт (рассеяние назад) с током I ~ (x = 0) = RI+ (x = 0). В отсутствии процессов рекомбинации T + R = 1. Предполагается также, что правый контакт идеальный, поглощает все входящие в него электроны. Результирующий ток, очевидно, будет I = (1- R)I+ (0) = TI+ (0). Так или иначе, в проводнике имеют место как прямые, так и обратные потоки электронов, и нам необходимо описать их пространственное распределение.
Определим обратное значение средней длины свободного пробега 1/X как вероятность (на единицу длины) обращения положительного потока электронов в отрицательный и наоборот. Именно исходя из этого определения, величину X и называют иногда средней длиной свободного обратного рассеяния (mean-free-pass for backscattering) или проще средней длиной свободного пробега. В диффузионном проводнике некоторая доля положительного потока в результате рассеяния назад обращается в отрицательный. В проводнике формируется отрицательный поток и его некоторая доля в результате рассеяния назад обращается, усиливая положительный поток. В результате градиент положительного тока складывается из двух величин:
0
0
dl+ (x)
I+ (x) +1 - (x)
(34)
^ Я Я
Пренебрегая процессами рекомбинации, результирующий ток
I = I+ (x) -1 - (x)
(35)
является постоянной величиной, так что градиент тока
dl+ (x)
dx
I
I
(36)
является константой. Другими словами, ток спадает линейно вдоль проводника:
I+ (x) = I+ (0) -- x .
I
(37)
Воспользуемся полученным уравнением для вычисления тока, входящего в правый контакт:
I+ (L) = I+ (0) - Ll = I+ (0) - [ I+ (L) -1 - (L)] L =
= I+ (0) -1+ (L)
(38)
где мы воспользовались уравнением баланса (35) и тем, что правый контакт идеальный (I - (Ь) = 0). Из последнего равенства находим
I+(L)=-
1
-I+ (0) = TI+ (0),
..........(39)
Я + Ь
если провести аналогичные рассуждения для впрыскивания электронов правым контактом, то получим аналогичное уравнение для электронов, входящих в левый контакт, а именно: I ~ (0) = Т' I ~ (Ь). Для однородного проводника
Т' = Т. Проводник под напряжением не является однородным, но нас интересует режим линейного отклика, так что вполне приемлемо положить Т'« Т. Окончательно, в предположении независимости друг от друга мод проводимости получаем искомое уравнение (19), связывающее коэффициент прохождения со средней длиной свободного пробега:
T (E) =-
1( E)
(40)
Я( Е) + Ь
Вывод уравнения (40) сделан в рамках простой модели, что никак не мешает успешному и широкому использованию его на практике. Важным моментом в проведенных рассуждениях является интерпретация 1 / Я как вероятности (на единицу длины) обращения потока частиц в обратном направлении в результате рассеяния. Именно поэтому саму длину Я, как уже упоминалось, часто называют средней длиной свободного рассеяния назад.
Установим связь между Я и временем рассеяния п. Пусть электрон совершает акт рассеяния в изотропном Ш проводнике. У него есть две возможности: рассеяться вперед и рассеяться назад. Только рассеяние назад существенно для
определения средней длины свободного пробега. Отсюда следует, что средняя длина рассеяния назад равна удвоенному значению средней длины рассеяния:
Я (Е) = 2Л(Е) = 2у(Е) т(Е) . (41)
Для проводника произвольной размерности средняя длина рассеяния назад дается выражением [14]
1(E) = 2
(К1):
(42)
где усреднение для 2D и 3D проводников ведется по углам. Для изотропных проводников
л
I2D (E) =- v(E)T(E) , 4
I3D (E) = 4 v(E)z(E) .
(43)
(44)
Для времени рассеяния п часто используется степенной закон
(45)
где значение показателя степени отражает тот или иной механизм рассеяния. Для акустических фононов в 3D проводниках с параболической дисперсией 5 = -1/2, для рассеяния на ионизированных примесях 5 = +3/2 [15].
Аналогичный степенной закон часто используется для средней длины рассеяния:
1(E) = 1
( E - Er
(46)
\ кТ
Для параболической зонной структуры у(Е) <х Е12, так что г = 5 +1/2. Для акустических фононов, соответственно, г = 0, а для рассеяния на ионизированных примесях г = 2 .
Возвращаясь к нашей первоначальной задаче -выводу уравнения (30) из (33), для 3D проводника, согласно (44), имеем:
4
Кь = з V (Й®) Трк (М = где, согласно (41),
^к (й®) тРк ) = АРк (й®) = так что окончательно
4
(47)
(48)
(49)
7. Плотность фононных состояний и фононные моды
Ранее было показано, что для 3D проводника число мод и плотность электронных состояний связаны соотношением (ф-ла (63) в [12]):
Ые1 (Е) = ЛМ»(Е) = Л-Ь+ (Е))Бз»(Е) . (50)
r
Мы хотим эту формулу переписать применительно к фононам. Сначала вспомним, что спиновое вырождение ^ = 2 входит в плотность состояний, так что в пересчете на один спин
(Е) = 2ДС (Е). Далее, для сферических зон в 3D
V (E))
Л, (E) 2
(51)
(52)
Сочетая (51) и (52) с (50) и переходя к фононам, получаем
h (V AhmV Mph(hm) = A-\-^^\2Dph(hm) =
.h
(53)
= A~vph(hm)Dph(hm)
Подставляя (49) и (53) в (33), получаем очередной промежуточный результат:
^\^ptl(ba,)vptl(ba,)Dptl(ba,)Wph(tia,)d(ha,)
о_
]Dph(hm)Wph(f,o,)d(ho,)
Cr
(54)
Умножая и деля (54) на
со
\уркфт)Оркфт)]¥ркфт)ё(Нт), (55)
о
окончательно получаем искомое уравнение (30) с коэффициентом пропорциональности между и С в виде произведения средней длины свободного пробега фононов
I фт)урк фт)ПрЬ (Нт) 1¥рЬ (Тгт) ¿(Па,) -
\урк (й а,)Пр, фа,) \¥рк (Па,) ¿(Па,)
о
и средней скорости фононов
ы
]vph (hm)Dph (hm) Wph (hm) d(hm)
со '
\Dph(hco)Wph(hco)d(hco)
(57)
усредненных должным образом.
Уравнение (30) часто используется для вычисления средней длины сводного пробега фононов из результатов измерений теплопроводности и теплоемкости. Для этого нужно знать среднюю скорость фононов, которая нередко предполагается равной продольной составляющей скорости звука. Приведенные выше выкладки точно определяют смысл усреднений при вычислении средней длины свободного пробега фононов и их средней скорости. Например, задавшись дисперсией фононов всегда можно вычислить их среднюю скорость по (57). Как правило, вычисленная таким
образом скорость сильно отличается от продольной составляющей скорости звука, что далее ведет к ошибочным результатам по средней длине свободного пробега фононов [16].
Рассмотрим плотность фононных состояний и число фононных мод с более общих позиций. Плотность состояний фононов определяется аналогично плотности состояний электронов, как число разрешенных фононных состояний на единицу энергии и единицу длины Ь для Ш проводников, площади поперечного сечения А для 2Б проводников и на единицу объема 0=ЬА для 3Б проводников. В случае электронов число разрешенных состояний со значением импульса меньшим, чем заданное значение р, для проводников размерности ё={ 1, 2, 3} [17]
N ( p )= 2-
L ■ W
4п L ■ A
(58)
^ И / р (И / р у 3 (И / р )3 ]
Аналогично, число разрешенных фононных состояний со значением волнового вектора к = р / й меньшим, чем заданное значение к,
NlD (к) =
N2D (к) =
N3D (к) =
2к 2л / L' лк2 4л2/ L2' 4лк3/3 8л3/ iL '
(59)
где единообразия ради далее удобно площадь 2D проводника обозначить как X2 = Ь ■ Ш, а объем для 3D проводника как Ь = 0 = Ь ■ А . Величину Ь условно назовем «объемом» проводника размерности ё.
Плотность фононных состояний, согласно определению выше, очевидно, такова:
Ad (к) = D2d (к) =
D3D(к) =
1 dNw (к) 1
L dк л
1 dN2D (к) к
L dк 2л
1 dN3D (к) к2
(60)
L ёк
2л2
Плотность фононных состояний часто выражают через групповую скорость фононов
vg (m) = dm / ёк,
(61)
а именно [18]:
1 Nd (к) 1 dN1 d (к) ёк D1D(m) = Т ~
1
D2D (Ю) =
D3D (m) =
L dm L ёк dm nvg (m) 1 dNw (к) 1 dNw (к) dh к(m) (62)
L dm L dk dm 2nvg (m) 1 dM3D(к) 1 Nd(к) к2(m)
L dm L dk dm 2л v (m)
0
Количество тепла, переносимого фононами, определяется выражением (5), в котором число состояний поляризации фононов учтено в числе фононных мод Мрк От). Перепишем это вьфажение через интеграл по частоте:
поперечном относительно направления движения фононов. Действительно, предполагая линейную зависимость а (к), что вполне обосновано для LA и TA акустических мод, имеем к = 2я/Я, так что уравнения (68) сводятся к
д=1-](Псо)Трк{со)Мрк{со){п1 -п2)с1ш . (63)
о
При переносе тепла по Ш проводнику Мрк (а)
есть просто единица, по 2D проводнику (а)
пропорционально ширине проводника, а по 3D проводнику - площади его поперечного сечения:
Мрк (а) = {Мш {а), WМгD (а), ЛМЗВ (а)}. (64) Для вычисления и интерпретации М2В (а) и Мзс (а) нам потребуются соответствующие выражения для потоков тепла
JQh = {JQh(1D), JQh(2D), JQh(3D)} = {Q, Q / w, Q / A},
(65)
связанных непосредственно с количеством переносимого тепла (33), а именно [18]:
jg{\D) = ]{hw)Tph{w)
о
J^(2D) = ](hw)Tph(w)
о
J^(3D) = ](hw)Tph(w)
vg (m)DiD (m)
(П - n2)dm,
1 2vg (m)
2 ж
1 vs (m)
2 2
D2d (m)
(n - n ) dm, (66)
D3D (m)
(n - n ) dm,
где плотности фононных состояний определяются уравнениями (62). Уравнивая потоки тепла (66) с количеством перенесенного тепла (63) согласно (65) и учитывая выражения для плотности фононных состояний (62), окончательно получаем:
Ml d (m) = ж [vg (m)] Dl d (m) = 1,
M2D (m) = ж
M3D (m) = ж
(m) ж
Vg (m)
D2 D (m) =
к (m)
(67)
D3 d (m) =
к 2(m) 4ж '
или иначе число фононных мод (64), определяющих количество перенесенного тепла (63)
1
Mph (m) =
M1D (m) = ж\уе (m)]
жv (m)
= 1,
WM2D (m) = Wж
AM3D (m) = A ж
2vg (m) ж
vg(m)'
к (m)
2^Vg (m) к 2(m)
2ж Vg (m)
(68)
что позволяет интерпретировать число мод (а)
как максимальное число де-бройлевских полуволн, укладывающихся в «сечении» проводника,
Mph (m) =
Mid (m) = 1, WM2D (m) = W / (X / 2), AM3B (m) = A /[(4/ ж)(Х/ 2)2
(69)
где множитель 4/ я появляется в связи с тем, что поперечное сечение 3D проводника предполагается круговым с площадью А, равной квадрату диаметра окружности, деленному на 4/я. Обратим также внимание на то, что выражения в квадратных скобках уравнений (68) есть ничто иное как усредненные по всем угловым переменным значения скорости фононов в направлении их движения (вдоль оси +х в нашем рассмотрении) [19].
8. Дебаевская модель
В случае электронов ширина зоны проводимости » кТ, так что заняты только состояния вблизи дна зоны проводимости и поэтому так неплохо работает модель эффективной массы с параболической дисперсией (рис. 1). Совсем другая ситуация в случае фононов: BW « кТ (рис. 4). Заняты состояния по всей зоне Бриллюэна. Акустические моды тем лучше аппроксимируются дебаевским линейным приближением (2), чем ближе волновые числа к центру зоны Бриллюэна (рис. 3 и 4).
Перепишем (2) в виде
hco = hvDq ,
(70)
где под дебаевской скоростью \в понимают среднюю скорость продольных и поперечных акустических мод. Тогда, по аналогии с электронами, плотность фононных состояний
2я (hvD)
(71)
где множитель 3 учитывает три поляризации акустических мод. Обратим внимание, что в учебниках плотность фононных состояний чаще
т—т -1 -3
выражается в Гц •м :
Dph(®)
Число мод на единицу площади поперечного сечения проводника вычисляем по (53):
Mph(ha>) =
3(hco)2 4 7i(hvD)2
(72)
Поскольку все состояния в зоне Бриллюэна при умеренных температурах стремятся быть заполненными, нужно проконтролировать, что действительно учтены все состояния. В кристалле с N / О атомами на единицу объема имеем 3Ы / О
состояний на единицу объема. Полное число состояний дается интегралом
J Dph фа) йфа) .
(73)
в котором верхний предел, определяемый дебаевской частотой, выбирается таким образом, чтобы получить правильное число состояний 3И / О, а именно:
hcoD = hvD
6 n2N Q
- T.
(74)
Дебаевская частота является предельной, выше которой фононные состояния уже не существуют или не учитываются. Это ограничение можно ввести и через дебаевский волновой вектор дв или через дебаевскую температуру Тв через (74). Для Т <к Тв заняты только состояния с (/—>(). для которых хорошо выполняется дебаевское приближение.
Теперь интегрирование в (24) можно вести только до дебаевской энергии:
тт2к2Тл? М Aha)
i — л— Kh WPh фю)йфа) (75)
3h
и брать число фононных мод по (72). Интеграл можно взять численно или аналитически с приемлемым выражением для средней длины свободного пробега фононов. Именно таким образом была впервые вычислена решеточная теплопроводность [20, 21]. Дальнейшее развитие теория и методы вычисления термоэлектрических транспортных коэффициентов получили в работах [14, 16, 22].
9. Рассеяние фононов
Фононы рассеиваются на незаряженных дефектах, примесных атомах, изотопах, на других фононах, на поверхностях и границах, на электронах. Фонон-фононное рассеивание обязано тому, что потенциальная энергия колеблющегося атома в кристалле не совсем гармоническая. Все более высокие порядки учитываются в виде потенциала рассеяния. Рассматривают две возможности фонон-фононного взаимодействия. В обычном процессе два фонона взаимодействуя порождают третий фонон с сохранением импульса и энергии:
(76)
+ Щ2 = Щ3
ы>х + ы>2 = ы>г
Поскольку суммарный импульс сохраняется, то такой процесс не может дать существенный вклад в перенос тепла.
В итк1арр/и-процессах рассеяния суммарный импульс не сохраняется, генерируемый фонон характеризуется импульсом за пределами первой зоны Бриллюэна. Это может быть как ангармоничное фонон-фононное, так и электрон-фононное взаимодействия. и-рассеяние является основным процессом в переносе тепла, особенно при высоких темпе-
ратурах. Ощутимый вклад дают также рассеяние на дефектахЮ и на границах^. Скорости всех трех процессов складываются, так что
1
1 1 1
- +-+ -
(77)
тифа) твфа) твфа)
или через соответствующие длины свободного пробега
1
1
1
1
, (78)
\jihco) Лвфа) 1вфа)
поскольку для 3D сферических дисперсий Я и т связаны уравнением (47).
Для всех механизмов рассеяния предложены выражения для скоростей [23]. Например, при рассеянии на дефектах
1
твфа)
■ ОС со
(79)
что напоминает рэлеевское рассеяние света на частицах пыли в атмосфере; высокочастотные фононы особенно чувствительны к рассеянию на малых дефектах. При рассеянии на границах и поверхностях
1
vph(hco)
(80)
тв фа) Ь
где Ь - наименьшая протяженность образца. Сложнее с ^/-рассеянием. Обычно используется выражение 1
тифа)
т3 2 -тп!ьт
■ ос 1 со е D
(81)
10. Температурная зависимость теплопроводности
Температурная зависимость решеточной теплопроводности к1 кремния представлена на рис. 8. Согласно уравнению (26), пропорциональна числу фононных мод (Мрк / А и среднему значению
длины свободного пробега фононов ^(Я)) .
Изменение теплопроводности с температурой качественно можно объяснить следующим образом. Из определения числа мод (25) можно показать, что при низких температурах
{Мрк) <* Т3; (Т ^ 0) (82)
быстрое увеличение числа заселенных мод с температурой вызывает рост теплопроводности. Важную роль играет рассеяние на границах. По мере роста температуры все больше продуцируется коротковолновых фононов (большие значения д). Возрастает роль рассеяния таких фононов на дефектах. По мере приближения температуры к дебаевской заселяются все фононные моды и дальнейшее увеличение температуры уже не увеличивает числа мод. Более того, высокие температуры порождают ^/-процессы рассеяния и теплопроводность начинает падать.
0
KL =
io4
ю1-
10° Ю1 Ю2 103
Т(К)
Рис. 8. Температурная зависимость решеточной теплопроводности кремния: экспериментальные точки по [24], расчетная кривая по [16]
11. Различие между решеточной теплопроводностью и электрической проводимостью
Уже подчеркивалось близкое подобие пары фононных транспортных уравнений (22) и (26) и пары электронных уравнений (28) и (29). Длины свободных пробегов электронов и фононов одного порядка. Почему же электрическая проводимость веществ меняется на много порядков, тогда как теплопроводность всего разнообразия веществ изменяется лишь на несколько порядков? Ответ кроется в различии между окнами проводимости для электронов (13) и фононов (15). Для тех и других рост температуры расширяет окна проводимости и тем самым увеличивает заселенность состояний. Для электронов, однако, все зависит от положения уровня Ферми, сдвигая который вверх или вниз можно менять проводимость на много порядков. Для фононов же ширина окна проводимости определяется только температурой.
Подчеркнем еще одно практически важное различие между описанием электронов и фононов, связанное с особенностями заселения состояний. В термоэлектрических устройствах Ер «Ес и окна проводимости электронов и фононов схожи. У электронов, однако, ширина зон очень большая и лишь небольшое число состояний вблизи дна зоны проводимости заселено. Хорошо работает параболическая дисперсия и легко получить аналитические решения. В сопоставимой ситуации у фононов ширина зон незначительная, при повышенных температурах практически все состояния в зоне Бриллюэна заселены и простые аналитические модели для вычисления теплопроводности не работают.
12. Квантование теплопроводности
По аналогии с квантованием электронной проводимости
О""" = 2- М (Е„ ) (83)
к
еще 30 лет назад Пендри [25] показал, что должны существовать квантовые ограничения на перенос
тепла. Действительно, в уравнении (17) при Т —>0 окно фононной проводимости Wh{tico) стремится к ¿-функции, так что в пределе
я2 k T
K =—^TPh(0)Mph(0) . (84)
Для массивных проводников А/^Дй®)—>0
при й® —> 0, однако, для наноразмерных проводников всегда есть конечное число фононных мод. В случае баллистического транспорта ^ = 1 и в этом случае можно ожидать что
я2 к T
Kl =—^MPh. (85)
Именно этот результат был подтвержден экспериментально на 4-хмодовом нанорезисторе при Т<0.8 К [13]; результаты измерений теплопроводности полностью согласуется с предсказаниями для одномерных баллистических резисторов [26-28].
Квант теплопроводности
g0 = — к T /3h (86)
есть максимальная энергия, переносимая в пересчете на одну фононную моду. Удивительно, что квант теплопроводности не зависит от статистики частиц, переносящих тепло, он универсален для фермионов, бозонов и любых аньенов/anyons [29-31].
13. Выводы
В настоящем обзоре подчеркивается, что концепции, используемые для описания электронного транспорта, с успехом могут быть перенесены на фононный транспорт. И в том и в другом случае подход Ландауэра, обобщенный впоследствии Даттой и Лундстромом, позволяет количественно описать транспортные процессы в проводниках любой размерности и при любых режимах транспорта - от диффузионного до баллистического. Решеточная теплопроводность и электронная проводимость описываются очень схожими по форме уравнениями. Есть, однако, два существенных отличия, связанных с физикой электронных и фононных процессов.
Первое отличие касается ширины зон. В случае электронов ширина зон намного превышает kT и при комнатной температуре заселены лишь уровни вблизи дна зоны проводимости. Ширина фононных дисперсий всего лишь порядка kT и при комнатной температуре все акустические моды в пределах всей зоны Бриллюэна практически заселены. В результате дебаевская модель для акустических мод далеко не столь успешно работает по сравнению с широкой применимостью параболической дисперсии для электронов.
Второе существенное различие связано с тем, что в случае электронов заселенность мод проводимости, контролируемая фермиевским окном проводимости, зависит не только от температуры, как это
имеет место в случае фононов, но и от положения уровня Ферми. Именно по этой причине электронная проводимость различных материалов может изменяться в очень широких пределах, тогда как теплопроводность изменяется в очень узких пределах всего на несколько порядков.
В основу настоящего обзора положены лекции Марка Лундстрома «Near-Equilibrium Transport: Fundamentals and Applications» [2], Суприе Датты «Fundamentals of Nanoelectronics, Part I: Basic Concepts» [1] и Тимоти Фишера «Thermal Energy at the Nanoscale» [18], прочитанных в 2011-2013 годах в рамках инициативы Purdue University/ nanoHUB-U [www.nanohub.org/u], а также наши статьи [32, 33].
14. Благодарности
Благодарю Н. Е. Кругляк за помощь в работе по подготовке рукописи к печати.
Литература
1. Supriyo Datta, Lessons from Nanoelectronics: A New Perspective on Transport [Electronic resource] / Hackensack, New Jersey: World Scientific Publishing Company, 2012. - Available at: www.nanohub.org/ courses/FoN1
2. Lundstrom, M. Near-Equilibrium Transport: Fundamentals and Applications [Electronic resource] / M. Lundstrom, C. JeongHackensack. - New Jersey: World Scientific Publishing Company. 2013. - Available at: www.nanohub.org/resources/11763
3. Kruglyak, Yu. Landauer-Datta-Lundstrom Generalized Transport Model for Nanoelectronics [Text] / Yu. Kruglyak // Journal of Nanoscience. - 2014. - Vol. 2015. -P. 15. doi: 10.1155/2014/725420
4. Kruglyak, Yu. A. A Generalized Landauer-Datta-Lundstrom Electron Transport Model [Text] / Yu. A. Kruglyak // Russian Journal of Physical Chemistry. - 2014. - Vol. 88, Issue 11. - P. 1826-1836. doi: 10.1134/s0036024414110119
5. Займан, Дж. Электроны и фононы. Теория явлений переноса в твердых телах [Текст] / Дж. Займан. -М: ИИЛ, 1962. - 488 с.
6. Займан, Дж. Принципы теории твердого тела [Текст] / Дж. Займан. - М: Высшая школа, 1974. - 472 с.
7. Киттель, Ч. Введение в физику твердого тела [Текст] / Ч. Киттель. - М: Наука: 1978. - 791 с.
8. Ашкрофт, Н. Физика твердого тела, тома 1 и 2 [Текст] / Н. Ашкрофт, Н. Мермин. - М: Мир: 1979. - 486 с.
9. Mohr, M. Phonon dispersion of graphite by inelastic x-ray scattering [Text] / M. Mohr, J. Maultzsch, E. Dobardzic, S. Reich, I. Milosevic, M. Damnjanovic, A. Bosak, M. Krisch, C. Thomsen // Physical Review B. - 2007. - Vol. 76, Issue 3. doi: 10.1103/physrevb.76.035439
10. Елецкий, А. В. Графен: методы получения и теплофизические свойства [Текст] / А. В. Елецкий, И. М. Искандарова, А. А. Книжник, Д. Н. Красиков // Успехи физических наук. - 2011. - Т. 181, № 3. -С. 227-258. doi: 10.3367/UFNr.0181.201103a.0233
11. Katsnelson, M. I. Graphene: Carbon in Two Dimensions [Text] / M. I. Katsnelson. - New York: Cambridge University Press, 2012. doi: 10.1017/cbo9781139031080
12. Кругляк, Ю. А. Обобщенная модель электронного транспорта Ландауэра-Датты-Лундстрома [Текст] / Ю. А. Кругляк // Nanosystems, Nanomaterials, Nanotech-nologies, 11, № 3: 519 - 549 (2013); Erratum: ibid, 12, № 2: 415 (2014).
13. Schwab, K. Measurement of the quantum of thermal conductance [Text] / K. Schwab, E. A. Henriksen, J. M. Worlock, M. L. Roukes // Nature. - 2000. - Vol. 404. - P. 974.
14. Jeong, C. On Landauer vs. Boltzmann and Full Band vs. Effective Mass Evaluation of Thermoelectric Transport Coefficients [Text] / C. Jeong, R. Kim, M. Luisier, S. Datta, M. Lundstrom // Journal of Applied Physics. - 2010. -Vol. 107, Issue 2. - P. 023707. doi: 10.1063/1.3291120
15. Lundstrom, M. Fundamentals of Carrier Transport [Text] / M. Lundstrom. - Cambridge UK: Cambridge University Press, 2012.
16. Jeong, C. Full Dispersion vs. Debye Model Evaluation of Lattice Thermal Conductivity with a Landauer Approach [Text] / C. Jeong, S. Datta, M. Lundstrom // Journal of Applied Physics. - 2011. - Vol. 109. - P. 073718/8. doi: 10.1063/1.3567111
17. Кругляк, Ю. А. Уроки наноэлектроники. 3. Электронная проводимость и моды проводимости в концепции «снизу - вверх» [Текст] / Ю. А. Кругляк, Н. Е. Кругляк // Физическое образование в вузах. - 2013. - Т. 19, № 3. -С. 99-110.
18. Timothy S. Fisher, Thermal Energy at the Nanoscale [Electronic resource] / Hackensack, New Jersey: World Scientific Publishing Company, 2013. - Available at: www.nanohub.org/courses/2
19. Кругляк, Ю. А. Уроки наноэлектроники. 2. Модель упругого резистора и новая формулировка закона Ома в концепции «снизу - вверх» [Текст] / Ю. А. Кругляк, Н. Е. Кругляк // Физическое образование в вузах. - 2013. - Т. 19, № 2. - С. 161-173.
20. Callaway, J. Model for lattice thermal conductivity at low temperatures [Text] / J. Callaway // Physical Review. -1959. - Vol. 113, Issue 4. - P. 1046-1051. doi: 10.1103/ physrev.113.1046
21. Holland, M. G. Analysis of lattice thermal conductivity [Text] / M. G. Holland // Physical Review. -1963. - Vol. 132, Issue 6. - P. 2461-2471. doi: 10.1103/ physrev.132.2461
22. Jeong, C. Thermal conductivity of bulk and thin-film silicon: a Landauer approach [Text] / C. Jeong, S. Datta, M. Lundstrom // Journal of Applied Physics. -2012. -Vol. 111. - P. 093708. doi: 10.1063/1.4710993
23. Gang, C. Nanoscale Energy Transport and Conversion: A Parallel Treatment of Electrons, Molecules, Phonons, and Photons [Text] / C. Gang. - New York: Oxford University Press, 2005. - 560 p.
24. Glassbrenner, C. J. Thermal Conductivity of Silicon and Germanium from 3°K to the Melting Point [Text] / C. J. Glassbrenner, G. A. Slack // Physical Review. - 1964. -Vol. 134, Issue 4A. - P. A1058-A1069. doi: 10.1103/ physrev.134.a1058
25. Pendry, J. B. Quantum limits to the flow of information and entropy [Text] / J. B. Pendry // Journal of Physics A: Mathematical and General. - 1983. - Vol. 16, Issue 10. - P. 2161-2171. doi: 10.1088/0305-4470/16/10/012
26. Angelescu, D. E. Heat transport in mesoscopic systems [Text] / D. E. Angelescu, M. C. Cross, M. L. Rou-kes // Superlattices and Microstructures. - 1998. -Vol. 23, Issue 3-4. - P. 673-689. doi: 10.1006/spmi.1997.0561
27. Rego, L. G. C. Quantized Thermal Conductance of Dielectric Quantum Wires [Text] / L. G. C. Rego, G. Kirczenow // Physical Review Letters. - 1998. - Vol. 81, Issue 1. - P. 232-235. doi: 10.1103/physrevlett.81.232
28. Blencowe, M. P. Quantum energy flow in mesoscopic dielectric structures [Text] / M. P. Blencowe // Physical Review B. - 1999. - Vol. 59, Issue 7. - P. 4992-4998. doi: 10.1103/physrevb.59.4992
29. Rego, L. G. C. Fractional exclusion statistics and the universal quantum of thermal conductance: A unifying approach [Text] / L. G. C. Rego, G. Kirczenow // Physical Review B. - 1999. - Vol. 59, Issue 20. -P. 13080-13086. doi: 10.1103/physrevb.59.13080
30. Krive, I. V. Transport properties of quasiparticles with fractional exclusion statistics [Text] / I. V. Krive, E. R. Mucciolo // Physical Review B. - 1999. - Vol. 60, Issue 3. - P. 1429-1432. doi: 10.1103/physrevb.60.1429
31. Caves, C. M. Quantum limits on bosonic communication rates [Text] / C. M. Caves, P .D. Drummond // Reviews of Modern Physics. - 1994. - Vol. 66, Issue 2. -P. 481-537. doi: 10.1103/revmodphys.66.481
32. Кругляк, Ю. O. Уроки наноелектрошки. Термоелектричш явища в концепцй «знизу - вгору» [Text] / Ю. O. Кругляк, Н. Ю. Кругляк, М. В. Стр1ха // Sensor Electronics Microsys. Tech. - 2013. - Т. 13, № 1. - С. 6-21.
33. Кругляк, Ю. А. Уроки наноэлектроники. 5. Перенос тепла фононами в концепции «снизу - вверх» [Текст] / Ю. А. Кругляк // Физическое образование в вузах. - 2014. - Т. 20, № 1. - С. 39-43.
References
1. Supriyo Datta, Lessons from Nanoelectronics: A New Perspective on Transport (2012). Hackensack, New Jersey: World Scientific Publishing Company. Available at: www.nanohub.org/courses/FoN1
2. Lundstrom, M., Jeong, C. (2013). Near-Equilibrium Transport: Fundamentals and Applications. Hackensack, New Jersey: World Scientific Publishing Company. Available at: www.nanohub.org/resources/11763
3. Kruglyak, Yu. (2014). Landauer-Datta-Lundstrom Generalized Transport Model for Nanoelectronics. Journal of Nanoscience, 2014, 15. doi: 10.1155/2014/725420
4. Kruglyak, Yu. A. (2014). A Generalized Landauer-Datta-Lundstrom Electron Transport Model. Russian Journal of Physical Chemistry, 88 (11), 1826-1836. doi: 10.1134/ s0036024414110119
5. Ziman, J. M. (1960). Electrons and phonons. The theory of transport phenomena in solids. Oxford at the Clarendon Press, Oxford, 488.
6. Ziman, J. M. (1964). Principles of the theory of solids. Cambridge University Press, Cambridge, 472.
7. Kittel, C. (1971). Introduction to solid state physics, 4th ed. John Wiley and Sons, New York, 791.
8. Ashcroft, N. W., Mermin, N. D. (1979). Solid State Physics (Philadelphia: Suanders College, 486.
9. Mohr, M., Maultzsch, J., Dobardzic, E., Reich, S., Milosevic, I., Damnjanovic, M., Bosak, A., Krisch, M., Thomsen, C. (2007). Phonon dispersion of graphite by inelastic x-ray scattering. Physical Review B, 76 (3), 035439/7. doi: 10.1103/physrevb.76.035439
10. Eletskii, A. V., Iskandarova, I. M., Knizhnik, A. A., Krasikov, D. N. (2011). Graphene: fabrication methods and thermophysical properties. Physics Uspekhi, 54, 227-258. doi: 10.3367/UFNe.0181.201103a.0233
11. Katsnelson, M. I. (2012). Graphene: Carbon in Two Dimensions. New York: Cambridge University Press. doi: 10.1017/cbo9781139031080
12. Kruglyak, Yu. A. (2013). The Generalized Landauer-Datta-Lunstrom Electron Transport Model. Nanosystems, Nanomaterials, Nanotechnologies, 11 (3), 519-549. Erratum: ibid, 12 (2), 415.
13. Schwab, K., Henriksen, E. A., Worlock, J. M., Roukes, M. L. (2000). Measurement of the quantum of thermal conductance. Nature, 404, 974.
14. Jeong, C., Kim, R., Luisier, M., Datta, S., Lundstrom, M. (2010). On Landauer vs. Boltzmann and Full
Band vs. Effective Mass Evaluation of Thermoelectric Transport Coefficients. Journal of Applied Physics, 107 (2), 023707. doi: 10.1063/1.3291120
15. Lundstrom, M. (2012). Fundamentals of Carrier Transport. Cambridge UK: Cambridge University Press.
16. Jeong, C., Datta, S., Lundstrom, M. (2011). Full Dispersion vs. Debye Model Evaluation of Lattice Thermal Conductivity with a Landauer Approach. Journal of Applied Physics, 109, 073718/8. doi: 10.1063/1.3567111
17. Kruglyak, Yu. A., Kruglyak, N. E. (2013). Lessons of nanoelectronics. 3. Electronic conductivity and conductivity modes by «bottom - up» approach. Physics in Higher Education, 19 (3), 99-110.
18. Fisher, T. S. (2013). Thermal Energy at the Nanoscale (Hackensack, New Jersey: World Scientific Publishing Company. Available at: www.nanohub.org/ courses/2
19. Kruglyak, Yu. A., Kruglyak, N. E. (2013). Lessons of nanoelectronics. 2. Elastic resistor model and new Ohm's law by «bottom - up» approach. Physics in Higher Education, 19 (2), 161-173.
20. Callaway, J. (1959). Model for lattice thermal conductivity at low temperatures. Physical Review, 113 (4), 1046-1015 doi: 10.1103/physrev.113.1046
21. Holland, M. G. (1963). Analysis of lattice thermal conductivity. Physical Review, 132 (6), 2461-2471. doi: 10.1103/physrev. 132.2461
22. Jeong, C., Datta, S., Lundstrom, M. (2012). Thermal conductivity of bulk and thin-film silicon: a Landauer approach. Journal of Applied Physics, 111, 093708. doi: 10.1063/1.4710993
23. Gang, C. (2005). Nanoscale Energy Transport and Conversion: A Parallel Treatment of Electrons, Molecules, Phonons, and Photons. New York: Oxford University Press, 560.
24. Glassbrenner, C. J., Slack, G. A. (1964). Thermal Conductivity of Silicon and Germanium from 3°K to the Melting Point. Physical Review, 134 (4A), A1058-A1069. doi: 10.1103/physrev.134.a1058
25. Pendry, J. B. (1983). Quantum limits to the flow of information and entropy. Journal of Physics A: Mathematical and General, 16 (10), 2161. doi: 10.1088/0305-4470/16/10/012
26. Angelescu, D. E., Cross, M. C., Roukes, M. L. (1998). Heat transport in mesoscopic systems, Superlattices and Microstructures, 23 (3-4), 673-689.
doi: 10.1006/spmi.1997.0561
27. Rego, L. G. C., Kirczenow, G. (1998). Quantized Thermal Conductance of Dielectric Quantum Wires. Physical Review Letters, 81 (1), 232-235. doi: 10.1103/ physrevlett.81.232
28. Blencowe, M. P. (1999). Quantum energy flow in mesoscopic dielectric structures. Physical Review B, 59 (7), 4992-4998. doi: 10.1103/physrevb.59.4992
29. Rego, L. G. C., Kirczenow, G. (1999). Fractional exclusion statistics and the universal quantum of thermal conductance: A unifying approach. Physical Review B, 59 (20), 13080-13086. doi: 10.1103/physrevb.59.13080
30. Krive, I. V., Mucciolo, E. R. (1999). Transport properties of quasiparticles with fractional exclusion statistics. Physical Review B, 60 (3), 1429-1432.
doi: 10.1103/physrevb.60.1429
31. Caves, C. M., Drummond, P. D. (1994). Quantum limits on bosonic communication rates, Reviews of Modern Physics, 66 (2), 481-537. doi: 10.1103/revmodphys.66.481
32. Kruglyak, Yu. A., Kruglyak, N. Yu., Strikha, M. V. (2013). Lessons of nanoelectronics. Thermoelectric phenomena in «bottom - up» approach, Sensor Electronics Microsys. Tech., 13 (1), 6-21.
33. Kruglyak, Yu. A. (2014). Lessons of ach. Physics in Higher Education, 20 (1), 39-43. nanoelectronics. 5. Phonon transport in «bottom - up» appro-
Рекомендовано до публтаци д-р ф1з.-мат. наук Глушков О.В.
Дата надходження рукопису 14.01.2015
Кругляк Юрий Алексеевич, доктор химических наук, профессор, кафедра информационных технологий, Одесский государственный экологический университет, ул. Львовская, 15, г. Одесса, Украина, 65016 E-mail: [email protected]
УДК 537.32
DOI: 10.15587/2313-8416.2015.36443
ГРАФЕН В ТРАНСПОРТНОЙ МОДЕЛИ ЛАНДАУЭРА-ДАТТЫ-ЛУНДСТРОМА © Ю. А. Кругляк
Обсуждаются такие свойства графена как плотность электронных состояний и носителей тока, число мод и максимальная проводимость, рассеяние и подвижность в графене, циклотронная частота и эффективная масса, плотность фононных состояний, сравнительный вклад электронов и фононов в теплопроводность графена. В справочных целях дается сводка термоэлектрических коэффициентов для графена в баллистическом и диффузионном режимах проводимости со степенным законом рассеяния
Ключевые слова: нанофизика, наноэлектроника, графен, число мод, максимальная проводимость, эффективная масса, фононные состояния, теплопроводность, термоэлектрические коэффициенты
There are discussed the following properties of grapheme such as the density of electronic states and current carriers, the number of modes and maximum conductivity, scattering and mobility in graphene, the cyclotron frequency and the effective mass, phonon density of states, the relative contribution of electrons and phonons in the thermal conductivity of graphene. For reference purposes a summary of thermoelectric coefficients for graphene in ballistic and diffusive conduction regimes with the power law of scattering is given Keywords: nanophysics, nanoelectronics, graphene, mode numbers, maximum conductivity, effective mass, phonon states, thermal conductivity, thermoelectric coefficients
1. Введение
Модель электронного транспорта Ландауэра-Датты-Лундстрома (ЛДЛ) в режиме линейного отклика сводится к паре транспортных уравнений [1]
d ( Ер / q ) dT 2
Jx = а—--- - S а- [А/м ]
dx dx
diEp/ф + sdT,
dx dx
dT
JQx =KJ, - к-
dx
■ = 1/ p = fa'(E)dE,
а'(Е ) =
h W
5E
(1)
d (E„ / q) dT r„ ,2 =F - K0 -T [Вт/м2 (2)
dx dx
или в обращенной форме
(3)
(4)
с транспортными коэффициентами для 2D проводников
(5)
2q М(Е) Я(Е) [-50,, (6)
S = -k 1 ff E-EPL I а'(E)dE = ^^ = £, (7) q аf ^ kT ) а T
к°=Т (к ] ]'(Е)йЕ=к+л8а • (8) Уравнения (1)-(8) справедливы для диффузионных резисторов любой размерности и любого масштаба, любого типа проводимости, включая биполярную проводимость. Уравнение (6) выписано для 2D проводника; это видно по множителю М(Е) / IV, где IV - ширина 2D проводника. При рассмотрении 3D проводников этот множитель нужно заменить на М(Е) / А, где А - площадь поперечного сечения 3D проводника, а для Ш проводника заменить просто на число мод проводимости М(Е). При рассмотрении транспорта во всех режимах от диффузионного до баллистического среднюю длину свободного пробега Л(Е) нужно заменить на кажущуюся длину по
уравнению 1 / Ларр = 1/ Л +1/ Ь, где Ь - длина проводника. Фактически в этих уравнениях нужно задаться только двумя параметрами, М(Е) и Л(Е). Число мод для проводников любой размерности с параболической дисперсией определяется форму-
