能带及计算方法.doc

第三章 能带的计算方法 周期场中的单电子波动方程除了少数几种简单的理想模型外，都只能用近似方法求解。目前，主要的近似方法有：准自由电子近似，紧束缚近似，原胞法，正交化平面波法，赝势法和法等。每一种近似方法都有其优点，也有其局限性，只能用于一定的情况。在这一章中简单介绍两种。 §3-1准自由电子近似法 在这种近似方法中假设原子的外层电子在晶体的周期性势场中运动，且势能的周期性变化部分很小，可作为微扰来处理。这种处理，电子的运动一方面和自由电子相近，另一方面又能反映出周期场中运动的电子所具有的周期性特征。这种方法较粗糙，适用于金属中的电子。 一．一维情况 设周期为a、长度为L的线状晶体沿x方向。电子波动方程为 （3-1） 式中， （为任意倒格矢）具有晶格的周期性，V0是电子在晶体中的平均势能。由于V(x)为实数，故有 令：W(x)为势函数中周期性变化部分，则 （3-2） 于是波函数可改写为 （3-3） 根据准自由电子近似的基本假设，W(x)很小，可当作微扰。从而可先求解无微扰的电子波动方程 （3-4） 其解为平面波 （3-5） 相应的能量谱值 （3-6） 这里，k是平面波的波矢量。在周期性边界条件下，k只能取断续值： ， 这些满足周期性边界条件的平面波彼此正交并归一化 （3-7） 当存在周期性变化的微扰W(x)时，波动方程的零级能量谱值为E0(k)。下面分两种情况讨论。 1．非简并情况。选择零级近似波函数为平面波，从而根据量子力学公式，能量一级修正项为 =0 （3-8） 故能量一级修正为零。进一步计算需考虑微扰矩阵元 （3-9） 故能量谱值的二级修正为 （3-10） 波函数的一级修正为 （3-11） 从而，考虑到二级近似后的能量为 （3-12） 考虑到一级近似后的波函数为 ]=， ] （3-13） 2．简并情况。当波矢量变化到使（3-10）式求和号中某一项的分母等于零或接近于零时，则该项所占的比例就会很大，不能再被认为是修正项了。这时，非简并化微扰理论就不再适用，需采用简并化微扰理论处理。在这种情况下，必须把能量彼此相近且矩阵元的平面波同时包含在零级近似波函数中。波矢量变化到使（3-10）式求和中第n项分母为零的条件为 ，即 （） 该条件正是确定布里渊边界的条件。当变化到布里渊边界附近时，零级近似波函数应该把波矢量为的和的平面波同时包含进去。即 （3-14） 如果忽略二级小量，则将零级近似波函数代入波动方程后，该式应近似地成立。于是有 （3-15） 先后用乘以上式并在L内积分，则有 （3-16） 上式为决定A和B的联立线性齐次代数方程组。要使A、B有不为零的解，其系数行列式必为零，即有 （3-17） 上式是关于能量E的久期方程，其解为 （3-18） 当k变化到布里渊区边界附近时，则存在 ，此时，（3-18）式可简化为 （3-19） 若令： ，，则由 和 得 （3-20-1） 和 （3-20-2） 其中，一支为上弯抛物线，另一支为下弯抛物线。在布里渊区边界上处，上弯抛物线的极小值为 （3-21） 下弯抛物线的极大值为 （3-22） 两者间能量间隙为。在此能量范围内，没有允许的能级存在。分别将上两式代入（3-16）式中。便有 。设，则。将此式代入（3-14）式，得波函数 上式括号中取正号得 取负号得 在布里渊区边界上，波函数为两个驻波，与相对应的驻波能量较高，与相对应的驻波