固体能带计算方法探讨.pdfVIP

第 ll卷 第 2期 丹 东 纺 专 学 报 Vo1． 1l No． 2 2004年 6月 JOURNAL OF DANDONG TEXTILE COLLEGE June2004 固体能带计算方法探讨 尤凤翔 龚善初 陈文涛 (揭阳学 院，广东 揭阳 522000) 摘 要：本文分析 了固体能带的计算方法，并用紧束缚方法采用Stater型和Gaussian型的最佳轨道。计算了闪锌矿结构 半导体化合物电子能带，并进行了分析。 关键词 ：固体能带；紧束缚方法；闪锌矿结构 中图分类号 ：0481．1 文献标识码 ：A 文章编号 ：1008—2174(2004)02—0075—03 固体能带理论无疑是凝 聚态物理 中最成功的理论之一， 然会在芯区引进振荡 的成分 ，这恰好能描写导 电电子的特 是固体 电子论的支柱 。固体的许 多基本物理性质 ，如振动 征。因此，这种同心层电子态正交 的平面波称为正交化平面 谱、磁有序、电导率、热导率、光学介电函数等，原则上都 波 (opw) 可由固体的能带理论阐明和解释，或需要对具体材料的能带 里 一 I= e：③ 结构有所了解 。因此 ，确定 固体电子能级 ，也就是能带的计 正交化平面波 … 一 算是十分重要的 ，( 、 1 对于三维的周期场中的单电子问题 ，通常采用各种近似 · ) 了 ““ 喾。 ) 方法求解 ，首先 ，选取某个具有布洛赫函数形式的完全集 久期方程为 ： 合，把晶体电子态的波函数用此函数集合展开 ，然后代入薜 定谔方程，确定展开式的系数所必须满足的久期方程 ，据此 detl IH I 一E I l=0 ④ ★+ I+t ★+★ ★+ 可求得能量本征值 ，再依照逐个本征值确定波函数展开的系 数。选择不同的函数集合 ，有不同的计算方法。借助快速大 OPW方法 已被广泛用于 Ⅳ (金刚石结构)，Ⅲ——V族 容量的电子计算机 ，使实际的计算更快捷。实际中，一些半 和 Ⅱ——Ⅵ族 (vr方闪锌矿结构和六角纤维锌矿结构)半导 径的经验方法也起相当重要的作用，它借助实验数据来确定 体的能带计算。图 1是 F，Herman等 (1)用 OPW方法得到 计算较困难的量。本文对这些计算方法进行阐述 ，并通过紧 的Ge(锗)和 si(硅)的能带结构 。计算时用 了 113个 束缚方法计算闪锌矿结构半导体化合物电子能带。 1 计算方法 1．1 平面波方法 ：平面波方法是计算方法中最简洁 ，物理 定义较明显的晶体电子能带计算方法，其方程形式如下 ： 【( ‘)一 (+ (K。’v一 )0(【‘)=。① 久期方程为上式的系数行列式必须等于零。即 detIf三一(七+七)一E(七)】sI．．七．+V(K．一K．)l=0 旧 二，，I 图1 Ge (al和 SiIbl的能带 ②式为无限阶的行列式，实际计算只能取有限阶的行列 - ．