Температура Дебая и теплоёмкость

Температура Дебая $\theta_D$ - это характерный энергетический масштаб колебаний кристаллической решётки, выраженный в кельвинах. Она появляется в модели Дебая для теплоёмкости твёрдых тел и разделяет два предельных режима: высокотемпературный, где работает классический закон Дюлонга–Пти $C_V = 3Nk_B$, и низкотемпературный, где $C_V$ падает как $T^3$. Сама $\theta_D$ - это параметр, который однозначно определяется упругими константами и плотностью материала; для свинца она около 105 К, для меди 343 К, для алмаза 2230 К. Ниже разбираем, как из спектра фононов получается формула $C_V(T/\theta_D)$, почему именно $T^3$, и как $\theta_D$ извлекают из эксперимента.

Колебания решётки и почему классика проваливается

В классической статфизике каждая степень свободы в гармоническом осцилляторе даёт $k_B T/2$ кинетической и $k_B T/2$ потенциальной энергии. У кристалла из $N$ атомов - $3N$ нормальных мод, значит полная тепловая энергия $U = 3Nk_B T$, а молярная теплоёмкость $C_V = 3R \approx 24{,}94$ Дж/(моль·К). Это закон Дюлонга–Пти (1819), и при комнатной температуре он действительно неплохо работает для большинства металлов. Но при понижении $T$ опытная $C_V$ падает почти до нуля, тогда как классика упорно предсказывает константу. Объяснение нашлось только в квантовой механике: колебания решётки квантуются, и при низких $T$ высокочастотные моды просто не успевают возбуждаться.

Модель Эйнштейна как первый шаг

В 1907 году Эйнштейн предложил считать, что все $3N$ мод колеблются с одной и той же частотой $\omega_E$. Тогда каждая мода - квантовый осциллятор со средней энергией $\hbar\omega_E/(e^{\hbar\omega_E/k_B T}-1)$, а теплоёмкость:

$C_V = 3Nk_B \left(\frac{\theta_E}{T}\right)^2 \frac{e^{\theta_E/T}}{(e^{\theta_E/T}-1)^2},\qquad \theta_E = \hbar\omega_E/k_B.$

При $T \gg \theta_E$ формула даёт $3Nk_B$ (Дюлонг–Пти), а при $T \to 0$ - экспоненциальное падение $\sim e^{-\theta_E/T}$. Эксперимент показал, что падение действительно есть, но оно степенное, а не экспоненциальное: модель Эйнштейна качественно правильна, но количественно занижает $C_V$ при низких $T$.

Модель Дебая: спектр частот с обрезанием

Дебай (1912) заменил одночастотный спектр на сплошной - такой же, как у звуковых волн в континууме, но обрезанный сверху. Плотность мод по частоте:

$g(\omega) = \frac{9N}{\omega_D^3}\,\omega^2,\qquad 0 \leq \omega \leq \omega_D.$

Множитель $\omega^2$ берётся из подсчёта плоских волн в фазовом пространстве $\sim k^2 dk$ при линейной дисперсии $\omega = c_s k$. Обрезание $\omega_D$ выбирается так, чтобы полное число мод равнялось $3N$:

$\int_0^{\omega_D} g(\omega)\,d\omega = 3N.$

Дебаевская частота $\omega_D$ связана с температурой Дебая через $\hbar\omega_D = k_B \theta_D$. Физически $\theta_D$ - это температура, при которой возбуждены все моды решётки вплоть до самых жёстких.

Формула теплоёмкости

Подставляя $g(\omega)$ в полную энергию $U = \int_0^{\omega_D} g(\omega) \hbar\omega / (e^{\hbar\omega/k_B T} - 1)\, d\omega$ и дифференцируя по $T$, получаем:

$C_V = 9Nk_B \left(\frac{T}{\theta_D}\right)^3 \int_0^{\theta_D/T} \frac{x^4 e^x}{(e^x - 1)^2}\,dx,$

где $x = \hbar\omega/k_B T$. Это и есть формула Дебая. Интеграл $D(\theta_D/T)$ называют функцией Дебая; в общем виде он не выражается через элементарные функции, но численно табулируется и быстро считается. Все материалы укладываются на одну универсальную кривую, если по оси абсцисс отложить $T/\theta_D$.

Низкотемпературный предел: закон $T^3$

