紧束缚应用实例.docx

应用实例（1）紧束缚模型求解六角格子的色散关系和态密度 对于六角格子，现以石墨烯为例。石墨烯是由碳原子以sp3杂化形成的二维六角蜂窝状结构，理论上认为石墨烯不可能二维存在。2004年，英国曼彻斯特大学的Geim首先用微机械剥离法制备了石墨烯，并与2010年获得诺贝尔物理学奖，引起人们的广泛关注。不再像正方、三角格子是简单格子，六角格子是由A、B两套子格套构而成的复式格子。对应的紧束缚模型哈密顿量为：H=-ti,ai+bi+hc-ti,ai+ai+hc-ti,bi+bi+hc-iai+ai-ibi+b (3-1)式中，i和j分别代表六角格子中最近邻和次近邻的格点，ai和ai+分别代表在a子格i格点电子的湮灭和产生算符。做傅里叶变换，可得紧束缚模型的哈密顿量为：H=-tk,eikak+bk+hc-2tk,coskak+ak+bk+bk-kak+ak+bk+bk (3-2) 可令上式中-teik=-tcosk+isink =xk+iyk-2tk,cosk=k (3-3) 对于六角格子，上式中的的坐标为0，1和（32，12），的坐标为3，0和（32，32）。则 xk=-tcos3kx2+ky2+cos3kx2-ky2+cosky (3-4) xy=-tsin3kx2+ky2+sin3kx2-ky2-sinky (3-5) k=-2tcos3kx2+ky2+cos3kx2-ky2+cos3kx (3-6) 最后哈密顿量可以简化为：H=kxk+iykak+bk+xk-iykbk+ak+(k-)ak+ak+bk+bk， (3-7) 可以进一步写成矩阵形式，即H=kak+,bk+k- xk+iyk xk-iyk k- akbk (3-8) 色散关系可通过下式得到：k-k xk+iykxk-iyk k-k =0 (3-9) 行列式可得k=k-xk2+yk2 (3-10) 将（3-4）-（3-6）式代入（3-10）式可得六角格子的色散关系为k=-2t2cos3ky2cos3kx2+cos3kx- t3+2cos3kx+4cos3ky2cos3kx2 (3-11)若只考虑电子在最近邻格点上的跳跃，色散关系为：k=t3+2cos3kx+4cos3ky2cos3kx2 (3-12) 以上结果是从二次量子化后的哈密顿量得到的，与从波函数角度出发理论计算的结果一致。若同时考虑电子在次近邻格点上的跳跃，则(3-11) 式中的t=0.2t。同时，在计算中都设=0。 用Matlab软件将六角晶格的色散关系表达式进行图形化，会考虑电子在近邻格点上的跳跃和考虑到电子在次近邻格点上的跳跃色散关系结果分别如图a和b所示。 六角格子的态密度表达式为：=1Nk-k=11Nk-12+2+-22+2 (3-13)式（3-13）中1和2分别对应于(3-11) 式中的k取正号和取负号的值，按照（3-13）式，用Fortran软件编程计算对应的态密度，然后用Origin作图，结果如图c和d所示。 由图7a可知，只考虑电子在最近邻格点上的跳跃，上、下两带是对称的，相交于第一布里渊区的6个顶点，也被称为Dirac点。Dirac点附近能带具有线性色散关系。能带对应的能量取值范围为-3t,3t，带宽为6t 。费米面刚好处于价带和导带相交的顶点处，可知，石墨烯是带隙为零的半导体。由图7c可以看出，在Dirac点附近，态密度具有线性关系，而且态密度关于E=0对称。如果考虑电子在次近邻格点上的跳跃，由图7d可知，上、下两带是不对称的，这是因为电子-空穴对称性被破坏。能带对应的取值范围为-4.2t,1.85t，带宽为6.05t。由图7d可以看出，价带和导带的态密度也不再关于Dirac点对称。这些结果与文献中的一致。 图七 六角格子的色散关系（2）紧束缚模型在DNA分子中的应用紧束缚模型方法的基础是紧束缚近似，即在分子或固体体系中，电子在某原子附近时，将主要受到该原子势场的作用，其它原子作用可以看作微扰，根据这个近似思想，紧束缚方法可以计算周期性晶格中的电子能带结构、能量、力及其它性质。紧束缚模型方法是紧束缚方法在简单晶格体系中的应用。该模型把晶胞抽象为格点，每个格点提供一个电子轨道，并且仅考虑邻近格点间的相互作用，紧束缚模型的哈密顿为：HTB=nncn+cn-tn,n+1cn+1+cn+cn+cn+1，式中c和c+分别为电子（或空穴）的湮灭算符和产生算符，表示个点的在位能，t为格点间电荷迁移积分。这一模型最初用于聚乙炔等导电聚合物分子，碳原子被视为格点，其上的2pz轨道为格点轨道，类似于量子化学中的休克尔方法。1979年，Su等人对聚乙炔的紧束缚模型做了改进，建立了著名的Su-Schrieffer-Heeger（SSH）模型。确定合适的在位能和迁移积分后，紧束缚模型方法可用于计算DNA、蛋白质等大生物分子的电子结构，研究其电荷输运性质。生物大分子中多电子体系的薛定谔方程为：H=E 其中，多电子的哈密顿为： H=i-12i2-a,iZarai+ij1rij+abZaZbRab 式中第一项为电子动能，第二项为核对电子的吸引能，第三项是电子相互作用能，第四项是核间排斥能，指标i和j对电子求和，a和b对原子核求和。生物大分子都是由一些固定的结构单元聚合而成的，例如蛋白质的结构单元是氨基酸，DNA的结构单元是核苷酸。根据这个特点，把结构单元视为格点，可将生物大分子抽象为格点模型。根据紧束缚近似的思想，电子主要受格点势场的作用，因而将局域在格点指代结构单元内。因此可将多电子哈密顿写为以下形式：H=i-12i2-i,LaLZarai+12i,LjL,ji1rij+La,bL,abZaZbRab+L,MaL,bMZaZbRab其中L、M指代格点。暂时略去对电子结构无影响的和排斥能项，上式可写为单电子算符加和的形式即：H=i-12i2+LVL(i) 其中， VLi=-aLZarai+12i,LjL,ji1rij 相应地，多电子薛定谔方程的求解转化为单电子的求解，h=其中为生物大分子的单电子波函数，单电子的哈密顿为h=-12i2+LVL，根据紧束缚近似，电子在格点附近时，主要受该格点势场 的影响，其他格点势场可作为微扰，求解方程应先求解无微扰时即孤立格点对应的单电子薛定谔方程，即-12i2+LVLl=l0l其中l是格点L方上的格点轨道，l0为相应的轨道能。生物大分子整体的电子轨道可展开为各格点轨道的线性组合，即=lcll代入h=可得以下线性方程：lclmhl=lml根据紧束缚近似，略去非近邻作用项，得mhmcm+m-1hmcm-1+mhm+1cm+1=cm该式与紧束缚模型哈密顿量HTB=nncn+cn-tn,n+1cn+1+cn+cn+cn+1相对应，可得在位能和迁移积分的计算公式分别为 n=nhn=n-122+LVLn tn,n+1=-nhn+1=-n-122+LVLn+1以上两式中的求和项L应走遍所有的格点，但是根据紧束缚近似，可只考虑近邻格点的影响。所以当计算格点n的在位能时，仅需考虑格点n及其两侧格点n-1与n+1的影响（线性分子），相应地在位能计算公式简化为nn-122+Vn+Vn-1+Vn+1n=n0+nVn-1+Vn+1n从该式可以看出格点的在位能并不等于格点上的轨道能量，因为还应该考虑至少两侧格点的影响，而对于格点n和n+1之间的迁移积分，仅需要考虑这两个格点势场的影响，即tn,n+1tn,n+1+12sn,n+1(n+n+1) 式中s为格点轨道间的重叠积分。而这种正交变换对格点在位能的影响很小，可忽略不计。这样，我们得出了紧束缚模型哈密顿中在位能和迁