.固体能带计算方法.pdfVIP

第 23卷 第 1期 贵州大学学报( 自然科学版) Vol. 23 No. 1 2006年 2月 Journal ofGuizhou Un iversity ( Natural Sciences) Fe . 2006 1000- 5269(2006) 01- 0068- 05 * * 闫万 , 谢 泉, 朱林山, 金石声 (贵州大学 计算机科学与工程学院, 贵阳 550025) , Si Ge, 。 ;; O48111 A 0 早期, 人们对于晶体中 子能带结构的了解只是定性的, 随着科技的发展情况就完全不同了。一方面 由于纯单晶的制备, 低温、微波等技术的发展, 人们可以通过回旋共振、磁阻、磁光、光谱等实验来直接测定 若干晶体的 子能带结构; 另一方面, 由于 子计算机的广泛应用, 在理论上也可以对能带结构进行更精 [ 1, 2] 确的计算。能带理论在材料的研究中已经成为非常有用的工具 。 对于三维的周期场中的单 子问题, 通常采用各种近似方法求解, 首先, 选取某个具有布洛赫函数形 式的完全集, 把晶体 子态的波函数用此函数集展开, 然后代入薜定谔方程, 确定展开式的系数所必须满 足的久期方程, 据此可求得能量本征值, 再依照逐个本征值确定波函数展开的系数。选择不同的函数集 合, 有不同的计算方法。借助快速大容量的 子计算机, 使实际的计算更快捷。本文将对这些常用的计算 方法进行论述, 并采用第一性原理分子动力学方法对 Si和 Ge的能带进行了计算。 1 111 (LCAO) [ 3] 紧束缚方法(Tight2Banding)又称为原子轨道的线性组合法, 简称 LCAO方法。第一次由 F. B loch在 1929年推出的, 其中心思想就是用原了轨道的线性组合来作为一组基函数, 由此而求解固体的薛定谔方 程。这个方法是基于这样的物理图像, 即认为固体中的 子态与其组成的自由原子态差别不大。按 Bloch 观点, 当晶体 子离某个原子较近时, 原子势场对 子有较强的束缚作用, 此时 子又有共有化运动, 它有 等同的几率出现在每个原子范围, 就是说, 子的共有化运动带有强烈的原子态的色彩。因此就用原子轨 道的线性组合来构成 子的共有化轨道, 即 Bloch 函数。 紧束缚方法在绝缘体的能带结构研究中是很成功的。由于原子轨道处在不同的格点上, 由它们组成 的基函数一般是非正交和。因此, 必然会遇到中心积分的计算问题, 而且本征方程形式也不简便。为了克 服这些困难, 人们道过 Slater2Koster参量法, W oAoHarr ison键轨道近似以及 OlcAo方法解决。 112 11211 平面波方法 平面波方法是计算方法中最简洁, 物理定义较明显的晶体 子能带计算方法, 其方程形式为 2 õ( k + K ) - E ( k) 8D + V(K c- K ) c ( k, K ) = 0 E n KKc n K * 收稿日期: 2005- 11- 11 基金项目: 贵州省教育厅重点基金( 05JJ002), 贵州大学人才引进基金资助项目( 04RCJJ001), 湖南省自然科学基金( 04JJ3030), 湖南 省教育厅青年基金资助项目( 03B006), 国家人事部留学归国人员择优资助项目( 200403) 。 作者简介: 闫万 ( 1978- ), 女, 贵州大学计算机科学与工程学院研究生, 研究方向: 能带计算。