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

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

2、第一个脚本是用来提取态密度数据的，脚本命令如下a=`head -6 DOSCAR|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

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、使用方法：bash gap.sh你会看到输出三行内容：HOMO=LUMO=能隙

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