dos处理及能带做图说明.doc

DOS处理： 在dos计算文件夹中直接使用命令split_dos 处理得到DOS0 DOS1 DOS2 DOS3……分别对应总态密度，各个原子的态密度（一一对应POSCAR原子顺序）。处理后的DOS0，DOS1….等文件的能量值是以费米能级作为能量参考点，DOS0的第一列数据是能量值，单位为eV,第二列数据是总态密度的值，单位为electrons. DOS1是第一个原子的分波态密度值，其中第一列数据是能量值，单位为eV。第二三四列分别对应s, p, d, 态的分波态密度值，单位为State/eV.atom,其他文件类似DOS1. Sum_dos(sum_dos_np)使用说明： Sum_dos使用于自旋计算，sum_dos_np对应非自旋计算。他们是把同一类原子的dos对应态密度值加起来！以GaN计算为例，若计算结果里有4个N 4个Ga（对应POSCAR里原子顺序为4个N在前，4个Ga在后），那么在dos计算出来split_dos处理后会出现DOS0 DOS1——DOS8，其中DOS1——DOS4对应的是4个N原子态密度，DOS5——DOS8对应是4个Ga态密度。若为自旋计算，利用命令sum_dos 0 1 4处理得到DOS.SUM.1.to.4; DOS.SUM.5.to.8两个文件，分别为4个N（Ga）的s p d 态各自相加后的总值（里面每一列对应4个s态的和，4个p态的和，4个d态的和及s+p+d总和）。非自旋计算利用sum_dos_np同样方式可得。 BAND画图处理：（以GaN为例） 布里渊区路径选择为G-K-M-G如图所示： 计算过程中两高对称点间插入30个点，首先测出G-K-M-G各段距离：在MS中选定一个线段左侧就能给出此线段距离如图所示： 选择G-K段后此线段变为黄色，左侧显示出长度0.206490.同理可得到K-M M-G各段长度分为0.103245，0.179004。 然后新建EXCLE文件（此过程相当于把三段长度分别扩大30倍）： 如图所示：，然后左键点击0不放松下拉到0.20649 选定1，2行，变为，然后鼠标放到此黑框右下角位置（鼠标变为细黑体十字）点击下拉30行如图： 第31行与30行同一个数如图：。然后第32行数为5.98821+0.103245（K-M线段长）=6.091455，接下来同理选定30，31行下拉到60行，如图所示： 61行与60行相同，以后重复上面步骤加上M-G线段值下拉到90行，如图： EXCLE里面第二列在30，31（60，61）行写入-10 10（对应画图的能量范围） 第三列1-90行都写入0（对应画出图的费米能级），如图所示： 最后就是origin画图：把bandplot处理得到的band.dat导入origin，选定第二列，右键选择Insert添加一列空列如图：一共添加两列空列， 然后选定前三列（一直按着Ctrl键用鼠标点前三列）如图 接下来把EXCLE里面的三列数字复制到此前三列中，如图： 接下来直接画图就可以直接得到能带图！