固体物理晶格振动与晶体的热力学函数
第三章 晶格振动与晶体的热力学函数
一、填空体
1. 若在三维空间中，晶体由N 个原胞组成,每个原胞有一个原子,则共有_ 3 N_个独立的 振动,_ N__个波矢, 3N_支格波。
2. 体积为V 的ZnS 晶体，如果晶胞的体积为Ω，则晶格振动的模式书为24N/Ω 。
3. 三维绝缘体晶体的低温比热Cv 与温度T 的关系为Cv~T 3。
4. 某三维晶体由N 个原胞组成，每个原胞内有3个原子。考虑晶体的晶格振动，其色散关系共有 9N 支，其中 3N 支声学波，包括 2N 支横声学波， 1N 支纵声学波；另有 6N
π
2L 。
二、基本概念 1. 声子
晶格振动的能量量子。
2.波恩-卡门条件
即周期性边界条件，设想在实际晶体外，仍然有无限多个相同的晶体相连接，各晶体中相对应的原子的运动情况都一样。 3.波矢密度
波矢空间单位体积内的波矢数目，三维时为
3
c
)2(V ，Vc 为晶体体积。 4. 模式密度
单位频率间隔内模式数目。 5.晶格振动。
答：由于晶体内原子间存在着相互作用，原子的振动就不是孤立的，而要以波的形式在晶体中传播，形成所谓格波，因此晶体可视为一个互相耦合的振动系统，这个系统的运动就叫晶格振动。 6.简谐近似 答：当原子在平衡位置附近作微小振动时，原子间的相互作用可以视为与位移成正比的虎克力，由此得出原子在其平衡位置附近做简谐振动。这个近似即称为简谐近似。 7.格波
答：晶格中的原子振动是以角频率为ω的平面波形式存在的，这种波就叫格波。 三、简答题
1. 试分析爱因斯坦模型和德拜模型的特点及局限性. 特点：
1）爱因斯坦模型假设晶体中所有原子都以相同的频率作振动； 2）德拜模型的基本思想是把格波作为弹性波来处理。 局限性：
1） 在爱因斯坦的假设下，解释了在甚低温时温度的变化趋势，但是不能解释为什么晶体热
熔随温度T 3的速度变化，这是因为，爱因斯坦模型只考虑了光学支格波，忽略了声学支格波，而在甚低温决定晶体热容的主要是长声学波。爱因斯坦模型过于简化。
2） 德拜模型不仅能够很好解释在甚低温时晶体热容随温度的变化趋势，同时得出了在甚低
温下，热容与T 3成正比的规律。但是德拜模型忽略了晶体的各向异性，即光学波和高频声学波对热容的贡献。
2. 长光学支格波与长声学支格波本质上有何差别?
答：长光学支格波的特征是每个原胞内的不同原子做相对振动, 振动频率较高, 它包含了晶格振动频率最高的振动模式. 长声学支格波的特征是原胞内的不同原子没有相对位移, 原胞做整体运动, 振动频率较低, 它包含了晶格振动频率最低的振动模式, 波速是一常数. 任何晶体都存在声学支格波, 但简单晶格(非复式格子)晶体不存在光学支格波.
3. 晶体中声子数目是否守恒? 答：频率为
的格波的(平均) 声子数为
,
即每一个格波的声子数都与温度有关, 因此, 晶体中声子数目不守恒, 它是温度的变量.
4. 温度一定，一个光学波的声子数目多呢, 还是声学波的声子数目多? 答：频率为 的格波的(平均) 声子数为
.
因为光学波的频率
比声学波的频率
高, (
)大于(
), 所以在
温度一定情况下, 一个光学波的声子数目少于一个声学波的声子数目.
5. 对同一个振动模式, 温度高时的声子数目多呢, 还是温度低时的声子数目多? 答：设温度TH>TL, 由于(
)小于(
), 所以温度高时的声子数目多
于温度低时的声子数目.
6. 高温时, 频率为 的格波的声子数目与温度有何关系? 答：温度很高时,
, 频率为
的格波的(平均) 声子数为
.
可见高温时, 格波的声子数目与温度近似成正比.
7. 长声学格波能否导致离子晶体的宏观极化? 答：长光学格波所以能导致离子晶体的宏观极化, 其根源是长光学格波使得原胞内不同的原子(正负离子)产生了相对位移. 长声学格波的特点是, 原胞内所有的原子没有相对位移. 因此, 长声学格波不能导致离子晶体的宏观极化.
8. 试定性给出一维单原子链中振动格波的相速度和群速度。 答：由一维单原子链的色散关系2
sin
2
qa
m
β
ω= 可求得一维单原子链中振动格波的相速度为2
/2sin
qa
qa m
a
q
p β
ω
υ== 群速度为
2cos qa
m qa dq d g βωυ==
9. 周期性边界条件的物理含义是什么？引入这个条件后导致什么结果？如果晶体是无限大，q 的取值将会怎样？
答：由于实际晶体的大小总是有限的，总存在边界，而显然边界上原子所处的环境与体内原子的不同，从而造成边界处原子的振动状态应该和内部原子有所差别。考虑到边界对内部原子振动状态的影响，波恩和卡门引入了周期性边界条件。其具体含义是设想在一长为Na 的
有限晶体边界之外，仍然有无穷多个相同的晶体，并且各块晶体内相对应的原子的运动情况一样，即第j个原子和第Nt+j个原子的运动情况一样，其中t ＝1，2，3…。
q只能取一些分立的不同值。如果晶体是无引入这个条件后，导致描写晶格振动状态的波矢
q的取值将趋于连续。
限大，波矢
10.下图表示一维双原子复式晶格振动的两支格波的色散关系。请简要分析并判断：在长波极限下，图中哪一条曲线反映了初基元胞内两个原子的质心振动？图中哪一条曲线反映了初基元胞内两个原子的相对振动？
答:
上半部分曲线表示光学支，光学支格波反映了晶体中分子内两个原子的相对振动；下半部分曲线表示声学支，声学支格波反映了晶体中分子的质心振动。
由N个原胞所组成的复式三维晶格，每个原胞内有r个原子，试问晶格振动时能得到多少支色散关系？其波矢的取值数和模式的取值数各为多少？
答：共有3r支色散关系，波矢取值数=原胞数N，模式取值数=晶体的总自由度数。
11.对于初基晶胞数为N的二维晶体，基元含有四个原子，声学支震动模式和光学支震动模式的数目各为多少？
答：2N，6N。
12.在三维晶体中，格波独立的点数N，格波个数，格波总支数，声学波支数分别等于多少？答：在三维晶格中，格波独立的点数是，格波个数有3Nn，格波总支数是3nN，对每个波矢q，有3支声学波，（3n-3）支光学波。
13.试述长光学波与长声学波的本质区别?
答：长光学支格波的特征是每个原胞内的不同原子做相对振动, 振动频率较高, 它包含了晶格振动频率最高的振动模式。长声学支格波的特征是原胞内的不同原子没有相对位移, 原胞做整体运动, 振动频率较低, 它包含了晶格振动频率最低的振动模式, 波速是一常数。任何晶体都存在声学支格波, 但简单晶格(非复式格子)晶体不存在光学支格波。
14. 长声学格波能否导致离子晶体的宏观极化?
答：长光学格波所以能导致离子晶体的宏观极化，其根源是长光学格波使得原胞内不同的原子(正负离子)产生了相对位移。长声学格波的特点是, 原胞内所有的原子没有相对位移. 因此，长声学格波不能导致离子晶体的宏观极化。
15.爱因斯坦模型在低温下与实验存在偏差的根源是什么?
1013, 属于答:按照爱因斯坦温度的定义, 爱因斯坦模型的格波的频率大约为Hz 光学支频率. 但光学格波在低温时对热容的贡献非常小, 低温下对热容贡献
大的主要是长声学格波. 也就是说爱因斯坦没考虑声学波对热容的贡献是爱因斯坦模型在低温下与实验存在偏差的根源。
16. 在甚低温下, 德拜模型为什么与实验相符? 答：在甚低温下, 不仅光学波得不到激发, 而且声子能量较大的短声学格波也未被激发,得到激发的只是声子能量较小的长声学格波. 长声学格波即弹性波. 德拜模型只考虑弹性波对热容的贡献. 因此, 在甚低温下, 德拜模型与事实相符, 自然与实验相符。
四、证明计算
1. 证明一维单原子链的运动方程，在长波近似下，可以化成弹性波方程，
)
2()(2222l
n l n u q a t
u m ++-=??β
观上的质点位移u ，从宏观上看，原子的位置可视为准连续的，原子的分离a l n )(+可视为准连续坐标x ，即
u
Ae Ae u t qx i t l n q i l n ===--++][])([ωω
于是（2）化成
22
22
2x u v t
u ??=?? 其中
m a
v β
=
2. 在一维双原子链中，如1>>m M ，求证
qa M sin 21β
ω=
)cos 21(222qa M m
m +=
βω
()()}]cos [(
12/122
22qa m M m M m M m
+++-++≈
}]cos 4)[(
12/122qa M m
m M m M m
++-+≈
β
}cos 42111{2qa M m
m
+
+≈β
}cos 1{22qa M m m +≈
β
qa M m m 22cos 12+=
∴βω)cos 21(22qa M m
m +≈β
220=-=M m A B ββ 故B ＝0， 重原子静止。
3.在一维无限长的简单晶格中，原子质量为M ，若只考虑近邻原子之间的相互作用，恢复力系数为β，试求格波的色散关系。
解：设原子的质量为 M ，第n 个原子对平衡位置的位移为un 第n+1和n-1个原子对平衡
??B 讨论当温度很高时，结果又会怎样？ 证明：按照量子理论，一个谐振子的能级是
ω
εη??
? ??+=n n 21
式中，ω为谐振子的角频率；n 取正整数。在热平衡条件下，谐振子
的平均能量为 ∑=n n
n P εε
式中
n
P 为谐振子处于能级
n ε的几率。若按玻耳兹曼统计计算，上式
写成
∑∑∞
=∞
=????????? ??+-??
???????
??+-??? ??+=00/2
1exp /21exp 21n B n B T k n T k n n ωωωεηηη
[]
∑∞
-/exp B
T k n n ωωηη在高温下，B ，有 ωωηηT k T k cth B B 22≈???? ?? 故得 T k B ≈ε
可见，在高温下，一个量子谐振子的平均能量与经典理论的结论相同。
5.在一维无限长的简单晶格中，若考虑原子间的长程作用力，第 n 个与第 n +m 或 n-m 个
原子间的恢复力系数为
m β，试求格波的色散关系。
解：设原子的质量为 M ，第n 个原子对平衡位置的位移为un 第n+m 和n-m 个原子对平衡位置的位移分别为un+m 与 un-m ，则第n+m 和n-m 个原子对第n 个原子的作用力为
)
2()()(,n m n m n m m n n m n m n m m n u u u u u u u f -+=---=-+-+βββ
第 n 个原子受力的总合为
∑∑∞
=-+∞=-+==1
1
,)
2(m n m n m n m m m n n u u u f F β
因此第 n 个原子的运动方程为
0z y x 解：2220z y x Cq Bq Aq ++=-ωωΘ
则 102
0202=-+-+-C
q B q A q z
y x ωωωωωω
这是q 空间的一个椭球面，其体积为abc π3
4
，而
2
/10A
a ω
ω-=
，2
/10B
b ω
ω-=
，2
/10C
c ω
ω-=
q 空间内的波矢密度()33
)2(2ππρV
L q =
??
? ??= ，故椭球内的总状态数N 为 ()2
/302
/13
1342ω
ωππ-??
?
???=ABC V N
所以
)
2
1sin()
21cos()21sin(21)21cos()21sin(212
2qa qa qa qa qa qa qa dq d m m m ωωω
ωω==
22
2221)]21(sin 1[21)21cos(21ωωωωω-=-==m m m a qa a qa a dq d 模式密度为
2
222
1
22
122)(ω
ωπωωπω-=
-=
m
m N
a L D
7. 已知一个频率为i ω的简谐振动在温度T 下的平均能量 1
21
/-+=
T k i
i i B i e ωωωεηηη 试用爱因斯坦模型求出由N 个原子组成的单原子晶体晶格振动的总能量，并求其
在高温和低温极限情况下的表达式。
解:由N 个原子组成的单原子晶体共有3N 个自由度，独立晶格振动方式数也等于
因斯坦模型下的零点振动能。 在低温极限下，x>>1，x x e e ≈-1，从(1)式得
T E B E B x B E e Nk k N
xe x T Nk E /32
3)21(3Θ--Θ+Θ=+=
8. 设晶格中每个振子的零点振动能为2ω
η，试用德拜模型求三维晶格的零点振
动能
解：状态密度
()()3
2
223v V V g ωπωωρ== 则
()ωωπωωωρεωωd v V d E D
D 32
20
002321η?
?==
D
D v V d v V ω
ωωπωωπ04320332163143ηη==? 4
3
2163D v V ωπη=
v
N V D 3
/126??? ??
=πωΘ D D N v V N v V E ωωππηη8961633
23
20=?=∴
9.设有三维间立方晶格，在德拜近似下计算比热，并论述在低温极限声子数目与T 3。
解：按照德拜模型, 晶体中的声子数目N’为
.
作变量代换
，
.
其中
是德拜温度. 高温时,
,
即高温时, 晶体中的声子数目与温度成正比. 低温时,
,
,
即低温时, 晶体中的声子数目与T 3成正比.
10. 有N 个相同原子组成的面积为S 的二维晶格，在德拜近似下计算比热，并论述在低温极限比热正比与2
T 。
证明：在k 到k dk +间的独立振动模式对应于平面中半径n 到n dn +间圆环的面积
2ndn π，且
()22
532222L s ndn kdk kdk d v ρ
ω
πρωωπππ===即则
()()2
3
3220
//2
22
22
333212121
m
D
D
B B x B B B B k T k T x D
D
d s k T s k T k T k T s d x dx
E E v e v e v e ωωωωρρρωωωωπππ????
? ?
????=
+==
---?
?
?
h h h h h h h h
20,(
)v s E
T E T C T T ?→∝∴=∝?3时，
11.有三维简单晶格，在德拜近似下计算比热，并论述在低温极限声子数目与T 3。 按照德拜模型, 晶体中的声子数目N '为
.
作变量代换
，
.
其中
是德拜温度. 高温时,
,
即高温时, 晶体中的声子数目与温度成正比. 低温时,
,
,
即低温时, 晶体中的声子数目与T 3成正比.
12.有N 个相同原子组成的体积为L 的一维晶格，在德拜近似下计算比热，并论述在低温极限比热正比与T 。.
13. 在一维无限长的简单晶格中，原子质量为M ，若只考虑近邻原子之间的相互作用，恢复力系数为β，试求格波的色散关系。
解：设原子的质量为 M ，第n 个原子对平衡位置的位移为un 第n+1和n-1个原子对平衡位置的位移分别为un+1与 un-1，则第n+m 和n-m 个原子对第n 个原子的作用力为
)
2()()(4111n n n n n n n u u u u u u u f -+=---=-+-+βββ
因此第 n 个原子的运动方程为
)
2(1122n n n n
u u u t d u d M -+=-+β
将格波的试解
)
(t qna i n Ae u ω-=
代入运动方程，得
计算色散关系为2
cq =ω的模式密度一维的模式密度。
解：一维情况下的q 空间中的等频面退化为两个等频的点，因此有 q 空间有两个等频点+q 和-q 。仿上面的方法可以得到：
1()2(2)(2)2()q dn L dq L g d Cq q ωωππω=
==?=??
15 对三维单原子点阵，计算德拜模型下的模式密度。 解：（ 解法一）
设横波和纵波具有相间波速v ，有
()()
()
()()
()
3
23
22323
23
2D
D
D
D K K K K
K K dK
g vK d dKK vK dKK vK ωδωπδωπδωπ<<<'=-=Ω-=
-?
??? (1)
令,z vK dz vdK ==，上式化为
()()2
23
3
2D D z dzz g z v
ωωδωπ<=
-?以上是就色散关系的一支求得的，考虑到一个波矢K 有三种偏振态，单原子点阵的色散关系有三支，纵波和横波有不同波速，总的模式密度应对各支求和，于是
()2
23
31122D l i g v v ωωπ??=
+ ???
式中
l
v 是纵波的波速，
i
v 是两个横波的波速．如果用v 表示纵波和横波的平均波速
333312l i v v v =+
德拜模式密度又可写为
()2
23
32D g v ωωπ=