晶格振动与晶体热力学性质热容

晶格振动与晶体热力学性质热容第1页，课件共47页，创作于2023年2月3、晶格振动模式密度晶格振动模式密度函数的定义表示在间隔内晶格振动模式的数目。constant确定了一个等频率面，那么在等频在q空间可计算如下：率面和之间的振动模式数目为第2页，课件共47页，创作于2023年2月首先计算N个波矢代表点在q空间的分布密度晶格振动模（格波）在q空间分布是均匀的：N很大，q值很密集，可认为是准连续的。由于q是限定在第一布里渊区的，而第一布里渊区在波矢（倒格子）空间的体积（倒格子原胞体积）为波矢q的数目等于N原胞（原子数）N个波矢代表点在q空间的分布密度为第3页，课件共47页，创作于2023年2月（频率为的等频率面间的体积）dqn表示沿ds面积元法线方向上的增量，因为沿法线方向频率的改变率第4页，课件共47页，创作于2023年2月得到模式密度的一般表达式知道了色散关系，便可由上式求得模式密度。对于具体的晶体，D(ω)的计算往往十分复杂，在一般讨论中，常采用简化的爱因斯坦模型及德拜模型。第5页，课件共47页，创作于2023年2月第六节晶格振动热容理论3.4.1热容理论本节主要内容:3.4.2爱因斯坦模型3.4.3德拜模型第6页，课件共47页，创作于2023年2月

引入声子概念后，研究晶格振动的热效应时，就可等效为研究3nN种声子组成的多粒子体系，在简谐近似下，这些声子是相互独立的，因而构成近独立子系。热容理论

热容量：一种物质每升高一摄氏度需要的能量或每降低一摄氏度释放的能量，被称为该物质的比热或热容量。定容比热定义：——固体的平均内能。

本节用统计理论的方法，讨论晶格振动的热容理论。第7页，课件共47页，创作于2023年2月固体的内能由两部分组成：

绝缘体：与温度有关的内能是晶格振动能量。

金属:与温度有关的内能由两部分组成，即晶格振动能量和价电子的运动能量。当温度不太低时，电子对比热的贡献远比晶格的贡献小，一般可以略去。本节只讨论晶格振动对比热的贡献。一部分内能与温度无关：例如，在简谐近似下，原子在平衡位置时的相互作用势能；

另一部分内能与温度有关。对比热有贡献的是依赖温度的内能。第8页，课件共47页，创作于2023年2月1、经典热容理论－－杜隆-帕替定律--------经典理论缺陷固体中的晶格振动的基本单元是谐振子。定容比热若固体中有N个原子，则有3N个简谐振动模，则总的平均能量

即定容比热是一个与温度和材料性质无关的常数。

根据经典统计理论的能量均分定理，在温度T时，每个自由度的平均能量是第9页，课件共47页，创作于2023年2月

高温和室温时这个理论结果与杜隆-帕替在1818年由实验发现的结果符合得很好；

低温时实验指出CV与温度T有关，即比热随温度降低的很快，当温度于绝对零度时，比热也趋于零。这个事实经典理论不能解释。为了解决经典理论的缺陷，爱因斯坦发展了普朗克的量子假说，第一次提出了量子的热容量理论。第10页，课件共47页，创作于2023年2月

2、量子热容理论简谐振动的能量本征值是量子化的，即频率为ωi的谐振子的振动能量为：其中代表零振动能，对比热没有贡献，略去不计，而将Ei写成：利用玻尔兹曼统计理论，在温度T时的单个谐振子的平均能量为：波尔兹曼权重x所有量子态求和第11页，课件共47页，创作于2023年2月第12页，课件共47页，创作于2023年2月

因此，在温度T时，频率为ωi的振动的平均能量为其中，表示温度为T时，振动模式为ωi的声子的平均数目。第13页，课件共47页，创作于2023年2月

把晶体看成一个热力学系统，晶体中有N个原子，每个原子有3个自由度；

在简谐近似下，各简正坐标Qi(i=1,2….3N)所代表的振动是相互独立的，因此因而可以认为这些振子构成近独立的子系；晶体有3N个正则频率，它们的统计平均能量应为：第14页，课件共47页，创作于2023年2月对实际晶体，晶格振动波矢q的代表点密集地均匀分布在布里渊区内，频率分布可以用一个积分函数表示，上式可改成积分形式计算。

设D(ω)dω表示角频率在ω和ω＋dω之间的格波数（即振动模式的数目）而且：模式密度D(ω)：单位频率区间的格波振动模式数目。又称角频率的分布函数。

ωm：最大的角频率，又称截止频率。

第15页，课件共47页，创作于2023年2月

只要知道晶格的模式密度D(ω)，就可以求出比热。平均能量可以写成：比热可写成：第16页，课件共47页，创作于2023年2月4、爱因斯坦模型

爱因斯坦模型的假设：固体中的原子都以同一频率ω振动，振动能量是量子化的。每一个原子都有三个振动自由度，每个振动自由度上有一个振子，固体中的N个原子可以看成3N个频率为ω的谐振子。根据以上假设，晶体的平均能量为：第17页，课件共47页，创作于2023年2月——爱因斯坦热容函数。第18页，课件共47页，创作于2023年2月则比热可简化为：令一个量子的能量等于一个经典振子的能量kBT，将用这种方法得到的T称为爱因斯坦温度，记为第19页，课件共47页，创作于2023年2月金刚石理论计算和实验结果比较

