体心立方晶格紧束缚近似能带结构的计算机模拟.doc

体心立方晶近似 肖瑞春陶松涛摘 要:利用MATLAB晶在紧束缚近似下s态进行计算机模拟布里渊区内等能面的清晰图像关键字：立方等能面 能带理论电子能量色散关系能等能面等能面的表达式比较复杂,几何结构很难想象.de Haas-van Alphen 效应的实验数据，绘制了一些金属的费米面的三维图形[1]，受到广泛关注。随着计算技术的发展，人们开发出了许多卓越的分析软件，便于研究涉及大量数据的复杂问题。借助于这软件，人们开始了晶格能MATLAB软件，对紧束缚近似下晶s态布里渊区内等能面清晰的图像紧束缚近似立方晶s态能带 紧束缚近似下i态原子能级形成的能带为 （1） 其中为孤立i的、是重叠积分s态个近邻相同，记为，有（1）式可得: （2） 其中a为晶格常数，为波矢量. 根据Bloch定理可推知，晶体的电子能带具有周期性，即只要研究清楚一个倒格子原胞（取简约布里渊区）的情况即可. 体心立方晶格的倒格子为面心立方格子（单胞边长）。根据布里渊区的界面方程 （是倒格矢） （3） 利用MATLAB可作出体心立方晶格体心立方晶格、、的特定的方向，可以表示为： （4） 这就是该方向的能带曲线方程。其中： 沿轴: 沿轴: 沿轴: 2利用MATLAB模拟能带的设计思路及步骤[] 2. 1等能面与能带曲线方程的无量纲化 表达式 （5） 其中 （6） （7） 对于（4）式，记 ，则有 （8） 于是，沿轴、沿轴、沿轴的能带曲线方程分别为： （9） 2. 2 1)构建三维数组，生成三维网格 [x,y,z]=meshgrid(linspace(-2,2,n)); 其中n为linspace生成一维数组元素个数，已事先定义 2）MATLAB里（*）为数组运算，数组元素的运算用点乘（.*）来实现。首先将（5）式转化MATLAB符号的表达式为： val=1-(cos(0.5*pi*x).*cos(0.5*pi*y).*cos(0.5*pi*z) Val即为()里的 3)限制图形绘制范围若不对x、y、z进行限制，以上等能面将在的在倒格子单胞里绘制，故要把布里渊区以外的图形切掉 可以通过逻辑索引的方式来实现，具体思路是：通过逻辑运算与关系运算产生一个索引数组（实际为一个由0、1组成的关系数组），然后利用这个索引数组来访问原数组，并进行重新赋值。 w=x+y-2>0|x-y-2>0|-x-y-2>0|-x+y-2>0|x+z-2>0|-x+z-2>0|y+z-2>0|-y+z-2>0|x-z-2>0|-x-z-2>0|y-z-2>0|-y-z-2>0; val(w)=NaN; 该命令以布里渊区边界面方程为判读语句产生一个索引数组w，并用w访问原数组将布里渊区以外的val赋值为NaN（非数），从而实现了三维图形的切割。 4)用isosurface、patch、set命令绘制等能面图像，并对图形属性进行设置绘制布里渊区轮廓在绘制前将布里渊区顶点坐标和各多边形片块的顶点索引并编成数组. 用以下绘图命令实现 patch('Vertices',Vm1,'Faces',Fm1,'Linestyle','-','LineWidth',2,'FaceColor','none'); 其中V1数组是布里渊区顶点坐标，Fm1为布里渊区顶点索引. 对视角、光照效果、坐标轴进行设置. 模拟绘图结果 s态能带的等能面在简约布里渊区的形态结构及其变化的清晰图像，如图2.1—2. 12所示. 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 图2. 紧束缚近似下体心立方晶格的等能面 图3 不同等能面的嵌套和二维视图 4. 几点讨论 从模拟结果可以看出 （1）简约布里渊区等能面在能量取值较小时与自由电子等能面相似 （2）随着电子能量增加，等能面的形状的正方体表面； （4）当等能面抵达里渊区里渊区 （5）接近能带顶部，等能面与半导体材料的费米面相差较大[5]. 实际上，这时公式（3）或公式（5）以不适用，通过修正能带公式，也可以得到与实际相符的结果，本文不做进一步分析； （6）利用MATLAB图形旋转功能可以使三维旋转，这就使得这些抽象的概能够形象、生动地展现出来 参考文献: [1] D. Shoenberg，The de Haas-van Alphen E