固体能带理论I2

1、4.2原子轨道线性组合紧束缚方法(tight-binding,TB)第一次由F.Bloch在1929年提出，其中心思想就是用原子轨道的线性组合(Linearcombinationofatomicorbitals,LCAO)来作为一组基函数，由此而求解固体的薛定调方程。这个方法是基于这样的物理图像，即认为固体中的电子态与其组成的自由原子差别不大。紧束缚方法在绝缘体的能带结构研究中是很成功的。由于原子轨道处于不同的格点上，由它们组成的基函数一般是非正交的。因此必然会遇到多中心积分的计算问题，而且本征方程形式也不简便。1紧束缚方法考虑固体中单电子的薛定调方程：Hnkr-茄乂Vr-nkr：尸EnJnk

2、r(4.2.1)式中哈密顿量的第一项是电子的动能，第二项是晶体势场；Ek是第n个能带且具有动量k的能级；中仙nknk描述固体中电子的波函数。at晶体势场可以表述为原子势场V(r)的线性叠加，即Vr='、Vatr一R-t.(4.2.2)l、工这里R是晶格矢量，t口是第l个原胞中第a个原子的位矢。TB方法的中心思想是利用原子轨道的线性组合作为基矢，即波函数中nk可用LCAO的基矢%来展开k")二，Anjjkr(4.2.3)这里的布洛赫函数r)由原子轨道线性组合：1jkrejr-Rl7工)(424)Nl,:.式中W(r-R-tu)第l个原胞中第a个原子的第j个轨道，N是单位体积的晶

3、格数目。值得注意的是，在同一格点上的原子轨道是相互正交的，但相邻原子间的轨道函数却一般是非正交的，因此%一般是非正交的。人)是线性组合参数，由解本征问题而得到。将方程(4.2.3)带入方程(4.2.1)并和弧(rX乍内积，得到£Anj他,kH©k)=Enk£Anj僧k|©k)(4.2.5)jj定义上式则可简化成这里Hj,j为哈密顿量的矩阵元,AjHj,j-EnkSj,j=0Sj,j为原子轨道交叠积分。(4.2.6)(4.2.7)为求展开参数%的非零解，需要解如下的本征方程:(4.2.8)detHj,j)=0上式即为TB方法的出发点及原始形式。在解本征方程

