Потік тепла, отже, не змінюється при N-процесі.
Такий простий доказ справедливий у відсутність дисперсії, але насправді результат виявляється абсолютно загальним. Один із способів показати це – досліджувати вплив N-процесів на «зміщений» фононний розподіл, що задається формулою(3.1)
У формулі (3.1) вектор u направлений так само, як і потік тепла, і представляє собою швидкість дрейфу. Розподіл відповідає рівноважному, але в системі координат, яка рухається з швидкістю u. Такий розподіл для одновимірного випадку показаний на рис. 3.1, так що з його асиметричного вигляду ясно, що існує потік тепла в праву сторону.
Вираз для швидкості зміни N(q) внаслідок N-процесів, отримане в теорії збурень, містить вірогідність процесу, при якому фонон цієї моди з'являється як кінцевий продукт взаємодії, а також вірогідність процесу, при якому відбувається зникнення фонона цієї моди і в результаті з'являються два інших фонона. Якщо розглянути процес, при якому q1 + q2 → q3, і відповідний зворотний процес, то швидкість зміни N(q1) визначається різницею між виразом від прямих процесів
і виразом від зворотних процесів
(такі вагові множники з'являються завжди в системах, що підкоряються статистиці Бозе – Ейнштейна).
Рис. 3.1. Одновимірний асиметричний розподіл фононів,
на який не впливають N-процеси.
Розподіл відповідає потоку тепла вправо. Воно описується формулою (2.1) при u = 0,1 ω/q, Т = 100К. По ординаті відкладено число фононів на одиничний інтервал q.
Якщо підставити формулу (3.1) для N(q) в ці два вирази і взяти їх різницю, то просте обчислення приводить до результату:
він визначає швидкість зміни N(q1). Якщо умови (2.1) справедливі при g = 0, тобто мають місце тільки N-процеси, то видно, що швидкість зміни N(q1) перетворюється в нуль. Для всіх пар значень q2 і q3, які можуть брати участь в N-процесах, виходить той же результат, як при процесах q1 → q2 + q3, що впливають на число фононів q1. Зміщений фононний розподіл, що визначається формулою (3.1) і відповідає ненульовому потоку тепла, таким чином, не змінюється внаслідок N-процесів.
Той же результат можна отримати і за допомогою варіаційного принципу. Для трифононних N-процесів чисельник перетвориться на вигляд:
(3.2)
Якщо записати
(3.3)
де u' – довільний вектор, то чисельник міститиме вираз (q1 + q2 – q3) • u', яке обертається на нуль, якщо q1 + q2 = q3. При такій пробній функції тепловий опір обертається на нуль, що, звичайно, відповідає мінімуму. Таким чином, за допомогою варіаційного принципу знову знаходимо, що у випадку тільки N-процесів теплопровідність нескінченна.
Легко показати, що для малих відхилень від рівноваги вираз для N(q), отримуване за допомогою формули (3.3), співпадає з формулою (3.1). З формули (2.1) маємо
для випадку, коли мало, можна написатиПідставляючи , знаходимо
Обидва ці вирази співпадають, якщо ħu = u'. Значення u і u' довільні й визначаються величиною потоку тепла; для заданого потоку тепла відхилення від рівноважного розподілу співпадають в обох випадках.
Не слід думати, проте, що оскільки N-процеси самі по собі не приводять до появи теплового опору, то ними можна зовсім нехтувати. Вони можуть істотно впливати на теплопровідність, якщо інтенсивності інших процесів розсіяння залежать від частоти; у такій ситуації N-процеси заважають модам, які розсіюються внаслідок цих процесів, «зноситися» потоком тепла. Багато зусиль було витрачено для того, щоб пояснити, як N-процеси спільно з процесами, які приводять до опору, визначають теплопровідність.
4. ОБЛІК НОРМАЛЬНИХ ПРОЦЕСІВ
4.1 Релаксаційний метод
У простій формі релаксаційного методу передбачається, що кожен механізм розсіяння характеризується часом релаксації, який для даної моди не залежить від населеності фононів у всіх інших модах. Якщо є декілька механізмів розсіяння, то швидкості
, , ... складаються і сумарний час релаксації τ(х), визначається виразомЦей вираз стає, звичайно, складнішим, якщо не користуватися спрощуючими припущеннями теорії Дебая про відсутність дисперсії і квадратичний закон для щільності мод: f(ω) ~ ω2.
Оскільки N-процеси самі по собі не приводять до встановлення рівноважного розподілу фононів, то вони не можуть входити в суму для τ(х) на тих же підставах, що і процеси, що ведуть до встановлення рівноваги (резистивне розсіяння). Проте ними не можна нехтувати, оскільки, перерозподіляючи енергію між модами, вони роблять відчутним для всіх мод наявність резистивних процесів розсіяння, залежних від частоти.
Для аналізу експериментальних даних по теплопровідності широко використовується розгляд Каллуея. Він припустив, що N-процеси переводять будь-який розподіл фононів, що відповідає деякому потоку тепла, в розподіл, визначуваний формулою (3.1), відповідний тому ж потоку тепла і далі вже не змінний внаслідок N-процесів. Час релаксації для таких процесів є τN (для простоти залежність часу релаксації від q, поляризації і температури не указується). Повна швидкість зміни N(q) дається тоді виразом
де у величину τR вносять внесок тільки процеси, що приводять до встановлення рівноважного розподілу N0(q). У моделі Дебая Каллуей, ввівши комбінований час релаксації
, отримав наступний вираз для теплопровідності:(4.1.1)
де
(4.1.1.а)
і(4.1.1.б)
Цей результат, як видно, знаходиться відповідно до твердження про те, що нормальні процеси впливають на теплопровідність, але трохи інакше, ніж чисто резистивні процеси. У формулу для ϰ1 нормальні процеси входять на тих же підставах, як і інші процеси, оскільки між ними не робиться ніякої відмінності у виразі для τС. Тому зазвичай вважається, що ϰ1 дає занижену оцінку теплопровідності, проте є другий член ϰ2, який декілька заповнює її «втрату».
Для пояснення експериментальних результатів часто необхідно користуватися повною формулою (4.1.1); обчислення, проте, дуже громіздкі, і корисно розглянути загальні результати, які виходять в трьох граничних випадках.
1) Випадок переважання резистивного розсіяння
Для кристала з великою кількістю дефектів всі моди сильно розсіваються унаслідок резистивних процесів; тоді для всіх мод τN >> τR , отже, τC ≈ τR. У такому разі ϰ2 << ϰ1 (якісно це можна зрозуміти, припустивши, що всі часи релаксації не залежать від частоти, тому при порівнянні ϰ1 та ϰ2 інтеграли скорочуються і ми маємо ϰ2/ ϰ1 = τR/τN << 1). Пізніше буде видно, що це порівняльно простий вираз придатний для аналізу експериментальних даних по теплопровідності не дуже ідеальних кристалів.
2) Випадок переважання N-процесів за наявності резистивного розсіяння
В цьому випадку час релаксації τC головним чином визначається N-процесами; тоді τR >> τN і τC ≈ τN. Звідси легко побачити, що ϰ2 >> ϰ1 (якісно це можна зрозуміти, припустивши незалежність часів релаксації від частоти, і отримати ϰ/ϰ1 = τR/τN >> 1). Для коефіцієнта теплопровідності тоді маємо
Перш за все дивно, що формула (4.1.2), яка визначає теплопровідність у разі переважання N-процесів, не містить τN. Проте N-процеси впливають на розподіл фононів і приводять його до форми (3.1). Коли N-процеси грають домінуючу роль, розподіл фононів стає «зміщеним» і не залежить від інтенсивності N-процесів. Тепловий опір виникає внаслідок резистивних процесів, що діють на цей розподіл.
Інший цікавий аспект формули (4.1.2) видно, якщо з її допомогою записати тепловий опір: