晶格振动与晶体的热学性质习题集.docx

第三章晶格振动与晶体的热学性质1 .什么是简谐近似？解：当原子在平衡位置附近作微小振动时，原子间的相互作用可以视为与位移成正比的虎克力，由此得出原子在其平衡位置附近做简谐振动。这个近似即称为简谐近似。2 .试定性给出维单原子链中振动格波的相速度和群速度对波矢的关系曲线，并简要说明其意义。解：由一维单原子链的色散关系3=,可求得一维单原子链中振动格波的相速度为2而其群速度为qa2由(1)式和(2)式可做出一维单原子链中振动格波的相速度和群速度对波矢的关系曲线如下列图3.1所示：Sinqcl上图中A=aj2,=-。曲线1代表%,=色2-,曲线2代表Nma1qmqaT由(1)式及结合上图3.1中可以看出，由于原子的不连续性，相速度不再是常数。但当g0时,为一常数。这是因为当波长很长时，一个波长范围含有假设干个原子，相邻原子的位相差很小，原子的不连续效应很小，格波接近与连续媒质中的弹性波。=”口,这说明波 V in由(2)式及结合上图3.1中可以看出，格波的群速度也不等于相速度。但当40,7=N表达出弹性波的特征，当夕处于第一布区边界上，即4=2时，O=O,而PP矢位于第一布里渊区边界上的格波不能在晶体中传播，实际上它是种驻波。3 .周期性边界条件的物理含义是什么？引入这个条件后导致什么结果？如果晶体是无限大，g的取值将会怎样？解：由于实际晶体的大小总是有限的，总存在边界，而显然边界上原子所处的环境与体内原子的不同，从而造成边界处原子的振动状态应该和内部原子有所差异。考虑到边界对内部原子振动状态的影响,波恩和卡门引入了周期性边界条件。其具体含义是设想在一长为的有限晶体边界之外，仍然有无穷多个相同的晶体，并且各块晶体内相对应的原子的运动情况一样，即第7个原子和第fN+/个原子的运动情况一样，其中f=l,2,30引入这个条件后，导致描写晶格振动状态的波矢g只能取一些分立的不同值。如果晶体是无限大，波矢g的取值将趋于连续。4 .什么叫声子？对于一给定的晶体，它是否拥有一定种类和定数目的声子？解：声子就是晶格振动中的简谐振子的能量量子，它是一种玻色子，服从玻色一爱因斯坦统计，即具有能量为力WZ（夕）的声子平均数为对于一给定的晶体，它所对应的声子种类和数目不是固定不变的，而是在一定的条件下发生变化。5.试比拟格波的量子声子与黑体辐射的量子光子；“声子气体”与真实理想气体有何相同之处和不同之处？解：格波的量子声子与黑体辐射的量子光子都是能量量子，都具有一定的能量和动量，但是声子在与其它粒子相互作用时，总能量守恒，但总动量却不一定守恒；而光子与其它粒子相互作用时，总能量和总动量却都是守恒的。“声子气体”与真实理想气体的相同之处是粒子之间都无相互作用，而不同之处是“声子气体”的粒子数目不守恒，但真实理想气体的粒子数目却是守恒的。6 .晶格比热容的爱因斯坦模型和德拜模型采用了什么简化假设？各取得了什么成就？各有什么局限性？为什么德拜模型在极低温度下能给出精确结果？解：我们知道晶体比热容的一般公式为由上式可以看出，在用量子理论求晶体比热容时，问题的关键在于如何求角频率的分布函数夕（口）。但是对于具体的晶体来讲，夕（口）的计算非常复杂。为此，在爱因斯坦模型中，假设晶体中所有的原子都以相同的频率振动，而在德拜模型中，那么以连续介质的弹性波来代表格波以求出p（0）的表达式。爱因斯坦模型取得的最大成就在于给出了当温度趋近于零时，比热容Cy亦趋近于零的结果，这是经典理论所不能得到的结果。其局限性在于模型给出的是比热容CV以指数形式趋近于零，快于实验给出的以趋近于零的结果。德拜模型取得的最大成就在于它给出了在极低温度下，比热和温度尸成比例，与实验结果相吻合。其局限性在于模型给出的德拜温度。应视为恒定值，适用于全部温度区间，但实际上在不同温度下，德拜温度0。是不同的。在极低温度下，并不是所有的格波都能被激发，而只有长声学波被激发，比照热容产生影响。而对于长声学波，晶格可以视为连续介质，长声学波具有弹性波的性质，因而德拜的模型的假设根本符合事实，所以能得出精确结果。7 .声子碰撞时的准动量守恒为什么不同于普通粒子碰撞时的动量守恒？U过程物理图像是什么？它违背了普遍的动量守恒定律吗？解：声子碰撞时，其前后的总动量不一定守恒，而是满足以下的关系式其中上式中的G“表示一倒格子矢量。对于G=0的情况，即有内q+人q2=力q3,在碰撞过程中声子的动量没有发生变化，这种情况称为正规过程，或N过程，N过程只是改变了动量的分布，而不影响热流的方向，它对热阻是没有奉献的。对于G“WO的情况，称为翻转过程或U过程，其物理图像可由下列图3.2来描述：图3.2U过程物理示意图在上图3.2中，q+q2是向“右”的，碰撞后q.,是向“左”的，从而破坏了热流的方向，所以U过程对热阻是有奉献的。U过程没有违背普遍的动量守恒定律，因为声子不是实物量子，所以其满足的是准动量守恒关系。8 .简要说明简谐近似下晶体不会发生热膨胀的物理原因；势能的非简谐项起了哪些作用？解：由于在简谐近似下，原子间相互作用能在平衡位置附近是对称的，随着温度升高，原子的总能量增高，但原子间的距离的平均值不会增大，因此，简谐近似不能解释热膨胀现象。势能的非简谐项在晶体的热传导和热膨胀中起了至关重要的作用。N个相同原子组成的一维单原子晶格格波的态密度可表示为P(G)=(tfl-)2。式中例是格波的最高频率。求证它的振动模总数恰好等于No解：由题意可知该晶格的振动模总数为(O=cq2和G=g-cq2,试导出它们的状态密度表达式。解：根据状态密度的定义式可知p)=Iim(1)>oAG其中加表示在GG+Ag间隔内晶格振动模式的数目。如果在q空间中，根据。(q)=ct"作出等频率面，那么在等频率面/和/+0之间的振动模式的数目就是4?。由于晶格振动模在q空间分布是均匀的，密度为V(2)3(V为晶体体积)，因此有Ie十A®=?dSdq(2)(2万)T将(2)式代入(1)式可得到状态密度的一般表达式为夕(G)=y7(3)(2;T)-JI/I(3)式中,g3)|表示沿法线方向频率的改变率。当G=Cq2时，将之代入(3)式可得当。=OO-C将之代入(3)式可得加,原子间距为/2,力常数交错为四，2的一维原子链振动的色散关系。当夕2=lA时，求在4=()和4=工处的火夕)，并粗略画出色散关系。在最近邻近似和简谐近似下，第2n和第(2n÷l)个原子的运动方程为wzl=A+-)一A区一x2z,-),月(1)fn¥=&2+2-+l)-鱼2%)at当月2=1OA时，上述方程组(1)可变为d2m-JTl=10AU2n+-)-(积-h-i)中(2)md*=B*2/2一兀2+1)-10用(!-)为求格波解，令|(2怨-砌"Ae(3)H(2n+D岑-M2,1÷=Be将(3)式代入(2)式，可导出线性方程组为-2)A-(Ov2+e-iqa/2)B=O,2m(4)-*(eiqa,2+10e'z2)A+(叫-2)B=0mm令2L=说，从A,8有非零解的系数行列式等于零的条件可得m(Il口；-2)2-说(IOe碗a+e-驷/2e碗/2+Ne-驷/2)=0由(5)式可解出d?=G：(U±j20cos4a+101)当夕=0时，cosqa=1,+=V220，M=O当g=工时，cosq=-l,+=V269o>co_=y2coa其色散关系曲线如下列图3.4所示：图3.4原子间的力常数不相噂的双原子链的晶格振动色散关系曲线12.如有一维布喇菲格子，第2个原子与第2+1个原子之间的力常数为6；而第2个原子与第2-1个原子的力常数为夕(1)写出这个格子振动的动力学方程；(2) 说明这种情况也有声学波和光学波；(3) 求q=0时，声学波和光学波的频率；TT(4) 求g=±(。为晶格常数)时，声学波和光学波的频率。解：(1)此题与(11)题根本相似，在最近邻近似和简谐近似下，同样可以写出第2和第2+1个原子的动力学方程为wzl=伙巧”+1-)一,(2l.-x2l,-)理(1)w-¾-=-+l)-U2n+l-x2n)(2)为求出方程组(1)的格波解，可令,Adg)Wex2n+_Be"(2"+i)缈-于是将(2)式代入(1)式，可导出线性方程组为/I三一2)A-g*a+2CB=O<2mmC)一(26碗+2«一碗)A+(i-Q2)B=ommm令红£=诙，殳=32=从A、B有非零解的系数行列式等于零的条件可得mmm(就-2)2-(6914+21cos2qa)=0(4)由(4)式可解出2=就±Jty14+0；+2说妖CQslqa(5)由此可知，的取值也有+和口_之分，即存在声学波和光学波(3)由(5)式可知当夕=0时，CoS2网=1,有声学波频率砧=y/;一(说+就),光学波频率g=y/;+3；+&；)(4)同样由(5)式可知当q=±时，cos2qa=-1,有2a声学波频率M=M-忸；-蜀,光学波频率+=M+Ma13.在一维双原子链中，如Mm>>l,(1)求证：(2)画出/与4的关系图(设Ms=IO).解：(1)在一维双原子链中，其第2个原子与第2+1个原子的运动方程为m=以工2T+2n+l-2/“)d(1)M="均+2-+)为解方程组（1）可令匕_ei(2n)qa-(jtx2n+1=&水2+IMe将（2）式代入（1）式可得出（-2）A-（"CoSqa）B=0<mM（3）-（j-cosa）A+（*-2）B=0从A、B有非零解，方程组13）的系数行列式等于零的条件出发，可得可解出得4(2+4±J(2+分-注&山2qa(4)MmNMmMm当(4)式中取“一”号时，有2伙M+Mz4Mtn.2一69.=-1-(17sn")2(5)mM(M+m)29.,Mm>>,：.(5)式中有理但以丝=幺Sin“。坐si112go=驯sin%<<lMmMinm(M+m)2M-M那么（5）式可简化为当（4）式中取“+”号时，有4Mm 21 ÷-cos qa(M-m)2(6)2(M+m)0(M-tn)5=1MtnMmVM11>>1,(6)式中有P(M+m)pM_p_(M-喻pM_2MmMmmMmMmm那么(6)式可简化为1.栏(YcoL（2）当Mm=10时，那么（4）式可化为此时，0与g的关系图，即色散关系图如下列图3.5所示：图3.5一维双域子链振动的色散关系曲线14 .在一维复式格子中，如果n=5l.67xK)-24g,mw=4,=.5Nm.求：(1)光学波频率的最大值、最小值及声学波频率的最大值；(2)相应的声子能量是多少eV?(3)这3种声子在300K时各有多少个？(4)如果用电磁波激发光频振动，要激发最大光学频率的声子所用的电磁波长在什么波段？解：(1)由于光学波频率的最大值和最小值的计算公式分别为：.Mm5×1.67×1024,.-24“工八工日上式中=6.68×1()-4g为约化质量m+M机/M+11/4+1所以有：而声学波频率的最大值的计算公式为：所以有：(2)相应的声子能量为：(3)由于声子属于玻色子，服从玻色一爱因斯坦统计，那么有(4)如用电磁波来激发光频振动，那么要激发最大光学频率的声子所用的电磁波长应满足如下关系式：15 .在一维双原子晶格振动的情况下，证明在布里渊区边界夕=±处，声学支格波中所有轻原子加静2a止，而光学支格波中所有重原子M静止。画出这时原子振动的图像。解：设第2个原子为轻原子，其质量为优，第2+1个原子为重原子，其质量为，那么它们的运动方程为d2x二以T+n+l-2)Mdjr=Kn+一+l)为解方程组(1)可令0_AeH(2")g-d'x2m+i=/g+m将(2)式代入(1)式可得出(-2)A-(-cosa)B=0<mm(3)-(*cosq)A+(工-2)B=0从A、B有非零解，方程组(3)的系数行列式等于零的条件出发，可得可解出得=(-+-)±J(g+2)2-4与”山2W(4)MmNMmMnI令夕=±爰，那么可求得声学支格波频率为力光学支格波频率为。+由方程组(3)可知，在声学支中，轻原子机与重原子M的振幅之比为由此可知，声学支格波中所有轻原子也静止。而在光学支中，重原子M与轻原子加的振幅之比为由此可知，光学支格波中所有重原子M静止。此时原子振动的图像如下列图3.6所示：图3.6(a)声学支格波原子振动图；(b)光学支格波原子振动图16.从一维双原子晶格色散关系出发，当加逐渐接近M和m=M时，在第一布里渊区中，晶格振动的色散关系如何变化？试与一维单原子链的色散关系比拟，并对结果进行讨论。解：一维双原子品格的色散关系为由上图可以看出，当阳逐渐接近”时，霍第一布里渊区边界，即g=±二处，声学波的频率开始2a增大，而光学波的频率那么开始减小，而当机二M时，那么声学波的频率和光学波的频率在g=±处相等，都等于J°而在一维单原子链中，其色散关系为。2="si112史，由此可见，在一维单原子链中只存在一支M2格波，其色散关系曲线与一维双原子链中的声学波的色散关系曲线根本相似，在其布里渊区边界，即q=±三处,其格波频率为。=2、但，是双原子链的格波在布里渊边界的频率值的2倍。a4VMoq一_N个原子细成，试用德拜底K证明格波的状态密度为2。ap()=-.式中例M为格波的截止频率。解：在德拜模型中，假设晶体的振动格波是连续介质的弹性波，即有色散关系那么格波的状态密度为V22又根据Jp()d=3N0式代入(3)式得FVI"2.CAraco3Nvp(4)由(4)式可得vP=Vn9N把(5)式代入(2)式即可得p()=-2-,试用德拜模型求一维、二维和三维晶体的总零点振动能。设原子总数为N,一维晶格长度为2二维晶格的面积为S,三维晶格的体积为V。解：(1)一维晶体的总零点振动能为：设表示角频率在33+d之间的格波数，而且p()d=N上式中：GW,是最大的角频率；N为晶体中的原子数。那么上述的总零点能可以写成:=-p()d考虑到维晶体中，其状态密度为:/、dZdZdqp(co)=-ddqd(3)由于德拜模型考虑的是长声学波的影响，而长声学波可以看成连续媒质弹性波。对于弹性波，一个波矢对应一个状态，那么有：JZdq(4)对于弹性波，=vpq,那么d-J-=vPdq(5)将(4)和(5)式代入(2)式，得:P(三)=(6)将(6)式代入(1)式，可得：G,”=将(6)式代入(2)式，可得维晶体的总零点振动能：(2)对于二维晶体来说，计算其总零点振动能根本方法与一维晶体的方法相似，只是对于(1)式要改p()d=2N(7)而对于二维晶体，其状态密度函数为：P(三)=-vp将(8)式代入(7)式可得：用J将(8)式代入(2)式可得二维晶体的总零点振动能为：(3)对于三维晶体来说，计算其总零点振动能根本方法与一维晶体的方法也根本相似，只是对于(1)式要改为：C(t)mJ。p()d=3N(9)而对于三维晶体，其状态密度函数为：(10)3V¾?22Vp将(10)式代入(9)式可得：3帆=LJ将(10)式代入(2)式可得三维晶体的总零点振动能为：19.应用德拜模型，计算一维、二维情况下晶格振动的状态密度、德拜温度、品格比热容。解：在德拜模型中，假设晶体的振动格波是连续介质的弹性波，即有色散关系co=vpq(1)(1)在一维情况下，品格振动的状态密度为p(,co)=-7-2=(2)2dc加Pdq上式中，L表示一维晶格的总长度。%又由关系式p()d=N(3)将式(3)代入式(2)可得1d=N,由此求得例”=JowPL-疳INUn于是德拜温度。力=,+f , (其中R = kJ晶体的比热容为=7or就VP0(2)在二维情况下，晶体振动的格波有2支，即支纵波和一支横波，在德拜模型中，假设纵波和横波的波速相等，都等于勺，即纵波和横波都有如下的色散关系先考率纵波，其状态密度为类似地可以写出横波的状态密度为p,(o)=2r2叫,加起来总的状态密度为/、/、/、SCO(4)0(G)=1(G)+G(G)=-啊又由关系式Jp()d = 2N (5)O将(4)式代入(5)式得'沼d=2N,由此可得©,“=%殳大叫IS于是得德拜温度为=也L=_lj网/TS)而晶体的比热容为A372q介©E的丁3A(其中R = kB空？f上JdX成R；I(/一1)XlO12NZm23O试计算金刚石的德拜温度。解：假设金刚石的原子振动的格波为一连续介质的弹性波，其波速为匕，=但二、回咚=IWlOFsPpV3.5×IO3而又金刚石的原子密度为n=旦%=3,5X1×6.02XlO23=1.756X1029个11Mc12×103由此可知金刚石的德拜温度为=2817K21.具有简单立方布喇菲格子的晶体，原子间距为2X10/°m,由于非线性相互作用，一个沿100方向传播，波矢大小为q=l3xlOK)m的声子同另一个波矢大小相等当沿110方向传播的声子相互作用，合成为第3个声子，试求合成后的声子波矢。解：易知简单立方格子的倒格子仍是一简单立方格子，其倒格基矢E、b2和b3互相垂直，长度为=O*/=3.14x10°mL第一布里渊区就是原点和六个近邻格点连线的垂直平分面围成的立方a2×10,°体。又因为由此可知q+q2落在第一布里渊区之外，即可知题所述两声子的碰撞过程是一个翻转过程或U过程，此时两声子的碰撞产生第三声子满足准动量守恒，即有力qI+方q2=方q3+方G(其中G“表示一倒格矢)为使q-3落在第一布里渊区里，取G“=3.14XK)K)i,那么有其大小为q.J=卜0.92×1010i+0.92×IO10j=1.3×1010me2B+0r*式中8为待定常数；厂为近邻原子间距。求该晶体的线膨胀系数。近邻原子的平均距离为3X10"°m°解：由平衡条件也Ddr=O,可得e22乌=0由此可得8 =4 在 °FW*%于是可求得C（八2!dr0 =1 /% 、 g一（H那么线膨胀系数为=5.4×10-5K-1