При $T \ll \theta_D$ верхний предел интеграла $\theta_D/T \to \infty$, и интеграл превращается в табличный:

$\int_0^\infty \frac{x^4 e^x}{(e^x-1)^2}\,dx = \frac{4\pi^4}{15}.$

Тогда теплоёмкость:

$C_V \approx \frac{12\pi^4}{5} Nk_B \left(\frac{T}{\theta_D}\right)^3 \approx 234\,Nk_B \left(\frac{T}{\theta_D}\right)^3.$

Это знаменитый закон $T^3$ Дебая. Физический смысл: при низкой температуре возбуждаются только длинноволновые акустические фононы; число тепловых мод в фазовом пространстве растёт как $T^3$, каждая даёт по $k_B T$, отсюда $U \propto T^4$ и $C_V \propto T^3$. Эксперимент подтверждает этот закон для изоляторов и неметаллов с впечатляющей точностью. Для металлов к фононному вкладу добавляется электронный $\gamma T$ от носителей заряда, и при очень низких $T$ зависимость $C_V/T = \gamma + \beta T^2$ даёт линейный график, по которому одновременно извлекают $\gamma$ и $\theta_D$.

Высокотемпературный предел: Дюлонг–Пти

При $T \gg \theta_D$ верхний предел интеграла $\theta_D/T \ll 1$, разложение подынтегральной функции $\frac{x^4 e^x}{(e^x-1)^2} \approx x^2$ даёт:

$\int_0^{\theta_D/T} x^2\,dx = \frac{1}{3}\left(\frac{\theta_D}{T}\right)^3,$

откуда $C_V \to 3Nk_B$. Это и есть классический закон Дюлонга–Пти. Условный порог «классики» - $T \gtrsim \theta_D$: для меди ($\theta_D = 343$ К) при комнатной температуре отклонение от $3Nk_B$ уже меньше 5%, а вот для алмаза ($\theta_D = 2230$ К) даже при $T = 300$ К теплоёмкость составляет лишь около четверти от $3Nk_B$.

Типичные значения $\theta_D$ для веществ

Чем легче атомы и жёстче связи, тем выше $\theta_D$. Несколько ориентиров:

Соотношение $\theta_D \propto c_s / a$, где $c_s$ - средняя скорость звука, $a$ - постоянная решётки. У алмаза высокая $c_s \approx 18000$ м/с и малый $a$ - отсюда экстремально высокая $\theta_D$.

Как $\theta_D$ извлекают из эксперимента

Способов несколько:

1. Из низкотемпературной теплоёмкости. Измеряют $C_V$ при $T \ll \theta_D$, строят $C_V/T^3$ как функцию $T$ (для диэлектриков) или $C_V/T$ как функцию $T^2$ (для металлов, чтобы отделить электронный вклад) и из наклона извлекают $\theta_D$. Это самый прямой метод.
2. Из скорости звука. $\theta_D = \frac{\hbar}{k_B}(6\pi^2 n)^{1/3} \bar c_s$, где $n$ - концентрация атомов, $\bar c_s$ - усреднённая по поляризациям скорость звука. Совпадение «упругого» и «калориметрического» значений - хороший тест применимости модели.
3. Из упругих констант. Для изотропного кристалла $\bar c_s$ вычисляется через модули объёмной и сдвиговой упругости.
4. Из фактора Дебая–Валлера в рентгеновской дифракции - амплитуда тепловых смещений атомов прямо связана с $\theta_D$.

Ограничения модели и реальные спектры

Дебаевская плотность $g(\omega) \propto \omega^2$ - это приближение для длинноволновых акустических фононов. Реальный спектр сложнее: помимо акустических ветвей есть оптические (для решёток с базисом из нескольких атомов), а на краях зоны Бриллюэна $\omega(k)$ загибается и появляются ван-Хововы сингулярности. Поэтому единое $\theta_D$ - фикция: подгоночное значение медленно «плывёт» с температурой. Часто различают $\theta_D^\text{th}$ (из теплоёмкости при разных $T$), $\theta_D^\text{el}$ (из упругости) и $\theta_D^\text{DW}$ (из дифракции) - они близки, но не совпадают точно. Для двумерных систем (графен) закон становится $C_V \propto T^2$, для одномерных - $C_V \propto T$. В сильно анизотропных слоистых кристаллах (висмут, графит) дебаевская модель плохо работает в промежуточной области и её часто заменяют комбинацией двух «эффективных» $\theta_D$.

