Теперь обратимся к вычислению названных интегралов. Нетрудно видеть, что эти интегралы можно записать и форме
J =  (5.17)
                                                     (5.17)           
где  — некоторая произвольная функция энергии. Проинтегрируем (5.17) по
частям:
— некоторая произвольная функция энергии. Проинтегрируем (5.17) по
частям:
J =  (5.18)
                                            (5.18)
Здесь  = g(ε),
 = g(ε),  .
.
Видно, что в (5.18) первый член
обращается в нуль, т.к. g(0)=0
и G(ε) g(ε) → 0 при ε →∞. Второй член в (5.18)
может быть вычислен благодаря следующему обстоятельству. Как отмечено выше, в
металлах  . Поэтому
при T ≠ 0 функция дf/дε имеет симметричный пик пиком шириной ~
. Поэтому
при T ≠ 0 функция дf/дε имеет симметричный пик пиком шириной ~ при
 при  (см. рис. 5.2), а при T = 0 ведет себя как δ-функция δ(ε -
(см. рис. 5.2), а при T = 0 ведет себя как δ-функция δ(ε -  ). При этом G можно разложить в ряд Тейлора вблизи
). При этом G можно разложить в ряд Тейлора вблизи  :
:
 
  
и представить интеграл J в виде ряда
J = - =
 = 
=-  (5.19)
                               (5.19)
Входящие в J интегралы можно привести к безразмерному виду
 =
 = =
=
= (kBT)n =
 =  .
.
Поскольку εF/kBT >> 0 и In ≈ , интеграл (5.19) можно представить в
виде ряда
, интеграл (5.19) можно представить в
виде ряда
J = -
-  (kBT)n In, In =
(kBT)n In, In = 
 
Для четных n подынтегральная функция антисимметрична и эти интегралы равны нулю I2n-1 = 0. Используя это представление для интеграла I2n , получим
I0 =  , I2n≠ 0, n = 0,
1, 2.   , I2 = π2/5.                 (5.20)
, I2n≠ 0, n = 0,
1, 2.   , I2 = π2/5.                 (5.20)
Из (5.89) следует
 .                                           
(5.21)
.                                           
(5.21)
(1/2) I2 =(1/2)
I2 =(1/2) I2 =(
I2 =( )εF
)εF  (5.22)
                                  (5.22)
Следовательно, для любой функции от энергии g(ε) с точностью до отношения (kBT/εF)2<< 1, получаем
 + O((kBT)4).                       (5.23)
+ O((kBT)4).                       (5.23)
Воспользуемся этим методом для расчета всех термодинамических функций. Сначала рассмотрим зависимость химического потенциала от температуры, т.е. зависимость μ(Т). Число электронных состояний, с использованием (5.22) можно записать как
 (5.24)
                        (5.24)
где  ,
,  (T)
- зависящая от температуры добавка к химическому потенциалу электронов при T = 0, равному
(T)
- зависящая от температуры добавка к химическому потенциалу электронов при T = 0, равному  
 
Учтем, что число электронов в кристалле N0 не зависит от T и равно:
 
                                     
Тогда, преобразуя (5.24), получаем  N0 =  =
 =
 (5.25)
 (5.25)
Из (5.25)
получаем зависящую от температуры добавку к химическому потенциалу  (T),
(T), 
 .                             
(5.26)
.                             
(5.26)
Оценим ее отношение к химическому
потенциалу электронов при T
= 0, равному  , поскольку
, поскольку  ν(
ν( )~
)~ 1/2,
1/2,  ~
~ .
.
 (5.27)
       (5.27)
Окончательно зависимость химического потенциала от температуры имеет вид
 ,=
,=  (5.28)
                          (5.28)
Вычислим внутреннюю энергию электронов при T ≠ 0, используя соотношение (5.22):
 =
= 
 .                                      
(5.29)
 .                                      
(5.29)
Зависимость внутренней энергии электронов от температуры имеет вид
 .                                            
(5.30)
.                                            
(5.30) 
Теплоемкость при постоянном объеме равна
 .                                                (5.31)
.                                                (5.31)
Обсудим получившееся выражение для
теплоемкости электронов в кристалле. Физический смысл полученного соотношения
ясен — только электроны вблизи энергии Ферми могут принимать участие в тепловом
возбуждении системы. Число таких электронов, как говорилось выше, порядка  , умноженное на плотность
состояний вблизи
, умноженное на плотность
состояний вблизи  . Каждый
электрон вносит в (5.102) вклад
. Каждый
электрон вносит в (5.102) вклад  . Важно отметить, что при получении (5. 102) нам не нужны были никакие данные о функции
плотности состояний
. Важно отметить, что при получении (5. 102) нам не нужны были никакие данные о функции
плотности состояний  .
Следовательно, теплоемкость для всех металлов является линейной функцией
температуры. Понятно также, что при малых температурах основной вклад в
теплоемкость твердого тела дают электроны, поскольку cVe~ Т, а для решеточной теплоемкости cVp~ Т5.
.
Следовательно, теплоемкость для всех металлов является линейной функцией
температуры. Понятно также, что при малых температурах основной вклад в
теплоемкость твердого тела дают электроны, поскольку cVe~ Т, а для решеточной теплоемкости cVp~ Т5. 
Рассмотрим теперь соотношение между cVe и cVpпри высоких температурах (Т>ΘD). Получим сначала для оценок формулу теплоемкости для газа свободных электронов. В этом случае, как нетрудно показать, химический потенциал мало отличается от энергии Ферми
 ≈
 ≈ .                                       
Теплоемкость электронной подсистемы растет линейно
.                                       
Теплоемкость электронной подсистемы растет линейно
 (5.32)
              (5.32)
Подставив в (5.32) характерные
значения пeи
 (например при пe ~1022 см-3,
 (например при пe ~1022 см-3,  7 эВ, T~0.03 эВ), получим
7 эВ, T~0.03 эВ), получим
 (5.33)
                                  (5.33)
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.