4、时，我们会遇到两个困难:(1)多中心积分。在计算矩阵元Hj,j和Sj时，会遇到多中心积分问题。例如：Hj,j中包含如下的两中心和三中心积分：僧,(-%W忒«N')(429)'Yr-t:Vatr-t：,r-t!：;(4.2.10)严格计算这些多中心积分是非常困难和繁琐的，因此人们通常忽略三中心积分而只考虑两中心积分。它们通常的形式为a'Vatr-tatr-t(4211)(2)本征方程的形式不简便，除了对角项外，非对角项也包含有Enko这是因为基矢函数是非正交的，使的非对角元5口为非零。在从头计算中，矩阵元Hj,j和Sj对给定的哈密顿量H和原子轨道要作严格的求解。

5、一般讲，这些矩阵元在实空间中收敛很慢，因此计算量相当大，几种近似方法由此产生，如Slater-Koster参数法和键轨道近似。2Slater-Koster参量方法从量子力学知道，原子轨道波函数可以写为(4.2.12)(4.2.13)对l=0,1,2,3,4的态，分别称为sP,(4.2.14)我们先讨论一下原子轨道nlm=RYm其中n为主量子数，l是表征角动量大小的量子数，而m二，1,，l-1,lRnl(r共轨道波函数的径向部分；Ylm为球谐函数，仅与角度有关。d,f,g,态。对s轨道，l=0,m=0,这个态是球对称的，其中球谐函数为.1丫10二五对p轨道，l=1,m=-1,0,1,相应的球谐函

6、数为丫1，=sinie'4丫1。(旦甲产(4.2.15)3sin二8二如果采用笛卡儿坐标，可以将球谐函数线性组合，将球谐函数表示为x,y,z的函数，这样轨道函数rnlmnlIxrry.rzr(4.2.16)分别记为Px，Py,Pz轨道。对p轨道，l=2,m=2,-1,0,1,2,在采用笛卡儿坐标中其函数形式为'yz/r2zx/r2，:LI242m(r)=;!5Rn2(r乂(4.2.17)一、4n-(x2.y2y(2r2)(3z2-r2)«2如2)现在以s和p轨道为例来讨论原子间的相互作用矩阵。原子A的和B的p)轨道之间的相互作用作如下标记(slHIs)三Vs灯(4.2

7、.18)j(s1Hp)=%"gH|Px)appQppyHPy)=pzHPz=Vppn上面是两原子间相互位矢d=品-ra平行于y轴的情况。对于一般情况，原子间相互作用可以分解成七的线性组合，其线性组合参数与位矢d的方向余弦l,m,n有关。下式中给出了用V表示的s,p,d轨道间相互作用矩阵元，取自文献Ess=sHs)Esx=(sHPx)=1Vspq-Exx=PxHPx=l2Vpp；：+1-l2Vpp二(4.2.19)Exy=(pxHPyj=lmVpp。ImVppnExz=PxHPz)=lnVpplnVppnJ.C.Slater和G.F.Koster在1954年的一篇经典文章中提出了一个非

8、常有价值的参量方法，他们建议将LCAO的哈密顿量矩阵元看成参量，其大小由布里渊区中心或边界上高对称k点的精确理论值或试验值拟合而得。为了避免多中心积分的计算，他们用Lwdin定理，从原子轨道来构成一组正交基函数。3键轨道模型在某些情况下，用沿键方向dj的杂化轨道要比一般的s,p轨道方便得多。这些所谓的杂化轨道，是由s及p轨道线性组合而成，称为sp杂化。容易证明，沿方向dj的杂化轨道可以写为r小2-1(4.2.20)111.dir.1hdjr=石7777%Rsr广7T7Rpr-14n3+731|dj|r|式中R(r网Rp(r的s,p轨道波函数的径向部分。由于同一原子的轨道波函数是正交归一化的，因

9、此杂化轨道波函数也是正交归一化的。参数"由正交归一化条件得到一1党didj上,、，、+31厂i,j=1,2,3,4(4.2.21)(3+入(3%)|di|dj'如果一个原子的四个键沿正四面体棱线分布，两个键之间的夹角等于109吃8',%=1(j=1,2,3,4),我们得到所谓的sp3杂化轨道，对IV族半导体，如金刚石，Si,Ge,a-Sn就是这种情况。即沿口11方向沿力彳彳方向(4.2.22),It.,t,h4)=2Is-Px)-Py)+|如果d2,d3和d4在同一平面上且相互之间的夹角为|.7=0120°,而di垂直于这个平面，则沿d方向的轨道是纯粹的p轨

10、道，其余三个轨道是所谓的如果d2,d3和d4相互垂直，则有i22SP2杂化轨道，石墨就是这种情况。(4.2.23)(4.2.24)三个轨道为纯粹的需要注意的是,P轨道，第四个轨道为s轨道。杂化轨道并不是本征态，它对应的能量称为杂化能，记作(4.2.25)这里(4.2.26)对极性半导体来讲,定义=(PxH|Px)=(PyHPy)=(pzHPz)A原子与B原子是不同的，所以有两种不同的杂化能，能量低的记为(4.2.27)耳，高的记为名：。V3=；h-f2为杂化极性能(hybridpolarenergy)。对IV族半导体，如金刚石,金属能Vi(4.2.28)Si,Ge,a-Sn来讲=0。定义一个V

11、i=_hj(4.2.29)显然，金属能决定于sp分裂能级的大小(Sp-己)。对极性半导体来讲，两种原子的金属能是不同的。不同原子间的哈密顿量矩阵元，称为杂化共价能(hybridcovalentenergy)V22；fNs二23Vsp二3Vpp4(4.2.30)将两组杂化轨道线性组合成键轨道(bondorbital):W)=Uihi)+u2h2J(4.2.31)假设两组杂化轨道是正交的而忽略了轨道交叠，即hih2=0(4.2.32)这样能量期待值为W|H|，uM-2uiu2V2城«运用能量极小原理，即对上式求极小，得到方程(叩)(4.2.33)(4.2.34)为了方便起见，定义平均杂化

12、能;方程(4.2.34)变成解方程得到其中的一个解称为成键能：相应的键轨道参数u1,u2由下式给出:在求U1和U2时，用了归一化条件而%是键极化率，定义为这样就可以构成所谓的成键轨道：另外一个解称为反键能：相应的键轨道参数U1,U2为相应的反键轨道为:Ui-V2U2=Eui2-V2U1-ahU2=EU2(4.2.35)-V3U1-V2U2=E-7Ui-V2U1V3U2=E-二U21ip(4.2.36)(4.2.37)(4.2.38),1：2(4.2.39)YpI(4.2.40)221+U2=1(4.2.41)(4.2.42)(4.2.43)(4.2.44)(4.2.45)4原子轨道正交化线性组

13、合方法在经马的TB方法中，交叠积分或者作为参数，或者被忽略掉。为了得到更精确的能带，交叠积分需要作计算。在上两节的TB方法中，我们只考虑了价层电子，没有考虑芯电子。下面介绍第一性原理的原子轨道正交化线性组合(orthogonalizedlinearcombinationsofatomicorbitals,OLCAO)方法。这个方法同用平面波作基函数的第一性原理方法相比，具有计算量较小的优点，因此可以处理较多原子体系。在局域密度近似(localdensityapproximation,LDA)下多电子体系的能量由电子数密度P(r)唯一确定，电子数密度对应的能量为体系的本基态能量。这样可以导出如下

14、的自洽方程：1.22，"rVeerVxc：个4rg：nr2:r卜nrocc(4.2.46)(4.2.47)这里VeN,Vee分别代表电子与原子核，电子与电子库仑相互作用；Vxc为有效局域势,电子数密度的求和是对所以占据态进行。在上式中已经用了约化单位A72m=1通常称为交换关联势。LDA假设交换关联部分可以写为Excr=dr!:,r；xc-'j：r(4.2.48)因此，交换关联势可以通过以下变分求得：Vxcri：31：!v.xcRr)(4.2.49)由于不知道acP(rg对不同系统的严格形式，Vxc(r)可有不同的近似。这里用X口的交换关联，即3一3T3以,)=展二5J。(4

15、.2.50)其中a是参数，对不同的原子，其值介于1:2/3之间。在OLCAO方法下的布洛赫函数由原子或类原子轨道。(r)组成，即这里i代表轨道的对称性,¥表示单位原胞里的不同原子，q是v原胞的丁原子的位矢,(4.2.51)Rv是晶格位矢。原子轨道的径向部分用高斯型轨道(Gaussian-typeorbitals,GTO)展开：ir-Cjn%xp-1/2丫皿(4.2.52).j其中n代表主量子数，l和m代表角量子数和磁量子数，它们统一用i来表示。上式的原子轨道可以包括芯轨道、价轨道和空轨道。例如Cu2O,最少基矢包括Cu的1s,2s,2px,2py,2pz,3s,3px,3py,3Pz

16、(芯轨道);Cu的4s,4px,4py,4pz,3dxy,3dyz,3dzx,3d*24皿2d(价轨道)；。的1s(芯轨道);。的2s,2px,2py,2Pz(价轨道)；也可以考虑Cu的5s,5px,5py,5pz,4dxy,4dyz,4dzx,4dx2T2,4d3z22及。的3s,3px,3py,3Pz空轨道，构成所谓的全基矢。在一般情况下，用最少基矢可以得到比较好的结果。全基矢用于精度要求高的计算，如总能计算或那些牵涉到高能量导带的物理量的计算。单电子的波函数可以由布洛赫函数展开：(4.2.52)':nr,Aikbik,ri,这里n代表能带编号。这样我们可以得到如下的本征方程：H&

17、quot;.kkE=0(4.2.53)其中(4.2.54)(4.2.55)(4.2.56)Hi即)=(b小,r)Hda(k,r=£exp(TkR'©(rTy)“。川(r)+Vxc(r)由(r%Rv)v.Si,jk=bik,rbjk,r二'exp：-ikRv.；rr-Rv分别为哈密顿量和交叠积分矩阵元。库仑势为Vcoulr=VeNrVeer总的晶格态密度及库仑势用高斯函数表示为cryr八，-AJ.a,AVcoulVXC(r)=ZVC(r-Ta),r户',5r-.a,;ar='Cjexp：；：?."2VCr=ZaJexp：；T2:Dje

18、xp1r2jVxFjexp-r2j(4.2.57)(4.2.58)(4.2.59)这里Za为原子A的原子序数自洽是这样来实现的：解单电子本征方程，计算电荷密度及势能；再将新的势能带入本征方程计算电荷密度及势能，直到输入和输出的势能差异小于预定值。5自旋轨道相互作用对于重原子，自旋-轨道相互作用会使本来简并的原子能级发生分裂。分裂。前面的讨论中没有考虑自旋-轨道相互作用。这个相互作用的哈密顿量为对固体来说，则是简并的能级Hso=工(|A,i,。,R)2AAIa5(A,j,b,RB,i,ff,R+d)2%52,仃',R+dIt(4.2.60)R,Crtii,j这里R是A原子晶格矢量，d=q

19、/4111是B晶格相对于A晶格之位矢，i和j代表s,px,py,pz轨道，仃代表自旋(向上或向下J),A,i,cr,R)表示位于R处且具有自旋仃的A原子轨道波函数。因为s轨道的角动量l=0,因此只有处于同一原子且自旋相反的两个p轨道发生耦合。133%是p能级的重整化原子自旋-轨道分裂大小。如果考虑自旋-轨道相互作用,动量l和自旋角动量s之和：j=ls轨道波函数|i，m)(m=-i,i)不再是本征态，因而用总角动量相对应的波函数1j=i2，m-j,j可以将总角动量态用轨道角动量态和自旋态表示。容易证明，(4.2.61)则总的角动量j为轨道角(4.2.62)j,mj更为方便。这里(4.2.63)3

