能带论4课件.pptVIP

* 显然，NFE和TBA用于计算可以与实验结果比较的实际能带是太粗糙了。已经发展了许多计算实际晶体能带的模型和方法，这些方法既需要比较深的量子力学基础，又需要大量繁琐的数学运算。 近代的能带计算多使用大型计算机，采用建立在密度泛函理论基础上的局域密度近似，但早期的几个模型均可用来做密度泛函计算，所以这里我们简要的定性地介绍一些曾获得一定成功的模型和方法。 参见：阎守胜《固体物理基础》3.4节; 李正中《固体理论》7章; Ashcroft《Solid State Physics》11章 冯端《凝聚态物理学》上册12章 Kittel 8版 9.3节p163-169 6.6 能带结构的计算方法 一. 原胞法（ Cellular Method）Wigner-Seitz 1933 缀加（增广）平面波法（Augemented Plane Wave Method ）Slater 1937 三. 格林函数法 (Korring1947,Kohn Rostoker 1954) 四. 正交平面波法（Orthogonalized Plane Wave Method ） Herring 1940 五. 赝势法（Pseudopotentials） Harrison 1966 六. 密度泛函理论（The density-function theory） Wolter kohn 1960 一. 原胞法（ Cellular Method） 这是能带计算最早使用的方法，1933 年由Wigner和Seitz引入。它曾成功地用于碱金属，特别是Na和K。第一章中使用的维格纳-赛兹原胞就是此时提出的。例如金属 Na 具有 bcc结构，它的 W-S 原胞是一个截角八面体，离子在其中心位置。晶体（110）晶面示意图如下页。 现在假定电子处在 A 原胞中，求解运动方程时，可以认为它只受到此原胞中离子势场的影响，其它原胞中离子势场对A原胞中电子的影响可以忽略不计。 只有求出一个原胞中的波函数就可以把整个晶体的问题解决了（平均地说，每个原胞都被另一个传导电子所占据，这些电子往往有屏蔽离子的作用，从而强烈地消弱了离子势场。） 原胞法示意图 这是一个近似图，并不准确。 根据Bloch定理，方程的解为： 其中表征周期部分的函数 u(r) 应该在原胞的两个对称点上 例如 p1和 p2 点上取值相同。 为了方便地处理周期性的要求， Wigner和Seitz用相同体积的WS球代替实际原胞，原胞内的势 场也看成球对称的，因而可以数值求解，在 k＝0 点得到的波 函数 划在图中。在带底附近，k 的另外一些点，波函数 近似表示为： 同样也可求出E0 ,和 Wigner和Seitz 用这种方法得到的能量去计算简单金属的结 合能，其结果令人满意地与实验一致。见陈洗书p340 这是原胞法求出的 曲线（实线），可以看出波函数在离子实内是振荡的，而一旦离开离子实部分，就基本是常数。波函数的这个常数部分几乎占原胞体积的90%，因此在晶体中波函数基本是一个平面波。电子在晶体中的运动基本是自由的，所以Na的导电电子是自由电子。 和原子波函数（虚线）相比，变平是由于加上边界条件产生的，而不是离子势场有什么特殊的性质，这个结果对以后的能带计算有启示。 实际势函数 近似势的等能面 这是 Ashcroft 书对原胞方法的图示p196 原胞法的球形近似，所得结果仅与原子体积有关，完全忽略了实际晶体结构的影响，实在过于简单。而使用实际的复杂的原胞形状，又给计算带来难以克服的困难，所以近来已经很少使用，改用下面几种方法。 二. 缀加（增广）平面波法 （Augemented Plane Wave Method ） 1937年Slater 从原胞法结果中发现晶体离子之间的广大区域里，晶体势几乎没有变化，因此采用了一种特殊的势重新求解。“muffin-tin” （糕模势）在离子实处就是自由离子势，在离子实外面严格是一个常数（见下页图）。因此波函 数为在 时是平面波： 而在 时，为原子波函数。r0 是离子实半径。后者可以通过求解自由原子的薛定谔方程得到，边界条件是在离子实的表面 r0 处和平面波连续衔接。 上面给出的函数，可能不满足Bloch定理，但可以通过线性组合来弥补： Ashcroft 书 p200对“The