- •I Статистическая физика 8
- •8 Неидеальные газы 91
- •Часть I
- •1.1 Статистический вес состояния
- •1.2 Микроканоническое распределение
- •1.3 Энтропия и температура
- •1.4 Работа при постоянной энтропии
- •1.5 Каноническое распределение
- •1.6 Усреднение по времени и по ансамблю
- •VimTfl2 лт
- •2.1 Распределение Максвелла
- •2.2 Термодинамические функции идеального газа
- •2.3 Теплоемкость двухатомного газа
- •2.4 Отсутствие магнетизма классического газа. Диамагнетизм Ландау
- •3.1 Некоторые условия равновесия
- •3.2 Большое каноническое распределение
- •3.3 Энтропия идеального газа в неравновесном состоянии
- •4.1 Степень ионизации газа
- •4.2 Закон действующих масс
- •4.3 Тепловой эффект реакции
- •5.1 Задача о компьютерном газе биллиардных шаров 5.1.1 Распределение шаров по скоростям
- •5.2.2 Шары в пенале
- •5.2.3 Циклические граничные условия
- •5.3 Статистическое моделирование
- •7.1 Идеальный ферми-газ
- •7.4 Теплоемкость вырожденного электронного газа
- •7.5 Полупроводники
- •7.6 Идеальный бозе-газ
- •7.7 Конденсация Бозе-Эйнштейна
- •7.8 Газ фотонов
- •7.9 Теплоемкость твердого тела
- •7.10 Тепловое расширение твердых тел
- •7.11 Силы притяжения между нейтральными проводниками (силы Казимира)
- •8.1 Вириальный коэффициент
- •8.3 Учет взаимодействия в плазме
- •8.4 Пространственная дисперсия диэлектрической проницаемости плазмы (в электростатических задачах)
- •9.1 Фазовые переходы первого рода
- •9.2 Фазовые переходы второго рода
- •9.3 Теория Вейсса
- •9.4 Модель Изинга
- •9.5 Теория Ландау фазовых переходов второго рода
- •9.6 Критические индексы
- •10.1 Квазистационарные флуктуации
- •10.2 Флуктуации числа частиц
- •10.3 Рассеяние света
- •10.4 Дублет Мандельштама — Бриллюена
- •10.5 Флуктуации параметра порядка
- •11.1 Корреляционная функция скоростей и диффузия
- •11.2 Корреляционная функция случайных сил
- •12.1 Флуктуации в электрических цепях
- •12.2 Спектральное разложение флуктуации
- •12.3 Пример применения формулы Найквиста
- •Часть II
- •14.1 Кинетическое уравнение
- •14.2 Ленгмюровские колебания плазмы
- •14.3 Затухание Ландау
- •14.4 Интеграл столкновений
- •14.5 Уравнения гидродинамики
- •14.7 Кинетическое уравнение для электронов в металле
- •14.8 Теплопроводность электронного газа
- •14.9 Термоэлектрические эффекты
- •14.10 Квантовое кинетическое уравнение
7.9 Теплоемкость твердого тела
Имеется в виду кристаллическая решетка, о теплоемкости электронного газа речь уже шла.
Движение атомов друг относительно друга можно рассматривать как малые колебания. Вспомним, что для цепочки атомов, была зависимость частот от номера нормального колебания. В кристалле ситуация аналогична.
Атомы в кристалле можно представлять набором осцилляторов с разными частотами. При температурах, высоких по сравнению с наибольшей из этих частот, на каждый осциллятор приходится теплоемкость с = 1, на N атомов — С = SN — это закон Дюлонга и Пти.
При низких температурах теплоемкость падает. Эйнштейн, впервые применивший квантовый подход к вычислению теплоёмкости кристаллической решётки, принял, что все частоты этих осцилляторов одинаковы. Тогда при малых Т теплоемкость ос е-^о/г Эксперимент дает ос Т3.
Дебай предположил, что спектр примерно такой же, как для колебаний упругого тела, ш = ик, где к — волновой вектор, и — скорость звука, некоторым образом усредненная (см. ниже). Полное же число нормальных колебаний было выбрано равным числу степеней свободы, что обеспечивает выполнение в модели закона Дюлонга и Пти.
При низких температурах осцилляторы с высокими частотами не возбуждаются, поэтому ограничение спектра не сказывается и получается такая же зависимость теплоемкости от температуры, как для фотонного газа С ос Т3.
В промежуточной области температур модель Дебая можно считать интерполяционной.
Удобно ввести циклические граничные условия и рассматривать в качестве нормальных колебаний бегущие волны. Нужно учитывать, что
81
возможны разные поляризации звуковых волн, в изотропном теле — продольные (I) и поперечные (t). Число колебаний
Vd?k
3Vu2du
~
^ (2тг)3«?
i=l,h,h V ' 1
где и — средняя скорость звука в решетке, которая определяется так:
3 _ 1 2
W"1 Uf Щ
Максимальная частота в дебаевской модели принимается одной и той же для продольных и поперечных волн и определяется из условия:
SV
= 3N
Обычно вводят дебаевскую температуру
TD = hwm = ( бтг — J Теперь легко записать энергию
Ни.
Е =
л Td/T ,
3VT4 Г /1
= 3N
'т
где
е* -
при при
— функция Дебая (ее значения имеются в таблицах, что имело смысл до распространения компьютеров).
При Г < TD, E = (9/8)iVTD+(37r4/5)iV(r4/^D)5 c = (12тг4/5)ЛГ(Т3/Т|
При Г > Гд, Е = 3NT, С = 3N.
Теория Дебая оказалась очень удачной.
82
Примеры значений дебаевской температуры: очень высокая — для алмаза (Td=1860K) и низкая — для свинца (Тд=88К).
Тепловое движение можно считать газом фононов, каждое нормальное колебание — сорт фононов, N^ — число фононов с волновым вектором к. Отвлекаясь от возможных явлений, связанных поверхностью, удобно принять циклические граничные условия для колебаний решетки. Тогда каждый фонон можно представлять себе как бегущую плоскую волну с определенными поляризацией и волновым вектором.
Истинный спектр колебаний известен из опытов по рассеянию нейтронов. При неупругом рассеянии проявляются отдельные фононы, и законы сохранения позволяют по потере импульса и энергии нейтроном узнать зависимость энергии фонона от его квазиимпульса. Спектр оказался не очень похож на дебаевский.
В выражении энергии есть слагаемое 9NTd/8 — вклад нулевых колебаний. Оно не проявляется в теплоемкости. Однако в других явлениях (рассеяние рентгеновских лучей, эффект Моссбауера и др.) эти колебания полностью проявляются.
ЗАДАЧИ
1. В модели Дебая найти средний квадрат амплитуды колебаний атома.
2.