Калашников, Н.П. Руководство к решению задач по физике Основы квантовой физики. Строение вещества. Атомная и ядерная физика
.pdfРис. 8
Число Z элементарных ячеек в единице объема кристалла Z = Zm / Vm или в общем случае
Z = ρ kn NμA ;
для кристалла, состоящего из одинаковых атомов (k =1),
Z= ρ NnμA .
●Параметр а кубической решетки
a = 3 nμ/ (kρNA ) .
Расстояние, d между соседними атомами в кубической решетке;
а) в гранецентрированной d = a / |
2 ; |
б) в объемно центрированной d = |
3 a / 2 . |
● Вращательная энергия двухатомной молекулы:
EJ = 2I2 J (J +1),
гдеI – моментинерциимолекулы, J – квантовоечисло(J = 0, 1, 2, …). ● Колебательная энергия двухатомной молекулы:
Ev = ω(v +1/ 2),
131
где ω – собственная частота колебаний молекулы, v = 0, 1, 2, … – квантовое число.
•Средняя энергия квантового гармонического осциллятора при температуреT :
ε = |
ω |
+ |
|
ω |
|
, |
(3.1.1) |
2 |
e |
ω/k T |
−1 |
||||
|
|
Б |
|
|
где ω – круговая частота колебаний осциллятора, kБ – постоянная Больцмана, T – термодинамическая температура.
• Внутренняя молярная энергия кристалла (энергия тепловых колебаний кристаллической решетки) в квантовой теории теплоемкости Эйнштейна:
|
|
U =U0 +3R |
θE |
|
|
|
|
|
, |
(3.1.2) |
|
|
|
eθE /T −1 |
|||
где |
θE = ω/ kБ |
– характеристическая температура Эйнштейна; |
|||
U0 |
= 3RθE / 2 – |
молярная энергия нулевых |
колебаний; R – |
||
универсальная газовая постоянная. |
|
• Молярная теплоемкость кристалла в квантовой теории теплоемкости Эйнштейна:
|
dU |
|
|
|
θ |
|
2 |
|
|
eθE /T |
|
|
|
||
cp = |
|
= 3R |
|
|
E |
|
|
|
|
. |
(3.1.3) |
||||
|
|
|
( |
|
) |
2 |
|||||||||
|
dT |
|
|
T |
|
|
|
|
|
||||||
|
|
|
|
eθE /T −1 |
|
|
|
||||||||
При высоких температурах (T |
|
|
θE ): |
|
|
|
|
|
|
||||||
|
|
|
|
1 |
θ |
|
|
2 |
|
|
|
|
|||
cp |
= 3R 1 |
− |
|
|
|
|
E |
|
|
+... . |
|
|
(3.1.4) |
||
|
|
|
|
|
|
||||||||||
|
|
|
|
12 |
T |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В пределе бесконечно высоких температур cp →3R (закон Дюлонга и Пти).
132
При низких температурах (T θE ): |
|
|
|
|
|
|||
cp |
|
θ |
2 |
−θ |
|
/T |
. |
|
= 3R |
|
E e |
|
E |
|
(3.1.5) |
||
|
|
T |
|
|
|
|
||
• В квантовой теории |
теплоемкости |
|
Дебая число dN |
собственных колебаний кристаллической решетки, приходящихся на интервал частот от ω до ω+ dω, определяется выражением:
dN = 9N |
ω2dω |
, |
(3.1.6) |
|
ω3 |
|
|
|
max |
|
|
где N — число атомов, ωmax — максимальная частота, ограничивающая спектр колебаний (т.н. дебаевская частота):
ωmax = v 3 6π2n . |
(3.1.7) |
Здесь n = N /V – концентрация атомов, v – средняя скорость волны колебаний в кристалле.
•Скорости продольных (v ) и поперечных (v ) звуковых волн
вкристалле определяются по формулам
v = |
E / ρ , v = |
G / ρ , |
(3.1.8) |
||||||||||
где ρ – плотность среды, |
а E и |
G |
– модули продольной и |
||||||||||
поперечной упругости, соответственно. Средняя |
скорость v |
||||||||||||
звуковой волны связана с v |
и v |
соотношением |
|
||||||||||
|
|
3 |
= |
|
2 |
|
+ |
1 |
. |
|
(3.1.9) |
||
|
|
v2 |
|
|
|
|
|
||||||
|
|
|
|
v2 |
v2 |
|
|
||||||
• Наименьшая длина волны, возбуждаемая в кристалле, |
|||||||||||||
λ |
|
= |
|
2πv |
≈ |
2 |
|
|
≈ 2d , |
(3.1.10) |
|||
min |
|
|
3 n |
||||||||||
|
|
|
ωmax |
|
|
где d – расстояние между соседними атомами в решетке.
• Молярная внутренняя энергия кристалла по Дебаю:
133
ω |
|
|
|
T |
|
3 |
θ |
/T |
|
x3 |
|
|
max |
|
dN =U0 |
|
D |
|
|
|
|||||
U = ∫ |
ε |
+9RT |
|
|
|
∫ |
e |
x |
|
dx . (3.1.11) |
||
0 |
|
|
|
θD |
|
0 |
|
−1 |
|
Здесь U0 = 98 RθD — молярная нулевая энергия кристалла по
Дебаю; θD = ωmax / kБ — характеристическая температура Дебая.
• Молярная теплоемкость кристаллической решетки по Дебаю:
|
|
|
|
T |
|
3 |
θD /T |
|
c |
|
= 3R 12 |
|
∫ |
||||
p |
|
θ |
|
|
|
|||
|
|
|
|
|||||
|
|
|
|
|
D |
|
0 |
При низких температурах (T
cp ≈ 125
|
x3 |
|
|
|
|
|
|
θD / T |
|
|
|
|
|
|
|
|
|
dx − |
3 |
|
|
. |
|
|
ex −1 |
eθD /T − |
|
||||||||
|
|
|
|
1 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
θD ): |
|
|
|
|
|
|
|
|
||
4 |
R |
|
|
|
T 3 |
|
|
|
|
||
π |
|
|
|
|
, |
|
|
|
|||
|
|
|
|
|
|
||||||
|
|
|
|
θD |
|
|
|
|
при высоких температурах (T θD ):
(3.1.12)
(3.1.13)
|
|
|
1 |
θ |
|
2 |
1 |
θ |
D |
4 |
|
||
cp |
≈ 3R 1 |
− |
|
|
|
D + |
|
|
|
|
+... . (3.1.14) |
||
|
|
|
|
|
|||||||||
|
|
|
20 |
T |
|
560 |
T |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В пределе бесконечно высоких температур cp |
стремится к 3R |
(закон Дюлонга и Пти). Сравнение приближенных формул с точными показано на рис. 9.
• Результаты для тепловых свойств кристалла выписаны в предположении, что его молекулы состоят из одного атома. Если в
молекулу входят na атомов, то всерезультаты надо умножить на na .
• Энергия ε и квазиимпульс p фонона связаны с круговой частотой ω колебаний обычными соотношениями
ε = ω, p = |
2π |
. |
(3.1.15) |
|
|||
|
λ |
|
134
Рис. 9 Молярная теплоемкость кристалла по Дебаю (сплошная линия)
в сравнении с приближенными формулами (3.1.13) (штриховая линия) и (3.1.14) (точечная линия). Видно, что первую формулу
можно применять с высокой степенью точности при T ≤ 0,1 θD , а вторую T ≥ 0,5 θD . Тонкой сплошной линией показана теплоемкость, рассчитанная в рамках теории Эйнштейна при
θ |
E |
= |
3 |
θ |
D |
(см. пример 3.1.1). Видно, что при T > 0,3 θ |
|
обе |
|
5 |
D |
||||||||
|
|
|
|
|
теории практически неразличимы.
3.1.2. Методические рекомендации по решению задач
Пример 3.1.1. Найти связь характеристических температур Эйнштейна и Дебая.
Решение. Как известно, в теории теплоемкости Дебая характеристическая температура выражается через параметры материала, в то время как в теории Эйнштейна характеристическая температура должна быть выбрана феноменологически, путем сравнения теоретической формулы с данными эксперимента. В то же время обе теории должны давать более или менее сходные результаты, по крайней мере, в области высоких температур, где
135
значения теплоемкости неплохо описываются законом Дюлонга и Пти.
Поэтому мы приравняем первые два члена высокотемпературных разложений (3.1.4) и (3.1.14) в обеих теориях:
|
|
|
|
|
|
|
1 |
θ |
|
2 |
|
|
|
|
|
|
|
|
1 |
θ |
D |
2 |
|
|
|
||||||
|
|
|
3R |
1 |
− |
|
|
|
|
E |
|
+... |
= 3R 1 |
− |
|
|
|
|
|
|
+... , |
|
(3.1.16) |
||||||||
|
|
|
12 |
|
|
20 |
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
откуда находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
θ |
E |
= |
|
3 |
θ |
D |
= |
0,775θ |
D |
. |
|
|
|
|
|
(3.1.17) |
||||
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Другой возможный путь установления связи между двумя |
|||||||||||||||||||||||||||||
теориями |
|
|
|
– |
|
|
|
|
приравнивание |
|
|
|
|
нулевых |
энергий: |
||||||||||||||||
(3 / 2)RθE = (9 / 8)RθD , |
откуда |
следует |
|
|
близкий |
результат: |
|||||||||||||||||||||||||
θ |
E |
= 3 |
θ |
D |
= |
0, 75θ |
D |
. Результаты теорий Эйнштейна и Дебая при |
|||||||||||||||||||||||
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
связи температур (3.1.17) можно сравнить на рис. 9. |
|
|
|||||||||||||||||||||||||||||
|
|
Пример |
3.1.2. |
Определить |
количество |
теплоты |
Q , |
необходимое для нагревания кристалла NaCl массой m = 20 г на T = 2 К, если нагревание происходит от температуры:
1) |
T1 = θD ; 2) T2 = 2 |
К. Температуру |
Дебая θD |
для NaCl |
||||
принять равной 320 К. |
|
|
|
Q подводимое для нагревания |
||||
Решение. Количество теплоты |
||||||||
тела от температуры Ti |
до Ti +ΔT , может быть вычислено по |
|||||||
формуле |
|
|
|
Ti +ΔT |
|
|
|
|
|
|
|
μ |
|
|
|
||
|
|
|
T∫ |
|
|
|
||
|
|
Q = |
m |
|
|
cp dT , |
(3.1.18) |
|
|
|
|
i |
|||||
|
|
|
|
|
|
|
|
|
где |
μ — молярная масса вещества, а |
cp — его |
молярная |
теплоемкость, которая в общем случае есть функция температуры, поэтому за знак интеграла ее выносить нельзя. Однако в первом случае изменение температуры мало по сравнению с ее начальным значением, так что теплоемкость можно считать постоянной:
136
cp (T )≈ cp (θD ). Для вычисления последней можно применить приближенную формулу (3.1.14), откуда находим
|
|
1 |
|
1 |
|
|
|
|
cp (θD )≈ 3R 1 |
− |
|
+ |
|
|
= 2,855R . |
(3.1.19) |
|
20 |
560 |
|||||||
|
|
|
|
|
|
Подставляя это значение в формулу (3.1.19), получим Подставляя это значение в формулу (3.1.19), получим
Q = 2,855R |
m |
T = 2,855 ×8,31× |
20 |
× 2 =16, 2 Дж. (3.1.20) |
|
μ |
23,0 + 35,5 |
||||
|
|
|
|||
Во втором случае нахождение |
Q облегчается тем, что можно |
воспользоваться предельным законом Дебая (3.1.13), в согласии с которым теплоемкость пропорциональна кубу абсолютной температуры. В этом случае теплоемкость сильно изменяется в пределах заданного интервала температур и ее нельзя выносить за знак интеграла в формуле (3.1.19). Получаем в этом случае
|
12m |
4 |
|
T2 +ΔT |
T 3 |
|
12m |
4 |
|
|
(T2 +ΔT )4 −T24 |
|
|||||
Q = |
|
π |
R |
|
|
|
dT = |
|
|
π |
R |
|
|
3 |
. (3.1.21) |
||
|
|
|
|
|
|
||||||||||||
|
5μ |
|
|
|
T∫ |
θD |
|
|
5μ |
|
|
|
4 |
θD |
|
||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя численные значения, находим |
|
|
|||||||||||||||
|
|
Q = |
12×20 |
π |
4 |
×8,31 |
44 −24 |
=1, 22 мДж. |
|||||||||
|
|
|
5(23,0 +35,5) |
|
4×3203 |
|
Пример 3.1.3. Вычислить дебаевскую частоту и дебаевскую температуру для железа, если скорости распространения продольных и поперечных колебаний равны 5,85 и 3,23 км/с соответственно.
Решение. Во-первых, найдем концентрацию атомов железа, используя таблицы в последнем разделе. Плотность железа
ρ = 7,87 106 |
г/м3. Молярная масса железа приведена в таблице |
Менделеева: |
μ =55,85 г/моль. Отсюда находим сначала |
молярный объем: Vm = μ / ρ = 7,10 10−6 м3/моль.
137
Теперь вычисляем концентрацию атомов в кристалле железа: n = NA /Vm =8, 48 1028 м-3.
Во-вторых, определим среднюю скорость звука в железе, используя соотношение (3.1.9): 3 / v2 = 2 / 32302 +1/ 58502 , откуда v = 3685 м/с.
Дебаевскую частоту находим теперь по формуле (3.1.7):
ωmax = 36853 6π2 ×8, 48 1028 = 6,31 1013 Гц.
И, наконец, вычисляем дебаевскую температуру:
θD = ωmax / kБ = 482 К.
Пример 3.1.4. Оценить давление фононов в меди при температуре T = θD , если θD = 320 К. Фононы рассматривать
как идеальный газ.
Решение. Для определения давления фононов воспользуемся уравнением Менделеева-Клапейрона, записанным в форме
pV = N |
p v |
= N |
pv |
, |
(3.1.22) |
|
3 |
|
3 |
|
|
где p и v — импульс и скорость фононов, а N |
— их число в |
системе объемом V . Мы заменили скалярное произведение на обычное, так как импульс фонона параллелен его скорости. В качестве v мы выбираем усредненную скорость звука, а импульс
фонона связан с его энергией ε |
соотношением |
p = ε/ v , так что |
||||
получаем |
N ε |
|
|
|
|
|
p = |
= |
W |
, |
(3.1.23) |
||
3V |
3V |
|||||
|
|
|
|
где W — полная кинетическая энергия фононов. Рассмотрим моль меди. Ее молярный объем находим с
использованием таблиц в последнем разделе, разделив молярную массу μ на плотность ρ :
V = |
μ |
= |
63,55 10−3 |
= 7,12 10−6 |
м3/моль. (3.1.24) |
|
ρ |
8,93 103 |
|||||
m |
|
|
|
138
Полную кинетическую энергию W фононов можно найти: 1) вычитая из молярной энергии U кристалла энергию U0
нулевых колебаний и 2) уменьшая результат в два раза (так как на кинетическую энергию фононов приходится ровно половина
полной энергии): W = (U −U0 )/ 2 . Получаем тогда
|
p = |
U −U0 |
. |
|
(3.1.25) |
||
|
|
|
|
||||
|
|
6V |
|
|
|
|
|
|
|
|
m |
|
|
|
|
При T = θD находим из (3.1.11): |
|
|
|
|
|||
1 |
x3dx |
|
|
|
|
||
U −U0 = 9RθD ∫0 |
|
= 9 |
×8,31×320 |
×0,0778 |
= |
||
exp(x)−1 |
|||||||
=1,86 103 Дж/моль |
|
(3.1.26) |
(интеграл в правой части взят численно). Находим окончательно
p = |
|
1,86 103 |
= 43,5 |
МПа. |
(3.1.27) |
|
6 |
×7,12 10−6 |
|||||
|
|
|
|
Пример 3.1.5. Определить среднее число фононов в моде упругих колебаний кристаллической решетки, для которой спектральная плотность числа фононов максимальна. Считать
выполненным условие T < 0,5θD .
Решение. Определим частоту моды, соответствующей максимальной спектральной плотности числа фононов. Функцию
спектральной плотности числа фононов f (ω) получим, перемножая спектральную плотность фононных мод D (ω) на среднее число фононов в моде < n >:
f (ω)= D(ω)< n >= |
3Vω2 |
|
1 |
. |
(3.1.28) |
|
|
||||
|
2π2vзв3 e ω/kБT −1 |
|
Для упрощения последующих преобразований произведем замену k ωT = x . После этого имеем
Б
139
|
|
f (x) |
|
3V (kБT )2 |
x2 |
x2 |
|
|
|
|||||||||
|
|
= |
|
|
|
|
|
|
= a |
|
, |
(3.1.29) |
||||||
|
|
2π2vзв3 |
2 |
|
ex |
−1 |
ex −1 |
|||||||||||
где a = |
3V (kБT )2 |
. При x → 0 |
и |
x →∞ эта |
неотрицательная |
|||||||||||||
|
|
|||||||||||||||||
|
2π2v3 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
зв |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
функция |
стремится |
к |
|
нулю. |
Следовательно, |
в промежутке |
||||||||||||
0 ≤ x < +∞ функция |
|
f (x) |
имеет максимум. |
Дифференцируя |
||||||||||||||
f (x) по x и приравнивая результат к нулю, получим |
|
|
|
|||||||||||||||
|
|
x (2ex − 2 − xex ) |
= 0 |
или |
x = 2 1−e−x . |
(3.1.30) |
||||||||||||
|
( |
|
|
) |
|
|||||||||||||
|
|
|
2 |
|
|
|
|
( |
|
|
) |
|
||||||
|
|
|
ex −1 |
|
|
|
|
|
|
|
|
Решая полученное трансцендентное уравнение методом последовательных приближений, находим xmax =1,59 , т.е. частота моды, соответствующей максимуму функции спектральной
плотности числа фононов f (ω): |
|
ωmax =1,59 kБT . |
(3.1.31) |
Тогда искомое среднее число фононов в этой моде вычисляется согласно соотношению
1 |
|
|
< n >= e ω/ KБT −1 |
= 0, 26 |
(3.1.32) |
Пример 3.1.6. Атом может находиться в любом из двух квантовых состояний с разностью энергий Е. Построить качественную зависимость удельной теплоемкости С ансамбля таких атомов от абсолютной температуры Т. Рассмотреть поведение теплоемкости в предельных случаях T → 0 и T →∞.
− |
E |
|
|
kБT , N1 N и N1 + N2 = N , |
|||
Решение. Поскольку N2 Ne |
|||
то, мы имеем |
|
|
140