20、32,2户旧卜312,2ipx，(4.2.64)px,ti3sp3作为基轨道，则相应的哈密顿用这些新的基函数组成布洛赫和，可以求得新的哈密顿矩阵。如果用量矩B16x16矩阵6简单情况如果晶体中原子之间的间距比较大，每个原子的势场对电子有较强的束缚作用，电子距离位于某格点Rm原子比较近时，将主要受到该原子势场V(r_RmF勺作用，电子的行为同孤立原子中的电子行为相似。孤立原子中电子的薛定调方程、本征波函数(r-Rm心口本征能量6作为晶体中布洛赫电子的零级近似。把其它原子势场的影响H'=V(r)_V(r_Rm那成是微扰作用。为了简单，假设晶体是由基元中只有一个原子的N个格点组成，环绕N个不

21、同的格点，将有N个类似的波函数，它们具有相同的能量名，问题变成了N重简并的微扰问题，微扰后的波函数是N个简并态波函数的线性组合：r='、'am：r-Rm(4.2.65)m此处求和式遍及所有的格点。因此这种处理方法也称为原子轨道线性组合法，简称LCAOo考虑到微扰后的波函数应具有布洛赫波函数的性质，可以取1 /2ikRmam=Nem(4.2.66)则微扰后的波函数为：kr=N/2eikRmr-Rm(4.2.67)m不难证明这个波函数具有布洛赫函数的性质：,krT=N1/2%'eikRm；rT-Rm=N，/2eik八'eik*Rm；Pr-Rm-T=eikT-；r(4

22、.2.68)mm一由此波函数可以得到N个不同的能量本征值，原来一个原子能级分裂成一个由N个不同能级组成的能带,由此可以得到原子能级和晶体中能带之间的相互联系。按照简并微扰的计算方法，能量一级近似为：(4.2.69)(4.2.70)kH'k=N-:eik/j-jH''mj其中叼三p(rRj、令Pj=RjRmkH'k-'ek：jdVr-；?jH':r除了在属于同一个原子上和在被连结起来的最近邻原子之间的积分外，忽略所有其它的积分。令dV:rH'：r-(4.2.71)dVr：jH':r=一(4.2.72)于是kH'ke'

23、k心(4.2.73)考虑到一级近似布洛赫电子的能量为E=；kH'k_一_'、,PjekL：j(4.2.74)代入紧束缚近似的能量表达式得到:紧束缚近似对原子的内层电子相当不错,常常用来研究s电子形成的s能带。例如简单立方晶格体,每个原子的最近邻原子有6个:(a,0,0),(0,a,0),(0,0,a)(-a,0,0),(0,-a,0),(0,0,同(4.2.75)Ek=；_-2coskxa-coskya-coskza在第一布里渊区的r点k=p0,0)、X点k=9,0,0抽R点kE;，；的能量分别为：Er=8ot6y,aaaa12工EX=ea2¥和ER=8-(/+6了。

24、因为¥>0,1点和R点分别对应带底和带顶，因此能带宽度为由勺大小主要决定于最近邻原子波函数的相互重叠，重叠愈多，形成的能带愈宽。外层价电子波函数交叠程度大，形成的能带宽，内层电子波函数交叠程度小，能带比较窄对于体心立方晶体，最近邻有aaa一.一I.2,2,2aaaI一.II222,a-2a-28个原子:、1.11.a-2a-2一一,a2a-2a2a-2-代人紧束缚近似的能量表达式得到:kxakyakza(4.2.76)(4.2.77)能带中能量的极大值在H点：±,0,0；,a0,2二c十,01,acc2二0,0,-a能带中能量极小值在第一布里渊区中心(0,00)：Em

25、in；-1-8(4.2.78)能带宽度为16工在能带底附近余弦函数展开至二次项得:J1i'k、-8/1-1212.a1性访.JL一2=Emina2k：ky2kz=Emin而k：k2kZ(4.2.79)能带底部电子的有效质量(4.2.80)由此可知，波函数交叠越大，有效质量越小。同样可以得到带顶电子的有效质量为:一2mt2a2(4.2.81)此有效质量小于零。如果此能带近于被电子填满，则在能带顶部有空穴存在，空穴的有效质量为mh=mt按照LCAO近似，Si的能带中的电子态是这出组合系数和能量本征值。成键态对应的8个布洛赫和的线性组合。确定波函数后代入薛定调方程解4个能带交叠在一起，形成S

26、i晶体的价带；反键态对应的能带交叠在一起，形成Si晶体的导带。例题4.2.1用紧束缚近似求简单立方、体心立方和面心立方晶体解：E-：、je4kGs态原子能级相又t应的能带ES(k)函数。简单立方有6个最近邻原子。a,0,00,a,00,0,a-a,0,00,-a,00,0,-aEsk=1e-xa在带底k=0ikxaJkyaikyaikzaikzaee_ee-e=；_：_2cosKacoskyacoskza在带顶k=二aEsk6Esk-；6因此能带宽度为12黑旦:三_12=_6a2k2222.2电子的有效质量为：me=-2 a2-kza在带顶附近k=-kx,-ky,aa电子的有效质量为：m；-2

27、2a2体心立方有8个最近邻原子aaa丐kx*ykz2,2,2kx-kykzI-kxky-kz-i-kx""ky""kz-e2-e2-e2-i-a-kxkykz_i：-kx-kykz：ky«|i：-kx-ky-kz：l+e2+e2+e2+e2(=;一二_ekyakya11_L-i22e+eKai2.kzaie2JJ-8coskacoskyacoska在带底k=0在带顶,0,0aEsk=8因此能带宽度为16Z在带底附近k0kxa8电子的有效质量为:J-22ha8_;-8;r+;1a2k2在带顶附近k=me=22a212Ji.-kx,ky，akzEs(kI-a九?1k旦iI8人(kya),8,一七a2k2电子的有效质量为：me-2<,22a面心立方有12个最近邻原子aaa2-1,-1,020,-1,-12-1,0,-1.s.Ek-;-"aaaaaa_e2k、ky.e2k、K4ygy咋2ky上I+e/长)+1(上)+e2(")+e21晨/)+e»*IJY4i丝Y上匕”虫iy占上)/上上y4&=o(一¥(，2+e2e2+e2+|e2+e2e2+e2+e2+e2e2+e2IA