VASP入门到精通[32]态密度和HOMO-LUMO

这一节小编写了个小脚本，提取DOS和HOMO-LUMO能级

工具/原料

vasplinux

方法/步骤

1、注意：这2个脚本是针对ISPIN=2的计算！！！！！！！！！！！！！小编写了两个脚本，一个名字叫data-dos一个名字叫gap.sh

2、第一个脚本是用来提取态密度数据的，脚本命令如下a=`head-6DOSCAR|tail咯悝滩镞-1|awk'{print$3}'`b=$((a敫苻匈酃+6))f=`awk'{if(NR==6)print$4}'DOSCAR`sed-n'7,'$b'p'DOSCAR>DOS.datawk'{print$1-'$f',$2}'DOS.dat>DOS-final.dat

3、注释：第一行：提取DOS数（就是取了多少个点）第二行：DOS的终止行，就是我们要提金泫险缋取a到b行的数据第三行：提取锗咱貉沲费米能级第四行：提取a到b行的DOS数据并写入DOS.dat（这时会生成一个DOS.dat文件）第五行：把提取的DOS.dat做费米修正，输出一个DOS-final.dat这个DOS-final.dat就是你要的DOS数据，放到origin就可以看到了

4、第一个脚本的使用方法：运行：bashdata-dos你会发现多了两个文件:DOS.dat和DOS-final.dat把DOS-final.dat放到origin就可以了

5、第二个脚本代码如下：m=`awk'{if($1<0&&(铮喘饿怊$2>0.05||$3>0.05))a=$1}END{print惺绅寨瞀a}'DOS-final.dat`n=`awk'{if($1>0&&($2>0.05||$3>0.05))print$1}'DOS-final.dat|head-1`echo“HOMO=$m”echo“LUMO=$n”echo"$n"-"$m"|bc

6、注意这个脚本就五行哦，小编有截图，怕在这里编辑就自动换行，所以截图给大家

7、注释：第一行提取homo第二行是LUMO第三四五行就输出HOMO，LUMO和能隙了

8、使用方法：bashgap.sh你会看到输出三行内容：HOMO=LUMO=能隙

9、有问题可以留言问小编哦！！