第20页，课件共47页，创作于2023年2月讨论：（1）在高温时:即在高温下，爱因斯坦近似过渡到经典的杜隆－帕替定律。第21页，课件共47页，创作于2023年2月（2）当温度T比爱因斯坦温度低很多时可以忽略比热公式分母中的1，则得到爱因斯坦模型的固体比热为：从上式可知，比热是随着温度指数下降的，这与很多固体在低温下的实验规律不符，而是更快地趋近于零。造成这一偏差的根源就在于爱因斯坦模型过于简单，它忽略了各格波对热容贡献的差异。第22页，课件共47页，创作于2023年2月按照爱因斯坦温度的定义可以估计出爱因斯坦频率，当相当于远红外光频率。频率为ω的一个格波的平均热振动能按照上式可以绘出格波的振动能与频率的关系曲线。

图中可以看出，格波的频率越高，热振动能越小。

爱因斯坦考虑的格波的频率很高，其热振动能很小，对比热的贡献本来就很小，当温度很低时，就更微不足道了。第23页，课件共47页，创作于2023年2月

爱因斯坦模型的单一频率格波实际上只适用于近似描写格波中的光学支，因为光学支一般频宽很窄，因而可以用一个固体频率描述。爱因斯坦模型实际忽略了频率较低的声学波对比热的贡献。而在低温时，声波对比热的贡献恰恰又是主要的。这就是为什么（16）式所示的比热随温度下降比实验结果更快的原因。本质上的原因：

当温度一定，频率越高的格波，其平均声子数越少。即频率高于的格波被“冻结”，对比热无贡献。第24页，课件共47页，创作于2023年2月5、德拜模型

德拜模型的基本假设：在三维晶体振动的能谱中忽略光学支对比热的贡献，将晶格视为连续介质，长声学波具有弹性波的性质。德拜引进了一个截止频率即德拜频率，以满足晶格振动的总自由度数（波的总数）限制条件即，不考虑频率超过德拜频率的高能量声子对固体比热的贡献。

三个声学支的色散关系简化为即一支纵波和二支偏掁方向不同的横波的波速相等。第25页，课件共47页，创作于2023年2月

下面我们先计算波矢q的频率分布函数：

在三维波矢空间中，N个波矢代表点在q空间的分布密度是

因此，在三维波矢空间中，波矢在ω到ω+dω两个等频面之间的振动模式数目为第26页，课件共47页，创作于2023年2月由于波的传播速度与波的传播方向q无关，在q空间等频面是球面，选用球坐标，所以第27页，课件共47页，创作于2023年2月根据上述模型，有：因此对于各向同性介质中的弹性波ω＝qvp，则可得频率在ω-ω+dω之间的模式数为：第28页，课件共47页，创作于2023年2月考虑到弹性波有三支格波，得出德拜近似的频率分布函数（模式密度）于是振动能量和比热分别为：第29页，课件共47页，创作于2023年2月截止频率可将模式密度（20）式代入振动模式的数目（7）式求出，即称截止频率为德拜频率，并记作对应是德拜温度。它是一个待定参数，由实验确定。令德拜温度定义为第30页，课件共47页，创作于2023年2月可得式中积分限第31页，课件共47页，创作于2023年2月在高温极限下:是小量。因此，比热的积分函数所以，高温比热即高温极限下，比热近似等于常数3NkB，与爱因斯坦模型的结果一致，也与杜隆－帕替定律相符。第32页，课件共47页，创作于2023年2月低温时:则（27）式积分上限可近似看作无穷大，将被积函数按二项式定理展成级数因此第33页，课件共47页，创作于2023年2月

德拜理论与很多实验事实符合。而且温度越低，近似越好。在低温下容易被激发的是长声学波振动，由于波长较长，晶体可看成连续介质，因而性质很象弹性波，这就是德拜近似取得成功的原因。这就是著名的德拜低温比热定律。所以第34页，课件共47页，创作于2023年2月例1、一维单原子布喇菲格子晶格振动的频率ω和波矢q的关系为其中m是原子质量，a是原子间距，β是原子间相互作用的力常数。1、按照ω和q和关系，求出晶格比热的表达式；2、给出高温、低温极限时比热随温度的变化关系。或者：3、按照德拜模型求出晶格比热的表达式；4、给出高温、低温极限时比热随温度的变化关系。第35页，课件共47页，创作于2023年2月先计算单位频率间隔的振动模式数（模式密度），即角频率的分布函数(1)：晶格振动的平均能为一维简单格子的色散关系dω区间对应两个同样大小的波矢区间dq。2π/a区间对应N=L/a个振动模式，单位波矢区间对应有L/2π个振动模式。dω范围包含的模式数为第36页，课件共47页，创作于2023年2月因此模式密度为由色散关系式得按色散关系可算出第37页，课件共47页，创作于2023年2月而频率为ω的谐振子的平均声子数目所以第38页，课件共47页，创作于2023年2月（2）高温极限为常数值。第39页，课件共47页，创作于2023年2月低温极限令即C正比于温度T。利用公式第40页，课件共47页，创作于2023年2月（3）：按德拜模型计算，弹性波处理由于得代入能量公式，所以第41页，课件共47页，创作于2023年2月（4）高温极限第42页，课件共47