dos处理及能带做图说明分析.docxVIP

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?里面的三列数字复制到此前三列中，如图： 接下来直接画图就可以直接得到能带图！