紧束缚应用实例.docxVIP

应用实例 （1）紧束缚模型求解六角格子的色散关系和态密度 对于六角格子，现以石墨烯为例。石墨烯是由碳原子以sp3杂化形成的二维六角蜂窝状结构，理论上认为石墨烯不可能二维存在。2004年，英国曼彻斯特大学的Geim首先用微机械剥离法制备了石墨烯，并与 不再像正方、三角格子是简单格子，六角格子是由A、B两套子格套构而成的复式格子。对应的紧束缚模型哈密顿量为： H=-ti,δa 式中，i和j分别代表六角格子中最近邻和次近邻的格点，ai和ai+分别代表在a H=-t μkak+ 可令上式中-tδe =xk+i 对于六角格子，上式中的δ的坐标为0，1和（±32，1 xk=-t xy=-tsin δk=-2t 最后哈密顿量可以简化为： H=kxk+i 可以进一步写成矩阵形式，即 H=kak+ 色散关系可通过下式得到： δk-μ-εk 行列式可得εk=δk-μ±xk2+yk2 (3-10) 将（3-4）-（ ε t3+2cos3 若只考虑电子在最近邻格点上的跳跃，色散关系为： εk=±t 以上结果是从二次量子化后的哈密顿量得到的，与从波函数角度出发理论计算的结果一致。若同时考虑电子在次近邻格点上的跳跃，则(3-11) 式中的t=0.2t。同时，在计算中都设 用Matlab软件将六角晶格的色散关系表达式进行图形化，会考虑电子在近邻格点上的跳跃和考虑到电子在次近邻格点上的跳跃色散关系结果分别如图a和b所示。 六角格子的态密度表达式为： ρε 式（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=nεncn+cn-tn,n+1cn+1+cn+cn+ 生物大分子中多电子体系的薛定谔方程为：HΨ=EΨ 其中，多电子的哈密顿为： H= 式中第一项为电子动能，第二项为核对电子的吸引能，第三项是电子相互作用能，第四项是核间排斥能，指标i和j对电子求和，a和b对原子核求和。生物大分子都是由一些固定的结构单元聚合而成的，例如蛋白质的结构单元是氨基酸，DNA的结构单元是核苷酸。根据这个特点，把结构单元视为格点，可将生物大分子抽象为格点模型。根据紧束缚近似的思想，电子主要受格点势场的作用，因而将局域在格点指代结构单元内。因此可将多电子哈密顿写为以下形式： H=i-12? H=i-12 相应地，多电子薛定谔方程的求解转化为单电子的求解，hψ=εψ其中ψ为生物大分子的单电子波函数，单电子的哈密顿为h=-12?i2+LVL，根据紧束缚近似，电子在格点附近时，主要受该格点势场 的影响，其他格点势场可作为微扰，求解方程应先求解无微扰时即孤立格点对应的单电子薛定谔方程，即-12?i2+ 生物大分子整体的电子轨道可展开为各格点轨道的线性组合，即ψ=lcl?l代入hψ=εψ可得以下线性方程：l εn 以上两式中的求和项L应走遍所有的格点，但是根据紧束缚近似，可只考虑近邻格点的影响。所以当计算格点n的在位能时，