晶格振动模式密度研究王晴晴;宫昊;程荣龙;葛立新【摘要】晶格振动模式密度(声子态密度)即单位频率间隔内的模式数,是反映声子在波矢空间分布疏密程度的物理量.为了准确地求出品格热容量随温度的变化关系,必须用较精确的办法计算出晶格振动的模式密度,进而掌握材料的热力学性质.一般教材中对该部分的讲解晦涩难懂,本文从品格振动的物理意义开始,分析说明并推导一维、二维、三维不同体系的晶格振动模式密度公式,进而求出德拜模型下不同体系晶格热容公式.【期刊名称】《大学物理》【年(卷),期】2018(037)008【总页数】4页(P4-7)【关键词】晶格振动;声子态密度;德拜模型;晶格热容【作者】王晴晴;宫昊;程荣龙;葛立新【作者单位】蚌埠学院理学院,安徽蚌埠233030;蚌埠学院理学院,安徽蚌埠233030;蚌埠学院理学院,安徽蚌埠233030;蚌埠学院理学院,安徽蚌埠233030【正文语种】中文【中图分类】O481晶格振动模式密度，也称为频谱密度、声子态密度、简正模式密度、频率分布函数等，是在解决晶格热容随温度的变化关系过程中提出来的，一般用g(ω)表示，黄昆先生的《固体物理学》[1]教材中给出的定义式为(1)其中dn表示在ω→ω+dω频率间隔内晶格振动模式的数目．最初对晶格振动模式密度的计算是为了准确求出晶格热容随温度的变化关系，而晶格热容的大小直接反映了固体材料的热力学性质．后来发现，计算晶格振动模式密度对研究晶体的某些电学性质和光学性质也非常有用．但是教材中对于晶格振动模式密度的计算这部分讲解晦涩难懂，特别是对于一般本科院校的学生很难完全理解如何计算和应用．本文将以通俗的语言从晶格振动的物理意义开始，引入一维、二维、三维不同体系格波分布密度和等频率面间的体积，然后举例讲解如何计算晶格振动模式密度，并获得计算通用公式，最后在德拜模型下推导不同体系的晶格热容公式．1 晶格振动模式晶体中原子的排列是具有周期性的，而其排列的具体形式即称为晶格，格子的位置代表了原子的平衡位置．晶体的一些宏观性质如导电性、导热性等就是因为晶体中电子和原子的运动造成的，而晶格振动即晶体中原子在其平衡位置附近的振动，温度越高振动越激烈，温度达到熔点附近则原子将脱离格点的束缚而发生相变．在简谐近似即温度不太高的小振动过程中，研究发现晶格振动中所有原子都以相同的振动频率参与振动，且振动的模式具有波的形式，一般把体系中所有原子一起参与的共同振动叫做一个振动模式．振动模式具有波的形式又叫做格波，格波的能量量子叫声子，声子是一种假想的能量量子．原则上只要找到了体系的振动模式或格波的运动方程，就能方便的求解薛定谔方程，解决相应的问题．振动模式具有波的形式，需要描述体系格波的特性，如格波方程、格波波长、格波波矢等参数．由于晶格具有周期性，只需讨论一个周期内的特性就可反映整个晶体的特性，因此引入了布里渊区的概念．第一布里渊区即倒格子空间的维格纳-赛兹原胞，是倒格子空间的最小周期性单元．在求解晶格振动模式密度的过程中比较关注的就是格波波矢q与频率ω之间的关系即色散关系式ω(q)(又称为晶格振动谱)[2]．研究发现一维单原子链一个波矢q对应一个频率ω，而二维、三维体系中一个波矢q对应多个频率ω，即一个波矢将对应多个振动模式．2 晶格振动模式密度的计算2.1 格波分布密度不管是一维、二维还是三维体系都存在边界问题，在求解所有原子的联立方程过程中，很难求出格波解的形式，因此所有体系都要在玻恩-卡曼边界条件下进行，即其中N=N1·N2·N3为总的原胞个数，α1、α2、α3为晶格的三个基矢，h1、h2、h3为一组整数．即波矢q不是连续的，而是q空间均匀分布的量子化的点，点的个数正好为h1·h2·h3个．格波分布密度就是单位体积中的量子化点的个数，所以不同体系的格波分布密度计算公式如下．1) 三维体系q空间的体积为格波分布密度为其中V0=α1·(α2×α3)为原胞的体积，V为晶体的体积．2) 二维体系q空间的面积为格波分布密度为其中A0=|α1×α2|为二维体系原胞的面积，A为二维晶面的面积．3) 一维体系q空间的长度为格波分布密度为其中α为一维体系原胞的间距，L为一维原子链的长度．2.2 晶格振动模式密度计算特例根据晶格振动模式密度的定义，其中dn表示ω→ω+dω频率间隔内晶格振动模式的数目，若能求出ω→ω+dω等频率面(ω(q)=常数对应的q空间中的面)间的体积，则dn=格波频率密度×等频率间体积1) 三维体系，以等频率面是球面为例．图1 三维体系等频率面等频率面间的体积为其中忽略(dq)2项和(dq)3项．ω→ω+dω频率间隔内晶格振动模式的数目为晶格振动模式密度为(2)若可以通过实验的方法获得体系的晶格振动谱ω(q)，就可直接求出g(ω)的值．2) 二维体系，以等频率线是圆形为例．等频率线间的面积为π(q+dq)2-πq2≈2πqdq其中忽略(dq)2项．ω→ω+dω频率间隔内晶格振动模式的数目为晶格振动模式密度为(3)3) 一维体系等频率点间距可直接为2dqω→ω+dω频率间隔内晶格振动模式的数目为晶格振动模式密度为(4)2.3 晶格振动模式密度一般表达式实际三维晶体的等频率面一般很复杂，不能用简单的球形来代表．若用dS表示等频率曲面上的面积元，dq为q空间频率ω→ω+dω的垂直距离，则表示第j个模式等频率面的梯度．(假设体系包含n个原子，则三维体系共有3n支格波，其中3支声学波，(3n-3)支光学波；二维体系共有2n支格波，其中2支声学波，(2n-2)支光学波；一维体系共有n支格波，其中1支声学波，(n-1)支光学波．)晶格振动模式密度的一般表达式如下：1) 三维体系：∬2) 二维体系：3) 一维体系：3 实验与理论计算方法求解晶格振动模式密度(声子态密度)的关键在于晶格振动谱的具体表达式．实验上一般采用Raman散射获得晶格振动谱：即将特定频率的光波照射待测样品表面，光子会与声子发生碰撞并散射，若能测定散射前后粒子的频率与波长的改变，就可以根据能量守恒和动量守恒公式确定声子的频率和波矢的关系．另外也可以采用布里渊散射、中子非弹性散射、X射线散射等实验手段确定晶格振动谱．对于晶格振动谱比较简单的体系，可以近似的描述等频率面的形状并求出等频率面间体积，直接用上文中的公式计算声子态密度．对于晶格振动谱复杂的体系，不能单纯的通过公式计算求解，可用Raman振动谱的积分面积和拟合峰高来定性描述声子态密度，若积分面积和拟合峰高比较大，说明声子态密度大，反之则声子态密度小[3]．理论模拟计算中，一般采用基于密度泛函理论的第一性原理计算分析[4-6]，如利用Materials Studio中的CASTEP程序软件包、Quantum-ESPRESSO软件的PWscf软件包均可以画出声子态密度曲线．4 德拜模型下由声子态密度求解晶格热容材料的热力学性质与其内部的晶格振动有很大关系，基于密度泛函理论的第一性原理可通过较小的计算量获得布里渊区内声子谱，从而求出相应的热力学参数．德拜模型提出将格波看成弹性介质波(长声学波)，而热容是原子的各种频率(0到极大值ωD，ωD为德拜频率)振动贡献的总和．即德拜模型下三维体系包括一支纵声学波、两支横声学波，二维体系包括一支纵声学波、一支横声学波，一维体系包括一支纵声学波．弹性波色散关系式可简单写成：纵波：ωL=vLq；横波：ωT=vTq(5)其中vL、vT分别为纵声学波和横声学波的传播速度.将式(5)分别带入式(2)、式(3)、式(4)可得声子态密度．三维体系：二维体系：一维体系：假设体系包含N个原胞，根据式(1)即声子态密度的定义式：1) 三维体系g(ω)dω=3N⟹(6)2) 二维体系g(ω)dω=2N⟹(7)3) 一维体系g(ω)dω=N⟹(8)根据晶格热容的量子理论，晶体内能和晶格热容公式分别为：(9)(10)分别将式(6-8)带入式(9)、式(10)计算得：1) 三维体系2) 二维体系3) 一维体系其中为德拜温度，表征了原子间结合力的一重要物理量，不同材料的德拜温度不同，熔点高，即材料原子健结合力强，则德拜温度越高[7,8]．根据德拜温度的定义它是一个跟温度无关的物理量，但实验证实德拜温度随着温度的变化而变化[9]，这因为德拜模型忽略了晶体的各向异性及光学波、高频声学波对热容的贡献[10]．实际晶格热容随温度的变化曲线直接反映了声子态密度的整体行为[11]．而热容是固体原子热运动在宏观性质上的最直接体现．因此研究和分析实际固体的声子态密度具有重大的意义．【相关文献】[1] 黄昆. 固体物理学[M]. 韩汝琦,改.北京：高等教育出版社, 1988.[2] 田强, 洪馥男. 具有在位势的一维双原子链晶格振动的色散关系[J].大学物理, 2006, 25(4):17-17.[3] 苏方宁, 邓再德. 声子态密度对掺稀土玻璃材料上转换发光强度的影响[J]. 玻璃与搪瓷, 2007,35(2):1-5.[4] 郭连权, 林琳, 马贺,等. Mg2Si的声子谱与力热性能的第一性原理计算[J]. 沈阳工业大学学报, 2015, 37(3):294-298.[5] 蒋文灿, 陈华, 张伟斌. TATB晶体声子谱及比热容的第一性原理研究[J]. 物理学报, 2016,65(12):216-224.[6] 付佳琦. 第一性原理研究纤锌矿结构AIN和InN的声子及热力学性质[D]. 内蒙古大学, 2015.[7] 程本培, 虞炳西. 用测定德拜温度的方法判断合金中原子间结合力的变化[J]. 物理测试,1991(5):6-11.[8] 文潮, 孙德玉, 关锦清,等. 用X射线衍射强度测定纳米金刚石的德拜特征温度和熔点[J]. 高压物理学报, 2003, 17(3):199-203.[9] 刘洋. 纳米材料德拜温度、体膨胀系数及热容的尺寸效应[D]. 吉林大学, 2008.[10] 高钦翔, 田强. 光频支格波对晶格热容的贡献[J]. 大学物理, 2002, 21(7):16-17.[11] 王晓旭, 赵琉涛, 成海霞,等. La3Co29-xFexSi4B10的择优占位、电子结构和晶格振动性质的理论研究[J]. 物理学报, 2016, 65(5):266-272.。
