一维双原子链的色散关系毕业论文

一维双原子链的色散关系 The dispersion relation of the one-dimensional diatomic chain 学生姓名: 吴刚 所在专业: 应用物理学 班 级: 物理091 指导教师: 陈冰泉，刘彦平 申请学位: 理学学士 论文提交日期: 2013-6-5 论文答辩日期: 2013-6-18 学位授予单位: 青岛理工大学 青岛理工大学毕业设计论文 摘 要 物理学中对晶格振动的研究一直是一个重要且有意义的课题。关于晶格振动的研究 通常建立在原子链的研究上。本文在介绍关于原子链研究基础理论的基础上,讨论了一 维单原子链晶格的色散关系,从一维双原子链的角度介绍了晶格的色散关系,然后在前 面讨论的基础上对三维晶格的色散关系进行了推导。关于原子链色散关系的研究,让我 们对于晶格振动有个更加清楚地认识。论文重点介绍了一维双原子链的色散关系,在公 式推导的基础上,作者完成计算机编程和模拟计算,得到色散关系的曲线。 关键词:晶格振动,一维单原子链,一维双原子链,三维晶格,色散关系 Abstract The study of lattice vibration has been an important and meaningful topic in physics. The investigation of lattice vibration is usually based on the study of atomic chain. With the introduction of the theoretical basis of atomic chain, this thesis discusses the dispersion relation of one-dimensional monatomic chain lattice, as well as the dispersion relation of one dimensional diatomic chain lattice. Based on the knowledge above, the equations for describing the dispersion relation of three dimensional lattice are then derived. The study of dispersion relations allows us to have a more clear understanding of lattice vibration. This thesis mainly presents the study and discussion of the dispersion relation of one dimensional diatomic chain. In addition to the equation derivation, we carry out programming and simulations for obtaining some important dispersion-relation curves. II 青岛理工大学毕业设计论文 目录 前 言 ............................................................................................................................ 1 第一章 理论基础 ............................................................................................................ 3 一维单原子链的色散关系 ..................................................................................... 6 第二章 2.1 建立振动模型 .................................................................................................... 6 2.2 建立振动方程并求解 .......................................................................................... 6 2.3 玻恩-卡曼条件 .................................................................................................. 8 w,q2.4.的函数关系 ............................................................................................ 10 第三章 一维双原子链的色散关系 ................................................................................... 13 3.1建立振动模型................................................................................................... 13 3.2 原子运动方程的求解 ........................................................................................ 13 3.3 周期性边界条件 .............................................................................................. 15 3.4 对于声学波和光学波的讨论 .............................................................................. 16 第四章 三维晶格振动的推导.......................................................................................... 21 4.1 一维多原子链问题的处理 ................................................................................. 21 4.2 建立三维模型和求解运动方程 ........................................................................... 21 ,q4.3 波矢的取值和范围 ........................................................................................ 23 4.4 理论上的计算 .................................................................................................. 25 第五章 结论和讨论 ....................................................................................................... 28 致 谢 .......................................................................................................................... 29 参考文献 ..................................................................................................................... 31 III 青岛理工大学毕业设计论文 前 言 讨论晶体结构时,我们把晶体内的原子看作是处于自己的平衡位置上固定不动的。但实际上,物质是在不断运动的,而量子力学告诉我们,即使达到绝对零度,仍然具有零点能的振动。我们一般常用晶体中的格点来表示原子的平衡位置,用晶格振动来表示原子在各个点附近的振动。晶格振动是研究固体宏观性质和微观性质的主要基础,人们最早对于晶格振动的一些研究是从晶体的热力学性质开始的。上世纪人们根据经典的统计规律来对杜隆-珀替经验定律进行说明，每摩尔固体都会有3N个振动自由度,按照能量的均分定律,每个自由度的平均热能是kT,那么摩尔热容量就为3Nk=3R，,这是把原子振动和热容量相互联系起来的一个非常重要的成就。但是,上世纪人们通过大量的实验研究居然发现杜隆-珀替定律仅仅只是在室温或者更高的温度下,才与固体基本上是一致的,可是在温度较低的情况下,伴随着温度的降低,固体的热容量也开始再不断降低。在这个问题的困扰下,爱因斯坦进一步对于普朗克的量子假说进行了发展,第一次得到量子的热容量理论,发现了热容量在低温的情况下下降,并在T趋于零时趋于零的结论,这项在量子理论的发展中占有重要地位的研究,对于原子振动也有重要的影响。量子理论的热容量值和经典不同,它与原子振动的具体频率有关,也推动了对固体原子振动进行了具体的研究。 研究晶格振动的意义远远不限于热学性质。晶格振动也是研究晶体的电学性质、光学性质、超电导性、结构相变、磁性………等一系列物理问题,晶格振动都起着很重要的作用。 1 青岛理工大学毕业设计论文 晶格振动虽然是一个十分复杂的多粒子问题,但是在一定条件下,依然可以在经典范畴求解,一维原子链的振动就是最典型的例子,它的振动简单可解,又能反映出晶格振动的基本特点。原子振动也不是孤立的,而是相互联系着的,在晶体内形成各种模式 [1,2]波,我们通过研究原子链振动得出的色散关系能够帮助我们了解晶格振动. 2 青岛理工大学毕业设计论文 第一章 理论基础 , Rn如果晶体中包含N个原子,我们把平衡位置设为,把偏离平衡处的位移矢量设为 ,,,,'R(t),R(t)，,(t)(t),nnnn,则得原子位置为。在处理一些小振动的问题时,我们常常选 , ,n用偏离平衡位置的量为宗量。用分量来表示位移矢量,所以N个原子的位移矢量总共 ,i会有3N个分量,写成，i=1,2,3……,3N，。N个原子体系所包含的势能也可以在平衡位置附近用泰勒级数展开: N3,V,，(),，高阶项VV,i00,,i,1i ，1.1， 下脚符号0说明是在平衡位置时才所具有的值。 那么设V0=0,并且 V,(),00,,i ，1.2， 忽略高于二阶的高阶项,我们得到: 23N,1V,(),,V,0ij,,2,,,i,j1ij ，1.3， ,i该系统的势能函数就保留至的二次方程,我们称它为简谐近似。处理小的振动问题时一般都采用简谐近似,至于一个具体的物理问题对于它是否可以采取简谐近似,还 3 青岛理工大学毕业设计论文 要看在简谐近似条件下我们推导得到的理论结果是否与实验结果互相一致。在有些物理的问题上就需要考虑到高阶项的作用,一般称为非谐作用。 为了使问题简化,引入所谓的简正坐标: Q,Q,？？，Q123N ,i简正坐标与原子的位移坐标之间通过如下形式的正交变换相互联系: 3N m,,aQ,iiijj,j1 ，1.4， 应用正则方程可以得到: 2，，Q，wQ,0ii i=1,2,……,3N ，1.5， 这是三个相互无关的方程,表明各个简正坐标描述独立的简谐振动,其中任意简正坐标的解为: Q,Asin(wt，,)i ，1.6， ww,2,viii是振动的圆频率,原子的位移坐标和简正坐标间存在着正交变换关系 Qj，1.4，。只考察某一个的振动时,公式，1.4，可以化为: ,ij,,Asin(wt，,)ijmi ，1.7) 这表明,一般讲,一个简正坐标并不是表示某一个原子的振动,而是表示整个晶体所有原子都参与的振动,而且它们的振动频率相同。我们常称体系中所有原子一起参与 [3]的共同振动称为振动模。 4 青岛理工大学毕业设计论文 ,,i,p,Qii按照一般方法把写成,得到谐振子方程的解: 1,,,,n，,w,,iii2,, (1.8， 2,,,w,,，，，，,Q,exp,H,nini,,i,2,, ，1.9， 系统的本证态为: N31,,E,n，,w,,,ii2,,,i1 (1.10) N3 ，，，，,Q,Q,？,Q,,Q,Nni123i,i1 (1.11) 5 青岛理工大学毕业设计论文 第二章 一维单原子链的色散关系 2.1 建立振动模型 假设一维无限长的单原子链,在平衡时相邻原子的距离为a，即原胞体积为a，,每个原胞内含一个原子,质量为m,原子限制在沿链的方向运动,偏离格点的位移用……,-2,Un-1,Un,……表示,如图1.1表示。 Un 图2.1 一维单原子链 2.2 建立振动方程并求解 假如只考虑近邻原子间存在互相作用,互相作用能可写成: 12()()va，,,va，,,，高阶项2 ，2.1， 2a,,其中表示对平衡距离的偏离。按照一般小振动近似互相作用能保留到项,即简谐近似,在这种近似下,相邻原子间的作用力为: 6 青岛理工大学毕业设计论文 dvF,,,,,,d, ，2.2， 它表明存在于相邻原子间的是正比于相对位移的弹性恢复力。 求解链的振动模,考察图2.1中第n个原子的运动方程,它受到左右两个近邻原子对它的作用力: ，，mu,,(u,u),,(u,u),,(u，u,2u)nn，1nnn,1n，1n,1n ，2.3， 每个原子都对应一个方程,若原子链有N个原子,则有N个方程,公式，2.3，实际上表示着N个联立的线性齐次方程。 验证方程具有下列“格波”形式的解: i(wt,naq)u,Aen (2.4) 其中w、A为常数,代入方程，2.3，,有 2i(wt,naq)i(wt,(n，1)aq)i(wt,(n,1)nq)i(wt,naq)m(iw)Ae,[Ae，Ae,2Ae], 241,,22w,[1,cosaq],sin(aq)mm2 (2.5) 公式(2.5)与n无关,表明N个联立的方程都归结为同一个方程。即是说,只要w和q之间满足公式，2.5，的关系,公式，2.4，就表示了联立方程的解。我们通常把w和q的关系称为色散关系。 式，2.4，与一般连续介质波 7 青岛理工大学毕业设计论文 x,i(wt,2)i(wt,qx),Ae,Ae (2.6) ,2q,,,有完全相似的形式。其中w是波的园频率,是波长,称为波数。区别是: na在连续介质波中x表示空间任意一点,而在，2.4，中只能取格点的位置,这是一系 aq列呈周期排列的点。相邻原子之间的位相差为,格波与连续介质波一个重要的区别 (2,)aq的涵义。如果把式，2.4，中的改变一个的整数倍,那么所有原子的在于波数q aq振动实际上没有任何不同。所以可限制在下面范围内: ,,，aq,, ，2.7， 或: ,,q,，,aa ，2.8， 这个范围以外的q值,并不能提供其他不同的波,这个q的取值范围常称为布里渊区。 2.3 玻恩-卡曼条件 前面所考虑的运动方程实际上只适用于无穷长的链,然而一个有限的链两端的原子显然和内部的原子有所不同。例如,在只有近邻原子作用时,最两端的两个原子只受到一个近邻原子的作用,所以它们具有与其他原子形式不同的运动方程。为了避免这种情况,玻恩-卡曼提出如图2.2所示的包含N个原胞的环状链作为一个有限链的模型,它包含有限数目的原子,而且保持所有原胞等价。以前的运动方程依旧适用如果N很大使得环的半径很大,沿环的运动依旧可以看成是直线运动。和以前的差距仅仅在于必须考 8 青岛理工大学毕业设计论文 虑链的循环性。即,原胞的标数n增加N,振动情况必须复原。参见格波解公式，2.4，的形式可以知道,这等于要求: 图2.2 一维链的玻恩-卡曼边界条件 ,i(Naq)e,1 (2.9) 或者 2,q,，hNa ，h是整数， (2.10， ,,NN,，,，aa22前面指出,q的取值范围是由到,所以式，2.9，中的h只能取到,一共有N个不同值。所以,N个原胞组成的链,q可以取N个不同的值,每一个q值对应着一个格波,共有N个不同的格波。N就是一维单原子链的自由度数,表明,现已取得了链的全部振动模。 玻恩-卡曼的模型相当于要求一个有限链头尾相衔接,起到了一个边界条件的作用,其实我们也注意到,用这个模型并没有改变运动方程的解,只是对解提出了一定的条件，2.9，.我们将它称为玻恩-卡曼条件,或者称为周期性边界条件. 9 青岛理工大学毕业设计论文 w,q2.4.的函数关系 ,w(q) 我们这里对式，2.5，做两点讨论。表面上看一个q对应两个频率,事实上, 的偶函数,所以只需要取式，2.5，的正根就足够了,因此式，2.5，可以写成: 由于w是q ,1w,2sinaqm2 ，2.11， ,,,，w,qaa图2.4画出了关系之间的函数曲线。由于格波的特性,q取值在到之间, 的允许值为这一区间中的均匀分布的N个点。 由于周期性边界条件q 我们令m=8.35×10-24,β=1.5N/m,并且为了更易显示图形我们令ɑ=1,利用matlab 软件把公式，2.11，编程,如下: 图2.3 一维单原子链色散关系编程 10 青岛理工大学毕业设计论文 w,q 图2.4 一维单原子链的函数关系 根据分析力学的原理,引入简正坐标求解链的振动模。我们可得每一个简正坐标,对应一个谐振子方程,波函数是以简正坐标为宗量的谐振子波函数,能量本征值为: 1,,,,n，,w,,nqq2,, ，2.12， 所以,由N个原子组成的一维单原子链,其振动模为N个格波,在简谐近似下格波是互相独立的,格波的振幅对应着系统的简正坐标,按量子理论每种简正振动的能级是量 ,wq子化的,能量激发的单元时。 ,wq当以上的结论用声子的“语言”描述,声子就是格波的量子,它的能量等于。 1,,n，,w,,qq2,,一个格波,也就是一个振动模,称为一种声子;当这种振动模处于本征态的 11 青岛理工大学毕业设计论文 nnqq时候,一般称有个声子,个声子数;当电子或者光子与晶格振动互相作用时,交 ,w,wqq换的能量为,假如电子从晶格中获得的能量,称为吸收了一个声子,如果点子 ,wq给了晶格的能量,称为发射一个声子。我们用声子的“语言”描述晶格的振动不但会使的表述简化,而且还有这深刻的理论意义。其实声子不是真实的粒子,我们成为“准粒子”,它反映的是晶格原子集体运动状态的激发单元。 12 青岛理工大学毕业设计论文 第三章 一维双原子链的色散关系 3.1建立振动模型 一维双原子链可以看做是最简单的复式晶格:每个原胞内包含2个不同的原子P和Q,如图3.1所示。在平衡时相邻原子距离用a表示,P、Q质量用m和M表示。原子限制在沿链 xxx2n,12n，12n的方向运动运动,偏离格点的位移用……,,,,……表示。 图3.1 一维双原子链的模型 3.2 原子运动方程的求解 仍假设只有近邻原子间存在互相作用能,互相作用能取简谐近似、类比一维单原子链的情况,可以得到原子的运动方程: ，，mx,,(2x,x,x),2n2n2n，12n,1 P原子: ，，Mx,,(2x,x,x),2n，12n，12n，22n Q原子: ，3.1， 这是两个典型的运动方程。当原子链包含N个原胞，即有N个P原子和N个Q原子，时,它事实上代表的是2N个方程的联立方程组。这个方程组有如下形式的格波解: 13 青岛理工大学毕业设计论文 i[wt,(2na)q]x,Ae2n i[wt,(2n，1)aq]x,Be2n，1 ，3.2， 把式，3.2，带入式，3.1，,除去共同的指数因子后,可以得到: 2,iaqiaq,mwA,(e，e)B,2A,, 2,iaqiaq,MwB,(e，e)A,2B,, ，3.3， 方程与n无关,表明所有联立方程对于格波形式的解式，3.2，都归结为同一对方称。 式，3.3，可看作是以A、B为未知数的线性齐次方程: 2(mw,2)A，2cosaqB,0,, 22cosaqA，(Mw,2)B,0,, ，3.4， 它的有解条件是: 2,,mw,22cosaqB 22,cosaqMw,2, 4222mMw,2(m，M)w，4sinaq,0,, = ，3.5， 22ww则可以看作是决定的方程,从而得到两个值 1,,2，mM4mM,2w2，2,w,,1,[1,sinaq]2,,2W,，mM(mM),,, ，3.6， 22ww，,把和代回公式，3.4，,可以求的相应的A和B的解: 2,Bmw2,，(),,，A2,cosaq 14 青岛理工大学毕业设计论文 2,Bmw2,,(),,,A2,cosaq ，3.7， 3.3 周期性边界条件 由格波解，3.2，可知道相邻原胞之间，原胞为2a，的位相差2aq,所谓相邻原胞之 2,间位相差应理解为相邻原胞P原子，或者Q原子，之间的位相差。所以,若把2aq改变的整数倍,所有原子的振动实际上没有任何不同,这说明q的取值范围限制在: ,,,2aq,，, ，3.8， 或者 ,,,,q,2a2a ，3.9， 范围内,这个范围就是一维双原子链的布里渊区。在这个范围内任意的波数q有两 ww，,个格波解,它们的频率为式，3.6，所给出和,和一般波的解一样,格波解可以有任意的振幅和位相,但是两种原子的振幅比和位相差事确定的,并由A、B的解决定,仍采用周期性边界条件: N(2aq),2,h , ，h=整数， 即: hq,2,2Na ，3.10， 15 青岛理工大学毕业设计论文 ,,NN,,，2a2a22 由于q的取值范围是由到,所以式，3.10，中的h取值只能是到,一共有N个不同取值。所以,由N个原胞组成的一维双原子链,q可以取N个不同的值,每个q对应两个解,总共加起来,共有2N个不同的格波,数目恰好等于链的自由度,这说明,我们已经得到链的全部振动模。 3.4 对于声学波和光学波的讨论 ww，,我们根据公式，3.6，运用画图软件我们能得到各个不同q的格波频率和,属ww，,于的格波我们成为光学波,属于的格波我们称为声学波。 对于一维双原子链我们令m=5×1.67×10-24g,M/m=4,β=1.5×10N/m,为了更清楚地显示图形,我们让a=1,编程过程如下图3.2中所示: [4] 图3.2 一维双原子链色散关系编程 16 青岛理工大学毕业设计论文 w,q在下图3.3中声学波和光学波,我们也称为的色散关系。 图3.3 声学波和光学波 q,0的长波在许多实际问题中具有特别重要的作用,光学波和声学波的命名也主 要是根据它们在长波极限的性质,下面重点讨论长波极限的问题: w, 错误!未找到引用源。?先讨论声学波的长波极限的情况,从公式，3.6，可 q,0w,0,以看出,当时,,正如图3.3所示。当q很小时: 4mM4mM22sinaq,(aq),,122(m，M)(m，M) ，3.11， 2q可以把式，3.6，中根式对展开得到: 17 青岛理工大学毕业设计论文 2,22w,(aq),m，M 或者: 2,w,aq,m，M ，3.12， 表明对于长声学波频率正比于波数正比于波数。长声学波就是把一维链看做连续介 q,0w,质时的弹性波,这也就是为什么称支为声学波的原因。对于长声学波,当时 w,0,,因此: B(),1A ，3.13， 表明在长声学波时,原胞中的两种原子的运动基本是完全一致的,振幅和位相都没 有什么差别。 q,0 错误!未找到引用源。?对于长光学波,当时,频率趋于下列有限值: 2,w,，mM m，M ，3.14， cosaq,1代入式，3.7，,并且令,我们能得到: Bm(),,，AM ，3.15， q,0当时,由式，3.2，可知,同一种原子具有相同的相位,所以每一种原子，P q,0原子或者Q原子，形成的格子象一个刚体一样整体的振动,公式，3.15，表明,在 18 青岛理工大学毕业设计论文 时,两种原子振动有完全相反的位相,长光学波的极限在实际上是P原子和Q两个格子的相对振动,振幅中保持他们的质心不变。 在离子晶体中的长光学波有特别重要的作用,因为,不同离子间的相对振动产生一定的电偶极矩,从而可以和电磁波互相作用。具体分析对于单原子过程的一级谱,电磁波只和电磁波互相作用。具体分析证明对于单原子过程的一级谱,电磁波只和电磁波相 w~q互作用,如果它们具有相同的频率就可以发生共振,在图3.4中,我们把光波的关系: w,cq(c,光速)000 和格波画在同一图中,代表光波的直线与光学波的图线的交点对应于它们共振的情 c0况。代表长声学波的直线的斜率c为弹性波速度,仅仅是光波直线斜率的约1/105,所以,在图中光波直线应当十分陡峻,在图上难以区分,图中把它画的较倾斜是为了便于 ，，w0q,0，辨认。这种情况表明,与光波共振的将是的长光学波。实际晶体的长光学波的 131410,10/s在的范围,对应于远红外的光波。离子晶体中的光学波的共振能够引起对 w,w，远红外光在附近的强烈吸收,这是红外光谱中的一个重要的效应。正是因为长光 w，学波的这种特点,的格波支称为光学波。 19 青岛理工大学毕业设计论文 图3.4 电磁波和光学波的共振 20 青岛理工大学毕业设计论文 第四章 三维晶格振动的推导 4.1 一维多原子链问题的处理 我们处理一维多原子链的步奏基本上分为了: ,,,,, 模型运动方程试探解色散关系波矢q的范围波矢q的取值 由前面的推导我们知道: 晶格振动的波矢数目 =晶体的原胞数N,格波振动频率数目=晶体的自由度数,格波的支数=原胞内原子的自由度数。 一维单原子链,设晶体有N个原胞。原胞内原子的自由度数=1,即有一只格波;晶体的自由度数=N,即频率数为N。 一维双原子链,设晶体有N个原胞。原胞内原子的自由度数=2,即2支格波;晶体的自由度数=2N,频率数为2N。 4.2 建立三维模型和求解运动方程 一维单原子和一维双原子链的模型已经比较全面的表现了晶格振动的基本特征,现在我们来对比一维单双原子链的方法来说明三维晶格的振动。 mm12设三维无限大的晶体,每个原胞中有n个原子,各原子的质量分别为,,mml(lll)3n123,……,,原胞以标志,表明它位于格点: ,,,,(),,,Rl,l，l，l112233 ，4.1， 21 青岛理工大学毕业设计论文 原胞中个原子的位置用: ,,,lll,,,,,,,,,,,,R,R,？？,R,,,,,,12n,,,,,, ，4.2， 表示,偏离格点的位移则写成: lll,,,,,,,,,,,,,,,,,,,,？,,,,,,,n12,,,,,, ，4.3， 和双原子链一样,可以写成一个典型原胞中的运动方程: l，，m,(),？？k,k ，4.4， ,其中k表明原胞中的各原子,取值1,2,3,……,n。表示原子的三个位移分量,方程右端是原子位移的线性齐次函数。方程解的形式和一维完全相似,可以表示成: ,l,,,,i[wt,Rq]l,,,,,k,,,,,Ae,k,,k,, ，4.5， ,q指数函数表示各种原子的振动都具有共同的平面波的形式,是其波数矢量,也即 ,qA(A,A,A)A(A,A,A)q11x1y1z22x2y2z是波数,的方向是指波传播的方向。,,……都可以是复数,表示个原子的位移分量的振幅和位相可以有区别。公式，4.5，实际上可以表示为三位晶格格波的一般形式。 A,A,AA,A,A1x1y1znxnynz同样可以证明公式，4.5，代入式，4.4，中,得到以,……,为未知数的3n个线性齐次联立方程 22 青岛理工大学毕业设计论文 ,q2mwA,C()A',,,kka',k'k,kk, ，4.6， 2w(j,1,2,？？，3n)jw它的有解条件是的一个3n次方程式,从而给出了3n个解。具 ,w,qA,A,？？,Aq,0j12n体分析表明,当时,有三个解,从而对这三个解趋于相同, -3，也就是说在长波极限整个原胞一起移动。这三个解实际上与弹性波相合。另外，3n 个解的长波极限描述n个格子之间的相对振动,并具有有限的振动频率。所以在三维晶 ,有3个声学波,，3n-3，个光学波,或者说有3支声学波,，3n-3，格中,对于一定的波矢q 支光学波。 ,q4.3 波矢的取值和范围 ,q在三维情况中,一样受到边界条件的限制,只能取得某些特定的值而不是任意的。 ,,,qqq我们常引入所谓的“空间”来表示边界条件所允许的值,即是说,我们常常把看 ,,qq作是空间的矢量,而边界条件允许的值将表示为这个空间中的点子。“空间”是以,,,,b,b,bq123倒矢量为基矢,也就是可写成: ,,,,q,xb，xb，xb112233 ，4.7， 的形式。仍采用玻恩-卡曼边界条件,在三维情况下为: ,,,,,,,,,RNR(，),()l11l,,,,,,,,,RNR(，),(),l22l,,,,,,,RN,,R(，),()l33l, ，4.8， ,,,N,N,N,,,,,123123其中为晶格基矢,为三个基矢方向的原胞数,显然有晶体总原胞 ,,,R,(R)N,NNNll123数。代表格点上原胞的位移，可以是k中任何一类原子的位移，。 23 青岛理工大学毕业设计论文 ,N,ii便捷一傲剑表示,沿着方向,原胞的标语增加,振动情况必须相同，i=1,2,3，。 公式，4.8，边界条件要求: ,,h,1,,qN,h2,x,1111,N1,,h,,2,,qN,h2,x,,2222N2,,h,,3qN,,h2,,x,3333,N3, ，4.9， h,h,h123其中为整数。代入式，4.7，可得: ,,,hhh,123q,b，b，b123NNN123 ，4.10， ,,qq代表空间均匀分布的点子,每个点子占据的空间体积为: ,,, bbb1312(，),，倒格子原胞的体积NNNN123 考虑到倒格子原胞的“体积”与正格子原胞的体积之间的关系,可以得到边界条件 ,,qq允许的在空间均匀分布的密度: 1NvV0,,,,,33(2)(2),,bbb312()，NNN123 分布密度= (4.11) V为晶体的体积。 ,q从晶格的振动考察,的作用只是在于确定不同原胞之间振动位相的联系。我们把 ,q的取值范围选为上述倒格子原胞并不是最方便的,我们通常是选在在第一布里渊区 24 青岛理工大学毕业设计论文 内。可以证明第一布里渊区的体积等于原胞的体积,且第一布里渊区具有环绕原点更为对称的优点。 V,3q(2,)的分布密度为,因此不同的q的总数应当是: 由于边界条件所允许的 V(倒格子的原胞体积)，,N3(2,) ,q和晶体中包含的原胞数目相同。对于每个有3个声学波,，3n-3，个光学波,所以不同的格波的总数: N(3，3n,3),3nN 正好等于晶格Nn个原子的自由度。这表明,上述的格波已经概括了晶格的全部振动模。 4.4 理论上的计算 ,,，，wqqj作为的函数称为晶格振动谱,或称为格波的色散关系。它可以通过实验的方法测量得到,也可以根据原子间相互作用力的模型从理论上进行计算。由理论与实验的比较重获得对相互作用力的认识。共价晶体、离子晶体、金属晶体、分子晶体等由于它们的原子间互相作用力有着不同的特点,因而在格波的谱上也有相应的特征。下面我 ,,qq们举几个典型的例子。三维晶格中是矢量,在作图时总是固定的方向,一般选典型 ,,，，wqqj的对称轴方向。分别画出沿不同方向时的变化。三维晶格还需要考虑原子位移 ,,q2方向与格波传播方向之间的关系,若沿着晶体的一个对称轴,晶体绕这个轴转,或,2,33、是对称操作时,这时格波可以分为纵波和横波,纵波原子位移平行于波的传播 25 青岛理工大学毕业设计论文 方向;横波原子位移垂直于波的传播方向,而且包括两个频率兼并的波。通常用TA表示横声学波,LA表示纵声学波,TO表示横光学波,LO表示纵光学波。 图4.1中给出了硅的格波谱,由于金刚石结构中每个原胞包含两个原子,因而存在 ,TO是两重兼并的，,可以看出长声学波极限纵波与横波有纵光学波和横声学波。，TA 不同的波速,长光学波极限,纵波横波有相同的频率。 图4.1 硅的格波谱 图4.2给出了GaAs的格波谱,GaAs是错误!未找到引用源。III—错误!未找到引用源。V族化合物,具有闪锌矿结构,它的格波谱与si很相似,主要区别在于q=0时纵光学波与横声学波的频率是不同的,这是离子性的反应,离子型越强,两个频率只差越大。 26 青岛理工大学毕业设计论文 图4.2 GaAs的格波谱 图4.3中给出了金属Pb的格波谱,由于它具有面心立方简单晶格结构,只有声学波。 ,,,qq，，wq在图中某些值附近函数出现扭转,这是因为对于这些值的格波与金属中电子之间耦合特别强的结果,科恩1959年曾预言与此相关的效应,故称为科恩异常。 图4.3 Pb的格波谱 27 青岛理工大学毕业设计论文 第五章 结论和讨论 晶格振动在在研究固体物理宏观微观性质方面有着重要的作用,可是晶格振动是一个多原子的复杂运动问题,不好直接研究说明。所以我们通过研究原子链的色散关系来形象说明晶格振动的一些具体性质和在固体物理学中的重要作用。通常一维双原子链的色散关系能更加形象的让我们了解晶格振动,可是我们首先要先了解一维单原子链的色散关系,只有这样我们才能更加容易一维双原子链的色散关系。最后我们在通过研究三维晶格的研究,把我们的研究方向由一维增加到多维,这样更贴近实际,能更好地了解晶格振动。 关于一维双原子的色散关系的论文终于完成。除了物理知识的提高和研究能力的进步外,作者获得的一个重要的收获是:无论做什么事情,只有真正用心地、去持之以恒地做,才能有成功。具体来说,如正在做的毕业论文研究,要将这个事情做好、做完整,需要锻炼培养以具备系统的思维方式和方法,并且要善于运用一切可以利用的资源来充实自己,对待要解决的问题要耐心、坚持不懈,对待问题要迎难而上,如果一个方法不行可以换个角度和去尝试。这一点,会帮助我们有信心有准备地进入新的研究领域。 28 青岛理工大学毕业设计论文 致 谢 四年的读书生活在这个季节即将划上一个句号,而于我的人生却只是一个逗号,我将面对又一次征程的开始。四年的求学生涯在师长、亲友的大力支持下,走得辛苦却也收获满囊,在论文即将付梓之际,思绪万千,心情久久不能平静。 伟人、名人为我所崇拜,可是我更急切地要把我的敬意和赞美献给我的导师。我不是你们最出色的学生,而你们却是我最尊敬的老师。 首先感谢指导我的毕业论文工作的陈冰泉老师、刘彦平老师。在这几个月里,得到了两位导师细心的指导,在这个过程中,使我学到了丰富的专业知识,并锻炼了独立进行科学研究工作的能力。老师渊博的学识、严谨的治学态度、敏锐的洞察力、活跃的学术思想、积极进取的敬业精神、雷厉风行的工作作风和为人坦诚的思想品德,使我受益非浅,终生难忘。同时也感谢所有的任课老师,谢谢你们的辛勤付出,是你们的敦敦教诲才让我获得扎实的理论基础。 授人以鱼不如授人以渔,置身其间,耳濡目染,潜移默化,使我不仅接受了全新的思想观念,树立了宏伟的学术目标,领会了基本的思考方式,从论文题目的选定到论文写作的指导,经由老师悉心的点拨,再经思考后的领悟,常常让我有“山重水复疑无路,柳暗花明又一村”。 感谢我的爸爸妈妈,焉得谖草,言树之背,养育之恩,无以回报,你们永远健康快乐是我最大的心愿。在论文即将完成之际,我的心情无法平静,从开始进入课题到论文的顺利完成,有多少可敬的师长、同学、朋友给了我无言的帮助,在这里请接受我诚挚谢意! 29 青岛理工大学毕业设计论文 感谢徐玉国,柴飞等同学在学习和研究上给予我的帮助和支持。 同时也感谢学院为我提供良好的做毕业设计的环境。 最后再一次感谢所有在毕业设计中曾经帮助过我的良师益友和同学,以及在设计中被我引用或参考的论著的作者. 30 青岛理工大学毕业设计论文 参考文献 [1] 黄昆 韩汝琦,半导体物理基础,科学出版社,2010 [2] 刘恩科 朱秉升 罗晋生等,半导体物理学,国防工业出版社，第四版，,2008 [3] 黄昆原著 韩汝琦改编,固体物理学,高等教育出版社，第一版，,1988 [4] Stephen J.Chapman,MATLAB Programming for Engineers,科学出版社,2003 31