dos处理及能带做图说明

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