Связь с другими моделями

В пределе $\theta_D/\theta_E$ → высокая модель Эйнштейна - это частный случай модели Дебая, в котором весь спектр сжимается в одну линию. Современная теория теплоёмкости использует прямой счёт $g(\omega)$ из ab initio расчётов фононных дисперсий (DFT-DFPT) и интегрирует $C_V = \int g(\omega) \cdot c_\text{осц}(\hbar\omega/k_B T)\,d\omega$ численно - без всяких подгоночных $\theta_D$. Тем не менее в инженерных задачах (термоизоляция, аккумуляция тепла, тепловые свойства керамик) модель Дебая по-прежнему удобный язык: одно число $\theta_D$ компактно описывает всю температурную зависимость.

Частые ошибки

• Путать $\theta_D$ и $\theta_E$. Эйнштейн - одна частота, экспоненциальное падение; Дебай - спектр $\omega^2$ с обрезанием, степенное падение $T^3$. Численно $\theta_D$ и $\theta_E$ для одного материала отличаются: обычно $\theta_E \approx 0{,}75\,\theta_D$.
• Применять закон $T^3$ во всём диапазоне. Он работает только при $T \lesssim \theta_D/10$. Между $\theta_D/10$ и $\theta_D$ обязательно интегрировать функцию Дебая, без сокращений.
• Забывать про электронный вклад в металлах. При $T \to 0$ у металлов $C_V = \gamma T + \beta T^3$. Если интерпретировать всё как фононный вклад, $\theta_D$ получится заниженной.
• Использовать $\theta_D$ из калориметрии в формулах для звука. Они близки, но не тождественны; для прецизионных расчётов берите $\theta_D$ из соответствующего метода.
• Применять модель Дебая к стёклам и аморфным телам без оговорок. У аморфных систем при $T \lesssim 1$ К появляется аномальный линейный по $T$ вклад от двухуровневых систем - модель Дебая его не описывает.

FAQ

Почему именно $T^3$, а не $T^2$ или $T^4$? Показатель степени равен размерности пространства акустических мод: длинноволновые фононы образуют 3D-«газ» в обратном пространстве. Число тепловых мод в шаре радиуса $k_T \propto T$ растёт как $T^3$, каждая даёт $k_B T$ энергии - отсюда $U \propto T^4$ и $C_V \propto T^3$. Для 2D-кристаллов аналогично получается $T^2$.

Чем модель Дебая лучше модели Эйнштейна? Эйнштейн при $T \to 0$ даёт экспоненциальное падение теплоёмкости, что противоречит эксперименту. Дебай учитывает длинноволновые акустические моды, которые остаются возбуждёнными даже при сколь угодно низких $T$, и даёт правильный закон $T^3$. Зато для оптических ветвей и узких пиков в спектре модель Эйнштейна нередко точнее.

Как связаны $\theta_D$ и скорость звука? $\theta_D = (\hbar/k_B)(6\pi^2 n)^{1/3} \bar c_s$. Чем выше скорость звука и плотность атомов, тем выше $\theta_D$. Поэтому алмаз с его жёсткими ковалентными связями и лёгкими атомами углерода - рекордсмен по $\theta_D$ среди объёмных материалов.

Коротко

Температура Дебая $\theta_D$ - единый параметр, описывающий теплоёмкость кристалла в рамках модели Дебая. Формула $C_V = 9Nk_B (T/\theta_D)^3 \int_0^{\theta_D/T} x^4 e^x/(e^x-1)^2\,dx$ даёт классический предел Дюлонга–Пти $C_V = 3Nk_B$ при $T \gg \theta_D$ и знаменитый закон $T^3$ Дебая при $T \ll \theta_D$, выводимый из спектра $g(\omega) \propto \omega^2$ с обрезанием на $\omega_D = k_B\theta_D/\hbar$. Типичные значения: 105 К у свинца, 343 К у меди, 2230 К у алмаза - чем легче атомы и жёстче связи, тем выше $\theta_D$. Модель Дебая остаётся лучшим аналитическим описанием фононной теплоёмкости и языком, на котором обсуждают тепловые свойства твёрдых тел.