56能带结构计算方法：.pdfVIP

5.6 能带结构的计算方法： 显然，NFE和TBA用于计算可以与实验结果比较的实际 能带是太粗糙了。已经发展了许多计算实际晶体能带的模型 和方法，这些方法既需要比较深的量子力学基础，又需要大 量繁琐的数学运算。近代的能带计算多使用大型计算机，采 用建立在密度泛函理论基础上的局域密度近似，但早期的几 个模型均可用来做密度泛函计算，所以这里我们简要的定性 地介绍一些曾获得一定成功的模型和方法。 参见：阎守胜《固体物理基础》3.4节; 李正中《固体理论》7章; Ashcroft 《Solid State Physics》11章 冯端《凝聚态物理学》上册12章 Kittel 8版9.3节p163-169 一. 原胞法（Cellular Method）Wigner-Seitz 1933 二.缀加（增广）平面波法（Augmented 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原胞中电子的影响可以忽略不计。  2 2    U r  r E k  r           k k  2m  只要求出一个原胞中的波函数就可以把整个晶体的问题解决 了 （平均地说，每个原胞都被一个传导电子所占据，这些电子 往往有屏蔽离子的作用，从而强烈地消弱了离子势场。） 这是一个近似图，并不准确。  r a   0 0 原胞法示意图 ik r  r e u r 根据Bloch定理，方程的解为： k   k   其中表征周期部分的函数u(r) 应该在原胞的两个对称点上 例如p 和p 点上取值相同。为了方便地处理周期性的要求， 1 2 Wigner和Seitz用相同体