Для объяснения первого из перечисленных фактов рассмотрим основы квантовой теории теплоемкости, в которой колебания атомов в кристалле описываются совокупностью квазичастиц – фононов, энергия которых ei подчиняется законам квантовой статистики.
Каждому кванту энергии упругой волны удобно сопоставить фонон с энергией
. Тепловая энергия
тела будет суммой энергий фононов
|
(6.4) |
Суммирование
в формуле (6.4) проводится по всем разрешенным значениям частот, заключенных в
первой зоне Бриллюэна. С учетом функции распределения фононов по частотам (частотный спектр фононов), тепловая энергия тела в
интегральной форме будет иметь вид
|
(6.5) |
где − среднее
значение энергии фонона.
Следовательно, для определения тепловой энергии кристалла, а затем и его теплоемкости необходимо знать функцию распределения D(w) и среднюю энергию тепловых колебаний атомов.
В гл 5 было показано (формула 5.61), что энергия гармонического
квантового осциллятора может быть представлена в виде двух слагаемых: энергии
нулевых колебаний атомов и энергии
. Нулевые колебания не несут тепловой энергии и без учета
теплового расширения (в гармоническом приближении) от температуры не зависят, а
второе слагаемое характеризует возбужденное состояние системы. квант энергии
возбужденного состояния называют фононом.
Число возбужденных фононов n, имеющих энергию , зависит от величины возбуждающей энергии. Если это
тепловая энергия, то каждый раз, когда температура возрастает на
, амплитуда колебаний
возрастает на
величину, определяемую из условия
и равную
|
(6.6) |
Энергия колебаний с частотой возрастает при
этом на величину энергии фонона.
Очевидно, что
при данной температуре число возбужденных фононов будет максимально у наиболее
низкочастотных колебаний. Это число убывает при увеличении частоты .
Заметим, что величина одного кванта энергии , т. е. минимальная величина энергии, необходимая для
возбуждения наиболее высокочастотных колебаний в кристаллической решетке, не
является малой величиной. Сравнивая
с энергией
классического осциллятора при некоторой температуре
, т. е. приравнивая
и подставляя табличные
величины и
, находим, что
.
Общая картина энергетического спектра колебаний кристаллической решетки,
определяемого формулой (5.61), изображена на рис. 6.2. Вертикальные линии
соответствуют разрешенным значениям частот. Разрешенные значения энергии
задаются точками пересечения наклонных прямых с вертикальными, а сами наклонные
прямые соответствуют различным значениям квантовых чисел . Расстояние между точками на вертикальных прямых
равно кванту энергии
.
|
Рис. 6.2.
Зависимость энергии колебаний от частоты |
При температуре, равной нулю, в спектре присутствуют только нулевые
колебания, значения которых даны точками на наклонной прямой . Возбужденным при Т1
квантам энергии отвечают точки, расположенные ниже горизонтальной пунктирной прямой
. При температуре
возбуждаются колебания
с частотами
, а более высокочастотные колебания в спектре
отсутствуют. При дальнейшем повышении температуры возбуждаются колебания с
более высокими частотами, и одновременно растет число квантов низкочастотных
колебаний. Из рис. 6.2 следует, что при некотором значении температуры T=T*
возбуждаются все возможные колебания системы с частотами от
до
. Дальнейшее повышение температуры не приводит к появлению
волн с новыми частотами
, а ведет лишь к увеличению амплитуды колебаний (числа
возбужденных квантов) с каждой частотой
.
При повышении температуры в первую очередь возбуждаются низкочастотные
колебания. Экспериментальным путем было установлено, что минимальная частота для кристалла
размерами приблизительно 1 см составляет около 105 Гц, т. е. на
восемь порядков меньше максимальной частоты. Тогда
, т. е. низкочастотные колебания возбуждаются уже при
температуре около 10-6 K.
Макроскопические тела представляют собой совокупность очень большого числа частиц, движущихся по законам классической или квантовой механики. В таких системах свойства подчиняются статистическим закономерностям.
Найдем среднее значение энергии фонона как гармонического квантового
осциллятора. Распределение фононов по состояниям при тепловом возбуждении в гармоническом
приближении подчиняются статистике
Больцмана. В гармоническом приближении рассматривается система
невзаимодействующих фононов, т. е. ее можно представить как идеальный фононный газ. Согласно
статистике Больцмана, вероятность нахождения осциллятора в n-м квантовом состоянии с энергией равна [59]
|
(6.7) |
Коэффициент определяется из условия нормировки
. Следовательно,
. Таким образом,
|
(6.8) |
В этом случае средняя энергия осциллятора при заданной
температуре будет равна сумме произведений возможных энергий осциллятора
на их вероятность
:
|
(6.9) |
Обозначим в формуле (6.9) , тогда прямым дифференцированием можно убедиться, что
|
(6.10) |
где . Найдем величину g, подставив в (6.9)
выражение для энергии осциллятора в данном состоянии
:
,
где − бесконечно убывающая геометрическая прогрессия со
знаменателем
и первым слагаемым
. Таким образом, учитывая, что
, получим
|
(6.11) |
Отсюда
|
(6.12) |
Подставляя выражение (6.12) в (6.10), получим
|
(6.13) |
Первое слагаемое в формуле (6.13) соответствует нулевой энергии
осциллятора, а второе слагаемое можно рассматривать как произведение энергии
фонона на среднее число
фононов
, находящихся в рассматриваемом квантовом состоянии
|
(6.14) |
где .
Значит фонон – это возбуждение кристалла над нулевым уровнем энергии, соответствующим нулевым колебаниям атомов, совершающимся при температуре абсолютного нуля.
Оценим количество атомов, находящихся на возбужденных уровнях при . Поскольку каждый атом в решетке совершает одновременно
колебания со всеми частотами
, возбужденными при данной температуре, то при
на первом
энергетическом уровне (
) для атомов, колеблющихся с частотой
, энергия колебаний
. Тогда
и вероятность нахождения атома в состоянии с
данным квантовым числом n будет
. Пользуясь этим соотношением, для различных
получим
значения вероятности
:
|
|
Количество
ат. % от общего числа |
1 |
0,232 |
23% |
2 |
0,0855 |
8,5% |
3 |
0,031 |
3,1% |
Следовательно, pn
показывает, что на каждом энергетическом уровне при температуре возбуждаются не все атомы:
с энергией первого возбужденного уровня, соответствующей частоте
, колеблется только приблизительно 23 % атомов, с энергией
второго – 8,5 % и третьего – 3,1 % от общего числа атомов кристалла. Таким
образом, при температуре
значительная часть
атомов совершает только нулевые колебания с частотой
.
Число атомов, возбужденных на первом энергетическом уровне при частоте , экспоненциально быстро уменьшается с понижением
температуры. Так, при
оно составляет около
12 %; при
− около 5 % от
общего числа атомов в элементарной ячейке; при
(несколько кельвинов)
это число составляет около
(N − число атомов в ячейке), т. е. в решетке практически
отсутствуют возбужденные на частоте
атомы.
Следовательно, температуру можно рассматривать
как граничную характеристическую температуру при достижении которой в
кристалле возбуждаются колебания со всеми возможными частотами. При увеличении температуры выше
число возбужденных на
частоте
атомов быстро
возрастает, так что при
возбуждаются колебания
с частотой
практически у всех
атомов. Характеристическую температуру обычно называют температурой Дебая и обозначают
(
=
).
Используя полученную выше зависимость для средней энергии фонона (6.14), запишем выражение для среднего значения энергии тепловых колебаний всей решетки
|
(6.15) |
Поскольку нулевые колебания тепловой энергии не несут, то в (6.15) отсутствует эта часть энергии.
Для расчета необходимо знать
функцию распределения фононов по частотам
. Однако даже для простой трехмерной структуры получить
аналитическое выражение для
очень сложно. Поэтому вычисление энергии
колебаний производится для конкретных моделей, в которых вводится предположение
о виде функции
. Существуют два основных приближения: Эйнштейна (1907) и
Дебая (1912).